求解微分方程

最新推荐文章于 2024-05-10 09:00:00 发布

云端暮雪

最新推荐文章于 2024-05-10 09:00:00 发布

阅读量9.9k

点赞数 16

分类专栏：数学基础文章标签：高等数学微分方程

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/yunduanmuxue/article/details/11924089

版权

数学基础专栏收录该内容

4 篇文章 3 订阅

订阅专栏

一般地，含有未知函数及未知函数的导数或微分的方程称为微分方程。微分方程中出现的未知函数的最高阶导数的阶数称为微分方程的阶。

按照不同的分类标准，微分方程可以分为线性或非线性，齐次或非齐次。

一般地，微分方程的不含有任意常数的解称为微分方程的特解，含有相互独立的任意常数，且任意常数的个数与微分方程阶数相等的解称为微分方程的通解（一般解）。

下面介绍微分方程的求解方法。

一、一阶微分方程

一阶微分方程具有如下一般形式：

$\frac{dy}{dx}=f(x,y)$

1.可分离变量

这类方程可以化为如下形式：

$\frac{dy}{dx}=f(x)g(y)$

设 $g(y)\neq0$ ，可以通过下式求解

$\int\frac{1}{g(y)}dy=\int f(x)dx$

如果，则易知也是方程的解。

2.齐次方程

形如

$\frac{dy}{dx}=f(\frac{y}{x})$

的一阶微分方程称为齐次微分方程，简称齐次方程。

通过变量替换，可以将这类方程化为可分离变量的方程来求解。

令

或

其中是新的未知函数，则有

$\frac{dy}{dx}=u+x\frac{du}{dx}$

将式代入原方程，得

$u+x\frac{du}{dx}=f(u)$

分离变量，两边积分，有

$\int \frac{du}{f(u)-u}=\int \frac{dx}{x}$

求出积分后，再将回代，便得到方程的解。

3.可化为齐次方程的微分方程

有些问题本身虽然不是齐次的，但是通过适当变换，可以化为齐次方程。

例如，对于形如

$\frac{dy}{dx}=f\left(\frac{a_1x+b_1y+c_1}{a_2x+b_2y+c_2} \right )$

的方程，先求出两条曲线

$a_1x+b_1y+c_1=0, a_2x+b_2+c_2=0$

的交点，于是作平移变换

$\begin{cases} X=x-x_0\\ Y=y-y_0 \end{cases}$

这时，原方程可以化为齐次方程

$\frac{dY}{dX}=f\left(\frac{a_1X+b_1Y}{a_2X+b_2Y} \right )$

4.一阶线性微分方程

形如

$\frac{dy}{dx}+P(x)y=Q(x)$

的方程称为一阶线性微分方程。

当 $Q(x)\equiv0$ 时，上式变为

$\frac{dy}{dx}+P(x)y=Q(x)$

这个方程称为一阶齐次线性微分方程。

一阶齐次线性微分方程是可分离变量的方程，由上面的方法可以得到方程的通解为

$y = Ce^{-\int P(x)dx}$

其中 $C(C=\pm e^{C_1})$ 是任意常数。

下面再讨论一阶齐次非线性方程的通解。

将方程变形为

$\frac{dy}{y}=\left[\frac{Q(x)}{y}-p(x) \right]dx$

两边积分，得

$ln|y|=\int \frac{Q(x)}{y}dx+\int P(x)dx$

若记 $\int \frac{Q(x)}{y}dx = v(x)$ ，则

$ln|y| = v(x) - \int P(x) dx$

即

$y = \pm e^{v(x)}e^{-\int P(x)dx}=u(x)e^{-\int P(x) dx)}$

与齐次方程的通解相比较，易见其表达形式一致，只需将换为函数。由此可以引出求解一阶非齐次线性微分方程的常数变异法。即在求出齐次方程通解后，将通解中的常数变异为待定函数，并设一阶非齐次方程的通解为

$y = u(x)e^{-\int P(x) dx}$

求导，得

$y' = u'e^{-\int P(x) dx} + u[-P(x)]e^{-\int P(x) dx}$

将和代入原方程，得

$u'(x) e^{-\int P(x) dx}=Q(x)$

积分，得

$u(x) = \int Q(x) e^{\int P(x) dx}dx+C$

从而得到一阶非齐次线性微分方程的通解为

$y = \left[\int Q(x)e^{\int P(x)dx}dx+C\right] e^{-\int P(x)dx}$

上述公式可以写成

