控制系统计算机仿真试题,中国矿业大学《控制系统计算机仿真》实验试题及仿真程序及结果...

中国矿业大学《控制系统计算机仿真》实验试题及仿真程序及结果

实验一 1、(教材P37页 例2.1)已知一阶微分方程：

y 30y,y(0) 1/3,

h 0.05

试编写程序，用欧拉法求 ，即 y30的值；并画出图形。 y(t 1.5)解：程序如下： clc clear

Tf=input('仿真时间Tf='); %输入仿真时间 h=input('计算步长h='); %输入仿真步长 x0=1/3; t=0; x=x0; %输入初始值 for i=1:Tf/h K1=-30*x0;

x1=x0+h*K1;

x=[x;x1]; %变量以向量形式保存 t=[t;t(i)+h]; %对应时刻以向量形式保存 x0=x1; end

[t,x] %以数据形式输出 plot(t,x) %以曲线形式输出 所绘图形如下：

0.5

0.4

0.3

0.2

0.1

-0.1

-0.2

00.20.40.60.811.21.41.6

1 2、2. (教材P79习题2.5)已知系统的状 x

态方程和输出方程为：

x式中 u(t)=1(t)。初始条件为： 2

x1(0)=x2(0)=x3(0)=0。取h=0.05， 3 x 试用RK4法求t=0.5时的y 值。

解：程序如下：

8 1 0 x1 0 19 0 1 x2 4u 12 0 0 10 x3

x1

y [1 0 0]x2