大学信号与系统实验MATLAB编程答案

2.已知频率特性函数为:

2(j?)3?(j?)?4 H(j?)?432(j?)?3(j?)?2(j?)?5(j?)?2求其幅频特性和相频特性。(必做) w = linspace(-5,5,200); b = [2 0 1 4]; a = [1 3 2 5 2]; H = freqs(b,a,w); subplot(2,1,1); plot(w,abs(H));

set(gca,'xtick',[-5 -4 -3 -2 -1 0 1 2 3 4 5]); set(gca,'ytick',[0 0.4 0.707 1]); xlabel('\\omega'); ylabel('|H(j\\omega)|'); title('幅频特性曲线'); subplot(2,1,2); plot(w,angle(H));

set(gca,'xtick',[-5 -4 -3 -2 -1 0 1 2 3 4 5]); xlabel('\\omega'); ylabel('\\phi(\\omega)'); title('相频特性曲线');

四、系统响应的仿真

1.设H(s)?s

(s?p1)(s?p2)设①p1=-2,p2=-30; ②p1=-2,p2=3

(1)针对极点参数①,画出系统零、极点分布图,判断该系统稳定性。 num=[1,0];

den=[1,32,60];

[r,poles,k]=residue(num,den); p=roots(den); z=roots(num);

plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on; title('零、极点分布');

针对极点参数②,画出系统零、极点分布图,判断该系统稳定性。 num=[1,0];

den=[1,-1,-6];

[r,poles,k]=residue(num,den); p=roots(den); z=roots(num);

plot(real(p),imag(p),'*');hold on; plot(real(z),imag(z),'o');grid on; title('零、极点分布');

(2)针对极点参数①,绘出系统的脉冲响应曲线,并观察t→∞时,脉冲响应变化趋势。

num=[1,0]; den=[1,32,60]; T=0:0.001:3;

y1=impulse(num,den,T); y2=step(num,den,T); U=sin(T);

y3=lsim(num,den,U,T); plot(T,y1);title('脉冲响应')

针对极点参数②,绘出系统的脉冲响应曲线,并观察t→∞时,脉冲响应变化趋势。

num=[1,0]; den=[1,-1,-6]; T=0:0.001:3;

y1=impulse(num,den,T); y2=step(num,den,T); U=sin(T);

y3=lsim(num,den,U,T); plot(T,y1);title(‘脉冲响应’)

联系客服:779662525#qq.com(#替换为@) 苏ICP备20003344号-4