$y =Ce^{-\int P(x)dx} + e^{-\int P(x)dx}\cdot \int Q(x)e^{\int P(x)dx}dx$

从中可以看出，一阶非齐次线性微分方程的通解是对应的齐次线性方程的通解与其本身的一个特解之和。这个结论对高阶非齐次线性方程亦成立。

5.伯努利方程

形如

$\frac{dy}{dx} + P(x)y = Q(x)y^n$

的方程称为伯努利方程。其中为常数，且 $n\neq 0,1$ 。

伯努利方程是一类非线性方程，但是通过适当的变换，就可以把它转化为线性方程。

方程两端同时除以，得

$\frac{1}{1-n}\cdot (y^{1-n})' + P(x)y^{1-n}=Q(x)$

于是，令 $z = y^{1-n}$ ，就得到关于变量的一阶线性方程

$\frac{dz}{dx} + (1-n)P(x)z = (1-n)Q(x)$

利用线性方程的求解方法求出通解后，再回代原变量，便可得到伯努利方程的通解

$y^{1-n} = e^{\int (1-n)P(x)dx} \left( \int Q(x)(1-n)e^{\int (1-n)P(x)dx}dx+C \right)$

二、特殊二阶微分方程

本节介绍三种可以通过化简来求解的二阶微分方程形式。

1. 型

这是最简单的二阶微分方程，求解方法是逐次积分，得到

$y = \int \left[ \int f(x)dx + C_1 \right ]dx + C_2$

这类方程的解法，可以推广到阶微分方程

$y^{(n)} = f(x)$

只要连续积分次，就可以得到这个方程含有个任意常数的通解。

2. 型

这类方程的特点是不显含未知函数，求解方法是：

令，则，原方程化为以为未知函数的一阶微分方程

设其通解为

$p = \varphi (x,C_1)$

然后回代变量，又可以得到一个一阶微分方程

$\frac{dy}{dx} = \varphi(x,C_1)$

对它进行积分，可以得到方程的通解

$y = \int \varphi(x,C_1) dx + C_2$

3. 型

这类方程的特点是不显含未知函数，求解方法是：

把暂时看成自变量，并作变换，于是，由符合函数求导法则，有

$y'' = \frac{dp}{dx}=\frac{dp}{dy}\cdot\frac{dy}{dx}=p\frac{dp}{dy}$

这样将原方程化为

$p\frac{dp}{dy} = f(y,p)$

这是一个关于变量的一阶微分方程，设它的通解为

$y' = p =\varphi(y,C_1)$

这是可分离变量的方程，对其积分得到原方程的通解

$\int \frac {dy}{\varphi(y,C_1)}=x + C_2$

三、二阶线性微分方程

二阶线性微分方程的一般形式是

$\frac{d^2y}{dx^2} + P(x)\frac{dy}{dx} + Q(x)y = f(x)$

函数成为方程的自由项。

当时，方程变为

$\frac{d^2y}{dx^2} + P(x)\frac{dy}{dx} + Q(x)y = f(x)$

后者称为二阶齐次线性微分方程，前者称为二阶非齐次线性微分方程。

下面介绍几个重要的定理。

定理一：如果函数与是方程 (2) 的两个解，则

也是方程 (2) 的解，其中是任意常数。

这个性质表明其次线性方程的解符合叠加定理。

补充一个定义：如果是定义在区间内的两个函数。如果存在两个不全为零的常数，使得在区间内恒有

$k_1y_1(x) + k_2y_2(x)\equiv 0$

则称这两个函数在区间内线性相关，否则称为线性无关。

定理二：如果函数与是方程 (2) 的两个线性无关的特解，则

也是方程 (2) 的通解，其中是任意常数。

定理三：设是方程 (1) 的一个特解，而是其对应的齐次方程 (2) 的通解，则

是二阶非齐次线性微分方程 (1) 的通解。

定理四：设分别是方程

的特解，则是方程的特解。

这个定理通常称为非齐次线性微分方程的解的叠加定理。

定理五：设是方程

的解，其中为实值函数，为纯虚数，则分别是方程

的解。

在方程 (1) 中，系数是随变化的，这类方程求解比较困难，下面介绍处理这类方程的两种方法。

1.降阶法

考虑二阶齐次线性方程

$\frac{d^2y}{dx^2} + P(x)\frac{dy}{dx} +Q(x)y = 0$

设是方程的一个已知的非零特解，作变量替换

其中为待定系数，求的一阶导数和二阶导数，得

