6-Advanced 常数变易法求解高阶线性微分方程

6-Advanced 常数变易法求解高阶线性微分方程

什么是常数变易法？

有以下一阶线性微分方程：
$$y^{\prime }+P(x)y=Q(x)\text{(1)}$$
其中，$P(x)>0$ 且 $Q(x)>0$。

解齐次方程

若解其对应的齐次方程：
$$y^{\prime }+P(x)y=0\text{(2)}$$
则易有：
$$y=C{e}^{-\int P(x)dx}(C\ne 0)$$
即为齐次方程的通解.

这时，我们可以用常数变易法来求非齐次方程(1)的通解，即将齐次方程(2)的通解中的常数，变易为一个关于未知函数的函数。变易之后，非齐次方程通解表示如下：
$$y=u(x)e^{-\int Pdx}\text{(3)}$$

解非齐次方程

于是将该通解形式代入原方程(1)，可以解得：
$$u(x)=\int Q(x)e^{\int Pdx}dx+C$$
将上式代入(3)式，即可解得：
$$y=e^{-\int Pdx}\left(\int Q(x)e^{\int Pdx}dx+C\right)$$
这就是所谓常数变易法。

不准确的直观理解

对于常数变易法，以前的理解是：

既然 $y=C{e}^{-\int Pdx}$ 可以使齐次方程 $y+P(x)y=0$ 成立，那么在其基础上增添一个函数，就应该使得该方程运算结果多出一个与自由项相关的余项 $Q(x)$，所以可以使用常数变易法。

这样的理解是基于表面形式做出的一个解释，然而还是不能够明确地说明这个方法的正当性与正确性。所以我们需要进一步探究其内在的原理。

常数变易法的原理

基本思想

容易理解，我们可以将任意函数表示成为两个函数之积，即：
$$y(x)=u(x)\cdot v(x)\text{(4)}$$
对 $y(x)$ 求导，得：
$$y^{\prime }(x)=u^{\prime }(x)v(x)+u(x)v^{\prime }(x)$$

实现过程

将 $y(x)=u(x)v(x)$ ，$y^{\prime }(x)=u^{\prime }(x)v(x)+u(x)v^{\prime }(x)$ 代入非齐次方程(1)，整理得到：
$$u^{\prime }(x)v(x)+u(x)\cdot [v^{\prime }(x)+P(x)v(x)]=Q(x)\text{(5)}$$
由解一阶线性微分方程的常用方法分离变量法容易想到，如果没有 $u(x)\cdot [v^{\prime }(x)+P(x)v(x)]$ 这一项，我们就可以简便地利用分离变量法进行计算。

现在单独考察 $u(x)\cdot [v^{\prime }(x)+P(x)v(x)]$ ，我们可知 $u(x)$ 是不确定的，不能用来确保$u(x)[v^{\prime }(x)+P(x)v(x)]\equiv 0$ 。

考虑另一个因式 $v^{\prime }(x)+P(x)v(x)=0$，显然$v(x)$ 是不确定的，故$v^{\prime }(x)+P(x)v(x)$这一项可以任意取值，假设
$$v^{\prime }(x)+P(x)v(x)\equiv 0\text{(6)}$$
观察式(6)，可以看到其形式与式(2)基本一致。

对式（6）进行求解，可以得到其通解形式：
$$v(x)=C_1\cdot e^{-\int P(x)dx}\text{(7)}$$
将所得通解代入（4），则
$$y(x)=u(x)\cdot C_1\cdot e^{-\int P(x)dx}\text{(8)}$$
将（7）式代入（5）式，得到：
$$u’(x)\cdot C_1\cdot e^{-\int P(x)dx}=Q(x)\text{(9)}$$
使用分离变量法，容易解得：
$$u(x)=\frac{1}{C_1}\int Q(x)e^{\int P(x)dx}dx+C_2\text{(10)}$$
将（7）和（9）同时代入式（4），则
$$y(x)=e^{-\int P(x)dx}\left(\int Q(x)e^{\int P(x)dx}dx+C_1{C}_2\right)\text{(11)}$$
令 $C=C_1{C}_2$, 则得原一阶线性微分方程的通解为：
$$y(x)=e^{-\int P(x)dx}\left(\int Q(x)e^{\int P(x)dx}dx+C\right)\text{(12)}$$

