方根法公式_计算方法（三）平方根法及其改进解线性方程组

和泉千晶，这位姐太……………… 遭不住遭不住。

一：概述

本篇文章介绍解线性方程组的平方根法及改进平方根法，适用范围为系数矩阵为正定Hermite矩阵(下称H阵)的线性方程组。这个方法的理论依据我觉得是来自Schur引理的H阵结构定理，从这个角度我们就可以理解课程上教授的那些摸不着头脑的机械的计算步骤背后的原理。

二：具体步骤

按照惯例，我们先举两个例子介绍一下两种方法是如何进行计算的，方便那些只想知道怎么算的同学。其实这个玩意和LU分解的计算方法完全一样

平方根法中的分解方法：

改进平方根中的分解方法：这张图片第一次编辑好之后鬼使神差地点了差号，气到吐血，第二遍就没了那么多耐心。

然后再大致解释一下

首先是平方根法：根据矩阵分析的理论，对于任意的正定H阵A，

$\exists$ 存在唯一的正线下三角矩阵L，使得

$A=LL^{T}$ ，把系数矩阵分解为三角矩阵后，我们就可以使用上次文章中的回代公式求解。

(1)分解系数矩阵得到正线下三角矩阵L

(2)使用回代公式得到X

已知

$LL^{T}X=b，设y=L^{T}X,则Ly=b。$ 使用两次回代公式，先求出y，再求出X即可。

然后是改进的平方根法：在平方根法中，对角元素的计算涉及开方，如何避免？我们可以对L做列变化，对

$L^{T}$ 做行变化，分别提出对角元素组成对角矩阵，再把两个对角阵相乘，即可避免开方运算。

直接从A矩阵得到以上两个三角矩阵的公式如下：

然后还是再使用两次回代公式，就可以得到方程组的解。

三：原理分析

为什么这么算呢？仔细一想，其实和上一篇的原理完全一样，当系数矩阵是正定H阵时，U矩阵就自然与L矩阵共轭转置。

四：算法实现(MATLAB)

改进的平方根法就不用再列出了，因为和上一篇的完全一样。

function[x]=pingfanggenfa(A,b)%平方根法

n=length(A);

for k=1:n

A(k,k)=sqrt(A(k,k));

A(k+1:n,k)=A(k+1:n,k)/A(k,k);

for j=k+1:n

A(j:n,j)=A(j:n,j)-A(j:n,j)*A(j,k);

end

for j=1:n-1

b(j)=b(j)/A(j,j);

b(j+1:n)=b(j+1:n)-b(j)*A(j+1:n,j);

end

for j=n:-1:2

b(j)=b(j)/A(j,j);

b(1:j-1)=b(1:j-1)-b(j)*A(1:j-1,j);

end

x=b

五：总结

这一篇说的两种方法，都要求系数矩阵是对称正定阵。这两种方法都可以看成是上一篇lu分解的特例，但是不同的是，可以利用矩阵的对称性试编程起来更方便，虽然这么做会损害程序的泛用性。

上一篇水法天天：计算方法(二)直接三角分解法解线性方程组zhuanlan.zhihu.com

方根法公式_计算方法（三）平方根法及其改进解线性方程组相关推荐

计算方法（三）平方根法及其改进解线性方程组
一:概述本篇文章介绍解线性方程组的平方根法及改进平方根法,适用范围为系数矩阵为正定Hermite矩阵(下称H阵)的线性方程组.这个方法的理论依据我觉得是来自Schur引理的H阵结构定理,从这个角度我 ...
方根法公式_方根的简易算法
1.1已知平方根整数部分,求平方根小数部分的算法: 任何一个数都可以写成某个整数的平方加上一个余数,如15可以写成32+6.8.9可以写成22+4.9.任何一个正数写成某个整数的平方加一个余数项的通式 ...
方根法公式_仓储管理笔记之库存分析法：ABC分析法、区域合并法......
导读国内有庞大的仓储物流从业人员队伍(根据中国物流与采购联合会的调查,2016年底我国物流从业人员5012万,是人员增速最快的行业),很多人只是想深入了解仓库从无序到有序,从源头开始应该如何管理.如 ...
方根法公式_初中数学根式运算法则公式
初中数学根式运算法则公式2020-10-27 15:17:20文/樊越很多同学都学习了根式,小编整理了一些根式运算法则,大家一起来看看吧. 根式运算法根式开方法则是根式的运算法则之一,算术根开n次 ...
方根法公式_层次分析法中用方根法计算权重在Excel中的具体操作
Exce l 表 A B C D E F G H I J K L M N 1 总目标子目标 1 子目标 2 子目标 3 子目标 4 M ij 几何平均数权重 W i AW i ...
方根法公式_(完整word版)层次分析法中用方根法计算权重在Excel中的具体操作
Exce l 表 A B C D E F G H I J K L M N 1 总目标子目标 1 子目标 2 子目标 3 子目标 4 M ij 几何平均数权重 W i AW i ...
米的换算单位和公式_小学三年级数学常用公式和单位换算，孩子复习宝典！
1 单位换算小技巧 1.将同类单位从大到小或者从小到大依次排列如长度单位:千米.米.分米.厘米 2.记住各单位之间的进率,高级单位X进率=低级单位:低级单位÷进率=高低单位,差多少位就进行多少次换算 ...
sin的傅里叶变换公式_第三篇傅里叶变换重要公式.pdf
第三章傅里叶变换重要概念与重要公式一.傅里叶级数 1.三角函数形式的傅里叶级数任何周期信号f (t ) 可以分解为 ∞ (1)f (t ) a + a cos n ωt +b sin n ωt ...
卡尔曼滤波五个公式_基于ROS的卡尔曼滤波姿态解算
前段时间由于项目关系需要实现基于卡尔曼滤波的姿态解算,也就是说融合加速度计.陀螺仪及磁罗盘进行AHRS姿态解算,得到机器人的姿态角. 本文的学习需要有一定的卡尔曼滤波器基础,可以参考白巧克力亦唯心的卡 ...

方根法公式_计算方法（三）平方根法及其改进解线性方程组

方根法公式_计算方法（三）平方根法及其改进解线性方程组相关推荐

最新文章

热门文章