上一篇文章讲述了Ax=0的解和矩阵A的零空间, 这里我们讨论Ax=b的解以及矩阵A的列空间。

Ax=0是肯定有解的,因为总存在x为全零向量,使得方程组成立。而Ax=b是不一定有解的,我们需要高斯消元来确定。我们还是利用上一篇讲述了Ax=0的解的矩阵A来举例说明:

我们可以得到上述方程组的增广矩阵(等式右侧不是全零向量,消元时值会改变,所以需要用增广矩阵)如下:

然后我们进行高斯消元可以得到:

从上面的矩阵可以看出,等式成立必须有

我们假设一个满足上面条件的b向量,例如:b=[1 5 1+5];并且令两个自由变量x2=0,x4=0,则我们将消元后的矩阵写成方程组的形式如下:

得到的解为:

Xc是这个方程组的一个特解,因为当X2,X4取不同的值时,会得到不同的特解。那么我们如何得到方程的同解呢?即怎样用一般形式来表示所有的特解?
求解Ax=b的过程:
1、求解特解Xc
2、求解Ax=0的解Xn
Ax=b的解就是特解Xc+Xn,证明如下:

Xc我们上面已经得到,Xn在上一篇文章中得到,则通解可以表示为:

至此,我们就得到了Ax=b的解。
通过上面的分析求解,我们知道当b满足下式时,方程组有解:

实际上,方程有解的条件是向量b属于矩阵A的列空间,即向量b可以表示为矩阵A的各列的线性组合。例如上面的例子:

方程的解就是矩阵A中各列前面的系数。
下面推广到更一般的情况,我们以矩阵A的不同情况来看解的结构(假设矩阵A为m*n的矩阵,秩为r):
1、r=n<m,即列满秩(所有列都有主元)
由于所有列都有主元,则自由变量的个数为0,矩阵A的零空间中只有零向量。Ax=b的解的个数为0个或者1个.
举例说明:

当b=[4 3 6 7]时,Ax=b的唯一解为x=[1 1]。
2、r=m<n,即行满秩(所有行都有主元)
由于所有行都有主元,消元后不会出现全为0的行,则Ax=b有无穷多解。且自由变量的个数为n-r,矩阵A的零空间中不只有零向量。
例如:

3、r=m=n,即列、行都满秩(矩阵可逆)
由于列、行都满秩,则具有列满秩,行满秩的一些性质:零空间只有零向量,方程总有解且解唯一。

4、r<m,r<n,非满秩矩阵

Ax=b有无穷多解或则没有解。
从上面的四种情况的讨论,我们可以总结如下:
如果想看一个线性方程组的解的情况,我们可以通过高斯消元法得到矩阵A的最简形式R,R的可能情况如下:

这四种情况分别对应的解的情况为:
1、唯一解或无解
2、无穷多解
3、唯一解
4、无解或无穷多解

原文:http://blog.csdn.net/tengweitw/article/details/40921003

作者:nineheadedbird

