矩阵——条件数与方程组的性态

Posted 炫云云

tags:

篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了矩阵——条件数与方程组的性态相关的知识,希望对你有一定的参考价值。

“病态”矩阵与方程

考虑线性方程组
( 2 6 2 6.00001 ) ( x 1 x 2 ) = ( 8 8.00001 ) \\left(\\begin{array}{cc} 2 & 6 \\\\ 2 & 6.00001 \\end{array}\\right)\\left(\\begin{array}{l} x_{1} \\\\ x_{2} \\end{array}\\right)=\\left(\\begin{array}{c} 8 \\\\ 8.00001 \\end{array}\\right) (2266.00001)(x1x2)=(88.00001)
它有准确解 x = ( 1 , 1 ) T x=(1,1)^{T} x=(1,1)T, 如果方程组右端项发生微小的变化 δ b = ( 0 , 0.00001 ) T \\delta \\boldsymbol{b}=(0,0.00001)^{T} δb=(0,0.00001)T
( 2 6 2 6.00001 ) ( x ~ 1 x ~ 2 ) = ( 8 8.00002 ) \\left(\\begin{array}{cc} 2 & 6 \\\\ 2 & 6.00001 \\end{array}\\right)\\left(\\begin{array}{l} \\tilde{x}_{1} \\\\ \\tilde{x}_{2} \\end{array}\\right)=\\left(\\begin{array}{c} 8 \\\\ 8.00002 \\end{array}\\right) (2266.00001)(x~1x~2)=(88.00002)
其解为 x ~ = ( − 2 , 2 ) T \\tilde{\\boldsymbol{x}}=(-2,2)^{T} x~=(2,2)T, 可以看出,
∥ x − x ~ ∥ ∞ ∥ x ∥ ∞ = ∥ ( 3 − 1 ) ∥ ∞ ∥ ( 1 1 ) ∥ ∞ = 3 , ∥ δ b ∥ ∞ ∥ b ∥ ∞ = ∥ ( 0 0.00001 ) ∥ ∞ ∥ ( 8 8.00001 ) ∥ ∞ = 0.00001 8.00001 ≈ 1 800000 , \\frac{\\|\\boldsymbol{x}-\\tilde{\\boldsymbol{x}}\\|_{\\infty}}{\\|\\boldsymbol{x}\\|_{\\infty}}=\\frac{\\left\\|\\left(\\begin{array}{c} 3 \\\\ -1 \\end{array}\\right)\\right\\|_{\\infty}}{\\left\\|\\left(\\begin{array}{c} 1 \\\\ 1 \\end{array}\\right)\\right\\|_{\\infty}}=3, \\\\ \\frac{\\|\\delta \\boldsymbol{b}\\|_{\\infty} }{\\|\\boldsymbol{b}\\|_{\\infty}}= \\frac{\\left\\|\\left(\\begin{array}{c} 0 \\\\ 0.00001 \\end{array}\\right)\\right\\|_{\\infty}}{{\\left\\|\\left(\\begin{array}{c} 8 \\\\ 8.00001 \\end{array}\\right)\\right\\|_{\\infty}}} =\\frac{0.00001}{8.00001} \\approx \\frac{1}{800000}, xxx~=(11)(31)=3,bδb=(88.00001)(00.00001)=8.000010.000018000001,
即解的相对误差是右端的相对误差的 2400000 倍.

算子范数: 矩阵——向量范数和矩阵范数

定 义 1 \\large\\color{magenta}{\\boxed{\\color{brown}{定义 1} }} 1 如果线性方程组 A x = b \\boldsymbol{A x}=\\boldsymbol{b} Ax=b 中, A \\boldsymbol{A} A b \\boldsymbol{b} b 的元素的微小变化, 就会引起方程组解的 巨大变化,则称方程组为“病态”方程组,矩阵 A \\boldsymbol{A} A 称为“病态”矩阵. 否则称方程组为“良 态”方程组,矩阵 A A A 称为“良态”矩阵.

我们需要一种能刻画矩阵和方程组“病态”标准的量. 设线性方程组 A x = b \\boldsymbol{A} \\boldsymbol{x}=\\boldsymbol{b} Ax=b 中, A \\boldsymbol{A} A 为 非奇异矩阵, x \\boldsymbol{x} x 为方程组的准确解. 考虑 b \\boldsymbol{b} b 有误差 δ b \\delta \\boldsymbol{b} δb, 其解为 x + δ x \\boldsymbol{x}+\\boldsymbol{\\delta x} x+δx ,即
A ( x + δ x ) = b + δ b (1) \\boldsymbol{A}(\\boldsymbol{x}+\\boldsymbol{\\delta} \\boldsymbol{x})=\\boldsymbol{b}+\\boldsymbol{\\delta} \\boldsymbol{b}\\tag{1} A(x+δx)=b+δb(1)
因为 x \\boldsymbol{x} x 为方程组的准确解, 故 (1)为
A δ x = δ b ,  即  δ x = A − 1 δ b , \\boldsymbol{A} \\delta \\boldsymbol{x}=\\boldsymbol{\\delta} \\boldsymbol{b}, \\text { 即 } \\boldsymbol{\\delta} \\boldsymbol{x}=\\boldsymbol{A}^{-1} \\boldsymbol{\\delta} \\boldsymbol{b}, Aδx=δb,以上是关于矩阵——条件数与方程组的性态的主要内容,如果未能解决你的问题,请参考以下文章

条件数:逆矩阵与线性方程组

N元线性方程组AX=B无解的充要条件是?

线性代数——矩阵

matlab中如何求解齐次线性方程组(代数矩阵)的非零解

线性代数,求解矩阵方程

非齐次线性方程组Ax=b的解都有哪些?