3.2.4 解对称正定矩阵方程组的平方根法

这篇具有很好参考价值的文章主要介绍了3.2.4 解对称正定矩阵方程组的平方根法。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

平方根法解对称正定方程组,数值计算方法,数值计算方法,Powered by 金山文档
平方根法解对称正定方程组,数值计算方法,数值计算方法,Powered by 金山文档

在工程技术问题中,常常需要求解系数矩阵是对称正定矩阵的线性代数方程组。对于这类方程组,若利用矩阵三角分解法求解,就可得到一个有效法平方根法,其设计原理。

定理3 若A为对称正定矩阵,则存在唯一分解

A=~L~L^(T) (3.28)

其中~L是对角元为正的下三角形矩阵(对称正定矩阵的这种分解称为楚列斯基(Cholesky)分解)。

证明 由矩阵三角分解基本原理,存在唯一杜利特尔分解A=LU.若以Ak,Lk,Uk,依次表示矩阵A,L,U的k阶顺序主子阵,则

det A = det (Lk,Uk) = det Lk • det Uk,

u11u2……ukk (k =1,2,--. , n).

因A 对称正定,det A, >0(4=1,2,•,几),从而飛u >0(i=1,2,..,n).手是A 又可进一步分解为 A=LDU。,其中

平方根法解对称正定方程组,数值计算方法,数值计算方法,Powered by 金山文档

因为A 是对称矩阵,又有

A=A'=(LDUO'=UE(DL')显然,U。是单位下三角形矩阵,DL”是上三角形矩阵. 由杜利特尔分解的唯一性知U。二L,故A=LDL"I另一方面,由山i>0知D又可分解成D=D'D'2,其中VuID12Vuzz53如令之=1D',则之 是对角元为正的下三角形短阵,且由(3.29)式可得分解式(3.28)• 由从上述分解过程可以看出,这种分解是唯一的。下面导出实现分解A=记元°的递推算式.设美油

𤠒𡒗𢝵 𥮳食爺

采用自左向右逐列计算待定数l,的计算过程.由矩阵乘法规则与相等条件,依次

可得到确定立 的第一列元的算式

4n

=Nan ,ln =as/ lni

(¡=2.3..

,几)

青Y=送(3.30)

以及在算出亡的第—列至第;-1列元后确定第i列元的算式 !

號武手,依湖

16-(0-82)%

2

節路食(3.31)

帮式

因此,用楚列斯基分解解对称正定矩阵方程组 AX=6的过程可归纳为

1°实现楚列斯基分解,即

(a)按算式(3.30)计算立的第一列元;

(b)对j=2,3,•,几,按算式(3.31)计算之的第j列元.

甘本2°求解三角形方程组 云Y=6,相应的递推算式是

1:=(6-246o1014 Ci=2.3..0)。含

3°求解三角形方程组 亡-区=了,相应的递推算式是

(3.32)

不路其目丽-楼p

14=101-24.5)/2

(¡=n-1,...,2,1).

上述求解对称正定矩阵方程组的方法称为平方根法:

例5 用平方根法解方程组

改进的平方

根法

2

-27

12,

1-2

-3

14J 83.

解该方程组的系数矩阵是对称正定矩阵,可用平方根法求解.按算

式(3.30)计算亡的第一列元得

41=2,

21=1,

2,=-1.

按算式(3.31)依次计算亡的第二列与第三列元得

2=1, 12=-2, 153=3.

2

0

=2

"

解讠Y=6得y1=5,92=0,93=3.

解立"X=飞得x=1,52=2,81=2.

用平方根法解系数矩阵是几阶对称正定矩阵的线性方程组,当几较大时约

需作口。次乘除法运算(是高斯消元法或杜利特尔分解法的一半).此外,平方根

法还具有数值稳定、存储量小(利用对称性只需用一维数组存放矩阵4对角线

及对角线以下元,并将算得的 工 元存放在A 对应元的位置上)等优点•但平方根法在计算亡的对角线上元时需要用到开方运算文章来源地址https://www.toymoban.com/news/detail-732531.html

改进的平方根算法:

平方根法解对称正定方程组,数值计算方法,数值计算方法,Powered by 金山文档
平方根法解对称正定方程组,数值计算方法,数值计算方法,Powered by 金山文档
平方根法解对称正定方程组,数值计算方法,数值计算方法,Powered by 金山文档