【线性代数】线性方程组的求解相关推荐

  1. 第2-2课:线性代数方程组的求解

    多元一次方程组,又称为线性代数方程组.在数值分析领域里有很多算法都会用到线性代数方程组的求解,比如三次样条曲线拟合时用到的插值算法.求解线性代数方程组可以用高斯消元法,高斯消元法是一种代数的方法,其主 ...

  2. 线性代数---线性方程组

    线性代数-线性方程组 常见题型的解题技巧 如果存在什么问题,欢迎批评指正!谢谢!

  3. 重磅!一文读懂线性方程组的求解方法

    目录 1.AAA为方阵 直接法 迭代法 2.AAA为非方阵且A∈Rm×n,m>nA\in R^{m\times n},m>nA∈Rm×n,m>n 2.1. r(A)=n<mr( ...

  4. matlab 病态方程组,数值分析(Hilbert矩阵)病态线性方程组的求解Matlab程序

    (Hilbert 矩阵)病态线性方程组的求解 理论分析表明,数值求解病态线性方程组很困难.考虑求解如下的线性方程组的求解Hx = b ,期中H 是Hilbert 矩阵,()ij n n H h ?=, ...

  5. 按照高等代数的传统解法编写c++程序实现N(N256)元的线性方程组的求解

    对于c++实现线性方程组的求解问题,前辈已经有太多的好算法.有一本书叫<C数值算法>,里面就详尽的论述了如何求解各种各样的线性方程组. 然而,前辈们的算法,无一例外几乎都是只得出一个解.面 ...

  6. 【计算方法】高斯消去法(线性方程组的求解)

    在本章,我们需要了解一些基本的线性代数基础知识.如非奇异矩阵,矩阵的初等变换等. 高斯消元法 我们以一个例子来讲解高斯消元法 那么,我们可以将它抽象为一个增广矩阵, 该矩阵由系数矩阵和y矩阵组成.我们 ...

  7. 线性方程组的求解(C++)

    1. 最佳求解方案 Most efficient way to solve a system of linear equations 求解形如 Ax=bAx=b 的最佳方式 将 A 分解为三角矩阵,A ...

  8. MATLAB数值分析学习笔记:线性代数方程组的求解和高斯消元法

    工程和科学计算的许多基本方程都是建立在守恒定律的基础之上的,比如质量守恒等,在数学上,可以建立起形如 [A]{x}={b} 的平衡方程.其中{x}表示各个分量在平衡时的取值,它们表示系统的状态或响应: ...

  9. 【复习笔记】线性代数——线性方程组

    一.同解方程组 若是同解方程组,则有: 1.的解满足,且的解满足(解出两个解相互代入都成立)定义 2.,且的解满足或 的解满足 (秩相同,单方满足) 3.,三秩相等 二.两个方程组的公共解 1.两个方 ...

  10. 线性代数:线性方程求解、矩阵的逆、线性组合、线性独立

    本文参考www.deeplearningbook.org一书第二章2.3 Identity and Inverse Matrices 2.4 Linear Dependence and Span 本文 ...

最新文章

  1. URL化 替换空格
  2. 在 Keil uVision4 MDK下配置开发STM32F103Z完整教程
  3. C++ Double Ended Queues(双向队列)
  4. 股骨截骨php钢板,股骨远端截骨(DFO)术前设计及手术步骤【附视频】
  5. 启明云端分享|SSD20X_烧录mac地址文档参考
  6. zabbix2.2升级到zabbix3.0.2
  7. 被清华免试录取的围棋天才,横扫60位围棋大师的最强AI,竟然都输给了高中生!?...
  8. JS判断一个数字是否为小数
  9. 现代软件工程 第一章 四则运算的实现--栈实现
  10. navicat连接mysql报错ssh_Navicat 连接SSH方式连接mysql报错 ssh:expected key exchange group packet...
  11. 在python中如何将两个list合并成一个list,不用for语句
  12. 正则表达式的学习使用
  13. Lazy Load Plugin for jQuery延迟加载测试成功
  14. 基于双向长短期记忆神经网络【biLSTM】模型的污染数据预测实战
  15. 车载系统华山论剑:Ali OS、Android、QNX孰优孰劣
  16. python print退格_python退格输入
  17. 手机无线上网何时进入战国时代?
  18. Maya vray 3S皮肤材质球设置与材质节点连接
  19. 心电图心电轴怎么计算_心电图心电轴计算表
  20. B2B商业模式以及交易模式

热门文章

  1. 电路图中电阻分类字母速记说明图文
  2. checkbox不全部选中,会跳出提示(这个是我记录的重点)
  3. 网页怎么看不到css,怎么查看网页的CSS代码?如何阅读?
  4. P2756 飞行员配对方案问题【二分匹配】【网络流Dinic、EK详解】【萌新初写】
  5. 怎么把知网的外文文献翻译成中文_知网查重中文译成英文可行吗?
  6. 怎样压缩ppt的大小?
  7. 形容谣言的四字词语_四字词语加解释大全
  8. 计算机桌面不同步,电脑时间不同步怎么回事 电脑时间不能自动更新如何修复...
  9. Java相关技术文档汇总
  10. duboo仿猫眼微服务架构—微服务入门