, 王继东1,2, 刘秀翀1, 孙鉴1 1. 东北大学 信息科学与工程学院, 辽宁 沈阳 110819;
2. 华北水利水电大学 电力学院, 河南 郑州 450045
收稿日期：2018-08-21
基金项目：国家自然科学基金资助项目(61473070, 61433004, 61627809)。
作者简介：王占山(1971-), 男, 辽宁抚顺人, 东北大学教授, 博士生导师。
摘要：针对具有复特征根的二阶线性系统的鲁棒H∞控制性能进行了研究, 特别是探讨了虚部特征根对H∞性能的作用, 揭示了当特征根的实部固定不变的情况下, 特征根的虚部增加到一定程度就会对H∞性能不起作用, 进而针对给定的系统矩阵, 在H∞性能的框架下分析了特征根实部和虚部之间的比例关系.同时, 比照经典二阶系统的部分成果, 探讨了经典二阶系统的阻尼比系数选取范围与H∞性能框架下复特征根参数之间的对照关系, 进而从时域的角度给出了阻尼比参数选取的一种合理解释, 并形成时频交互综合分析, 为极点配置控制理论中的期望极点的选取提供一种理论参考.
关键词：H∞范数H∞性能H∞控制干扰阻尼比时频交互综合分析
Impact of Complex Eigenvalues on Their Robust H∞ Control Performance in Second-Order Systems
WANG Zhan-shan1
, WANG Ji-dong1,2, LIU Xiu-chong1, SUN Jian1 1. School of Information Science & Engineering, Northeastern University, Shenyang 110819, China;
2. School of Electrical Engineering, North China University of Water Resources and Electric Power, Zhengzhou 450045, China
Corresponding author: WANG Zhan-shan, E-mail:wangzhanshan@ise.neu.edu.cn
Abstract: Robust H∞ performance is studied for a class of second-order linear systems with complex eigenvalues, especially the impact of imaginary part of complex eigenvalues on the H∞ performance is studied. It is found that when the real part of the complex eigenvalues is fixed, the impact of the imaginary part of the complex eigenvalues on the H∞ performance becomes little when the value of the imaginary part exceeds a threshold. Based on this observation, a fixed proportional relation between the imaginary part and the real part of the complex eigenvalues is found for a given system matrix under the framework of H∞ performance. Referring to the classical control theory of second-order systems in frequency domain, the relationship between the selection ranges of damping ratio and the complex eigenvalue is further discussed in the framework of H∞ performance, in which a reasonable explanation on the selection ranges of damping ratio is presented. This analysis procedure forms the time-frequency domain interaction synthesis method, which provides a theoretical guidance for the selection and determination of the desired closed-loop poles in the control theory based on the pole placement method.
Key words: H∞ normH∞ performanceH∞ controldisturbancedamping ratiotime-frequency domain interaction synthesis method
H∞控制理论主要是用来解决系统存在干扰和不确定情况下的优化控制问题, 进而实现系统外界干扰或不确定性对系统的期望输出的影响达到最小.H∞方法最初是从频域入手进行的研究, 并在1989年转换到时域, 进而开启了基于状态空间方程描述的受控系统的鲁棒H∞控制研究[1-6].H∞控制实现的是静态向量或矩阵的H∞范数的动态化, 并能确保动态系统对干扰具有一定的抑制性能.围绕如何设计H∞控制器及其相关优化求解方法成为控制界研究的重点内容之一[7-9].实际上, 从不同分析观点来看, H∞控制既属于鲁棒控制, 也属于优化控制范畴.不论属于哪类控制模式, 最终影响系统性态的还是闭环系统的特征根的分布.这样, 实现H∞性能的控制器是满足一定性能指标的集合, 而特征值也将在一定的域内变动.H∞控制器的设计方法, 大体上可以分为基于极点配置的控制方法和基于有界实引理的优化方法[5-6].基于极点配置方法是H∞控制早期的主要工作, 它是从系统结构配置上解决系统抗扰性问题, 但这类研究只是针对期望的给定极点进行的控制器设计, 并没有具体考虑期望极点如何选取以及期望极点的实部和虚部对H∞性能的作用[10-11].另一种是根据干扰到期望输出的H∞范数最小的设计原则, 基于有界实引理方法实现H∞控制器的设计[12-14].这类控制器的设计是基于优化原理进行的, 并由此得到控制器存在的充分条件.这类方法的特点是直接对闭环系统的系统矩阵进行分析, 并求得满足期望性能的一组控制器, 但没有考虑H∞控制器作用后是怎样改变闭环系统的特征值分布的, 以及系统矩阵的特征值、特别是复特征值的实部和虚部对H∞性能的影响如何.作为研究系统特性最基本的特征值法, 应该在H∞控制中体现特征值不同分量的作用, 进而为H∞控制器的设计, 特别是低能耗、低阶控制器的设计提供评价和设计标准.
一般来讲, 系统的实部和虚部都很大的话, 从能量角度来看, 都需要很大的投入.控制器装置作为能量系统的一部分, 相应也会消耗大量的能量, 进而在达到控制效果的前提下, 不论减少实部还是虚部的大小都将是有益的[15-16].具体减少哪部分作用, 这首先需要由理论科技工作者给出解答, 并再由工程技术人员结合具体的问题进行具体运用.
二阶系统是控制系统中应用最广泛、最具代表性的系统之一.许多高阶系统在一定的条件下(如在工作点附件采用的线性化分析方法)常常近似地作为二阶系统来研究, 二阶系统的分析方法也是分析高阶系统的基础[17].同时鉴于在频域内对二阶系统已经有着成熟的研究结果, 如何实现频域和时域分析方法相对照相补充并形成时频交互综合分析, 特别是在鲁棒干扰抑制方面的认识进行统一, 这方面的研究还不多见.进而, 本文针对一类线性二阶系统, 首次从系统特征根的实部和虚部角度来研究H∞控制是通过怎样的系统特征值的配置来达到干扰抑制的目的, 并由此建立与经典频域下的二阶系统理论部分成果之间的相互联系.
本文一是揭示了H∞性能下的实部与虚部之间的比例关系, 说明了在实部固定情况下, 过大增加虚部值对系统性能提升没有任何作用; 二是探究了H∞性能下的二阶系统的阻尼比问题, 验证了经典控制理论中阻尼比选在0.4~0.8的合理性, 并由此为基于极点配置的控制理论中如何确定期望闭环极点提供了基于H∞性能的时频交互综合分析方法.
1 H∞控制基础H∞控制理论是20世纪80年代开始兴起的一门“新”的现代控制理论, 主要研究抑制干扰和不确定性问题, 也就是针对一个给定的性能指标, 使得期望输出对干扰的H∞范数小于一个极小值, 进而抑制干扰对期望输出的影响[1-6].H∞控制问题的基本结构如图 1所示.
图 1(Fig. 1)
 | 图 1 H∞控制基本框图Fig.1 Basic block diagram of H∞ control |
