自动实验2预习报告leeee.docx

1、自动实验2预习报告leeee2010-2011 学年第 1 学期 院 别: 控制工程学院 课程名称: 自动控制原理 实验名称: PID控制特性的实验研究 实验教室: 6111 指导教师: 周蓉 小组成员（姓名，学号）: 程凤 2008071048 李同学2008071059 顾新武 2008071060 刘虎 2008071061 实验日期： 2010 年 月 日 评 分： 一、实验目的1、学习并掌握利用MATLAB编程平台进行控制系统复数域和时域仿真的方法；2、通过仿真实验，学习并掌握应用根轨迹分析系统性能及根据系统性能选择系统参数的方法；3、通过仿真实验研究，总结PID控制规律及参数变化对

2、系统性能影响的规律。二、实验任务及要求（一）实验任务针对如图所示系统，设计实验及仿真程序，研究在控制器分别采用比例（P）、比例积分（PI）、比例微分（PD）及比例积分微分（PID）控制规律和控制器参数（Kp、KI、KD）不同取值时，控制系统根轨迹和阶跃响应的变化，总结PID控制规律及参数变化对系统性能、系统根轨迹、系统阶跃响应影响的规律。具体实验内容如下：1、比例（P）控制，设计参数Kp使得系统处于过阻尼、临界阻尼、欠阻尼三种状态，并在根轨迹图上选择三种阻尼情况的Kp值，同时绘制对应的阶跃响应曲线，确定三种情况下系统性能指标随参数Kp的变化情况。总结比例（P）控制的规律。2、比例积分（PI）控

3、制，设计参数Kp、KI使得由控制器引入的开环零点分别处于：1）被控对象两个极点的左侧；2）被控对象两个极点之间；3）被控对象两个极点的右侧（不进入右半平面）。分别绘制三种情况下的根轨迹图，在根轨迹图上确定主导极点及控制器的相应参数；通过绘制对应的系统阶跃响应曲线，确定三种情况下系统性能指标随参数Kp和KI的变化情况。总结比例积分（PI）控制的规律。3、比例微分（PD）控制，设计参数Kp、KD使得由控制器引入的开环零点分别处于：1）被控对象两个极点的左侧；2）被控对象两个极点之间；3）被控对象两个极点的右侧（不进入右半平面）。分别绘制三种情况下的根轨迹图，在根轨迹图上确定主导极点及控制器的相应参

4、数；通过绘制对应的系统阶跃响应曲线，确定三种情况下系统性能指标随参数Kp和KD的变化情况。总结比例积分（PD）控制的规律。 4、比例积分微分（PID）控制，设计参数Kp、KI、KD使得由控制器引入的两个开环零点分别处于：1）实轴上：固定一个开环零点在被控对象两个开环极点的左侧，使另一个开环零点在被控对象的两个极点的左侧、之间、右侧（不进入右半平面）移动。分别绘制三种情况下的根轨迹图，在根轨迹图上确定主导极点及控制器的相应参数；通过绘制对应的系统阶跃响应曲线，确定三种情况下系统性能指标随参数Kp、KI和KD的变化情况。2）复平面上：分别固定两个共轭开环零点的实部（或虚部），让虚部（或实部）处于三

5、个不同位置，绘制根轨迹图并观察其变化；在根轨迹图上选择主导极点，确定相应的控制器参数；通过绘制对应的系统阶跃响应曲线，确定六种情况下系统性能指标随参数Kp、KI和KD的变化情况。综合以上两类结果，总结比例积分微分（PID）控制的规律。（二）实验要求1、合理选择P、PI、PD、PID控制器参数，使开环系统极零点分布满足实验内容中的要求。通过绘图展示不同控制规律和参数变化对系统性能的影响。根轨迹图可以单独绘制，按照不同控制规律、不同参数将阶跃响应绘制于同一幅面中。2、通过根轨迹图确定主导极点及参数值，根据阶跃响应曲线确定系统性能指标并列表进行比较，总结控制器控制规律及参数变化对系统特性、系统根轨迹

6、影响的规律。3、总结在一定控制系统性能指标要求下，根据系统根轨迹图和阶跃响应选择PID控制规律和参数的规则。4、全部采用MATLAB平台编程完成。三、实验方案设计（含实验参数选择、控制器选择、仿真程序等）1、比例（P）控制，设计参数Kp使得系统处于过阻尼、临界阻尼、欠阻尼三种状态，并在根轨迹图上选择三种阻尼情况的Kp值，同时绘制对应的阶跃响应曲线，确定三种情况下系统性能指标随参数Kp的变化情况。总结比例（P）控制的规律。实验系统模型： 控制器结构应为:绘制根轨迹图，选择实验参数：（1）、过阻尼实验参数选择：序号（组）y1y2y3y4y5y6Kp0.511.522.53实验代码：t=0:0.01

