数值分析·学习 | 平方根法和追赶法matlab实现

这篇具有很好参考价值的文章主要介绍了数值分析·学习 | 平方根法和追赶法matlab实现。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

目录

一、前言:

二、算法描述:

三、实现代码:

1、平方根法:

2、改进的平方根法:

3、追赶法:

四、总结:


一、前言:

个人学习内容分享


二、算法描述:

平方根法:

        如果A为n阶对称正定矩阵,则存在一个实的非奇异下三角矩阵L,使,当限定L的对角元素为正时,这种分解是唯一的。

        下面我们用直接分解方法来确定计算L元素的递推公式,因为

                

        其中.由矩阵乘法及(当j<k时),得

                数值分析·学习 | 平方根法和追赶法matlab实现

        于是得到对于求解Ax=b的平方根法计算公式:

        对于j=1,2,...,n

  

数值分析·学习 | 平方根法和追赶法matlab实现

        求解Ax=b,即求解两个三角形方程组:

         ①Ly=b,求y;     ②LTx=y,求x.

数值分析·学习 | 平方根法和追赶法matlab实现

        其中:

                        

所以

改进的平方根法:

        为了避免开方,我们下面用定理9的分解式A=

即 

        由矩阵乘法,并注意,得

数值分析·学习 | 平方根法和追赶法matlab实现

        于是得到计算L得元及D得对角元素公式:

        对于i=1,2,...,n.

        为了避免重复计算,我们引进tij=lijdj,由上式得到按行计算L,T元素得公式:d1=a11

对于i=2,3,...,n.

        

        

        

        计算出T=LD的第i行元素tij(j=1,2,...,i-1)后,存放在A的第i行相应位置,然后再计算L的第i行元素,存放在A的第i行.D的对角线元素存放在A的相应位置。例如

        ​​​​​​​        

        对称正定矩阵A按分解和按分解计算量差不多,但分解不需要开方计算。 

        求解Ly=b,DLTx=y计算公式:

数值分析·学习 | 平方根法和追赶法matlab实现

追赶法:

对于A为三对角矩阵时,可分解成以下形式

即可得到

        数值分析·学习 | 平方根法和追赶法matlab实现

从而求解Ax=f等价于解两个三角形方程组:

        ①Ly=f,求y,②Ux=y,求x.

进而追赶法公式为:

(1)计算{Bi}的递推公式

        

        

(2)解Ly=f

        

        

(3)解Ux=y

        

        数值分析·学习 | 平方根法和追赶法matlab实现文章来源地址https://www.toymoban.com/news/detail-408516.html

三、实现代码:

1、平方根法:

function X=squareRoot(A,b)
%%功能:平方根算法
n=length(A);
y=b;
X=y;
%变换矩阵
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,k)*A(j,k);
    end
end
%求解矩阵
for k=1:n
    y(k)=(b(k)- sum(A(k,1:k-1)*y(1:k-1)))/A(k,k);
end
A=A';
for k=n:-1:1
    X(k)=(y(k)-sum(A(k,k+1:n)*X(k+1:n)))/A(k,k);
end

2、改进的平方根法:

function X=ImprovedSquareRoot(A,b)
%%功能:改进的平方根公式
n=length(A);
v=b;
y=b;
X=b;
%变换矩阵
for j=1:n
    for k=1:j-1
        v(k)=A(j,k)*A(k,k);
    end
    A(j,j)=A(j,j)-sum(A(j,1:j-1)*v(1:j-1));
    for i=j+1:n
        A(i,j)=(A(i,j)-sum(A(i,k)*v(1:j-1)))/A(j,j);
    end
end
%求解矩阵
D=diag(diag(A));
A=A-D+diag(diag(ones(3)));
for k=1:n
    y(k)=(b(k)-sum(A(k,1:k-1)*y(1:k-1)))/A(k,k);
end
d=diag(D);
X(n)=y(n)/d(n);
for k=n-1:-1:1
    X(k)=y(k)/d(k)-sum(A(k+1:n,k).*X(k+1:n));
end
end

3、追赶法:

function X=CatchUpMethod(a,b,c,f)
%%功能:追赶法
n=length(b);
B=f;
B(1)=c(1)/b(1);
for k=2:n-1
    B(k)=c(k)/(b(k)-a(k)*B(k-1));
end
y(1)=f(1)/b(1);
for k=2:n
    y(k)=(f(k)-a(k-1)*y(k-1))/(b(k)-a(k-1)*B(k-1));
end
X(n)=y(n);
for k=n-1:-1:1
    X(k)=y(k)-B(k)*X(k+1);
end
end 

四、总结:

到了这里,关于数值分析·学习 | 平方根法和追赶法matlab实现的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

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

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

