内容发布更新时间 : 2024/12/26 9:18:01星期一 下面是文章的全部内容请认真阅读。

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(‘脉冲响应’)