计控实验报告.docx

1、计控实验报告 数字PID 控制器的设计 直流闭环调速实验一、实验目的：1. 理解晶闸管直流单闭环调速系统的数学模型和工作原理；2. 掌握PID 控制器参数对控制系统性能的影响；3. 能够运用MATLAB/Simulink 软件对控制系统进行正确建模并对模块进行正确的参数设置；4. 掌握计算机控制仿真结果的分析方法。二、实验工具：MATLAB 软件（7.0 版本）。三、实验内容：已知晶闸管直流单闭环调速系统的转速控制器为PID 控制器，如图1所示。试运用MATLAB 软件对调速系统的P、I、D 控制作用进行分析。四、实验步骤：（一）模拟PID 控制作用分析：运用MATLAB 软件对调速系统的P、

2、I、D 控制作用进行分析。（1）比例控制作用分析为分析纯比例控制的作用，考察当Td = 0, Ti = , Kp =15时对系统阶跃响应的影响。1) MATLAB实验程序如下：G1=tf(1,0.017 1);G2=tf(1,0.075 0);G12=feedback(G1*G2,1);G3=tf(44,0.00167 1);G4=tf(1,0.1925);G=G12*G3*G4;Kp=1:1:5;for i=1:length(Kp)Gc=feedback(Kp(i)*G,0.01178);step(Gc),hold onendaxis(0 0.2 0 130); legend(kp=1,kp

3、=2,kp=3,kp=4,kp=5);2) 实验结果如图2图2 P 控制阶跃响应曲线3) 结果分析观察比例控制作用的阶跃响应曲线知，在比例控制作用下，系统的响应速度与比例系数成正比，比例系数越大，比例作用越强，响应速度越快；同时，可以看出比例系数越大系统的稳态误差越小；但是，比例作用会加大超调量，比例系数越大，系统的超调量越大。（2）积分控制作用分析保持Kp =1不变，考察Ti = 0.03 0.07 时对系统阶跃响应的影响。1) MATLAB 程序如下：G1=tf(1,0.017 1);G2=tf(1,0.075 0);G12=feedback(G1*G2,1);G3=tf(44,0.001

4、67 1);G4=tf(1,0.1925);G=G12*G3*G4;Kp=1;Ti=0.03:0.01:0.07;for i=1:length(Ti)Gc=tf(Kp*Ti(i) 1,Ti(i) 0);Gcc=feedback(G*Gc,0.01178);step(Gcc),hold onendlegend(Ti=0.03,Ti=0.04,Ti=0.05,Ti=0.06,Ti=0.07)2) 实验结果如图3图 3 I 控制阶跃响应曲线3) 结果分析 比例积分控制作用下的单位阶跃响应可以看出，在开环比例系数不变的情况下，积分作用提高了系统的型别，积分系数变大，系统的响应速度变慢。（3）微分控制作

5、用分析为分析微分控制的作用，保持Kp = 0.01 ，Ti = 0.01, 考察Td = 12 84 时对系统阶跃响应的影响。1) Matlab程序G1=tf(1,0.017 1);G2=tf(1,0.075 0);G12=feedback(G1*G2,1);G3=tf(44,0.00167 1);G4=tf(1,0.1925);G=G12*G3*G4;Kp=0.01;Ti=0.01;Td=12:36:84;for i=1:length(Td)Gc=tf(Kp*Ti*Td(i) Ti 1,Ti 0)Gcc=feedback(G*Gc,0.01178) ; step(Gcc),hold onen

6、dlegend(Td=12,Td=48,Td=84);2) 仿真结果如图4图4 D 控制阶跃响应曲线3) 结果分析又系统单位阶跃响应曲线可知，在微分作用下，微分系数越小，系统的响应速度越快，同时系统的振荡减小，更有利于系统的动态性能的改善。（二）数字PID 控制作用分析：仿照上述过程，独立完成PID 离散化仿真程序编写及结果分析。（1）比例控制作用下，取采样时间0.001 秒。ts=0.001;G1=tf(1,0.017 1);G2=tf(1,0.075 0);G12=feedback(G1*G2,1);G3=tf(44,0.00167 1);G4=tf(1,0.1925);G=G12*G3*

