常见一阶微分方程@可降阶的二阶微分方程@Bernoulli方程

可用变量代换法求解的一阶微分方程

可分离变量的方程

  • d y d x = f ( x ) g ( y ) g ( y ) ≠ 0 \frac{\mathrm{d}y}{\mathrm{d}x}=f(x)g(y)\quad g(y)\neq{0} dxdy=f(x)g(y)g(y)=0可以分离变量为 d y g ( y ) = f ( x ) d x \frac{\mathrm{d}y}{g(y)}=f(x)\mathrm{d}x g(y)dy=f(x)dx

  • 两边积分: ∫ d y g ( y ) = ∫ f ( x ) d x + C \int{\frac{\mathrm{d}y}{g(y)}}=\int{f(x)\mathrm{d}x}+C g(y)dy=f(x)dx+C

齐次方程

  • 形如或可化为 d y d x = f ( y x ) \frac{\mathrm{d}y}{\mathrm{d}x}=f(\frac{y}{x}) dxdy=f(xy), ( x ≠ 0 ) (x\neq{0}) (x=0)(0)的微分方程称为齐次方程
    • 令函数 u = y x u=\frac{y}{x} u=xy(1),则 y = x u y=xu y=xu,(2)对其两侧求关于 x x x的导数: d d x y \frac{\mathrm{d}}{\mathrm{d}x}y dxdy= d d x ( x u ) \frac{\mathrm{d}}{\mathrm{d}x}(xu) dxd(xu),即 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= u + x d u d x {u}+x{\frac{\mathrm{d}u}{\mathrm{d}x}} u+xdxdu(3)
      • 其中式(1)用一个新的未知函数 u u u代替了原未知函数 y y y,
      • 由于 y y y是关于 x x x得函数,所以由式(1)可知,函数 u u u也是关于 x x x的函数
      • 对(1)变形得到(2)式,再对(2)式两边求导,得到(3)式
    • 将(1)代入(0),得 d y d x = f ( u ) \frac{\mathrm{d}y}{\mathrm{d}x}=f(u) dxdy=f(u)(4),将(4)代入(3),得 f ( u ) = u + x d u d x f(u)=u+x\frac{\mathrm{d}u}{\mathrm{d}x} f(u)=u+xdxdu(5),这就将原微分方程中的 y y y u u u代替
    • 再对式(5)进行分离变量,得 f ( u ) − u d u = x d x \frac{f(u)-u}{\mathrm{d}u}=\frac{x}{\mathrm{d}x} duf(u)u=dxx(6),两边取倒数,得 d u f ( u ) − u = d x x \frac{\mathrm{d}u}{f(u)-u}=\frac{\mathrm{d}x}{x} f(u)udu=xdx,两边积分: ∫ d u f ( u ) − u \int\frac{\mathrm{d}u}{f(u)-u} f(u)udu= ∫ d x x + C \int{\frac{\mathrm{d}x}{x}}+C xdx+C= ln ⁡ ∣ x ∣ + C \ln|x|+C lnx+C(7)
    • 设(7)左边的积分结果为 Φ ( u ) \Phi(u) Φ(u),代回式(1),得微分方程(0)得通解 Φ ( y x ) \Phi(\frac{y}{x}) Φ(xy)= ln ⁡ ∣ x ∣ + C \ln|x|+C lnx+C

