Discussion :

[narcisseb]

Bonjour,
j'ai 25 images ciné_IRM qui représente le mouvement du cœur pendant le cycle cardiaque, à chaque fois j'obtiens la courbe de variation de l'intensité au cours du temps en suivant la variation de l'intensité du pixel au cours de ces 25 images
Ci-joint un exemple d'une courbe de variation de 'intensité au cours du temps pour un pixel donnée.
maintenant j'ai besoin de faire l'ajustement de cette courbe par une sinusoïde en utilisant la méthode des moindres carrées
la fonction de cette sinusoïde est: A1 sin (wt+φ1) + A2 sin (2wt + φ2)

j'ai essayé avec la fonction lsqcurvefit mais ca n'a pas marché
j'ai besoin d'aide svp

[narcisseb]

Bonjour
voici le code de la fonction de fitting :

Code :

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

1
2
3
4
5
function y =f(param,t)
a1= param(1);
a2=param(2)
w=0.2;
y = a1*sin(w*t + a2)  ;

et le code principal:

Code :

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

1
2
3
4
5
6
7
8
9
10
11
t=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25];
A= [910  773  656  548  477    422   414    412   398   382   366   352    380    514    630  666    811    923     1030    1052   1007   982   1023    1049  1041];
figure(2),
c=plot(t,A(t));
a0 = [500  120 ]
options = optimset('lsqcurvefit');
options = optimset(options,'LevenbergMarquardt' ,'on');
options.MaxFunEvals = 100;
options.MaxIter = 6;
asol = lsqcurvefit (@f, a0, t, A,[],[],options);
plot (t, A, '*', t, f(asol, t))

ce code me donne de mauvais résultats Ci-joint le résultat obtenu
SVP j'ai besoin de vos aides ?

[Invité]

Bonjour,

L'erreur vient du modèle que tu as choisi : a1*sin(w*t + a2).
Tu considères que ta courbe oscille autour de 0, or elle le fait plutôt aux environs de 800.
Il te faut au moins ajouter un offset : a1*sin(w*t + a2) + a3

[narcisseb]

je vous remercie énormément pour votre réponse
le problème est résolu grâce à votre remarque,
ci-joint la courbe obtenue.
permettez moi de poser une autre question,
Pour chaque pixel, j'ai une courbe de variation de l'intensité au cours du temps,
donc je dois faire le lissage pour toutes les courbes mais le problème que les paramètres de a0 varient à chaque fois

Code :

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

a0 = [500  0  800]

, est ce qu'il y'a un moyen pour que mon code soit valable pour toutes les courbes ?