7、G4;Kp=1:1:5;for i=1:length(Kp)Gc=feedback(Kp(i)*G,0.01178);dsys=c2d(Gc,ts,z);step(dsys),hold onendaxis(0 0.2 0 130); legend(kp=1,kp=2,kp=3,kp=4,kp=5);仿真结果如图5图5 数字P 控制阶跃响应曲线（2）比例积分控制作用下，取采样时间0.001 秒。ts=0.001;G1=tf(1,0.017 1);G2=tf(1,0.075 0);G12=feedback(G1*G2,1);G3=tf(44,0.00167 1);G4=tf(1,0.1925);G

8、=G12*G3*G4;Kp=1;Ti=0.03:0.01:0.07;for i=1:length(Ti)Gc=tf(Kp*Ti(i) 1,Ti(i) 0);Gcc=feedback(G*Gc,0.01178);dsys=c2d(Gcc,ts,z);step(Gcc),hold onendlegend(Ti=0.03,Ti=0.04,Ti=0.05,Ti=0.06,Ti=0.07)仿真结果如图6图6 数字I 控制阶跃响应曲线（3）比例积分微分控制作用下，取采样时间0.05 秒。clcclearfigure(1);grid onhold oncolor=r,g,b,y,mts=0.001;G1=t

9、f(1,0.017 1);G2=tf(1,0.075 0);G12=feedback(G1*G2,1);G3=tf(44,0.00167 1);G4=tf(1,0.1925);G=G12*G3*G4;Ti=0.01;Td=0.084;kp1=1:1:5for i=1:length(kp1)kp=kp1(i);dsys=c2d(G,ts,z);num,den=tfdata(dsys,v);u1=0.0;u2=0.0;u3=0.0;y1=0;y2=0;y3=0;x=0 0 0;error1=0;for k=1:1:400 time(k)=k*ts; rin(k)=1; u(k)=kp*x(1);%

