给出矩阵 m 的特征多项式.
CharacteristicPolynomial[{m,a},x]
给出关于 a 的广义特征多项式.
CharacteristicPolynomial
给出矩阵 m 的特征多项式.
CharacteristicPolynomial[{m,a},x]
给出关于 a 的广义特征多项式.
更多信息
- m 必须是一个方阵.
- 它可以包含数值或符号项.
- CharacteristicPolynomial[m,x] 实际上等价于 Det[m-id x],其中 id 是适当尺寸的单位矩阵. »
- CharacteristicPolynomial[{m,a},x] 实际上是 Det[m-a x]. »
范例
打开所有单元 关闭所有单元基本范例 (3)
CharacteristicPolynomial[{{1, 2}, {3, 4}}, x]Plot[%, {x, -5, 10}]m = (| | |
| - | - |
| a | b |
| c | d |);CharacteristicPolynomial[m, x]Det[m - x IdentityMatrix[2]]CharacteristicPolynomial[IdentityMatrix[3], λ]CharacteristicPolynomial[ConstantArray[0, {10, 10}], λ]范围 (17)
基本用法 (7)
CharacteristicPolynomial[{{1.1, 2.2, 3.25}, {0.76, 4.6, 5}, {0.1, 0.1, 6.1}}, x]CharacteristicPolynomial[Table[N[1 / (i + j + 1), 20], {i, 3}, {j, 3}], x]CharacteristicPolynomial[{{1.2 + I, 3 - 2 I, 3 π}, {-0.2, 5I, 2}, {1, 2.3, E}}, x]CharacteristicPolynomial[{{(1/3), (1/2), (3/5)}, {(1/2), (4/5), 1}, {(3/5), 1, (9/7)}}, x]Plot[%, {x, -1.5, 2.5}]m = RandomReal[{1, 9}, {100, 100}];CharacteristicPolynomial[m, x];//Timingℱ = FiniteField[43, 2];
CharacteristicPolynomial[{{ℱ[12], ℱ[23], ℱ[34]}, {ℱ[45], ℱ[56], ℱ[67]}, {ℱ[78], ℱ[89], ℱ[90]}}, x]含有 CenteredInterval 对象的矩阵的特征多项式:
(m = Map[CenteredInterval, RandomReal[{-10, 10}, {3, 3}, WorkingPrecision -> 10], {2}])//MatrixFormp = CharacteristicPolynomial[m, x]ranrep[e_CenteredInterval] := e["Center"] + RandomInteger[{-1000, 1000}] / 1000 e["Radius"]
(mrep = Map[ranrep, m, {2}])//MatrixFormMapThread[IntervalMemberQ, {CoefficientList[p, x], CoefficientList[CharacteristicPolynomial[mrep, x], x]}]广义特征值 (4)
m1 = {{1, 2}, {5, 4}};m2 = {{4, 3}, {6, 4}};CharacteristicPolynomial[{m1, m2}, x]Det[m1 - x m2]a = {{1., 1.5, 2.}, {3.1, 2., 2.9}, {3., 2., 1.}};b = {{1.3, .5, 1.1}, {0., 1.5, 2.3}, {1., 0., 1.}};CharacteristicPolynomial[{a, b}, x]a = {{1, 1, 1}, {1, 0, 1}, {0, 0, 1}};b = {{0, 1, 1}, {0, 1, 1}, {1, 0, 0}};CharacteristicPolynomial[{a, b}, x]Eigenvalues[{a, b}]CharacteristicPolynomial[N[{a, b}, 20], y]a = {{x, 1 + x}, {1 - x, x}};b = {{1, 1}, {1, 2x}};
CharacteristicPolynomial[{a, b}, y]特殊矩阵 (6)
SparseArray[{{1, 3} -> 2, {2, 2} -> 3, {3, 1} -> 1, {4, 2} -> 5}, {4, 4}]CharacteristicPolynomial[%, x]SparseArray[{{x_, y_} /; Abs[x - y] < 3 -> 1}, {10, 10}]CharacteristicPolynomial[%, x]SymmetrizedArray[{{1, 1} -> 2, {1, 2} -> 1}, {2, 2}, Symmetric[All]]CharacteristicPolynomial[%, x]QuantityArray[{{1, 2}, {3, 4}}, {"Meters", "Meters"}]CharacteristicPolynomial[%, x]CharacteristicPolynomial[IdentityMatrix[12], λ]Factor[%]HilbertMatrix[n] 的特征多项式:
CharacteristicPolynomial[HilbertMatrix[5], λ]JordanMatrix[λ,n] 的特征多项式为
:
CharacteristicPolynomial[JordanMatrix[λ, 4], x] == (λ - x)^4//SimplifyCharacteristicPolynomial[JordanMatrix[λ, 5], x] == (λ - x)^5//SimplifyCompanionMatrix[{c0,c1,…,cn}] 的最小多项式:
CompanionMatrix[{c0, c1, c2, c3}]//MatrixFormCharacteristicPolynomial[%, x]应用 (6)
m = (| | | |
| -- | --- | --- |
| -6 | 28 | 21 |
| 4 | -15 | -12 |
| -8 | a | 25 |)cp = CharacteristicPolynomial[m, x]roots = SolveValues[cp == 0, x]roots /. a -> 32roots /. a -> 31roots /. a -> 33{Plot[Table[cp, {a, {32, 31, 33}}]//Evaluate, {x, -3, 6}, PlotLegends -> LineLegend[{a == 32, a == 31, a == 33}], ImageSize -> Small],
Plot[cp /. a -> 32, {x, 0.5, 2.5}, PlotLegends -> LineLegend[{a == 32}], ImageSize -> Small]}//Rowm = {{6, 5, 8, 1, 4}, {8, 0, 2, 7, 1}, {7, 4, 0, 8, 3}, {2, 3, 1, 1, 6}, {0, 4, 8, 2, 3}};cp = CharacteristicPolynomial[m, x]cp /. x -> 0Times @@ SolveValues[cp == 0, x] // FullSimplify与用 Det 直接计算的结果相比较:
Det[m]m = {{4, 3, 8, 9}, {8, 1, 0, 7}, {1, 2, 6, 7}, {2, 5, 1, 2}};cp = CharacteristicPolynomial[m, x](-1)^Length[m] - 1Coefficient[cp, x^Length[m] - 1]Total[SolveValues[cp == 0, x]]//FullSimplify与用 Tr 直接计算的结果相比较:
Tr[m]m = {{1, 2, 3}, {4, 5, 6}, {7, 8, 9}};SolveValues[CharacteristicPolynomial[m, x] == 0, x]与用 Eigenvalues 直接计算的结果相比较:
Eigenvalues[m]a = {{3, 9, 4}, {3, 7, 1}, {8, 1, 8}};cp = CharacteristicPolynomial[a, x]CharacteristicPolynomial[Transpose[a], x] == cpλ = SolveValues[cp == 0, x]Subscript[v, λ] = Join[NullSpace[a - λ[[1]]IdentityMatrix[3]], NullSpace[a - λ[[2]]IdentityMatrix[3]], NullSpace[a - λ[[3]]IdentityMatrix[3]]]//FullSimplifyEigensystem 给出相同的结果,尽管它按绝对值对特征值进行排序:
Eigensystem[a]Subsuperscript[v, λ, * ] = Join[NullSpace[a - λ[[1]]IdentityMatrix[3]], NullSpace[a - λ[[2]]IdentityMatrix[3]], NullSpace[a - λ[[3]]IdentityMatrix[3]]]//FullSimplify{MatrixPlot[Subscript[v, λ]], MatrixPlot[Subsuperscript[v, λ, * ]]}a = {{1., 1.5, 2.}, {3.1, 2., 2.9}, {3., 2., 1.}};b = {{1.3, .5, 1.1}, {0., 1.5, 2.3}, {1., 0., 1.}};CharacteristicPolynomial[{a, b}, x]λ = SolveValues[% == 0, x]Subscript[v, λ] = Join[NullSpace[a - λ[[1]]b], NullSpace[a - λ[[2]]b], NullSpace[a - λ[[3]]b]]//FullSimplify与用 Eigensystem 直接计算的结果相比较:
Eigensystem[{a, b}]属性和关系 (10)
特征多项式等价于 Det[m - id x]:
m = RandomInteger[9, {10, 10}];cp = CharacteristicPolynomial[m, x]% == Det[m - IdentityMatrix[Length[m]] x]广义的特征多项式等价于 Det[m-a x]:
{m, a} = RandomInteger[9, {2, 5, 5}];CharacteristicPolynomial[{m, a}, x]% == Det[m - a x]矩阵是其特征多项式的一个根 (Cayley–Hamilton 定理 [更多信息...]):
m = RandomInteger[10, {10, 10}];
cp = CharacteristicPolynomial[m, x];
MatrixPolynomialValue[cp, m, x] == ConstantArray[0, {10, 10}]MatrixPolynomialValue[cp, m, x]//AbsoluteTimingm = RandomReal[1, {10, 10}];Product[(v - x), {v, Eigenvalues[m]}] - CharacteristicPolynomial[m, x]//Simplify//Chop特征多项式的根的和是矩阵的迹 (Tr):
m = RandomReal[1, {10, 10}];roots = SolveValues[CharacteristicPolynomial[m, x] == 0, x];Total[roots] == Tr[m]同样,根的积是行列式 (Det):
Times @@ roots == Det[m]m = RandomReal[1, {10, 10}];CharacteristicPolynomial[m, x] == CharacteristicPolynomial[Transpose[m], x]t1 = (| | | |
| - | - | - |
| 1 | 0 | 0 |
| 0 | 1 | 0 |
| 0 | 0 | 1 |);t2 = (| | | |
| - | - | - |
| 1 | 4 | 0 |
| 0 | 1 | 3 |
| 0 | 0 | 1 |);t3 = (| | | |
| - | - | - |
| 1 | 0 | 0 |
| 5 | 1 | 0 |
| 6 | 7 | 1 |);CharacteristicPolynomial[t1, x] == CharacteristicPolynomial[t2, x] == CharacteristicPolynomial[t3, x]cl = RandomInteger[9, 10];
p = cl.x ^ Range[0, 9] + x ^ 10CharacteristicPolynomial[CompanionMatrix[cl], x] == pMatrixMinimalPolynomial[m,λ] 整除 CharacteristicPolynomial[m,λ],其商为一个(可能为常数的)多项式:
mat = (| | | |
| - | -- | -- |
| 0 | -1 | 1 |
| 1 | 2 | -1 |
| 1 | 1 | 0 |);
PolynomialQuotientRemainder[CharacteristicPolynomial[mat, t], MatrixMinimalPolynomial[mat, t], t]a = (| | | | |
| -- | -- | - | - |
| -6 | 4 | 0 | 9 |
| -3 | 0 | 1 | 6 |
| -1 | -2 | 1 | 0 |
| -4 | 4 | 0 | 7 |);
χ = CharacteristicPolynomial[a, λ];
μ = MatrixMinimalPolynomial[a, λ];
Roots[μ == 0, λ] == DeleteDuplicates[Roots[χ == 0, λ]]PolynomialRemainder[μ ^ LCM@@(Max@Cases[#, λ ^ k_ :> k, -1]& /@ {μ, χ}), χ, λ]历史
2003年引入 (5.0) | 在以下年份被更新:2007 (6.0) ▪ 2024 (14.0)
文本
Wolfram Research (2003),CharacteristicPolynomial,Wolfram 语言函数,https://reference.wolfram.com/language/ref/CharacteristicPolynomial.html (更新于 2024 年).
CMS
Wolfram 语言. 2003. "CharacteristicPolynomial." Wolfram 语言与系统参考资料中心. Wolfram Research. 最新版本 2024. https://reference.wolfram.com/language/ref/CharacteristicPolynomial.html.
APA
Wolfram 语言. (2003). CharacteristicPolynomial. Wolfram 语言与系统参考资料中心. 追溯自 https://reference.wolfram.com/language/ref/CharacteristicPolynomial.html 年
BibTeX
@misc{reference.wolfram_2026_characteristicpolynomial, author="Wolfram Research", title="{CharacteristicPolynomial}", year="2024", howpublished="\url{https://reference.wolfram.com/language/ref/CharacteristicPolynomial.html}", note=[Accessed: 13-September-2026]}
BibLaTeX
@online{reference.wolfram_2026_characteristicpolynomial, organization={Wolfram Research}, title={CharacteristicPolynomial}, year={2024}, url={https://reference.wolfram.com/language/ref/CharacteristicPolynomial.html}, note=[Accessed: 13-September-2026]}