二阶常系数非齐次线性微分方程的特解公式

本方法只有一个公式，和算子法比减少记忆负担，和待定系数法相比减少计算量。
二阶常系数非齐次线性微分方程：
$$y^{\prime \prime }+p{y}^{\prime }+qy=f(x)(1)$$
其对应的齐次微分方程为：
$$y^{\prime \prime }+p{y}^{\prime }+qy=0(2)$$
设$y_1(x),y_2(x)$是方程(2)有两个线性无关的特解，则方程(1)有特解:
$$y^{\ast }=y_1\int \frac{f(x)}{y_2{\left(\frac{y_1}{y_2}\right)}^{\prime }}\mathrm{d}x+y_2\int \frac{f(x)}{y_1{\left(\frac{y_2}{y_1}\right)}^{\prime }}\mathrm{d}x$$
与一阶线性微分方程求解公式类似，这里不定积分不用+C

例题1，$y^{\prime \prime }+2y^{\prime }+y=(3x+2)e^{-x}$

首先解出两个齐次特解：
$$r^2+2r+1=0\Rightarrow r_{1,2}=-1\Rightarrow y_1=x{e}^{-x},y_2=e^{-x}$$
再代公式解出非齐次特解：
$$\begin{array}{rl}& y^{\ast }=x{e}^{-x}\int \frac{(3x+2)e^{-x}}{e^{-x}(x{)}^{\prime }}\mathrm{d}x+e^{-x}\int \frac{(3x+2)e^{-x}}{x{e}^{-x}{\left(\frac{1}{x}\right)}^{\prime }}\mathrm{d}x\\ & =x{e}^{-x}\int (3x+2)\mathrm{d}x-e^{-x}\int \left(3x^2+2x\right)\mathrm{d}x\\ & =x{e}^{-x}\left(\frac{3}{2}{x}^2+2x\right)-e^{-x}\left(x^3+x^2\right)\\ & =\frac{1}{2}{x}^3{e}^{-x}+x^2{e}^{-x}\end{array}$$

例题2，[2007数1,2]求$y^{\prime \prime }-4y^{\prime }+3y=2{\mathrm{e}}^{2x}$的通解

首先解出两个齐次特解：
$$r^2-4r+3=0⟹r_1=1,r_2=3⟹y_1={\mathrm{e}}^x,y_2={\mathrm{e}}^{3x}$$
再代公式解出非齐次特解：
$$\begin{array}{rl}{y}^{\ast }& ={\mathrm{e}}^x\int \frac{2{\mathrm{e}}^{2x}}{{\mathrm{e}}^{3x}{\left(\frac{{\mathrm{e}}^x}{{\mathrm{e}}^{3x}}\right)}^{\prime }}\mathrm{d}x+{\mathrm{e}}^{3x}\int \frac{2{\mathrm{e}}^{2x}}{{\mathrm{e}}^x{\left(\frac{{\mathrm{e}}^{3x}}{{\mathrm{e}}^x}\right)}^{\prime }}\mathrm{d}x\\ & ={\mathrm{e}}^x\int \frac{2{\mathrm{e}}^{2x}}{-2{\mathrm{e}}^{-2x}\cdot {\mathrm{e}}^{3x}}\mathrm{d}x+{\mathrm{e}}^{3x}\int \frac{2{\mathrm{e}}^{2x}}{2{\mathrm{e}}^x\cdot {\mathrm{e}}^{2x}}\mathrm{d}x\\ & =-{\mathrm{e}}^x\int {\mathrm{e}}^x\mathrm{d}x+{\mathrm{e}}^{3x}\int {\mathrm{e}}^{-x}\mathrm{d}x\\ & =-2{\mathrm{e}}^{2x}\end{array}$$
因此该方程的通解为：
$$y=c_1{\mathrm{e}}^{3x}+c_2{\mathrm{e}}^x-2{\mathrm{e}}^{2x}$$

例题3，[2010数1]求$y^{\prime \prime }-3y^{\prime }+2y=2x{\mathrm{e}}^x$的通解

首先解出两个齐次特解：
$$r^2-3r+2=0⟹r_1=1,r_2=2⟹y_1={\mathrm{e}}^x,y_2={\mathrm{e}}^{2x}$$

