《计算机仿真》上机实验报告.pdfVIP

《计算机仿真》

上机实验报告

姓名：______________

学号：—

专业：测控技术与仪器

班级：⑵班_________

实验常微分方程的求解及系统数学模型的转换

.实验目的

通过实验熟悉计算机仿真中常用到的Matlab指令的使用方法，掌握常微分

方程求解指令和模型表示及转换指令，为进步从事有关仿真设计和研究工作打

下基础。

二.实验设备

个人计算机，Matlab软件。

三.实验准备

预习本实验有关内容(如教材第2、3、5章中的相应指令说明和例题)，编

写本次仿真练习题的相应程序。

四.实验内容

1.Matlab中常微分方程求解指令的使用

题目：请用MATLAB的0DE45算法分别求解下列二个方程。要求：1.编写出Matlab

仿真程序；2.画出方程解的图形并对图形进行简要分析；3.分析下列二个方程的

关系，

e[x=^x2\x=x2

=-

k)iU=i

1、仿真程序

functionfl=f1(t,x)

fl=-x2;

,,

[t,x]=ode45(fl,[0,20],[1]);

plot(t,x);

grid

functionf2=f2(t,x)

f2=x.2;

[t,x]=ode45Cf2\[0,20],[-1])；

figure(2);plot(t,x；;

grid

2.方程解的图形图形进行简要分析

3.二个方程的关系

两函数关于X轴对称

题目二：下面方程组用在人口动力学中，可以表达为单化的捕食者-被捕食者

模式（例如，狐狸和兔子）。其中阳表示被捕食者，当表示捕食者。如果被捕

食者有无限的食物，并且不会出现捕食者。于是有王=王，则这个式子是以指数

形式增长的。大量的被捕食者将会使捕食者的数量增长；同样，越来越少的捕食

者会使被捕食者的数量增长。而且，人口数量也会增长。请分别调用0DE45、0DE23

算法求解下面方程组。要求编写出Matlab仿真程序、画出方程组解的图形并对图

形进行分析和比较。

Xi=Xi-01XiX2+001r

002004r

X?=F+X1占+

耳(0)=30

%(0)=20

MATLAB仿真程序:

functionf3=f3(t,x)

f3=[x(l)-0.1*x(l)*x(2)+0.01*t;-x(2)+0.02*x(1)*x(2)+0.0*t]

,

[t,x]=odc5(f3,,[0,20],[30;20]);

plot(t,x);

arid

方程组图形:

2.Matlab中模型表示及模型转换指令的使用

题目三：若给定系统的的传递函数为

32

…、6.v+12.v+65+10

G(s)=—7；；

2

/+2?+35+h-+l

请用MATLAB编程求解其系统的极零点模型。

MATLAB仿真程序:

num=[612610];

den=[l2311];

sys=zpk(num,den,k)

程序运行结果:

num=：612610];