(3)分析导致上述两种情况下
取值范围差异的原因。
答:
连续系统比离散系统稳定性好,加入采样开关以后,采样周期越大,离散系统系统稳定性越差,能使系统稳定的K的范围越小。
2.时域特性分析
令
(1)假设不考虑采样开关和零阶保持器的影响,即看作一连续系统,观察其单位阶跃响应,记录上升时间、超调量、调节时间、峰值时间等一系列的时域性能指标;
G1=tf([20],[110]);
step(feedback(G1,1));
由图数据:
上升时间=0.254s
超调量=70.2%
调节时间=7.82s
峰值时间=0.702s
(2)考虑采样开关和零阶保持器的影响,观察其单位阶跃响应,记录上升时间、超调量、调节时间、峰值时间等一系列的时域性能指标;
G1=tf([20],[110]);
G=c2d(G1,0.01,'zoh');
Q=step(feedback(G,1));
[num,den]=tfdata(Q);
dstep(num,den)
上升时间=24.7*T=0.247s
超调量=72.8%
调节时间=857*T=8.57s
峰值时间=71*T=0.71s
(3)分析其时域性能指标的差异及产生原因。
由于采样开关和零阶保持器的存在,使得离散系统的时域响应与连续系统相比发生变化,稳定性相对降低,动态性能相对变差。
3.频域特性分析
(1)假设不考虑采样开关和零阶保持器的影响,即看作一连续系统,绘制其频率特性响应;
G1=tf([1],[110]);
bode(G1)
(2)考虑采样开关和零阶保持器的影响,绘制其频率特性响应;;
w=logspace(-1,10);
dbode([1],[110],0.01,w)
(2)讨论上述两种情况下频率特性响应的区别和联系。
G1=tf([1],[110]);
bode(G1)
holdon
w=logspace(-1,10);
dbode([1],[110],0.01,w)
频率特性对比:
1、频率特性是w的周期函数,当wT沿着单位圆每转一周时,频率特性周期性重复一次,这是连续系统没有的。
2、幅频特性是w的偶函数,相频特性是w的奇函数,连续系统也有这个特性。
3、离散环节频率特性形状与连续系统频率特性形状有较大差别,特别是当采样周期较大以及频率较高时,由于混叠,使频率特性形状有较大变化,主要表现有:
高频时会出现多个峰值;
可能出现正相位;
仅在较小的采样周期或低频段与连续系统频率特性相接近。
南京邮电大学自动化学院
实验报告
课程名称:
计算机控制系统
实验名称:
数字PID控制
所在专业:
自动化
学生姓名:
王站
班级学号:
B11050107
任课教师:
程艳云
2013/2014学年第二学期
实验二数字PID控制
一、实验目的
1.研究PID控制器的参数对系统稳定性及过渡过程的影响。
2.研究采样周期T对系统特性的影响。
3.研究I型系统及系统的稳定误差。
二、实验仪器
1.EL-AT-III型计算机控制系统实验箱一台
2.PC计算机一台
三、实验内容
1.系统结构图如3-1图。
图3-1系统结构图
图中Gc(s)=Kp(1+Ki/s+Kds)
Gh(s)=(1-e-TS)/s
Gp1(s)=5/((0.5s+1)(0.1s+1))
Gp2(s)=1/(s(0.1s+1))
2.开环系统(被控制对象)的模拟电路图如图3-2和图3-3,其中图3-2对应GP1(s),图3-3对应Gp2(s)。
图3-2开环系统结构图1图3-3开环系统结构图2
3.被控对象GP1(s)为“0型”系统,采用PI控制或PID控制,可系统变为“I型”系统,被控对象Gp2(s)为“I型”系统,采用PI控制或PID控制可使系统变成“II型”系统。
4.当r(t)=1(t)时(实际是方波),研究其过渡过程。
5.PI调节器及PID调节器的增益
Gc(s)=Kp(1+K1/s)
=KpK1((1/k1)s+1)/s
=K(Tis+1)/s
式中K=KpKi,Ti=(1/K1)
不难看出PI调节器的增益K=KpKi,因此在改变Ki时,同时改变了闭环增益K,如果不想改变K,则应相应改变Kp。
采用PID调节器相同。
6.“II型”系统要注意稳定性。
对于Gp2(s),若采用PI调节器控制,其开环传递函数为
G(s)=Gc(s)·Gp2(s)
=K(Tis+1)/s·1/s(0.1s+1)
为使用环系统稳定,应满足Ti>0.1,即K1<10
7.PID递推算法如果PID调节器输入信号为e(t),其输送信号为u(t),则离散的递推算法如下:
u(k)=u(k-1)+q0e(k)+q1e(k-1)+q2e(k-2)
其中q0=Kp(1+KiT+(Kd/T))
q1=-Kp(1+(2Kd/T))
q2=Kp(Kd/T)
T--采样周期
四、实验步骤
1.连接被测量典型环节的模拟电路(图3-2)。
电路的输入U1接A/D、D/A卡的DA1输出,电路的输出U2接A/D、D/A卡的AD1输入。
检查无误后接通电源。
2.启动计算机,双击桌面“计算机控制实验”快捷方式,运行软件。
3.测试计算机与实验箱的通信是否正常,通信正常继续。
如通信不正常查找原因使通信正常后才可以继续进行实验。
4.在实验项目的下拉列表中选择实验三[数字PID控制],鼠标单击鼠标单击
按钮,弹出实验课题参数设置窗口。
5.输入参数Kp,Ki,Kd(参考值Kp=1,Ki=0.02,kd=1)。
6.参数设置完成点击确认后观察响应曲线。
若不满意,改变Kp,Ki,Kd的数值和
与其相对应的性能指标p、ts的数值。
7.取满意的Kp,Ki,Kd值,观查有无稳态误差。
8.断开电源,连接被测量典型环节的模拟电路(图3-3)。
电路的输入U1接A/D、D/A卡的DA1输出,电路的输出U2接A/D、D/A卡的AD1输入,将纯积分电容的两端连在模拟开关上。
检查无误后接通电源。
9.重复4-7步骤。
10.计算Kp,Ki,Kd取不同的数值时对应的p、ts的数值,测量系统的阶跃响应曲线及时域性能指标,记入表中:
实验结果
参数
δ%
Ts
阶跃响应曲线
Kp
Ki
Kd
1
0
0
0
0.623
如图1
1
0.02
1
13.8%
0.654
如图2
1
0.1
0
63.1%
5.8
如图3
1
0.05
0
35.2%
1.392
如图4
1
0.02
0
15.9%
0.821
如图5
1
0.02
2
11.6%
0.696
如图6
实验结果
参数
δ%
Ts
阶跃响应曲线
Kp
Ki
Kd
1
0
0
0
0.623
如图1
1
0.05
0
95.3%
1.715
如图7
1
0.03
0
63.1%
1.569
如图8
1
0.02
0
31%
0.840
如图9
1
0.02
1
22.4%
0.738
如图10
图1
图2
图3
图4
图5
图6
图7
图8
图9
图10
五、实验报告
1.画出所做实验的模拟电路图。
答:
所做实验的模拟电路图如图3-2和图3-3所示.
2.当被控对象为Gp1(s)时取过渡过程为最满意时的Kp,Ki,Kd,画出校正后的Bode图,查出相稳定裕量和穿越频率c。
答:
取kp=1,ki=0.01,kd=1
Gc(s)=1+0.02/s+s
Gp1(s)=5/((0.5s+1)(0.1s+1))
根据MATLAB指令:
G=tf([550.1],[0.050.610]);
bode(G)
由伯德图可以得到:
相稳定裕量=96deg,穿越频率c=100rad/s
3.总结一种有效的选择Kp,Ki,Kd方法,以最快的速度获得满意的参数。
答:
根据给定的参考值kp,ki,kd,运用控制变量法,每次保证其中两个值不变,改变另外的一个值,然后画出阶跃响应曲线,通过对比每次得到的曲线,找出较为理想的曲线,即可获得较为满意的参数kp,ki,kd。
六、预习要求
1.熟悉PID控制器系统的组成。
2.熟悉PID控制器的参数对系统稳定性的影响。
七、PID软件流程图
图中ek为误差,ek1为上一次的误差,ek2为误差的累积和,uk是控制量
南京邮电大学自动化学院
实验报告
课程名称:
计算机控制系统
实验名称:
最小拍控制系统
所在专业:
自动化
学生姓名:
王站
班级学号:
B11050107
任课教师:
程艳云
2013/2014学年第二学期
实验三:
最小拍控制系统
一、实验目的:
1.建立计算机最小拍控制系统的一般概念;
2.掌握有纹波最小拍控制器的设计方法
3.观察无纹波最小拍控制器的设计方法;
4.了解最小拍控制器的优缺点;
5.掌握最小拍控制系统的改进方法。
二、实验内容:
图1
采样周期为T=0.1s
1.针对图一所示的计算机控制系统,考虑输入为单位速度信号时,进行计算机控制算法D(Z)设计,编程实现最小拍有纹波系统;
第一步:
求广义脉冲传递函数
Gs=tf([5],[1,1,0]);
Gz=c2d(Gs,0.1,'zoh');//求解广义对象的脉冲传递函数
Transferfunction:
0.02419z+0.02339
----------------------
z^2-1.905z+0.9048
Samplingtime:
0.1
第二步:
确定闭环脉冲传递函数
由广义对象的闭环脉冲传递函数得延迟因子为
,则
应包含该延迟因子
[num,den]=c2dm([5],[1,1,0],0.1,'zoh')
num=
00.02420.0234
den=
1.0000-1.90480.9048
tf2zpk(num)
ans=
-0.9672
由此可得:
G(z)没有单位圆上或圆外的零点
tf2zpk(den)
ans=
1.0000
0.9048
由此可得:
G(z)包含一个单位圆上极点
根据物理可实现条件和稳定性条件,
所以
Qez=tf([1-21],[1,0,0],0.1);
Qz=1-Qez
Transferfunction:
2z-1
-------
z^2
Samplingtime:
0.1
Dz=1/Gz*Qz/(1-Qz)
Transferfunction:
2z^5-4.81z^4+3.715z^3-0.9048z^2
----------------------------------------------------
0.02419z^5-0.02498z^4-0.0226z^3+0.02339z^2
Samplingtime:
0.1
Qz=minreal(Dz*Gz/(1+Dz*Gz));
Q=2;
t=0:
0.1:
Q;u=t;plot(0:
0.1:
Q,u,'r*');holdon
yt=lsim(Qz,u,t,0);plot(0:
0.1:
Q,yt);
2.讨论纹波的生成原因,编程实现最小拍无纹波系统;
纹波产生原因:
控制信号u(t)的变化引起输出的波动
根据物理可实现条件和稳定性条件以及无纹波条件,
Qez=tf([1-1.2666-0.46680.7334],[1,0,0,0],0.1);
Qz=1-Qez
Transferfunction:
1.267z^2+0.4668z-0.7334
-----------------------------
z^3
Samplingtime:
0.1
设计最小拍无纹波系统控制器
Dz=1/Gz*Qz/(1-Qz)
Transferfunction:
1.267z^7-1.946z^6-0.4765z^5+1.819z^4-0.6636z^3
---------------------------------------------------------------------
0.02419z^7-0.007241z^6-0.04092z^5+0.006818z^4+0.01716z^3
Samplingtime:
0.1
Qz=minreal(Dz*Gz/(1+Dz*Gz));
Q=2;
t=0:
0.1:
Q;u=t;plot(0:
0.1:
Q,u,'r*');holdon
yt=lsim(Qz,u,t,0);plot(0:
0.1:
Q,yt);
3.讨论最小拍系统的特点,采取惯性因子法对最小拍控制器加以改进,并研究惯性因子对系统性能的影响。
最小拍系统的特点:
系统对应于典型的输入信号具有最快的响应速度,被控量在最短的时间达到设定值。
局限性:
1)仅适应与一种类型的输入信号
2)系统参数变化引起系统极点位置产生偏移,系统动态性能明显变坏。
3)采样周期的限制。
惯性因子法:
基本思想:
牺牲有限拍的性质为代价,换取系统对不同输入类型性能皆能获得比较满意的控制效果。
方法:
引进惯性因子,改进系统的闭环脉冲传递函数,使系统对多种信号的输入信号都有满意的响应。
定义
相当于增加单位圆内实轴上的极点。
南京邮电大学自动化学院
实验报告
课程名称:
计算机控制系统
实验名称:
炉温控制实验
所在专业:
自动化
学生姓名:
王站
班级学号:
B11050107
任课教师:
程艳云
2013/2014学年第二学期
实验四炉温控制实验
一、实验目的
1.了解温度控制系统的特点。
2.研究采样周期T对系统特性的影响。
3.研究大时间常数系统PID控制器的参数的整定方法。
二、实验仪器
1.EL-AT-III型计算机控制系统实验箱一台
2.PC计算机一台
3.炉温控制实验对象一台
三、炉温控制的基本原理
1.
系统结构图示于图7-1。
图7-1系统结构图
图中Gc(s)=Kp(1+Ki/s+Kds)
Gh(s)=(1-e-TS)/s
Gp(s)=1/(Ts+1)
2.系统的基本工作原理
整个炉温控制系统由两大部分组成,第一部分由计算机和A/D&D/A卡组成,主要完成温度采集、PID运算、产生控制可控硅的触发脉冲,第二部分由传感器信号放大,同步脉冲形成,以及触发脉冲放大等组成。
炉温控制的基本原理是:
改变可控硅的导通角即改变电热炉加热丝两端的有效电压,有效电压的可在0~140V内变化。
可控硅的导通角为0~5CH。
温度传感是通过一只热敏电阻及其放大电路组成的,温度越高其输出电压越小。
外部LED灯的亮灭表示可控硅的导通与闭合的占空比时间,如果炉温温度低于设定值则可控硅导通,系统加热,否则系统停止加热,炉温自然冷却到设定值。
第二部分电路原理图见附录一。
3.PID递推算法:
如果PID调节器输入信号为e(t),其输送信号为u(t),则离散的递推算法如下:
Uk=Kpek+Kiek2+Kd(ek-ek-1),其中ek2是误差累积和。
四、实验内容:
1.设定炉子的温度在一恒定值。
2.调整P、I、D各参数观察对其有何影响。
五、实验步骤
1.启动计算机,双击桌面“计算机控制实验”快捷方式,运行软件。
2.测试计算机与实验箱的通信是否正常,通信正常继续。
如通信不正常查找原因使通
信正常后才可以继续进行实验。
3.20芯的扁平电缆连接实验箱和炉温控制对象,检查无误后,接通实验箱和炉温控
制的电源。
开环控制
4.在实验项目的下拉列表中选择实验七[七、炉温控制],鼠标单击
按钮,弹出实验课题参数设置对话框。
在参数设置对话框中设置相应的实验参数后鼠标单击确认等待屏幕的显示区显示实验结果。
测量系统响应时间Ts和超调量p。
5.重复步骤4,改变参数设置,观测波形的变化,记入下表:
性能指标
占空比
阶跃响应曲线
δ%
Tp(秒)
Ts(秒)
20%
如图1
无
无
81.220
40%
如图2
无
无
75.130
60%
如图3
无
无
67.365
70%
如图4
无
无
63.114
图1
图2
图3
图4
闭环控制
6.在实验项目的下拉列表中选择实验七[七、炉温控制]鼠标单击
按钮,弹出实验
课题参数设置对话框,选择PID,在参数设置窗口设置炉温控制对象的给定温度以及
Ki、Kp、Kd值,点击确认在观察窗口观测系统响应曲线。
测量系统响应时间Ts和超
调量p。
7.重复步骤6,改变PID参数,观测波形的变化,记入下表中:
性能指标
参数
阶跃响应曲线
δ%
Tp(秒)
Ts(秒)
Kp
Ki
Kd
1.00
0.02
1.00
如图5
86%
146.597
440.721
1.00
0.40
1.00
如图6
45%
230.886
468.051
1.00
0.30
1.00
如图7
52%
47.319
267.680
1.00
0.50
1.00
如图8
54%
205.515
440.257
图5图6
图7图8
六、实验报告
1、过渡过程为最满意时的Kp, Ki, Kd及其响应曲线:
最满意时的Kp=0.5, Ki=3.5, Kd=6.5.响应曲线如图5。
2、温度控制系统有效的选择Kp, Ki, Kd方法:
本实验采用的是控制变量法来做的,即在初始数值的基础上改变(变大或者变小都行)一个数值,观察炉温曲线的变化趋势,并初步得出这个变量对于曲线哪一部分有影响。
对其中每一个可控变量应用此方法,最后根据需要改变所有变量值,这样得出的曲线应该就能更加切合需要。
3、实验中问题分析:
总的来说比例环节KP越大,快速性增强,振荡加剧;积分环节是用以消除静差;微分环节:
用以改善系统的振荡性。
七、实验总结
通过这次试验,我真正接触了PID调节原理在实际中的应用,大大加深了自己对PID调节原理的理解和认识。
PID控制原理的优点在于能够在控制过程中根据预先设定好的控制规律不停地自动调节控制量以使被控系统朝着设定的平衡状态过度,最后达到控制范围精度内的稳定的动态平衡状态。
比例系数Kp决定控制作用的强弱,Kp加大时可减少系统的稳态误差,提高系统的动态响应速度,但Kp过大会引起振荡或导致系统不稳定。
积分项Ki/S用于消除系统的稳态误差,但积分项会使动态过程变慢,增大系统的超调量,使系统的稳定性变好。
微分KdS的控制作用与偏差的变化速度有关。
微分能产生超前的校正作用,有助于减少超调和振荡,并能减少调整时间,从而改善系统的动态性能。