主题
Search

稳定多项式


一个 实多项式 P 被称为稳定的,如果它的所有 都位于 左半平面 中。“稳定”一词用于描述这种多项式,因为在线性伺服机构理论中,系统表现出 e^(st) 形式的无外力时变运动,其中 s 是某个 实多项式 P(s)=0 的根。因此,一个系统在机械上是稳定的 当且仅当 P 是一个稳定多项式。

多项式 x+a 是稳定的 当且仅当 a>0,且 不可约多项式 x^2+ax+b 是稳定的 当且仅当 ab 都大于零。劳斯-赫尔维茨定理 可用于确定多项式是否稳定。

给定两个实多项式 PQ,如果 PQ 是稳定的,那么它们的乘积 PQ 也是稳定的,反之亦然(Séroul 2000,第 280 页)。因此,稳定实多项式的系数要么都为正,要么都为负(尽管这不是一个 充分 条件,如反例 x^3+x^2+x+1 所示)。此外,稳定多项式在 x>=0 时的值永远不为零,并且与多项式的系数具有相同的符号。

无需预先知道多项式的根,就可以使用 Strelitz (1977) 的以下定理来判断多项式是否稳定。设 A=x^n+a_(n-1)x^(n-1)+...+a_0 为一个实多项式,其根为 alpha_1,...,alpha_n,构建 B=x^m+b_(m-1)x^(m-1)+...+b_0 为首一实多项式,其次数为 m=n(n-1)/2,根为 alpha_i+alpha_j,其中 1<=i<j<=n。那么 A 是稳定的 当且仅当 AB 的所有系数都为正数(Séroul 2000,第 281 页)。

例如,给定三阶多项式 A=x^3+ax^2+bx+c,根和多项式 B 由下式给出

 B=x^3+2ax^2+(a^2+b)x+(ab-c).
(1)

通过求解要求 AB 的每个系数都大于零的不等式,可以得到 A 稳定的条件为 a>0b>00<c<ab

同样,对于四阶多项式 A=x^4+ax^3+bx^2+cx+d,根和多项式为

 x^6+3ax^5+(3a^2+2b)x^4+(a^3+4ab)x^3+(2a^2b+b^2+ac-4d)x^2+(ab^2+a^2c-4ad)x+(abc-c^2-a^2d),
(2)

因此,A 稳定的条件可以解析为 a>0b>00<c<ab0<d<(abc-c^2)/a^2

五阶多项式为

 x^(10)+4ax^9+(6a^2+3b)x^8+(4a^3+9ab+c)x^7+(a^4+9a^2b+3b^2+4ac-3d)x^6+(3a^3b+6ab^2+5a^2c+2bc-5ad-11e)x^5+(3a^2b^2+b^3+2a^3c+6abc-c^2-2a^2d-2bd-22ae)x^4+(ab^3+4a^2bc+b^2c-4cd-16a^2e-4be)x^3+(2ab^2c+a^2c^2-bc^2+a^2bd+b^2d-3acd-4d^2-4a^3e-9abe+7ce)x^2+(abc^2-c^3+ab^2d-4ad^2-4a^2be-b^2e+4ace+4de)x+(abcd-c^2d-a^2d^2-ab^2e+bce+2ade-e^2).
(3)

以下 Wolfram 语言 代码计算根和多项式 B 以及从系数获得的不等式

  RootSumPolynomial[r_List, x_]:=Module[
      {n = Length[r], i, j},
      RootReduce@Collect[Expand[
        Times@@((x - #)&/@Flatten[
          Table[r[[i]] + r[[j]], {i, n},
            {j, i+1, n}]])
      ], x]
  ]
  RootSumPolynomial[p_?PolynomialQ, x_]:=
    RootSumPolynomial[RootList[p, x], x]
  RootList[p_?PolynomialQ, x_]:=
    x /. {ToRules[Roots[p==0, x,
      Cubics -> False, Quartics -> False
    ]]}
  RootSumInequalities[p_?PolynomialQ, x_]:=
    And @@ (# > 0& /@
      Flatten[CoefficientList[#, x]& /@
        {RootSumPolynomial[p, x], p}])

而以下代码将三次情况下的不等式简化为最小集合

  Resolve[Exists[x, Element[(a | b | c | x), Reals],
    RootSumInequalities[x^3 + a x^2 + b x + c, x]
  ], {a, b, c}]

另请参阅

左半平面, 劳斯-赫尔维茨定理

使用 Wolfram|Alpha 探索

参考文献

Séroul, R. “稳定多项式。” Programming for Mathematicians. 第 10.13 节。柏林:Springer-Verlag,第 280-286 页,2000 年。Strelitz, S. “关于劳斯-赫尔维茨问题。” Amer. Math. Monthly 84, 542-544, 1977 年。Tóth, J.; Szili, L.; 和 Zachár, A. “多项式的稳定性。” Mathematica Educ. Res. 7, 5-12, 1998 年。

在 Wolfram|Alpha 上被引用

稳定多项式

请引用为

Weisstein, Eric W. “稳定多项式。” 来自 MathWorld--Wolfram Web 资源。 https://mathworld.net.cn/StablePolynomial.html

学科分类