code matlab

[Leond95]

Bonsoir à tous,

J'ai écris un programme matlab qui permet de tracer le diagramme de bode en gain et en phase, sans utilisé la commande bode mais quand j’exécute le programme, la phase est constante (0°), est ce que quelqu'un peut corrigé mon programme.

Cdlt

w=logspace(-3,3,120);
h=@(s) 1./(s.^2+2.*s+1);
z = 20*log(abs(h(w*i)));

y = angle(h(w*i))*180/pi

subplot(2,1,1)
semilogx(w,z)
grid on
subplot(2,1,2)
semilogx(w,y)
grid on

-----

[Antoane]

Bonjour,

Probablement la variable i ne renvoie pas sqrt(-1) mais a gardé une valeur réelle prise au cours d'un traitement précédent.
Cela peut se solutionner en ajoutant un "clear i" ou un "i=sqrt(-1);" en début de programme ou, mieux, en remplacant "i" par "1i" ou "ij" dans le programme.

[Leond95]

merci Antoane ça a marché