线性方程组的直接解法

已于 2023-09-28 17:25:18 修改
阅读量90 收藏
点赞数
文章标签： 算法 机器学习 人工智能
于 2023-09-28 10:51:06 首次发布
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/Cherry_csc/article/details/133376439
版权

高斯消去法

消元法的基本思想：利用第一个方程消去后面n-1个方程中的第一个未知数，再利用第二个方程消去后面n-2个方程的第二个未知数，以此类推，直到把系数矩阵变成上三角阵的形式

然后从最后一个方程开始，解出x_n,代入前一个方程解出x_{n-1}，以此类推。

顺序消元和回代

 算法：Gauss顺序消元算法：

 

 计算复杂度

顺序消元能进行至第k步充要条件

定理：高斯顺序消元法能进行到第k步\Leftrightarrow前面k个顺序主子式都不为零

 矩阵的三角分解

由消元过程的矩阵形式及上述定理，可得到矩阵三角分解定理

存在性：高斯顺序消元已给出

唯一性：归纳+反证法

 三角分解常用语需解多个右端向量方程组情况。解一次三角分解，多次回代，回代复杂度O(n^2)低于分解复杂度O(n^3)

顺序消元法存在问题

列主元高斯消去法

 

 列主元高斯消去法算法

用矩阵语言来描述上述算法，可写成：

 

 证明

 Gauss-Jordan消去法

消元时同时消去对角线上下元素，把A化成对角阵，可以计算A^{-1},但是计算量会增大到\frac{n^3}{2}+O(n^2)

 矩阵三角分解法——Gauss消去法的变形

Doolittle三角分解（LU,LR分解）

利用矩阵的三角分解来求解Ax=b,假设A的顺序主子式都不为零，将A分解为A=LU,再解LUx=b，分为两步：

Ly=b

Ux=y

选主元LU分解

 对称阵的三角分解

 

证明：已经证过A可唯一分解为LU +A^T=A

 什么样的对称阵有Cholesky分解

 非正定-->D矩阵对角元存在负值，无法开根

使用Cholesky分解计算量：\frac{n^3}{6}+O(n^2)

需不需要选主元？

 算法缺陷
——需要计算很多平方根，平方根的计算远大于乘除法。

改进：用A=LDL^T分解来计算

 

带状矩阵的三角分解法

 

 

证明：归纳法

 

变形1——分块

 变形2——循环三对角阵

