АЧХ и ФЧХ

function plottanh
 
kzv1 = tf([1 1 -2], [0.12 0.16 3.04  1]) % обьявление передаточной функции
kzv2=tf([1 1 -2], [0.12 0.16 3.04  1],0.01)
 
num=[1 1 -2];
den =[0.12 0.16 3.04  1];
 zplane (num,den),grid 
 
 
 title(' карта нулей и полюсов'),
  figure
 
[r,p,c]=residuez(num, den)% цифровая 
[num, den]=residuez(r,p,c)
[r,p,c]=residue(num, den) %аналогавая
[num, den]=residue(r,p,c)  
 
k1 = c2d(kzv1,0.01,'tustin')%перевод в дискретную форму
k2 = d2c(kzv2,'tustin')%переводит дискретную систему в непрерывную
 
step(kzv1),grid %характнристики аналоговой системы
 figure
impulse(kzv1),grid 
 figure
nyquist(kzv1),grid 
 figure
bode(kzv1),grid 
 figure
 
 
t = 0:0.01:40; u = sin(t);
ssys=tf(kzv1)
lsim(ssys,u,t);grid
 figure
 
 
num2=[ 0.04157  - 0.04116   - 0.04158     0.04116   ];
den2=[1 - 2.984  2.971    - 0.9868 ];
 
 
num3=[     0.04116 - 0.04158 - 0.04116  0.04157   ];
den3=[  - 0.9868  2.971 - 2.984 1];
 
 
 
 b=num3;
a =den3;
 
 N=500;
Fs=200000;
 
 [h,nT]=impz(b,a,N,Fs);%характнристики дискретной системы
 stem(nT,h),grid
 title('impz'),
 
 figure
 
 w=0:pi/10:pi;
H=freqz( num,den,w);
 
 
f=0:0.01:500;
Fs=1000;
H=freqz( num,den,f,Fs);
A=abs(H);
FI=angle(H);
plot(f,A),grid
 
 gtext('АЧХ')
hold off
 
  title('AЧХ'),
 xlabel('частота'),
 ylabel('амплитуда'),
 figure
 
 plot(f,FI),grid 
  gtext('ФЧХ')
 title('ФЧХ'),
 xlabel('частота'),
 ylabel('фаза'),
 
 
 
 
end