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));