可齐次化的方程

  • 形如 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= f ( a x + b y + c a 1 x + b 1 y + c 1 ) f(\frac{ax+by+c}{a_1x+b_1y+c_1}) f(a1x+b1y+c1ax+by+c)(0)的方程可以通过变换化为齐次方程
  • 分三种情形讨论
    1. c = c 1 = 0 c=c_1=0 c=c1=0,方程(0)就是齐次方程
    2. a a 1 = b b 1 \frac{a}{a_1}=\frac{b}{b_1} a1a=b1b时,设此比值为 λ \lambda λ,则方程(0)可以表示为 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= f ( λ ( a 1 x + b 1 y ) + c a 1 x + b 1 y + c 1 ) f(\frac{\lambda(a_1x+b_1y)+c}{a_1x+b_1y+c_1}) f(a1x+b1y+c1λ(a1x+b1y)+c)(0-1)
      • 引入新变量(函数) u = a 1 x + b 1 y u=a_1x+b_1y u=a1x+b1y(1),对其两边关于 x x x求导: d u d x \frac{\mathrm{d}u}{\mathrm{d}x} dxdu= a 1 + b 1 d y d x a_1+b_1\frac{\mathrm{d}y}{\mathrm{d}x} a1+b1dxdy(2),变形可得 d y d x \frac{\mathrm{d}y}{\mathrm{d}{x}} dxdy= 1 b 1 ( d u d x − a 1 ) \frac{1}{b_1}(\frac{\mathrm{d}u}{\mathrm{d}x}-a_1) b11(dxdua1)(2-1)
      • 将(1)代入(0-1),得 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= f ( λ u + c u + c 1 ) f(\frac{\lambda{u}+c}{u+c_1}) f(u+c1λu+c)(3);比较(2-1),(3),得 1 b 1 ( d u d x − a 1 ) \frac{1}{b_1}(\frac{\mathrm{d}u}{\mathrm{d}x}-a_1) b11(dxdua1)= f ( λ u + c u + c 1 ) f(\frac{\lambda{u}+c}{u+c_1}) f(u+c1λu+c)(4),这是可分离变量( u , x u,x u,x)的方程,问题得到转化
        • 将(4)分离变量,可得 b 1 f ( λ u + c u + c 1 ) + a 1 b_1f(\frac{\lambda{u}+c}{u+c_1})+a_1 b1f(u+c1λu+c)+a1= d u d x \frac{\mathrm{d}u}{\mathrm{d}x} dxdu,即 d x \mathrm{d}x dx= d u b 1 f ( λ u + c u + c 1 ) + a 1 \frac{\mathrm{d}u}{b_1f(\frac{\lambda{u}+c}{u+c_1})+a_1} b1f(u+c1λu+c)+a1du
    3. a a 1 ≠ b b 1 \frac{a}{a_1}\neq{\frac{b}{b_1}} a1a=b1b,且 c , c 1 ≠ 0 c,c_1\neq{0} c,c1=0
      • 借助平面直线方程性质来求解
      • 方程 a x + b y + c = 0 ax+by+c=0 ax+by+c=0; a 1 x + b 1 y + c = 0 a_1x+b_1y+c=0 a1x+b1y+c=0表示平面上的两条直线 l 1 , l 2 l_1,l_2 l1,l2
        • 设两直线交点为 P ( α , β ) P(\alpha,\beta) P(α,β),即 α , β \alpha,\beta α,β满足 a α + b β + c = 0 a\alpha+b\beta+c=0 aα+bβ+c=0, a 1 α + b 1 β + c 1 = 0 a_1\alpha+b_1\beta+c_1=0 a1α+b1β+c1=0
        • 由直线的一般式方程的性质可知, c = 0 c=0 c=0表示直线过原点(0,0)
        • 显然 l 1 , l 2 l_1,l_2 l1,l2的均不过原点,其交点 P P P不可能为原点,即 α , β \alpha,\beta α,β不同时为0( α ≠ 0 \alpha\neq{0} α=0 β ≠ 0 \beta\neq{0} β=0),否则就是情形(1)的直线方程表示
      • 为了解决情形3,参考情形1,可以将坐标系( x O y xOy xOy)平移到 P P P,这就将问题转化为情形(1)
      • 而平移后的坐标记为 X O ′ Y XO'Y XOY其中 O ′ O' O就是 P P P,即新坐标为 X P Y XPY XPY
        • 同一位置的点的新旧坐标转换公式(5): X = x − α X=x-\alpha X=xα; Y = y − β Y=y-\beta Y=yβ
          • 或表示为 x = X + α x=X+\alpha x=X+α; y = Y + β y=Y+\beta y=Y+β
        • l 1 , l 2 l_1,l_2 l1,l2分别表示为 a X + b Y = 0 aX+bY=0 aX+bY=0, a 1 X + b 1 Y = 0 a_1X+b_1Y=0 a1X+b1Y=0;
        • 在新坐标系上,方程(0)表示为 d Y d X \frac{\mathrm{d}Y}{\mathrm{d}X} dXdY= f ( a X + b Y a 1 X + b 1 Y ) f(\frac{aX+bY}{a_1X+b_1Y}) f(a1X+b1YaX+bY)(6)
        • 这就将问题转化为情形(1)(齐次化)
      • 解出(6)的通解后,将(5)代入得到方程(0)的通解

Bernoulli方程

  • 方程 y ′ + p ( x ) y = q ( x ) y n y'+p(x)y=q(x)y^{n} y+p(x)y=q(x)yn(1), ( n ≠ 0 , n ≠ 1 ) (n\neq{0},n\neq{1}) (n=0,n=1)称为Bernoulli方程
