matlab频域分析

MATLAB频域分析实验1. 给定系统 (a) 画出系统极零图； 由题知，该系统零点为Z1=0，极点为P1=-(8/10)j和P2=(8/10)j。由matlab作出极零图如下：代码如下：b=0,-0.2;a=1,0,0.8;zplane(b,a);(b) 求出并绘出系统幅频响应与相频响应； 由题知，用ejw替换z后得到；利用matlab做出幅频响应和相频响应如下：代码如下：clear all;b=0,-0.2;a=1,0,0.8;H,w=freqz(b,a,512,1);Hr=abs(H);Hphase=angle(H)*180/3.1415;Hphase=unwrap(Hphase);subplot(221);plot(2*w,Hr);grid on;xlabel( omega/pi);ylabel(幅频响应 );subplot(222);plot(w,Hphase);grid on;xlabel( omega/pi);ylabel(相频响应/);(c) 求出并绘出系统的单位样值响应；由(b) 中结果，利用matlab中的iztrans命令求得,在matlab中做出单位样值响应如下：代码如下：clear all;b=0,-0.2;a=1,0,0.8;N=50;h,t=impz(b,a,N);subplot(223);stem(t,h,.);grid;xlabel( n);ylabel(单位样值响应 h(n) );(d) 令x(n) = u(n), 求出并绘出系统的单位阶跃响应.由题知，由卷积定理得G(z)=X(z)*H(z),则，利用matlab进行逆变换后得，利用matlab做出图形如下：代码如下：clear all; b=0,-0.2;a=1,0,0.8;x=ones(100);t=1:100;y=filter(b,a,x);subplot(224);plot(t,y);grid on;xlabel( n);ylabel(阶跃响应 );2. 某离散系统用以下差分方程表示.(a) 画出系统极零图； 将原式两边进行z变换得，由于初始状态为零状态，将上式整理后得，则可得系统函数，极点为P1=0.4，P2=-0.6，零点为Z1=0。由matlab作出极零图如下：代码如下：b=0,2;a=1,0.2,-0.24;zplane(b,a);(b) 求出并绘出系统幅频响应与相频响应； 用ejw接替z得，利用matlab作出频率响应和相频响应如下：代码如下：clear all;b=0,2;a=1,0.2,-0.24;H,w=freqz(b,a,512,1);Hr=abs(H);Hphase=angle(H)*180/3.1415;Hphase=unwrap(Hphase);subplot(221);plot(2*w,Hr);grid on;xlabel( omega/pi);ylabel(幅频响应 );subplot(222);plot(w,Hphase);grid on;xlabel( omega/pi);ylabel(相频响应/);(c) 求出并绘出系统的单位样值响应；利用matlab求逆变换得，用matlab作出单位样值响应如下代码如下：clear all;b=0,2;a=1,0.2,-0.24;N=50;h,t=impz(b,a,N);subplot(223);stem(t,h,.);grid;xlabel( n);ylabel(单位样值响应 h(n) );(d) 已知输入x(n) = (-1)n u(n), 求出并绘出系统的零状态响应.零状态响应，由卷积定理得Rzs(z)=X(z)*H(z),求得，利用matlab求逆变换得零状态响应为，再由matlab作出其图形如下：代码如下：n=5