实验一控制新版系统的稳定性分析.docx

试验一 控制系统稳定性分 班级：光伏2班 姓名：王永强 学号：120XX09067 试验一 控制系统稳定性分析 一、试验目标 1、研究高阶系统稳定性，验证稳定判据正确性； 2、了解系统增益改变对系统稳定性影响； 3、观察系统结构和稳态误差之间关系。 二、试验任务 1、稳定性分析 欲判定系统稳定性，只要求出系统闭环极点即可，而系统闭环极点就是闭环传输函数分母多项式根，能够利用MATLAB中tf2zp函数求出系统零极点，或利用root函数求分母多项式根来确定系统闭环极点，从而判定系统稳定性。 （1）已知单位负反馈控制系统开环传输函数为，用MATLAB编写程序来判定闭环系统稳定性，并绘制闭环系统零极点图。 在MATLAB命令窗口写入程序代码以下： z=-2.5 p=[0,-0.5,-0.7,-3] k=1 Go=zpk(z,p,k) Gc=feedback(Go,1) Gctf=tf(Gc) dc=Gctf.den dens=ploy2str(dc{1},'s') 运行结果以下： Gctf = s + 2.5 --------------------------------------- s^4 + 4.2 s^3 + 3.95 s^2 + 2.05 s + 2.5 Continuous-time transfer function. dens是系统特征多项式，接着输入以下MATLAB程序代码： den=[1,4.2,3.95,1.25,0.5] p=roots(den) 运行结果以下： p = -3.0058 + 0.0000i -1.0000 + 0.0000i -0.0971 + 0.3961i -0.0971 - 0.3961i p为特征多项式dens根，即为系统闭环极点，全部闭环极点全部是负实部，所以闭环系统是稳定。 下面绘制系统零极点图，MATLAB程序代码以下： z=-2.5 p=[0,-0.5,-0.7,-3] k=1 Go=zpk(z,p,k) Gc=feedback(Go,1) Gctf=tf(Gc) [z,p,k]=zpkdata(Gctf,'v') pzmap(Gctf) Grid 运行结果以下： z = -2.5000 p = -3.0297 + 0.0000i -1.3319 + 0.0000i 0.0808 + 0.7829i 0.0808 - 0.7829i k =1 输出零极点分布图 （2）已知单位负反馈控制系统开环传输函数为，当取=1，10，100用MATLAB编写程序来判定闭环系统稳定性。 只要将（1）代码中k值变为1，10，100，即可得到系统闭环极点，从而判定系统稳定性，并讨论系统增益k改变对系统稳定性影响。 1.当K=1时 在MATLAB命令窗口写入程序代码以下： z=-2.5 p=[0,-0.5,-0.7,-3] k=1 Go=zpk(z,p,k) Gc=feedback(Go,1) Gctf=tf(Gc) dc=Gctf.den dens=ploy2str(dc{1},'s') 运行结果以下： Gctf = s + 2.5 --------------------------------------- s^4 + 4.2 s^3 + 3.95 s^2 + 2.05 s + 2.5 Continuous-time transfer function. 接着输入以下MATLAB程序代码： den=[1,4.2,3.95,1.25,0.5] p=roots(den) 运行结果以下： p = -3.0058 + 0.0000i -1.0000 + 0.0000i -0.0971 + 0.3961i -0.0971 - 0.3961i 下面绘制系统零极点图，MATLAB程序代码以下： z=-2.5 p=[0,-0.5,-0.7,-3] k=1 Go=zpk(z,p,k) Gc=feedback(Go,1) Gctf=tf(Gc) [z,p,k]=zpkdata(Gctf,'v') pzmap(Gctf) Grid 运行结果以下： z = -2.5000 p = -3.0297 + 0.0000i -1.3319 +