2021-05-07 10:20:42

This commit is contained in:
wizardforcel
2021-05-07 10:20:42 +08:00
parent c411fa2a99
commit 957b754226
52 changed files with 2953 additions and 22 deletions

151
docs/linalg/chapter07.md Normal file
View File

@@ -0,0 +1,151 @@
# 第七讲:求解$Ax=0$,主变量,特解
举例:$3 \times 4$矩阵
$
A=
\begin{bmatrix}
1 & 2 & 2 & 2\\
2 & 4 & 6 & 8\\
3 & 6 & 8 & 10\\
\end{bmatrix}
$,求$Ax=0$的特解:
找出主变量pivot variable
$$
A=
\begin{bmatrix}
1 & 2 & 2 & 2\\
2 & 4 & 6 & 8\\
3 & 6 & 8 & 10\\
\end{bmatrix}
\underrightarrow{消元}
\begin{bmatrix}
\underline{1} & 2 & 2 & 2\\
0 & 0 & \underline{2} & 4\\
0 & 0 & 0 & 0\\
\end{bmatrix}
=U
$$
主变量pivot variable下划线元素的个数为2即矩阵$A$的秩rank为2即$r=2$。
主变量所在的列为主列pivot column其余列为自由列free column
自由列中的变量为自由变量free variable自由变量的个数为$n-r=4-2=2$。
通常,给自由列变量赋值,去求主列变量的值。如,令$x_2=1, x_4=0$求得特解
$x=c_1\begin{bmatrix}-2\\1\\0\\0\\\end{bmatrix}$
再令$x_2=0, x_4=1$求得特解
$x=c_2\begin{bmatrix}2\\0\\-2\\1\\\end{bmatrix}$。
该例还能进一步简化,即将$U$矩阵化简为$R$矩阵Reduced row echelon form即简化行阶梯形式。
在简化行阶梯形式中,主元上下的元素都是$0$
$$
U=
\begin{bmatrix}
\underline{1} & 2 & 2 & 2\\
0 & 0 & \underline{2} & 4\\
0 & 0 & 0 & 0\\
\end{bmatrix}
\underrightarrow{化简}
\begin{bmatrix}
\underline{1} & 2 & 0 & -2\\
0 & 0 & \underline{1} & 2\\
0 & 0 & 0 & 0\\
\end{bmatrix}
=R
$$
将$R$矩阵中的主变量放在一起,自由变量放在一起(列交换),得到
$$
R=
\begin{bmatrix}
\underline{1} & 2 & 0 & -2\\
0 & 0 & \underline{1} & 2\\
0 & 0 & 0 & 0\\
\end{bmatrix}
\underrightarrow{列交换}
\left[
\begin{array}{c c | c c}
1 & 0 & 2 & -2\\
0 & 1 & 0 & 2\\
\hline
0 & 0 & 0 & 0\\
\end{array}
\right]
=
\begin{bmatrix}
I & F \\
0 & 0 \\
\end{bmatrix}
\textrm{,其中}I\textrm{为单位矩阵,}F\textrm{为自由变量组成的矩阵}
$$
计算零空间矩阵$N$nullspace matrix其列为特解有$RN=0$。
$$
x_{pivot}=-Fx_{free} \\
\begin{bmatrix}
I & F \\
\end{bmatrix}
\begin{bmatrix}
x_{pivot} \\
x_{free} \\
\end{bmatrix}=0 \\
N=\begin{bmatrix}
-F \\
I \\
\end{bmatrix}
$$
在本例中
$
N=
\begin{bmatrix}
-2 & 2 \\
0 & -2 \\
1 & 0 \\
0 & 1 \\
\end{bmatrix}
$,与上面求得的两个$x$特解一致。
另一个例子,矩阵
$
A=
\begin{bmatrix}
1 & 2 & 3 \\
2 & 4 & 6 \\
2 & 6 & 8 \\
2 & 8 & 10 \\
\end{bmatrix}
\underrightarrow{消元}
\begin{bmatrix}
1 & 2 & 3 \\
0 & 2 & 2 \\
0 & 0 & 0 \\
0 & 0 & 0 \\
\end{bmatrix}
\underrightarrow{化简}
\begin{bmatrix}
1 & 0 & 1 \\
0 & 1 & 1 \\
0 & 0 & 0 \\
0 & 0 & 0 \\
\end{bmatrix}
=R
$
矩阵的秩仍为$r=2$,有$2$个主变量,$1$个自由变量。
同上一例,取自由变量为$x_3=1$,求得特解
$
x=c
\begin{bmatrix}
-1 \\
-1 \\
1 \\
\end{bmatrix}
$