一阶常微分方程

这篇具有很好参考价值的文章主要介绍了一阶常微分方程。希望对大家有所帮助。如果存在错误或未考虑完全的地方,请大家不吝赐教,您也可以点击"举报违法"按钮提交疑问。

前言

第一次写博客记录自己的学习过程,写的不好希望大家斧正。

讲的更多的是常微分方程的解法的理解,也是我在学习中遇到各种证明的关键点,希望通过记录博客深化对于证明的理解,建立起数学思维,而不是知其然而不知其所以然。因此阅读中需要读者有一定的基础,与实践相结合,

2.1 分离变量和变量代换法

2.1.1 分离变量法

一阶常微分方程
我们很容易对上式进行变量分离,并计算;需要注意的是:分离变量时f(x)或者g(y)做分母时不为零,此时若有f(x0) = 0 、g(y0) = 0,则方程可能还存在着特殊解y = y0、x = x0。

2.1.2 变量代换法

1.齐次方程

一阶常微分方程

2.类似于两直线方程相除的形式

一阶常微分方程
该方程求解的时候分为3种情况:

  • c1、c2为零的情况:直接就是齐次方程
  • 两直线方程平行的情况,且c1、c2不同时为0:令u = a1 * x + b1 * y,然后进行变量分离
  • 两直线有焦点的情况:求出交点,然后将交点移动到原点的位置,此时有新的变量X、Y,并且方程可以化成其次方程

将求解分为这三类主要就是将方程向齐次方程或者可分离变量方程的形式靠拢。

3. 第三类

一阶常微分方程

2.2 常数变易法

常数变易法是专门针对一阶线性方程、高阶线性方程和线性方程组总结的一种特定的方法;这里先讲一阶的,后面会有相应高阶方程应用该方法进行求解。

2.2.1 线性方程的通解公式

一阶常微分方程
上图中方程是一阶非齐次的常微分方程,对于其齐次形式,即 q ( x ) = 0 q(x) = 0 q(x)=0,我们很容易可以解出:

y = C e ∫ p ( x ) d x . y = Ce^{\int{p(x)dx}}. y=Cep(x)dx.

那么对于非齐次形式的方程我们猜想解为: y = C ( x ) e ∫ p ( x ) d x y = C(x)e^{\int{p(x)dx}} y=C(x)ep(x)dx
后面进行验证,有两种方法:

  • y = C ( x ) e ∫ p ( x ) d x y = C(x)e^{\int{p(x)dx}} y=C(x)ep(x)dx带入到原方程中,解出 C ( x ) C(x) C(x)
  • 方程 y = C ( x ) e ∫ p ( x ) d x y = C(x)e^{\int{p(x)dx}} y=C(x)ep(x)dx两边对 x x x求导,然后求解(高阶常微分方程的常数变易法与此类似)

结果均为:
y = e ∫ p ( x ) d x ( C + ∫ q ( x ) e ∫ − p ( x ) d x d x ) y = e^{\int{p(x)dx}}(C + \int{q(x)e^{\int{-p(x)dx}}dx}) y=ep(x)dx(C+q(x)ep(x)dxdx)

2.2.2 伯努利方程

方程如下:
d y d x + p ( x ) y = q ( x ) y n ( n ≠ 0 , 1 ) \frac{\mathrm{d} y}{\mathrm{d}x} + p(x)y = q(x)y^n \quad(n \neq 0,1) dxdy+p(x)y=q(x)yn(n=0,1)
这种类型的方程同样也化成常数变易的形式,将方程两端同时除以 y n y^n yn,化为如下格式:
y − n d y d x + p ( x ) y 1 − n = q ( x ) y^{-n}\frac{\mathrm{d} y}{\mathrm{d}x} + p(x)y^{1-n} = q(x) yndxdy+p(x)y1n=q(x)
此时令 z = y 1 − n z = y^{1-n} z=y1n,常数变易法即可解出y。

2.3 凑全微分法与积分因子法

2.3.1 恰当方程

