搜索
中国科学技术大学信息学院电子科学与技术系MadebySheldong-1-信号与系统实验报告6学号：姓名：信息科学技术学院电子科学与技术系中国科学技术大学信息学院电子科学与技术系MadebySheldong-2-一、实验目的1.熟悉离散信号Z变换的原理及性质2.熟悉常见信号的Z变换3.了解正/反Z变换的MATLAB实现方法4.了解离散信号的Z变换与其对应的理想抽样信号的傅氏变换和拉氏变换之间的关系5.了解利用MATLAB实现离散系统的频率特性分析的方法二、实验原理1.正/反Z变换Z变换分析法是分析离散时间信号与系统的重要手段。如果以时间间隔sT对连续时间信号f(t)进行理想抽样，那么，所得的理想抽样信号()ft为：()()*()()*()TsskftfttfttkT理想抽样信号()ft的双边拉普拉斯变换F(s)为：()()*()()sksTstsskkFsfttkTedtfkTe若令()()sfkTfk，sTsze，那么()ft的双边拉普拉斯变换F(s)为：()()()sTskzekFsfkzFz则离散信号f(k)的Z变换定义为：()()kkFzfkz从上面关于Z变换的推导过程中可知，离散信号f(k)的Z变换F(z)与其对应的理想抽样信号()ft的拉氏变换F(s)之间存在以下关系：()()sTszeFsFz同理，可以推出离散信号f(k)的Z变换F(z)和它对应的理想抽样信号()ft的傅里叶变换之间的关系为()()jTszeFjFz如果已知信号的Z变换F(z)，要求出所对应的原离散序列f(k)，就需要进行反Z变换，反Z变换的定义为：11()()2kfkFzzdzj中国科学技术大学信息学院电子科学与技术系MadebySheldong-3-其中，C为包围1()kFzz的所有极点的闭合积分路线。MATLAB中专门对信号进行正反Z变换的函数ztrans()和itrans()F=ztrans(f)对f(n)进行Z变换，其结果为F(z)F=ztrans(f,v)对f(n)进行Z变换，其结果为F(v)F=ztrans(f,u,v)对f(u)进行Z变换，其结果为F(v)f=itrans(F)对F(z)进行Z反变换，其结果为f(n)f=itrans(F,u)对F(z)进行Z反变换，其结果为f(u)f=itrans(F,v,u)对F(v)进行Z反变换，其结果为f(u)注意事项在调用函数ztran()及iztran()之前，要用syms命令对所有需要用到的变量（如t,u,v,w）等进行说明，即要将这些变量说明成符号变量。例①用MATLAB求出离散序列()(0.5)()kfkk的Z变换。MATLAB程序如下：symskzf=0.5^k;%定义离散信号Fz=ztrans(f)%对离散信号进行Z变换运行结果如下：Fz=2*z/(2*z-1)例②已知一离散信号的Z变换式为2()21zFzz，求出它所对应的离散信号f(k)。MATLAB程序如下：symskzFz=2*z/(2*z-1);%定义Z变换表达式fk=iztrans(Fz,k)%求反Z变换运行结果如下：fk=(1/2)^k中国科学技术大学信息学院电子科学与技术系MadebySheldong-4-2.离散系统的频率特性同连续系统的系统函数H(s)类似，离散系统的系统函数H(z)也反映了系统本身固有的特性。对于离散系统来说，如果把其系统函数H(z)中的复变量z换成jTe，那么所得的函数()jTHe就是此离散系统的频率响应特性，即离散时间系统的频率响应为：()()()()jTjTjTjzeHeHeeHz其中，()jTHe称为离散系统的幅频特性，()称为系统的相频特性。同连续系统一样，离散时间系统的幅频特性也是频率的偶函数，相频特性也是频率的齐函数。由于jTe是频率的周期函数，所以离散系统的频率响应特性也是频率的周期函数，其周期为2T，或者频率周期为2TT。实际上，这就是抽样系统的抽样频率，而其中的T则是系统的抽样周期。频率响应呈现周期性是离散系统特性区别于连续系统特性的重要特点。因此，只要分析()jTHe在2T范围内的情况，便可分析出系统的整个频率特性。鉴于离散系统频率响应特性的特点，为了表示方便起见，我们通常将其中的T用一个变量来代替，即令jze代入系统函数H(z)中，用()jHe函数来表示离散系统的频率响应特性。相应地，用()jHe表示幅频特性，而相频特性仍用()来表示。应该特别注意的是，虽然这里的变量仍然称为频率变量，但是它已经不是原来意义上的角频率概念，而实际上是表示角度的概念。我们称之为数字角频率。它与原来角频率的关系为：T。也就是说，根据离散系统的系统函数H(z)，令其中的jze，并且代入0～2范围内不同的频率值（实际上是角度值），就可以逐个计算出不同频率时的响应，求出离散系统的频率响应特性。再利用离散系统频率特性的周期性特点（周期为2），求出系统的整个频率特性。离散系统的幅频特性曲线和相频特性曲线能够直观地反映出系统对不同频率的输入序列的处理情况。在函数()jHe随的变换关系中，在=0附近，反映了系统对输入信号低频部分的处理情况，而在附近，则反映了系统对输入信号高频部分的处理情况。一般来说，分析离散系统频率响应特性就要绘制频率响应曲线，而这是相当麻烦的。虽然可以通过几何矢量法来定性画出频率响应特性曲线，但一般来说这也是很麻烦的。值得庆幸的是，MATLAB为我们提供了专门用于求解离散系统频率响应的函数freqz()，其调用格式如下：中国科学技术大学信息学院电子科学与技术系MadebySheldong-5-[H,w]=freqz(B,A,N)其中，B和A分别是表示待分析的离散系统的系统函数的分子，分母多项式的向量，N为正整数，返回向量H则包含了离散系统频率响应函数()jHe在0~范围内的N个频率等分点的值。向量则包含0~范围内的N个频率等分点。在默认情况下N=512。[H,w]=freqqz(B,A,N,'whole')其中，B，A和N的意义同上，而返回向量H包含了频率响应函数()jHe在0~2范围内N个频率等分点的值。由于调用freqz()函数只能求出离散系统频率响应的数值，不能直接绘制曲线图，因此，我们可以先用freqz()函数求出系统频率响应的值，然后再利用MATLAB的abs()和angle()函数以及plot()命令，即可绘制出系统在0~或0~2范围内的幅频特性和相频特性曲线。例①若离散系统的系统函数为0.5()zHzz，请用MATLAB计算0~频率范围内10个等分点的频率响应()jHe的样值。MATLAB程序如下：A=[10];%分母多项式系数向量B=[1-0.5];%分子多项式系数向量[H,w]=freqz(B,A,10)%求出对应0~范围内10个频率点的频率响应样值例②用MATLAB计算前面离散系统在0~2频率范围内200个频率等分点的频率响应值，并绘出相应的幅频特性和相频特性曲线。MATLAB程序如下：A=[10];B=[1-0.5];%[H,w]=freqz(B,A,200);[H,w]=freqz(B,A,200,'whole');%求出对应0~2范围内200个频率点的频率响%应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)%画出幅频特性曲线subplot(2,1,2);plot(w,HX)%画出相频特性曲线运行结果如下：中国科学技术大学信息学院电子科学与技术系MadebySheldong-6-运行结果分析：从该系统的幅频特性曲线可以看出，该系统呈高通特性，是一阶高通滤波器。三、实验内容1．求出下列离散序列的Z变换【实验步骤】：1.通过符号函数构造fk2.利用库函数ztrans进行Z变换④4()(1)()(5)fkkkkk【理论值计算】：Fz=2z^-2+6z^-3+12z^-4【程序如下】：clear;symskz;Fk=k*(k-1)*(heaviside(k)-heaviside(k-5));%定义FkFz=ztrans(Fk,k,z)%使用库函数计算Z变换【matlab结果】：(-25*z^3+3*z^2+43*z-25)/(2*z^5*(z-1)^3)-z/(z-1)^2+(5*z^2+2*z-5)/(2*z^5*(z-1)^2)+(z*(z+1))/(z-1)^3【差异分析】：matlab2018中heasvide（0）=0.5，这导致了Z变换的时候离散点的数值出现了偏差，这导致了仿真结果的偏差2．已知下列单边离散序列的z变换表达式，求其对应的原离散序列。【实验步骤】：1．通过符号函数构造fz2.利用库函数iztrans进行反Z变换中国科学技术大学信息学院电子科学与技术系MadebySheldong-7-④24(1)()(1)(2)(3)zzzFzzzz【理论值计算】：Fk=((7*2^k)/15-(-1)^k/6+(7*(-3)^k)/10)*Heaviside(k)【程序如下】：clear;symskz;Fz=z*(z*z+z+1)/((z+1)*(z-2)*(z+3));%定义FzFk=iztrans(Fz,z,k)%用库函数求Fk【matlab结果】：Fk=(7*2^k)/15-(-1)^k/6+(7*(-3)^k)/10【差异分析】：在matlab的iztrans中由于是有限数值计算，所以没有计算k0的数值，这导致在总结Fk的时候没有加heaviside函数，这导致了最后结果的偏差。3.已知离散系统的系统函数H(z)如下，请绘出系统的幅频和相频特性曲线，并说明系统的作用【实验步骤】：1.构造矢量A、B描述差分方程2.利用库函数freqz进行频域分析①122344()()()zHzzz【零极点分析】：零点是z=-1，极点是z=-1/2和-2/3【程序如下】：clear;A=[11/2+2/31/3];B=[044];%[H,w]=freqz(B,A,200);[H,w]=freqz(B,A,200,'whole');%求出对应范围内200个频率点的频率响应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)%画出幅频特性曲线xlabel(w);ylabel(abs(H));subplot(2,1,2);plot(w,HX)%画出相频特性曲线中国科学技术大学信息学院电子科学与技术系MadebySheldong-8-xlabel(w);ylabel(angle(H));【画图如下】：【系统分析】：对相频和幅频分析可知，该系统对中低频有比较好的响应，而在对高频有很高的抑制作用。②221()0.81zHzz【零极点分析】：零点是z=+-1，极点是z=+-0.9i【程序如下】：clear;A=[100.81];B=[10-1];%[H,w]=freqz(B,A,200);[H,w]=freqz(B,A,200,'whole');%求出对应范围内200个频率点的频率响应样值HF=abs(H);%求出幅频特性值HX=angle(H);%求出相频特性值subplot(2,1,1);plot(w,HF)%画出幅频特性曲线xlabel(w);ylabel(abs(H));subplot(2,1,2);plot(w,HX)%画出相频特性曲线xlabel(w);ylabel(angle(H));中国科学技术大学信息学院电子科学与技术系MadebySheldong-9-【画图如下】：【系统分析】：对相频和幅频分析可知，该系统对高频有比较好的响应，而在对低频有很高的抑制作用。4.已知描述离散系