实验4Z变换及系统频域特性的MATLAB实现.docxVIP

小实验4 Z变换和系统频域特性的MATLAB实现实验目的学习通过Z变换来分析离散系统的频率响应，并用MATLAB实现。加深对系统的零、极点分布概念的理解。2.实例分析2.1通过Z变换分析求解系统的冲激响应已知用线性常系数差分方程：表示的线性时不变系统其系统函数为：（3-1）上式为两个关于的多项式之比，即为有理分式。同时，式可以表示成部分分式的形式：则可以通过所熟悉的常见序列的Z变换形式求得的Z反变换，从而求得系统的冲激响应函数。MATLAB提供了一个内部函数residuez( )，来计算有理多项式的留数和直接项，residuez( )函数有几种调用方式：（1）[R,p,C]=residuez(b,a)，在已知以分子行向量b和分母行向量a下，得到列向量R含有留数，列向量p是极点位置，行向量C包含直接项；（2）[b,a] = residuez(R,p,C)，将部分分式展开式转换到分子行向量b和分母行向量a.MATLAB还提供了一个内部函数impz(b,a,N)，在已知分子行向量b和分母行向量a下，计算N点的单位冲激响应。例3.1求系统：的单位冲激响应。解：由上述基本原理和MATLAB提供的函数，将下列指令编辑到“exe3impz.m”文件中，可以得到所求系统的冲激响应。% exe3impz.m h(n)求解b=[0.001836, 0.007344, 0.011016, 0.007374, 0.001836];a=[1, -3.0544, 3.8291, -2.2925, 0.55075];[h,n] = impz(b,a,40);stem(n,h,.);ylabel(h[n]);grid on运行“exe3impz.m”文件将产生如图3-1所示的序列。图4-1 系统冲激响应例3.2求的Z反变换。解：由上述基本原理和MATLAB提供的函数，将下列指令编辑到“exe3zinver.m”文件中，来求Z反变换。% exe3zinver.mz反变换求解b=[0,1,0];a=[3,-4,1];[R,p,C]=residuez(b,a)[b,a]=residuez(R,p,C)运行结果如下：R = 0.5000 -0.5000p = 1.0000 0.3333C = 0b = -0.0000 0.3333 0a = 1.0000 -1.3333 0.3333因此得到因式分解后的，所以Z反变换的结果为。2.2通过Z变换分析离散系统的频率响应和零、极点分布将式（3-1）表示的的分子、分母进行因式分解，可采用根的形式表示多项式，即其中，为分子多项式的根，称为系统函数的零点，为分母多项式的根，称为系统函数的极点，A为比例。这样系统函数就表示成了零极点的形式，零极点在Z平面的位置刻画了系统很重要的特性，可以通过系统函数零极点位置估算出系统函数的频率响应，进而判断系统的滤波特性，这是一种非常实用的方法，也称为频率响应的几何确定法。根据频率响应的定义，将代入上式，系统频率响应为：MATLAB提供了一个内部函数zplane(b,a)，在已知系统函数的分子行向量b和分母行向量a下，画出极点和零点。MATLAB还提供了一个内部函数freqz( )来求系统的频率响应，freqz( )函数有几种调用方式：（1）[H,w]=freqz(b,a,N)，在已知分子行向量b和分母行向量a情况下得到N点的频率向量w和N点的系统复频率响应向量H，频率响应在单位圆的上半圆的N个等分点上求值；（2）[H,w]=freqz(b,a,N,’whole’)，在已知分子行向量b和分母行向量a 情况下得到N点的频率向量w和N点的系统复频率响应向量H，频率响应用环绕整个单位圆的N个等分点求值；（3）H=freqz(b,a,w)，在已知分子行向量b和分母行向量a情况下得到在频率向量w上的频率响应向量H。例3.3画出因果系统的零极点图和频率响应曲线，并求出脉冲响应。解：由上述的基本原理和MATLAB提供的函数，将下列指令编辑到“exe3freqz.m”文件中，即可求得系统的零极点图和频率响应。% exe3freqz.m系统求解b=[1,0];a=[1,-0.9];subplot(2,2,1);zplane(b,a);[H,w]=freqz(b,a,100);magH=abs(H);phaH=angle(H);subplot(2,2,3);plot(w/pi,magH);grid onxlabel(Frequency in Pi uints);ylabel(Magnitude);title(Magnitude response);subplot(2,2,4);plot(w/pi,phaH/pi);xlabel(Freque