再代公式解出非齐次特解：
$$\begin{array}{rl}{y}^{\ast }& ={\mathrm{e}}^x\int \frac{2x{\mathrm{e}}^x}{{\mathrm{e}}^{2x}{\left(\frac{{\mathrm{e}}^x}{{\mathrm{e}}^{2x}}\right)}^{\prime }}\mathrm{d}x+{\mathrm{e}}^{2x}\int \frac{2x{\mathrm{e}}^x}{{\mathrm{e}}^x{\left(\frac{{\mathrm{e}}^{2x}}{{\mathrm{e}}^x}\right)}^{\prime }}\mathrm{d}x\\ & ={\mathrm{e}}^x\int \frac{2x{\mathrm{e}}^x}{-{\mathrm{e}}^x}\mathrm{d}x+{\mathrm{e}}^{2x}\int \frac{2x{\mathrm{e}}^x}{{\mathrm{e}}^{2x}}\mathrm{d}x\\ & =-{\mathrm{e}}^x\int 2x\mathrm{d}x+2{\mathrm{e}}^{2x}\int x{\mathrm{e}}^{-x}\mathrm{d}x\\ & =-x^2{\mathrm{e}}^x-2{\mathrm{e}}^{2x}\left(x{\mathrm{e}}^{-x}+{\mathrm{e}}^{-x}\right)\\ & =-x(x+2){\mathrm{e}}^x+2{\mathrm{e}}^x\end{array}$$
因此该方程的通解为：
$$y=c_1{\mathrm{e}}^x+c_2{\mathrm{e}}^{2x}-x(x+2){\mathrm{e}}^x$$

例题4，[同济第七版高数例题]求$y^{\prime \prime }+y=x\cos 2x$的特解