其中:w是干扰; y是系统输出; u是控制输入; x是系统期望输出; G(s)是系统结构矩阵或者是传递函数矩阵; K(s)是反馈控制增益.控制问题就是设计适当的控制器, 使得干扰w对期望输出x的影响最小.
H∞控制的理论基础是H∞范数.向量的H∞范数是指向量诸元素中最大的那个元素, 传递函数的H∞范数是指输入到输出传递函数的波特图中增益的最大值, 反映了输出信号与输入信号能量之比的最大值.如果输入是外界干扰, 显然希望能够通过加入H∞控制器使得输出对输入的放大倍数(传递函数的H∞范数)越小越好[1-6].传递函数的H∞范数刻画了输入到期望输出影响的一种量度, 称之为H∞性能, 它是矩阵H∞范数概念在控制系统理论中的一种延拓.静态向量或矩阵的H∞范数与动态系统的H∞性能具有异曲同工之处, 进而静态H∞范数的动态化是H∞控制的实现基础.由于难以精确求解H∞范数, 进而逼近论方法被提出来, 并结合近似寻优思想, 在干扰抑制方面界实原理被提出来[12-14], 并将频域空间解析求解H∞问题转化为时域状态空间的约束优化问题, 即将求解频域H∞范数的问题转化为求解时域内某期望输出满足‖ x(t)‖2≤γ‖w(t)‖2条件下的界实引理的数值优化问题, 并由此求得最小正标量γ, 即H∞范数.
在图 1中, 如果将传递函数G(s)看作是一个孤立的系统, 则w是外部干扰, u是外部控制输入, K(s)是反馈增益矩阵.如果将传递函数G(s)和反馈增益矩阵K(s)看作是一个整体, 则w是内部干扰, u是内部调节输入.不论是内部反馈控制还是外部反馈控制, 最终的表征都是改变系统的自身结构特征, 即特征值的分布特性.状态空间方法更能够反映系统结构特性, 特别是特征根的实部和虚部分布问题.这样, 为简明起见, 上述H∞控制问题可以用如下的线性状态空间来近似描述：
 | (1) |
