當前位置：首頁 > 编程资源 > 编程问答 >内容正文

编程问答

4.1 简单方程的解

發布時間：2023/12/20 编程问答 30 豆豆

生活随笔收集整理的這篇文章主要介紹了 4.1 简单方程的解小編覺得挺不錯的,現在分享給大家,幫大家做個參考.

前面章節介紹了向量組和矩陣理論，利用這些理論可以解決線性方程 $Amnx=bA_{mn}\mathbf{x}=\mathbf{b}$ 解的存在性和唯一性問題。向量組理論如下：向量 $b\mathbf{b}$ 能被矩陣 $A$ 向量組表示時，則存在解，否則不存在解。矩陣 $A$ 向量組的極大無關組是基時，必存在解。存在解時，如果矩陣 $A$ 向量組是無關組，則解唯一，否則解無窮多。矩陣理論如下： $(A,\mathbf{b})=rank A$ 時存在解，否則不存在解。 $r a n k A = m$ 必存在解。存在解時，如果 $r a n k A = n$ 解唯一，否則解無窮多。

雖然理論上完美地解決了方程解的存在性和唯一性問題，但實際中，如何計算出解，求矩陣的秩，得到向量組的極大無關組等具體問題，都需要通過算法得到，高斯消元法即是這樣的算法，能解決這些具體問題。

4.1 簡單方程的解

先以大家最熟悉的可逆方陣介紹方程的解，此時解存在且唯一。方程 $Annx=bA_{nn}\mathbf{x}=\mathbf{b}$ 按矩陣運算規則展開，即可得熟悉的 $n$ 元一次方程。
$Ax=[a11a12?a1na21a22?a2n?an1an2?ann][x1x2?xn]=[a11x1+a12x2+?+a1nxna21x1+a22x2+?+a2nxn?an1x1+an2x2+?+annxn]=[b1b2?bn]A\mathbf{x}= \left[ \begin{matrix} a_{11} & a_{12} & \cdots & a_{1n}\\ a_{21} & a_{22} & \cdots & a_{2n}\\ \vdots \\ a_{n1} & a_{n2} & \cdots & a_{nn} \end{matrix} \right] \left[ \begin{matrix} x_1 \\ x_2 \\ \vdots \\ x_n \end{matrix} \right]= \left[ \begin{matrix} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n\\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n\\ \vdots \\ a_{n1}x_1 + a_{n2}x_2 + \cdots + a_{nn}x_n \end{matrix} \right]= \left[ \begin{matrix} b_1 \\ b_2 \\ \vdots \\ b_n \end{matrix} \right]$
得 $n$ 元一次方程：
$[a11x1+a12x2+?+a1nxn=b1a21x1+a22x2+?+a2nxn=b2?an1x1+an2x2+?+annxn=bn]\left[ \begin{matrix} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n=b_1\\ a_{21}x_1 + a_{22}x_2 + \cdots + a_{2n}x_n=b_2\\ \vdots \\ a_{n1}x_1 + a_{n2}x_2 + \cdots + a_{nn}x_n=b_n \end{matrix} \right]$

上述記號中矩陣 $A=[ a_{ij}]$ 其中 $a_{ij}$ 是矩陣第 $i$ 行第 $j$ 列的元素。

一般來說，當矩陣 $A$ 是任意的，方程解 $x\mathbf{x}$ 與 $A$ 和 $b\mathbf{b}$ 的數值關系十分復雜，但是 $A$ 是特殊矩陣時，很容易得到解。前面已經介紹了三種特殊矩陣，單位陣、對角陣和正交矩陣都十分容易得到解。

單位陣時，方程 $Ex=bE\mathbf{x}=\mathbf{b}$ ，解就是向量 $b\mathbf{b}$ ，即 $x=b\mathbf{x}=\mathbf{b}$ ，是最簡單的情況。

對角陣時， $D =diag(d_1,d_2,...,d_n)$ ，當對角線元素都不等于0時，對角矩陣可逆，逆為： $D?1=diag(1d1,1d2,...,1dn)D^{-1} =diag(\frac{1}{d_1},\frac{1}{d_2},...,\frac{1}{d_n})$ ，則方程 $Dx=bD\mathbf{x}=\mathbf{b}$ 解為 $x=D?1b\mathbf{x}=D^{-1}\mathbf{b}$ 即 $x_i = b_i/d_i$ 是次最簡單的情況。

正交矩陣時，方程 $Qx=bQ\mathbf{x}=\mathbf{b}$ 解為 $x=QTb\mathbf{x}=Q^T\mathbf{b}$ 是較簡單的情況。

還有一種比較簡單的情況，矩陣為三角矩陣：上三角矩陣 $U$ 和下三角矩陣 $L$ 。

