matlab系统响应例子

)()(')(2)('2)()1(tftftytyty1、求零输入响应,绘仿真波形图。eq='D2y+4*Dy+3*y=0';%齐次方程求零输入响应%cond='y(0)=2,Dy(0)=1';yzi=dsolve(eq,cond);ezplot(yzi,[-1,8]);gridontitle('零输入响应');eq='D2y+4*Dy+4*y=0';%齐次方程求零输入响应%cond='y(0)=2,Dy(0)=1';yzi=dsolve(eq,cond);ezplot(yzi,[-1,8]);gridontitle('零输入响应');2、求冲击响应和阶跃响应,绘图t=0:0.001:8;sys=tf([1,1],[1,2,2]);ht=impulse(sys,t);%冲击响应gt=step(sys,t);%阶跃响应subplot(2,1,1)plot(t,ht),gridonxlabel('Time(sec)'),ylabel('h(t)');title('冲击响应');subplot(2,1,2)plot(t,gt),gridonxlabel('Time(sec)'),ylabel('g(t)');title('阶跃响应');)(')(2)('2)()2(tftytytyt=0:0.001:8;sys=tf([1],[1,2,2]);ht=impulse(sys,t);%冲击响应gt=step(sys,t);%阶跃响应subplot(2,1,1)plot(t,ht),gridonxlabel('Time(sec)'),ylabel('h(t)');title('冲击响应');subplot(2,1,2)plot(t,gt),gridonxlabel('Time(sec)'),ylabel('g(t)');title('阶跃响应');【例5】已知LTI系统的微分方程为，求频率响应w=-3*pi:0.01:3*pi;b=[137];a=[11085];H=freqs(b,a,w);subplot(2,1,1);plot(w,abs(H)),gridonxlabel('\omega(rad/s)'),ylabel('|H(\omega)|');title('H(w)的幅频特性')subplot(2,1,2);plot(w,angle(H)),gridonxlabel('\omega(rad/s)'),ylabel('\phi(\omega)');title('H(w)的相频特性'))(6)(6)('5)(tftytyty)()(tuetft求零输入响应ts=0;te=5;dt=0.01;sys=tf([6],[1,5,6]);t=ts:dt:te;f=6*exp(-t).*u(t);y=lsim(sys,f,t);plot(t,y),gridonxlabel('Time(sec)'),ylabel('y(t)');title('零状态响应');