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
| clear all; clc;
% Question 1
n=[0 400];
d=[0 1 12 400];
f=tf(n,d) % FT = n / d
K=dcgain(f) %Gain statique
Zeros=zero(f) %Zeros de la FT 'f'
Poles=pole(f) %Pole de la FT 'f'
[r,p,k]=residue(n,d) % r= coeff a,b,c de la decomp' en element simple
% p= Poles
% k= facteur de gain
F=zpk(f) %forme factorisée ou Evans ou 'zpk'
[Wn,Z] = damp(f) % Wn pulsation propre et Z coefficient d'amortissement
W=Wn(1);
z=Z(1);
disp('Forme Generale de la reponse indicielle')
t=0:0.1:10;
q=(1-((exp((-z*W)*t))/sqrt(1-z^2))*cos(W*t*sqrt(1-z^2)-asin(z)));
plot(t,q);
figure(1); %Reponse a un echelon
step(f);grid; |
Partager