伯努利方程求解步骤
  • 原式化为 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)(2),确定出其中的 p ( x ) , q ( x ) , n p(x),q(x),n p(x),q(x),n三要素
  • z = y 1 − n z=y^{1-n} z=y1n(3),列出关于函数 z ( x ) z(x) z(x)的一阶线性微分方程(标准形式): z ′ + ( 1 − n ) p ( x ) z = ( 1 − n ) q ( x ) z'+(1-n)p(x)z=(1-n)q(x) z+(1n)p(x)z=(1n)q(x)(3-0)
  • 代入线性微分方程的通解公式,再代回 y y y,得出 y = y ( x ) y=y(x) y=y(x)(3-1)(得原微分方程的通解
分析和推导
  • n = 0 n=0 n=0 n = 1 n=1 n=1时,方程(1)退化为普通线性微分方程;
  • n ≠ 0 , n ≠ 1 n\neq{0},n\neq{1} n=0,n=1时,此方程虽然不是线性的,但可以通过变量代换,化为线性方程
    • 通常以 y n y^{n} yn除方程(1)两边,得方程(2)
    • 式(2)左端第一项 y − n d y d x y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x} yndxdy(3-2) y 1 − n y^{1-n} y1n的导数 d y 1 − n d x \frac{\mathrm{d}y^{1-n}}{\mathrm{d}x} dxdy1n= ( 1 − n ) y − n ⋅ d y d x (1-n)y^{-n}\cdot\frac{\mathrm{d}y}{\mathrm{d}x} (1n)yndxdy(4)仅差一个常数因子 ( 1 − n ) (1-n) (1n),变形式(4),可得 y − n d y d x y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x} yndxdy= 1 1 − n d y 1 − n d x \frac{1}{1-n}{\frac{\mathrm{d}y^{1-n}}{\mathrm{d}x}} 1n1dxdy1n(4-1)
      • Note: y y y x x x函数, y 1 − n y^{1-n} y1n x x x求导是复合函数求导: ( y 1 − n ) x ′ (y^{1-n})_{x}' (y1n)x= ( 1 − n ) y − n ⋅ y ′ (1-n)y^{-n}\cdot{y'} (1n)yny
    • 因此引入新的因变量 z = y 1 − n z=y^{1-n} z=y1n,即式(3),从而式(4-1)改写为 y − n d y d x y^{-n}\frac{\mathrm{d}y}{\mathrm{d}x} yndxdy= 1 1 − n d z d x \frac{1}{1-n}{\frac{\mathrm{d}z}{\mathrm{d}x}} 1n1dxdz(4-2)
      • 由式(4-2),式(2)写成: 1 1 − n d z d x + p ( x ) z \frac{1}{1-n}\frac{\mathrm{d}z}{\mathrm{d}x}+p(x)z 1n1dxdz+p(x)z= q ( x ) q(x) q(x)
      • d z d x {\frac{\mathrm{d}z}{\mathrm{d}x}} dxdz+ ( 1 − n ) p ( x ) z (1-n)p(x)z (1n)p(x)z= ( 1 − n ) q ( x ) (1-n)q(x) (1n)q(x)(5)
    • 方程(5)显然是一阶线性方程,这就将Bernoulli方程转换为一阶线性方程问题
  • 求一阶常系数微分方程 y ′ − y + 2 x y y'-y+2\frac{x}{y} yy+2yx= 0 0 0(1)
    • 这是一个伯努利方程,将其化为标准形: y ′ − y y'-y yy= − 2 x ⋅ y − 1 -2x\cdot{y^{-1}} 2xy1(2)
    • p ( x ) p(x) p(x)= − 1 -1 1; q ( x ) = − 2 x q(x)=-2x q(x)=2x; n = − 1 n=-1 n=1,即 1 − n = 2 1-n=2 1n=2
    • z = y 1 − n z=y^{1-n} z=y1n= y 2 y^2 y2(3)
    • 由Bernoulli方程的解公式得 z z z的一阶线性微分方程: z ′ + 2 ( − 1 ) z = 2 ( − 2 x ) z'+2(-1)z=2(-2x) z+2(1)z=2(2x),即 z ′ − 2 z = − 4 x z'-2z=-4x z2z=4x(4)
    • 代入一阶线性微分方程解公式
      • g ( x ) g(x) g(x)= ∫ − 2 d x \int{-2}\mathrm{d}x 2dx= − 2 x -2x 2x
      • 从而 z z z= e 2 x ( C + ∫ − 4 x ⋅ e − 2 x d x ) e^{2x}(C+\int{-4x}\cdot{e^{-2x}}\mathrm{d}x) e2x(C+4xe2xdx)= C ( e 2 x + 2 x + 1 ) C(e^{2x}+2x+1) C(e2x+2x+1)(5),由分部积分可以算出该结果
    • 回代(3)到(5),得 y 2 y^2 y2= C ( e 2 x + 2 x + 1 ) C(e^{2x}+2x+1) C(e2x+2x+1),(6)这就是原方程的隐式通解

可降阶的二阶微分方程

  • 二阶及以上的微分方程称为高阶微分方程
  • 某些高阶微分方程可以通过代换法转化为较低阶的方程

