Méthode itérative par MATLAB
Bonjour,
je me trouve dans un impasse depuis quelque jours,, j'arrive pas à résoudre un problème que je dois très vite exposé.
Mon problème est d’estimer les coordonnées d'un hypocentre/épicentre d'un évènement sismique par méthode itérative en prenant l'exemple de stations rapprochées.
on m'a imposer le développement en série de Taylor pour linéarisation de ce problème que je le présente tel que:
T= Distance/Vitesse + T0
Tiobs = Di/Vi ; i=1à6.
Di = sqrt(((x1-X).^2)+((y1-Y).^2))
Di/V ; V=5km/s.
donc la formulation est:
Tiobs= sqrt(((x1-X).^2)+((y1-Y).^2))/5.
par méthode de moindre carrées ça marche le plus normalement du monde,, mais la méthode itérative,,,,????
voici l'algorithme que j'utilise et qu'est apparemment ça ne marche pas;
%calcule initial
t1=sqrt(((x1-X).^2)+((y1-Y).^2))/5;
m(1)=t1; m(2)=x1; m(3)=y1;
% le nombre d'étiration
for i=1:6
for i=1:length(Y)
diffx=(X(i)-m(2));
Diffy=(Y(i)-m(3));
D0(i)=sqrt(diffx^2+ Diffy^2);
d(i)=m(1)-D0(i)/5-m(1);
G(i,3) = -(Diffy/D0(i))/5; G(i,2)= -(diffx/D0(i))/5; G(i,1)= 1;
end
% sollution par moindre carré
A=G*G';
B=d*G';
B=A*m;
dm=A\B
end
svppppppppppppppppp quelqu'un pourrais m'aider là dessus?????
