线性方程组

线性方程组

设存在线性方程组
{ a 1 , 1 x 1 + a 1 , 2 x 2 + ⋯ + a 1 , n x n = b 1 a 2 , 1 x 1 + a 2 , 2 x 2 + ⋯ + a 2 , n x n = b 2 ⋮ ⋮ a m , 1 x 1 + a m , 2 x 2 + ⋯ + a m , n x n = b m \left.\left\{\begin{array}{l}a_{1,1}x_1+a_{1,2}x_2+\cdots+a_{1,n}x_n=b_1\\a_{2,1}x_1+a_{2,2}x_2+\cdots+a_{2,n}x_n=b_2\\\vdots\quad\vdots\\a_{m,1}x_1+a_{m,2}x_2+\cdots+a_{m,n}x_n=b_m\end{array}\right.\right. a1,1x1+a1,2x2++a1,nxn=b1a2,1x1+a2,2x2++a2,nxn=b2am,1x1+am,2x2++am,nxn=bm
为了简化书写写成矩阵形式
( a 11 a 12 ⋯ a 1 n b 1 a 21 a 22 ⋯ a 2 n b 2 ⋮ ⋮ ⋱ ⋮ ⋮ a m 1 a m 2 ⋯ a m n b m ) \left.\left(\begin{array}{llll|l}a_{11}&a_{12}&\cdots&a_{1n}&b_{1}\\a_{21}&a_{22}&\cdots&a_{2n}&b_{2}\\\vdots&\vdots&\ddots&\vdots&\vdots\\a_{m1}&a_{m2}&\cdots&a_{mn}&b_{m}\end{array}\right.\right) a11a21am1a12a22am2a1na2namnb1b2bm
所以行数为方程组的个数,列数为未知数的个数

高斯消元法

主元不能为0,主元要消去主元下方的元素

高斯消元法时间复杂度需要 n 3 3 + n 2 − n 3 \frac{n^3}3+n^2-\frac n3 3n3+n23n次乘除法以及 n 3 3 + n 2 2 − 5 n 6 \frac{n^3}3+\frac{n^2}2-\frac{5n}6 3n3+2n265n次加减法

Gauss-Jordan method

Gauss-Jordan 方法在高斯消元法的基础上增加了两个规则

  1. 主元必须是1
  2. 主元上方和下方的所有的项都应被消去


( a 11 a 12 ⋯ a 1 n b 1 a 21 a 22 ⋯ a 2 n b 2 ⋮ ⋮ ⋱ ⋮ ⋮ a n 1 a n 2 ⋯ a n n b n ) Gauss-Jordan method → ( 1 0 ⋯ 0 s 1 0 1 ⋯ 0 s 2 ⋮ ⋮ ⋱ ⋮ ⋮ 0 0 ⋯ 1 s n ) \left.\left(\begin{array}{llll|l}a_{11}&a_{12}&\cdots&a_{1n}&b_{1}\\a_{21}&a_{22}&\cdots&a_{2n}&b_{2}\\\vdots&\vdots&\ddots&\vdots&\vdots\\a_{n1}&a_{n2}&\cdots&a_{nn}&b_{n}\end{array}\right.\right) \underrightarrow{\text{Gauss-Jordan method}} \left.\left(\begin{array}{cccc|c}1&0&\cdots&0&s_1\\0&1&\cdots&0&s_2\\\vdots&\vdots&\ddots&\vdots&\vdots\\0&0&\cdots&1&s_n\end{array}\right.\right) a11a21an1a12a22an2a1na2nannb1b2bn Gauss-Jordan method 100010001s1s2sn
其中Gauss-Jordan method的时间复杂度需要 n 3 2 + n 2 2 \frac{n^3}2+\frac{n^2}{2} 2n3+2n2次乘除法以及 n 3 3 − n 2 \frac{n^3}3-\frac{n}2 3n32n次加减法

image-20231115103808538

部分主元法

在选择主元法过程中,在候选主元位置所在的列选择绝对值最大的数字作为主元,若不在主元位置则交换位置使其在主元的位置上

在部分主元法中只涉及到了行交换

完全主元法

在完全主元法过程中,在候选主元位置选择整个系数矩阵中绝对值最大的数字作为主元,若不在主元位置则交换位置使其在主元的位置上

在部分主元法中不只涉及到了行交换而且涉及到了列交换

修改后的高斯消元法

因为高斯消元法存在主元无法选择的可能,所以提出了修改


( 1 2 1 3 3 2 4 0 4 4 1 2 3 5 5 2 4 0 4 7 ) ⟶ ( 1 2 1 3 3 0 0 − 2 − 2 − 2 0 0 2 2 2 0 0 − 2 − 2 1 ) \begin{pmatrix}1&2&1&3&3\\2&4&0&4&4\\1&2&3&5&5\\2&4&0&4&7\end{pmatrix}\longrightarrow\begin{pmatrix}1&2&1&3&3\\0&0&-2&-2&-2\\0&0&2&2&2\\0&0&-2&-2&1\end{pmatrix} 12122424103034543457 10002000122232223221

矩阵的秩=矩阵的主元的个数=行阶梯型的非零行的个数=矩阵基本列的个数

其中矩阵的基本列卫矩阵主元所在的列

行阶梯型

主元下的元素为0,主元左侧的元素为0

行最简型

首先是行阶梯型,然后主元所在的列只有主元不为0,且主元为1

所以非主元所在的列可以由左边的主元所在的列线性组合而成,即
E ∗ k = μ 1 E ∗ b 1 + μ 2 E ∗ b 2 + ⋯ + μ j E ∗ b j = μ 1 ( 1 0 ⋮ 0 ⋮ 0 ) + μ 2 ( 0 1 ⋮ 0 ⋮ 0 ) + ⋯ + μ j ( 0 0 ⋮ 1 ⋮ 0 ) = ( μ 1 μ 2 ⋮ μ j ⋮ 0 ) \begin{aligned} \mathbf{E}_{*k}& \begin{aligned}=\mu_1\mathbf{E}_{*b_1}+\mu_2\mathbf{E}_{*b_2}+\cdots+\mu_j\mathbf{E}_{*b_j}\end{aligned} \\ &=\mu_1\begin{pmatrix}1\\0\\\vdots\\0\\\vdots\\0\end{pmatrix}+\mu_2\begin{pmatrix}0\\1\\\vdots\\0\\\vdots\\0\end{pmatrix}+\cdots+\mu_j\begin{pmatrix}0\\0\\\vdots\\1\\\vdots\\0\end{pmatrix}=\begin{pmatrix}\mu_1\\\mu_2\\\vdots\\\mu_j\\\vdots\\0\end{pmatrix} \end{aligned} Ek=μ1Eb1+μ2Eb2++μjEbj=μ1 1000 +μ2 0100 ++μj 0010 = μ1μ2μj0
E ∗ k \mathbf{E}_{*k} Ek为非主元列, E ∗ b i \mathbf{E}_{*b_i} Ebi E ∗ k \mathbf{E}_{*k} Ek左边的主元所在的列

a 11 x 1 + a 12 x 2 + ⋯ + a 1 n x n = 0 , a 21 x 1 + a 22 x 2 + ⋯ + a 2 n x n = 0 , ⋮ a m 1 x 1 + a m 2 x 2 + ⋯ + a m n x n = 0. \begin{aligned}a_{11}x_1+a_{12}x_2+\cdots+a_{1n}x_n&=0,\\a_{21}x_1+a_{22}x_2+\cdots+a_{2n}x_n&=0,\\\vdots\\a_{m1}x_1+a_{m2}x_2+\cdots+a_{mn}x_n&=0.\end{aligned} a11x1+a12x2++a1nxna21x1+a22x2++a2nxnam1x1+am2x2++amnxn=0,=0,=0.
像这样的方程组被称为齐次方程组(homogeneous Systems)

当等号右边存在非零值时,被称为非齐次方程组(nonhomogeneous Systems)

x 1 = x 2 = ⋯ = x n = 0 x_1=x_2=\dots=x_n=0 x1=x2==xn=0,则被称为平凡解

基本列位置的未知数称为基本变量(basic variables),对应于非基本列位置的未知数称为自由变量(free variables)

其中矩阵 A A A为齐次方程组时,其中 r a n k ( A ) = r rank(A)=r rank(A)=r,通解为 x = x f 1 h 1 + x f 2 h 2 + ⋯ + x f n − r h n − r x=x_{f1}h_1+x_{f2}h_2+\cdots+x_{fn-r}h_{n-r} x=xf1h1+xf2h2++xfnrhnr,其中 x f 1 , x f 2 , … , x f n − r x_{f 1}, x_{f 2}, \ldots, x_{f n-r} xf1,xf2,,xfnr为自由变量。

其中矩阵 A A A为非齐次方程组时,其中 r a n k ( A ) = r < n rank(A)=r<n rank(A)=r<n,通解为 x = x f 1 h 1 + x f 2 h 2 + ⋯ + x f n − r h n − r + ξ x=x_{f1}h_1+x_{f2}h_2+\cdots+x_{fn-r}h_{n-r}+\xi x=xf1h1+xf2h2++xfnrhnr+ξ,其中 ξ \xi ξ为A的一个特解, x f 1 h 1 + x f 2 h 2 + ⋯ + x f n − r h n − r x_{f1}h_1+x_{f2}h_2+\cdots+x_{fn-r}h_{n-r} xf1h1+xf2h2++xfnrhnr为A为齐次方程组时候的通解。

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/194065.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

windows排除故障工具RELI、PSR

RELI 在windows10中&#xff0c;如果按某些书上或网上所属&#xff0c;点击开始按钮&#xff0c;然后再搜书欧昂立输入RELI&#xff0c;在开始菜单中可以看到Reliable History&#xff08;查看可靠性历史记录&#xff09;蓝色的图标。 可惜&#xff0c;我试验后没能找到。 这…

年货FPS大作,艾尔莎EA B450M-E和你玩转《使命召唤20》

说到动视旗下的《使命召唤》系列&#xff0c;相信大家都不陌生&#xff0c;它以出色爽快的游戏体验以及精良的画面著称&#xff0c;而且每年一部的更新节奏也是如今为数不多的“年货”游戏之一了。时至今日&#xff0c;该系列已经来到了第20部作品&#xff0c;也就是《使命召唤…

Sectigo证书——全球SSL证书第一品牌

随着互联网技术的发展和网络安全意识的提升&#xff0c;SSL证书已经成为确保网络数据传输安全的重要工具。而在众多SSL证书品牌中&#xff0c;Sectigo无疑是其中的佼佼者&#xff0c;被誉为“全球证书第一品牌”。 Sectigo是一个有着超过二十年行业经验的老牌认证机构&#xff…

矩阵的QR分解

矩阵的QR分解 GramSchmidt 设存在 B { x 1 , x 2 , … , x n } \mathcal{B}\left\{\mathbf{x}_{1},\mathbf{x}_{2},\ldots,\mathbf{x}_{n}\right\} B{x1​,x2​,…,xn​}在施密特正交化过程中 q 1 x 1 ∣ ∣ x 1 ∣ ∣ q_1\frac{x_1}{||x_1||} q1​∣∣x1​∣∣x1​​ q k …

心怀祖国放眼世界 爱国人士华国中应邀参加美国旧金山2023(APEC)峰会

据相关媒体美国旧金山报道:2023亚太经合组织&#xff08;APEC&#xff09;领导人非正式会议将于11月15日至17日在美国旧金山召开。11月11日&#xff0c;本年度APEC高级财政官员和部长会晤在旧金山率先启动&#xff0c;APEC CEO峰会将于11月14日至16日开幕。著名爱国人士、亚太一…

基于鼠群算法优化概率神经网络PNN的分类预测 - 附代码

基于鼠群算法优化概率神经网络PNN的分类预测 - 附代码 文章目录 基于鼠群算法优化概率神经网络PNN的分类预测 - 附代码1.PNN网络概述2.变压器故障诊街系统相关背景2.1 模型建立 3.基于鼠群优化的PNN网络5.测试结果6.参考文献7.Matlab代码 摘要&#xff1a;针对PNN神经网络的光滑…

量化交易:借助talib使用技术分析指标

什么是技术分析&#xff1f; 所谓股票的技术分析&#xff0c;是相对于基本面分析而言的。基本分析法着重于对一般经济情况以及各个公司的经营管理状况、行业动态等因素进行分析&#xff0c;以此来研究股票的价值&#xff0c;衡量股价的高低。而技术分析则是透过图表或技术指标…

资深测试总结,现在软件测试有未来吗?“你“的底气在哪里?

目录&#xff1a;导读 前言一、Python编程入门到精通二、接口自动化项目实战三、Web自动化项目实战四、App自动化项目实战五、一线大厂简历六、测试开发DevOps体系七、常用自动化测试工具八、JMeter性能测试九、总结&#xff08;尾部小惊喜&#xff09; 前言 1、为什么会有 “…

图像分类(四) 全面解读复现GoogleNet_InceptionV1-V4

论文解读 InceptionV1 前言 论文题目: Going Deeper with Convolutions Googlenet论文原文地址:https://arxiv.org/pdf/1409.4842.pdf 之前看过VGG的论文&#xff08;VGG精读直达&#xff09;。当时VGG获得了 2014 ILSVRC 图像分类的第二名&#xff0c;今天来看一下第一名…

盘点54个Python实用工具源码Python爱好者不容错过

盘点54个Python实用工具源码Python爱好者不容错过 学习知识费力气&#xff0c;收集整理更不易。 知识付费甚欢喜&#xff0c;为咱码农谋福利。 链接&#xff1a;https://pan.baidu.com/s/1OXyEh-Yy3JI90jvn6d6wRw?pwd8888 提取码&#xff1a;8888 项目名称 7z辅助破解工…

【win32_001】win32命名规、缩写、窗口

整数类型 bool类型 使用注意&#xff1a; 一般bool 的false0&#xff1b;true1 | 2 | …|n false是为0&#xff0c;true是非零 不建议这样用&#xff1a; if (result TRUE) // Wrong! 因为result不一定只返回1&#xff08;true&#xff09;&#xff0c;当返回2时&#xff0c…

【C刷题】day7

&#x1f3a5; 个人主页&#xff1a;深鱼~&#x1f525;收录专栏&#xff1a;【C】每日一练&#x1f304;欢迎 &#x1f44d;点赞✍评论⭐收藏 一、选择题 1、以下对C语言函数的有关描述中&#xff0c;正确的有【多选】&#xff08; &#xff09; A: 在C语言中&#xff0c;一…