Discussion :

[math09]

bonjour,
je travaille avec matlab ma question est:
j'ai une variable x=-1:1 dans une matrice A d'ordre 9x9 A=[1 x 5;2*x 3 6;1 3 je veux tracer les valeurs propres de A en fct de x
le code utilise est

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
23
24
25
26
27
28
29
30
31
32
33
34
x=-1:1;
 
for k = 1:numel(x)
 
 
b1 = 14*x(k) + 4*(1-abs(x(k)))
b2 = 18*x(k)-20*(1-abs(x(k)))
b12 =(8.3666)*x(k)-6*(8.3666)*(1-abs(x(k)))
b3 = 14*x(k) +4*(1-abs(x(k)))
b23 = (8.3666)*x(k)-6*(8.3666)*(1-abs(x(k)))
b4 = -21*x(k)-17*(1-abs(x(k)))
b45 = 5*(2.6458)*x+3*(2.6458)*(1-abs(x(k)))
b5 = 9*x(k) + 87*(1-abs(x(k)))
b6 = -21*x(k)-17*(1-abs(x(k)))
b67 = 5*(2.6458)*x+3*(2.6458)*(1-abs(x(k)))
b7 = 9*x(k) + 87*(1-abs(x(k)))
b8 = -11*x(k)+22*(1-abs(x(k)))
b98 = 15*x(k)+42*(1-abs(x(k)))
b9 = -11*x(k)+22*(1-abs(x(k)));
A=[b1 b12 0 0 0 0 0 0 0;
b12 b2 0 0 0 0 0 0 0;
0 b23 b3 0 0 0 0 0 0;
0 0 0 b4 b45 0 0 0 0;
0 0 0 b45 b5 0 0 0 0;
0 0 0 0 b6 b67 0 0 0;
0 0 0 0 b67 b7 0 0 0;
0 0 0 0 0 0 0 b8 b98;
0 0 0 0 0 0 0 b98 b9];
 
d(:,k)=eig(A)
 
end
 
plot(x,d)

mais ca marche pas
ou est le probleme
merci

[christophe_halgand]

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
23
24
25
26
27
28
29
30
31
32
33
34
x=-1:1;

for k = 1:numel(x)


b1 = 14*x(k) + 4*(1-abs(x(k)));
b2 = 18*x(k)-20*(1-abs(x(k)));
b12 =(8.3666)*x(k)-6*(8.3666)*(1-abs(x(k)));
b3 = 14*x(k) +4*(1-abs(x(k)));
b23 = (8.3666)*x(k)-6*(8.3666)*(1-abs(x(k)));
b4 = -21*x(k)-17*(1-abs(x(k)));
b45 = 5*(2.6458)*x(k)+3*(2.6458)*(1-abs(x(k)));
b5 = 9*x(k) + 87*(1-abs(x(k)));
b6 = -21*x(k)-17*(1-abs(x(k)));
b67 = 5*(2.6458)*x(k)+3*(2.6458)*(1-abs(x(k)));
b7 = 9*x(k) + 87*(1-abs(x(k)));
b8 = -11*x(k)+22*(1-abs(x(k)));
b98 = 15*x(k)+42*(1-abs(x(k)));
b9 = -11*x(k)+22*(1-abs(x(k)));
A=[b1 b12 0 0 0 0 0 0 0;...
b12 b2 0 0 0 0 0 0 0;...
0 b23 b3 0 0 0 0 0 0;...
0 0 0 b4 b45 0 0 0 0;...
0 0 0 b45 b5 0 0 0 0;...
0 0 0 0 b6 b67 0 0 0;...
0 0 0 0 b67 b7 0 0 0;...
0 0 0 0 0 0 0 b8 b98;...
0 0 0 0 0 0 0 b98 b9];

d(:,k)=eig(A)

end

plot(x,d)