利用特征值与特征向量求解耦合微分方程组的疑问

我来帮你一步步梳理这个问题，看看哪里可能出了偏差。首先，我们先明确你给出的方程组：

$$\frac{d}{dt}\begin{bmatrix}
y \
z
\end{bmatrix} = \begin{bmatrix} -2g & -\Delta \ \Delta & 0 \end{bmatrix}\begin{bmatrix}y \ z\end{bmatrix}, \quad \text{初始条件 } \begin{bmatrix}y(0) \ z(0)\end{bmatrix} = \begin{bmatrix}0 \ 1\end{bmatrix}$$

教程给出的解是 $z(t) = \cos(2\pi t)e^{-gt}$，但你算不出来，那我们先从特征值计算开始，一步步推导正确的解，再对比教程的结果找差异。

第一步：计算系统矩阵的特征值 #

设系统矩阵为 $A = \begin{bmatrix} -2g & -\Delta \ \Delta & 0 \end{bmatrix}$，特征值 $\lambda$ 满足特征方程 $\det(A - \lambda I) = 0$：
$$
\det\begin{bmatrix} -2g - \lambda & -\Delta \ \Delta & -\lambda \end{bmatrix} = 0
$$
展开行列式：
$$
(-2g - \lambda)(-\lambda) - (-\Delta)(\Delta) = \lambda^2 + 2g\lambda + \Delta^2 = 0
$$
解这个二次方程，得到特征值：
$$
\lambda = -g \pm i\sqrt{\Delta^2 - g^2} \quad (\text{假设 } \Delta > g \text{，此时特征值为复数，对应振荡解})
$$
记 $\omega = \sqrt{\Delta^2 - g^2}$，则特征值可写成 $\lambda_1 = -g + i\omega$，$\lambda_2 = -g - i\omega$。

第二步：求对应的特征向量 #

对 $\lambda_1 = -g + i\omega$，代入 $(A - \lambda_1 I)\vec{v}1 = 0$：
$$
\begin{bmatrix} -g - i\omega & -\Delta \ \Delta & g - i\omega \end{bmatrix}\begin{bmatrix}v{11} \ v_{12}\end{bmatrix} = 0
$$
由第一行方程可得：$(-g - i\omega)v_{11} = \Delta v_{12}$，取 $v_{11} = \Delta$，则 $v_{12} = -g - i\omega$，即特征向量 $\vec{v}_1 = \begin{bmatrix}\Delta \ -g - i\omega\end{bmatrix}$。

同理，$\lambda_2$ 对应的特征向量 $\vec{v}_2 = \begin{bmatrix}\Delta \ -g + i\omega\end{bmatrix}$。

第三步：构造通解并代入初始条件 #

线性齐次方程组的通解为：
$$
\begin{bmatrix}y(t) \ z(t)\end{bmatrix} = C_1 e^{\lambda_1 t}\vec{v}_1 + C_2 e^{\lambda_2 t}\vec{v}_2
$$
代入初始条件 $t=0$ 时 $\begin{bmatrix}0 \ 1\end{bmatrix}$：

1. $C_1\Delta + C_2\Delta = 0 \implies C_2 = -C_1$
2. $C_1(-g - i\omega) + C_2(-g + i\omega) = 1$

将 $C_2 = -C_1$ 代入第二个方程：
$$
C_1(-g - i\omega) - C_1(-g + i\omega) = 1 \implies -2i\omega C_1 = 1 \implies C_1 = \frac{i}{2\omega}, \quad C_2 = -\frac{i}{2\omega}
$$

第四步：推导 $z(t)$ 的表达式 #

提取通解的第二个分量（即 $z(t)$）：
$$
z(t) = C_1(-g - i\omega)e^{(-g+i\omega)t} + C_2(-g + i\omega)e^{(-g-i\omega)t}
$$
代入 $C_1$ 和 $C_2$，并利用欧拉公式 $e^{i\theta} = \cos\theta + i\sin\theta$ 化简，最终会得到：
$$
z(t) = e^{-gt}\left( \cos(\omega t) + \frac{g}{\omega}\sin(\omega t) \right)
$$
其中 $\omega = \sqrt{\Delta^2 - g^2}$。

和教程结果的差异分析 #

教程给出的解是 $z(t) = \cos(2\pi t)e^{-gt}$，这意味着我们推导结果中的 $\frac{g}{\omega}\sin(\omega t)$ 项消失了，且 $\omega = 2\pi$。要满足这个情况，有两种可能：

1. 矩阵参数或形式有误：如果系统矩阵是 $\begin{bmatrix} -g & -\Delta \ \Delta & -g \end{bmatrix}$（注意右下角是 $-g$ 而非 $0$），那么特征值会是 $\lambda = -g \pm i\Delta$，代入初始条件后，$z(t)$ 会直接化简为 $\cos(\Delta t)e^{-gt}$，当 $\Delta=2\pi$ 时就和教程结果完全一致。
2. 初始条件或参数设定不同：如果 $g=0$，那么 $\omega=\Delta$，此时 $z(t)=\cos(\Delta t)$，但会缺少 $e^{-gt}$ 项；如果教程的初始条件不是 $\begin{bmatrix}0 \1\end{bmatrix}$，也可能得到纯余弦形式的解，但结合 $e^{-gt}$ 项来看，第一种矩阵形式错误的可能性更大。

建议的排查方向 #

• 核对教程中的系统矩阵是否和你写的一致，大概率是矩阵元素的抄写错误；
• 检查自己计算特征值、特征向量时的符号是否出错，尤其是行列式展开和特征向量求解步骤；
• 可以尝试把教程的解代入你给出的方程组，验证是否满足，这样能快速判断解和方程组是否匹配。

备注：内容来源于stack exchange，提问作者LieAlgebraGuy1999