7、:3;Kp1=0.5;Kp2=1;Kp3=1.5;Kp4=2;Kp5=2.5;Kp6=3;num1=Kp1;num2=Kp2;num3=Kp3;num4=Kp4;num5=Kp5;num6=Kp6;den1=1 8 12+Kp1;den2=1 8 12+Kp2;den3=1 8 12+Kp3;den4=1 8 12+Kp4;den5=1 8 12+Kp5;den6=1 8 12+Kp6;step(num1,den1,t);hold on;step(num2,den2,t);step(num3,den3,t);step(num4,den4,t);step(num5,den5,t);step(nu

8、m6,den6,t);grid;gtext(Kp1=0.5);gtext(Kp2=1);gtext(Kp3=1.5);gtext(Kp4=2);gtext(Kp5=2.5);gtext(Kp6=3);系统阶跃响应曲线：系统性能指标随参数Kp的变化情况:序号（组）y1y2y3y4y5y6Tr (s)0.9090.9480.991.0841.081.14Tp (s)333333Ts (s)1.61.681.761.851.952.05P.O.000000（2）、临界阻尼实验参数选择：序号（组）yKp4实验代码：t=0:0.01:10;Kp=4;num=Kp;den=1 8 Kp+12;y,x,t=

9、step(num,den,t);plot(t,y);legend(y);xlabel(t);ylabel(y(t);grid;title(Kp=4)系统阶跃响应曲线：系统性能指标:Tr (s)Tp (s)Ts (s)P.O.ess0.843.51.4600（3）、欠阻尼实验参数选择：序号（组）y1y2y3y4y5y6Kp51015202530实验代码：t=0:0.01:3;Kp1=5;Kp2=10;Kp3=15;Kp4=20;Kp5=25;Kp6=30;num1=Kp1;num2=Kp2;num3=Kp3;num4=Kp4;num5=Kp5;num6=Kp6;den1=1 8 12+Kp1;d

10、en2=1 8 12+Kp2;den3=1 8 12+Kp3;den4=1 8 12+Kp4;den5=1 8 12+Kp5;den6=1 8 12+Kp6;step(num1,den1,t);hold on;step(num2,den2,t);step(num3,den3,t);step(num4,den4,t);step(num5,den5,t);step(num6,den6,t);grid;gtext(Kp1=5);gtext(Kp2=10);gtext(Kp3=15);gtext(Kp4=20);gtext(Kp5=25);gtext(Kp6=30);系统阶跃响应曲线：g0=tf(1,

11、conv(1 2,1 6);figure(1);rlocus(g0);kp=2.68,4,31.5,78;hold onfor i=1:length(kp) gc=kp(i) g=g0*gc sys=feedback(g,1) figure(2) subplot(2,2,i) step(sys)end2、比例积分（PI）控制，设计参数Kp、KI使得由控制器引入的开环零点分别处于：（1）被控对象两个极点的左侧；（2）被控对象两个极点之间；（3）被控对象两个极点的右侧（不进入右半平面）。分别绘制三种情况下的根轨迹图，在根轨迹图上确定主导极点及控制器的相应参数；通过绘制对应的系统阶跃响应曲线，确定三

12、种情况下系统性能指标随参数Kp和KI的变化情况。总结比例积分（PI）控制的规律。实验系统模型： 控制器结构应为：（1）由控制器引入的开环零点分别处于被控对象两个极点的左侧；绘制根轨迹图，选择实验参数：Kp=1;Ki=10;p=1 Ki/Kp;q=1 8 12 0;sys=tf(p,q);rlocus(sys);grid;rlocfind(sys)实验参数选择(固定Ki/Kp=10,Ki随Kp变化取值，y3的Kp值为临界阻尼点的K值)：y1y2y3y4y5y6Kp0.40.50.551510Ki455.51050100实验代码：t=0:0.01:16;Kd1=2.128; Kp1= 0.152;

13、Ki1=6.08; Kd2=6.916; Kp2=0.494; Ki2=19.76;Kd3=14; Kp3=1; Ki3=40;Kd4=27.86; Kp4= 1.99; Ki4=79.6;num1=Kd1 Kp1 Ki1;num2=Kd2 Kp2 Ki2;num3=Kd3 Kp3 Ki3;num4=Kd4 Kp4 Ki4;den1=1 Kd1+8 12+Kp1 Ki1;den2=1 Kd2+8 12+Kp2 Ki2;den3=1 Kd3+8 12+Kp3 Ki3;den4=1 Kd4+8 12+Kp4 Ki4;step(num1,den1);grid;holdon;step(num2,den

14、2);grid;hold on;step(num3,den3);grid;hold on;step(num4,den4);grid;hold ontitle(Kp,Kd;Ki)=(2.18,0.152,6.08),(6.916,0.494,19.76),(14,1,40),(27.86,1.99,79.6);gtext(Kp=2.18,Kd=0.152,Ki=6.08);gtext(Kp=6.916,Kd=0.494,Ki=19.76);gtext(Kp=14,Kd=1,Ki=40)gtext(Kp=27.86,Kd=1.99,Ki=79.6)系统阶跃响应曲线：系统性能指标随参数Kp、Ki的变

15、化情况:序号（组）y1y2y3y4Tr (s)3.811.891.191.19Tp (s)6.6 9.152.772.49P.O.1.51 1921.817.8ess1.021.191.221.18（2）、由控制器引入的开环零点分别处于被控对象两个极点之间：绘制根轨迹图，选择实验参数：Kp=1;Ki=4;p=1 Ki/Kp;q=1 8 12 0;sys=tf(p,q);rlocus(sys);grid;rlocfind(sys)实验参数选择(固定Ki/Kp=4,Ki随Kp变化取值，y3的Kp值为临界阻尼点的K值)：y1y2y3y4y5y6Kp0.611.673510Ki2.446.681220

16、40实验代码：t=0:0.01:25;Kp1=0.6;Kp2=1;Kp3=1.67;Kp4=3;Kp5=5;Kp6=10;Ki1=2.4;Ki2=4;Ki3=6.68;Ki4=12;Ki5=20;Ki6=40;num1=Kp1 Ki1;num2=Kp2 Ki2;num3=Kp3 Ki3;num4=Kp4 Ki4;num5=Kp5 Ki5;num6=Kp6 Ki6;den1=1 8 12+Kp1 Ki1;den2=1 8 12+Kp2 Ki2;den3=1 8 12+Kp3 Ki3;den4=1 8 12+Kp4 Ki4;den5=1 8 12+Kp5 Ki5;den6=1 8 12+Kp6 K

17、i6;step(num1,den1,t);holdon;step(num2,den2,t);step(num3,den3,t);step(num4,den4,t);step(num5,den5,t);step(num6,den6,t);grid;gtext(Kp1=0.6 ki1=2.4);gtext(Kp2=1 ki2=4);gtext(Kp3=1.67 ki3=6.68);gtext(Kp4=3 ki=12);gtext(Kp5=5 ki5=20);gtext(Kp6=10 ki6=40);系统阶跃响应曲线：系统性能指标随参数Kp、Ki的变化情况:序号（组）y1y2y3y4y5y6Tr (

18、s)10.15.73.111.590.9660.542Tp (s)2525253.332.041.23Ts (s)18.210.35.383.093.66 2.93P.O. 0002.198.2618.9ess0001.021.081.19（3）、由控制器引入的开环零点分别处于被控对象两个极点的右侧（不进入右半平面）：绘制根轨迹图，选择实验参数：Kp=1;Ki=1;p=1 Ki/Kp;q=1 8 12 0;sys=tf(p,q);rlocus(sys);grid;rlocfind(sys)实验参数选择(固定Ki/Kp=1,Ki随Kp变化取值，y3的Kp值为临界阻尼点的K值)：y1y2y3y4y

19、5y6Kp455.37101520Ki455.37101520实验代码：t=0:0.01:20;Kp1=4;Kp2=5;Kp3=5.37;Kp4=10;Kp5=15;Kp6=20;Ki1=4;Ki2=5;Ki3=5.37;Ki4=10;Ki5=15;Ki6=20;num1=Kp1 Ki1;num2=Kp2 Ki2;num3=Kp3 Ki3;num4=Kp4 Ki4;num5=Kp5 Ki5;num6=Kp6 Ki6;den1=1 8 12+Kp1 Ki1;den2=1 8 12+Kp2 Ki2;den3=1 8 12+Kp3 Ki3;den4=1 8 12+Kp4 Ki4;den5=1 8 1

20、2+Kp5 Ki5;den6=1 8 12+Kp6 Ki6;step(num1,den1,t);hold on;step(num2,den2,t);step(num3,den3,t);step(num4,den4,t);step(num5,den5,t);step(num6,den6,t); grid;gtext(Kp1=4 ki1=4);gtext(Kp2=5 ki2=5);gtext(Kp3=5.37 ki3=5.37);gtext(Kp4=10 ki=10);gtext(Kp5=15 ki5=15);gtext(Kp6=20 ki6=20);系统阶跃响应曲线：系统性能指标随参数Kp、Ki

21、的变化情况:序号（组）y1y2y3y4y5y6Tr (s)7.015.665.292.91.911.34Tp (s)202020202020Ts (s)12.910.59.915.944.43.6P.O.000000ess000000g0=tf(1,conv(1 2,1 6);gc=tf(1 4,1 0);figure(1);rlocus(g0*gc);kp=1.47,15.4,38.6;hold on;for i=1:length(kp) g=g0*gc*kp(i) sys=feedback(g,1) figure(2) subplot(2,2,i) step(sys)endhold off

22、gtext(kp=1.47);gtext(kp=15.4);gtext(kp=38.6)3、比例微分（PD）控制，设计参数Kp、KD使得由控制器引入的开环零点分别处于：（1）被控对象两个极点的左侧；（2）被控对象两个极点之间；（3）被控对象两个极点的右侧（不进入右半平面）。分别绘制三种情况下的根轨迹图，在根轨迹图上确定主导极点及控制器的相应参数；通过绘制对应的系统阶跃响应曲线，确定三种情况下系统性能指标随参数Kp和KD的变化情况。总结比例积分（PD）控制的规律。实验系统模型： 控制器结构应为：Gc(s)=Kp+Kd*s=Kd(s+Kp/Kd)（1）由控制器引入的开环零点分别处于被控对象两个极点

23、的左侧；绘制根轨迹图，选择实验参数：Kp=10;Kd=1;p=1 Kp/Kd;q=1 8 12;sys=tf(p,q);rlocus(sys);grid;rlocfind(sys)实验参数选择(固定Kp/Kd=10,Kp随Kd变化取值，y3的Kd值为临界阻尼点的K值)：y1y2y3y4y5y6Kp456.8101520Kd0.40.50.6811.52实验代码：t=0:0.01:3;Kp1=4;Kp2=5;Kp3=6.8;Kp4=10;Kp5=15;Kp6=20;Kd1=0.4;Kd2=0.5;Kd3=0.68;Kd4=1;Kd5=1.5;Kd6=2;num1=Kd1 Kp1;num2=Kd2

24、 Kp2;num3=Kd3 Kp3;num4=Kd4 Kp4;num5=Kd5 Kp5;num6=Kd6 Kp6;den1=1 Kd1+8 12+Kp1;den2=1 Kd2+8 12+Kp2;den3=1 Kd3+8 12+Kp3;den4=1 Kd4+8 12+Kp4;den5=1 Kd5+8 12+Kp5;den6=1 Kd6+8 12+Kp6;step(num1,den1,t);hold on;step(num2,den2,t);step(num3,den3,t);step(num4,den4,t);step(num5,den5,t);step(num6,den6,t);grid;gt

25、ext(Kp1=4 ki1=0.4);gtext(Kp2=5 ki2=0.5);gtext(Kp3=6.8 ki3=0.68);gtext(Kp4=10 ki=1);gtext(Kp5=15 ki5=1.5);gtext(Kp6=20 ki6=2);系统阶跃响应曲线：系统性能指标随参数Kp、Kd的变化情况:序号（组）y1y2y3y4y5y6Tr (s)0.4270.5030.6180.7260.8070.859Tp (s)33331.311Ts (s)1.481.381.221.010.7940.625P.O.0000.002880.1120.375ess0.3620.3620.3620.45

26、50.5560.627（2）由控制器引入的开环零点分别处于被控对象两个极点之间；绘制根轨迹图，选择实验参数：Kp=4;Kd=1;p=1 Kp/Kd;q=1 8 12;sys=tf(p,q);rlocus(sys);grid;rlocfind(sys)实验参数选择：y1y2y3y4y5y6Kp41220283644Kd1357911实验代码：t=0:0.01:2;Kp1=4;Kp2=12;Kp3=20;Kp4=28;Kp5=36;Kp6=44;Kd1=1;Kd2=3;Kd3=5;Kd4=7;Kd5=9;Kd6=11;num1=Kd1 Kp1;num2=Kd2 Kp2;num3=Kd3 Kp3;n

27、um4=Kd4 Kp4;num5=Kd5 Kp5;num6=Kd6 Kp6;den1=1 Kd1+8 12+Kp1;den2=1 Kd2+8 12+Kp2;den3=1 Kd3+8 12+Kp3;den4=1 Kd4+8 12+Kp4;den5=1 Kd5+8 12+Kp5;den6=1 Kd6+8 12+Kp6;step(num1,den1,t);hold on;step(num2,den2,t);step(num3,den3,t);step(num4,den4,t);step(num5,den5,t);step(num6,den6,t);grid;gtext(Kp1=4 kd1=0.41);gtext(Kp2=12 kd2=3);gtext(Kp3=20 kd3=5);gtext(Kp4=28 kd4=7);gtext(Kp5=36 kd5=9);gtext(Kp6=44 kd6=11);系统阶跃响应曲线：系统性能指标随参数Kp、Kd的变化情况:序号（组）y1y2y3y4y5y6Tr (s)0.7330.4830.3540.2760.225