plot avec une boucle for

bonsoir,
Je vais essayer d'expliquer mon problème clairement.
Je dois résoudre une EDO pour caractériser le comportement mécanique d'un métal subissant des chargement cyclique, on obtient un cycle d'hytérésis.
Mon système est définit comme ceci:

Code:

---

1 2 3 4 5 6 7 8 9 10 11

function dy=f_iso(t,y) global Q b K n E sigma0 epoint if abs(y(1))-Q*(1-exp(-b*y(2)))-sigma0<0 dy(2)=0 else dy(2)=((abs(y(1))-Q*(1-exp(-b*y(2)))-sigma0)/K)^n end dy(3)=sign(y(1))*dy(2); dy(1)=E*(epoint-dy(3)); dy=dy'; end

---

Puis ma fonction principale est:

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14

global Q b K n E sigma0 epoint b=5; Q=-150 ; E=140000; sigma0=200; K=800; n=6; epoint=0.001; % 1ere branche [Tn,N]=ode45('f_iso',[0 10],[0 0 0]); sigma1=N(:,1); ep1=N(:,3); e1=ep1+sigma1./E; plot(e1,sigma1)

---

Ce qui me donne la première branche.
Ensuite je veux tracer 9 autres branches sachant que à partir de la deuxième ça devient cyclique.
Pour cela je fais une boucle for :

Code:

---

1 2 3 4 5 6 7 8 9

L=length(N) for w = 0:10 epoint=-epoint [Tn,N]=ode45('f_iso',[10*(2*w+1) (30+20*w)],[N(L,:)]); sigma2=N(:,1); ep2=N(:,3); e2=ep2+sigma2./E; L=length(N) end

---

A chaque itération Ti=Ti+20 et Tf=Tf+20 et comme condition initiale je prends la valeur du dernier vecteur précédent .
Mais je ne sais pas comment ployer pour faire en sorte d'avoir toutes les courbes sigma=f(e) sur un meme grave avec la boucle for...
De plus avec cette boucle, les valeurs vont être écraser à chaque itération donc ça va pas marcher...

merci

PLus de précisions.
La boucle marche très bien fait , mais lorsque je plot, je n'ai que le tracé de la dernière itération cependant moi je voudrais tracer toutes les courbes mais je ne sais pas comment faire...

Salut,
si j'ai bien compris tu cherches à sauvegarder e2 et sigma2 à chaque itération de ta boucle :

Code:

---

1 2 3 4 5 6 7 8 9

sauvegardeE2 = zeros(L,11); sauvegardeSigma2 = zeros(L,11) for w = 0:10 %[...] L=length(N); sauvegardeE2(1:L,w+1) = e2; sauvegardeSigma2(1:L,w+1) = sigma2; end plot(sauvegardeE2,sauvegardeSigma2);

---

Bonjour,

Pourquoi ne pas tout simplement utiliser un hold on après le premier plot et rajouter un plot de e2 et sigma 2 dans la boucle ?

Code:

---

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 26 27 28 29

global Q b K n E sigma0 epoint b=5; Q=-150 ; E=140000; sigma0=200; K=800; n=6; epoint=0.001; % 1ere branche [Tn,N]=ode45('f_iso',[0 10],[0 0 0]); sigma1=N(:,1); ep1=N(:,3); e1=ep1+sigma1./E; plot(e1,sigma1) hold on L=length(N); for w = 0:10 epoint=-epoint; [Tn,N]=ode45('f_iso',[10*(2*w+1) (30+20*w)],[N(L,:)]); sigma2=N(:,1); ep2=N(:,3); e2=ep2+sigma2./E; L=length(N); plot(e2,sigma2) end

---

Salut,

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12

L=length(N) for w = 0:10 epoint=-epoint [Tn,N]=ode45('f_iso',[10*(2*w+1) (30+20*w)],[N(L,:)]); sigma2=N(:,1); ep2=N(:,3); e2=ep2+sigma2./E; L=length(N) figure % ouvre une nouvelle figure et sera donc appelé à chaque itération de ta boucle % Ici tu peux donc tracer ce que tu veux dans ta nouvelle figure end

---

cordialement,