自控课程设计1

  • 格式:doc
  • 大小:234.50 KB
  • 文档页数:12

课程设计报告( 2006-- 2007年度第一学期)名称:自动控制理论课程设计题目:院系:自动化系班级:自动化046学号:200403020128学生姓名:张金营指导教师:于希宁孙建平刘鑫屏孙海蓉魏乐黄宇设计周数:1周成绩:日期:年月日一、目的与要求本课程为《自动控制理论A》的课程设计,是课堂的深化。

设置《自动控制理论A》课程设计的目的是使MATLAB成为自动化专业学生的基本技能,熟悉MATLAB这一解决具体工程问题的标准软件,能熟练地应用MATLAB软件解决控制理论中的复杂和工程实际问题,并给以后的模糊控制理论、最优控制理论和多变量控制理论等的学习奠定良好的基础。

作为自动化专业的学生很有必要学会应用这一强大的工具,并掌握利用MATLAB对控制理论内容进行分析和研究的技能,以达到加深对课堂上所讲内容理解的目的。

通过使用这一软件工具把学生从繁琐枯燥的计算负担中解脱出来,而把更多的精力用到思考本质问题和研究解决实际生产问题上去。

通过此次计算机辅助设计,掌握以下基本技能:1.能用MATLAB软件分析复杂控制系统的性能。

2.能用MATLAB软件定性或定量设计控制系统以满足具体的性能指标要求。

3.能灵活应用MATLAB的CONTROL SYSTEM 工具箱和SIMULINK仿真软件,分析系统的性能。

二、设计正文2.1.数学模型题目设计:2.1.1 某一以微分方程描述系统的传递函数,其微分方程描述如下:(3)(2)(2)(1)111048y y y u u u++=++求其传递函数。

解:建立模型的MA TLAB代码如下:num=[1,4,8];%分子多项式den=[1,11,11,10];%分母多项式G=tf(num,den)%建立传递函数模型程序运行结果如下:Transfer function:s^2 + 4 s + 8------------------------s^3 + 11 s^2 + 11 s + 102.1.2 已知某系统的传递函数2()3248111110ss sGs s s++=+++,求其分子分母多项式并绘制零极点图。

解:MATLAB程序代码如下:num=[1,4,8]; %分子多项式den=[1,11,11,10]; %分母多项式G=tf(num,den) %建立传递函数模型[tt,ff]=tfdata(G,'v')%提取传递函数的分子分母多项式[z,p,k]=tf2zp(num,den)%提取传递函数的零极点与增益pzmap(G)%绘制零极点图grid on %给图形加网格程序运行结果如下:Transfer function:s^2 + 4 s + 8------------------------s^3 + 11 s^2 + 11 s + 10tt =0 1 4 8ff =1 11 11 10z =-2.0000 + 2.0000i-2.0000 - 2.0000ip =-10.0000-0.5000 + 0.8660i-0.5000 - 0.8660ik =1输出图形如下(图一):图一2.2. 时域分析题目设计:2.2.1 稳定性分析:典型二阶输入如下:2()222nsn nGs sωξωω=++式中,nω是自然振荡频率(无阻尼振荡频率),ξ是阻尼系数。

要求绘制当ξ=0.5,nω分别为2、4、6、8、10、12时系统的单位阶跃响应。

