计算方法(三)平方根法及其改进解线性方程组

一:概述

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

二:具体步骤

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

平方根法中的分解方法:
在这里插入图片描述

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

然后再大致解释一下

首先是平方根法:根据矩阵分析的理论,对于任意的正定H阵A, [公式] 存在唯一的正线下三角矩阵L,使得 [公式] ,把系数矩阵分解为三角矩阵后,我们就可以使用上次文章中的回代公式求解。

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

在这里插入图片描述
在这里插入图片描述
(2)使用回代公式得到X

已知[公式]使用两次回代公式,先求出y,再求出X即可。
在这里插入图片描述

然后是改进的平方根法:在平方根法中,对角元素的计算涉及开方,如何避免?我们可以对L做列变化,对 [公式] 做行变化,分别提出对角元素组成对角矩阵,再把两个对角阵相乘,即可避免开方运算。
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
直接从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
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分解的特例,但是不同的是,可以利用矩阵的对称性试编程起来更方便,虽然这么做会损害程序的泛用性。

  • 13
    点赞
  • 95
    收藏
    觉得还不错? 一键收藏
  • 1
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论 1
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值