Discussion :

[aymenvictoire]

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 :

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
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

[FLB]

Bonjour,
Ci dessous un exemple de solution :

Code :

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

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);