当前位置: 首页 > article >正文

矩阵的可解性:关于Ax=b的研究

上一篇文章讲解了如何求解 A x = 0 Ax=0 Ax=0,得到 A A A的零空间。

类似的,我们今天学习的是如何求解 A x = b Ax=b Ax=b,并以此加强你对线性代数中,代数与空间的理解。


同样的,我们举与上一次一样的例子,矩阵 A A A为:

[ 1 2 2 2 2 4 6 8 3 6 8 10 ] \left[ \begin{matrix} 1 & 2 & 2 &2\\ 2 & 4 & 6&8\\ 3 & 6& 8 &10 \end{matrix} \right] 1232462682810

关于这个矩阵的详细分析与消元过程在上一篇文章讨论过,这里就不再赘述。


首先,我们将 b b b增广到矩阵 A A A中,得到如下矩阵:

∣ 1 2 2 2 b 1 2 4 6 8 b 2 3 6 8 10 b 3 ∣ \left| \begin{array}{lccc|c} {1}&{2}&{2}&{2} &{b_1}\\ {2}&{4}&{6}&{8} &{b_2}\\ {3}&{6}&{8}&{10}&{b_3} \end{array} \right| 1232462682810b1b2b3

经消元处理,能得到如下矩阵:

∣ 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 ∣ \left| \begin{array}{lccc|c} {1}&{2}&{2}&{2} &{b_1}\\ {0}&{0}&{2}&{4} &{b_2-2\times b_1}\\ {0}&{0}&{0}&{0}&{b_3-b_2-b_1} \end{array} \right| 100200220240b1b22×b1b3b2b1


在继续进行下一步操作前,让我们想一想这个问题: A x = b Ax=b Ax=b在何时有解?

观察消元过后的第三行,不难发现, b b b的元素应该满足 b 3 − b 2 − b 1 = 0 b_3-b_2-b_1=0 b3b2b1=0,这样才能使矩阵第三行成立。对这个结论进行拓展,不难想到,当 b b b在矩阵 A A A的列空间内时,方程有解。明白这一点也会对我们接下来的操作有指导意义。


如同我们求零空间的方法,我们利用消元过后的自由列能快速得到一个关于 A x = b Ax=b Ax=b的特殊解。

具体到这道题上,我们可以看到 A 1 , 1 与 A 2 , 3 A_{1,1}与A_{2,3} A1,1A2,3为主元。因为自由列的变量可以取任意值,为求计算方便,我们一般取其为0,即 x 2 = 0 , x 4 = 0 x_2=0,x_4=0 x2=0,x4=0

那么此时的方程就变为了这样:

x 1 + 2 x 3 = b 1 2 × x 2 = b 2 − 2 b 1 x_1+2x_3=b_1 \\ 2 \times x_2 = b_2-2b_1 x1+2x3=b12×x2=b22b1

因为 b 1 , b 2 b_1,b_2 b1,b2为参数,所以现在我们就求得了特解 x p a r t i c u l a r , 即 x p x_{particular},即x_p xparticular,xp


又一次同样的,我们采用求零空间时的方法,利用特解来求得所有的解,而这里也会用上零空间 N N N,设其中任意的元素为 n n n吧。

那么,我们有:

A x p = b A n = 0 Ax_p=b \\ An = 0 Axp=bAn=0

不难发现, A ( x p + n ) = b A(x_p+n)=b A(xp+n)=b,即特解加上零空间的和后得到的向量同样是方程的解。不妨猜想,特解加上零空间即使所有的解。前面证明了充分性,下面证明必要性:

x x x为一个任意的方程的解,有
A x = b A x p = b → A ( x − x P ) = 0 Ax=b \\ Ax_p=b \\ \rightarrow A(x-x_P)=0 Ax=bAxp=bA(xxP)=0
换言之 n + x p = x n+x_p=x n+xp=x
证得必要性成立。

所以,我们得到了 A x = b Ax=b Ax=b的解,即为其特解加上 A A A的零空间。

此时,再来想象一下,零空间是经过原点的向量空间,那么 A x = b Ax=b Ax=b的解就应是将零空间向特解的方向平移过去所得。要注意的是,其解并不包含原点,所以不是向量空间。



http://www.kler.cn/news/366776.html

相关文章:

  • 【Python爬虫实战】Selenium自动化网页操作入门指南
  • WASM 使用说明23事(RUST实现)
  • 基于neo4j的新冠治疗和新冠患者轨迹的知识图谱问答系统
  • C++ [项目] 愤怒的小鸟
  • 【云原生】Kubernets1.29部署StorageClass-NFS作为存储类,动态创建pvc(已存在NFS服务端)
  • 分页列表缓存
  • WeakHashMap详解
  • 谷粒商城の订单服务分布式事务
  • 京准电钟:NTP网络授时服务器应用航管自控系统
  • jmeter学习(6)逻辑控制器
  • 【Python】Whoosh:全流程自建搜索引擎
  • 【c++篇】:从基础到实践--c++内存管理技巧与模版编程基础
  • 免费的CMS系统有哪些?
  • windows中的tracert命令
  • Apache Paimon介绍
  • IIS下FTP服务器的PASV端口范围修改方法
  • SQL 干货 | SQL 反连接
  • OBOO鸥柏丨液晶拼接大屏分布式基本管理系统架构显示技术曝光
  • python subproces模块
  • 安装Maven配置以及构建Maven项目(2023idea)
  • 智慧铜矿厂综合管控平台 | 图扑软件
  • 计算机网络的主要知识点小结
  • 反向代理服务器---NGINX
  • 【PHP】在ThinkPHP6中Swoole与FPM的简单性能测试对比
  • FullRankMatrix
  • 【Vulnhub靶场】DC-4