2 问题描述考虑二阶系统(1), 现重写如下:
 | (2) |
所考虑的问题是, 在干扰w的作用下, 如何获得状态x的最大暂态上界.在数学意义上, 这一问题等价为具有状态输出的系统(1)的H∞性能估计问题[12-13].
为了凸显系统特征值对于H∞性能估计的作用, 考虑到系统矩阵具有复特征值, 则依据矩阵理论, 存在可逆矩阵H, 使得D=－HAH－1, 其中
 |
 | (3) |
 | (4) |
需说明的是, 系统矩阵A的特征值是与参数α值一一对应的, 尽管γ和Γ正交矩阵的选取不同, 但是不会改变系统矩阵A的表示形式, 进而不会影响系统的H∞性能.相反, 选取不同的参数α值, 具有相同特征值的系统矩阵A存在多种表示形式, 但此时的系统H∞性能可能会不一样.因此, 针对给定的系统矩阵A, 本文这里的变换不会影响系统的H∞性能.
3 仿真算例下面针对系统(3), 利用界实引理(4)来分析不同情况下的系统特征根的实部和虚部对H∞性能的影响.针对系统特征根的实部λ固定, 虚部v从0逐渐增大的情况, 系统的H∞性能以及稳态值如表 1, 表 2所示.当α选取不同值时的H∞性能及其阻尼比, 见表 3.
表 1(Table 1)
 表 1 H∞范数估计(α=2, λ=2)Table 1 Estimation of H∞ norm(α=2, λ=2)
