- 1、本文档共11页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
PAGE
PAGE 1
实验六 基于MATLAB控制系统的Nyquist图及其稳定性分析
一、实验目的
1、熟练掌握使用MATLAB命令绘制控制系统Nyquist图的方法。
2、能够分析控制系统Nyquist图的基本规律。
3、加深理解控制系统乃奎斯特稳定性判据的实际应用。
4、学会利用奈氏图设计控制系统。
二、实验原理
奈奎斯特稳定性判据(又称奈氏判据)
反馈控制系统稳定的充分必要条件是当从变到时,开环系统的奈氏曲线不穿过点且逆时针包围临界点点的圈数R等于开环传递函数的正实部极点数。奈奎斯特稳定性判据是利用系统开环频率特性来判断闭环系统稳定性的一个判据,便于研究当系统结构参数改变时对系统稳定性的影响。
1、对于开环稳定的系统,闭环系统稳定的充分必要条件是:开环系统的奈氏曲线不包围点。反之,则闭环系统是不稳定的。
2、对于开环不稳定的系统,有个开环极点位于右半平面,则闭环系统稳定的充分必要条件是:当从变到时,开环系统的奈氏曲线逆时针包围点次。
三、实验内容
1、绘制控制系统Nyquist图
例1、系统开环传递函数,绘制其Nyquist图。
M-file
clc
clear all
den=[10];
num=[1 2 10];
sys=tf(den,num)
nyquist(sys);
2、根据奈氏曲线判定系统的稳定性
例2、已知
绘制Nyquist图,判定系统的稳定性。
M-file
clc
clear
den=[0.5];
num=[1 2 1 0.5];
sys=tf(den,num);
nyquist(sys)
roots(num)
ans =
-1.5652
-0.2174 + 0.5217i
-0.2174 - 0.5217i
【分析】由于系统奈氏曲线没有包围且远离(-1,j 0)点,且p=0,因此系统闭环稳定。
四、实验能力要求
1、熟练使用MATLAB绘制控制系统Nyquist曲线的方法,掌握函数nyquist ( )的三种调用格式,并灵活运用。
2、学会处理奈氏图形,使曲线完全显示ω从-∞变化至+∞的形状。
3、熟练应用奈氏稳定判据,根据Nyquist图分析控制系统的稳定性。
4、改变系统开环增益或零极点,观察系统Nyquist图发生的变化以及系统稳定性的影响。
实验七 基于MATLAB控制系统的伯德图及其频域分析
一、实验目的
1、熟练掌握运用MATLAB命令绘制控制系统伯德图的方法。
2、了解系统伯德图的一般规律及其频域指标的获取方法。
3、熟练掌握运用伯德图分析控制系统稳定性的方法。
二、实验原理
对数频率稳定性判据的内容为:
闭环系统稳定的充分必要条件是当从零变化到时,时,在开环系统对数幅频特性曲线的频段内,相频特性穿越的次数为。其中 ,为正穿越次数,为负穿越次数,为开环传递函数的正实部极点数。
1、相角裕度
对于闭环稳定系统,如果开环相频特性再滞后γ度,则系统将变为临界稳定。当 γ 0时,相角裕度为正,闭环系统稳定。当 γ = 0 时,表示奈氏曲线恰好通过点,系统处于临界稳定状态。当 γ 0 时,相角裕度为负,闭环系统不稳定。
2、幅值裕度
对于闭环稳定系统,如果系统开环幅频特性再增大h倍,则系统将变为临界稳定状态。当h (dB) 0时 ,闭环系统稳定。当h (dB) = 0时,系统处于临界稳定状态。当h (dB) 0 ,闭环系统不稳定。
三、实验内容
1、绘制连续系统的伯德图
例1、已知控制系统开环传递函数,绘制其Bode图。
M-file
clc
clear
den=[10];
num=[1 2 10];
sys=tf(den,num)
bode(sys);
clc
clear
den=[10];
num=[1 2 10];
sys=tf(den,num)
margin(sys);
2、系统对数频率稳定性分析
r = 0 -10.0000 -2.0000例2
r =
0
-10.0000
-2.0000
令K=1时,根据跟轨迹可知K=12时临界增益,则
M-file
clc
clear
k=1;
den=[k];
num=conv([1 0],conv([0.5 1],[0.1 1]));
sys=tf(den,num);
margin(sys);
r=roo
文档评论(0)