非齐次微分方程转化为线性系统及含参数a的微分方程流求解咨询
没问题,我来帮你理清楚这个非齐次二阶微分方程怎么转化为一阶线性系统,再一步步分析不同a值下的解(也就是你说的“流”)。
首先,你熟悉齐次方程$x'' + ax' + bx=0$的转化方法,那带常数项的非齐次方程其实只是多了一步处理非齐次项的操作。我们依然用状态变量替换来降阶:
- 设第一个状态变量 $x_1 = x$(就是原函数)
- 设第二个状态变量 $x_2 = x'$(原函数的一阶导数)
接下来写出每个状态变量的导数:
- $x_1' = x' = x_2$,这个和齐次情况完全一样
- $x_2' = x''$,而原方程变形后是$x'' = -ax' -1$,代入$x'=x_2$,就得到 $x_2' = -a x_2 -1$
这样就得到了一个非齐次一阶线性系统,写成矩阵形式就是:
$$
\begin{pmatrix} x_1' \ x_2' \end{pmatrix} = \begin{pmatrix} 0 & 1 \ 0 & -a \end{pmatrix} \begin{pmatrix} x_1 \ x_2 \end{pmatrix} + \begin{pmatrix} 0 \ -1 \end{pmatrix}
$$
这里的矩阵$A = \begin{pmatrix}0&1\0&-a\end{pmatrix}$就是系统的系数矩阵,后面的$\begin{pmatrix}0\-1\end{pmatrix}$就是非齐次项的常数向量。
接下来我们可以通过分析系数矩阵的特征值,结合非齐次系统的解法来求通解,也就是你要的“流”:
分情况讨论参数a的取值
情况1:$a \neq 0$
先求系数矩阵$A$的特征值,解特征方程$\det(A - \lambda I) = 0$:
$$
\det\begin{pmatrix} -\lambda & 1 \ 0 & -a - \lambda \end{pmatrix} = \lambda(\lambda + a) = 0
$$
得到两个不同的特征值:$\lambda_1=0$,$\lambda_2=-a$,矩阵可对角化。
- 齐次系统的通解:对应两个特征值的特征向量分别是$\begin{pmatrix}1\0\end{pmatrix}$(对应$\lambda=0$)和$\begin{pmatrix}1\-a\end{pmatrix}$(对应$\lambda=-a$),所以齐次解为:
$$
\mathbf{y}_h(t) = C_1 \begin{pmatrix}1\0\end{pmatrix} + C_2 \begin{pmatrix}1\-a\end{pmatrix} e^{-a t}
$$ - 非齐次系统的特解:因为非齐次项是常数,而$\lambda=0$是特征值,常数特解不存在,所以假设特解为$\mathbf{y}_p(t) = \begin{pmatrix} -t/a \ -1/a \end{pmatrix}$(代入验证满足系统方程)。
把齐次解和特解合并,再转化回原变量$x(t)=x_1(t)$,就得到原方程的通解:
$$
x(t) = C_1 + C_2 e^{-a t} - \frac{t}{a}
$$
情况2:$a = 0$
此时原方程简化为$x'' + 1 = 0$,对应的线性系统是:
$$
\begin{pmatrix} x_1' \ x_2' \end{pmatrix} = \begin{pmatrix} 0 & 1 \ 0 & 0 \end{pmatrix} \begin{pmatrix} x_1 \ x_2 \end{pmatrix} + \begin{pmatrix} 0 \ -1 \end{pmatrix}
$$
系数矩阵$A$是Jordan块,特征值为$\lambda=0$(二重)。
- 齐次系统的通解:对应Jordan块的解为$\mathbf{y}_h(t) = C_1 \begin{pmatrix}1\0\end{pmatrix} + C_2 \begin{pmatrix}t\1\end{pmatrix}$
- 非齐次系统的特解:假设特解为$\mathbf{y}_p(t) = \begin{pmatrix} -t^2/2 \ -t \end{pmatrix}$,代入验证满足系统方程。
合并后转化回原变量,得到通解:
$$
x(t) = C_1 + C_2 t - \frac{t^2}{2}
$$
这样不管a取什么值,我们都得到了方程的通解,也就是对应不同初始条件下的“流”啦。
备注:内容来源于stack exchange,提问作者Kajice

