中国矿业大学《控制系统计算机仿真》实验试题及仿真程序及结果
实验一 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