1、自 动 控 制 原 理课 程 设 计（理工类）课程名称： 自动控制原理 专业班级： 08 自动化（1）班 学生学号： 0804110601 学生姓名： 丁丽华 所属院部： 机电工程学院 指导教师： 陈丽换 2009 2010 学年 第二学期金陵科技学院教务处制金 陵 科 技 学 院自动控制原理课程设计任务书课程序号 32 课程编号 04184500实践序号 10 设计名称 自动控制原理课程设计适用年级、专业 08 自动化 时间 1 周一、 设计目的：1、 了解控制系统设计的一般方法、步骤。2、 掌握对系统进行稳定性分析、稳态误差分析以及动态特性分析的方法。3、 掌握利用 MATLAB 对控制理
2、论内容进行分析和研究的技能。4、 提高分析问题解决问题的能力。二、 设计内容与要求：设计内容：1、阅读有关资料。2、对系统进行稳定性分析、稳态误差分析以及动态特性分析。3、绘制根轨迹图、Bode 图、Nyquist 图。4、设计校正系统，满足工作要求。设计条件： 1、已知单位负反馈系统被控制对象的传递函数为 0KG(S)1)(2设计要求：1、能用 MATLAB 解复杂的自动控制理论题目。2、能用 MATLAB 设计控制系统以满足具体的性能指标。3、能灵活应用 MATLAB 的 CONTROL SYSTEM 工具箱和 SIMULINK仿 真软件，分析系统的性能。设计题目: ，试用频率法设计串联滞
3、后超前校正装置，使系统0KG(S)1)(2的相角裕量 ，静态速度误差系数 ,截止频率不低于451vK0s1.5rads设计步骤：1、静态速度误差系数 ，即当 S0 时， =10，解得1vK0s0G()S)(2K0=20s-1。即被控对象的开环传递函数：G(S)= 。)(122、滞后校正器的传递函数为：G C1(S)= TSb1根据题目要求，取校正后系统的截止频率 WC=1.5rad/s,先试取 b=0.105，编写求滞后校正器的传递函数的 MATLAB 的程序如下：wc=1.5;k0=20;n1=1;d1=conv(conv(1 0,1 1),1 2);b=0.105;T=1/(0.1*wc)
4、;B=b*T;Gc1=tf(B 1,T 1)将程序输入 MATLAB Command Window 后，并按回车，Command Window 出现如下代数式： 13.67s由式可知：b=0.105,T=63.33。3、 求超前校正器的传递函数，而已知串联有滞后校正器的传递函数为：G(S)Gc1(S)= )2(10S13.67s根据校正后系统的传递函数，编写求超前校正器的传递函数的 MATLAB 程序，其中调用了求超前校正器传递函数的函数 leadc()，leadc.m 保存在 matlab6.5work文件夹下，leadc.m 编制如下：function Gc=leadc(key,sope,
5、vars)% MATLAB FUNCTION PROGRAM leadc.m%if key=1gama=vars(1);gama1=gama+5;mag,phase,w=bode(sope);mu,pu=bode(sope,w);gam=gama1*pi/180;alpha=(1-sin(gam)/(1+sin(gam);adb=20*log10(mu);am=10*log10(alpha);wc=spline(adb,w,am);T=1/(wc*sqrt(alpha);alphat=alpha*T;Gc=tf(T 1,alphat 1);elseif key=2wc=vars(1);num=
6、sope.num1;den=sope.den1;na=polyval(num,j*wc);da=polyval(den,j*wc);g=na/da;g1=abs(g);h=20*log10(g1);a=10(h/10);wm=wc;T=1/(wm*(a)(1/2);alphat=a*T;Gc=tf(T 1,alphat 1);elseif key=3gama=vars(1);wc=vars(2);gama1=gama+5;num=sope.num1;den=sope.den1;ngv=polyval(num,j*wc);dgv=polyval(den,j*wc);g=ngv/dgv;theta
7、g=angle(g);thetag_d=thetag*180/pi;mg=abs(g);gama_rad=gama1*pi/180;z=(1+mg*cos(gama_rad-thetag)/(-wc*mg*sin(gama_rad-thetag);p=(cos(gama_rad-thetag)+mg)/(wc*sin(gama_rad-thetag);nc=z,1;dc=p,1;Gc=tf(nc,dc);End在 Command Window 中编写下列程序：n1=conv(0 20,6.667 1);d1=conv(1 0,1 1);d2=conv(1 2,63.33 1);d3=conv(
8、d1,d2);sope=tf(n1,d3);wc=1.5;Gc=leadc(2,sope,wc)写完后，按回车，出现如下代码：，即超前传递函数为 Gc2(S)= = ,可得12087.3s TSa112087.3sa=10.21,T=0.2087。故校正后的开环系统总传递函数为：G(S) Gc1(S) Gc2(S)= )2(10S13.67s12087.3s验证校正后的闭环系统的性能指标,画出 bode 图：由图可知：剪切频率为 Wc=1.54rad/s，相角裕量为 = 45，符合设计要求。4、 用 MATLAB 求校正前后系统的特征根： 20先写出校正前系统单位负反馈传递函数：（S）= S3
9、+3*S2+2S+20故在 Commmand Window 中写入：p=1 3 2 20;roots(p) 回车后得到：ans =-3.8371 0.4186 + 2.2443i0.4186 - 2.2443i由于校正前系统单位负反馈得特征方程有右半平面的根，故校正前的闭环系统不稳定。写出校正后系统单位负反馈传递函数：11.5290S2+7.1480S+0.8128（S）= S5+7.8058*S4+16.4931*S3+21.3676*S2+7.2994*S+ 0.8128 故在 Commmand Window 中写入 p=1.00007.8058 16.4931 21.3676 7.299
10、4 0.8128;roots(p)回车后得到：p =1.0000 7.8058 16.4931 21.3676 7.2994 0.8128ans =-5.4570 -0.9582 + 1.3694i-0.9582 - 1.3694i-0.2163 + 0.0810i-0.2163 - 0.0810i 由于校正后系统单位负反馈得特征方程没有右半平面的根，故校正后的闭环系统稳定。5、 用 MATLAB 作出系统校正前和校正后的单位脉冲响应曲线，单位阶跃响应曲线，单位斜坡响应曲线，并求出系统校正前和校正后的动态性能指标 %、tr、tp、ts 以及稳态误差的值： 校正前的单位阶跃响应：编写的程序如下：
11、n1=20;d1=conv(conv(1 0,1 1),1 2);s1=tf(n1,d1);sope=s1;sys=feedback(sope,1);step(sys)y,t=step(sys);得到的单位阶跃响应曲线图： 校正后的单位阶跃响应：编写的程序如下：n1=20;d1=conv(conv(1 0,1 1),1 2);s1=tf(n1,d1);s2=tf(6.667 1,63.33 1);s3=tf(2.13 1,0.2087 1);sope=s1*s2*s3;sys=feedback(sope,1);step(sys)y,t=step(sys);得到的单位阶跃响应曲线图： 校正前的单
12、位冲击响应：编写的程序如下：n1=20;d1=conv(conv(1 0,1 1),1 2);s1=tf(n1,d1);sope=s1;sys=feedback(sope,1);impulse(sys)y,t=impulse(sys);得到的单位冲击响应曲线图： 校正后的单位冲击响应：编写的程序如下：n1=20;d1=conv(conv(1 0,1 1),1 2);s1=tf(n1,d1);s2=tf(6.667 1,63.33 1);s3=tf(2.13 1,0.2087 1);sope=s1*s2*s3;sys=feedback(sope,1);impulse(sys)y,t=impuls
13、e(sys);得到的单位冲击响应曲线图： 校正前的单位斜坡响应：在 Simulink 窗口里菜单方式下的单位斜坡响应的动态结构图如下：在 Command Window 中输入 plot(tout,dimout),回车后得到的波形图为： 校正后的单位斜坡响应：在 Simulink 窗口里菜单方式下的单位斜坡响应的动态结构图如下：在 Command Window 中输入 plot(tout,bimout),回车后得到的波形图为：由于校正前的系统为不稳定系统，故不讨论其动态性能指标；而校正后的系统是稳定的系统，其超调量 %和峰值时间可由下列图得到：超调量 %为 22%.上升时间 tr 为：0.806
14、s由下图可得到校正后系统的调节时间 ts（5%）:调节时间 ts（5%）=6.03s由下图可得到校正后系统的调节时间 ts（2%）：调节时间 ts(2%)=14.4s右下图可以得到校正后系统的稳态误差值：稳态误差 ess=1-1=0。单位脉冲、阶跃、斜坡响应曲线的关系是：单位脉冲响应的积分是单位阶跃响应曲线，单位阶跃响应的积分是单位斜坡响应。6、 绘制系统校正前与校正后的根轨迹，并求其分离点、汇合点及虚轴交点的坐标和相应的增益 K*值，得出系统稳定时增益 K*的变化范围：求校正前的根轨迹的程序如下：num=1;den=1 3 2 0;rlocus(num,den,k)校正前的根轨迹的分离点和汇合点由下图得：分离点 d=-0.423增益 K*=0.385校正前的根轨迹虚轴的交点由下图得：与虚轴的交点 w=+1.41、-1.41，增益 K*=6.00校正前系统稳定时增益 K*的变化范围是 K*0,故系统稳定性能很好。并且校正后的系统的相角裕量 、静态速度误差系数 Kv、剪切频率 Wc 均符合题目设计的要求。故系统校正到此全部完成。小组成员信息：班级：08 级自动化（1）丁丽华 0804110601王凯 0804110603王建英 0804110604尹燕彬 0804110602自动控制原理 程 鹏主编基于 MATLAB7.x 的系统分析与设计控制系统 （第二版）楼顺天 主编
