一阶线性微分方程、条件多阶线性微分方程通解公式

一阶线性微分方程: 

一阶线性微分方程标准形式为:

                                                \frac{\mathrm{d} y}{\mathrm{d} x}+P\left ( x \right )y=Q\left ( x \right )

其中P(x)、Q(x)均为连续函数,改方程通解为:y=e^{-\int P(x)dx}\left ( \int Q(x)e^{\int P(x)dx}dx+c \right )

下面是推导过程:

首先我们都知道假如知道了f(x)的导数方程,f(x)只用通过积分就可以得到,即

\because y^{(1)}=f(x)

\therefore f(x)=\int ydx+c

所以我们要得到f(x),就要知道其导数形式,但大多数情况我们不能直接看出,所以需要我们将没见过的形式转换成我们知道的能积分出来的形式。标准形式为y^{'}+P\left ( x \right )y=Q\left ( x \right ),很容易发现这与我们学过的(uv)^{'}=u^{'}v+uv^{'}十分相似,下面的工作就是找到u和v。

寻找过程中我们需要用到两个重要概念:(e^{})^{'}=e^{x};

(\int f(x)dx)^{'}=f(x);

将两者结合可以达到这样的效果:(e^{\int f(x)dx})^{'}=e^{\int f(x)dx}f(x)

将标准式两端同乘e^{\int P(x)dx}得到:y^{'}e^{\int P(x)dx}+P\left ( x \right )ye^{\int P(x)dx}=Q\left ( x \right )e^{\int P(x)dx}

即得到我们想要的:(e^{\int P(x)dx}y)^{'}=y^{'}e^{\int P(x)dx}+P\left ( x \right )ye^{\int P(x)dx}

(ye^{\int P(x)dx})^{'}=Q\left ( x \right )e^{\int P(x)dx}

积分得到:ye^{\int P(x)dx}=\int Q\left ( x \right )e^{\int P(x)dx}dx+c

左右同除e^{\int P(x)dx}e^{\int P(x)dx}不为0)得到:y=e^{-\int P(x)dx}\left ( \int Q(x)e^{\int P(x)dx}dx+c \right )

即得到通解,而c的值确定需要根据题目所给条件确定。

多阶线性微分方程(只包括最高阶和次高阶)

理解了一阶线性微分方程通解推导以后,多阶线性微分方程分析也是同样的方法。

以二阶为例,二阶标准形式为A(x)y^{''}+B(x)y^{'}=Q\left ( x \right ) ,

u=y^{'},则y^{''}=u^{'}

原式=A(x)u^{'}+B(x)u=Q\left ( x \right )

u^{'}+\frac{B(x)}{A(x)}u=\frac{Q(x)}{A(x)}

u^{'}e^{\int \frac{B(x)}{A(x)}dx}+\frac{B(x)}{A(x)}ue^{\int\frac{B(x)}{A(x)}dx}=\frac{Q\left ( x \right )}{A(x)}e^{\int \frac{B(x)}{A(x)}dx}

ue^{\int \frac{B(x)}{A(x)}dx}=\int \frac{Q\left ( x \right )}{A(x)}e^{\int \frac{B(x)}{A(x)}dx}dx+c1

u=e^{-\int \frac{B(x)}{A(x)}dx}(\int \frac{Q\left ( x \right )}{A(x)}e^{\int \frac{B(x)}{A(x)}dx}dx+c1)

y=\int udx+c2=\int e^{-\int \frac{B(x)}{A(x)}dx}(\int \frac{Q\left ( x \right )}{A(x)}e^{\int \frac{B(x)}{A(x)}dx}dx+c1)dx+c2

往往需要两个初始条件解得c1、c2。

例题

最近在准备数学竞赛,选了一道竞赛题。

解:\frac{\mathrm{d} y}{\mathrm{d} x}=\frac{\mathrm{d} y}{\mathrm{d} t}\frac{\mathrm{d} t}{\mathrm{d} x}=\frac{\varphi (t)^{'}}{2+2t}

\frac{\mathrm{d} \frac{\mathrm{d} y}{\mathrm{d} x}}{\mathrm{d} t}\frac{\mathrm{d} t}{\mathrm{d} x}=\frac{\varphi (t)^{''}(2+2t)-2\varphi (t)^{'}}{(2+2t)^{3}}=\frac{3}{4(1+t)}

解得:\varphi (t)^{''}-\frac{\varphi (t)^{'}}{1+t}=3(1+t)

\varphi (t)^{'}=u,则\varphi (t)^{''}=u^{'}

u^{'}-\frac{u}{1+t}=3(1+t)

u=e^{\int \frac{1}{1+t}dt}(\int 3(1+t)e^{-\int \frac{1}{1+t}dt}dt+c1)=(1+t)(3t+c1)

由题意得,t=1时,\varphi (1)=\frac{3}{2e}    (1)

\varphi (1)^{'}=\frac{2}{e}      (2)

将(1)式带入u得c1=\frac{1}{e}-3

\varphi (t)=\int udt+c2

将c2带入\varphi (t)得c2=2

\therefore \varphi (t)=t^{3}+\frac{t^{2}}{2e}+(\frac{1}{e}-3)t+2   (t>-1)

  • 1
    点赞
  • 2
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值