矩阵方程求解方法

矩阵方程求解方法

本文所述的矩阵方程是指形如Ax=b的方程,其中A是一个mxn的矩阵,称为方程的系数矩阵。x和b是mx1的矩阵。特别的,当b=0时,这种方程又称为其次方程。本文将讨论这种矩阵的有解条件和求解方法。

矩阵方程的有解条件

为了解释矩阵方程的有解条件,我们首先要熟悉一些概念。

一个矩阵方程的增广矩阵是系数矩阵A和b并在一起构成的矩阵,记作(A,b)。 假定 , ,则矩阵方程的增广矩阵就是

矩阵的秩定义为其行向量中极大线性无关组中包含向量的个数,等价的说法是,矩阵的秩是r,则矩阵通过行列初等变换,变换成左上角是一个r阶单位矩阵,其他都是0的矩阵。矩阵A的秩记作r(A),其中r是英文单词rank的缩写。

有了这两个基本概念,我们就可以准确描述矩阵方程的有解条件了:矩阵方程Ax=b的有解条件是矩阵A的秩等于增广矩阵(A,b)的秩,也就是r(A)=r(A,b)。

证明很简单,既然矩阵A的秩是r,那么肯定可以找到两个可逆的矩阵P,Q,满足

--1)

其中Ir表示r阶单位矩阵。

应用到原来的方程,可以得到:

--2)

我们把Q-1x当作一个未知的变量,PAQ当作系数,这就构成一个新的矩阵方程。而这个矩阵方程的左侧系数除了前r行是有1的之外,其余行是0。为了它有解,Pb的后m-r行必须也是0。这样(A,b)的秩必然是r。

必须注意到Q-1是可逆的,因此以Q-1x为未知变量的方程有解意味着以x为未知变量的原方程也是有解的。

矩阵方程的解

对于矩阵方程Ax=b,如果满足r(A)=r(A,b),则矩阵方程是有解的。为了求它的解,我们首先把矩阵方程通过行列初等变换变化成前文2)式的形式,代入1)式后得到:

--3)

其中Q-1x和Pb是一个列向量,我们可以把它们分割成rx1和(n-r)x1的两个矩阵,分别记作x’1和x’2,及b’1和b’2。则很显然我们可以得到:

--4)

很显然,b’2必须为0,因为展开后b’2等于0 x’1 +0 x’2 =0

而由4式可以看出,x’1= b’1,x’2可以为任意向量。

所以方程最后的解为:

--5)

从解的形式可以看出解空间有如下特性:

1. 方程Ax=b的解空间的秩是n=r(A)

2. 如果A是满秩的,则方程的解唯一。

矩阵方程求解方法

本文所述的矩阵方程是指形如Ax=b的方程,其中A是一个mxn的矩阵,称为方程的系数矩阵。x和b是mx1的矩阵。特别的,当b=0时,这种方程又称为其次方程。本文将讨论这种矩阵的有解条件和求解方法。

矩阵方程的有解条件

为了解释矩阵方程的有解条件,我们首先要熟悉一些概念。

一个矩阵方程的增广矩阵是系数矩阵A和b并在一起构成的矩阵,记作(A,b)。 假定 , ,则矩阵方程的增广矩阵就是

矩阵的秩定义为其行向量中极大线性无关组中包含向量的个数,等价的说法是,矩阵的秩是r,则矩阵通过行列初等变换,变换成左上角是一个r阶单位矩阵,其他都是0的矩阵。矩阵A的秩记作r(A),其中r是英文单词rank的缩写。

有了这两个基本概念,我们就可以准确描述矩阵方程的有解条件了:矩阵方程Ax=b的有解条件是矩阵A的秩等于增广矩阵(A,b)的秩,也就是r(A)=r(A,b)。

证明很简单,既然矩阵A的秩是r,那么肯定可以找到两个可逆的矩阵P,Q,满足

--1)

其中Ir表示r阶单位矩阵。

应用到原来的方程,可以得到:

--2)

我们把Q-1x当作一个未知的变量,PAQ当作系数,这就构成一个新的矩阵方程。而这个矩阵方程的左侧系数除了前r行是有1的之外,其余行是0。为了它有解,Pb的后m-r行必须也是0。这样(A,b)的秩必然是r。

必须注意到Q-1是可逆的,因此以Q-1x为未知变量的方程有解意味着以x为未知变量的原方程也是有解的。

矩阵方程的解

对于矩阵方程Ax=b,如果满足r(A)=r(A,b),则矩阵方程是有解的。为了求它的解,我们首先把矩阵方程通过行列初等变换变化成前文2)式的形式,代入1)式后得到:

--3)

