Graphe à partir d'une matrice variable

bonjour à tous et à toutes

j'ai un système des matrices que je doit retire une valeur pour le présenter sur une graphe en fonction de x

Code:

---

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22

clear;clc;clf; Fnl=0 ts=pi nu=0.5   for x=0:10   u=0.8 a=u^2 b=x.^2 A0=(a-b)./(a-b-(a*b)*(1+nu)+b.^2) a0=b./(a-b-(a*b)*(1+nu)+b.^2) d=sqrt((1+(1+nu)*(a))^2-4*a) w1=sqrt(((1+(1+nu)*a)-d)./(2*b)) w2=sqrt(((1+(1+nu)*a)+d)./(2*b)) T1=((b.*w1.^2)-1)./(b.*w1.^2) T2=((b.*w2.^2)-1)./(b.*w2.^2) ts=pi H0=[0,1,0,1;w1,0,w2,0;0,T1,0,T2;T1*w1,0,T2*w2,0] H1=inv([a0+l,0,-l,0;0,-a0-l,0,-l;A0,0,-1,0;0,-A0,0,-1]) end H= H0*H1

---

comment sortir de cette boucle et obtenir une courbe de l’élément H(3,2) en fonction de x

merci

Bonjour,
Ci dessous un exemple de solution :

Code:

---

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

tabX = 0:10; H32 = zeros(1,11); for idX=1:numel(tabX) x = tabX(idX); %[...] H0=[0,1,0,1;w1,0,w2,0;0,T1,0,T2;T1*w1,0,T2*w2,0]; H1=inv([a0+l,0,-l,0;0,-a0-l,0,-l;A0,0,-1,0;0,-A0,0,-1]); H= H0*H1; H32(x+1) = H(3,2); end figure; plot(tabX,H32);

---