类型0

  • y ′ ′ = f ( x ) y''=f(x) y′′=f(x)(0)型的微分方程
    • 方程右端仅为 x x x的函数,作2次积分即可
    • 第1次: y ′ y' y= ∫ f ( x ) d x + C 1 \int{f(x)\mathrm{d}x}+C_1 f(x)dx+C1(1)
    • 第2次: y y y= ∫ ( ∫ f ( x ) d x ) d x + C 1 x + C 2 \int(\int{f(x)}\mathrm{d}x)\mathrm{d}x+C_1x+C_2 (f(x)dx)dx+C1x+C2(2)
  • Note:直接作2次积分: ∫ ( ∫ f ( x ) d x ) d x + C \int{(\int{f(x)\mathrm{d}x})\mathrm{d}x}+C (f(x)dx)dx+C会少一个自由变量
    • 因此,要们分开作2次积分,要么直接用公式(2)

类型1

  • y ′ ′ = f ( x , y ′ ) y''=f(x,y') y′′=f(x,y)(1)型的微分方程
    • 该方程右端不含有未知函数 y y y
    • y ′ = p y'=p y=p( y y y的一阶导数代换为某个函数 p p p)
    • 那么 y ′ ′ = d p d x = p ′ y''=\frac{\mathrm{d}p}{\mathrm{d}x}=p' y′′=dxdp=p;而方程(1)改写成 p ′ = f ( x , p ) p'=f(x,p) p=f(x,p)(2),这就将二阶方程以代换的方法降阶为一阶方程
    • 设(2)的通解为 p = ϕ ( x , C 1 ) p=\phi(x,C_1) p=ϕ(x,C1),即 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= ϕ ( x , C 1 ) \phi(x,C_1) ϕ(x,C1)(3),这是一个关于函数 y y y的一阶微分方程,继续解这个方程
    • 对(3)两边积分,得通解: y = ∫ ϕ ( x , C 1 ) d x + C 2 y=\int\phi(x,C_1)\mathrm{d}x+C_2 y=ϕ(x,C1)dx+C2(4)

类型2

  • y ′ ′ = f ( y , y ′ ) y''=f(y,y') y′′=f(y,y)(1)型的微分方程
    • 方程右端不明显得含有自变量 x x x
    • 同样令 y ′ = p y'=p y=p,利用复合函数求导法,将 y ′ ′ y'' y′′表示为 p p p关于 y y y的导数: d p d y \frac{\mathrm{d}p}{\mathrm{d}y} dydp(微商的角度描述导数)
      • y ′ ′ y'' y′′= d p d x \frac{\mathrm{d}p}{\mathrm{d}x} dxdp= d p d x \frac{\mathrm{d}p}{\mathrm{d}x} dxdp d y d y \frac{\mathrm{d}y}{\mathrm{d}y} dydy= d p d y \frac{\mathrm{d}p}{\mathrm{d}y} dydp d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= d p d y p \frac{\mathrm{d}p}{\mathrm{d}y}p dydpp(2),即公式 y ′ ′ y'' y′′= d p d y p \frac{\mathrm{d}p}{\mathrm{d}y}p dydpp,此公式是核心,配项手法是重点
    • 从而方程(1)化为 p d p d y p\frac{\mathrm{d}p}{\mathrm{d}y} pdydp= f ( y , p ) f(y,p) f(y,p)(3),这是一个一阶微分方程
    • 设(3)的通解为 y ′ = p = ϕ ( y , C 1 ) y'=p=\phi(y,C_1) y=p=ϕ(y,C1)(4),即 d y d x \frac{\mathrm{d}y}{\mathrm{d}x} dxdy= ϕ ( y , C 1 ) \phi(y,C_1) ϕ(y,C1)分离变量得 d y ϕ ( y , C 1 ) = d x \frac{\mathrm{d}y}{\phi(y,C_1)}=\mathrm{d}x ϕ(y,C1)dy=dx(5),
    • 两边积分,得方程的 ∫ d y ϕ ( y , C 1 ) \int{\frac{\mathrm{d}y}{\phi(y,C_1)}} ϕ(y,C1)dy= x + C 2 x+C_2 x+C2,可求得通解
posted @   xuchaoxin1375  阅读(106)  评论(0编辑  收藏  举报  
相关博文:
阅读排行:
· 分享4款.NET开源、免费、实用的商城系统
· 全程不用写代码,我用AI程序员写了一个飞机大战
· MongoDB 8.0这个新功能碉堡了,比商业数据库还牛
· 记一次.NET内存居高不下排查解决与启示
· 白话解读 Dapr 1.15:你的「微服务管家」又秀新绝活了
历史上的今天:
2022-12-17 math_基本导数公式@积分公式@部分推导@数学工具@公式识别@计算推演工具
2022-12-17 算法导论答案网Solutions to Introduction to Algorithms Third Edition_开源免费完整
点击右上角即可分享
微信分享提示