当方程
M ( x , y ) d x + N ( x , y ) d y = 0 , ( 3 − 0 ) \quad M(x,y)\mathrm{d}x + N(x,y)\mathrm{d}y = 0\quad,\qquad(3-0) M(x,y)dx+N(x,y)dy=0,(30)
等号左边可以表示为某二元函数的全微分 M d x + N d y = d u \quad M\mathrm{d}x + N\mathrm{d}y = \mathrm{d}u\quad Mdx+Ndy=du时,称此一阶方程为恰当方程。

  • 定理:当 ∂ M ∂ x ≡ ∂ N ∂ y \quad\frac{\partial M}{\partial x} \equiv \frac{\partial N}{\partial y}\quad xMyN时,方程时恰当方程。
    积分时先对N或者M偏积分,然后对剩下的未积分完全的变量在进行积分:
    一阶常微分方程
    证明过程:
    ∂ u ∂ y = ∂ ∂ y u = ∂ ∂ y ∫ x 0 x M ( x , y ) d x + ∂ ψ ( y ) ∂ y \frac{\partial{u}}{\partial{y}} =\frac{\partial}{\partial{y}}u= \frac{\partial}{\partial{y}}\int_{x0}^x{M(x,y)dx} + \frac{\partial{\psi(y)}}{\partial{y}} yu=yu=yx0xM(x,y)dx+yψ(y) = ∫ x 0 x ∂ ∂ y M ( x , y ) d x + ψ ′ ( y ) =\int_{x0}^x{\frac{\partial}{\partial{y}}M(x,y)dx} + \psi{'}(y) =x0xyM(x,y)dx+ψ(y) = ∫ x 0 x ∂ ∂ x N ( x , y ) d x + ψ ′ ( y ) =\int_{x0}^x{\frac{\partial}{\partial{x}}N(x,y)dx} + \psi{'}(y) =x0xxN(x,y)dx+ψ(y) = N ( x , y ) − N ( x 0 , y ) + ψ ′ ( y ) =N(x,y) - N(x0,y) + \psi{'}(y) =N(x,y)N(x0,y)+ψ(y)
    容易得到: N ( x 0 , y ) = ψ ′ ( y ) N(x0,y) = \psi{'}(y) N(x0,y)=ψ(y),进而求出$\psi(y) = $。
    积分为:
    u ≡ ∫ x 0 x M ( x , y ) d x + ∫ y 0 y N ( x 0 , y ) d y = C . u \equiv \int_{x0}^xM(x,y)dx + \int_{y0}^yN(x0,y)dy = C. ux0xM(x,y)dx+y0yN(x0,y)dy=C.

2.3.2 凑全微分法

一阶常微分方程

2.3.3 积分因子法

1.常规情况

∂ M ∂ x ≠ ∂ N ∂ y \quad\frac{\partial M}{\partial x} \neq \frac{\partial N}{\partial y}\quad xM=yN时,方程不是恰当的,不能表示为一个二元函数的全微分。我们希望找到一个函数 μ ( x , y ) \mu(x,y) μ(x,y),使得:
d u ( x , y ) ≡ μ ( x , y ) [ M ( x , y ) d x + N ( x , y ) d y ] du(x,y) \equiv \mu(x,y)[M(x,y)dx + N(x,y)dy] du(x,y)μ(x,y)[M(x,y)dx+N(x,y)dy]
其中 μ ( , x , y ) \mu(,x,y) μ(,x,y)为方程的积分因子。此时方程是恰当的,有以下推导:
∂ ( μ M ) ∂ y = ∂ ( μ N ) ∂ x \frac{\partial(\mu M)}{\partial y} = \frac{\partial(\mu N)}{\partial x} y(μM)=x(μN) μ ( ∂ M ∂ y − ∂ N ∂ x ) = N ∂ μ ∂ x − M ∂ μ ∂ y , ( 3 − 1 ) \mu(\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}) = N\frac{\partial \mu}{\partial x} - M\frac{\partial \mu}{\partial y} ,\quad(3-1) μ(yMxN)=NxμMyμ,(31)
从一般的偏微分方程求解 μ ( x , y ) \mu(x,y) μ(x,y)会比较困难,先不做考虑;这里我们假设 μ = μ ( x ) 或 μ = μ ( y ) \mu = \mu(x) 或\mu = \mu(y) μ=μ(x)μ=μ(y),此时分离变量法即可求出 μ \mu μ
μ ( ∂ M ∂ y − ∂ N ∂ x ) N = ∂ μ ∂ x , ( 3 − 2 ) \frac{\mu(\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x})}{N} = \frac{\partial \mu}{\partial x},\qquad(3-2) Nμ(yMxN)=xμ,(32) 积 分 得 μ ( x ) = e ∫ ∂ M ∂ y − ∂ N ∂ x N , ( 3 − 3 ) 积分得\mu(x) = e^{\int{{\frac{\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}}{N}}}},\qquad(3-3) μ(x)=eNyMxN,(33)
μ ( ∂ M ∂ y − ∂ N ∂ x ) − M = ∂ μ ∂ y , ( 3 − 4 ) \frac{\mu(\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x})}{-M} = \frac{\partial \mu}{\partial y},\qquad(3-4) Mμ(yMxN)=yμ,(34) 积 分 得 μ ( y ) = e ∫ ∂ M ∂ y − ∂ N ∂ x − M , ( 3 − 5 ) 积分得\mu(y) = e^{\int{{\frac{\frac{\partial M}{\partial y} - \frac{\partial N}{\partial x}}{-M}}}},\qquad(3-5) μ(y)=eMyMxN,(35)
在实际计算中应该注意寻找(3-2)、(3-4),即可解决这类问题。

