控制系统数字仿真与CAD第4章仿真题答案.docVIP

《控制系统数字仿真与CAD》 第四章上机实验报告 院系：电信工程学院 班级：11级控制工程基地1班 姓名： 张亚东 学号： 时间： 2014 年 5 月 27日 例 4-1 求阶跃响应的数值解。 解： 先建立一个名为sp4-1.m的子程序如下： b=b/a(1);a=a/a(1);A=a(2:n+1); A=[rot90(rot90(eye(n-1,n)));-fliplr(A)]; B=[zeros(1,n-1),1]; m1=length(b); C=[fliplr(b),zeros(1,n-m1)]; Ab=A-B*C*V; X=X0;y=0;t=T0; N=round(Tf-T0)/h; for i=1:N K1=Ab*X+B*R; K2=Ab*(X+h*K1/2)+B*R; K3=Ab*(X+h*K2/2)+B*R; K4=Ab*(X+h*K3)+B*R; X=X+h*(K1+2*K2+2*K3+K4)/6; y=[y,C*X]; t=[t,t(i)+h]; end [t,y] plot(t,y) 然后对上面的子程序进行调用，改变开环放大系数k，得到k取不同值时输出响应的曲线如下： k=2; a=conv([1 0 0],conv([0.25 1],[0.25 1])) a