Bonjour a tous,

je veux ajuster mes valeurs expérimentaux a un model théorique en utilisant la méthode de moindre carrée sous matlab (le but du jeu est la comparaison des d#ts erreurs de mes valeurs exp. par rapport a chaque courbes théorique et enfin choisir celle qui possède une R a peu prés égale à 1.
j'ai déjà essaie mais je trouve pas les résultats attendus,est ce que vous pouvez me corriger mon programme SVP ,merci d'avance
Voila le code:

Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
clear all
close all
 
 
x=[ 0   0.5      1       1.5      2        2.5      3        3.5     4         4.5      5        5.5      6        6.5      7        7.5      8        8.5     9         9.5      10       10.5       11     11.5       12      12.5      13       13.5      14     14.5       15      15.5     16      16.5     17      17.5      18      18.5       19      19.5       20] ;
y1=[0 8.52e-84 1.95e-57 1.82e-39 4.00e-28 2.22e-21 8.27e-18 2.99e-16 1.12e-15 1.97e-15 2.70e-15 3.26e-15 3.68e-15 3.97e-15 4.16e-15 4.28e-15 4.36e-15 4.36e-15 4.40e-15 4.43e-15 4.45e-15 4.45e-15 4.46e-15 4.47e-15 4.47e-15 4.47e-15 4.47e-15 4.47e-15 4.47e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 ];
w= [0 5.87e-84 1.20e-57 3.06e-40 2.59e-28 1.10e-21 7.21e-18 2.21e-16 1.00e-15 2.07e-15 2.51e-15 3.17e-15 3.63e-15 4.03e-15 4.18e-15 4.25e-15 4.38e-15 4.38e-15 4.38e-15 4.38e-15 4.40e-15 4.40e-15 4.41e-15 4.42e-15 4.42e-15 4.42e-15 4.42e-15 4.42e-15 4.42e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15];
a= [0 1.82e-86 3.37e-66 2.12e-42 3.16e-31 1.23e-21 9.48e-18 1.90e-16 1.37e-15 2.50e-15 2.80e-15 3.20e-15 3.70e-15 4.10e-15 4.19e-15 4.30e-15 4.40e-15 4.40e-15 4.40e-15 4.40e-15 4.41e-15 4.42e-15 4.42e-15 4.42e-15 4.42e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15 4.44e-15];
e1=[0 8.87e-85 1.68e-58 5.26e-40 4.56e-28 3.56e-21 8.89e-18 3.25e-16 1.22e-15 1.95e-15 2.65e-15 3.32e-15 3.75e-15 4.05e-15 4.15e-15 4.25e-15 4.25e-15 4.35e-15 4.45e-15 4.45e-15 4.45e-15 4.48e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 4.50e-15 ];
e2=[0 7.45e-85 2.05e-57 2.26e-39 5.56e-29 2.89e-21 7.95e-18 3.99e-16 1.08e-15 1.58e-15 2.25e-15 3.24e-15 3.55e-15 4.00e-15 4.18e-15 4.26e-15 4.38e-15 4.38e-15 4.42e-15 4.42e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.43e-15 4.44e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 4.45e-15 ];
e3=[0 8.86e-83 2.20e-67 0.26e-39 4.23e-40 2.23e-21 8.23e-18 2.25e-16 1.25e-15 1.92e-15 2.90e-15 3.16e-15 3.65e-15 4.10e-15 4.20e-15 4.30e-15 4.39e-15 4.37e-15 4.42e-15 4.42e-15 4.47e-15 4.47e-15 4.48e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.49e-15 4.50e-15 4.48e-15 4.49e-15 4.51e-15 4.51e-15 4.50e-15 4.50e-15 4.50e-15 4.51e-15 4.49e-15 4.49e-15 4.50e-15 ];
e4=[0 6.58e-84 1.42e-57 2.44e-40 3.65e-28 1.36e-21 9.35e-18 1.45e-16 0.89e-15 2.03e-15 2.65e-15 3.28e-15 3.70e-15 3.98e-15 4.145e-15 4.29e-15 4.34e-15 4.34e-15 4.41e-15 4.44e-15 4.46e-15 4.46e-15 4.47e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 4.48e-15 ];
 
 
 
plot(x,y1,x,w,x,a,x,e1,x,e2,x,e3,x,e4,'*') %les * sont les val exp 
R=y1-w;
R1=y1-a;
R2=y1-e1;
R3=y1-e1;
R4=y1-e1;