CN104199293A - 基于最大灵敏度指标的稳定分数阶pid参数优化方法 - Google Patents
基于最大灵敏度指标的稳定分数阶pid参数优化方法 Download PDFInfo
- Publication number
- CN104199293A CN104199293A CN201410395817.3A CN201410395817A CN104199293A CN 104199293 A CN104199293 A CN 104199293A CN 201410395817 A CN201410395817 A CN 201410395817A CN 104199293 A CN104199293 A CN 104199293A
- Authority
- CN
- China
- Prior art keywords
- parameter
- fractional order
- order pid
- omega
- pid controller
- 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.)
- Pending
Links
Landscapes
- Feedback Control In General (AREA)
Abstract
本发明公开了一种基于最大灵敏度指标的稳定分数阶PID参数优化方法,包括步骤:S1、获取分数阶PID控制器与被控对象组成的闭环***的参数稳定域,分数阶PID控制器的参数包括:kp、ki、kd、λ以及μ;S2、获取满足一预设的最大灵敏度指标的分数阶PID控制器的参数解集;S3、由参数稳定域以及参数解集的交集获得满足最大灵敏度的参数稳定解;S4、根据优化目标确定权系数大小;S5、在参数解集中选择一组最优的解;其中,优化目标的函数为:f(kp,ki,kd)=ξ1σ+ξ2ts,最优的解包括:使优化目标的函数值最小,且最优解属于参数解集,σ为超调量,ts为调节时间,ξ1,ξ2为权系数,并满足ξ1+ξ2=1。
Description
技术领域
本发明涉及计算、推算技术领域,特别涉及一种基于最大灵敏度指标的稳定分数阶PID参数最优整定方法。
背景技术
分数阶微积分作为一种数学分析工具,和整数阶微积分相比有许多优点。分数阶微积分是指微分,积分的阶次不再局限于整数,可以是任意的或者是分数。由于微分,积分阶次范围的扩大,分数阶微积分的描述能力比整数阶微积分更大。当描述对象本身是分数阶***时,分数阶微积分通常比整数阶微积分具有更加普遍的意义。
分数阶控制理论随着分数阶微积分与控制理论的结合而形成。随后分数阶控制理论应用于PID控制器参数整定,这给PID控制注入了新的活力。分数阶PID控制器是一种新型控制器,除了kp,ki,kd三个参数,它多了两个可调参数λ和μ,使得控制器设计及参数整定的难度加大,但同时却也使控制器参数调节更为灵活,控制效果更好。
然而大部分研究集中在对分数阶PID控制器参数整定,对于参数优化也是及其重要的,特别当对于控制器跟踪性能有要求时,优化就更为重要。
发明内容
本发明针对现有技术存在的上述不足,提供了一种基于最大灵敏度指标的稳定分数阶PID参数优化方法,本发明通过以下技术方案实现:
一种基于最大灵敏度指标的稳定分数阶PID参数优化方法,包括步骤:
S1、获取分数阶PID控制器与被控对象组成的闭环***的参数稳定域,分数阶PID控制器的参数包括:kp、ki、kd、λ以及μ;
S2、获取满足一预设的最大灵敏度指标的分数阶PID控制器的参数解集;
S3、由参数稳定域以及参数解集的交集获得满足最大灵敏度的参数稳定解;
S4、根据优化目标确定权系数大小;
S5、在参数解集中选择一组最优的解;
其中,优化目标的函数为:f(kp,ki,kd)=ξ1σ+ξ2ts,最优的解包括:使优化目标的函数值最小,且最优解属于参数解集,其中,σ为超调量,ts为调节时间,ξ1,ξ2为权系数,并满足ξ1+ξ2=1。
较佳的,采用D分割原理获取分数阶PID控制器的参数的稳定域,稳定域为由实根边界、无穷根边界以及复根边界为边界的区域。
较佳的,获取分数阶PID控制器的参数稳定域包括:
βn>…>β1>β0≥0,αn>…>α1>α0≥0,其闭环传递函数:
将分数阶PID控制器和被控对象代入闭环传递函数,得到特征多项式:Ψ(s)=sλD(s)+(kpsλ+ki+kdsμ+λ)N(s)e-Ls,对于一组参数(kp,ki,kd,λ,μ)能使特征多项式的值为0,则分数阶PID控制器的输入输出稳定。
较佳的,获得满足最大灵敏度的参数稳定解,包括:
将s=jω,和ejx=cosx+jsinx代入特征多项式, 得到一组KP,Ki关于ω的方程,在[ωmin,ωmax]内以一定步长取ω的值,得到KP,Ki的参数稳定解的解集。
较佳的,最大灵敏度指标为:MT≤γT,γT≥1,其中, γT为预设的最大灵敏度指标。
本发明可以得到满足灵敏度约束的分数阶PID控制器,并且控制器能兼顾多个输入跟踪性能,该方法是一种有效的分数阶PID控制器参数整定方法。
附图说明
图1所示的是分数阶PID控制器与被控对象组成的闭环***的结构图;
图2所示的是最大灵敏度指标几何表示示意图;
图3所示的是分数阶PID控制器的稳定域示意图;
图4所示的是稳定域内随机得到三个参数的单位阶跃响应示意图;
图5所示的是满足最大灵敏度指标的分数阶PID控制器参数解集示意图;
图6所示的是ξ1=0.4,ξ2=0.6时的阶跃响应示意图;
图7所示的是ξ1=1,ξ2=0时的阶跃响应示意图;
图8所示的是ξ1=0,ξ2=1时的阶跃响应示意图。
具体实施方式
以下将结合本发明的附图,对本发明实施例中的技术方案进行清楚、完整的描述和讨论,显然,这里所描述的仅仅是本发明的一部分实例,并不是全部的实例,基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例,都属于本发明的保护范围。
为了便于对本发明实施例的理解,下面将结合附图以具体实施例为例作进一步的解释说明,且各个实施例不构成对本发明实施例的限定。
本发明针对的分数阶PID控制器与被控对象组成的闭环***的结构如图1所示:
式中0<λ,μ<2,C(s)代表控制器;
其中,G(s)为被控对象,βn>…>β1>β0≥0,αn>…>α1>α0≥0,
其闭环传递函数为
将式(1)和式(2)代入(3),闭环***的特征多项式写作:
Ψ(s)=sλD(s)+(kpsλ+ki+kdsμ+λ)N(s)e-Ls (4)
对于一组参数(kp,ki,kd,λ,μ),若它使特征方程式Ψ(s)=0的根都有复实部,则该子***是输入输出稳定的。
由D分割原理(出处:Serdar E H.An Algorithm for Stabilization ofFractional-Order Time Delay Systems Using Fractional Order PID Controllers[J].IEEE Trans on Automatic Control,2007,52(4):1964-1969.),可将(kp,ki,kd,λ,μ)所构成的参数空间Φ分割成RRB,CRB,IRB为边界的区域D,区域D包含所有使子***稳定的点,D分割的边界定义如下:
其中,表示实根边界,表示无穷根边界,表示复根边界。
将s=jω代入特征多项式(4),可得到特征方程如下:
Ψ(jω)=(jω)λD(jω)+(kd(jω)λ+μ+kp(jω)λ+ki)N(jω)e-Ljω (7)
根据D分割原理,对给定被控***,分数阶PIλDμ控制器的参数空间(kp,ki,kd,λ,μ)被可能的平面(a),可能的平面(b)和曲面(c)分割成若干区域。
(a)实根边界(RRB)
根据式(6)的第一个式子,可以得到实根边界为:ki=0。
(b)无穷根边界(IRB)
(c)复根边界(CRB)根据和欧拉公式ejx=cosx+jsinx,代入特征方程(7)可得
其中:
通过上式(9)和(10),对于给定的参数kd,λ,μ,当ω从0到∞变化时,可以在(kp,ki)平面得到复根边界(CRB)。
最大灵敏度指标如图2表示:
引理1设F(s)=NF(s)/DF(s)为一稳定正则实(或复)有理数,且DF(s)的阶数为α次幂系数。定义:
则对于一给定的γ>0,当且仅当同时满足下面两个条件时,||F(s)||∞<γ:1)|cα|<γ|dα|;2)Hurwize稳定,
其中:cα和dα分别为NF和DF的s的α次幂系数。
由灵敏度约束指标和引理1,可得引理1中对应的NF和DF分别为:
NF(s)=N(s)e-Ls (18)
DF(s)=sλD(s)+(kpsλ+ki+kdsλ+μ)N(s)e-Ls (19)
由条件2)可得式
求取满足最大灵敏度指标的参数解集:
其中:f1(ω),f2(ω),f3(ω),f4(ω),f5(ω),f6(ω),h1(ω),h2(ω)同前
通过上式(23),对于给定的参数kd,λ,μ,当ω从0到∞变化时,可以得到最大灵敏度指标下的分数阶PID控制器参数解集Qθ。
定义优化函数为:f(kp,ki,kd)=ξ1σ+ξ2ts (20)
参数优化问题为:
min(f(kp,ki,kd)) (21)
s.t.(kp,ki,kd)∈Φ
其中:σ为超调量,ts为调节时间,ξ1,ξ2为权系数,并满足ξ1+ξ2=1。通过改变两个权系数的比重,可以达到不同的优化目的。
参数优化步骤:
(1)求取分数阶PID控制器参数稳定域;
(2)在稳定域内求得满足灵敏度约束指标下的参数解集Φ;
(3)根据实际优化要求确定权值ξ1和ξ2的大小;
(4)在参数解集Φ中通过求解参数优化问题(21)得到一组最优的解。
仿真采用被控对象函数如下式所示
设计分数阶PID控制器,且满足灵敏度指标:MT≤γT,γT=1.6。
根据分数阶PID参数整定步骤可以确定可以求得分数阶PID控制器参数为λ=1.5,μ=0.7,kd=1时得到最大稳定。如图3所示:
在稳定域内随机取三个点,点A:kp=8.2534,ki=2.0227,点B:kp=11.1401,ki=6.9027点C:kp=19.5514,ki=10.0121,分别得到单位阶跃响应的仿真图,如图4所示,从图中可以看出,***都是稳定的。
再根据灵敏度指标约束下的控制器参数解集求取策略,得到λ=1.5,μ=0.7,kd=1时满足灵敏度约束时控制器参数解集如图5所示。
参数优化:
(1)在满足灵敏度指标的分数阶PID控制器参数解集中根据本文的优化方法对参数进行优化,选取ξ1=0.4,ξ2=0.6。得到最优参数解为kp=3.9544,ki=0.6423。再在满足灵敏度约束的参数解集中随机取一个点kp=15.0453,ki=13.3374。单位阶跃响应仿真如图6所示。从图中可以看出,最优参数解下的闭环***阶跃响应跟踪性能明显优于优化前。
(2)当ξ1=1,ξ2=0时,此时得到优化后的参数解为kp=3.0452,ki=0.3373。阶跃响应如图7所示。
(3)当ξ1=0,ξ2=1时得到优化后的参数解为kp=10.0453,ki=0.3374,此时的阶跃响如图8所示。
本发明对分数阶PID控制器的参数整定进行了研究,给出了一种基于最大灵敏度指标的分数阶PID控制器参数最优整定方法。该方法先根据D分割法求取分数阶PID控制器参数稳定域,再得到满足灵敏度约束指标下的参数解集,最后由超调量和调节时间作为优化指标函数对参数进行优化得到最优参数解。最后结果表明该方法可以得到满足灵敏度约束的分数阶PID控制器,并且控制器能兼顾多个输入跟踪性能,该方法是一种有效的分数阶PID控制器参数整定方法。
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。
Claims (5)
1.一种基于最大灵敏度指标的稳定分数阶PID参数优化方法,其特征在于,包括步骤:
S1、获取分数阶PID控制器与被控对象组成的闭环***的参数稳定域,所述分数阶PID控制器的参数包括:kp、ki、kd、λ以及μ;
S2、获取满足一预设的最大灵敏度指标的所述分数阶PID控制器的参数解集;
S3、由所述参数稳定域以及所述参数解集的交集获得满足最大灵敏度的参数稳定解;
S4、根据优化目标确定权系数大小;
S5、在参数解集中选择一组最优的解;
其中,所述优化目标的函数为:f(kp,ki,kd)=ξ1σ+ξ2ts,所述最优的解包括:使所述优化目标的函数值最小,且所述最优解属于所述参数解集,其中,σ为超调量,ts为调节时间,ξ1,ξ2为权系数,并满足ξ1+ξ2=1。
2.根据权利要求1所述的基于最大灵敏度指标的稳定分数阶PID参数优化方法,其特征在于,采用D分割原理获取分数阶PID控制器的参数的稳定域,所述稳定域为由实根边界、无穷根边界以及复根边界为边界的区域。
3.根据权利要求1所述的基于最大灵敏度指标的稳定分数阶PID参数优化方法,其特征在于,获取分数阶PID控制器的参数稳定域包括:
βn>…>β1>β0≥0,αn>…>α1>α0≥0,其闭环传递函数:
将分数阶PID控制器和被控对象代入闭环传递函数,得到特征多项式:Ψ(s)=sλD(s)+(kpsλ+ki+kdsμ+λ)N(s)e-Ls,对于一组参数(kp,ki,kd,λ,μ)能使特征多项式的值为0,则所述分数阶PID控制器的输入输出稳定。
4.根据权利要求3所述的基于最大灵敏度指标的稳定分数阶PID参数优化方法,其特征在于,所述获得满足最大灵敏度的参数稳定解,包括:
将s=jω,和ejx=cosx+jsinx代入所述特征多项式, 得到一组KP,Ki关于ω的方程,在[ωmin,ωmax]内以一定步长取ω的值,得到KP,Ki的参数稳定解的解集。
5.根据权利要求1所述的基于最大灵敏度指标的稳定分数阶PID参数优化方法,其特征在于,最大灵敏度指标为:MT≤γT,γT≥1,其中, γT为所述预设的最大灵敏度指标。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410395817.3A CN104199293A (zh) | 2014-08-12 | 2014-08-12 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
CN201510475140.9A CN105045093B (zh) | 2014-08-12 | 2015-08-05 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410395817.3A CN104199293A (zh) | 2014-08-12 | 2014-08-12 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN104199293A true CN104199293A (zh) | 2014-12-10 |
Family
ID=52084595
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410395817.3A Pending CN104199293A (zh) | 2014-08-12 | 2014-08-12 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
CN201510475140.9A Expired - Fee Related CN105045093B (zh) | 2014-08-12 | 2015-08-05 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
Family Applications After (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510475140.9A Expired - Fee Related CN105045093B (zh) | 2014-08-12 | 2015-08-05 | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (2) | CN104199293A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106200619A (zh) * | 2016-08-15 | 2016-12-07 | 浙江大学 | 附带控制器输出约束的pi控制回路性能评估方法 |
CN106338913A (zh) * | 2016-11-04 | 2017-01-18 | 河北省科学院应用数学研究所 | 基于相位裕量和截止频率的分数阶pid控制器设计方法 |
CN106681150A (zh) * | 2017-01-19 | 2017-05-17 | 华侨大学 | 基于闭环参考模型的分数阶pid控制器参数优化整定方法 |
CN106802637A (zh) * | 2015-11-26 | 2017-06-06 | 发那科株式会社 | 根据主轴负载控制进给速度的数值控制装置 |
CN108121207A (zh) * | 2017-12-21 | 2018-06-05 | 浙江工业大学 | 同时满足跟踪性和鲁棒性的燃气锅炉最优分数阶pid控制方法 |
CN110119084A (zh) * | 2018-02-03 | 2019-08-13 | 苏州电器科学研究院股份有限公司 | 一种分数阶pid的无轴承电机驱动***参数整定方法 |
CN113156824A (zh) * | 2021-04-29 | 2021-07-23 | 苏州科技大学 | 一种基于神经网络的分数阶内模pid控制器 |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107045280A (zh) * | 2017-03-24 | 2017-08-15 | 清华大学 | 分数阶pid控制器的有理化实现方法 |
CN107657125B (zh) * | 2017-09-30 | 2021-07-20 | 中冶赛迪技术研究中心有限公司 | 基于联合仿真的控制***设计与pid参数整定方法及*** |
CN108536138B (zh) * | 2017-12-29 | 2021-04-20 | 哈尔滨航天恒星数据***科技有限公司 | 农机自动驾驶***的转向控制参数整定方法及*** |
CN110090728B (zh) * | 2019-01-24 | 2020-11-06 | 济南大学 | 用于控制水泥生料立磨中喂料量的方法、装置及设备 |
CN109799698B (zh) * | 2019-01-30 | 2020-07-14 | 上海交通大学 | 时滞视觉伺服***最优pi参数优化方法及*** |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04118702A (ja) * | 1990-09-10 | 1992-04-20 | Toshiba Corp | 制御パラメータ決定方法 |
JPH0713632A (ja) * | 1993-06-25 | 1995-01-17 | Kubota Corp | 水道用バルブの制御用コントローラ |
JP4474555B2 (ja) * | 2007-01-23 | 2010-06-09 | 国立大学法人広島大学 | Pid制御装置及びpid制御方法 |
JP5108554B2 (ja) * | 2008-02-20 | 2012-12-26 | アズビル株式会社 | Pidパラメータ調整装置および調整方法 |
CN103116270A (zh) * | 2013-01-30 | 2013-05-22 | 上海交通大学 | 基于灵敏度约束的稳定PIλDμ控制器参数整定方法 |
CN103454914B (zh) * | 2013-08-23 | 2016-04-13 | 上海交通大学 | 一种多指标约束的分数阶PIλDμ控制器整定方法 |
-
2014
- 2014-08-12 CN CN201410395817.3A patent/CN104199293A/zh active Pending
-
2015
- 2015-08-05 CN CN201510475140.9A patent/CN105045093B/zh not_active Expired - Fee Related
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106802637A (zh) * | 2015-11-26 | 2017-06-06 | 发那科株式会社 | 根据主轴负载控制进给速度的数值控制装置 |
CN106802637B (zh) * | 2015-11-26 | 2019-03-01 | 发那科株式会社 | 根据主轴负载控制进给速度的数值控制装置 |
US10488849B2 (en) | 2015-11-26 | 2019-11-26 | Fanuc Corporation | Numerical controller for controlling feed rate based on spindle load |
CN106200619A (zh) * | 2016-08-15 | 2016-12-07 | 浙江大学 | 附带控制器输出约束的pi控制回路性能评估方法 |
CN106200619B (zh) * | 2016-08-15 | 2018-11-20 | 浙江大学 | 附带控制器输出约束的pi控制回路性能评估方法 |
CN106338913A (zh) * | 2016-11-04 | 2017-01-18 | 河北省科学院应用数学研究所 | 基于相位裕量和截止频率的分数阶pid控制器设计方法 |
CN106681150A (zh) * | 2017-01-19 | 2017-05-17 | 华侨大学 | 基于闭环参考模型的分数阶pid控制器参数优化整定方法 |
CN108121207A (zh) * | 2017-12-21 | 2018-06-05 | 浙江工业大学 | 同时满足跟踪性和鲁棒性的燃气锅炉最优分数阶pid控制方法 |
CN110119084A (zh) * | 2018-02-03 | 2019-08-13 | 苏州电器科学研究院股份有限公司 | 一种分数阶pid的无轴承电机驱动***参数整定方法 |
CN113156824A (zh) * | 2021-04-29 | 2021-07-23 | 苏州科技大学 | 一种基于神经网络的分数阶内模pid控制器 |
Also Published As
Publication number | Publication date |
---|---|
CN105045093B (zh) | 2018-02-09 |
CN105045093A (zh) | 2015-11-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104199293A (zh) | 基于最大灵敏度指标的稳定分数阶pid参数优化方法 | |
Zhao et al. | Control design for the SISO system with the unknown order and the unknown relative degree | |
CN103293956A (zh) | 一种被控对象参数不确定***的分数阶pid控制器整定方法 | |
CN103454914B (zh) | 一种多指标约束的分数阶PIλDμ控制器整定方法 | |
CN104089270A (zh) | 一种发电机组锅炉负荷控制优化调整试验方法 | |
Adnan et al. | Self-tuning fuzzy PID controller for electro-hydraulic cylinder | |
CN104199437A (zh) | 基于区域极点指标的分数阶PIλDμ控制器的参数优化方法 | |
Roy et al. | Design of optimal PIλDδ controller for speed control of DC motor using constrained particle swarm optimization | |
CN102591203B (zh) | 一种伺服电机的基于微分器的直接神经网络控制方法 | |
CN107807522A (zh) | 水下机器人轨迹跟踪反步控制方法 | |
Romasevych et al. | PI-controller tuning optimization via PSO-based technique | |
Outanoute et al. | Synthesis of an Optimal Dynamic Regulator Based on Linear Quadratic Gaussian (LQG) for the Control of the Relative Humidity Under Experimental Greenhouse. | |
Balandin et al. | Design of generalized discrete-time H∞-optimal control over finite and infinite intervals | |
CN105320173A (zh) | 一种基于混沌优化免疫模糊pid的汽提温度控制方法 | |
CN111367178A (zh) | 一种船舶自动舵的自适应控制装置及方法 | |
CN104155876A (zh) | 一种pid控制器的分离实现方法 | |
Mitkowski et al. | Tuning of the half-order robust PID controller dedicated to oriented PV system | |
Marquez-Rubio et al. | Stabilization region of PD controller for unstable first order process with time delay | |
Pahasa et al. | Model predictive control-based wind turbine blade pitch angle control for alleviation of frequency fluctuation in a smart grid | |
Tang et al. | Model reference adaptive control based on Lyapunov stability theory | |
KR20140077720A (ko) | 학습 기반 온실 제어 방법 및 제어 시스템 | |
Arıcan et al. | Suboptimal control of a 3 dof helicopter with state dependent riccati equations | |
Gopinath | Design of fractional order controllers for first order plus time delay systems | |
Díaz-Rodríguez | Modern design of classical controllers: continuous-time first order controllers | |
Pahuja et al. | Design of Robust Control for Pertubed Systems Using QFT |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20141210 |