10、+ts/Ti*kp*x(3)+Td*kp*x(2)/ts; yout(k)=-den(2)*y1-den(3)*y2-den(4)*y3+num(2)*u1+num(3)*u2+num(4)*u3; error(k)=rin(k)-yout(k)*0.01178; u3=u2;u2=u1;u1=u(k); y3=y2;y2=y1;y1=yout(k); x(1)=error(k); x(2)=(error(k)-error1); x(3)=x(3)+error(k); error1=error(k);endplot(time,yout,color(i);endlegend(kp=1,kp=2,

11、kp=3,kp=4,kp=5); 仿真结果如图7图7 数字PID 控制作用仿真结果图结果分析：数字PID与模拟PID的性能相同。在PID控制作用下，在微分系数和积分系数不变的情况下，比例系数越大，系统的响应速度越快，同时系统稳态误差越小，但可以看出，比例系数增大会使系统的振荡加剧。实验三 复杂规律数字控制器的设计-DALIN 算法用于温度控制实验一、实验目的：1理解具有纯滞后的温度控制系统的数学模型和工作特点；2理解DALIN （大林）算法工作原理并掌握控制器参数对系统性能的影响；3能够运用MATLAB/Simulink 软件对控制系统进行正确建模并对模块进行正确的参数设置；4掌握计算机控制仿

12、真结果的分析方法。二、实验工具：MATLAB 软件（7.0版本）。三、实验内容：1实验原理图1 单回路温度控制系统图1 中：D(z)为数字调节器传递函数；A/D环节近似为一采样开关；H0(s)为D/A 环节传递函数，采用零阶保持器；G0(s)为温度被控对象函数，次环节可近似认为是一个纯滞后一阶惯性环节，所以系统是一个具有零阶保持器的单变量调节系统。根据大林算法的控制原理，将期望的闭环响应设计成一阶惯性加纯迟延，然后反过来得到能满足这种闭环响应的控制器。设图1所示单回路温度控制系统中，G0(s)为被控对象，D(z)为大林控制器，闭环控制系统传递函数为则有：由此我们可以知道，如果能事先设定系统的闭

13、环相应，则可得控制D(z)。2实验设计假定被控对象为：采样时间为0.5s ，期望的系统闭环传递函数为：，其中0.15s 为校正后闭环系统的时间常数。为便于比较，运用MATLAB软件使用大林算法及PID（KP=1.0,KI=0.50, KD=0.10）算法分别对一阶惯性纯滞后温度系统进行控制，并分析两种算法的优缺点。四、实验步骤：运用MATLAB 软件使用大林算法及PID算法分别对温度控制系统进行控制，并且分析各自的优缺点。MATLAB 程序设计具体步骤如下1被控对象离散化： num=1 den=0.4 1 G=tf(num,den)%被控对象传递函数 G.iodelay=0.76;dsys1=

14、c2d(G,ts,z);%广义被控对象脉冲传递函数2期望闭环系统离散化：2期望闭环系统离散化：%整个闭环离散化num1,den1=tfdata(dsys1,v);%期望闭环传递函数sys2=tf(1,0.15 1,inputdelay,0.76)dsys2=c2d(sys2,ts,z);%期望闭环脉冲传递函数 3设计DALIN 控制算法：%由公式，设计大林控制器，得num1,den1=tfdata(dsys1,v);%期望闭环传递函数sys2=tf(1,0.15 1,inputdelay,0.76)dsys2=c2d(sys2,ts,z);%期望闭环脉冲传递函数dsys=1/dsys1*dsy

15、s2/(1-dsys2);%大林算法控制器脉冲传递函数num,den=tfdata(dsys,v);得 Num=0.7981 -0.0909 -0.0454 0.0017 0 0 0 0 ；Den=0 0.4512 0.2301 -0.3782 -0.2712 -0.0335 0.0016 0；4与普通PID 控制进行比较： %使用大林算法，由den 和num 的值与公式,经计算推出u(k)的表达式 yout(k)=-den1(2)*y1+num1(2)*u2+num1(3)*u3;%被控对象输出 %返回大林算法中的参数 Num=0.7981 -0.0909 -0.0454 0.0017 0

16、0 0 0 ； Den=0 0.4512 0.2301 -0.3782 -0.2712 -0.0335 0.0016 0； %使用普通PID 算法，位置式数字PID 控制算法公式并KP=1.0,KI=0.50,KD=0.105、程序达林算法程序：ts=0.5;num=1den=0.4 1G=tf(num,den)%被控对象传递函数G.iodelay=0.76;dsys1=c2d(G,ts,z);%广义被控对象脉冲传递函数num1,den1=tfdata(dsys1,v);%期望闭环传递函数sys2=tf(1,0.15 1,inputdelay,0.76)dsys2=c2d(sys2,ts,z)

17、;%期望闭环脉冲传递函数dsys=1/dsys1*dsys2/(1-dsys2);%大林算法控制器脉冲传递函数num,den=tfdata(dsys,v);u1=0.0;u2=0.0;u3=0.0;u4=0;u5=0;y1=0;error1=0;error2=0;error3=0;for k=1:1:50 time(k)=k*ts; rin(k)=1.0; yout(k)=-den1(2)*y1+num1(2)*u2+num1(3)*u3;%被控对象输出 error(k)=rin(k)-yout(k); u(k)=(num(1)*error(k)+num(2)*error1+num(3)*er

18、ror2+num(4)*error3.%控制器输出 -den(3)*u1-den(4)*u2-den(5)*u3-den(6)*u4-den(7)*u5)/den(2);u5=u4;u4=u3;u3=u2;u2=u1;u1=u(k);y1=yout(k);error3=error2;error2=error1;error1=error(k);endsubplot(2,1,1);plot(time,rin,b,time,yout,r);xlabel(time(s);ylabel(rin,yout);title(DALIN method);grid on%#PID控制器输出#PID程序：u1=0.

19、0;u2=0.0;u3=0.0;u4=0;u5=0;y1=0;error1=0;error2=0;error3=0;ei=0;for k=1:1:50 time(k)=k*ts; rin(k)=1.0; yout(k)=-den1(2)*y1+num1(2)*u2+num1(3)*u3; error(k)=rin(k)-yout(k); ei=ei+error(k)*ts; u(k)=1.0*error(k)+0.1*(error(k)-error1)/ts+0.50*ei;%PID控制器输出 y1=yout(k); error3=error2;error2=error1;error1=err

20、or(k); u3=u2;u2=u1;u1=u(k); y1=yout(k); error3=error2;error2=error1;error1=error(k);endsubplot(2,1,2)plot(time,rin,b,time,yout,g);xlabel(time(s);ylabel(rin,yout);title(PID method);grid on六、实验结果及分析：1采用大林控制算法，画出阶跃响应曲线,仿真曲线如图上2。2. 采用普通PID控制算法，画出阶跃响应曲线仿真曲线如图2下。图2 达林算法和PID算法系统输出曲线3分析采用大林控制算法与PID 控制算法相比较的

21、控制效果。 从系统的响应曲线可知，采用大林算法控制器的系统响应速度比采用PID控制器的系统的响应速度快，而且系统输出无振荡。由此可见，大林算法比PID的控制效果更好。 最少拍计算机控制系统的设计一、实验目的：1学习并掌握有纹波最少拍控制器的设计和Simulink 实现方法；2研究最少拍控制系统对三种典型输入的适应性及输出采样点间的纹波；3学习并掌握最少拍无纹波控制器的设计和Simulink实现方法；4研究输出采样点间的纹波消除方法以及最少拍无纹波控制系统对三种典型输入的适应性。二、实验工具：MATLAB 软件（7.0版本）。三、实验内容：1实验原理最少拍设计，是指系统在典型输入信号（如阶跃信号

22、、速度信号、加速度信号等）作用下，经过最少拍（有限拍）使系统输出的系统稳态误差为零。因此，最少拍控制系统也称最少拍无差系统或最少拍随动系统，它实质上是时间最优控制系统，系统的性能指标就是系统调节时间最短或尽可能短，即对闭环Z传递函数要求快速性和准确性。下面以一个具体实例介绍最少拍系统的设计和仿真。考虑图1所示的采样数字控制系统，被控对象的脉冲传递函数为图1 最少拍采样数字控制系统设采样周期T=1s ，首先求取广义被控对象的脉冲传递函数：广义被控对象我们知道，最少拍系统是按照指定的输入形式设计的，输入形式不同，数字控制器也不同。因此，对三种不同的输入信号分别进行考虑：1 单位阶跃信号： 计算可得

23、到最少拍数字控制器为检验误差序列：由误差的变换函数得知，所设计的系统当k1后，e（k）=0就是说，一拍以后，系统输出等于输入，设计正确。 单位速度信号： 原理同上，我们可以得到： 检验误差： 从E(z)看出，按单位速度输入设计的系统，当k大于等于2之后，即二拍之后，误差e（k）=0，满足题目要求。2 单位加速度信号： 可知，按加速度输入信号设计的系统当k大于等于，即三拍之后，误差e（k）=0。2实验设计将所得结果分别用Matlab中的Simulink 工具箱进行仿真，并将输入、输出和误差三条曲线放置在同一图像内，比较三种情况下的跟随特性。四、实验步骤：在三种输入（单位阶跃/速度/加速度）分别作

24、用下，运用Simulink 对其控制结果进行仿真。 单位阶跃信号： 系统Simulink 仿真模型框图如图2 所示：可得输入、输出和误差三条曲线结果图如下所示：将其数据存为矩阵形式，命名为y，在Matlab 命令窗口输入： plot(y(:,1),y(:,2:4) grid on ,legend( 输入,误差,输出) 可得更加清晰的仿真结果如图3所示：图2 单位阶跃信号输入时最少拍控制系统图3 单位阶跃信号输入时系统的仿真结果 单位速度信号： 控制系统Simulink 框图如图4 所示：图4 单位速度信号输入时最少拍控制系统在Matlab 命令窗口输入 plot(tout(:,1),y2(:,

25、2:4) grid on ,legend( 输入,误差,输出) 仿真结果为图5 单位速度信号输入时系统的仿真结果 单位加速度信号： 控制系统Simulink 框图如图6 所示：图6 单位加速度信号输入时最少拍控制系统在Matlab 命令窗口输入 plot(tout(:,1),y2(:,2:4) grid on ,legend( 输入,误差,输出) 仿真结果图如图7 所示。五、实验结果分析：由上面的仿真结果图可知，按最少拍控制系统设计出来的闭环系统，在有限拍后进入稳态，这时闭环系统输出在采样时间精确的跟踪输入信号。如单位阶跃信号在一拍后，单位速度信号在两拍后，单位加速度信号则在三拍之后。图7 单位加速度信号输入时系统的仿真结果然而，进一步研究可以发现虽然在采样时刻系统输出与所跟踪的参考输入一致，但在两个采样时刻之间，系统的输出存在着纹波或振荡。例如单位阶跃信号在一拍后的稳态响应仍有许多振荡。这种纹波不仅影响系统的控制性能，产生过大的超调和持续振荡，而且还增加了系统功率损耗和机械磨损。因此我们需要设计无纹波最少拍计算机控制系统。