Cherry_csc
关注
  • 0
    点赞
  • 踩
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
3 线性方程组直接解法.ppt
10-31
67页ppt，目录如下： §2.0 概述 §2.1 消去法 一、Gauss消去法 二、Gauss列主元消去法 §2.2 直接三角分解法 一、消去法与系数矩阵的三角分解 二、Doolittle分解
线性方程组直接解法_furnituren7x_线性方程组解法_
09-29
1、求解线性方程组的Gauss消去法、LU分解法等直接解法。2、编制程序实现Gauss消去法和LU分解法，体会功能模块结构化程序设计方法的优点；3、通过三对角形线性方程组的求解，稀疏线性方程组解法的特点
线性方程组的直接解法_线性方程组的直接解法_
09-28
线性方程组的直接解法:高斯消元法、LU分解法
线性方程组直接解法相关论文研读1
08-03
线性方程组直接解法相关论文研读1
解线性方程组的直接解法
12-17
第三章 解线性方程组的直接解法
Java 7大排序
m0_74270127的博客
05-10 706
插入排序、希尔排序、选择排序、堆排序、冒泡排序、快速排序、归并排序
算法基础01一快速排序，归并排序，二分
m0_74340589的博客
05-06 146
假设两个有效序列 两个指针指向开头 新数组来存答案。1.快速 排序 基于分治。
1080：余数相同问题
如果没有特别说明，c++都是基于GCC-4.9.2运行的
05-11 222
【代码】：1080：余数相同问题
找第二大的数（分治法）
2201_75525567的博客
05-11 210
这段代码通过递归的方式，将数组分为两半，并分别找出左半部分和右半部分的两个最大数，然后比较这四个数的大小，得到第二大的数。这段代码是用来找出数组中第二大的数的。
Golang | Leetcode Golang题解之第67题二进制求和
weixin_66442839的博客
05-04 663
Golang | Leetcode Golang题解之第67题二进制求和
StarryCoding 算法小白周赛2 题解与代码（含视频题解）
qq_29495615的博客
05-04 973
因为他会按照“拜访时间安排表”的顺序拜访每一位好友，而无法成功拜访某一好友时，必定是出现靠后拜访的时间早于靠前拜访的，或者两次拜访的时间相同，所以问题的本质就是判断数组是否。），应该删除这个数还是前一个数，因为我们想要的是让后面的数的合法值域更大，所以我们需要最小化留下的那个数，也就是说，如果靠后的数。如果某个数比前一个数小，那么肯定不行，我们需要考虑删去这个元素。如果大于，则最好不要更改，给后面的元素更多的空间。由于数的范围很小，假如这个数合法，我们直接从此数。,能保证加的数量是最少的，并且。
从零学算法6
m0_53256503的博客
05-11 691
之后，你的输出需要从左往右逐行读取，产生出一个新的字符串，比如：“PAHNAPLSIIGYIR”。将一个给定字符串 s 根据给定的行数 numRows ，以从上往下、从左到右进行 Z 字形排列。输入：s = “PAYPALISHIRING”, numRows = 3。输入：s = “PAYPALISHIRING”, numRows = 4。输入：s = “A”, numRows = 1。输出：“PAHNAPLSIIGYIR”输出：“PINALSIGYAHRPI”
【算法基础】第六章：贪心
m0_65787507的博客
05-07 890
给定 N个闭区间 [a,b]，请你在数轴上。输出选择的点的最小数量。位于区间端点上的点也算作区间内。① 所有区间按右端点从小到大排序② 遍历每一个区间，如果当前区间的左与前一个区间的右有交集，则只需要一个点就可以覆盖掉两个区间。
AtCoder Beginner Contest 353
最新发布
2301_78527293的博客
05-11 635
abc
代码随想录算法训练营第四十六天
weixin_41902984的博客
05-08 247
昨天媳妇发烧啦，和他一顿激烈讨论该不该母乳，目前一切良好，题目在单位做完了，但是有个地方想在家里测试一下。139.单词拆分要注意先背包后物品，并且物品是变化的，不是直接用wordDict这里想试个东西，总感觉s.substr(j-i)的长度不太对。这里真的用的好巧啊，例子用applepenapple。5apple8pen13apple想说的就是，正好dp[x]中的x，是下一个起始位，所以没有那么多的±1来对齐的问题。最后直接看dp[s.size()]即可。
算法练习day4
rempter的博客
05-08 908
算法练习第四天，是链表方面的练习
解双曲型非线性方程的Harden-Yee算法(TVD格式)
程序员的博客
05-11 119
解双曲型非线性方程的Harden-Yee算法。先贴代码，教程后面有空再写。
【算法】并查集
Eristic0618的博客
05-08 878
算法学习笔记之并查集
[算法面试]_01_L1和L2正则化，为什么L1正则化更容易导致稀疏？
qq_41915623的博客
05-11 390
L1和L2都是是用来控制模型复杂度和防止过拟合的常见手段，其中L1正则化是加上参数的绝对值的和，而L2正则化是加上参数的平方和。
matlab做线性方程组的直接解法
06-08
在 中，可以使用“\”操作符进行线性方程组的直接解法，包括高斯消元法和 LU 分解法。 1. 高斯消元法： 假设有一个系数矩阵 A 和一个常数向量 b，可以使用 MATLAB 的“\”操作符进行高斯消元法求解线性方程组。具体的步骤如下： a. 构造增广矩阵 [A, b]。 b. 使用“\”操作符求解线性方程组，即 x = A \ b。 c. 得到解向量 x。 例如，对于下面的线性方程组： $$ \begin{cases} 2x+y+z=5 \\ 4x+6y+2z=18 \\ -2x+7y+2z=7 \end{cases} $$ 可以使用以下 MATLAB 代码求解： ``` A = [2 1 1; 4 6 2; -2 7 2]; b = [5; 18; 7]; x = A \ b; ``` 运行结果为： ``` x = 1.0000 2.0000 -10.0000 ``` 即解向量为 (1, 2, -10)。 2. LU 分解法： 使用 MATLAB 的“lu”函数进行 LU 分解，然后使用“\”操作符求解线性方程组。具体的步骤如下： a. 对系数矩阵 A 进行 LU 分解，即 [L, U, P] = lu(A)。 b. 将方程组转化为 Ly = Pb 和 Ux = y 两个方程组。 c. 使用“\”操作符求解方程组，即 y = L \ (P*b)，x = U \ y。 d. 得到解向量 x。 例如，对于下面的线性方程组： $$ \begin{cases} 2x+y+z=5 \\ 4x+6y+2z=18 \\ -2x+7y+2z=7 \end{cases} $$ 可以使用以下 MATLAB 代码求解： ``` A = [2 1 1; 4 6 2; -2 7 2]; b = [5; 18; 7]; [L, U, P] = lu(A); y = L \ (P * b); x = U \ y; ``` 运行结果为： ``` x = 1.0000 2.0000 -10.0000 ``` 即解向量为 (1, 2, -10)。 以上就是在 MATLAB 中使用“\”操作符进行线性方程组的直接解法，包括高斯消元法和 LU 分解法的实现步骤。

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

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

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

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

抵扣说明：

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

余额充值