一、增广矩阵

 假设我们要求解方程$Ax=b$,其中矩阵$A$和$b$如下所示:

$A = \left[\begin{array}{llll}{1} & {2} & {2} & {2} \\ {2} & {4} & {6} & {8} \\ {3} & {6} & {8} & {10}\end{array}\right]$

$b = \left[\begin{array}{llll}{b_1}\\ {b_2}\\{b_3}\end{array}\right]$

 我们可以在矩阵$A$右侧加上一列b,得到:

$\left[\begin{array}{llll}{1} & {2} & {2} & {2} & {b_1} \\ {2} & {4} & {6} & {8}& {b_2} \\ {3} & {6} & {8} & {10}& {b_3}\end{array}\right]$

此即为增广矩阵:即A和b一块二考虑

二、可解性

 方程$Ax=b$可解的条件是:当$b$属于$A$的列空间时

 换一种说法:如果矩阵$A$各行的线性组合得到零行,那么$b$分量执行相同的操作也必须得到零向量(参考上面的矩阵$A$)

三、求$Ax=b$所有解

 第一步:求一个特解

  将所有自由变量设为0

  我们将增广矩阵进行消元后:

$\begin{array}{l}{\left[\begin{array}{}{1} & {2} & {2} & {2} & {b_{1}} \\ {0} & {0} & {2} & {4} & {b_{2}-2 b_{1}} \\ {0} & {0} & {0} & {0} &{b_{3}-b_{2}-b_{1}}\end{array}\right]} \\ {} \end{array}$

  方程有解的条件便是:$O=b_{3}-b_{2}-b_{1}$,假设$b=\left[\begin{array}{l}{1} \\ {5} \\ {6}\end{array}\right]$,则消元矩阵为:

$\begin{array}{l}{\left[\begin{array}{}{1} & {2} & {2} & {2} & {1} \\ {0} & {0} & {2} & {4} & {3} \\ {0} & {0} & {0} & {0} &{0}\end{array}\right]} \\ {}\end{array}$

  然后求出主变量:

   针对上面的例子,$x_2, x_4 = 0$,则求主变量过程简化为:

$x_{1}+2 x_{3}=1$
$2 x_{3}=3$

   于是一个特解为:$x_p = \left[\begin{array}{c}{-2} \\ {0} \\ {3/2} \\ {0}\end{array}\right]$

  还记得07)节中的零空间吗?可以去看看,于是我们的通解为:即特解和零空间

$X = \left[\begin{array}{c}{-2} \\ {0} \\ {3/2} \\ {0}\end{array}\right]+C_{1}\left[\begin{array}{c}{-2} \\ {1}\\ {0} \\ {0}\end{array}\right]+ C_{2}\left[\begin{array}{c}{2} \\ {0} \\ {-2} \\ {1}\end{array}\right]$

转载于:https://www.cnblogs.com/always-fight/p/11412421.html

08-求解Ax=b:可解性和解的结构相关推荐

  1. MIT线性代数:8.求解Ax=b:可解性和解的结构

    1.Ax=b 最后消元得到: 所以可以得出有解的条件是: 1.1Ax=b可解性 (1)从列空间看:b必须满足是A各列的线性组合,当然这也是方程组有解的要求. (2)如果A的各行经过线性组合得到零行,那 ...

  2. 介绍求解AX=b:可解性与解的结构

    前面用高斯消元法或矩阵LU分解求解线性方程组的解,主要是针对有唯一解(矩阵A可逆)的情况,下面进一步介绍线性方程组有多个解的情况下,解的求解. 转载于:https://www.cnblogs.com/ ...

  3. Matlab求解AX=XB(手眼标定用)

    描述 使用matlab求解AX=XB的问题 可以正确求出X的值,我已经成功验证了 至于说X具体代表什么,要结合你自己推导出的公式AX=XB之中,A.X.B分别代表什么 代码 main.m 包含我做手眼 ...

  4. Matlab --- 求解Ax=b时的反斜杠“\“,backslash

    一直很喜欢两个老爷爷,一个是MIT主讲线性代数的Gilbert Strang,另一个就是Matlab的首席数学科学家Cleve Moler.记得作者在几年前,下载了他个人主页上的两本书.一本叫< ...

  5. opencv求解AX=0

    在Ax=0时候是不能用cvsolve来接函数的,但是可以利用一下函数来求解 [cpp] view plaincopyprint? #include<cv.h> #include<io ...

  6. 求解ax + by = c 这类方程

    基础知识: 1.对于任意的ax+by=c, 如果我们知道有一组解x0, y0; 那么 x1 = x0+kb'(b'=b/gcd(a,b)), y1 = y0-ka'(a'=a/gcd(a,b)); 求 ...

  7. 矩阵论 - 7 - 求解Ax=0:主变量、特解

    求解Ax=0:主变量.特解 求零空间(Nullspace) 矩阵 \(A\) 的零空间即满足 \(Ax=0\) 的所有构成 \(x\) 的向量空间. 对于矩阵 \(A\) 进行"行操作&qu ...

  8. matlab 求解 Ax=B 时所用算法

    x = A\B; x = mldivide(A, B); matlab 在这里的求解与严格的数学意义是不同的, 如果 A 接近奇异,matlab 仍会给出合理的结果,但也会提示警告信息: 如果 A 为 ...

  9. 扩展欧几里得求解ax+by=c的特殊解(模板)

    转自:https://www.cnblogs.com/033000-/p/10040183.html template<class T> void exgcd(T a,T b,T & ...

最新文章

  1. 正点原子探索者原理图_正点原子【STM32-F407探索者】第六章 跑马灯实验
  2. 浅谈 instanceof 和 typeof 的实现原理
  3. 为什么居住的地方离上班的地方很远
  4. 【算法】漫画:如何找到链表的倒数第n个结点?
  5. RedHat Linux AS4 LAMP经典网站搭建实例
  6. 利用VSCode阅读OpenFOAM源代码及其调试Debug【终极总结篇】
  7. 设置状态栏和标题栏的样式
  8. 深搜DFS\广搜BFS 图初步入门
  9. 一个交易平台源码,全源无接口
  10. 一秒点击手机屏幕次数_抓住夏天的尾巴|与自动点击评论器邂逅一次
  11. 比尔·盖茨:美国优先的世界观使我担心
  12. FRR BGP 协议分析 3 --- FSM 状态机
  13. serialVersionUID快捷键
  14. 根据城市首字母进行分类,区分多音字,获取城市首字母
  15. pythonai人脸识别_AI的强大!用Python实现一个简单的人脸识别
  16. C语言练习---杨辉三角
  17. matlab的程序设计心得和体会,程序设计心得体会(多篇).doc
  18. sim7600ce 拨号上网测试_SIM7600CE应用程序调试流程
  19. 微信小程序有哪些?401~500
  20. 大数据应用在医疗行业的5个经典案例

热门文章

  1. [C/CPP系列知识] C++中extern “C” name mangling -- Name Mangling and extern “C” in C++
  2. 黑马程序员_语句结构及数组
  3. 总结必须学习的10项.NET技术
  4. Spring在SSH中的角色和作用
  5. client_hello_cb、get_session_cb、servername_cb、cert_cb
  6. 2018-04-08椭圆曲线测试程序
  7. 简书优秀IT专栏作者推荐
  8. css字体倾斜角度_css如何实现渐变效果?css背景色渐变与文字渐变效果的实现(代码实例)...
  9. 地图整饰-框架与格网
  10. SpringCloud系列之Nacos+Dubbo+Seata应用篇