追赶法求解方程组备忘

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

本篇内容为数值分析中,用追赶法求解方程组的方法,备忘如下:

1. 原理部分

追赶法求解的矩阵格式一般如下:

a1 c1  0  0
b2 a2 c2  0
0  b3 a3 c3
0  0  b4 a4

如果矩阵 A 存在 doolittle 分解,则计算步骤:

  1. 首先需要对矩阵进行 LU 分解,得到两个三对角矩阵 L 和 U。
  2. 然后依次求解 Ly=b 和 Ux=y 两个线性方程组即可得到方程组的解。

L 和 U 的格式如下

      1  0   0  0            q1 c1  0  0
L =   p2 1   0  0   U =      0  q2 c2  0
      0  p3  1  0            0  0  q3 c3
      0  0   p4 1            0  0   0 q4

可以看出,L 对角线元素均为 1;U 中 C1、C2、C3 等都是照抄下来。
优势也会把这个矩阵合并化简成如下格式

q1 c1  0   0 
p2 q2  c2  0
0  p3  q3  c3
0  0   p4  q4

计算规则/步骤为:

  • q1 = a1
  • pi = bi/q(i-1)
  • qi = ai - pi*c(i-1)

更通俗的描述为:

  • q1 为 a1,照抄下来即可
  • pi 等于 bi/q(i-1) 即,“我自己” ➗ “我上面(U 矩阵中对应位置的元素)”
  • qi 等于 ai - pi*c(i-1) 即,原矩阵该位置 - “我左边(L 矩阵中)” ✖️ “我上面(U 矩阵中对应位置)”

2. 例题解答

设系数矩阵如下,则求解过程为:
追赶法解例题及答案,笔记,数值分析,追赶法解矩阵
解答过程如下
追赶法解例题及答案,笔记,数值分析,追赶法解矩阵

参考资料汇总:文章来源地址https://www.toymoban.com/news/detail-766479.html

  • 思路参考
  • 例题矩阵参考

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

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

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

相关文章

  • MATLAB:方程组的求解

    综合实例应用:方程组的求解 无论工程应用问题,还是数学计算问题,方程组都是解决问题转化的重要途径之一,将复杂问题转化为简单的方程组矩阵求解问题。 利用矩阵分解来求解线性方程组,是工程计算中最常用的计算。 LU分解法是先将系数矩阵A进行LU分解,得到LU=P

    2024年01月19日
    浏览(38)
  • 数值分析——线性方程组求解

    清理磁盘的时候偶然发现大二下数值分析的实验作业还在,本着在丢弃之前可以放在网上以备不时之需的原则,我便发了上来。 分别用直接法、Jacobi迭代法、Gauss-Seidel迭代法求解下列线性方程组AX = b,其中A为五对角矩阵(n=20),b是除第一个分量是1外,其他分量都是0的列向量

    2024年02月05日
    浏览(43)
  • MATLAB 之 线性方程组求解

    在 MATLAB 中,关于线性方程组的解法一般分为两类:一类是直接法,就是在没有舍入误差的情况下,通过有限步的矩阵初等运算来求得方程组的解;另一类是迭代法,就是先给定一个解的初始值,然后按照一定的迭代算法进行逐步逼近,求出更精确的近似解。 线性方程组的直

    2024年02月08日
    浏览(51)
  • 10种基于MATLAB的方程组求解方法

              线性方程组的求解包括直接法和迭代法,其中迭代法包括传统的高斯消元法,最速下降法,牛顿法,雅克比迭代法,共轭梯度法,以及智能启发式算法求解法和神经网络学习算法,传统算法可以相互组合改进,智能仿生启发式算法包括粒子群算法,遗传算法,模拟

    2024年02月02日
    浏览(42)
  • matlab使用教程(6)—线性方程组的求解

            进行科学计算时,最重要的一个问题是对联立线性方程组求解。在矩阵表示法中,常见问题采用以下形式:给定两个矩阵 A 和 b,是否存在一个唯一矩阵 x 使 Ax = b 或 xA = b?         考虑一维示例具有指导意义。例如,方程         7x = 21         是否具

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

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

    2024年02月12日
    浏览(39)
  • 40.利用欧拉法求解微分方程组(matlab程序)

    1. 简述        求解微分方程的时候,如果不能将求出结果的表达式,则可以对利用数值积分对微分方程求解,获取数值解。欧拉方法是最简单的一种数值解法。前面介绍过MATLAB实例讲解欧拉法求解微分方程,今天实例讲解欧拉法求解一阶微分方程组。 本文理论部分来自知乎

    2024年02月14日
    浏览(58)
  • 解线性方程组(一)——克拉默法则求解(C++)

    解线性方程组最基础的方法就是使用克拉默法则,需要注意的是,该方程组必须是线性方程组。 假设有方程组如下: { a 11 x 1 + a 12 x 2 + ⋯ + a 1 n x n = b 1 a 21 x 1 + a 22 x 2 + ⋯ + a 2 n x n = b 2 ⋯ ⋯ ⋯ a n 1 x 1 + a n 2 x 2 + ⋯ + a n n x n = b n begin{cases} a_{11}x_1+a_{12}x_2+cdots+a_{1n}x_n=b_1\\\\

    2024年02月20日
    浏览(38)
  • 线性代数代码实现(七)求解线性方程组(C++)

    前言:         上次博客,我写了一篇关于定义矩阵除法并且代码的文章。矩阵除法或许用处不大,不过在那一篇文章中,我认为比较好的一点是告诉了大家一种计算方法,即:若矩阵  已知且可逆,矩阵  已知,并且  ,求解矩阵 B 。我认为这种初等行变换的方法还是挺

    2023年04月23日
    浏览(48)
  • 牛顿(Newton)迭代法求解非线性方程以及方程组的Matlab实现

    必做题目比较简单,写得有些随意,主要还是第二个拓展题目的难度比较高 传入题设数据有: 另附运行截图  

    2024年02月11日
    浏览(49)

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

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

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

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

二维码1

领取红包

二维码2

领红包