- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
自动控制原理实验报告(线性系统的频域分析)
物电学院 电气12(1)班 徐楠 一、实验目的
1.掌握用MATLAB语句绘制各种频域曲线。
2.掌握控制系统的频域分析方法。
二、基础知识及MATLAB函数
频域分析法是应用频域特性研究控制系统的一种经典方法。它是通过研究系统对正弦信号下的稳态和动态响应特性来分析系统的。采用这种方法可直观的表达出系统的频率特性,分析方法比较简单,物理概念明确。
1.频率曲线主要包括三种:Nyquist图、Bode图和Nichols图。
1)Nyquist图的绘制与分析
MATLAB中绘制系统Nyquist图的函数调用格式为:
nyquist(num,den) 频率响应w的范围由软件自动设定
nyquist(num,den,w) 频率响应w的范围由人工设定
[Re,Im]= nyquist(num,den) 返回奈氏曲线的实部和虚部向量,不作图
2)Bode图的绘制与分析
系统的Bode图又称为系统频率特性的对数坐标图。Bode图有两张图,分别绘制开环频率特性的幅值和相位与角频率的关系曲线,称为对数幅频特性曲线和对数相频特性曲线。
MATLAB中绘制系统Bode图的函数调用格式为:
bode(num,den) 频率响应w的范围由软件自动设定
bode(num,den,w) 频率响应w的范围由人工设定
[mag,phase,w]=bode(num,den,w) 指定幅值范围和相角范围的伯德图
2.幅值裕量和相位裕量
幅值裕量和相位裕量是衡量控制系统相对稳定性的重要指标,需要经过复杂的运算求取。应用MATLAB功能指令可以方便地求解幅值裕量和相位裕量。
其MATLAB调用格式为:
[Gm,Pm,Wcg,Wcp]=margin(num,den)
其中,Gm,Pm分别为系统的幅值裕量和相位裕量,而Wcg,Wcp分别为幅值裕量和相位裕量处相应的频率值。
另外,还可以先作bode图,再在图上标注幅值裕量Gm和对应的频率Wcg,相位裕量Pm和对应的频率Wcp。其函数调用格式为:
margin(num,den)
三、实验内容
1.典型二阶系统
绘制出,,0.3,0.5,0.8,2的bode图,记录并分析对系统bode图的影响。
num=36;
den1=[1 1.2 36];
den2=[1 3.6 36];
den3=[1 6 36];
den4=[1 9.6 36];
den5=[1 24 36];
w=logspace(-2,3,100);
bode(num,den1,w)
grid
hold on
bode(num,den2,w)
bode(num,den3,w)
bode(num,den3,w)
bode(num,den4,w)
bode(num,den5,w)
结果分析:
从图中可看出ζ越小,中频段振荡越剧烈。该二阶系统是典型的振荡环节,当0ζ/2时, ωr(谐振频率), Mr(谐振峰值)均为ζ的减函数。ζ越小ωr,Mr越大。振荡幅度越大,超调量越大,过程越不平稳且系统响应速度越慢.当/2ζ1时。A(ω)单调减小,此时无谐振峰值和谐振频率,过程较平稳,即随着ζ(阻尼比)的增大,bode图在穿越频率处的尖峰越明显,此处用渐近线代替时误差越大。
2.系统的开环传递函数为
(1)
(2)
(3)
绘制系统的Nyquist曲线、Bode图,说明系统的稳定性,并通过绘制阶跃响应曲线验证。
bode图
num=10;
den=[5 26 5 0 0];
w=logspace(-2,3,100);
bode(num,den,w)
grid
Nyquist图
num=10;
den=[5 26 5 0 0];
w=logspace(-2,3,100);
nyquist(num,den)
v=[-0.5,0.5,-0.5,0.5]
axis(v )
grid
局部放大图:
阶跃响应曲线:
G=tf(10,[5,26,5,0,0]);
G_c=feedback(G,1);
t=0:1:100;
step(G_c,t)
bode图
num=100;
den=[1 8 17 10];
w=logspace(-2,3,100);
bode(num,den,w)
grid
Nyquist图:
um=100;
den=[1 8 17 10];
w=logspace(-2,3,100);
nyquist(num,den)
v=[-0.5,0.5,-0.5,0.5]
axis(v )
grid
局部放大图:
阶跃响应曲线:
num=100;
den=[1 8 17 10];
w=logspace(-2,3,100);
G=tf(nu
文档评论(0)