matlab z变换

1:  clear all;

2:  close all;

3:  clc;

4:  %  syms n k z;

5:  % f=1/n;         %定义离散信号

6:  %  F=ztrans(f)              %z变换

7:  %  pretty(F);

8:  %  fk=iztrans(F,k)

9:  %  pretty(fk)

10:  %   

11:  % b=[0 0 10 0];%分子的系数数组 

12:  % a=[1 -5 8 -4]; %分母的系数数组 

13:  % zplane(b,a)% 使用zplane函数绘制如下系统的零极点分布图 

14:  % [z,p,k]=tf2zp(b,a)%使用tf2zp求出系统函数的零、极点和增益运算结果

15:  % [b,a] = zp2tf(z,p,k) %根据求出的零、极点和增益，然后自学使用zp2tf还原出 分子和分母的系数 

16:  % [sos,g]=tf2sos(b,a)%使用tf2sos将系统函数分解成一系列二阶子系统的级联形式

17:  % [b,a]=sos2tf(sos,g)%根据求出的一系列二阶子系统，使用sos2tf还原出)(zH分子和分母的系数

18:  

19:  b=[0 1 1 1]; %分子的系数数组 

20:  a=[1 2 2 2]; %分母的系数数组

21:  n=(0:500)*pi/500; %在pi范围内取501个采样点 

22:  [h,w]=freqz(b,a,n);%求系统的频率响应 

23:  subplot(2,1,1),plot(n/pi,abs(h));

24:  grid %作系统的幅度频响图 

25:  axis([0,1,1.1*min(abs(h)),1.1*max(abs(h))]); 

26:  ylabel('幅度'); 

27:  subplot(2,1,2),plot(n/pi,angle(h));

28:  grid %作系统的相位频响图 

29:  axis([0,1,1.1*min(angle(h)),1.1*max(angle(h))]);

30:  ylabel('相位');

31:  xlabel('以pi为单位的频率'); 

32:  

33:  figure;plot(w,abs(h));