定義上三角矩陣 $U$ 是對角線左下角元素都是 $0$ ，即 $a_{ij}=0, i>j$ ，非零元素都在矩陣右上角。下三角矩陣 $L$ 是對角線右上角元素都是0，即 $a_{ij}=0, i < j$ ，非零元素都在矩陣左下角。
$\left[ \begin{matrix} 1 & 2 & 3\\ 0 & 4 & 5\\ 0 & 0 & 6 \end{matrix} \right]， L = \left[ \begin{matrix} 1 & 0 & 0\\ 2 & 4 & 0\\ 3 & 5 & 6 \end{matrix} \right]， A = \left[ \begin{matrix} 1 & 1 & 0\\ 2 & 4 & 0\\ 3 & 5 & 6 \end{matrix} \right]。$
上面矩陣 $A$ 不是三角陣。注意，上三角矩陣右上角元素也可以為 $0$ 。

三角矩陣能快速求得方程的解，以上三角矩陣為例，此時方程為：
$[a11x1+a12x2+?+a1nxn=b10x1+a22x2+?+a2nxn=b2?0x1+0x2+?+annxn=bn]\left[ \begin{matrix} a_{11}x_1 + a_{12}x_2 + \cdots + a_{1n}x_n=b_1\\ 0x_1 + a_{22}x_2 + \cdots + a_{2n}x_n=b_2\\ \vdots \\ 0x_1 + 0x_2 + \cdots + a_{nn}x_n=b_n \end{matrix} \right]$

根據最后一個方程，可得 $x_n=b_n/a_{nn}$ ，根據倒數第二個方程 $0x1+0x2+?+an?1,n?1xn?1+an?1,nxn=bn?10x_1 + 0x_2 + \cdots + a_{n-1,n-1}x_{n-1} + a_{n-1,n}x_n=b_{n-1}$ 可得 $x_{n-1}=(b_{n-1} - a_{n-1,n}x_n)/a_{n-1,n-1}$ ，同理，按倒推法，可得倒數第 $i + 1$ 個方程的解為
$xn?i=(bn?i?∑j=i+1nan?i,jxj)/an?i,n?ix_{n-i}=(b_{n-i} - \sum^n_{j=i+1} a_{n-i,j}x_j)/a_{n-i,n-i}$

同理，下三角矩陣用同樣方法可求得解，只是不需要倒推，而是順推。由于我們習慣上三角矩陣，所以后面只以上三角矩陣為例。根據公式 $(5)$ ，每個解分量需要除以對角線元素 $a_{ii}$ ，只要其不為零，則方程必有解且唯一。

重要性質 三角矩陣對角線元素 $a_{ii}$ 均不等于 $0$ ，則方程存在解且唯一，三角矩陣可逆。

重要性質 上三角矩陣的逆矩陣是上三角矩陣。

證：令上三角矩陣為 $U=[ a_{ij}],a_{ij}=0, i>j$ ，其逆矩陣為 $U^{-1}=[ b_{ij}]$ 。即 $b_{ij}=0, i>j$ 。根據 $UU^{-1} = E$ ， $e_{ij}], e_{ij}=1,i=j;0,i \ne j$ 根據矩陣乘法第 $4$ 種計算方式，即矩陣 $U$ 第 $i$ 行向量和矩陣 $U^{-1}$ 第 $j$ 列向量的內積為 $(i, j)$ 元素
$ai1b1j+ai2b2j+?+ainbnj=eija_{i1}b_{1j}+a_{i2}b_{2j}+\cdots+a_{in}b_{nj}=e_{ij}$
令 $i = n, j < i$ 根據 $a_{ij}=0, i>j$ ，得 $0b1j+0b2j+?+annbnj=00b_{1j}+0b_{2j}+\cdots+a_{nn}b_{nj}=0$ 得 $b_{nj}=0$ ，即逆矩陣最后一行均為 $0$ ，除了對角線元素外。同理令 $i = n ? 1, j < i$ 根據 $a_{ij}=0, i>j$ ，得 $0b1j+0b2j+?+an?1,n?1bn?1,j+an?1,nbnj=00b_{1j}+0b_{2j}+\cdots+a_{n-1,n-1}b_{n-1,j}+a_{n-1,n}b_{nj}=0$ 得 $b_{n-1,j}=0$ ，即逆矩陣倒數第二行前 $n ? 2$ 個元素均為 $0$ 。依次類推，可得逆矩陣倒數第 $i$ 行前 $n ? i$ 個元素均為 $0$ 。

這個性質很好，上三角矩陣左下角 $0$ 元素，逆矩陣對應位置也是 $0$ 元素， $0$ 元素保持了位置不變。但必須指出，上三角矩陣右上角 $0$ 元素，逆矩陣對應位置一般不是 $0$ 元素。

逆矩陣元素可采用待定系數法快速方便地計算，讀者可自行推導。

總結

以上是生活随笔為你收集整理的4.1 简单方程的解的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。

上一篇： 3.5 矩阵 $4$ 个空间和方程 $A
下一篇： 4.2 矩阵分解概念