$\frac{dy}{dx} = y_1\frac{du}{dx} + u\frac{dy_1}{dx},\frac{d^2y}{dx^2}=y_1\frac{d^2u}{dx^2}+2\frac{du}{dx}\frac{dy_1}{dx}+u\frac{d^2y_1}{dx^2}$

将他们代入方程，得

$y_1\frac{d^2y}{dx^2}+\left(2\frac{dy_1}{dx}+P(x)y_1 \right )\frac{du}{dx}+\left(\frac{d^2y}{dx^2}+P(x)\frac{dy_1}{dx}+Q(x)y_1 \right )u = 0$

各项系数是的已知函数，因为是原方程的解，所以，其中的系数为 0 。

故上式可以化为

$y_1\frac{d^2y}{dx^2}+\left(2\frac{dy_1}{dx}+P(x)y_1 \right )\frac{du}{dx} = 0$

再做变量替换 $\frac{du}{dx} = z$ ，得

$y_1\frac{dz}{dx}+\left(2\frac{dy_1}{dx}+P(x)y_1 \right )z=0$

分离变量

$\frac{1}{z}dz=-\left[\frac{2}{y_1}\frac{dy_1}{dx}+P(x) \right ]dx$

两边积分，得其通解

$z=\frac{C_2}{y^2_1}e^{-\int P(x)dx}$

其中为任意常数。

对 $\frac{du}{dx} = z$ 积分，得

$u=C_2\int\frac{1}{y^2_1}e^{-\int P(x)dx}dx+C_1$

其中为任意常数。

代回原变量，就得到原方程的通解

$y=y_1\left[C_1 +C_2\int\frac{1}{y^2_1}e^{-\int P(x)dx}dx \right ]$

这个公式称为二阶线性微分方程的刘维尔公式。

综上所述，对于二阶齐次线性方程，如果已知其一个非零特解，作变量替换 $y=y_1\int zdx$ 就可将其姜维一阶齐次线性方程，从而求得通解。

对于二阶非齐次线性方程，若已知其对应的齐次方程的一个特解，作同样的变量替换（因为这种变换并不影响方程的右端），也能使非齐次方程降一阶。

2.常数变异法

设有二阶非齐次线性方程

$\frac{d^2y}{dx^2}+P(x)\frac{dy}{dx}+Q(x)y=f(x)$

其中在某区间上连续，如果其对应的齐次方程

$\frac{d^2y}{dx^2} + P(x)\frac{dy}{dx} +Q(x)y = 0$

的通解已经求得，可通过如下方法求得其通解。

设非齐次方程具有形如

的特解，对求导数，得

把特解代入原方程，可以得到确定的一个方程，因为这里有两个未知函数，所以还需要添加另一个条件，为计算发辫，我们补充如下条件：，这样

代入原方程，并注意到是齐次方程的解，经整理得

与补充条件联立，得方程组

$\begin{cases} y_1u_1'+y_2u_2'=0\\ y_1'u_1'+y_2'u_2'=f(x) \end{cases}$

因为线性无关，即 $\frac{y_2}{y_1}\neq$ 常数，所以

$\left(\frac{y_2}{y_1} \right )'=\frac{y_1y_2'-y_2y_1'}{y_1^2}\neq0$

设则有 $w(x))\neq0$ ，所以上述方程有唯一解，解得

$\begin{cases} u_1'=\frac{-y_2f(x)}{y_1y_2'+y_2y_1'}=\frac{-y_2f(x)}{w(x)}\\ u_2'=\frac{y_1f(x)}{y_1y_2'+y_2y_1'}=\frac{y_1f(x)}{w(x)} \end{cases}$

积分并取一个原函数，得

$\begin{cases} u_1=\int\frac{-y_2f(x)}{w(x)}dx\\ u_2'=\int\frac{y_1f(x)}{w(x)}dx \end{cases}$

于是，所求特解为

$y^*=y_1\int\frac{-y_2f(x)}{w(x)}dx+y_2\int\frac{y_1f(x)}{w(x)}dx$

所以，所求方程的通解为

$y=Y+y^*=C_1y_1+C_2y_2+y_1\int\frac{-y_2f(x)}{w(x)}dx+y_2\int\frac{y_1f(x)}{w(x)}dx$

关注

16
点赞
踩
50

收藏

觉得还不错? 一键收藏
2
评论
求解微分方程

一、一阶微分方程具有如下形式：1.可分离变量具有如下形式：设，可以得到如果，则易知也是方程的解。
复制链接

扫一扫

专栏目录

评论 2

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

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