耦合谐振子二阶微分方程转化为一阶矩阵形式的相关技术咨询

我来帮你理清楚耦合谐振子从二阶微分方程组转化为一阶矩阵形式的整个过程哈，步骤很清晰，咱们一步步来：

首先，先看你给出的原始二阶耦合微分方程组：
$$
\left{\begin{array}{rcrcr}
x_1'' & = & -(k_1 + k_2),x_1 & + & k_2,x_2
\[1mm]
x_2'' & = & k_2,x_1 & - & (k_1 + k_2),x_2
\end{array}\right.
$$

要把二阶方程转化为一阶，核心方法就是引入新的状态变量降阶：
我们令：

• $y_1 = x_1'$（即$x_1$的一阶导数）
• $y_2 = x_2'$（即$x_2$的一阶导数）

这样一来，原来的二阶方程组就可以拆成四个一阶微分方程：
$$
\begin{aligned}
x_1' & = y_1, \
y_1' & = -(k_1 + k_2)x_1 + k_2x_2, \
x_2' & = y_2, \
y_2' & = k_2x_1 - (k_1 + k_2)x_2.
\end{aligned}
$$

接下来就是把这组一阶方程写成矩阵形式 $X' = AX$，首先要定义状态向量 $X$：
我们取状态向量为 $X = \begin{bmatrix} x_1 \ y_1 \ x_2 \ y_2 \end{bmatrix}$，那么它的一阶导数 $X' = \begin{bmatrix} x_1' \ y_1' \ x_2' \ y_2' \end{bmatrix}$。

把四个一阶方程对应到矩阵乘法的形式，就能得到系数矩阵 $A$：
$$A =
\begin{pmatrix}
0 & 1 & 0 & 0 \
-(k_1 + k_2) & 0 & k_2 & 0 \
0 & 0 & 0 & 1 \
k_2 & 0 & -(k_1 + k_2) & 0
\end{pmatrix}$$

咱们可以简单验证一下：把 $A$ 和 $X$ 做矩阵乘法，得到的结果正好是 $\begin{bmatrix} x_1' \ y_1' \ x_2' \ y_2' \end{bmatrix} = X'$，完全匹配咱们的一阶方程组。

你提到的“当计算th...”看起来内容没写完，如果是要计算矩阵的特征值、特征向量或者后续的解的形式，也可以补充完整问题，咱们再继续讨论哦！

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