MATLAB基础与控制系统CAD仿真练习1(11)

控制系统计算机辅助设计,matlab仿真的一些小练习

3-5．分别使用解微分方程方法、控制工具箱、simulink求解具有如下闭环传递函数的系统的阶跃响应。

10

(s) 4

32

s 8s 36s 40s 10

解答：

（1）解微分方程方法：

首先将 (s)转化为状态方程，在Command Window中输入命令行： num=[10];

den=[1 8 36 40 10];

[A B C D]=tf2ss(num,den)

得到结果：

A =

-8 -36 -40 -10 1 0 0 0 0 1 0 0 0 0 1 0 B = 1 0 0 0 C =

0 0 0 10 D = 0

从而得到状态方程：

8 36

0dx=10100y=[

40 101

00 x+0 u

00 0

01 0

] x

建立M-file：ztfc.m

function dx=ztfc(t,x)

dx=[-8*x(1)-36*x(2)-40*x(3)-10*x(4)+1;x(1);x(2);x(3)];

又建立M-file:zuoye251.m（调用ztfc.m）

[t,x]=ode45('ztfc',[0,18],[0;0;0;0]); y=10*x(:,4);

plot(t,y),grid on;