到了这里,关于3.2.4 解对称正定矩阵方程组的平方根法的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 排列矩阵和三角矩阵——Matlab解线性方程组(2)

    目录 前言 一、排列矩阵是什么? 二、三角形矩阵 总结         上一篇文章讲了线性方程组的高斯消元法 。本文是一个辅助概念,讲解上文得到的P矩阵和L与U矩阵所代表的排列矩阵和上三角矩阵。         排列矩阵(permutation matrix)是单位矩阵经过行列交换而得到的新矩

    2024年02月07日
    浏览(50)
  • LA@2@1@线性方程组和简单矩阵方程有解判定定理

    线性方程组有解判定 线性方程组 A x = b Abold{x}=bold{b} A x = b 有解的 充分必要条件是它的系数矩阵A和增广矩阵 ( A , b ) (A,bold{b}) ( A , b ) 具有相同的秩 R ( A ) = R ( A , b ) R(A)=R(A,bold{b}) R ( A ) = R ( A , b ) ,记 r = R ( A ) = R ( A , b ) r=R(A)=R(A,bold{b}) r = R ( A ) = R ( A , b ) : 若 r = n r=n r = n 有

    2024年02月12日
    浏览(41)
  • Markdown:常用公式、行列式、矩阵、方程组等

        当前整理出来的皆为实际使用过的,欢迎大佬路过补充说明或者指正错误点。无用请轻喷。 1.1 常用公式符号 1.1.1 上下标 显示效果 公式代码 描述 x y x^y x y $x^y$ 或 $x^{y}$ 上标,若独显一个上标直接用 ^ ,若需要实现: x x + y x^{x+y} x x + y ,则用 {} 即可 x y x_y x y ​ $

    2024年02月05日
    浏览(42)
  • 线性方程组系数矩阵的秩与解的个数的关系

    齐次方程组: A x = 0 Ax=0 A x = 0 系数矩阵 A n × n A_{n×n} A n × n ​ 的秩 解的个数 满秩: r ( A ) = n r(A)=n r ( A ) = n 仅有零解 不满秩: r ( A ) = r n r(A)=rn r ( A ) = r n 有无穷多解 注: 齐次线性方程 A x = 0 Ax=0 A x = 0 一定有解. 当 r ( A ) = r n r(A)=rn r ( A ) = r n 时, 基础解系 (线性无关的

    2024年02月01日
    浏览(48)
  • 【算法竞赛模板】求解线性方程组是否有解(求解矩阵的秩)

        在实际运用中需判断线性方程组有无解,可以通过矩阵运算判断线性方程组是否有解 线性方程组有无解总结: 矩阵求解秩流程:    所以:当我们遇到题目问线性方程组是否有解时,只需求解系数矩阵的秩与增广矩阵的秩的关系 。我们可以通过分别求系数矩阵与增

    2024年02月12日
    浏览(39)
  • 【线代】矩阵的秩和线性方程组的解的情况

    行最简型矩阵 :(也可以叫做行最简阶梯型矩阵,或者行简化阶梯型矩阵),其特点是:非零行的首非零元为1,且这些非零元所在的列的其它元素都为0。所谓的行最简的意思就是对应的方程组是“最简单的”,就是说,对应的方程组,最多只需要移项就行了,不再需要其他任何

    2024年01月19日
    浏览(41)
  • 线性代数——线性方程组和矩阵(Linear and Matrices)

    1.Identify which of the following equations are linear equations: (判断哪些是线性方程) 只有(4)是,一般形式如下 特征:每一项都是一次的,也不代幂什么的 线性方程组(System of linear equations) ai,j是系数(i代表是第几个方程里,j是代表在方程里的第几个),b1是右端项,xj是未

    2023年04月08日
    浏览(43)
  • latex 常用数学符号( 二项式系数、矩阵、数组、方程与方程组、条件定义、括号、括号尺寸、字体)

    类型 符号 LaTeX 二项式系数 ( n k ) binom{n}{k} ( k n ​ ) binom{n}{k} 小型二项式系数 ( n k ) tbinom{n}{k} ( k n ​ ) tbinom{n}{k} 大型二项式系数 ( n k ) dbinom{n}{k} ( k n ​ ) dbinom{n}{k} x y z v begin{matrix}x y \\\\z vend{matrix} x z ​ y v ​ ∣ x y z v ∣ begin{vmatrix}x y \\\\z vend{vmatrix} ​ x z ​ y v ​ ​ ∥

    2024年02月06日
    浏览(38)
  • 数值分析第二次作业-求解系数矩阵为Hilbert 矩阵的线性方程组

    现要求解系数矩阵由16 阶Hilbert 方程组构成的线性方程组,右端项为  即要求解方程组Ax = b,其中 A=A0,b=b0  分别用高斯-赛德尔方法、最速下降法、共轭梯度法求解如下。 2.1. 高斯-赛德尔方法     2.2. 最速下降法    2.3. 共轭梯度法  在最速下降法中,搜索方向p取的是函数减

    2024年02月05日
    浏览(82)
  • 【考研数学】矩阵、向量与线性方程组解的关系梳理与讨论

    两个原因让我想写这篇文章,一是做矩阵题目的时候就发现这三货经常绑在一起,让人想去探寻其中奥秘;另一就是今天学了向量组的秩,让我想起来了之前遗留下来的一个问题:到底存不存在系数矩阵的秩和增广矩阵的秩之差比 1 大的情况?可能这个问题有点抽象,不过看

    2024年02月11日
    浏览(62)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包