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