微分方程:无法凑出的结果

怎么导都不变

我们学习到过一个神奇的函数:y=exy=e^x:这个函数不仅处处可导,而且导完还是自己,就像这样:

f(x)=f(x)f\left(x\right)=f^\prime\left(x\right)

以上就是一个微分方程,代表着在方程中表达了函数与其导数之间的关系。由于f(x)=exf\left(x\right)=e^x是以上方程的解,所以以后的微分方程中,我们的结果大多都是u(x) eg(x)u\left(x\right)\cdot\ e^{g\left(x\right)}这个形式。

我先来,这个我会

在之前的学习中,我们已经掌握了一种微分方程的求解办法:我们把dxdxdydy放在等号两边,让dxdx那边只剩xx的项,另一边只有yy项,两边求导就完成了。

它们非常简单,比如下面这个等式:

dydx=2xy\frac{dy}{dx}=2xy

我们只需要将其变成:

1y dy=2x dx\frac{1}{y}\ dy=2x\ dx

两边反求导:

lny+C1=x2+C2\ln{y}+C_1=x^2+C_2

结果就是:

y=ex2+Cy=e^{x^2+C}

请注意,这不是正确结果!我们希望最终结果的CC并不出现在次幂上。因为eCe^C也是一个常数,所以我们真正的标准结果是:

y=Cex2y=Ce^{x^2}

齐次微分方程

热身结束!接下来我们面对的微分方程,将不能使用上面的方法获得答案。或许凑一凑其实还是能凑出来?

面对除开yy^\prime每项xxyy的次数和都相等的方程,这类方程的求解是:将每项都变成只包含(ynxm)q(\frac{y^n}{x^m})^q的形式。说起来可能难以理解,我们不如用个例子来展示一下。比如以下这个微分方程:

y2xyy+x2y=0y^2-xy\cdot y^\prime+x^2\cdot y^\prime=0

这里,我们将每项除以x2x^2,变成:

y2x2yxy+y=0\frac{y^2}{x^2}-\frac{y}{x}\cdot y^\prime+y^\prime=0

此时,我们就组成了只包含yx\frac{y}{x}的形式。那么这里,我们需要令u=yxu=\frac{y}{x}来简化:

u2uy+y=0u^2-uy^\prime+y^\prime=0

化简一下,就是:

y=u2u1y^\prime=\frac{u^2}{u-1}

唔,仿佛到了这里我们还是无从下手呢。原因很简单,因为我们并没有yyuu的任何关系。所以这时,我们需要使用y=uxy=ux来替换yy

y=(ux)=ux+ux=xdudx+uy^\prime=(ux)^\prime=u^\prime x+ux^\prime=x\cdot \frac{du}{dx}+u

带入之后,我们的等式就成为了关于uuxx的微分方程:

u+xdudx=u2u1u+x\cdot \frac{du}{dx}=\frac{u^2}{u-1}

化简一下,得到:

(11u)du=1xdx(1-\frac{1}{u})du=\frac{1}{x}dx

两边反求导,得到:

ulnu=lnx+Cu-\ln{u}=\ln{x}+C

两个ln没法处理,我们把两边都变成ee次幂,得到:

euu=Cx\frac{e^u}{u}=C\cdot x

这时候将uu替换成yx\frac{y}{x},得到:

eyxy=C\frac{e^\frac{y}{x}}{y}=C

化简得到答案:

y=Ceyxy=Ce^{\frac{y}{x}}

一阶线性微分方程

好,接下来是我们的重点了:我们需要求解形如y+P(x)y=Q(x)y^\prime+P\left(x\right)y=Q\left(x\right)的方程。叫做一阶的原因,是y\textbf{y}及其各阶导数都是一次的

此类方程的特征在于,当Q(x)=0Q\left(x\right)=0时,方程y+P(x)y=0y^\prime+P\left(x\right)y=0是齐次的。我们可以令Q(x)=0Q\left(x\right)=0来获取该方程的特解,以此作为我们解题的突破口。

比如我们来看出现在我们MTH019期末考试的微分方程题,两边求导会得到:

y=2x+yy^\prime=2x+y

让我们把它变成我们上面的形式:

yy=2xy^\prime-y=2x

我们可以得到,P(x)=1P(x)=-1Q(x)=2xQ(x)=2x。接下来,我们需要令Q(x)=0Q\left(x\right)=0,来解这个一阶齐次微分方程

yy=0y^\prime-y=0

非常简单!相信你一眼就看出来了,我们得到的特解是:

y=Cexy=Ce^x

好,我们现在需要细化结果中的CC是什么。因为在原方程中,CC是会受xx影响的,而不是一个常量。我们把这里的CC重命名为C(x)C(x),带入原式后求导。

y=C(x)ex+exC(x)y^\prime=C^\prime(x)e^x+e^xC(x)

好,我们现在需要把上面的这个式子和y=C(x)exy=C(x)e^x一起带入原式:

C(x)ex+exC(x)C(x)ex=2xC^\prime(x)e^x+e^xC(x)-C(x)e^x=2x

此时等式中的C(x)C(x)一定可以被约掉。整理一下,就是:

C(x)=2xexC^\prime(x)=\frac{2x}{e^x}

然后反求导,得到C(x)C(x)这个积分你一定会算吧):

C(x)=2(x+1)ex+CC(x)=\frac{-2(x+1)}{e^x}+C

好,让我们带入y=C(x)exy=C(x)e^x,得到:

y=Cex2x2y=Ce^x-2x-2

大功告成!是不是非常简单?

以上是解答题需要如何书写具体步骤,在其他题你完全可以直接使用这个公式直接得出结果:

y=CeP(x)dx+eP(x)dxQ(x)eP(x)dxdxy=Ce^{\int P\left(x\right)dx}+e^{-\int P\left(x\right)dx}\cdot\int{Q\left(x\right)e^{\int P\left(x\right)dx}}dx

其中,我们一阶齐次微分方程的特解是:

y=CeP(x)dxy=Ce^{\int P\left(x\right)dx}

高阶微分方程

唔,我们现在还剩一种微分方程还没有学习,面对多阶的微分方程,目前只有三种题型是我们能够下手的:第一种是y(n)=f(x)y^{(n)}=f(x)型,第二种是y=f(x,y)y^{\prime\prime}=f(x,y^\prime)型,第三种是y=f(y,y)y^{\prime\prime}=f(y,y^\prime)型。

对于y(n)=f(x)y^{(n)}=f(x)型,对等式两边进行nn次积分就能获得答案。需要注意答案中包含多个CC,需写为C1C_1,C2C_2的形式。

对于y=f(x,y)y^{\prime\prime}=f(x,y^\prime)型,由于等式中不包含yy,我们令y=py^\prime=p,这样y=py^{\prime\prime}=p^\prime以此来减少一层求导,得出ppxx的等式之后,再积分一次就能获得答案啦。

对于y=f(y,y)y^{\prime\prime}=f(y,y^\prime)型,这里我们也令y=py^\prime=p,但我们令y=pdpdyy^{\prime\prime}=p\frac{dp}{dy},就比如这个微分方程:

yy(y)2=0yy^{\prime\prime}-(y^\prime)^2=0

我们改写之后就成为了:

ypdpdyp2=0yp\frac{dp}{dy}-p^2=0

两边同时除以pp,整理一下得到:

1pdp=1ydy\frac{1}{p}dp=\frac{1}{y}dy

两边积分得:

y=p=C1yy^\prime=p=C_1y

所以结果为:

y=C2eC1xy=C_2e^{C_1x}