控制系统数字仿真与CAD（第3版）教学课件作者张晓华第四章2课件.docVIP

4-2设典型闭环结构控制系统如图4-47所示，当阶跃输入幅值 时，用sp4_1.m求取输出的响应。 解：用sp4_1.m求解过程如下： 在MATLAB语言环境下，输入以下命令语句 a=[0.016 0.864 3.27 3.42 1]; b=[30 25]; X0=[0 0 0 0]; %系统状态向量初值为零 V=2; %反馈系数 n=4; T0=0;Tf=10; h=0.01;R=20 ; %仿真步长h=0.01，阶跃输入幅值 sp4_1 %调用sp4_1.m函数 plot(t,y) 运行结果为： 附：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