推广到高阶线性微分方程

这是一阶线性微分方程的通解形式，其中：
$$u(x)=\frac{1}{C_1}\int Q(x)e^{\int P(x)dx}dx+C_2\text{(10)}$$

$$v(x)=C_1\cdot e^{-\int P(x)dx}\text{(7)}$$

可以根据具体问题中的 $P(x)$ 和 $Q(x)$ 进行相应的代入，得到具体的解。

多项式的高阶微分形式与展开

现在有一般n阶线性微分方程

给定方程如下：
$$P_0(x)y+P_1(x)y^{\prime }+P_2(x)y^{\prime \prime }+\dots +P_n(x)y^{(n)}=Q(x)(10)$$
由前述有，$y(x)$ 可以表示为 $y(x)=u(x)v(x)$。

两函数乘积的高阶微分形式

现在我们考察两函数乘积的高阶微分形式。比较二项式展开定理，我们不难发现，对 $y=uv$ 的高阶微分具有类似的形式。比如：
$$(uv{)}^{\prime }=v^{\prime }u+u{v}^{\prime }(uv{)}^{\prime \prime }=v^{\prime \prime }v+2u{v}^{\prime \prime }+u{v}^{\prime \prime }$$
从原理上来看，展开多项式的每一项都应有高阶微分，而这些高阶微分分别分配在 $u$、$v$ 上。对于多项式的每一项，相当于任选 $n$ 个微分算子作用于 $u$，则另有 $n$ 个微分算子作用于 $v$，与二项式展开定理本质相同，所以展开形式也应相同。

则有式（11）：
$$(uv{)}^{(n)}=C_n^0{u}^{(n)}v+C_n^1{u}^{(n-1)}{v}^{(1)}+C_n^2{u}^{(n-2)}{v}^{(2)}\dots +C_n^{n-1}{u}^{(1)}{v}^{(n-1)}+C_n^{n}u{v}^{(n)}(11)$$

对原方程进行转化

将这个一般形式代回式(10)，假设将 $u$ 作为主要研究对象（以 $v$ 为主要研究对象亦可，二者地位相同），则按 $u$ 的导数降阶排列多项式：
$$M_{n-1}(x)u^{(n)}+M_{n-2}(x)u^{(n-1)}+\dots +M_0(x)u^{\prime }+(P_0(x)v+P_1(x)v^{\prime }+P_2(x)v^{\prime \prime }+\dots +P_n(x)v^{(n)})u=Q(x)(12)$$
其中，$M_i(x)$ 为关于 $x$ 的多项式。

按一阶情况下的原理，可以令多项式 $P_0(x)v+P_1(x)v^{\prime }+P_2(x)v^{\prime \prime }+\dots +P_n(x)v^{(n)}$ 消去$u$项，解 $v$ 即为解式$10$对应的齐次线性微分方程：$P_0(x)v+P_1(x)v^{\prime }+P_2(x)v^{\prime \prime }+\dots +P_n(x)v^{(n)}=0$。

则剩下的式子为：
$$M_{n-1}(x)u^{(n)}+M_{n-2}(x)u^{(n-1)}+\dots +M_0(x)u^{\prime }=Q(x)$$
令$\alpha (x)=u^{\prime }(x)$，则上式化为：
$$M_{n-1}(x){\alpha }^{(n-1)}+M_{n-2}(x){\alpha }^{(n-2)}+\dots +M_0(x)\alpha =Q(x)$$

结论

比较式$(12)、(13)$，可以看到：通过常数变易法，成功地把求解一个n阶线性微分非齐次方程的问题，变成求解一个对应的$n$阶线性微分齐次方程和一个$(n-1)$阶线性微分非齐次方程的问题。