Discussion :

[louftansa]

Bonjour à tous,

Ceci est un post que j'ai déja publié mais qui a été supprimé à cause des problèmes de serveur. Je dispose de deux matrices comprenant chacune une colonne pour les coordonnées des points et une autre pour leurs images respectives. Pour chaque matrice, j'ai tracé la courbe correspondante. Maintenant je voudrais calculer l'erreur entre ces deux courbes. Quelqu'un aurait il une idée de comment faire ??? . Mes fichiers de données sont en pieces jointes ci-dessous.

Merci de vos réponses.

[mr_samurai]

Salut,

Tu peux voir -ici-

++

[casier]

il a jamais dit qu'il s'agissait de variables aléatoires, pourquoi aller chercher une erreur quadratique ?

Si tu veux lla différence entre tes deux tracés, il suffit simplement de faire la différence des colonnes d'images :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
t=A(:,1);   % vecteur temps, pour le tracé
erreur=A(:,2)-B(:,2); % vecteur d'erreur, ajoute un abs(...) éventuellement
plot(t,erreur)

[louftansa]

Effectivement, il ne s'agit pas de variables aléatoires. En ce qui concerne le fait de pouvoir faire la différence entre les deux colonnes image, je ne pense pas que ce soit la bonne solution car je n'ai pas les mêmes colonnes coordonnées pour mes 2 courbes...

[Hanz Peter Mariol]

"L'erreur quadratique d'un estimateur par rapport à un modèle est l'espérance du carré des différences entre les observations et les valeurs prédites par le modèle."
L'erreur quadratique sur un couple de point est le carré de l'opération précisée par mon précédent. Une fois ce carré effectué, tu fais la moyenne sur tous les points et tu as l'erreur quadratique moyenne (ou variance, il me semble, en probabilité).
Si tu prends la racine carré de ta moyenne quadratique, tu as ce qu'on appel l'écart type en proba. Mais dans ton cas ce sera donc l'écart moyen entre un point de la courbe A et un de la courbe B.
Si j'ai compris le problème et la question, j'espère y avoir répondu ^^


Bonne journée

[Hanz Peter Mariol]

Je viens de voir ton dernier post...
Après meme si t'as pas les meme coordéonnées, tu peux "bidouiller". Si les echelles sont les memes, tu fait des approximations (linéaires ou autres) pour obtenir tes valeurs à des coordonnées identiques.

[mr_samurai]

salut casier,

Je ne vois pas la relation entre l'erreur quadratique et les variables aléatoires . D'ailleurs, c'est l'erreur quadratique qui est utiliser (minimiser) dans la régression linéaire (Quadratique = Moindre Carré) ...

Bref, ce que tu propose n'est pas une évaluation d'un erreur de proximité, car ce n'est pas une 'Distance' sur l'espace concerné.

[Hanz Peter Mariol]

En fait la premiere réponse partait du postulat que les abscisses des points pris pour les deux courbes étaient les memes, d'ou sont rejet de la méthode qui fut appelé à tord "erreur quadratique".
Mais sinon, c'est bien les moindre coarré qu'il faut utiliser^^

[mr_samurai]

rebonjour,

Voici un exemple d'interpolation et de calcul d'erreur :

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
7
8
9
10
11
12
13
14
 
% Exemple
xx1 = 1:0.1:5;
xx2 = 1:0.13:5; % Pas Different!
 
% Calcul normal
A = sin(xx1); 
 
% Interpolation 
Atmp = sin(xx2);
Atest = interp1(xx2,Atmp,xx1,'cubic');
 
% Error : erreur d'interpolation
Err = mean((A-Atest).^2)

++ bonne chance

[casier]

Ah oui si les axes de temps ne sont pas les même là faut bidouiller un peu, interpolation etc., comme le propose mr_samurai

salut casier,

Je ne vois pas la relation entre l'erreur quadratique et les variables aléatoires .

Ben tu le renvoies sur un lien wiki qui parle d'erreur d'erreur quadratique moyenne, il s'agit bien de variables aléatoires (l'espérence d'une variable non aléatoire ça n'a pas trop d'intérêt ) (mais bien sur, c'est les moindres carrés qui sont utilisés pour le faire...)

A ce moment là parlons de moindre carrés et on est tous d'accord

[Hanz Peter Mariol]

Gloablement qu'est ce que l'espérance quand tu l'as calcul? C'est une moyenne... Et si tu regarde la fonction matlab utilisée par monkey d luffy(mr_samurai) : "Err = mean(...)". "Mean", qui en anglais signifie moyenne.
Ne te borne pas à une manière de décrire les choses, si la boulangère veut t'offrir ses miches, attends toi à avoir du pain...

[louftansa]

Merci de vos réponses. La proposition de Samurai marche bien. A plus...

[casier]

Gloablement qu'est ce que l'espérance quand tu l'as calcul? C'est une moyenne... Et si tu regarde la fonction matlab utilisée par monkey d luffy(mr_samurai) : "Err = mean(...)". "Mean", qui en anglais signifie moyenne.
Ne te borne pas à une manière de décrire les choses, si la boulangère veut t'offrir ses miches, attends toi à avoir du pain...

Heuu, ouais... Merci du conseil, mais je vais m'en tenir à la vision des choses (au passage je ne me suis bloqué nulle part moi, j'ai juste dit que le lien wiki n'allait pas beaucoup aider louftansa, je pense pas m'être trompé !). Une espérance c'est une espérance, une moyenne c'est une moyenne, si c'était la même chose on les aurait pas appelé différemment hein... (inutile de me sortir des théorèmes qui stipulent que patati et donc patata, je les connais moi aussi)

Parler dans ce topic d'espérance est un non sens, la moyenne qui est faite dans le code n'a rien à voir avec une quelconque espérance de quoi que ce soit (ah si... ce serait une approximation de l'E d'une certaine VA inconnue dont les valeurs données par louftansa seraient plusieurs tirages... la belle affaire...), d'ailleurs mr_samurai n'a jamais évoqué l'espérance dans ses messages, à partir de là je vois pas pourquoi le débat repart là-dessus... en tout cas moi j'en reste là

[Hanz Peter Mariol]

Certes, l'espérance est un outil des probabilités, cependant je pense qu'y trouver une certaine analogie avec un moyenne est loin d'etre un non sens...
Même si théoriquement, si l'on est strict et rigoureux(je n'ai pas dit rigide...), ce n'est pas exactement le bon terme, il me parait aisé de comprendre la solution avec l'expression "erreur quadratique moyenne".
Bref peu importe, notre camarade a la solution et nous avons deux manières différentes d'appréhender les choses.
Si jamais quelqu'un voulait continuer d'échanger des idées, n'hésitez pas à MP (ne floodons pas le topic), ce sera toujours enrichissant, Merci.
Bonne journée

[tatou006]

merciiiii infiniment mr_samurai