matlab频域分析

MATLAB频域分析实验1.给定系统8.02.02zzzH(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)求出并绘出系统幅频响应与相频响应；由题知，)8.0)(8.0(2.0)(jzjzzzH，用ejw替换z后得到)8.0)(8.0(2.0)(jejeeeHjwjwjwjw；利用matlab做出幅频响应和相频响应如下：代码如下：clearall;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);gridon;xlabel('\omega/\pi');ylabel('幅频响应');subplot(222);plot(w,Hphase);gridon;xlabel('\omega/\pi');ylabel('相频响应/°');(c)求出并绘出系统的单位样值响应；由(b)中结果，利用matlab中的iztrans命令求得)()21sin(8.08.041-)(nunnhn,在matlab中做出单位样值响应如下：代码如下：clearall;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),求出并绘出系统的单位阶跃响应.由题知，)()()(nhnxng，由卷积定理得G(z)=X(z)*H(z),则)8.0)(1(2.08.02.0*1)(*)()(222zzzzzzzzHzXzG，利用matlab进行逆变换后得)(91)()21cos(8.091)()21sin(0.8452)(nununnunngnn，利用matlab做出图形如下：代码如下：clearall;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);gridon;xlabel('n');ylabel('阶跃响应');2.某离散系统用以下差分方程表示1224.012.0nxnxnynyny.(a)画出系统极零图；将原式两边进行z变换得)1()()())2()1()((24.0))1()((2.0)(1122xzXzzXzyzzyzYzzyzYzzY，由于初始状态为零状态，将上式整理后得)(2)(24.0)(2.0)(1-2zXzzYzzYzzY，则可得系统函数24.02.02)()()(2zzzzXzYzH，极点为P1=0.4，P2=-0.6，零点为Z1=0。由matlab作出极零图如下：代码如下：b=[0,2];a=[1,0.2,-0.24];zplane(b,a);(b)求出并绘出系统幅频响应与相频响应；用ejw接替z得24.02.02)()()(2jwjwjwjwjwjweeeeXeYeH，利用matlab作出频率响应和相频响应如下：代码如下：clearall;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);gridon;xlabel('\omega/\pi');ylabel('幅频响应');subplot(222);plot(w,Hphase);gridon;xlabel('\omega/\pi');ylabel('相频响应/°');(c)求出并绘出系统的单位样值响应；利用matlab求逆变换得)(4.02)()6.0(2)(nununhnn，用matlab作出单位样值响应如下代码如下：clearall;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)nu(n),求出并绘出系统的零状态响应.零状态响应)()()(nhnxnrzs，由卷积定理得Rzs(z)=X(z)*H(z),求得)24.02.0)(1(224.02.021)(222zzzzzzzzzzRzs，利用matlab求逆变换得零状态响应为)()52(74)()53(3)()1(725)(nunununrnnnzs，再由matlab作出其图形如下：代码如下：n=50;stem(n,(4*(2/5)^n)/7-(25*(-1)^n)/7+3*(-3/5)^n,'.');xlabel('n');ylabel('零状态响应rzs(n)');