计算机仿真技术实验报告.docx
- 文档编号:5411117
- 上传时间:2022-12-16
- 格式:DOCX
- 页数:22
- 大小:321.79KB
计算机仿真技术实验报告.docx
《计算机仿真技术实验报告.docx》由会员分享,可在线阅读,更多相关《计算机仿真技术实验报告.docx(22页珍藏版)》请在冰豆网上搜索。
计算机仿真技术实验报告
计算机仿真技术实验报告
学院:
专业:
班级:
姓名:
学号:
实验一 常微分方程的求解及系统数学模型的转换
一.实验目的
通过实验熟悉计算机仿真中常用到的Matlab指令的使用方法,掌握常微分方程求解指令和模型表示及转换指令,为进一步从事有关仿真设计和研究工作打下基础。
二.实验设备
个人计算机,Matlab软件。
三.实验准备
预习本实验有关内容(如教材第2、3、5章中的相应指令说明和例题),编写本次仿真练习题的相应程序。
四.实验内容
1.Matlab中常微分方程求解指令的使用
题目一:
请用MATLAB的ODE45算法分别求解下列二个方程。
要求:
1.编写出Matlab仿真程序;2.画出方程解的图形并对图形进行简要分析;3.分析下列二个方程的关系。
1.2.
程序代码
方程一:
functionfun=fun1(t,x)
fun=-x^2;
[t,x]=ode45('fun1',[0,1],1);plot(t,x,'o',t,x,'-');xlabel('timet0=0,tt=20');ylabel('xvaluesx(0)=1');grid
执行结果
使用ode45
方程二:
functionfun=fun2(t,x)
fun=x^2;
[t,x]=ode45('fun2',[0,1],-1);plot(t,x,'o',t,x,'-');xlabel('timet0=0,tt=1');ylabel('xvaluesx(0)=-1');grid
执行结果
使用ode45
题目二:
下面方程组用在人口动力学中,可以表达为单一化的捕食者-被捕食者模式(例如,狐狸和兔子)。
其中
表示被捕食者,
表示捕食者。
如果被捕食者有无限的食物,并且不会出现捕食者。
于是有
,则这个式子是以指数形式增长的。
大量的被捕食者将会使捕食者的数量增长;同样,越来越少的捕食者会使被捕食者的数量增长。
而且,人口数量也会增长。
请分别调用ODE45、ODE23算法求解下面方程组。
要求编写出Matlab仿真程序、画出方程组解的图形并对图形进行分析和比较。
程序代码
functionfun=fun3(t,x)
fun=[x
(1)-0.1*x
(1)*x
(2)+0.01*t;-x
(2)+0.02*x
(1)*x
(2)+0.04*t];
[t,x]=ode45('fun3',[0,20],[30,20]);
plot(t,x);xlabel('timet0=0,tt=20');ylabel('xvaluesx1(0)=30,x2(0)=20');grid
执行结果
[t,x]=ode23('fun3',[0,20],[30,20]);
plot(t,x);xlabel('timet0=0,tt=20');ylabel('xvaluesx1(0)=30,x2(0)=20');grid
2.Matlab中模型表示及模型转换指令的使用
题目三:
若给定系统的的传递函数为
请用MATLAB编程求解其系统的极零点模型。
程序代码
num=[612610];
den=[12311];
sys=tf(num,den);
H=zpk(sys)
执行结果:
Zero/pole/gain:
6(s+1.929)(s^2+0.07058s+0.8638)
-------------------------------------------------
(s^2+0.08663s+0.413)(s^2+1.913s+2.421)
题目四:
习题2.4
系统传递函数为
,求其对角标准型实现。
程序代码:
num=[145];
den=[16116];
sys=tf(num,den);
canon(sys,'modal')
执行结果:
a=
x1x2x3
x1-300
x20-20
x300-1
b=
u1
x1-15.11
x2-22.98
x38.958
c=
x1x2x3
y1-0.066190.043520.1116
d=
u1
y10
Continuous-timemodel.
题目五:
习题5.8
已知某离散系统
,采样周期
,先需要在
下作数字仿真,求该系统的数字仿真模型。
分步执行:
D=zpk(0,[0.10.2],1,0.1)
Zero/pole/gain:
z
---------------
(z-0.1)(z-0.2)
Samplingtime:
0.1
D2=d2d(D,0.02)
Zero/pole/gain:
0.278(z-0.4926)
--------------------
(z-0.631)(z-0.7248)
Samplingtime:
0.02
直接执行:
D=zpk(0,[0.10.2],1,0.1);
D2=d2d(D,0.02)
Zero/pole/gain:
0.278(z-0.4926)
--------------------
(z-0.631)(z-0.7248)
Samplingtime:
0.02
五.总结与体会
通过本次实验了解了计算机仿真中常用到的Matlab指令的使用方法,掌握了常微分方程求解指令和模型表示及转换指令,为进一步学习仿真设计打下了基础。
实验二Matlab优化工具箱的使用
一.实验目的
通过上机操作熟悉Matlab优化工具箱的主要功能及其使用方法,掌握优化工具箱中常用函数的功能和语法,并利用其进行极值运算、求解线性和非线性问题等,为进一步的仿真设计和研究打下基础。
二.实验设备
个人计算机,Matlab软件。
三.实验准备
预习本实验有关内容(如教材第6章中的相应指令说明和例题),编写本次仿真练习题的相应程序。
四.实验内容
1.应用Matlab优化工具箱求解优化问题
例题6.6~6.10,选做2题,要求自行修改方程系数,比较运行结果。
6.6(原)求解非线性方程
初始值为x0=[-5-5]。
原程序
functionF=myfun(x)
F=[2*x
(1)-x
(2)-exp(-x
(1));
-x
(1)+2*x
(2)-exp(-x
(2))];
x0=[-5,-5];
>>options=optimset('Display','iter');
>>[x,fval]=fsolve(@myfun,x0,options)
NormofFirst-orderTrust-region
IterationFunc-countf(x)stepoptimalityradius
1347071.22.29e+0041
2612003.415.75e+0031
393147.0211.47e+0031
412854.45213881
515239.52711071
61867.0412130.81
72116.704219.051
8242.4278812.261
9270.0.0.2062.5
10307.03149e-0060.0.002942.5
11333.29525e-0130.6.36e-0072.5
Optimizationterminatedsuccessfully:
First-orderoptimalityislessthanoptions.TolFun.
x=
0.56710.5671
fval=
1.0e-006*
-0.4059
-0.4059
改变后
functionF=myfun(x)
F=[3*x
(1)-x
(2)-exp(-x
(1));
-x
(1)+3*x
(2)-exp(-x
(2))];
x0=[-5,-5];
options=optimset('Display','iter');
[x,fval]=fsolve(@myfun,x0,options)
NormofFirst-orderTrust-region
IterationFunc-countf(x)stepoptimalityradius
1350189.52.38e+0041
2613370.616.15e+0031
393741.711.65e+0031
4121109.0314661
515344.01811411
618105.246145.91
72126.6064115.11
8242.654413.521
9270.0.0.182.5
10309.03107e-0080.0.2.5
11339.44087e-0180.5.87e-0092.5
Optimizationterminatedsuccessfully:
First-orderoptimalityislessthanoptions.TolFun.
x=
0.35170.3517
fval=
1.0e-008*
-0.2173
-0.2173
例题6.8
functionf=fun(x)
f=exp(x
(1))*(4*x
(1)^2+2*x
(2)^2+4*x
(1)*x
(2)+2*x
(2)+1);
x0=[-11];
options=[];
[x,options]=fminu('fun',x0,options)
Warning:
ThefunctionFMINUisobsoleteandhasbeenreplacedby
FMINUNC.FMINUwillberemovedinafuturereleaseofthe
OptimizationToolbox.UpdateyourcodetocallFMINUNCto
suppressthiswarningmessage.See"Convertingyourcodeto
Version2Syntax"intheOptimizationToolboxUser'sGuide
formoreinformation.
>InC:
\MATLAB6p5\toolbox\optim\fminu.matline48
x=
0.5000-1.0000
options=
Columns1through11
00.00010.00010.00000000.0000036.00000
Columns12through18
00200.000000.00000.10001.0000
改变后:
functionf=fun(x)
f=exp(x
(1))*(5*x
(1)^2+3*x
(2)^2+2*x
(1)*x
(2)+4*x
(2)+1);
x0=[-11];
options=[];
[x,options]=fminu('fun',x0,options)
Warning:
ThefunctionFMINUisobsoleteandhasbeenreplacedby
FMINUNC.FMINUwillberemovedinafuturereleaseofthe
OptimizationToolbox.UpdateyourcodetocallFMINUNCto
suppressthiswarningmessage.See"Convertingyourcodeto
Version2Syntax"intheOptimizationToolboxUser'sGuide
formoreinformation.
>InC:
\MATLAB6p5\toolbox\optim\fminu.matline48
x=
0.1878-0.7292
options=
Columns1through11
00.00010.00010.0000000-0.5057038.00000
Columns12through18
00200.000000.00000.10001.0005
2.应用Matlab优化工具箱求解极值问题
已知函数f(x)=10*exp(-x)*cos(x),求函数的极值。
(1)x∈[2,5]时,求函数的最小值,并画出函数的曲线。
X=fmin('10*exp(-x)*cos(x)',2,5)
X=2.3562
f='10*exp(-x)*cos(x)';
fplot(f,[2,5])
Y=-0.6702
(2)x∈[3,9]时,求函数的最大值,并画出函数的曲线。
X=fmin('-10*exp(-x)*cos(x)',3,9)
X=5.4978
f='10*exp(-x)*cos(x)';
fplot(f,[3,9])
Y=0.0290
五.总结与体会
通过上机操作熟悉Matlab优化工具箱的主要功能及其使用方法,掌握优化工具箱中常用函数的功能和语法,并利用其进行极值运算、求解线性和非线性问题等,为进一步的仿真设计和研究打下基础。
实验三利用Matlab和Simulink进行系统仿真设计
一.实验目的
通过实验对一个汽车运动控制系统进行实际设计与仿真,掌握控制系统性能的分析和仿真处理过程,熟悉用Matlab和Simulink进行系统仿真的基本方法。
二.实验设备
个人计算机,Matlab软件。
三.实验准备
预习本实验相关说明,复习PID控制器的原理和作用,明确汽车运动控制系统问题的描述及其模型表示,编写本次仿真练习的相应程序。
四.实验说明
本实验是对一个汽车运动控制系统进行实际设计与仿真,其方法是先对汽车运动控制系统进行建摸,然后对其进行PID控制器的设计,建立了汽车运动控制系统的模型后,可采用Matlab和Simulink对控制系统进行仿真设计。
注意:
设计系统的控制器之前要观察该系统的开环阶跃响应,采用阶跃响应函数step()来实现,如果系统不能满足所要求达到的设计性能指标,需要加上合适的控制器。
然后再按照仿真结果进行PID控制器参数的调整,使控制器能够满足系统设计所要求达到的性能指标。
五.实验内容
1.问题的描述
如下图所示的汽车运动控制系统,设该系统中汽车车轮的转动惯量可以忽略不计,并且假定汽车受到的摩擦阻力大小与汽车的运动速度成正比,摩擦阻力的方向与汽车运动的方向相反,这样,该汽车运动控制系统可简化为一个简单的质量阻尼系统。
根据牛顿运动定律,质量阻尼系统的动态数学模型可表示为:
系统的参数设定为:
汽车质量m=1000kg,
比例系数b=50N·s/m,
汽车的驱动力u=500N。
根据控制系统的设计要求,当汽车的驱动力为500N时,汽车将在5秒内达到10m/s的最大速度。
由于该系统为简单的运动控制系统,因此将系统设计成10%的最大超调量和2%的稳态误差。
这样,该汽车运动控制系统的性能指标可以设定为:
上升时间:
tr<5s;
最大超调量:
σ%<10%;
稳态误差:
essp<2%。
2.系统的模型表示
假定系统的初始条件为零,则该系统的Laplace变换式为:
即
则该系统的传递函数为:
如果用Matlab语言表示该系统的传递函数模型,相应的程序代码如下:
num=1;den=[1000,50];sys=tf(num,den)
同时,系统的数学模型也可写成如下的状态方程形式:
如果用Matlab语言表示该系统状态空间模型,相应的程序代码如下:
A=-50/1000;B=1/1000;C=1;D=0;sys=ss(A,B,C,D)
3.系统的仿真设计
●利用Matlab进行仿真设计
I.求系统的开环阶跃响应
在Matlab命令窗口输入相应的程序代码,得出该系统的模型后,接着输入下面的指令:
step(u*sys)
可得到该系统的开环阶跃响应曲线,如下图所示:
从图上可看出该系统不能满足系统设计所要求达到的性能指标,需要加上合适的控制器。
II.PID控制器的设计
PID控制器的传递函数为:
在PID控制中,比例(P)、积分(I)、微分(D)这三种控制所起的作用是不同的(请注意在实验总结中进行归纳)。
下面分别讨论其设计过程。
(1)比例(P)控制器的设计
增加比例控制器之后闭环系统的传递函数为:
由于比例控制器可以改变系统的上升时间,现在假定Kp=100,观察一下系统的阶跃响应。
在MATLAB命令窗口输入指令:
num=100;den=[1000,150];sys=tf(num,den);step(500*sys)
由此仿真结果,分析系统的稳态值是否满足设计要求,系统的稳态误差和上升时间能不能满足设计要求?
上升时间为40s远远大于5s,不能满足设计要求,稳态误差为(500-333)/500远远大于2%,因此系统不满足设计要求
若减小汽车的驱动力为10N,重新进行仿真,仿真结果为:
num=100;den=[1000,150];sys=tf(num,den);step(10*sys)
如果所设计的比例控制器仍不能满足系统的稳态误差和上升时间的设计要求,则可以通过提高控制器的比例增益系数来改善系统的输出。
例如把比例增益系数Kp从100提高到10000重新计算该系统的阶跃响应,结果为:
程序:
num=10000;den=[1000,10050];sys=tf(num,den);step(10*sys)
此时系统的稳态误差接近为零,系统上升时间也降到了0.5s以下。
这样做虽然满足了系统性能要求,但实际上该控制过程在现实中难以实现。
因此,引入比例积分(PI)控制器来对系统进行调节。
(2)比例积分(PI)控制器的设计
采用比例积分控制的系统闭环传递函数可表示为:
增加积分环节的目的是减小系统的稳态误差,假设比例系数Kp=600,积分系数KI=1,编写相应的MATLAB程序代码如下:
num=[600,1];den=[1000,650,1];sys=(num,den)
可以调节控制器的比例和积分系数来满足系统的性能要求。
例如选择比例系数KP=800,积分系数KI=40时,可得系统阶跃响应曲线为:
num=[800,40];den=[1000,850,40];sys=(num,den);step(500*sys)
可见,此时的控制系统已经能够满足系统要求达到的性能指标设计要求。
但此控制器无微分项,而对于有些实际控制系统往往需要设计完整的PID控制器,以便同时满足系统的动态和稳态性能要求。
(3)比例积分微分(PID)控制器的设计
采用PID控制的系统闭环传递函数为:
假设该控制器的比例系数KP=1,积分系数KI=1,微分系数KD=1,编写MATLAB程序代码如下:
num=[1,1,1];den=[1001,51,1];sys=(num,den)
运行上述程序,并且调整PID控制器的控制参数,直到控制器满足系统设计的性能指标要求为止。
num=[10,650,50];den=[1010,700,50];sys=tf(num,den);step(500*sys)
最后,选择KP=650,KI=50,KD=10,此时系统的阶跃响应曲线如下:
从图中可以看出该系统能够满足设计的总体性能要求。
●利用Simulink进行仿真设计
I.求系统的开环阶跃响应
利用Simulink建立系统阶跃响应模型,如下图所示。
双击Step模块,设置模块属性:
跳变时间为0;初始值为0;终止值为10;采样时间为0。
单击◢按钮开始仿真,双击Scope模块,可得系统阶跃响应曲线。
分析此时系统的性能指标是否满足设计要求?
若不能则应加上合适的控制器。
答:
不满足
II.PID控制器的设计
在Simulink的模型窗口建立一个包含PID控制器的闭环系统阶跃响应模型,如图所示:
分别双击Kp、Ki、Kd模块设定比例、积分、微分系数,点击◢按钮开始仿真,双击Scope模块,观察系统的阶跃响应曲线,直到满足要求为止。
最终选取KP=680,KI=50,KD=15,此时控制器能满足系统设计所要求达到的性能指标,Simulink仿真的汽车运动PID控制系统的阶跃响应曲线如下:
六.总结与体会
实验总结:
(总结实验中遇到的问题及解决办法,分析PID控制器三种控制所起的作用等。
)
通过实验对一个汽车运动控制系统进行实际设计与仿真,掌握控制系统性能的分析和仿真处理过程,熟悉用Matlab和Simulink进行系统仿真的基本方法。
本实验是对一个汽车运动控制系统进行实际设计与仿真,其方法是先对汽车运动控制系统进行建摸,然后对其进行PID控制器的设计,建立了汽车运动控制系统的模型后,可采用Matlab和Simulink对控制系统进行仿真设计。
设计系统的控制器之前要观察该系统的开环阶跃响应,采用阶跃响应函数step()来实现,如果系统不能满足所要求达到的设计性能指标,需要加上合适的控制器。
然后再按照仿真结果进行PID控制器参数的调整,使控制器能够满足系统设计所要求达到的性能指标。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算机仿真 技术 实验 报告