状态方程simulink仿真_控制系统设计与仿真作业与复习资料

仿真的三次作业和三次上机实验报告，百度云链接: http://pan.baidu.com/s/1eRsWZ2i 密码: 9s8v注意不要直接抄袭喔！如果老师发现所有学生的作业都一样的话你们懂的

当时考试的时候总结的，考题基本【或者肯定】就从这里出……

---

采用四阶龙格库塔法求如下二阶系统的阶跃响应。

解：

T=10;%仿真时长b=0.5;w=5; %%%%%%%%%%%采用矩阵方法的四阶龙格库塔法%%%%%dt=0.05;%积分步长N=T/dt;%积分总步数Y1=[0;0];YY=Y1;A=[0 1;-w^2 -2*b*w];%%%%状态观测矩阵%%%%B=[0;1];C=[1 0];U=w^2;for i=1:N;    K1=A*YY+B*U;    K2=A*(YY+dt/2*K1)+B*U;    K3=A*(YY+dt/2*K2)+B*U;    K4=A*(YY+dt*K3)+B*U;    YY=YY+dt/6*(K1+2*K2+2*K3+K4);    Y1=[Y1 YY];end;figureplot([0:dt:T],C*Y1);grid on

求脉冲响应：

%把U定义在for中for i=1:N;if i < 1/dt        U = w^2;    else        U = 0;    end    K1=A*YY+B*U;    K2=A*(YY+dt/2*K1)+B*U;    K3=A*(YY+dt/2*K2)+B*U;    K4=A*(YY+dt*K3)+B*U;    YY=YY+dt/6*(K1+2*K2+2*K3+K4);    Y1=[Y1 YY];end;

ode45解：

function r=second_order(t,x);b=0.5;w=5;r=[x(2);-2*b*w*x(2)-w^2*x(1)+w^2];endt=[0 T];x0=[0 0];[t y]=ode45('second_order',t,x0);figure;hold onplot(t,y(:,1));

---

采用四阶龙格库塔法求高阶系统单位阶跃响应曲线的数值解，并从数值解中得到系统的上升时间、峰值时间、调节时间等时域指标。

解：

t=40;%仿真时长b=0.5;w=10;T = 5; %%%%%%%%%%%采用矩阵方法的四阶龙格库塔法%%%%%dt=0.05;%积分步长N=