首先解出两个齐次特解：
$$r^2+1=0⟹r_{1,2}=\pm i⟹y_1=\cos x,y_2=\sin x$$
再代公式解出非齐次特解：
$$\begin{array}{rl}{y}^{\ast }& =\cos x\int \frac{x\cos 2x}{\sin x\cdot {\left(\frac{\cos x}{\sin x}\right)}^{\prime }}\mathrm{d}x+\sin x\int \frac{x\cos 2x}{\cos x\cdot {\left(\frac{\sin x}{\cos x}\right)}^{\prime }}\mathrm{d}x\\ & =-\cos x\int x\sin x\cos 2x\mathrm{d}x+\sin x\int x\cos x\cos 2x\mathrm{d}x\\ & =-\frac{1}{2}\cos x\int x(\sin 3x-\sin x)\mathrm{d}x+\frac{1}{2}\sin x\int x(\cos 3x+\cos x)\mathrm{d}x\\ & =-\frac{1}{2}\cos x\cdot \left(-\frac{1}{3}x\cos 3x+\frac{1}{9}\sin 3x+x\cos x-\sin x\right)+\frac{1}{2}\sin x\cdot \left(\frac{1}{3}x\sin 3x+\frac{1}{9}\cos 3x+x\sin x+\cos x\right)\\ & =\frac{1}{6}x\cos x\cos 3x-\frac{1}{18}\cos x\sin 3x-\frac{1}{2}x{\cos }^{2}x+\frac{1}{2}\sin x\cos x+\frac{1}{6}x\sin x\sin 3x+\frac{1}{18}\sin x\cos 3x+\frac{1}{2}x{\sin }^{2}x+\frac{1}{2}\sin x\cos x\\ & =\frac{1}{6}x\cos 2x-\frac{1}{18}\sin 2x-\frac{1}{2}x\cos 2x+\frac{1}{2}\sin 2x\\ & =-\frac{1}{3}x\cos 2x+\frac{4}{9}\sin 2x\end{array}$$
看起来计算量很大，不过很多同类项可以快速合并，所以实际上仍然要比待定系数法计算量小，我们可以对比一下待定系数法的计算量。注意到$2i$不是特征方程的根，所以可设特解为：
$$y^{\ast }=(ax+b)\cos 2x+(cx+d)\sin 2x$$
代入方程得：
$$(-3ax-3b+4c)\cos 2x-(3cx+3d+4a)\sin 2x=x\cos 2x\text{. }$$
比较系数得：
$$\{\begin{array}{l}-3a=1\\ -3b+4c=0\\ -3c=0\\ -3d-4a=0\end{array}$$
由此解得：
$$a=-\frac{1}{3},b=0,c=0,d=\frac{4}{9}\text{. }$$
于是特解为：
$$y^{\ast }=-\frac{1}{3}x\cos 2x+\frac{4}{9}\sin 2x\text{. }$$
其中代入方程这里要求一个二阶导，计算量非常大

例题5，[同济第七版高数例题]求$y^{\prime \prime }-y={\mathrm{e}}^x\cos 2x$的特解

首先解出两个齐次特解：
$$r^2-1=0⟹r_{1,2}=\pm 1⟹y_1={\mathrm{e}}^x,y_2={\mathrm{e}}^{-x}$$
再代公式解出非齐次特解：
$$\begin{array}{rl}{y}^{\ast }& ={\mathrm{e}}^x\int \frac{{\mathrm{e}}^x\cos 2x}{{\mathrm{e}}^{-x}\cdot {\left(\frac{{\mathrm{e}}^x}{{\mathrm{e}}^{-x}}\right)}^{\prime }}\mathrm{d}x+{\mathrm{e}}^{-x}\int \frac{{\mathrm{e}}^x\cos 2x}{{\mathrm{e}}^x\cdot {\left(\frac{{\mathrm{e}}^{-x}}{{\mathrm{e}}^x}\right)}^{\prime }}\mathrm{d}x\\ & ={\mathrm{e}}^x\int \frac{{\mathrm{e}}^x\cos 2x}{{\mathrm{e}}^{-x}\cdot 2{\mathrm{e}}^{2x}}\mathrm{d}x+{\mathrm{e}}^{-x}\int \frac{{\mathrm{e}}^x\cos 2x}{-2{\mathrm{e}}^x{\mathrm{e}}^{-2x}}\mathrm{d}x\\ & =\frac{1}{2}{\mathrm{e}}^x\int \cos 2x\mathrm{d}x-\frac{1}{2}{\mathrm{e}}^{-x}\int {\mathrm{e}}^{2x}\cos 2x\mathrm{d}x\\ & =\frac{1}{4}{\mathrm{e}}^x\sin 2x-\frac{1}{4}{\mathrm{e}}^{-x}\int {\mathrm{e}}^{2x}\mathrm{d}(\sin 2x)\\ & =\frac{1}{2}{\mathrm{e}}^{-x}\int \sin 2x\cdot {\mathrm{e}}^{2x}\mathrm{d}x\\ & =\frac{1}{8}{\mathrm{e}}^x(\sin 2x-\cos 2x)\end{array}$$
其中最后一个等号可以套公式，具体公式及其推导如下：
$$\begin{array}{rl}& I_1=\int {\mathrm{e}}^{ax}\sin bx\mathrm{d}x,I_2=\int {\mathrm{e}}^{ax}\cos bx\mathrm{d}x\\ & I_1=-\frac{1}{b}\int {\mathrm{e}}^{ax}\mathrm{d}\cos bx=-\frac{1}{b}\left({\mathrm{e}}^{ax}\cos bx-a\int {\mathrm{e}}^{ax}\cos bx\mathrm{d}x\right)=\frac{a}{b}{I}_2-\frac{1}{b}{\mathrm{e}}^{ax}\cos bx\\ & I_2=\frac{1}{b}\int {\mathrm{e}}^{ax}\mathrm{d}\sin bx=\frac{1}{b}\left({\mathrm{e}}^{ax}\sin bx-a\int {\mathrm{e}}^{ax}\sin bx\mathrm{d}x\right)=\frac{1}{b}{\mathrm{e}}^{ax}\sin bx-\frac{a}{b}{I}_1\\ & \{\begin{array}{l}b{I}_1-a{I}_2=-{\mathrm{e}}^{ax}\cos bx\\ a{I}_2+b{I}_1={\mathrm{e}}^{ax}\sin bx\end{array}⟹\{\begin{array}{l}{I}_1=\frac{{\mathrm{e}}^{ax}}{a^2+b^2}(a\sin bx-b\cos bx)+C\\ I_2=\frac{{\mathrm{e}}^{ax}}{a^2+b^2}(b\sin bx+a\cos bx)+C\end{array}\end{array}$$
看起来计算量很大，但实际上仍然要比待定系数法计算量小，我们可以对比一下待定系数法的计算量。注意到$1+2i$不是特征方程的根，所以可设特解为：
$$y^{\ast }={\mathrm{e}}^x(a\cos 2x+b\sin 2x)$$
求一阶导和二阶导：
$$\begin{array}{rl}& y^{\ast \prime }={\mathrm{e}}^x[(a+2b)\cos 2x+(-2a+b)\sin 2x],\\ & y^{\ast \prime \prime }={\mathrm{e}}^x[(-3a+4b)\cos 2x+(-4a-3b)\sin 2x].\end{array}$$
代入方程得：
$$4{\mathrm{e}}^x[(-a+b)\cos 2x-(a+b)\sin 2x]={\mathrm{e}}^x\cos 2x,$$
比较系数得：
$$\{\begin{array}{l}-a+b=\frac{1}{4},\\ a+b=0,\end{array}⟹\{\begin{array}{l}a=-\frac{1}{8},\\ b=\frac{1}{8}.\end{array}$$

于是特解为：
$$y^{\ast }=\frac{1}{8}{\mathrm{e}}^x(\sin 2x-\cos 2x)$$
其中代入方程这里的计算量非常大

参考文献：

[1]数学群管理员魏念辉给出此公式(与文献3本质上相同，但形式上更简单)
[2]魏佬给出了云哥的证明过程(与文献3基本一致)
[3]朱德刚. (2010). 二阶常系数非齐次线性微分方程的特解公式. 高等数学研究, 13(3), 2.