MATLAB课程设计.docx
- 文档编号:25660633
- 上传时间:2023-06-11
- 格式:DOCX
- 页数:10
- 大小:86.17KB
MATLAB课程设计.docx
《MATLAB课程设计.docx》由会员分享,可在线阅读,更多相关《MATLAB课程设计.docx(10页珍藏版)》请在冰豆网上搜索。
MATLAB课程设计
一、设计任务
已知单位负反馈系统被控对象的传递函数为
使用MATLAB工具箱,设计系统的超前校正器
,要求:
(1)在斜坡信号
作用下,系统的稳态误差
;
(2)校正后系统的相位裕度
范围为:
;
(3)绘制系统校正前后的Bode图和阶跃响应曲线,并作简要分析。
二、设计目的
1.通过课程设计使学生更进一步掌握自动控制原理课程的有关知识,加深对内涵的理解,提高解决实际问题的能力。
2.理解自动控制原理中的关于开环传递函数,闭环传递函数的概念以及二者之间的区别和联系。
3.理解在自动控制系统中对不同的系统选用不同的校正方式,以保证得到最佳的系统。
4.理解在校正过程中的静态速度误差系数,相角裕度,截止频率,超前(滞后)角频率,分度系数,时间常数等参数。
5.学习MATLAB在自动控制中的应用,会利用MATLAB提供的函数求出所需要得到的实验结果。
6.从总体上把握对系统进行校正的思路,能够将理论操作联系实际、运用于实际。
三、设计方案
所谓的校正,就是在系统中加入一些其参数可以改变的机构或装置,使系统的整个特性发生变化,从而满足给定的各项性能指标。
如果性能指标以单位阶跃响应的峰值时间、调节时间、超调量、阻尼比、稳态误差等时域特征量给出时,一般采用时域法校正;如果性能指标以系统的相角裕度、谐振峰值、闭环带宽、静态误差系数等频域特征量给出,如本题,一般采用频率法校正。
在频域内进行系统设计,是一种间接而又简单的设计方法,它虽然以伯德图的形式给出非严格意义上的系统动态性能,但却能方便的根据频域指标确定校正装置的参数,特别是对已校正系统的高频特性有要求时,采用频域校正法较其他方法更为方便。
一般来说,开环频率特性的低频段表征了闭环系统的稳态性能;开环频率特性的中频段表征了闭环系统的动态性能;高频段表征了闭环系统地复杂性和噪声抑制性能。
因此,用频域校正法设计控制系统的实质,就是在系统中加入频率特性形状合适的校正装置,使开环频率特性形状变成所期望的形状:
低频段增益充分大,以保证稳态误差的要求;中频段对数幅频特性斜率一般为-20dB/dec,并占据充分宽的频带,以保证具备适当的相角裕度;高频段增益尽快减小,以削弱噪声影响,若系统原有部分高频段已经符合该种要求,则校正时可保持高频段形状不变,以简化校正装置形式。
常用的校正形式有串联超前校正、串联滞后校正、串联滞后-超前校正。
每种方法都有不同的适用范围,应当根据实际要求恰当的选择,由于本题要求采用串联超前校正,下面将着重介绍这种方法。
3.1设计思路
频域法中的超前校正是利用校正装置的超前相位在穿越频率处对系统进行相位补偿,以提高系统的相位稳定裕量,同时也提高了穿越频率值,从而改善系统的稳定性快速性。
串联超前校正主要适用于稳定精度不需要改变(即稳态性能较好),动态性能不佳,而穿越频率附近相位变化平稳的系统。
3.2参数的一般设计步骤
(1)根据给定的系统性能指标,确定开环增益K。
(2)利用步骤1中求的的K绘制未校正系统的伯德图。
(3)在伯德图上量取未校正系统的相位裕度和幅值裕度,并计算为使相位裕度达到给定的指标所需补偿的超前相角
其中
为给定的相位裕度指标,
为未校正系统的相位裕度,
为附加的角度。
(4)取
,从而求出
求出a。
(5)取未校正系统的幅值为-10lga(dB)时的频率作为校正后系统的截止频率
。
(6)由
计算出参数T,并写出超前校正传递函数。
(7)检验指标:
绘制系统校正后的伯德图,检验是否满足给定的性能指标。
当系统仍不能满足要求时增大
值,从步骤3开始重新计算。
3.3超前校正的设计
1)由稳态误差要求,求系统开环增益K。
由于要求系统跟随r(t)=2t的斜坡输入时,稳态误差
所以,系统开环增益K
取K=1000
则待校正的系统的开环传递函数为:
2)绘制系统校正前的Bode图和单位阶跃响应曲线,并计算校正前系统的增益裕量Gm,相位裕量Pm。
绘制系统校正前Bode图的MATLAB程序如下:
s=tf('s');
G=1000/(s*(0.1*s+1)*(0.001*s+1));
margin(G)
上述程序运行后,得到系统校正前的Bode图,如图3-1所示。
图3-1系统校正前的Bode图
放大后可知,穿越频率及截止频率如下图3-2
图3-2穿越频率及截止频率
该被控对象构成单位负反馈系统后,绘制该系统校正前的单位阶跃响应曲线的MATLAB程序如下:
s=tf('s');
G=1000/(s*(0.1*s+1)*(0.001*s+1));
sys=feedback(G,1);
step(sys)
上述程序运行后,得到系统校正前的单位阶跃响应曲线,如图3-3所示。
图3-3系统校正前的单位阶跃响应曲线
由图3-1可知,该系统校正前的增益裕量Gm=0.0864dB,相位裕量Pm=0.0584°,显然不满足设计要求,系统超调过大,需要对系统进行超前校正。
3)求超前校正器的传递函数。
校正后系统的相位裕度范围为
,取Pm=50°,按超前校正器设计步骤,编写的MATLAB程序如下:
s=tf('s');
G=1000/(s*(0.1*s+1)*(0.001*s+1));
[mag,phase,w]=bode(G);
[Gm,Pm]=margin(G);
QWPm=50;
FIm=QWPm-Pm+5;
FIm=FIm*pi/180;
alfa=(1-sin(FIm))/(1+sin(FIm));
adb=20*log10(mag);am=10*log10(alfa);
wc=spline(adb,w,am);
T=1/(wc*sqrt(alfa));
alfat=alfa*T;
Gc=tf([T1],[alfat1])
上述程序部分语句注释:
[mag,phase,w]=bode(G);该句只计算Bode图上多个频率点w对应的幅值mag和相位phase,三个变量w,mag,phase均为有多个元素的矢量,该句不绘制Bode图。
[Gm,Pm]=margin(G);该句只计算Bode图的增益裕量Gm和相位裕量Pm,不绘制Bode图。
Gc=tf([T1],[alfat1]);该句是以[T1]为分子,[alfat1]为分母构成传递函数Gc。
上述程序运行后,得到超前校正器的传递函数为
4)校正后系统的开环传递函数为Gc(s)G(s),根据该开环传递函数绘制Bode图,其MATLAB程序如下:
s=tf('s');
G=1000/(s*(0.1*s+1)*(0.001*s+1));Gc=(0.01794*s+1)/(0.00179*s+1);
margin(Gc*G)
上述程序运行后,校正后系统的Bode图如图3-4所示。
图3-4校正后系统的Bode图
由图3-4可知,校正后系统的相位裕量Pm=48.2°,满足设计要求。
校正后系统的单位阶跃响应曲线如图3-5所示,系统经过超前校正后是稳定的,其MATLAB程序如下:
s=tf('s');
G=1000/(s*(0.1*s+1)*(0.001*s+1));
Gc=(0.01794*s+1)/(0.00179*s+1);
step(feedback(Gc*G,1))
图3-5校正后系统的单位阶跃响应曲线
心得体会
通过这次对控制系统的滞后校正的设计与分析,让我对串联超前校正环节有了更清晰的认识,加深了对课本知识的理解,而且让我更进一步熟悉了相关的MATLAB软件的基本编程方法和使用方法。
在这次课程设计的过程中,从整体思路的构建到具体每一步的实现,过程并不是一帆风顺的,通过复习课本知识以及查阅图书馆有关资料确定了整体思路,然后在稿纸上演算确定校正网络参数,最后运用MATLAB软件编程验证,作图。
在word编辑和运用MATLAB软件时遇到了一些问题,通过上网查询或者请教同学都得到了解决,因而设计的过程中也提高了我对word和matlab软件运用能力。
本次课程设计的核心之一是MATLAB软件的使用,通过软件编程,我对MATLAB的语言和应用有了更进一步的了解,特别熟悉了一些对自动控制很实用的命令,伯德图绘制函数bode()等。
利用MATLAB对控制系统进行频域分析,大大简化了计算和绘图步骤,是一款很实用的软件,今后利用课余时间也可以拓展一下自己在MATLAB软件编程的能力。
这次课程设计完成后,我体会到了学习自动控制原理,不仅要掌握书本上的内容,还要灵活思考,善于变换,在提出问题、分析问题、解决问题的过程中不断提高自己分析和解决实际问题的能力。
要把理论知识与实践相结合起来,从而提高自己的实际动手能力和独立思考的能力。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 课程设计