2.4 一阶隐式微分方程与参数表示

一阶隐式微分方程通常表示: F ( x , y , y ′ ) = 0 F(x,y,y')=0 F(x,y,y)=0需要解决的四种类型:

  • y = f ( x , y ′ ) y=f(x,y') y=f(x,y)
  • x = f ( y , y ′ ) x=f(y,y') x=f(y,y)
  • F ( x , y ′ ) = 0 F(x,y')=0 F(x,y)=0
  • F ( y , y ′ ) = 0 F(y,y')=0 F(y,y)=0

类型一: y = f ( x , d y d x ) y=f(x,\frac{dy}{dx}) y=f(x,dxdy)

  • p = d y d x p=\frac{dy}{dx} p=dxdy于是,变成 y = f ( x , p ) y=f(x,p) y=f(x,p)
  • 不难想到,可以两边对 x x x求偏导,于是变成 p = ∂ f ∂ x + ∂ f ∂ p d p d x p=\frac{\partial f}{\partial x}+\frac{\partial f}{\partial p} \frac{d p}{dx} p=xf+pfdxdp
  • (1)求解,得到 p = φ ( x , c ) p=\varphi(x,c) p=φ(x,c)代入原式,有 y = f ( x , φ ( x , c ) ) y=f(x,\varphi(x,c)) y=f(x,φ(x,c))
    • (2)若求出的解为 x = ψ ( p , c ) x=\psi(p,c) x=ψ(p,c)则得到的参数形式通解为 { x = ψ ( p , c ) y = f ( ψ ( p , c ) , p ) \begin{cases}x=\psi(p,c)\\y=f(\psi(p,c),p)\end{cases} {x=ψ(p,c)y=f(ψ(p,c),p)
    • (3)若得到的通解为 Φ ( x , p , c ) = 0 \Phi(x,p,c)=0 Φ(x,p,c)=0则得到额通解为 { Φ ( x , p , c ) = 0 y = f ( x , p ) \begin{cases}\Phi(x,p,c)=0\\y=f(x,p)\end{cases} {Φ(x,p,c)=0y=f(x,p)
  • (2)与(3)中 p 为 参 数 , c 为 任 意 常 数 p为参数,c为任意常数 p,c

类型二: x = f ( y , d y d x ) x=f(y,\frac{dy}{dx}) x=f(y,dxdy)

  • 依然令 p = d y d x p=\frac{dy}{dx} p=dxdy得到 1 p = ∂ f ∂ x + ∂ f ∂ p d p d x \frac1p=\frac{\partial f}{\partial x}+\frac{\partial f}{\partial p} \frac{d p}{dx} p1=xf+pfdxdp
  • 求得 Φ ( y , p , c ) = 0 \Phi(y,p,c)=0 Φ(y,p,c)=0于是通解为 { x = f ( y , p ) Φ ( y , p , c ) = 0 \begin{cases}x=f(y,p)\\\Phi(y,p,c)=0\end{cases} {x=f(y,p)Φ(y,p,c)=0

一个题

求方程 ( d y d x ) 3 + 2 x d y d x − y = 0 \left(\frac{d y}{d x}\right)^{3}+2 x \frac{d y}{d x}-y=0 (dxdy)3+2xdxdyy=0的解

  • p = d y d x p=\frac{dy}{dx} p=dxdy,有 y = p 3 + 2 x p y=p^{3}+2 x p y=p3+2xp两边对 x x x求导,有 p = 3 p 2 d p d x + 2 x d p d x + 2 p p=3 p^{2} \frac{d p}{d x}+2 x \frac{d p}{d x}+2 p p=3p2dxdp+2xdxdp+2p
  • 经过求解,最终得到 p ≠ 0 p\ne 0 p=0时,有 3 p 4 4 + x p 2 = c \frac{3 p^{4}}{4}+x p^{2}=c 43p4+xp2=c发现可以求解出 x x x,即 x = c − 3 4 p 4 p 2 x=\frac{c-\frac{3}{4} p^{4}}{p^{2}} x=p2c43p4
  • 代入原方程中得到 y = p 3 + 2 ( c − 3 4 p 4 ) p y=p^{3}+\frac{2\left(c-\frac{3}{4} p^{4}\right)}{p} y=p3+p2(c43p4)
  • { x = c p 2 − 3 4 p 2 y = 2 c p − 1 2 p 3 p ≠ 0 \left\{\begin{array}{l} {x=\frac{c}{p^{2}}-\frac{3}{4} p^{2}} \\ {y=\frac{2 c}{p}-\frac{1}{2} p^{3}} \end{array} \quad p \neq 0\right. {x=p2c43p2y=p2c21p3p=0当然 p = 0 p=0 p=0时的情况也要考虑:当 p = 0 p=0 p=0时, y = 0 y=0 y=0也是方程的解

类型三:不显含y的方程 F ( x , y ′ ) = 0 F(x,y')=0 F(x,y)=0

  • 依然令 p = d y d x p=\frac{dy}{dx} p=dxdy(不过这里不是像之前那样对 x x x求导,这里要让 d y = p ⋅ d x dy=p\cdot dx dy=pdx)
  • 这里从几何的角度思考: F ( x , p ) = 0 F(x,p)=0 F(x,p)=0可以看做是 O x p Oxp Oxp平面上的一条曲线,将曲线表示为参数形式 { x = φ ( t ) y = ψ ( t ) \begin{cases}x=\varphi(t)\\y=\psi(t)\end{cases} {x=φ(t)y=ψ(t)
  • 再由上面的 d y = p ⋅ d x dy=p\cdot dx dy=pdx,得到 d y = ψ ( t ) φ ′ ( t ) d t \mathrm{d} y=\psi(t) \varphi^{\prime}(t) \mathrm{d} t dy=ψ(t)φ(t)dt两边积分 y = ∫ ψ ( t ) φ ′ ( t ) d t + c y=\int \psi(t) \varphi^{\prime}(t) \mathrm{d} t+c y=ψ(t)φ(t)dt+c
  • 所以通解为 { x = φ ( t ) y = ∫ ψ ( t ) φ ′ ( t ) d t + c \left\{\begin{array}{l} {x=\varphi(t)} \\ {y=\int \psi(t) \varphi^{\prime}(t) \mathrm{d} t+c} \end{array}\right. {x=φ(t)y=ψ(t)φ(t)dt+c

类型四:不显含x的方程 F ( y , y ′ ) = 0 F(y,y')=0 F(y,y)=0

  • 同类型三做法相似,令 p = d y d x p=\frac{dy}{dx} p=dxdy,引入参数 t t t { y = φ ( t ) p = ψ ( t ) \left\{\begin{array}{l} {y=\varphi(t)} \\ {p=\psi(t)} \end{array}\right. {y=φ(t)p=ψ(t)
  • d y = p   d x dy=p\ dx dy=p dx得到 φ ′ ( t ) d t = ψ ( t ) d x \varphi^{\prime}(t) \mathrm{d} t=\psi(t) \mathrm{d} x φ(t)dt=ψ(t)dx于是!!可以解出 d x = φ ′ ( t ) ψ ( t ) d t \mathrm{d} x=\frac{\varphi^{\prime}(t)}{\psi(t)} \mathrm{d} t dx=ψ(t)φ(t)dt x = ∫ φ ′ ( t ) ψ ( t ) d t + c x=\int \frac{\varphi^{\prime}(t)}{\psi(t)} \mathrm{d} t+c x=ψ(t)φ(t)dt+c
  • 最终解为 { x = ∫ φ ′ ( t ) ψ ( t ) d t + c y = φ ( t ) \left\{\begin{array}{l} {x=\int \frac{\varphi^{\prime}(t)}{\psi(t)} \mathrm{d} t+c} \\ {y=\varphi(t)} \end{array}\right. {x=ψ(t)φ(t)dt+cy=φ(t)
    • 还要注意 p = 0 p=0 p=0的情况,即 F ( y , 0 ) = 0 F(y,0)=0 F(y,0)=0,此时 y = k y=k y=k,则 y = k y=k y=k也是方程的解

今天戴十八给我带了她妈妈做滴鸭脖,结果她自己忍不住吃了一半?!最终还是在把对方喂胖的道路上同归于尽了
在这里插入图片描述

  • 9
    点赞
  • 30
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
MATLAB可以用于求解一微分方程。根据引用,MATLAB可以求解三种类型的一微分方程:显微分方程、线性微分方程和完全微分方程。对于显微分方程,可以直接给出解析解。对于线性微分方程和完全微分方程,可以利用数值方法进行求解。 对于线性微分方程和完全微分方程,可以使用MATLAB中的ode45函数进行求解。这个函数采用常微分方程的初始条件和微分方程的表达作为输入,并返回方程的数值解。ode45函数使用的是龙格-库塔法进行数值求解,可以得到较高的精度。 另外,根据引用,如果已知具体的微分方程表达和边界条件,可以使用MATLAB的ode45函数或其他适用的函数来求解一微分方程。 综上所述,MATLAB提供了丰富的工具和函数来求解一微分方程,可以根据具体的问题选择合适的函数进行求解。<span class="em">1</span><span class="em">2</span><span class="em">3</span> #### 引用[.reference_title] - *1* *2* *3* [MATLAB-常微分方程求解](https://blog.csdn.net/weixin_56691527/article/details/128581996)[target="_blank" data-report-click={"spm":"1018.2226.3001.9630","extra":{"utm_source":"vip_chatgpt_common_search_pc_result","utm_medium":"distribute.pc_search_result.none-task-cask-2~all~insert_cask~default-1-null.142^v93^chatsearchT3_1"}}] [.reference_item style="max-width: 100%"] [ .reference_list ]
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

universe_1207

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值