MATLAB程序如下:syms sw=[2:2:12];kosai=0.5;figure(1)hold onfor wn=wnum=wn^2;den=[1,2*kosai*wn,wn^2];step(num,den);endhold offgrid on;title('单位阶跃响应曲线')xlabel('时间')ylabel('振幅')程序运行结果如下:图二2.2.2 稳态性能分析已知下图示的单位负反馈控制系统的传递函数为)25(10)s s G s s +=+(,试计算当输入为单位阶跃信号时系统的稳态误差。

解:Simulink 模型如图所示:信号源选定“step ”,连好模型进行仿真,仿真结束后,双击示波器,输出波形如图所示,它是输入为单位阶跃时系统的输出误差:误差波形图2.3. 运用MA TLAB 对根轨迹进行分析2.3.1设一系统的开环传递函数为)()(1)(0.5)s s kG H s s s =++(,绘制该系统根轨迹图并对系统判稳。

解:MATLAB 程序为:num=[1];den=conv([1 1 0],[0.5 1]); G=tf(num,den); rlocus(G)[k,ploles]=rlocfind(G)程序执行后,可在图形窗口根轨迹图中显示十字行光标,当十字行光标指向根轨迹与纵坐标交点时,对应的开环增益与极点为: k =2.8528ploles =-2.9728 -0.0136 + 1.3853i -0.0136 - 1.3853i 系统根轨迹图为:当系统参数()03K →变动时,根轨迹均在S 平面纵坐标的左侧,对应的系统闭环是稳定的。

一旦根轨迹穿越纵坐标到达其右侧,对应的K>3,那么系统闭环就不稳定。

当在根轨迹实轴上的区段时,对应着系统闭环阶跃响应无超调,系统闭环稳定。

当不在根轨迹实轴上的区段时,系统闭环特征方程出现了共轭复根,意味着系统闭环阶跃响应有超调但系统闭环还是稳定的。

当系统的K=1.2、3、4时,绘制系统闭环单位阶跃响应曲线。

MATLAB 程序如下:t=0:0.1:30; num1=[1];den=conv(conv([1 0],[1 1]),[0.5 1]); num=1.2*num1; G1=tf(num,den); sys=feedback(G1,1); figure(1); step(sys,t); hold on num=3*num1; G2=tf(num,den); sys=feedback(G2,1); figure(2); step(sys,t); hold on num=4*num1' G3=tf(num,den); sys=feedback(G3,1); figure(3); step(sys,t)绘制图形如下,由图可见,当K=1.2时,系统单位阶跃响应曲线为衰减振荡(如下图A );当K=3时,系统的单位阶跃响应曲线为等幅振荡(如下图B ),对应着系统临界稳定;当K=4时,系统的单位阶跃响应曲线为发散的振荡(如下图C ),对应系统不稳定。

这完全验证了以上的结论。

图A 图B 图C2.3.2 绘制参数a 的根轨迹,使用MATLAB 画系统的根轨迹,该系统的开环传递函数如下:)()5()(1)(3)(12)s s s a G H s s s +=+++(,其中212a ≤≤解:MATLAB 程序如下:den=conv(conv([1 1],[1,3]),[1 12]);k=5;% 定义数组存储结果clpoles=[];param=[];for a=2:10num=[0 0 k k*a];clpoly=num+den;%计算闭环极点clp=roots(clpoly);clpoles=[clpoles;clp'];param=[param;a];end%打印a和极点表格disp([param,clpoles])plot(clpoles,'*')axis equal;%调整绘制区域axis([-4 0 -2 2])程序运行结果如下:2.0000 -11.4658 -3.3291 -1.20513.0000 -11.5249 -3.0000 -1.47514.0000 -11.5826 -2.4174 -2.00005.0000 -11.6388 -2.1806 - 0.6972i -2.1806 + 0.6972i6.0000 -11.6938 -2.1531 - 1.0041i -2.1531 + 1.0041i7.0000 -11.7475 -2.1262 - 1.2341i -2.1262 + 1.2341i8.0000 -11.8001 -2.1000 - 1.4251i -2.1000 + 1.4251i9.0000 -11.8516 -2.0742 - 1.5913i -2.0742 + 1.5913i10.0000 -11.9020 -2.0490 - 1.7399i -2.0490 + 1.7399i根轨迹如图所示:2.4用MA TLAB 对频域进行分析 2.4.1 Nyquist 曲线的绘制已知一个典型的一阶环节传递函数:()531s G s =+,试绘制该环节的Nyquist 图。

解:MATLAB 程序如下:num=5; den=[3 1] G=tf(num,den) nyquist(G) grid响应曲线如图所示:2.4.2 Bode 图行的绘制:MATLAB 代码如下: w=[0,logspace(-2,2,200)]wn=0.7kosai=[0.1,0.4,1.0,1.6,2.0] for j=1:5sys=tf([wn*wn],[1,2*kosai(j)*wn,wn*wn]) bode(sys,w) hold on endgtext('kosai=0.1') gtext('kosai=0.4') gtext('kosai=1.0') gtext('kosai=1.6') gtext('kosai=2.0')Bode 图形如下:2.4.3某控制系统的开环传递函数为)12.0()1(32)()(++=s s s H G S S ,试求系统的频域性能指标c w 、γ。

解:MATLAB 代码如下:num=3*[1 1];den=conv([1 0 0],[0.2 1]); GH=tf(num,den); margin(GH) 程序运行结果如下:Bode图系统性能指标如下:剪切频率:w=2.78rad/s 相角稳定裕度: =41.1 c2.5控制系统的校正以及校正前后的性能分析2.5.1 根轨迹法超前校正以及校正前后的性能分析2.5.2 频域法超前校正以及校正前后的性能分析2.5.3 频域法滞后校正以及校正前后的性能分析。