2. 一阶微分方程可以拆分时

先引入定理:

  • μ ( x , y ) \mu(x,y) μ(x,y)是方程的积分因子,且 u = u ( x , y ) u = u(x,y) u=u(x,y),则 μ ⋅ g ( u ) \mu ·g(u) μg(u)也是方程(3-0)的一个积分因子,简单推导一下充分性:
    已 知 : d u = μ M d x + μ N d y , 已知:du = \mu Mdx + \mu Ndy, du=μMdx+μNdy, 故 只 需 证 明 : ∂ [ μ ⋅ g ( u ) M ] ∂ y = ∂ [ μ ⋅ g ( u ) N ] ∂ x , 故只需证明:\frac{\partial [\mu·g(u) M]}{\partial y}= \frac{\partial [\mu·g(u) N]}{\partial x}, :y[μg(u)M]=x[μg(u)N] ∂ [ μ ⋅ g ( u ) M ] ∂ y = ∂ ( μ M ) ∂ y g + ∂ g ∂ y μ M = ∂ ( μ M ) ∂ y g + ∂ g ∂ u ( ∂ u ∂ y μ M ) , ( 1 ) \frac{\partial [\mu·g(u) M]}{\partial y} = \frac{\partial (\mu M)}{\partial y} g + \frac{\partial g}{\partial y}\mu M = \frac{\partial (\mu M)}{\partial y} g + \frac{\partial g}{\partial u}(\frac{\partial u}{\partial y}\mu M ) ,(1) y[μg(u)M]=y(μM)g+ygμM=y(μM)g+ug(yuμM),(1) ∂ [ μ ⋅ g ( u ) N ] ∂ x = ∂ ( μ N ) ∂ x g + ∂ g ∂ x μ N = ∂ ( μ N ) ∂ x g + ∂ g ∂ u ( ∂ u ∂ x μ N ) , ( 2 ) \frac{\partial [\mu·g(u) N]}{\partial x} = \frac{\partial (\mu N)}{\partial x} g + \frac{\partial g}{\partial x}\mu N = \frac{\partial (\mu N)}{\partial x} g + \frac{\partial g}{\partial u}(\frac{\partial u}{\partial x}\mu N),(2) x[μg(u)N]=x(μN)g+xgμN=x(μN)g+ug(xuμN),(2) ( 2 ) 、 ( 1 ) 两 式 中 对 应 项 相 等 , 充 分 性 得 证 。 (2)、(1)两式中对应项相等,充分性得证。 (2)(1)
  • 当方程(3-0)左边可以拆分成两项:
    ( M 1 d x + N 1 d y ) + ( M 2 d x + N 2 d y ) = 0 , (M_{1}dx +N_{1}dy) +(M_{2}dx +N_{2}dy) = 0, (M1dx+N1dy)+(M2dx+N2dy)=0
    μ 1 、 μ 2 \mu_{1}、\mu_{2} μ1μ2分别为积分因子时,如果能够找到共同的积分因子:
    μ ≡ μ 1 g 1 ( u 1 ) = μ 2 g 2 ( u 2 ) , \mu \equiv\mu_{1}g_{1}(u_{1}) = \mu_{2}g_{2}(u_{2}), μμ1g1(u1)=μ2g2(u2),
    μ \mu μ就是其共同的积分因子。
    拆分成两部分确实可以降低寻找积分因子的难度。

2.4 引入参数法

在处理一些隐式方程时常常借助引入参数的方法解决问题,常见的4种类型的隐式方程如下:
一阶常微分方程

2.4.1 可解出未知函数(或自变量)的方程

1. y = f ( x , d y d x ) y = f(x,\frac{dy}{dx}) y=f(x,dxdy)

