14 线性方程组之齐次与非齐次方程组

查看大图

齐次方程组一定有零解，非齐次方程组要先判断是否有解。有解时，通解常写成特解加齐次解。

查看大图

我会先看右侧是不是零向量，再用秩判断解空间。不要把两类方程组混着套结论。

在上一篇中，我们学习了线性方程组的高斯消元法，这是一种有效的求解线性方程组的方法。本篇将讨论线性方程组的两种基本类型：齐次方程组和非齐次方程组。理解这两种方程组的结构和解法对于深入学习线性代数及其在人工智能中的应用至关重要。

1. 齐次方程组

定义

一个线性方程组如果右边的常数项全为零，则称为齐次方程组。它可以表示为：

查看大图

比较齐次与非齐次方程组时，先看常数项、零解、自由变量、秩和通解结构。

$$A\mathbf{x} = \mathbf{0}$$

其中，$A$ 是系数矩阵，$\mathbf{x}$ 是未知变量向量，$\mathbf{0}$ 是零向量。

例子

考虑以下齐次方程组：

$$\begin{align*} 2x + 3y - z &= 0 \\ 4x + 6y - 2z &= 0 \\ -x + y + 2z &= 0 \end{align*}$$

这个方程组可以用矩阵形式表示为 $A\mathbf{x} = \mathbf{0}$，其中：

$$A = \begin{pmatrix} 2 & 3 & -1 \\ 4 & 6 & -2 \\ -1 & 1 & 2 \end{pmatrix}, \quad \mathbf{x} = \begin{pmatrix} x \\ y \\ z \end{pmatrix}, \quad \mathbf{0} = \begin{pmatrix} 0 \\ 0 \\ 0 \end{pmatrix}$$

解的性质

齐次方程组的解总是包括零解，即 $\mathbf{x} = \mathbf{0}$。如果该方程组有非零解，则称为非平凡解。通过高斯消元法，可以求得此方程组的解的形式。当自由变量的数量大于0时，方程组将有无穷多解。

Python 示例

以下是使用 NumPy 求解齐次方程组的示例代码：

import numpy as np

# 定义系数矩阵 A
A = np.array([[2, 3, -1],
              [4, 6, -2],
              [-1, 1, 2]])

# 使用 NumPy 的线性代数库求解 Ax = 0
# 因为是齐次方程组，我们只需要求出零解与非平凡解
u, s, vt = np.linalg.svd(A)  # 奇异值分解
    
# 非零奇异值对应的列就是解的基础
rank = np.sum(s > 1e-10)  # 发现有几个非零的奇异值
null_space_dim = A.shape[1] - rank  # 求解空域的维度

print(f"齐次方程组的解的维度为 {null_space_dim}")

2. 非齐次方程组

定义

如果一个线性方程组的常数项不全为零，则称为非齐次方程组。其形式为：

查看大图

《线性方程组之齐次与非齐次方程组》可以按“场景、概念、动作、结果”来读。先把这四件事对齐，再回到正文里的参数、代码或流程。

$$A\mathbf{x} = \mathbf{b}$$

其中，$\mathbf{b}$ 是非零向量。

例子

考虑以下非齐次方程组：

$$\begin{align*} 2x + 3y - z &= 5 \\ 4x + 6y - 2z &= 10 \\ -x + y + 2z &= 3 \end{align*}$$

用矩阵表示为：

$$A = \begin{pmatrix} 2 & 3 & -1 \\ 4 & 6 & -2 \\ -1 & 1 & 2 \end{pmatrix}, \quad \mathbf{x} = \begin{pmatrix} x \\ y \\ z \end{pmatrix}, \quad \mathbf{b} = \begin{pmatrix} 5 \\ 10 \\ 3 \end{pmatrix}$$

解的性质

非齐次方程组可能具有唯一解、无解或无穷多解。根据方程的阶数和自由变量的数量，我们可以使用高斯消元法或者其他数值方法求解。

Python 示例

以下是使用 NumPy 求解非齐次方程组的示例代码：

import numpy as np

# 定义系数矩阵 A 和常数项向量 b
A = np.array([[2, 3, -1],
              [4, 6, -2],
              [-1, 1, 2]])
b = np.array([5, 10, 3])

# 使用 NumPy 的线性代数库求解 Ax = b
x = np.linalg.solve(A, b)

print(f"非齐次方程组的解为: x = {x}")

查看大图

学完《线性方程组之齐次与非齐次方程组》后，不妨换一个自己的场景试一次，重点观察输入、处理和输出是否能对应起来。

查看大图

如果想把《线性方程组之齐次与非齐次方程组》用到自己的任务里，可以先缩小场景，只验证一个最关键的判断点。

结语

在本篇中，我们介绍了齐次方程组和非齐次方程组的定义、解的性质，以及如何利用Python实现求解。下篇我们将继续讨论特征值与特征向量，欢迎继续关注！