matlab仿真实例Word格式.docx
- 文档编号:21959015
- 上传时间:2023-02-01
- 格式:DOCX
- 页数:13
- 大小:198.88KB
matlab仿真实例Word格式.docx
《matlab仿真实例Word格式.docx》由会员分享,可在线阅读,更多相关《matlab仿真实例Word格式.docx(13页珍藏版)》请在冰豆网上搜索。
零极点图'
)
由图可知:
在S右半平面有极点,因此可知系统是不稳定的。
2、用MATLAB求出G(s)=(s^2+2*s+2)/(s^4+7*s^3+5*s+2)的极点。
程序及结果:
sys=tf([1,2,2],[1,7,3,5,2]);
p=pole(sys)
p=
-6.6553
0.0327+0.8555i
0.0327-0.8555i
-0.4100
(二)阶跃响应
1、二阶系统G(s)=10/s^2+2*s+10
1)键入程序,观察并记录单位阶跃响应曲线:
sys=tf(10,[1,2,10]);
step(sys);
G(s)=10/s^2+2*s+10单位阶跃响应曲线'
2)计算系统闭环跟、阻尼比、无阻尼振荡频率,并记录:
程序及结果:
-1.0000+3.0000i
-1.0000-3.0000i
[wn,z]=damp(sys)
wn=
3.1623
z=
0.3162
3)记录实际测取的峰值大小,峰值时间和过渡过程时间,并填表:
实际值理论值
峰值Cmax1.35s
峰值时间tp1.05s
过渡时间+5%3.54s
ts+2%3.18s
4)修改参数,分别实现ξ=1和ξ=2的响应曲线,并记录:
ξ=1:
>
zeta=1;
wn=sqrtm(10)
sys=tf(10,[1,2*wn*zeta,10]);
step(sys)
ξ=1响应曲线'
)
(图见下页)
ξ=2:
zeta=2;
ξ=2响应曲线'
(曲线见下页)
5)修改参数,分别实现wn1=wn/2和wn2=2*wn的响应曲线,并记录:
wn1=w0/2:
zeta=2/(wn*2)
zeta=
wn1=wn/2
wn1=
1.5811
sys=tf((wn1)^2,[1,2*wn1*zeta,(wn1)^2]);
wn1=wn/2响应曲线'
)(曲线见下页)
wn2=2*wn:
wn2=2*wn
wn2=
6.3246
sys=tf((wn2)^2,[1,2*wn2*zeta,(wn2)^2]);
wn2=wn*2响应曲线'
2、作出以下系统的阶跃响应曲线,并与原系统响应曲线进行比较,作出相应实验分析结果。
1)G1(s)=(2*s+10)/(s^2+2*s+10)
step(sys)
holdon
sys=tf([2,10],[1,2,10]);
gtext('
\leftarrowG(s)'
);
\leftarrowG1(s)'
G1(s)与G(s)阶跃响应曲线'
实验分析结果:
G1(s)与原系统响应曲线相比,峰值增加,峰值时间、
上升时间、调节时间提前,最终稳定值相等。
2)G2(s)=(s^2+0.5*s+10)/(s^2+2*s+10)
sys=tf([1,0.5,10],[1,2,10]);
gtext('
\leftarrowG2(s)'
G2(s)与G(s)阶跃响应曲线'
G2(s)与原系统响应曲线相比,峰值减小,峰值时间增加、
上升时间减小、调节时间增加,最终稳定值相等。
3)G3(s)=(s^2+0.5*s)/(s^2+2*s+10)
holdon
sys=tf([1,0.5],[1,2,10]);
\leftarrowG3(s)'
G3(s)与G(s)阶跃响应曲线'
G3(s)与原系统响应曲线相比,峰值减小,峰值时间减小、
上升时间减小、调节时间增加,最终稳定值不相等。
4)G4(s)=s/(s^2+2*s+10)
step(sys)>
\leftarrowG4(s)'
holdon>
G4(s)与G(s)阶跃响应
sys=tf(1,[1,2,10]);
曲线
G4(s)与原系统响应曲线相比,峰值减小,峰值时间、上升时间、
调节时间都相等,最终稳定值不相等。
、单位阶跃响应:
3
C(s)/R(s)=25/(s^2+4*s+25)求该系统单位阶跃响应曲线,并在所得图形上加网格和标题:
sys=tf(25,[1,4,25]);
gridon;
C(s)/R(s)=25/(s^2+4*s+25)单位阶跃响应曲线'
(图见下页图一)
(二)系统动态特性分析
用MATLAB求二阶系统G(s)=120/(s^2+12*s+120)和
G(s)=0.01/(s^2+0.002*s+0.01)的峰值时间tp上升时间tr调整时间ts超调量
σ%。
G(s)=120/(s^2+12*s+120):
sys=tf(120,[1,12,120]);
G(s)=120/(s^2+12*s+120)单位阶跃响应曲线'
(曲线见下页图二)
峰值时间tp=0.34s上升时间tr=0.158s调整时间ts=0.532s超
调量σ%=12.8%
G(s)=0.01/(s^2+0.002*s+0.01)
sys=tf(0.01,[1,0.002,0.01]);
title('
G(s)=0.01/(s^2+0.002*s+0.01)单位阶跃响应曲线'
(图见图三)
峰值时间tp=32s上升时间tr=10.3s调整时间ts=3.9e+003s
超调量σ%=96.7%
图一
图二
图三
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 仿真 实例
![提示](https://static.bdocx.com/images/bang_tan.gif)