CN105955029B - 一种保鲁棒性的pid控制参数优化方法 - Google Patents
一种保鲁棒性的pid控制参数优化方法 Download PDFInfo
- Publication number
- CN105955029B CN105955029B CN201610395162.9A CN201610395162A CN105955029B CN 105955029 B CN105955029 B CN 105955029B CN 201610395162 A CN201610395162 A CN 201610395162A CN 105955029 B CN105955029 B CN 105955029B
- Authority
- CN
- China
- Prior art keywords
- particle
- controller
- algorithm
- optimization
- robustness
- 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.)
- Active
Links
- 238000005457 optimization Methods 0.000 title claims abstract description 67
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 68
- 239000002245 particle Substances 0.000 claims abstract description 59
- 238000013461 design Methods 0.000 claims description 19
- 230000004044 response Effects 0.000 claims description 18
- 238000004458 analytical method Methods 0.000 claims description 10
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 3
- 238000003786 synthesis reaction Methods 0.000 claims description 3
- 230000006870 function Effects 0.000 description 41
- 238000012360 testing method Methods 0.000 description 12
- 238000010586 diagram Methods 0.000 description 6
- 230000003044 adaptive effect Effects 0.000 description 3
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 238000013459 approach Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- RZVHIXYEVGDQDX-UHFFFAOYSA-N 9,10-anthraquinone Chemical compound C1=CC=C2C(=O)C3=CC=CC=C3C(=O)C2=C1 RZVHIXYEVGDQDX-UHFFFAOYSA-N 0.000 description 1
- BHELIUBJHYAEDK-OAIUPTLZSA-N Aspoxicillin Chemical compound C1([C@H](C(=O)N[C@@H]2C(N3[C@H](C(C)(C)S[C@@H]32)C(O)=O)=O)NC(=O)[C@H](N)CC(=O)NC)=CC=C(O)C=C1 BHELIUBJHYAEDK-OAIUPTLZSA-N 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000019771 cognition Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000002156 mixing Methods 0.000 description 1
- 238000011056 performance test Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000000452 restraining effect Effects 0.000 description 1
- 238000013432 robust analysis Methods 0.000 description 1
- 238000010845 search algorithm Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000010200 validation analysis 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
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
本发明公开一种保鲁棒性的PID控制参数优化方法,先采用约束粒子群优化算法,将控制器的性能指标和控制器所要满足的约束条件进行综合,形成带约束条件的性能指标函数,从而将约束优化问题转化为一般的非约束优化问题,再采用标准粒子群算法获得满足约束条件的最优PID控制参数。此种优化方法可搜索最优的PID控制器参数,使得所得到的控制器满足μ综合鲁棒性指标,并进一步将该PID控制器参数优化方法应用于导弹PID控制器。
Description
技术领域
本发明属于自动控制技术领域,特别涉及一种保鲁棒性的PID控制参数优化方法。
背景技术
导弹的控制器设计是导弹实现自动攻击目标的关键技术之一,控制器设计的好坏直接关系到导弹的飞行品质和脱靶量指标,是导弹飞行控制系统设计的核心内容。目前的导弹机绝大多数采用PID控制器,随着导弹控制技术的发展,如第五代空对空格斗导弹等,一个日益突出的要求是需要所设计的控制器满足多变量鲁棒性指标。
传统的导弹PID控制律设计方法是单通道设计多通道验证的方法,并且所采用的单通道设计方法只能采用幅值裕度、相角裕度等传统单通道鲁棒性指标,且传统鲁棒方法设计的控制器存在阶数过高难以实现的弱点,而经典PID控制方法实现简单但又无法在设计过程中保证多变量系统的鲁棒性能。近年来,采用非线性优化算法来解决控制器优化设计问题的研究已逐步进入应用阶段。
发明内容
本发明的目的,在于提供一种保鲁棒性的PID控制参数优化方法,其可搜索最优的PID控制器参数,使得所得到的控制器满足μ综合鲁棒性指标,并进一步将该PID控制器参数优化方法应用于导弹PID控制器。
为了达成上述目的,本发明的解决方案是:
一种保鲁棒性的PID控制参数优化方法,先采用约束粒子群优化算法,将控制器的性能指标和控制器所要满足的约束条件进行综合,形成带约束条件的性能指标函数,从而将约束优化问题转化为一般的非约束优化问题,再采用标准粒子群算法获得满足约束条件的最优PID控制参数。
上述优化方法具体包括如下步骤:
步骤A、确定粒子种群数量N和优化最大代数kmax;
步骤B、在粒子的取值范围内初始化粒子种群每一个粒子的位置和速度;
步骤C、采用μ分析方法计算鲁棒性指标,同时计算时域性能指标作为约束条件并进行性能指标函数的综合;
步骤D、判断当前最优值是否满足性能指标停止条件或达到最大优化代数,如满足二者之一即停止算法并转到步骤F;
步骤E、更新每个粒子的位置并转到步骤C;
步骤F、输出最优粒子的位置和对应的性能指标函数值。
上述步骤A中,通过下式将约束优化问题转化为非约束优化问题:
其中fm(x)定义为:
其中,f(x)为原目标函数,hmax(x):=max([h1(x),h2(x),....,hn(x)])表示当前所有约束条件中违反量最大的约束条件。
上述优化方法应用的μ-PID控制器设计描述为:给定目标函数
搜索其最小值并满足由时域响应指标组成的约束条件;其中u=K(s,x)y表示控制器,x∈Rm为控制器参数,s为复数域自变量。
上述μ-PID控制器的样例导弹在15000米高度、2.8马赫、40度攻角下的特征工作点纵向通道状态空间模型,其状态方程为:
输出方程为:
式中:u、w分别为纵向和法向速度,α为攻角,q为俯仰角速率,Az为法向加速度,δe为等效升降舵偏角。
采用上述方案后,本发明利用一种改进的粒子群优化算法—简便约束粒子群优化算法搜索最优的PID控制器参数,并进一步将该PID控制器参数优化方法应用于导弹的保鲁棒性PID控制器。所述简便约束粒子群优化算法的基本思想是采用目标函数替换的方法将约束优化问题转化为非约束优化问题,具有简便易用的优点。利用本发明PID控制器参数优化方法对导弹PID控制器进行优化,所得到的导弹PID控制器能够在满足传统时域指标的同时使得多变量鲁棒性指标---μ值最小,即具有最佳的鲁棒性能。。
相比现有技术,本发明具有以下有益效果:
(1)所得到的控制器在满足时域指标的同时可以保证鲁棒性指标—μ值最小;
(2)对约束条件进行综合时没有引入任何新的算法参数;
(3)相比于现有约束PSO算法,本发明的综合算法更加简单;
(4)对原目标函数没有连续可微等要求;
(5)具有优良的搜索性能。
附图说明
图1是标准带不确定性反馈控制系统结构示意图;
图2是本发明中进行鲁棒性分析是采用的μ分析结构示意图;
图3是导弹Raytheon驾驶仪的俯仰/偏航通道PID控制器结构示意图;
图4是采用罚函数法PSO算法对样例导弹纵向通道PID控制器参数进行优化时的性能指标函数收敛性曲线;
图5是采用ALPSO算法对样例导弹纵向通道PID控制器参数进行优化时的性能指标函数收敛性曲线;
图6是采用本发明对样例导弹纵向通道PID控制器参数进行优化时的性能指标函数收敛性曲线;
图7是对本发明结果进行验证的样例导弹闭环系统仿真结构示意图;
图8是获得的PID控制器鲁棒特性曲线图;
图9是获得的PID控制器对标称系统的阶跃响应图;
图10是获得的PID控制器对摄动系统的阶跃响应图。
具体实施方式
以下将结合附图,对本发明的技术方案进行详细说明。
本发明提供一种保鲁棒性的PID控制参数优化方法,它是先采用约束粒子群优化算法,将控制器的性能指标和控制器所要满足的约束条件进行综合,形成新的带约束条件的性能指标函数,从而将约束优化问题转化为一般的非约束优化问题,再采用标准粒子群算法获得满足约束条件的最优PID控制参数,所述优化方法具体包括如下步骤:
步骤A、确定粒子种群数量N和优化最大代数kmax;
步骤B、在粒子的取值范围内初始化粒子种群每一个粒子的位置和速度;
步骤C、采用μ分析方法计算鲁棒性指标,同时计算时域性能指标作为约束条件并进行性能指标函数的综合;
步骤D、判断当前最优值是否满足性能指标停止条件或达到最大优化代数,如满足二者之一即停止算法并转到步骤F;
步骤E、更新每个粒子的位置并转到步骤C;
步骤F、输出最优粒子的位置和对应的性能指标函数值。
以下将对本发明中的几个技术特征进行详细说明。
1、PSO算法
粒子群优化算法是由美国社会心理学家James Kennedy和电气工程师RussellEberhart在1995年共同提出的一种进化计算技术,它是一种通过模拟鸟类群体行为进行建模与仿真研究而发展起来的一种基于群体协作的随机搜索算法,属于群体智能搜索算法的一种,也称为微粒群算法。该算法原理是从一组随机的初始值出发,采用基于邻域的搜索技术,通过迭代来寻找最优解,在每一次迭代中,粒子通过个体极值(粒子本身找到的最优解)和群体极值(种群目前找到的最优解)来更新自己。由于粒子群优化算法其算法本身结构简单,并在多种复杂非凸优化问题的应用中表现优异,近年来引起了越来越多的关注。
1.1标准PSO算法
在PSO算法中,每一个优化问题的解作为一只单个的鸟,或称为粒子,粒子的属性包括其位置x和速度v。每个粒子在多维空间中依据个体经验和群体中其他粒子的经验来动态调整自身移动的轨迹与速度。在第k次迭代中,第i个粒子xi=(xi,1,xi,2,...,xi,m)∈Rm依据以下公式对位置和速度进行调整:
其中,为粒子i在第k次迭代时第d维位置,xmax,d,xmin,d为粒子在第d维的最远和最近位置,粒子的运动被设定为不能超过此边界;为粒子i在第k次迭代时第d维速度,vmax,d为粒子在第d维的最大速度,粒子的运动被设定为不能超过此最大速度;为粒子i个体当前在第d维位置的最优值;为粒子群体当前在第d维位置的最优值;c1为粒子个体认知加速常数,c2为群体认知加速常数,参数c1和c2分别代表了个体自身行为和群体行为对个体影响的大小;是介于[0,1]之间的随机数;ω为惯性系数,代表了粒子的运动惯性,按照运动惯性随进化代数逐渐减小的原则,其表达式为:
式中:ωmax和ωmin分别为ω的最大值与最小值;kmax为最大进化代数。
1.2本发明涉及的约束PSO算法
粒子群优化算法最初是作为一种非约束优化算法被提出的,然而在大量的工程实践中,有相当多的应用属于要求满足一定不等式约束条件下的优化问题。典型的约束PSO算法有自适应惩罚函数粒子群优化算法和ALPSO算法,然而这两种算法均有一些如前所述不足之处。为了克服这些不足,本发明提出了一种新型简便易用的约束PSO算法。
一般约束优化问题可以描述为:
其中,函数h(x):=[h1(x),h2(x),....,hn(x)],(Rm→Rn)表示约束条件,F表示所有可行的区域,在本发明中,假设F为非空集合。
本发明采用如下一种新颖简便的方法将式(6)所描述的约束优化问题转化为式(7)所示的非约束优化问题,在该过程中不需要增加任何新的算法参数。
其中fm(x)定义为:
其中,f(x)为原目标函数(性能指标函数),hmax(x):=max([h1(x),h2(x),....,hn(x)])表示当前所有约束条件中违反量最大的约束条件。在自适应罚函数法中,每一步都需要考虑所有的约束条件,显然是没有必要的。取arctan(f(x))的形式可以使得函数在(-∞,+∞)上为单调增函数,同时,当粒子未处于可行域时(约束条件不满足的区域),即hmax(x)≥0时,考虑融合约束条件时采用指数函数的形式更加有利于在算法的初始阶段扩大搜索范围而在算法末段通过细致搜索从而增大得到最优值的概率。为了使得粒子能够进入可行域,即当hmax(xi k+1)<hmax(xi k)时,能满足fm(xi k+1)<fm(xi k),因此这里把底数取为另外,当处于可行域时,即hmax(x)<0,只考虑原目标函数,且减去使优化算法保持以下特性:当依据式(1)-式(4)进化到下一代,即使此时出现的情形,因为fm(xi k)<0<1≤fm(xi k+1),与仍将处于可行域中。可见,依据上述原理,所有违反约束条件的粒子都具有自动回到可行域的趋势。
显然,当式(8)所示的fm(x)取得最优值时,原目标函数f(x)也取得最优值,且所有的约束条件均得到满足。此时,原约束优化问题就转化为了非约束优化问题,可以采用1.1节的标准PSO算法进行寻优。
本文提出的算法最大的优势在于其简便性,它不需要像ALPSO算法那样要求目标函数连续可微,也不需要像自适应罚函数法那样增加新的需调参数。在后续内容中,该算法将与ALPSO算法、罚函数法进行性能比较。
2、算法性能测试
目前,非线性优化算法的优劣尚不能从理论上进行证明,通常只能通过测试函数的实测结果来进行比较。为了验证所提出的约束粒子群优化算法的性能,在本部分内容中采用了多个标准测试函数来对其性能进行全面测试与验证,其中包含24个标准测试函数,我们选择了其中带不等式约束的单目标优化问题共12个测试函数,这些测试函数在约束优化算法的性能测试中被广泛采用。
表1三种约束优化算法测试结果
测试程序在Matlab7.8环境下运行,采用一台Intel i5CPU PC机,内存为4G字节。每个任务运行30次,算法停止的条件是达到最大进化代数1000。同时为了与前述自适应罚函数法PSO算法、ALPSO算法进行对比,对这两种算法也进行了同样的测试。在进行比较分析时,粒子群算法参数取为:
●粒子数量:30
●粒子维数:依测试函数
●粒子最大移动速度:Vmax,d=xmax,d/2
●学习因子:c1=2,c2=2
●惯性系数最大与最小值:Wmax=0.9,Wmin=0.4
●最大搜索代数:1000
测试结果如表1所示,显然,本发明提出的简便约束PSO算法在处理带约束的优化问题中具有优越的性能,对测试函数寻优的成功率明显高于自适应罚函数PSO算法和ALPSO算法,并且得到的最终优化结果及其分布特性也优于上述两种算法。
3、导弹μ-PID控制器设计
3.1理论分析
结构奇异值μ(Structured Singular Value)是一种线性代数工具,用来定量地表征结构化不确定性对线性动态系统稳定性及性能的影响。相比于H∞方法,μ综合方法可以同时分析控制系统的鲁棒稳定性和鲁棒性能。
3.1.1μ分析方法
对任意多输入不确定线性闭环系统,可以从图1看出。其中M代表系统传递函数矩阵,由控制器和控制对象构成。
代表模型不确定性,即作为控制对象的数学模型与实际对象存在的差异,且有:
其中, 表示Δ的最大奇异值。
若系统不确定性Δ满足式(10),则系统传递函数矩阵M的结构奇异值定义为:
直接按照式(11)来求解结构奇异值μ是非常困难的,通常采用的是一种逼近方法。首先计算μ值的上下界,当上下界之差足够小时,取其为近似值,该方法称为“D-K”迭代法,由Doyle在1985年提出。
3.1.2μ-PID控制器设计问题
μ分析问题可由图2表示,在此结构中,M代表广义控制对象以及控制器,w1表示外部输入信号,z2和w2分别表示模型不确定性的输入和输出信号,z1表示控制输出信号。此时系统∑(s,x)可分解为如下形式:
则在本发明中,基于约束优化的μ-PID控制器设计可以描述为:给定目标函数
搜索其最小值并满足由时域响应指标组成的约束条件。其中式(13)中的K(s,x)表示控制器,x∈Rm为控制器参数,s为复数域自变量。
3.2导弹控制设计与仿真
样例导弹采用基于多工作点线性化模型的控制律设计,六自由度非线性模型验证的方法,因此本发明中控制律设计仍以线性化模型为设计对象。以样例导弹在15000米高度、2.8马赫、40度攻角下的特征工作点纵向通道状态空间模型为例,其状态方程为:
输出方程为:
式中:u、w分别为纵向和法向速度,α为攻角,q为俯仰角速率,Az为法向加速度,δe为等效升降舵偏角。
导弹的定结构控制器可以有各种形式,如姿态控制器、迎角控制器和过载控制器等,每种控制器依据控制对象本身的特性和控制需求还可以细分为更多的回路控制结构。不失一般性,本文以经典Raytheon驾驶仪控制结构为例,其纵向控制结构如图(3)所示。
样例导弹纵向通道所需满足的时域性能指标如表2所示:
表2纵向通道阶跃响应性能指标
由于控制参数均为正值,为了缩小搜索空间的范围,将待寻优的控制参数转换为:
(xe1 xe2 xe3 xe4)=(log10iy log10k5 log10k4 log10k0)
此时,根据手工设计的经验,待寻优参数的范围可以设置为:
Δfe:={(xe1,xe2,xe3,xe4)∈R4:-2<xei<2,i=1,2,3,4},样例导弹定结构鲁棒控制问题即为在满足上述阶跃响应时域性能指标的约束条件下寻找优化参数使得目标函数μ取得最小值。
由于导弹的μ-PID控制器的设计问题相比于测试函数而言要复杂得多,计算量也大很多,依据多轮设计与仿真的经验,此时简便约束PSO算法的参数可以取为:
●粒子数量:20
●粒子维数:4
●粒子最大移动速度:Vmax,d=xmax,d/2i=1...4,
●学习因子:c1=2,c2=2
●惯性系数最大与最小值:Wmax=0.9,Wmin=0.4
●最大搜索代数:20
●每个任务运行次数:50
●终止条件:达到最大代数
采用如上所述性能指标函数和约束条件,基于三种不同PSO算法的优化结果如表3-表5所示:
表3罚函数PSO算法优化结果
表4ALPSO算法优化结果
表5本发明算法优化结果
由表3-表5可知,罚函数PSO算法所得到的三种μ值均最大,μ最小值大于1,表明该优化方法所得到的控制器鲁棒性较差;ALPSO算法总体性能次之;本发明提出的简便约束PSO算法显然具有最佳的性能,并且传统的幅值和相角裕度指标也验证了这一结果。基于上述三种PSO算法的目标函数收敛性曲线如图4-图6所示。
由图4-图6明显可见,相比于罚函数PSO算法和ALPSO算法,本发明提出的简便约束PSO算法在设计样例导弹纵向通道控制律过程中具有最好的统计收敛特性。
在简便约束PSO算法得到的50次结果中,选取最优值结果:iy=0.97,K5=0.2,K4=2.44,K0=0.41作为控制参数,并且在图7所示的俯仰通道闭环系统中,模型不确定性函数和噪声干扰的大小分别取为:模型不确定性权函数Win:2·(s+3.2)/(s+160);加速度计噪声:1%;角速率陀螺噪声:0.1%;
对闭环控制系统的频域和时域性能分析如下:
(1)频域性能分析
a.鲁棒稳定性分析
图7所示的俯仰通道闭环系统中,令外输入为0,仅考虑输入端乘型不确定性Δin作用下z1的响应。图8所示从d1到z1的传递函数的结构奇异值上界曲线。由图中可见结构奇异值均在0.3以下,表明系统具有良好的鲁棒稳定性。
b.鲁棒性能分析
考察在外输入和输入端乘型不确定性同时作用下系统的性能。如图8下图所示,结构奇异值均在1以下,表明系统的鲁棒性能也满足指标要求。
(2)时域性能分析
a.标称系统单位阶跃响应
图7所示的俯仰通道闭环系统中,令Δin=0且扰动d2=0,仅考虑俯仰通道过载命令Azg=1作用下系统的阶跃响应。从图9可见,系统的标称性能满足表2所示的时域性能指标。
b.摄动系统单位阶跃响应
在外输入和输入端乘型不确定性同时作用下系统的阶跃响应如图10所示,可知系统在摄动下的性能仍然满足表2所示的时域性能指标,从而印证了频域分析的结果。
根据上述实例可以看出本发明所提出的简便约束粒子群优化算法用于PID控制器参数优化的优越性,本发明方法尤其适合于导弹鲁棒PID控制器设计,使得所设计的无人机PID控制律能够在满足传统时域指标的同时能保证多变量鲁棒性指标---μ值最小,即具有最佳的鲁棒性能。
以上实施例仅为说明本发明的技术思想,不能以此限定本发明的保护范围,凡是按照本发明提出的技术思想,在技术方案基础上所做的任何改动,均落入本发明保护范围之内。
Claims (4)
1.一种保鲁棒性的PID控制参数优化方法,其特征在于:先采用约束粒子群优化算法,将控制器的性能指标和控制器所要满足的约束条件进行综合,形成带约束条件的性能指标函数,从而将约束优化问题转化为一般的非约束优化问题,再采用标准粒子群算法获得满足约束条件的最优PID控制参数;
所述优化方法具体包括如下步骤:
步骤A、确定粒子种群数量N和优化最大代数kmax;
步骤B、在粒子的取值范围内初始化粒子种群每一个粒子的位置和速度;
步骤C、采用μ分析方法计算鲁棒性指标,同时计算时域性能指标作为约束条件并进行性能指标函数的综合;
步骤D、判断当前最优值是否满足性能指标停止条件或达到最大优化代数,如满足二者之一即停止算法并转到步骤F;
步骤E、更新每个粒子的位置并转到步骤C;
步骤F、输出最优粒子的位置和对应的性能指标函数值。
2.如权利要求1所述的一种保鲁棒性的PID控制参数优化方法,其特征在于:所述步骤A中,通过下式将约束优化问题转化为非约束优化问题:
其中fm(x)定义为:
其中,f(x)为原目标函数,hmax(x):=max([h1(x),h2(x),....,hn(x)])表示当前所有约束条件中违反量最大的约束条件。
3.如权利要求1所述的一种保鲁棒性的PID控制参数优化方法,其特征在于:所述优化方法应用的μ-PID控制器设计描述为:给定目标函数
搜索其最小值并满足由时域性能指标组成的约束条件;其中μ=K(s,x)表示控制器,x∈Rm为控制器参数,s为复数域自变量。
4.如权利要求3所述的一种保鲁棒性的PID控制参数优化方法,其特征在于:所述μ-PID控制器的样例导弹在15000米高度、2.8马赫、40度攻角下的特征工作点纵向通道状态空间模型,其状态方程为:
输出方程为:
式中:u、w分别为纵向和法向速度,α为攻角,q为俯仰角速率,Az为法向加速度,δe为等效升降舵偏角。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610395162.9A CN105955029B (zh) | 2016-06-06 | 2016-06-06 | 一种保鲁棒性的pid控制参数优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610395162.9A CN105955029B (zh) | 2016-06-06 | 2016-06-06 | 一种保鲁棒性的pid控制参数优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105955029A CN105955029A (zh) | 2016-09-21 |
CN105955029B true CN105955029B (zh) | 2019-03-29 |
Family
ID=56907772
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610395162.9A Active CN105955029B (zh) | 2016-06-06 | 2016-06-06 | 一种保鲁棒性的pid控制参数优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105955029B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105947184B (zh) * | 2016-05-31 | 2018-07-10 | 中国航空工业集团公司西安飞机设计研究所 | 一种迎角控制器的稳定性计算方法 |
CN108345216B (zh) * | 2018-01-12 | 2021-10-26 | 中国科学院理化技术研究所 | 一种基于多目标粒子群算法的磁悬浮轴承鲁棒控制器构造方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551642A (zh) * | 2009-04-08 | 2009-10-07 | 南京航空航天大学 | 用于无人飞机控制律参数自动优化的改进粒子群算法 |
CN103645636A (zh) * | 2013-11-25 | 2014-03-19 | 南京航空航天大学 | 一种pid控制器参数优化方法 |
CN103972908A (zh) * | 2014-05-23 | 2014-08-06 | 国家电网公司 | 一种基于自适应混沌粒子群算法的多目标无功优化方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8606424B2 (en) * | 2011-04-05 | 2013-12-10 | King Fahd University Of Petroleum And Minerals | Particle swarm optimization system and method for microgrids |
-
2016
- 2016-06-06 CN CN201610395162.9A patent/CN105955029B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551642A (zh) * | 2009-04-08 | 2009-10-07 | 南京航空航天大学 | 用于无人飞机控制律参数自动优化的改进粒子群算法 |
CN103645636A (zh) * | 2013-11-25 | 2014-03-19 | 南京航空航天大学 | 一种pid控制器参数优化方法 |
CN103972908A (zh) * | 2014-05-23 | 2014-08-06 | 国家电网公司 | 一种基于自适应混沌粒子群算法的多目标无功优化方法 |
Non-Patent Citations (1)
Title |
---|
基于约束粒子群优化的导弹H∞-PID控制器设计;张民 等;《弹箭与制导学报》;20160228;第36卷(第1期);第1-4、9页 |
Also Published As
Publication number | Publication date |
---|---|
CN105955029A (zh) | 2016-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103645636B (zh) | 一种pid控制器参数优化方法 | |
Gao et al. | Stabilization of nonlinear systems under variable sampling: a fuzzy control approach | |
Tao et al. | Multiple model predictive control for large envelope flight of hypersonic vehicle systems | |
CN110376882A (zh) | 基于有限时间扩张状态观测器的预定性能控制方法 | |
Wang et al. | A restart univariate estimation of distribution algorithm: sampling under mixed Gaussian and Lévy probability distribution | |
CN105182743A (zh) | 一种基于鲁棒h无穷的变增益解耦控制方法 | |
CN106502255B (zh) | 一种舰载机自动着舰控制系统的设计方法和控制方法 | |
Rahimi et al. | Designing and simulation for vertical moving control of UAV system using PID, LQR and Fuzzy Logic | |
CN105955029B (zh) | 一种保鲁棒性的pid控制参数优化方法 | |
CN112819303A (zh) | 基于pce代理模型的飞行器追踪效能评估方法及系统 | |
CN109085848A (zh) | 空空导弹直接力/气动力有限时间抗饱和控制方法 | |
CN102566446B (zh) | 基于线性模型组的无人直升机全包线数学模型构建方法 | |
CN109543284B (zh) | 基于Kriging空间插值的火星大气进入段最优制导方法 | |
Zhang et al. | Motion state recognition and trajectory prediction of hypersonic glide vehicle based on deep learning | |
CN114611416A (zh) | 导弹非线性非定常气动特性ls-svm建模方法 | |
Yuan-bin | Genetic algorithm’s application for optimization of PID parameters in dynamic positioning vessel | |
CN109033487A (zh) | 一种基于蒙特卡洛仿真的飞行器总体概率设计方法 | |
CN106382853B (zh) | 一种带终端弹道倾角和攻角约束奇异摄动次优制导律 | |
CN107194155A (zh) | 一种基于小数据集和贝叶斯网络的威胁评估建模方法 | |
Zheng et al. | An effective genetic algorithm to VDA with discontinuous “on-off” switches | |
Xu et al. | Parameter selection for particle swarm optimization based on stochastic multi-objective optimization | |
Xiao-ping | Aircraft evasive maneuver trajectory optimization based on QPSO | |
Na et al. | Modified particle swarm optimization based algorithm for BP neural network for measuring aircraft remaining fuel volume | |
Gan et al. | Wavelet neural network aerodynamic modeling from flight data based on pso algorithm with information sharing and velocity disturbance | |
Wang et al. | Research on performance of artificial bee colony algorithm based on benchmark test function |
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 |