其中Q-1x和Pb是一个列向量,我们可以把它们分割成rx1和(n-r)x1的两个矩阵,分别记作x’1和x’2,及b’1和b’2。则很显然我们可以得到:

--4)

很显然,b’2必须为0,因为展开后b’2等于0 x’1 +0 x’2 =0

而由4式可以看出,x’1= b’1,x’2可以为任意向量。

所以方程最后的解为:

--5)

从解的形式可以看出解空间有如下特性:

1. 方程Ax=b的解空间的秩是n=r(A)

2. 如果A是满秩的,则方程的解唯一。


相关内容

  • 预处理共轭梯度法求解线性方程组
  • [摘 要]针对共轭梯度法求解线性方程组,提出一种预处理思想.基于次思想,首先给出预处理矩阵,然后求解预处理线性方程组,再使用共轭梯度法求解.最后通过几个数值试验,与直接使用共轭梯度法求解线性方程组相比较,本文的方法提高了收敛速度. [关键词]线性方程组,预处理,共轭梯度法 中图分类号:E911 文献 ...

  • 求解三对角线性方程组两类并行算法的特点
  •  一、概述   三对角线性方程组的求解是许多科学和工程计算中最重要也是最基本的问题之一。在核物理、流体力学、油藏工程、石油地震数据处理及数值天气预报等许多领域的大规模科学工程和数值处理中都会遇到三对角系统的求解问题。很多三对角线性方程组的算法可以直接推广到求解块三对角及带状线性方程组。由于在理论和实 ...

  • 浅谈线性方程组的求解及其应用
  • 浅谈线性方程组的求解及其应用 [摘要] 线性代数是代数学的一个重要组成部分,广泛应用于现代科学的许多 分支.其核心问题之一就是线性方程组的求解问题.本文先简要介绍了线性方程 组求解的历史,然后给出线性方程组解的结构.重点介绍了解线性方程组的几种 方法:消元法,克拉默法则求解线性方程组的方法.最后介绍 ...

  • 现代控制理论
  • 现代控制理论知识结构 (研究的顺序符合控制理论的一般步骤:系统建模.模型求解.提出性能指标.系统设计校正.) 一.系统的状态空间建模 1. 从2个角度4种方法建立状态空间表达式. 2. 状态空间表达式的线性变换 3. 由状态空间表达式反推传递函数 4. 建立离散系统的状态空间表达式 5. 需要注意的 ...

  • 线性方程组的直接解法
  • 线性方程组的直接解法 §1线性方程组的条件数 我们曾对一个矩阵计算问题的病态性这个概念作了概略的定性说明.现在我们来考虑如何对线性方程组的病态程度做出定量估计. 设ARnn是非奇异的,bR.我们来考虑线性方程组 Axb, (1) n 之解x对数据A和b的微小扰动的敏感程度.为此,考虑如下的含 ...

  • 材料成型计算机模拟(纯手工打造)
  • 一.名词解释 1计算机模拟的概念:根据实际体系在计算机上进行模拟实验,通过将模拟结果与实际体系的实验数据进行比较,可以检验模型的准确性,也可以检验由模型导出的解析理论作为所作的简化近似是否成功.1 2材料设计是指(主要包含三个方面的含义):理论计算→预报→组分.结构和性能:理论设计→订做→新材料:按 ...

  • 计算方法实验报告_列主元高斯消去法
  • 计 算 方 法 实 验 报 告 动力与机械学院 08级自动化3班 唐禹 [1**********]78 2010.11.08 实验一 :使用列主元消去法求解线性方程组 列主元消去法是在Gauss消去法的基础上改进而得的一种比较快速和合理的求解线性方程组的方法.它的主要思路是通过对每次消元过程中主元的 ...

  • 实验八matlab支路电流法.回路分析法求解电压和电流(1)
  • 实验八 Matlab 支路电流法.回路分析法求解电压和电流 一.实验目的 1.加深理解电阻电路的分析方法,并求解出电压和电流. 2.验证支路电流法.回路分析法,加深对支路电流法.回路分析法的理解. 3.掌握Matlab中矩阵和数组的运算. 4.学习分析Matlab中简单的数值计算. 二.实验原理 电 ...

  • 线性方程组的直接法和迭代法
  • 线性方程组的直接法 直接法就是经过有限步算术运算,无需迭代可直接求得方程组精确解的方法. 线性方程组迭代法 迭代法就是用某种极限过程去逐步逼近线性方程组精确解的方法.该方法具有对计算机的存贮单元需求少,程序设计简单.原始系数矩阵在计算过程中不变等优点,是求解大型稀疏矩阵方程组的重要方法.迭代法不是用 ...