常微分方程2

作者: 夏羽兮 | 来源:发表于2020-06-09 10:33 被阅读0次

10-Python 科学计算_scipy 篇
常/偏微分方程求解
【常微分方程】$1.1 微分方程及其解的定义
常微分方程2
常微分方程2
几道微分方程题
二阶变系数线性常微分方程的求解
前向欧拉方程
常微分方程
常微分方程

一阶微分方程的解法：

（一）变量分离方程与变量变换

1.1形如 $\frac{dy}{dx} =f(x)\varphi (x)$ 的方程，称为变量分离方程，这里 $f(x)$ , $\varphi (y)$ 分别是 $x$ , $y$ 的连续函数。

（1）如果 $\varphi （y）≠0$ ,则 $\frac{dy}{\varphi(y)} =f(x)dx$ ,两边进行积分，得到 $\int_{}^{} \frac{dy}{\varphi (y)} =\int_{}^{} f(x)dx+c$ .

（2）如果 $\varphi （y）=0$ ,但如果存在 $y_{0}$ 使 $\varphi （y_{0} ）=0$ ,则直接验证 $y=y_{0}$ 也是方程的解。

例1

1.2可化为变量分离类型的方程

形如 $\frac{dy}{dx} =g(\frac{y}{x} )$ 的方程，称为齐次线性微分方程，这里 $g(u)$ 是 $u$ 的连续函数。

作变量变换 $u=\frac{y}{x}$ ，即 $y=ux$

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

代入原方程，则方程变为 $x\frac{du}{dx} +u=g(u)$

整理后得， $\frac{du}{dx} =\frac{g (u)-u}{x}$ ，是一个变量分离方程

代回原变量，即为方程的解。

例2

形如 $\frac{dy}{dx} =\frac{a_{1}x+b_{1}y+c_{1} }{a_{2}x+b_{2}y+c_{2} }$ 的方程也可经变量变换化为变量分离方程，这里 $a_{1}$ , $a_{2}$ , $b_{1}$ , $b_{2}$ , $c_{1}$ , $c_{2}$ 均为常数。

分3种情形讨论：

（1) $\frac{a_{1} }{a_{2} } =\frac{b_{1} }{b_{2} } =\frac{c_{1} }{c_{2} } =k$ (常数)情形。

这时方程化为 $\frac{dy}{dx} =k$ ，有通解 $y=kx+c$ ，其中 $c$ 为任意常数。

（2) $\frac{a_{1} }{a_{2} } =\frac{b_{1} }{b_{2} } =k\neq \frac{c_{1} }{c_{2} }$ 情形。

令 $u=a_{2}x+b_{2} y$ ，这时有 $\frac{du}{dx} =a_{2} +b_{2} \frac{dy}{dx} =a_{2} +b_{2} \frac{ku+c_{1} }{u+c_{2} }$ 是变量分离方程。

（3） $\frac{a_{1} }{a_{2} } \neq \frac{b_{1} }{b_{2} }$ 情形。

令 $a_{1} x+b_{1} y+c_{1} =0$ 且 $a_{2} x+b_{2} y+c_{2} =0$ ，可得两直线的交点 $（\alpha ，\beta ）$

再令 $X=x-\alpha$ 且 $Y=y-\beta$

则上式化为 $a_{1} X+b_{1} Y=0$ 且 $a_{2} X+b_{2} Y=0$

从而原方程变为 $\frac{dY}{dX} =\frac{a_{1}X+b_{1}Y }{a_{2}X+b_{2}Y } =g(\frac{Y}{X} )$

求解上述变量分离方程，代回原变量可得原方程的解。

（注：若方程中 $c_{1} =c_{2} =0$ ，直接取变换 $u=\frac{y}{x}$ 即可）

例3

10-Python 科学计算_scipy 篇
课程概要：1、非线性方程组2、数值积分3、常微分方程组 1、非线性方程组 2、数值积分（定积分） 3、常微分方程组
常/偏微分方程求解
常微分方程 1、四阶经典步长Runge-Kutta法求解高阶常微分方程 2、用ode15s函数、ode23s函数和...
【常微分方程】$1.1 微分方程及其解的定义
1.常微分方程、偏微分方程概念常微分方程：未知函数是一元函数。偏微分方程：未知函数是多元函数，微分方程...
常微分方程2
常微分方程2
一阶微分方程的解法：（一）变量分离方程与变量变换 1.1形如的方程，称为变量分离方程，这里,分别是,的连续函数。...
几道微分方程题
13. 求欧拉方程：令，原方程化为即，。此方程为二阶常系数齐次线性常微分方程。其特征方程为，所以，所以...
二阶变系数线性常微分方程的求解
之前总结过一次常见的常微分方程的解法，不过里面并没有二阶变系数线性常微分方程的求解方法，现在补充如下。对于形式的...
前向欧拉方程
欧拉方法是一种一阶数值方法，用以对给定初值的常微分方程（即初值问题）求解。它是一种解决数值常微分方程的最基本的一类...
常微分方程
import numpy as np import matplotlib.pyplot as plt from s...
常微分方程
一、一阶线性微分方程一般形式：当时，则方程为齐次方程，反之为非齐次。通解：通解公式：特解：特解是将通解中的根...