d y d x = p \frac{dy}{dx} = p dxdy=p,带入到上式中得 p = ∂ f ∂ x + ∂ f ∂ p ∂ p ∂ x p=\frac{\partial f}{\partial x} + \frac{\partial f}{\partial p} \frac{\partial p}{\partial x} p=xf+pfxp,这是关于x、p的导数 d p d x \frac{dp}{dx} dxdp可解出的一阶微分方程,利用前面的方法就可以可以解决。为方便理解,具体例子如下:

2. x = f ( x , d y d x ) x =f(x,\frac{dy}{dx}) x=f(x,dxdy)

d y d x = p \frac{dy}{dx} = p dxdy=p,关于y求导数,带入到上式中得 1 p = ∂ f ∂ y + ∂ f ∂ p ∂ p ∂ y \frac1p=\frac{\partial f}{\partial y} + \frac{\partial f}{\partial p} \frac{\partial p}{\partial y} p1=yf+pfyp,这是关于y、p的导数 d p d y \frac{dp}{dy} dydp可解出的一阶微分方程,利用前面的方法就可以可以解决。为方便理解,具体例子如下:
一阶常微分方程
一阶常微分方程

2.4.2 不显含未知函数(或自变量)的方程

1. 不显含y的方程 F ( x , y ′ = 0 ) F(x,y{'} = 0) F(x,y=0)

对于上面的方程,令 p = y ′ p = y{'} p=y,有限制条件 d y = p d x dy =pdx dy=pdx;若将x-p看作空间中的曲线,
参数方程为 x = φ ( t ) 、 p = ψ ( t ) x = \varphi(t)、p =\psi(t) x=φ(t)p=ψ(t),可得:
d y = ψ ( t ) φ ′ ( t ) d t , dy =\psi(t)\varphi{'}(t)dt, dy=ψ(t)φ(t)dt,
不难得出参数形式的方程解:
x = φ ( t ) , y = ∫ ψ ( t ) φ ′ ( t ) d t + C x = \varphi(t), y =\int\psi(t)\varphi{'}(t)dt +C x=φ(t)y=ψ(t)φ(t)dt+C
该解法的关键是将微分方程拆解成合适的参数形式,如下面例子:
一阶常微分方程
一阶常微分方程

2. 不显含x的方程 F ( y , y ′ = 0 ) F(y,y{'} = 0) F(y,y=0)

对于上面的方程,令 p = y ′ p = y{'} p=y,有限制条件 d y = p d x dy =pdx dy=pdx;若将y-p看作空间中的曲线,
参数方程为 y = φ ( t ) 、 p = ψ ( t ) y = \varphi(t)、p =\psi(t) y=φ(t)p=ψ(t),可得:
φ ′ ( t ) d t = ψ ( t ) d x , \varphi{'}(t)dt =\psi(t)dx, φ(t)dt=ψ(t)dx,
不难得出参数形式的方程解:
y = φ ( t ) , x = ∫ φ ′ ( t ) ψ ( t ) d t + C y = \varphi(t), x =\int\frac{\varphi{'}(t)}{\psi(t)}dt +C y=φ(t)x=ψ(t)φ(t)dt+C
该解法的关键是将微分方程拆解成合适的参数形式,如下面例子:
一阶常微分方程
一阶常微分方程文章来源地址https://www.toymoban.com/news/detail-446914.html

到了这里,关于一阶常微分方程的文章就介绍完了。如果您还想了解更多内容,请在右上角搜索TOY模板网以前的文章或继续浏览下面的相关文章,希望大家以后多多支持TOY模板网!

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处: 如若内容造成侵权/违法违规/事实不符,请点击违法举报进行投诉反馈,一经查实,立即删除!

领支付宝红包 赞助服务器费用

相关文章

  • 0702可分类变量的微分方程-微分方程

    本节至第四节我们学习的都是一阶微分方程 ​ y ′ = f ( x , y ) y^{\\\'}=f(x,y) y ′ = f ( x , y ) (2-1) 一阶微分方程对称形式 p ( x , y ) d x + Q ( x , y ) d y = 0 ( 2 − 2 ) p(x,y)dx+Q(x,y)dy=0qquad (2-2) p ( x , y ) d x + Q ( x , y ) d y = 0 ( 2 − 2 ) 若以x为自变量,y为因变量,则 d y d x = − P ( x , y ) Q (

    2024年02月04日
    浏览(52)
  • 常微分方程建模R包ecode(一)——构建常微分方程系统

    常微分方程在诸多研究领域中有着广泛应用,本文希望向大家介绍笔者于近期开发的R包 ecode ,该包 采用简洁易懂的语法帮助大家在R环境中构建常微分方程 ,并便利地调用R图形接口,研究常微分方程系统的相速矢量场、平衡点、稳定点等解析性质,或进行数值模拟,进行敏

    2024年02月16日
    浏览(43)
  • 【数学建模】常微分,偏微分方程

    普通边界   已知t0时刻的初值    ode45()  龙格-库塔法 一阶,高阶都一样 如下: s(1) = y , s(2)=y\\\'  s(3) = x , s(4)=x\\\'   分段边界 非匿名函数    手写改进的ode45()函数代码 复杂边界值(即已知初始值,也知道末尾值),用bvp4c()函数 1. pdepe()函数 椭圆-抛物线型 控制方程  左边界

    2024年02月09日
    浏览(45)
  • 微分方程应用——笔记整理

    首先,根据正常思路走,化简得到式子:    不难发现,设  后面得出该方程的通解: 这里要注意什么等于这个通解   --- z 又因为该曲线过点  所以可以求出c为3 该题虽然简单,但是要注意几个问题,该定义域在第一象限,还有就是求解中变量的代换

    2024年02月11日
    浏览(40)
  • 高等数学(微分方程)

    x y ′ ′ ′ + ( y ′ ) 3 + y 4 xy\\\'\\\'\\\'+(y\\\')^3+y^4 x y ′′′ + ( y ′ ) 3 + y 4 quad quad 三阶 y ′ = 2 x y\\\'=2x y ′ = 2 x quad quad quad quad quad quad 一阶 d y = 2 x d x dy=2xdx d y = 2 x d x quad quad quad quad 一阶 ( y ′ ′ ) 5 + 2 y ′ = 3 (y\\\'\\\')^5+2y\\\'=3 ( y ′′ ) 5 + 2 y ′ = 3 quad quad quad 二阶 quad 例1: 已知

    2024年02月10日
    浏览(48)
  • matlab解微分方程

    f=@(变量) 表达式; x1为2 3 4 5;x2为3 4 5 6的情况下求解函数f的值 用“dsolve” step1: 申明自变量和因变量 syms y(x) step2:编程 得到: step1: 申明自变量和因变量 syms y(x) step2:编程 得到 step1.写函数文件 step2.主函数 相当于定义了一个新向量y,然后列 匿名函数 ,方程的 左边都是一阶

    2024年02月13日
    浏览(81)
  • MATLAB-常微分方程求解

    MATLAB中可以用来求解常微分方程(组)的函数有ode23、 ode23s、 ode23t、 ode23tb 、ode45、ode15s和odel13等,见下表。它们的具体调用方法类似,为了方便后面的描述, 在后面的介绍中将使用solver统一代替它们。 函数的具体调用方法如下。 [T,Y] =solver( odefun, tspan,y0) [T,Y] = solver( odefun,

    2024年02月02日
    浏览(50)
  • 【数学建模】常微分方程

    博客园解释 https://www.cnblogs.com/docnan/p/8126460.html https://www.cnblogs.com/hanxi/archive/2011/12/02/2272597.html https://www.cnblogs.com/b0ttle/p/ODEaid.html matlab求解常微分方程 https://www.cnblogs.com/xxfx/p/12460628.html https://www.cnblogs.com/SunChuangYu/p/13415439.html https://www.cnblogs.com/tensory/p/6590783.html 高等数学-常微分

    2024年02月16日
    浏览(44)
  • Fortran 微分方程求解 --ODEPACK

    最近涉及到使用Fortran对微分方程求解,我们知道MATLAB已有内置的函数,比如ode家族,ode15s,对应着不同的求解办法。通过查看odepack的官方文档,我尝试使用了dlsode求解刚性和非刚性常微分方程组。 首先是github网址:https://github.com/jacobwilliams/odepack 具体使用办法: 1.我使用的

    2024年02月11日
    浏览(44)
  • 【scipy 基础】--积分和微分方程

    对于手工计算来说,积分计算是非常困难的,对于一些简单的函数,我们可以直接通过已知的积分公式来求解,但在更多的情况下,原函数并没有简单的表达式,因此确定积分的反函数变得非常困难。 另外,相对于微分运算来说,积分运算则具有更多的多样性,包括不同的积

    2024年02月05日
    浏览(41)

觉得文章有用就打赏一下文章作者

支付宝扫一扫打赏

博客赞助

微信扫一扫打赏

请作者喝杯咖啡吧~博客赞助

支付宝扫一扫领取红包,优惠每天领

二维码1

领取红包

二维码2

领红包