| 表 1 H∞范数估计(α=2, λ=2) Table 1 Estimation of H∞ norm(α=2, λ=2) | |||||||||||||||||||||||||||||||||||||||||||||
表 2(Table 2)
表 2 H∞范数估计(α=2, λ=4)Table 2 Estimation of H∞ norm(α=2, λ=4)
| 表 2 H∞范数估计(α=2, λ=4) Table 2 Estimation of H∞ norm(α=2, λ=4) | |||||||||||||||||||||||||||||||||||||||||||||
表 3(Table 3)
表 3 实部固定时H∞范数对应的虚部实部之最优比例值及相应的阻尼比Table 3 Optimal proportion value and damping ratio of the imaginary part and the real part corresponding to H∞ norm when the real part is fixed
| 表 3 实部固定时H∞范数对应的虚部实部之最优比例值及相应的阻尼比 Table 3 Optimal proportion value and damping ratio of the imaginary part and the real part corresponding to H∞ norm when the real part is fixed |
针对稳定系统矩阵A具有复特征根的情况, 根据自动控制理论[16-17]可知, 实部和虚部的作用各不相同, 实部增大, 稳定裕度增大; 虚部增大, 则振荡环节作用增强, 进而调节时间相应增加.实部和虚部的适当配合可实现快速稳定的动态响应, 不能仅仅强调系统的稳定性能.从表 1, 表 2可见, 当α和实部都给定的情况下, 系统的H∞性能随着虚部不断变化.但是当实部和虚部的比值达到某个界限时, H∞性能就保持不变了.这给控制理论研究提供如下启示:
1) 鲁棒控制器的介入能保证系统稳定的同时, 也能满足一定的鲁棒性能, 这一鲁棒性能是通过改变系统特征值的分布来实现的.适当的虚部值可以对鲁棒性能起到改善作用, 过大的虚部值将对H∞性能没有任何贡献.
2) 随着实部的增大, H∞性能或者H∞范数逐渐减小, 进而对干扰的抑制能力越来越强.这是与系统的稳定裕度增大相适应的, 越稳定的系统抗扰性越强.
3) 将受控系统当作发电系统, 从平均功率理论来看, 虚部往往代表的是无功出力, 进而增大虚部作用, 意味着无功功率增加, 系统进相深度增加, 不利于系统性能提升.将受控系统当作用电系统, 对于恒定的视在功率, 增加虚部无功的作用意味着用电系统的功率因素将降低, 不利于系统的节能降耗.因此, 从消耗能量角度来看, 鲁棒H∞控制器作用后的闭环系统的特征根的虚部值不宜过大.
4) 在频域的极点配置中, 期望性能极点的选取可参照时域H∞性能下的复特征根的分布情况, 形成时频交互综合分析, 进而为极点配置理论提供借鉴和补充作用.即针对给定的极点实部先分析其H∞性能指标, 由此来确定虚部的位置, 并考虑到抗扰性或阻尼比的权衡, 进一步确定虚部值, 并由此确定期望的闭环极点位置.这一时频交互综合分析, 将对传统控制理论中的期望极点的选取给出了一种合理解释, 而不是简单的一句依赖专家经验进行选择, 进而为经典的极点配置控制方法的有效运用提供更加有力的检验依据.
5) 从抗扰性的角度来看, 表 1~表 2中的阻尼比选择在0.5~0.9就可满足H∞性能要求.对照经典控制理论中无干扰时的控制器设计通常将阻尼比限制在0.4~0.8 [16-17], 显然具有一定的合理性.尽管此时的H∞性能在阻尼比大于0.8时就已实现, 但无干扰时的阻尼比选择在0.4~0.8已经具备足够的抗扰性了.由此可见经典控制理论中的二阶系统的设计已经考虑了外界干扰的潜在影响并具有足够的鲁棒性, 体现了系统内部结构特性以及外部专家经验的综衡问题.与传统的基于幅频裕度角度来解释传统二阶控制系统具有鲁棒性的特点不同, 本文是基于时域H∞性能的角度通过对比相应的阻尼比的选取范围来验证其鲁棒性的.再结合表 3可见, 对应H∞性能的闭环系统的最优阻尼比都是大于0.707的, 因此, 在干扰抑制控制系统设计中, 如果采用极点配置方法的话, 可将阻尼比相应提高一些.
6) 由表 3可知, 即使系统矩阵具有相同的特征根实部, 但是不同的α对应不同的系统矩阵表示形式, 进而具有不同的H∞性能以及不同的虚实部比值, 相应的阻尼比也不同.进而, 在追求控制系统卓越品质和友好交互的前提下, 精细设计好每一个环节, 如复特征根的实虚部值以及相应的变换矩阵, 最终实现系统矩阵的约简以及控制器的设计, 都是非常有意义的.
7) 从表 3可见, 复特征根的虚实部的比值最大为1.当比值趋近于1时, 所对应的阻尼比将是0.707.可见, 在H∞控制中, 对应H∞性能的闭环系统的最优阻尼比都是大于0.707的, 这是有干扰作用和无干扰作用下的控制效果的区别.
当特征根的虚部v固定, 则实部λ逐渐增大的仿真结果如表 4, 表 5所示.从表 4, 表 5可见, 随着实部的增加, H∞范数呈现减小趋势, 亦即说明, 随着系统的稳定裕度增大, 系统的抗扰能力逐渐增强, 以致当系统达到深度稳定时, 干扰将对系统没有什么影响.
表 4(Table 4)
表 4 H∞范数估计(α=2, v=2)Table 4 Estimation of H∞ norm(α=2, v=2)
| 表 4 H∞范数估计(α=2, v=2) Table 4 Estimation of H∞ norm(α=2, v=2) |
表 5(Table 5)
表 5 H∞范数估计(α=2, v=4)Table 5 Estimation of H∞ norm(α=2, v=4)
| 表 5 H∞范数估计(α=2, v=4) Table 5 Estimation of H∞ norm(α=2, v=4) |
4 结论针对时域二阶系统的H∞性能问题, 采用系统变换的方法, 探究了系统复特征根对于H∞性能的作用, 特别是研究了系统的虚部对于系统H∞性能的作用.由此得出, 在系统的H∞动态性能中, 系统的实部起着重要的作用, 系统的虚部不能超越实部的幅值, 进而保证系统的稳定特性和收敛特性, 这是系统对于外界干扰具有抑制作用的根本.虚部可以起到系统调节性能, 比如增加快速性等暂态性能, 是对系统品质的一种补充.与经典频域二阶系统理论进行了比较研究, 发现了阻尼比与H∞性能之间的某些关系, 并对于经典控制理论中阻尼比的参数选取范围进行了讨论, 从时域的H∞性能角度给出了一些启示, 进而对经典教材中的阻尼比在0.4~0.8之间的取值给出了解释.此外, 对于二阶系统的基于极点配置的控制器设计, 可以采用时域的H∞性能分析方法来确定期望闭环极点, 并形成时频交互综合分析实现控制器综合设计及性能验证.
参考文献
| [1] | Zames G. Feedback and optimal sensitivity:model reference transformations, multiplicative seminorms, and approximate inverses[J].IEEE Transactions on Automatic Control, 1981, 26(2): 301–320.DOI:10.1109/TAC.1981.1102603 |
| [2] | Francis B, Zames G. On H∞-optimal sensitivity theory for SISO feedback systems[J].IEEE Transactions on Automatic Control, 1984, 29(1): 9–16.DOI:10.1109/TAC.1984.1103357 |
| [3] | Doyle J C, Glover K. State-space solutions to standard H2 and H∞ control problems[J].IEEE Transactions on Automatic Control, 1989, 34(8): 831–847.DOI:10.1109/9.29425 |
| [4] | Gahinet P. Explicit controller formulas for LMI-based H∞ synthesis[J].Automatica, 1996, 32(7): 1007–1014.DOI:10.1016/0005-1098(96)00033-7 |
| [5] | 解学书, 钟宜生. H∞控制理论[M]. 北京: 清华大学出版社, 1994. ( Xie Xue-shu, Zhong Yi-sheng. H∞ control theory[M]. Beijing: Tsinghua University Press, 1994.) |
| [6] | 黄琳. 稳定性与鲁棒性理论基础[M]. 北京: 科学出版社, 2003. ( Huang Lin. Fundamentals of stability and robustness theory[M]. Beijing: Science Press, 2003.) |
| [7] | Yang R N, Wang Y Z. Finite-time stability analysis and H∞ control for a class of nonlinear time-delay Hamiltonian systems[J].Automatica, 2013, 49(2): 390–401.DOI:10.1016/j.automatica.2012.11.034 |
| [8] | Ma Y C, Lim C C, Li S N, et al. H-infinity control in coherent wind measurement lidar signal processing[J].Circuits Systems and Signal Processing, 2018, 37(7): 2824–2836.DOI:10.1007/s00034-017-0710-7 |
| [9] | 邢双云, 张庆灵, 朱宝彦. 随机奇异T-S模糊系统的限时间鲁棒H∞控制[J].东北大学学报(自然科学版), 2014, 35(11): 1521–1524. ( Xing Shuang-yun, Zhang Qing-ling, Zhu Bao-yan. Robust finite-time H∞ control for stochastic singular T-S fuzzy systems[J].Journal of Northeastern University(Natural Science), 2014, 35(11): 1521–1524.) |
| [10] | Chilali M, Gahinet P. H-infinity with pole placement constraints:an LMI approach[J].IEEE Transactions on Automatic Control, 1996, 41(3): 358–367.DOI:10.1109/9.486637 |
| [11] | Sarjas A, Svecko R, Chowdhury A. An H∞ modified robust disturbance observer design for mechanical positioning systems[J].International Journal of Control, Automation and Systems, 2015, 13(3): 575–586.DOI:10.1007/s12555-013-0531-9 |
| [12] | Iwasaki T, Skelton R E. All controllers for the general H∞ control problem:LMI existence conditions and state space formulas[J].Automatica, 1994, 30(8): 1307–1317.DOI:10.1016/0005-1098(94)90110-4 |
| [13] | Akos L. On the bounded real lemma[J].Systems and Control Letters, 2002, 45(4): 339–346. |
| [14] | 王锡范. 现代电力系统分析[M]. 北京: 科学出版社, 2003. ( Wang Xi-fan. Modern power system analysis[M]. Beijing: Science Press, 2003.) |
| [15] | 陈珩. 电力系统稳态分析[M]. 3版. 北京: 中国电力出版社, 2007. ( Chen Heng. Power system steady state snalysis[M]. 3rd ed. Beijing: China Electric Power Press, 2007.) |
| [16] | 陈启宗, 王纪文. 线性系统理论与设计[M]. 北京: 科学出版社, 1988. ( Chen Qi-zong, Wang Ji-wen. Linear system theory and design[M]. Beijing: Science Press, 1988.) |
| [17] | Kailath T. Linear systems[M]. Englewood Cliffs: Prentice-Hall, 1980. |
