CN105045233A - 火电厂热工***中基于时间量度的pid控制器的优化设计方法 - Google Patents
火电厂热工***中基于时间量度的pid控制器的优化设计方法 Download PDFInfo
- Publication number
- CN105045233A CN105045233A CN201510404996.7A CN201510404996A CN105045233A CN 105045233 A CN105045233 A CN 105045233A CN 201510404996 A CN201510404996 A CN 201510404996A CN 105045233 A CN105045233 A CN 105045233A
- Authority
- CN
- China
- Prior art keywords
- pid
- time
- pid controller
- parameter
- controlled device
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 58
- 238000013461 design Methods 0.000 title claims abstract description 41
- 230000010354 integration Effects 0.000 title abstract description 4
- 230000004069 differentiation Effects 0.000 title abstract 2
- 238000013178 mathematical model Methods 0.000 claims abstract description 9
- 238000005457 optimization Methods 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 9
- 230000009471 action Effects 0.000 claims description 7
- 238000013016 damping Methods 0.000 claims description 4
- 230000009699 differential effect Effects 0.000 claims description 4
- 238000005070 sampling Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 2
- 230000008859 change Effects 0.000 claims description 2
- 230000003247 decreasing effect Effects 0.000 claims description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 2
- 238000012544 monitoring process Methods 0.000 claims 1
- 230000000694 effects Effects 0.000 abstract description 6
- 238000012546 transfer Methods 0.000 abstract description 2
- 230000004044 response Effects 0.000 description 11
- 238000010586 diagram Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000011156 evaluation Methods 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 2
- 230000009514 concussion Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229910052760 oxygen Inorganic materials 0.000 description 2
- 239000001301 oxygen Substances 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 206010000372 Accident at work Diseases 0.000 description 1
- 101100502522 Mus musculus Fcor gene Proteins 0.000 description 1
- 235000019504 cigarettes Nutrition 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- YHXISWVBGDMDLQ-UHFFFAOYSA-N moclobemide Chemical compound C1=CC(Cl)=CC=C1C(=O)NCCN1CCOCC1 YHXISWVBGDMDLQ-UHFFFAOYSA-N 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/418—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM]
- G05B19/41885—Total factory control, i.e. centrally controlling a plurality of machines, e.g. direct or distributed numerical control [DNC], flexible manufacturing systems [FMS], integrated manufacturing systems [IMS] or computer integrated manufacturing [CIM] characterised by modeling, simulation of the manufacturing system
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02P—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN THE PRODUCTION OR PROCESSING OF GOODS
- Y02P90/00—Enabling technologies with a potential contribution to greenhouse gas [GHG] emissions mitigation
- Y02P90/02—Total factory control, e.g. smart factories, flexible manufacturing systems [FMS] or integrated manufacturing systems [IMS]
Landscapes
- Engineering & Computer Science (AREA)
- Feedback Control In General (AREA)
- Manufacturing & Machinery (AREA)
- General Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
Abstract
一种火电厂热工***中基于时间量度的PID控制器的优化设计方法,在不完全知道被控对象数学模型的情况下,PID控制器的传递函数为:GPID(s)=Kf(Kp+Ki/s+Kds/(Tds+1);设计整定第一步是设置PID控制器Kp,Ti,Kd的参数范围,应用黄金分割法确定Kp,Ti,Kd的数值;设计整定第二步是观察在第一步设置PID控制器Kp,Ti,Kd的参数下,***闭环输出,然后按照黄金分割法进行2—3次参数整定即可获得满意的控制效果。本发明能够较好较快的完成调试过程,并且能够获得满意的控制品质。
Description
技术领域
本发明涉及一种在火电厂热工***中,基于时间量度的无模型PID控制器的优化整定设计方法,设计中应用了黄金分割优化方法,与其它应用黄金分割优化方法设计PID控制器参数不同的是,本发明设计方法不需要被控制对象的数学模型,以及增加对前置系数Kf的设计,这是本发明设计方法的主要特点。它应用于火电厂热工***控制,属于自动控制技术领域。
背景技术
在工业过程控中,PID控制器是应用最广泛的控制器。通常,工程师设计PID控制器,需要知道被控对象的精确模型,但是在实际工程中,精确模型往往很难获得。在无模型的情况下,设计和整定PID控制器参数是一项复杂的调试过程,如果现场的工程技术人员不了解被控对象的特性,不适当的调试方法可能使***发散,引发生产事故,这种调试过程是一项费时费力的工作,还不一定能够获得满意的控制品质。本专利申请发明和设计了基于***时间量度慨念的无模型PID控制器的设计和整定方法,它能够较好较快的完成调试过程,并且能够获得满意的控制品质。程序中应用了黄金分割优化方法,与其它应用黄金分割优化方法设计PID控制器参数不同的是,本专利申请设计方法不需要被控制对象的数学模型,以及增加通过对前置系数Kf的设计,这是本专利申请设计方法的主要特点。
发明内容
现代控制***中的控制器都是应用计算机实现的,每个控制***都有一定的控制周期(或者采样周期)ts,ts是与计算机的主频率相关的。
定义:一个模拟量控制***在其控制(运行)周期ts确定的情况下,控制指令(一般是阶跃信号)发出后,***在没有超过控制指令值的情况下,快速平稳完成(达到)控制指令值的时间T=nts。其中T是***无超调过渡过程时间,n是控制***的时间量度。
本专利申请发明和设计的基于***时间量度慨念的无模型PID控制器的设计和整定方法,就是根据控制***的控制周期(或者采样周期)ts而出发的。
基于时间量度的PID控制器参数整定设计步骤
步骤一:在被控对象不包含不稳定的极点和纯滞后时间,由PID控制器构建的控制***中,PID控制器的积分时间Ti≈(1/3~1/6)nts其中n是控制***的时间量度。设被控对象为G(s),如果经过PID补偿后的开环***Q(s)=G(s)GPID(s)≈1/Tis,则闭环***H(s)=1/(Tis+1),这个***完成控制指令(阶跃信号)的时间T≈6Ti。
如果经过PID补偿后的开环***(大多数情况下,开环***的传递函数可以与此等价):
Q(s)=G(s)GPID(s)=1/(cs2+ds)=1/(Tis(Tis/4+1))
由于***无超调,阻尼系数ζ=1,则c=1/w2,d=2/w,Ti=d。闭环***H(s)=1/(cs2+ds+1),完成控制指令的时间T≈7/w,则Ti≈(2/7)T。
式中符号说明如下:T≈7/w,w是***无阻尼震动频率。
这是本发明专利设计PID控制器中积分时间的根据。
传统的PID控制器的传递函数为:GPID(s)=Kp+Ki/s+Kds
步骤二:本发明改进的PID控制器的传递函数为:
GPID(s)=Kf(Kp+Ki/s+Kds/(Tds+1))(1)
其中Kf,Td的作用是非常重要的,选择合适的Kf可以对大纯滞后***的PID控制器设计取得事半功倍的作用。
Kf=f(τ/tp)=f(h),h=τ/tp;其中,τ是***纯滞后时间,tp是在无纯滞后时间τ的情况下,***的过渡过程时间,Kf是关于h的单调递减函数,0<Kf≤1,(当τ=0时,Kf=1)。本发明通过大量工程数据和曲线拟合的方法得到下面的公式:
Kf=f(τ/tp)
Kf=f(h)=c0+c1h+c2h2+c3h3+c4h4+…,(2)
或者
实践证明,当系数h=τ/tp增大时,为了保证控制***的稳定和品质,必须减小Kf。
上面计算Kf公式的非线性,恰好补偿了被控制对象G(s)=k0e-τs/(as2+bs+1)=e-τsG1(s)中非最小相位部分e-τs对最小相位部分k/(as2+bs+1)的影响。
把传统PID控制器中微分Kds改进为Kds/(Tds+1),能够有效改进微分的冲击作用,通过Td的设置能够增加微分作用的时间,减少微分带来的干扰。
步骤三:在不完全知道被控对象数学模型的情况下,若设控制***的时间尺度为n,依据上述命题,可以按照下面方法整定PID控制器参数。
实际工程应用中,在分散控制***(DCS)控制周期ts确定的情况下,这种PID控制器参数整定方法不需要知道被控对象的数学模型,只需大致确定***的1个基本参数,即时间尺度n。***完成控制指令的时间T一般在现场也容易得到(T=nts),这反映了控制工程师对被控***特性的基本了解和认识,例如在火力发电厂对汽温、水位、汽压、负荷等***中,T是不同的。
如果不能够估计出n,就只能根据***的控制周期(或者采样周期)ts来进行优化设计。
(1)如果能够估计出n,设置积分时间Ti=nits,其中ni=(1/3~1/6)n。如果不能够估计出n,设置积分时间Ti=(60~200)ts,初步设定PID参数如下:
PID控制器前置系数Kf=1,比例系数Kp=(0.3~1.5);,微分系数Kd=0,微分时间Td=(1~5)ts。一般的,为了防止***输出出现大的波动,积分时间需要选择大一些,取Ti=200ts,在Kp,Ti,Kd这组参数范围内去控制***,能够保证被控量不会出现大的波动。应用黄金分割方法,可以在这组参数范围内确定一个Kp,Ti,Kd的值。
(2)应用上述第一步设定的PID参数,如果被控对象输出比较慢,应用黄金分割法减少积分时间Ti,增加比例系数Kp(如果能够估计出n,设置Ti=nts/(4.5~)6;)或者增大Kf,Kf=(1.1~1.7)。
应用上述第一步设定的PID参数,如果被控对象的输出太快且波动大,说明被控对象包含纯滞后时间环节,或者被控对象的阻尼系数ζ<0.76,可以应用黄金分割法增加积分时间Ti,减少比例系数Kp,(如果能够估计出n,设置Ti=nts/(2.5~3.5))。如果被控对象的输出仍有波动,但波动有所减少,可以加进微分Kds/(Tds+1)作用(Kd=Ti/4,Td=(1~5)ts);如果加进微分作用后,被控对象输出仍然有波动,可以减小Kf,取Kf=0.1~0.7。
(3)如果经过(1),(2)步设计仍然不能够获得满意的控制品质,需要回到第一步重新设置Kp,Ti,Kd的取值范围,重复(1)和(2)步的设计直到获得满意的控制品质。
按照上述(1),(2)和(3)的步骤设计整定PID控制器参数,这种基于时间尺度的PID参数整定方法,在第(2)步中,改变Kf等价于同时改变了其它3个参数Kp,Ti,Kd,这是一种简单有效的方法,一般2~4次即可以获得比较满意的控制品质。
优点及功效:本发明一种火电厂热工***中基于时间量度的PID控制器的优化设计方法,其优点是基于***时间量度慨念的无模型PID控制器的设计和整定方法,它能够较好较快的完成调试过程,并且能够获得满意的控制品质。程序中应用了黄金分割优化方法,与其它应用黄金分割优化方法设计PID控制器参数不同的是,本专利申请设计方法不需要被控制对象的数学模型,以及增加通过对前置系数Kf的设计,这是本专利申请设计方法的主要特点。
附图说明
图1典型的PID控制***示意图。
图2第一次设计整定的PID参数控制***闭环输出响应示意图。
图3第二次设计整定的PID参数控制***闭环输出响应示意图。
图4PID控制器中P,I,D的控制作用示意图。
图5a、b是设计的程序界面示意图。
图6是计算程序方框图。
具体实施方式:
由本发明设计的整定的PID控制器需要在用户的分散式控制***(DCS)上组态实现,然后进行实时运行控制。也可以在工控机上实施。
按照发明内容中的设计方法,结合PID控制器中比例P,积分I和微分D的作用(见说明书附图4),本发明设计的整定的PID控制器参数程序如下图所示:
程序中应用了黄金分割(0.618)优化方法,与其它应用黄金分割优化方法设计PID控制器参数不同的是,本专利设计方法不需要被控制对象的数学模型,以及增加前置系数Kf的设计,这是本发明设计方法的主要特点。
所谓黄金分割(0.618)优化方法程序是:
function[xo,fo]=Opt_Golden(f,a,b,TolX,TolFun,k)
%%%%黄金搜索算法求在区间[a,b]上的最优化解
%f为目标函数,TolX为x阈值,TolFun为函数阈值,k为迭代次数
r=(sqrt(5)-1)/2=0.618;%r为黄金分割点值,
h=b-a;%区间宽度
rh=r*h;
%%%取两点c、d,并计算相应的函数值fc和fd
c=b-rh;
d=a+rh;
fc=feval(f,c);
fd=feval(f,d);
%%%算法第二步判断是否停止迭代
%%%%算法第三步,进行新一轮迭代
上面程序Opt_Golden(f,a,b,TolX,TolFun,k)中
性能指标f按照下面的分析的计算方法确定
由扰动分析可知,当***的扰动为确定性扰动时,其控制***的性能评价是无法直接应用最小方差指标的。通常在确定性扰动下,***的性能分为动态性能和稳态性能,从研究领域来看又分为时域性能与频域性能。时域的稳态性能通常采用稳态误差。时域的动态性能指标通常如下:
上升时间:单位阶跃响应从稳态值的10%上升到稳态值的90%,所需要的时间。对于***单位阶跃响应单调变化的***如此,而对有震荡的***,可以被定义为从初始值到第一次到达稳态的时间。上升时间,反应的是***的响应速度。
峰值时间:对于有震荡的***,其单位阶跃响应穿越终值达到第一个峰值所用的时间。
调节时间:单位阶跃响应首次达到稳态值的5%(2%)的误差带并保持的时间。
超调量:当***的单位阶跃有超调时,其单位阶跃响应的最大峰值减去稳态值的差与稳态值之比的百分数。
在频域中也有一系列的指标来反应***的性能,如反映***稳态精度的零频值,反映***超***况的谐振峰值,反映***快速响应特性与抗干扰滤波特性的带宽,反映***稳态特性的相位裕度和幅值裕度等。
误差积分性能指标
在上节中介绍的各种指标是单项的性能指标,对于一些综合性能指标包括有误差的积分形式有各种形式,常用的有以下几种:
1)误差积分(IE)
2)绝对误差指标(IAE)
3)平方误差积分(ISE)
4)时间与误差绝对值积分(ITAE)
上述各种指标的估计的是针对整个过程的情况但侧重点各有不同。在实际性能评价中更为简洁实用的误差平和(SSE):
基于此可以定义***的确定性性能指标:
性能指标的统一形式
不论回路中存在的是确定性扰动还是随机性扰动,B.Huang将其统一到H2范数意义下的性能指标:
当***仅有随机性扰动时,假定***的设定值为0,即当随机性扰动αt≠0,确定性扰动时的性能评价如下:
其闭环H2范数定义如下:
其中:Gcl为控制***闭环传递函数
随机性扰动到控制误差的传递函数写成多项式的形式:
此时***的随机性能指标定义如下:
ηs越接近于1,认为***的性能越好。从式右边部分可以得出其随机性能指标的计算方法。对于最小方差的估计可以通过FCOR算法或者ARIMA建模来完成。
当***仅有设定值扰动时,对其进行性能评价和其在仅有随机性扰动情况下的评价类似,从控制误差入手,定义其跟踪性能指标。此时,αt=0式易得:
此时,其闭环H2范数做如下定义:
此时的跟踪性能(确定性性能)指标[44]定义如下:
从该指标的定义容易看出,在死区过后***越快的完全跟踪上阶跃信号,则ηd越接近于1,***的跟踪性能越好。
其它性能指标
控制***的性能评价指标除了最小方差,还有许多其他指标。
首先基于历史数据的指标是容易想到的而且实用的。运行人员通过自己的观察和经验发现某一控制***,在某一段时间内的控制效果较为理想,那么把这一时期内的运行数据作为基础,按照某一算法计算此时的值(如方差)作为对该***评价的基准。于此可得基于历史数据的性能指标为:
但是可以肯定的是该指标越接近于1则认为***性能的性能较理想。其缺点是该指标值不像最小方差评价指标限定在[0,1]之间,其范围是不定的。但是其优点也是很明显的,该指标克服了最小方差性能指标在某些情况下使得控制作用太激烈而无法实现的缺点,因为基于历史数据的性能值一定是可达的,而且计算是方便的,且易于理解。
其次,广义的最小方差(GeneralisedMinimumVariance,GMV)也是一个常用的指标。在计算MV指标时,计算的是***输出yt或者误差et的方差最下的控制律,而此控制律常常导致***的控制作用过激而无法在实际中实现最小方差控制。为克服这一缺点,GMV通过控制量ut与误差et进行加权构造新变量,控制律通过将该新变量方差最小化实现。通过对新变量的加权系数的选择,对***的控制作用进行考虑避免控制过程过激,从而使得取得的控制律在实际中可用。所设计的程序界面如图5a、图5b,计算程序方框图见图6。
本发明设计把人工经验(增加前置系数Kf和对积分时间的估计)与黄金分割(0.618)优化方法相结合,设计了下面的控制***。
说明书附图1是在火电厂热工过程中设计的一个典型的PID控制***方框图。
在给煤量和引风量不变化的情况下,被控对象为烟汽含氧量和送风量的之间的传递函数为:
G(s)=1.3e-5s/(54s2+14s+1);
氧量给定值为4.3(%),控制周期ts=0.2(s);
按照上面的计算程序,优化设计过程如下:
1.第一步参数设计:取Kf=1;ti=80ts=16;kp=0.62;kd=0.0;td=1;
说明书附图2是应用本发明设计方法,第一步设计整定的PID参数控制***闭环输出响应。
2.从图2看出,第一次设计整定的PID参数输出有超调,第二次参数设计:取Kf=0.55;其它参数不变,控制***闭环输出响应如图3所表示。第二次参数设计值控制效果很好。
图4是PID控制器中P,I,D和PID的控制作用。
图5a、b是设计的程序界面示意图,图6是计算程序方框图。
Claims (1)
1.火电厂热工***中基于时间量度的PID控制器的优化设计方法,其特征在于:该方法具体步骤如下:
步骤一:在被控对象不包含不稳定的极点和纯滞后时间,由PID控制器构建的控制***中,PID控制器的积分时间Ti≈(1/3~1/6)nts其中n是控制***的时间量度;设被控对象为G(s),如果经过PID补偿后的开环***Q(s)=G(s)GPID(s)≈1/Tis,则闭环***H(s)=1/(Tis+1),这个***完成控制指令即阶跃信号的时间T≈6Ti;
如果经过PID补偿后的开环***为
Q(s)=G(s)GPID(s)=1/(cs2+ds)=1/(Tis(Tis/4+1))
由于***无超调,阻尼系数ζ=1,则c=1/w2,d=2/w,Ti=d;闭环***H(s)=1/(cs2+ds+1),完成控制指令的时间T≈7/w,则Ti≈(2/7)T;
式中符号说明如下:T≈7/w,w是***无阻尼震动频率;
这是设计PID控制器中积分时间的根据,而传统的PID控制器的传递函数为:
GPID(s)=Kp+Ki/s+Kds;
步骤二:改进的PID控制器的传递函数为:
GPID(s)=Kf(Kp+Ki/s+Kds/(Tds+1))(1)
其中Kf,Td的作用很重要,选择合适的Kf对大纯滞后***的PID控制器设计很关键;
Kf=f(τ/tp)=f(h),h=τ/tp;其中,τ是***纯滞后时间,tp是在无纯滞后时间τ的情况下,***的过渡过程时间,Kf是关于h的单调递减函数,0<Kf≤1,当τ=0时,Kf=1;通过大量工程数据和曲线拟合的方法得到下面的公式:
Kf=f(τ/tp)
Kf=f(h)=c0+c1h+c2h2+c3h3+c4h4+…,(2)
或者
实践证明,当系数h=τ/tp增大时,为了保证控制***的稳定和品质,必须减小Kf;
上面计算Kf公式的非线性,恰好补偿了被控制对象G(s)=k0e-τs/(as2+bs+1)=e-τsG1(s)中非最小相位部分e-τs对最小相位部分k/(as2+bs+1)的影响;
把传统PID控制器中微分Kds改进为Kds/(Tds+1),能够有效改进微分的冲击作用,通过Td的设置能够增加微分作用的时间,减少微分带来的干扰;
步骤三:在不完全知道被控对象数学模型的情况下,若设控制***的时间尺度为n,依据上述命题,按照下面方法整定PID控制器参数;
实际工程应用中,在分散控制***DCS控制周期ts确定的情况下,这种PID控制器参数整定方法不需要知道被控对象的数学模型,只需大致确定***的1个基本参数,即时间尺度n;***完成控制指令的时间T一般在现场也容易得到T=nts,这反映了控制工程师对被控***特性的基本了解和认识,在火力发电厂对汽温、水位、汽压、负荷***中,T是不同的;
如果不能够估计出n,就只能根据***的控制周期或者采样周期ts来进行优化设计;
(1)如果能够估计出n,设置积分时间Ti=nits,其中ni=(1/3~1/6)n;如果不能够估计出n,设置积分时间Ti=(60~200)ts,初步设定PID参数如下:
PID控制器前置系数Kf=1,比例系数Kp=(0.3~1.5);,微分系数Kd=0,微分时间Td=(1~5)ts;一般的,为了防止***输出出现大的波动,积分时间需要选择大一些,取Ti=200ts,在Kp,Ti,Kd这组参数范围内去控制***,能够保证被控量不会出现大的波动;应用黄金分割方法,在这组参数范围内确定一个Kp,Ti,Kd的值;
(2)应用上述第一步设定的PID参数,如果被控对象输出比较慢,应用黄金分割法减少积分时间Ti,增加比例系数Kp,如果能够估计出n,设置Ti=nts/(4.5~6,)或者增大Kf,Kf=(1.1~1.7);
应用上述第一步设定的PID参数,如果被控对象的输出太快且波动大,说明被控对象包含纯滞后时间环节,或者被控对象的阻尼系数ζ<0.76,应用黄金分割法增加积分时间Ti,减少比例系数Kp,如果能够估计出n,设置Ti=nts/(2.5~3.5));如果被控对象的输出仍有波动,但波动有所减少,加进微分Kds/(Tds+1)作用Kd=Ti/4,Td=(1~5)ts;如果加进微分作用后,被控对象输出仍然有波动,减小Kf,取Kf=0.1~0.7;
(3)如果经过(1),(2)步设计仍然不能够获得满意的控制品质,需要回到第一步重新设置Kp,Ti,Kd的取值范围,重复(1)和(2)步的设计直到获得满意的控制品质;
按照上述(1),(2)和(3)步的设计整定PID控制器参数,这种基于时间尺度的PID参数整定方法,在第(2)步中,改变Kf等价于同时改变了其它3个参数Kp,Ti,Kd,这是一种简单有效的方法,一般2~4次即能获得满意的控制品质。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510404996.7A CN105045233B (zh) | 2015-07-10 | 2015-07-10 | 火电厂热工***中基于时间量度的pid控制器的优化设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510404996.7A CN105045233B (zh) | 2015-07-10 | 2015-07-10 | 火电厂热工***中基于时间量度的pid控制器的优化设计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105045233A true CN105045233A (zh) | 2015-11-11 |
CN105045233B CN105045233B (zh) | 2017-09-19 |
Family
ID=54451837
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510404996.7A Active CN105045233B (zh) | 2015-07-10 | 2015-07-10 | 火电厂热工***中基于时间量度的pid控制器的优化设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105045233B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105487375A (zh) * | 2015-12-31 | 2016-04-13 | 中国船舶重工集团公司第七一九研究所 | 一种离散pid控制器参数整定方法 |
CN105700353A (zh) * | 2016-01-30 | 2016-06-22 | 河南城建学院 | 一种基于差分进化算法的pid控制器参数优化整定方法 |
CN106054610A (zh) * | 2016-06-23 | 2016-10-26 | 东南大学 | 一种优化的pi控制器参数工程整定方法 |
CN106094510A (zh) * | 2016-06-30 | 2016-11-09 | 电子科技大学 | 一种基于干扰补偿器的pid参数调节方法 |
CN107361760A (zh) * | 2017-07-12 | 2017-11-21 | 中国科学院上海微***与信息技术研究所 | 心磁图仪诊断***及应用该***的心磁图仪 |
CN110083877A (zh) * | 2019-03-29 | 2019-08-02 | 北京国电龙源环保工程有限公司 | 一种波动热工***传递函数建模方法 |
CN113110034A (zh) * | 2021-05-10 | 2021-07-13 | 常州市新港热电有限公司 | 一种基于dcs的引风机模糊pid控制*** |
CN113139291A (zh) * | 2021-04-23 | 2021-07-20 | 广东电网有限责任公司电力科学研究院 | 一种获取被控过程最优滑动窗滤波模型的方法及装置 |
CN113325692A (zh) * | 2021-04-28 | 2021-08-31 | 东南大学 | 基于近邻等效的pid控制器拉回式整定方法 |
CN113708952A (zh) * | 2020-10-22 | 2021-11-26 | 天翼智慧家庭科技有限公司 | 用于保障web服务连接稳定性的方法和装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10302585A1 (de) * | 2003-01-22 | 2004-08-12 | Endress + Hauser Wetzer Gmbh + Co Kg | Verfahren zum Einstellen eines Reglers |
US20050065621A1 (en) * | 2000-06-20 | 2005-03-24 | Danyang Liu | Methods of designing optimal linear controllers |
CN101069902A (zh) * | 2007-06-13 | 2007-11-14 | 湖南大学 | 基于数值模拟的高强铝、镁合金等温挤压方法 |
JP2012190364A (ja) * | 2011-03-11 | 2012-10-04 | Hiroshima Univ | 制御パラメータ調整装置及び制御パラメータ調整方法 |
CN103439880A (zh) * | 2013-08-30 | 2013-12-11 | 上海电力学院 | 基于mcp标准传递函数的pid参数整定方法 |
CN104360595A (zh) * | 2014-11-19 | 2015-02-18 | 上海电力学院 | 无迟延惯性过程的mcp-pid参数整定方法 |
-
2015
- 2015-07-10 CN CN201510404996.7A patent/CN105045233B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050065621A1 (en) * | 2000-06-20 | 2005-03-24 | Danyang Liu | Methods of designing optimal linear controllers |
DE10302585A1 (de) * | 2003-01-22 | 2004-08-12 | Endress + Hauser Wetzer Gmbh + Co Kg | Verfahren zum Einstellen eines Reglers |
CN101069902A (zh) * | 2007-06-13 | 2007-11-14 | 湖南大学 | 基于数值模拟的高强铝、镁合金等温挤压方法 |
JP2012190364A (ja) * | 2011-03-11 | 2012-10-04 | Hiroshima Univ | 制御パラメータ調整装置及び制御パラメータ調整方法 |
CN103439880A (zh) * | 2013-08-30 | 2013-12-11 | 上海电力学院 | 基于mcp标准传递函数的pid参数整定方法 |
CN104360595A (zh) * | 2014-11-19 | 2015-02-18 | 上海电力学院 | 无迟延惯性过程的mcp-pid参数整定方法 |
Non-Patent Citations (5)
Title |
---|
何芝强: "PID控制器参数整定方法及其应用研究", 《中国优秀博硕士学位论文全文数据库(硕士)工程科技Ⅱ辑》 * |
刘聪: "锅炉快速控制原型开发及PID参数自整定方法研究", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑》 * |
刘邦利 等: "《电厂化学仪表及自动化》", 28 February 2015, 中国电力出版社 * |
朱晓青: "《过程检测控制技术与应用》", 30 June 2002, 冶金工业出版社 * |
李国林: "PID控制器参数整定技术研究与优化设计", 《中国优秀硕士学位论文全文数据库 信息科技辑》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105487375B (zh) * | 2015-12-31 | 2018-06-22 | 中国船舶重工集团公司第七一九研究所 | 一种离散pid控制器参数整定方法 |
CN105487375A (zh) * | 2015-12-31 | 2016-04-13 | 中国船舶重工集团公司第七一九研究所 | 一种离散pid控制器参数整定方法 |
CN105700353A (zh) * | 2016-01-30 | 2016-06-22 | 河南城建学院 | 一种基于差分进化算法的pid控制器参数优化整定方法 |
CN106054610A (zh) * | 2016-06-23 | 2016-10-26 | 东南大学 | 一种优化的pi控制器参数工程整定方法 |
CN106054610B (zh) * | 2016-06-23 | 2019-01-25 | 东南大学 | 一种优化的pi控制器参数工程整定方法 |
CN106094510A (zh) * | 2016-06-30 | 2016-11-09 | 电子科技大学 | 一种基于干扰补偿器的pid参数调节方法 |
CN106094510B (zh) * | 2016-06-30 | 2019-11-05 | 电子科技大学 | 一种基于干扰补偿器的pid参数调节方法 |
CN107361760A (zh) * | 2017-07-12 | 2017-11-21 | 中国科学院上海微***与信息技术研究所 | 心磁图仪诊断***及应用该***的心磁图仪 |
CN110083877A (zh) * | 2019-03-29 | 2019-08-02 | 北京国电龙源环保工程有限公司 | 一种波动热工***传递函数建模方法 |
CN113708952A (zh) * | 2020-10-22 | 2021-11-26 | 天翼智慧家庭科技有限公司 | 用于保障web服务连接稳定性的方法和装置 |
CN113708952B (zh) * | 2020-10-22 | 2024-06-07 | 天翼数字生活科技有限公司 | 用于保障web服务连接稳定性的方法和装置 |
CN113139291A (zh) * | 2021-04-23 | 2021-07-20 | 广东电网有限责任公司电力科学研究院 | 一种获取被控过程最优滑动窗滤波模型的方法及装置 |
CN113325692A (zh) * | 2021-04-28 | 2021-08-31 | 东南大学 | 基于近邻等效的pid控制器拉回式整定方法 |
CN113325692B (zh) * | 2021-04-28 | 2022-12-27 | 东南大学 | 基于近邻等效的pid控制器拉回式整定方法 |
CN113110034A (zh) * | 2021-05-10 | 2021-07-13 | 常州市新港热电有限公司 | 一种基于dcs的引风机模糊pid控制*** |
Also Published As
Publication number | Publication date |
---|---|
CN105045233B (zh) | 2017-09-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105045233A (zh) | 火电厂热工***中基于时间量度的pid控制器的优化设计方法 | |
JP5005695B2 (ja) | 適応多変数mpcコントローラ | |
CN103558755B (zh) | 分数阶积分pid控制器整定和自整定方法 | |
CN103676651B (zh) | 基于状态观测模型的锅炉汽温预测控制方法 | |
CN103406909B (zh) | 一种机械臂***的跟踪控制设备及方法 | |
CN113031565B (zh) | 一种火电机组一次调频功率响应性能预测方法及*** | |
CN107450310B (zh) | 一种连续过程模型的设定值激励闭环辨识方法 | |
CN106530118A (zh) | 多区互联电力***负荷频率的约束gpc优化控制方法 | |
Priya et al. | Particle swarm optimisation applied to real time control of spherical tank system | |
CN101126917A (zh) | 槽式反应器基于多核支持向量机的非线性预测控制方法 | |
Chertovskikh et al. | An adaptive PID controller with an online auto-tuning by a pretrained neural network | |
Anusha et al. | Comparison of tuning methods of PID controller | |
Folorunso et al. | Control of a two layered coupled tank: Application of IMC, IMC-PI and pole-placement PI controllers | |
CN104808486A (zh) | 基于模糊ts模型的压电陶瓷执行器的预测控制方法和装置 | |
Yang et al. | Application of fuzzy neural network PID algorithm in oil pump control | |
Maasoumy et al. | Comparison of control strategies for energy efficient building HVAC systems | |
CN105607486A (zh) | 一种通用的针对一系列化工积分时滞过程的控制方法 | |
Škarda et al. | Bode-like control loop performance index evaluated for a class of fractional-order processes | |
Zhang et al. | Nonlinear model predictive control of ultra-supercritical once through boiler-turbine unit | |
Grelewicz et al. | Practical Verification of the Advanced Control Algorithms Based on the Virtual Commissioning Methodology-A Case Study | |
Campo | PID control design | |
Ionesi et al. | On-line parameter and state estimation of an air handling unit model: experimental results using the modulating function method | |
Dulău et al. | Conventional control vs. robust control on heat-exchangers | |
Ruusu et al. | Sliding mode SISO control of model parameters for implicit dynamic feedback estimation of industrial tracking simulation systems | |
Aravind et al. | Reliability of pid controller for a real time sensor interfaced temperature process |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |