基于Matlab的零输入和零状态响应

)?Р及初始状态。输入函数。Р先求出系统的冲激响应,写出其特征方程Р Р求出其特征根为p和p,及相应的留数r,r;则冲激响应为Р Р输入y(t)可用输入u(t)与冲激响应h(t)的卷积求得。РMATLAB程序:Рa=[1,3,6];Рn=length(a)-1;РY0=[0,10];Рb=[1,6,8];Рdt=0.001;te=9;Рt=0:dt:te;Рu=sin(2*t)+cos(5*t);Рp=roots(a); %求出齐次解РV=rot90(vander(p)); %生成2×2矩阵并旋转90度Рc=V\Y0'; %求出待定系数r1,r2的值Рy1=zeros(1,length(t)); %建立一个1×length(t)的矩阵Рfor k=1:n y1=y1+c(k)*exp(p(k)*t);end %求得冲击响应Рdt=te/(length(t)-1); [r,p,k]=residue(b,a);Рh=r(1)*exp(p(1)*t)+r(2)*exp(p(1)*t);Рy2=conv(u,h)*dt;Рsubplot(2,1,1),plot(t,y1),gridРxlabel('t'); ylabel('y1');title('零输入响应');Рsubplot(2,1,2),plot(t,y2(1:length(t)));gridРxlabel('t');ylabel('y2');title('零状态响应');Р参考文献:信号与线性系统分析(第四版) 主编吴大正Р五、总结与体会Р通过这次信号的课程设计,我和组员加强了对matlab软件的学习,更加清楚了信号系统中函数卷积和系统零输入、零响应状态的公式及其算法,同时在编程过程了发现了许多问题,并通过不断的修正,从而使我们对matlab软件的编程思想有了更进一步的认识与了解,总之通过这次信号系统的课程设计,收获颇丰。