相关文章

  • Python算法例4 求平方根

    实现int sqrt(int x)函数,计算并返回x的平方根。 sqrt(3)=1;sqrt(4)=2;sqrt(5)=2;sqrt(17)=4。 要实现计算整数x的平方根函数sqrt(x),可以使用二分查找法。 首先,我们定义一个变量left = 0用来表示搜索区间的左边界,以及一个变量right = x用来表示搜索区间的右边界。初

    2024年02月05日
    浏览(39)
  • leetcode69 x 的平方根

    题目变形为找到 f ( x ) = x 2 − c = 0 f(x)=x^2-c=0 f ( x ) = x 2 − c = 0 的根,其中 x x x 是非负整数。由于 f ( 0 ) = − c ≤ 0 , f ( c ) = c 2 − c ≥ 0 f(0)=-cle0,f(c)=c^2-cge0 f ( 0 ) = − c ≤ 0 , f ( c ) = c 2 − c ≥ 0 ,则 [ 0 , c ] [0,c] [ 0 , c ] 之间必然存在一个根,使用二分法。 但是由于计算

    2024年02月22日
    浏览(42)
  • leetcode69---x 的平方根

    大家好,我是大唐,刚刷完了几道经典的leetcode题,今天给大家分享一道leetcode上面的二分查找经典题型---x 的平方根,我们往下看。 给你一个非负整数  x  ,计算并返回  x  的  算术平方根  。 由于返回类型是整数,结果只保留  整数部分  ,小数部分将被  舍去 。 注意

    2024年03月18日
    浏览(43)
  • C语言—求平方根(sqrt函数)

            在数学当中,我们知道了平方根。那么在C语言当中求一个数的平方根是如何实现的呢?今天我们就来讲解。  sqrt()函数为库函数,所以要包含对应的头文件,这个头文件包含了sqrt()函数的定义  下图中,x为要计算平方根的参数,sqrt()函数返回的是x的平方根,返回值

    2024年01月22日
    浏览(36)
  • FPGA verilog 简单的平方根求法

    用下面的平方根求法不需要乘法,只需简单的移位就能实现。 原理参照论文 A New Non-Restoring Square Root Algorithm and Its VLSI Implementations

    2024年02月04日
    浏览(35)
  • Java习题之实现平方根(sqrt)函数

    目录 前言 二分查找 牛顿迭代法 总结 🎁博主介绍:博客名为tq02,已学C语言、JavaSE,目前学了MySQL和JavaWeb 🎥学习专栏:  C语言        JavaSE      MySQL基础 🎄博主链接:tq02的博客_CSDN博客-C语言,Java,MySQL领域博主         可使用java.lang.Math类的 sqrt(double)方法 求平方根。

    2024年02月15日
    浏览(44)
  • LeetCode每日一题——x 的平方根

    乍一看题目只需要算一个数的平方根,根据我们之前学的C语言我们能很快的想到使用sqrt,pow这类的math.h库函数,但是题目要求我们不能使用,那么我们便可以使用我们的数学思想,将给的整数拆成两个一样的数相乘。 代码实现: 运行结果:   PS:看到这里了,码字不易,给

    2024年03月23日
    浏览(47)
  • 【算法专题突破】二分查找 - x 的平方根(18)

    目录 1. 题目解析 2. 算法原理 3. 代码编写 写在最后: 题目链接:69. x 的平方根 - 力扣(LeetCode) 这道题就是求算数平方根, 要注意的点是他只需要保留整数部分,小数部分会舍去 我们确定好一个区间 1 ~ x,数字 x 的算数平方根一定在这里面, 最简单的思路就是用暴力解法

    2024年02月07日
    浏览(39)
  • 力扣:69. x 的平方根(Python3)

    给你一个非负整数  x  ,计算并返回  x  的  算术平方根  。 由于返回类型是整数,结果只保留  整数部分  ,小数部分将被  舍去 。 注意: 不允许使用任何内置指数函数和算符,例如  pow(x, 0.5)  或者  x ** 0.5  。 来源:力扣(LeetCode) 链接:力扣(LeetCode)官网 - 全球

    2024年02月12日
    浏览(53)
  • 算法leetcode|69. x 的平方根(rust重拳出击)

    给你一个非负整数 x ,计算并返回 x 的 算术平方根 。 由于返回类型是整数,结果只保留 整数部分 ,小数部分将被 舍去 。 注意 :不允许使用任何内置指数函数和算符,例如 pow(x, 0.5) 或者 x ** 0.5 。 0 = x = 2 31 - 1 面对这道算法题目,二当家的再次陷入了沉思。 要开平方,但

    2024年02月13日
    浏览(33)

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

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

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

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

二维码1

领取红包

二维码2

领红包