08-求解Ax=b:可解性和解的结构
一、增广矩阵
假设我们要求解方程$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:可解性和解的结构相关推荐
- MIT线性代数:8.求解Ax=b:可解性和解的结构
1.Ax=b 最后消元得到: 所以可以得出有解的条件是: 1.1Ax=b可解性 (1)从列空间看:b必须满足是A各列的线性组合,当然这也是方程组有解的要求. (2)如果A的各行经过线性组合得到零行,那 ...
- 介绍求解AX=b:可解性与解的结构
前面用高斯消元法或矩阵LU分解求解线性方程组的解,主要是针对有唯一解(矩阵A可逆)的情况,下面进一步介绍线性方程组有多个解的情况下,解的求解. 转载于:https://www.cnblogs.com/ ...
- Matlab求解AX=XB(手眼标定用)
描述 使用matlab求解AX=XB的问题 可以正确求出X的值,我已经成功验证了 至于说X具体代表什么,要结合你自己推导出的公式AX=XB之中,A.X.B分别代表什么 代码 main.m 包含我做手眼 ...
- Matlab --- 求解Ax=b时的反斜杠“\“,backslash
一直很喜欢两个老爷爷,一个是MIT主讲线性代数的Gilbert Strang,另一个就是Matlab的首席数学科学家Cleve Moler.记得作者在几年前,下载了他个人主页上的两本书.一本叫< ...
- opencv求解AX=0
在Ax=0时候是不能用cvsolve来接函数的,但是可以利用一下函数来求解 [cpp] view plaincopyprint? #include<cv.h> #include<io ...
- 求解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 - 求解Ax=0:主变量、特解
求解Ax=0:主变量.特解 求零空间(Nullspace) 矩阵 \(A\) 的零空间即满足 \(Ax=0\) 的所有构成 \(x\) 的向量空间. 对于矩阵 \(A\) 进行"行操作&qu ...
- matlab 求解 Ax=B 时所用算法
x = A\B; x = mldivide(A, B); matlab 在这里的求解与严格的数学意义是不同的, 如果 A 接近奇异,matlab 仍会给出合理的结果,但也会提示警告信息: 如果 A 为 ...
- 扩展欧几里得求解ax+by=c的特殊解(模板)
转自:https://www.cnblogs.com/033000-/p/10040183.html template<class T> void exgcd(T a,T b,T & ...
最新文章
- 正点原子探索者原理图_正点原子【STM32-F407探索者】第六章 跑马灯实验
- 浅谈 instanceof 和 typeof 的实现原理
- 为什么居住的地方离上班的地方很远
- 【算法】漫画:如何找到链表的倒数第n个结点?
- RedHat Linux AS4 LAMP经典网站搭建实例
- 利用VSCode阅读OpenFOAM源代码及其调试Debug【终极总结篇】
- 设置状态栏和标题栏的样式
- 深搜DFS\广搜BFS 图初步入门
- 一个交易平台源码,全源无接口
- 一秒点击手机屏幕次数_抓住夏天的尾巴|与自动点击评论器邂逅一次
- 比尔·盖茨:美国优先的世界观使我担心
- FRR BGP 协议分析 3 --- FSM 状态机
- serialVersionUID快捷键
- 根据城市首字母进行分类,区分多音字,获取城市首字母
- pythonai人脸识别_AI的强大!用Python实现一个简单的人脸识别
- C语言练习---杨辉三角
- matlab的程序设计心得和体会,程序设计心得体会(多篇).doc
- sim7600ce 拨号上网测试_SIM7600CE应用程序调试流程
- 微信小程序有哪些?401~500
- 大数据应用在医疗行业的5个经典案例
热门文章
- [C/CPP系列知识] C++中extern “C” name mangling -- Name Mangling and extern “C” in C++
- 黑马程序员_语句结构及数组
- 总结必须学习的10项.NET技术
- Spring在SSH中的角色和作用
- client_hello_cb、get_session_cb、servername_cb、cert_cb
- 2018-04-08椭圆曲线测试程序
- 简书优秀IT专栏作者推荐
- css字体倾斜角度_css如何实现渐变效果?css背景色渐变与文字渐变效果的实现(代码实例)...
- 地图整饰-框架与格网
- SpringCloud系列之Nacos+Dubbo+Seata应用篇