当前位置: 首页 > news >正文

【高等数学】第七章 微分方程——第四节 一阶线性微分方程

上一节【高等数学】第七章 微分方程——第三节 齐次方程
总目录【高等数学】 目录

文章目录

  • 1. 线性方程
  • 2. 伯努利方程

1. 线性方程

  • 方程
    dydx+P(x)y=Q(x)\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = Q(x)dxdy+P(x)y=Q(x)
    叫做一阶线性微分方程,因为它对于未知函数yyy及其导数是一次方程. 如果Q(x)≡0Q(x) \equiv 0Q(x)0,那么方程是齐次
    如果Q(x)≢0Q(x) \not\equiv 0Q(x)0,那么方程是非齐次的.
  • 解法
    Q(x)=0Q(x)=0Q(x)=0,写出对应的齐次线性方程
    dydx+P(x)y=0\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = 0dxdy+P(x)y=0
    分离变量后,两边积分得
    ln⁡∣y∣=−∫P(x)dx+C1\ln \vert y \vert = -\int P(x) \mathrm{d}x + C_1lny=P(x)dx+C1
    y=Ce−∫P(x)dxy = C\mathrm{e}^{-\int P(x) \mathrm{d}x}y=CeP(x)dx
    使用常数变易法,求非齐次线性方程的通解
    u=Cu=Cu=Cy=ue−∫P(x)dxy = u\mathrm{e}^{-\int P(x) \mathrm{d}x}y=ueP(x)dx
    dydx=u′e−∫P(x)dx−uP(x)e−∫P(x)dx\displaystyle\dfrac{\mathrm{d}y}{\mathrm{d}x} = u'\mathrm{e}^{-\int P(x) \mathrm{d}x} - uP(x)\mathrm{e}^{-\int P(x) \mathrm{d}x}dxdy=ueP(x)dxuP(x)eP(x)dx
    代入dydx+P(x)y=Q(x)\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y = Q(x)dxdy+P(x)y=Q(x)
    u′e−∫P(x)dx−uP(x)e−∫P(x)dx+P(x)ue−∫P(x)dx=Q(x)u'\mathrm{e}^{-\int P(x) \mathrm{d}x} - uP(x)\mathrm{e}^{-\int P(x) \mathrm{d}x} + P(x)u\mathrm{e}^{-\int P(x) \mathrm{d}x} = Q(x)ueP(x)dxuP(x)eP(x)dx+P(x)ueP(x)dx=Q(x)
    u′=Q(x)e−∫P(x)dxu'=Q(x)\mathrm{e}^{-\int P(x) \mathrm{d}x}u=Q(x)eP(x)dx
    u=∫Q(x)e−∫P(x)dxdx+C\displaystyle u=\int Q(x)\mathrm{e}^{-\int P(x) \mathrm{d}x} \mathrm{d}x+ Cu=Q(x)eP(x)dxdx+C
    因此非齐次线性方程的通解
    y=e−∫P(x)dx(∫Q(x)e∫P(x)dxdx+C)=Ce−∫P(x)dx+e−∫P(x)dx∫Q(x)e∫P(x)dxdx\begin{aligned}y &= \mathrm{e}^{-\int P(x) \mathrm{d}x} \left( \int Q(x) \mathrm{e}^{\int P(x) \mathrm{d}x} \mathrm{d}x + C \right)\\ &= C\mathrm{e}^{-\int P(x) \mathrm{d}x} + \mathrm{e}^{-\int P(x) \mathrm{d}x} \int Q(x) \mathrm{e}^{\int P(x) \mathrm{d}x} \mathrm{d}x\end{aligned}y=eP(x)dx(Q(x)eP(x)dxdx+C)=CeP(x)dx+eP(x)dxQ(x)eP(x)dxdx
    右边是齐次线性方程的通解,左边是非齐次线性方程的一个特解

2. 伯努利方程

  • 方程
    dydx+P(x)y=Q(x)yn(n≠0,1)\dfrac{\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)
    叫做伯努利(Bernoulli)方程.
    n=0n = 0n=0n=1n = 1n=1时,这是线性微分方程.
    n≠0n \neq 0n=0n≠1n \neq 1n=1时,这方程不是线性的,但是通过变量的代换,便可把它化为线性的.
  • 解法
    两边同除以yny^nyn
    y−ndydx+P(x)y1−n=Q(x).y^{-n}\dfrac{\mathrm{d}y}{\mathrm{d}x} + P(x)y^{1 - n} = Q(x).yndxdy+P(x)y1n=Q(x).
    z=y1−nz=y^{1-n}z=y1ndzdx=(1−n)y−ndydx\dfrac{\mathrm{d}z}{\mathrm{d}x}=(1-n)y^{-n}\dfrac{\mathrm{d}y}{\mathrm{d}x}dxdz=(1n)yndxdy,则可得线性方程
    dzdx+(1−n)P(x)z=(1−n)Q(x)\dfrac{\mathrm{d}z}{\mathrm{d}x} + (1-n)P(x)z = (1-n)Q(x)dxdz+(1n)P(x)z=(1n)Q(x)
    求得通解后,再将z=y1−nz=y^{1-n}z=y1n回代即可得伯努利方程通解

下一节
总目录【高等数学】 目录

http://www.dtcms.com/a/307056.html

相关文章:

  • LNN+XGBoost:优化多层供应链订购:缓解牛鞭效应
  • C++STL系列之bitset
  • Git——分布式版本控制系统
  • #C语言——学习攻略:深挖指针路线(四)--字符指针变量,数组指针变量,二维数组传参的本质,函数指针变量,函数指针数组
  • ConvertX:自托管的在线文件转换器,支持1000+种格式!
  • Linux系统编程Day1-- 免费云服务器获取以及登录操作
  • CH347使用笔记:CH347作为FPGA下载器的几种方式
  • Maven 配置阿里云镜像加速
  • huggingface是什么?2025-07-30
  • Mac 上配置jdk 环境变量
  • 2. Agent与 React流程
  • 【LY88】双系统指南及避坑
  • Python 的 match-case
  • 从映射到共生:元宇宙、物联网与AI的智能融合生态图谱
  • (LeetCode 面试经典 150 题) 141. 环形链表(快慢指针)
  • HPCtoolkit的下载使用
  • Oracle11g数据库迁移达梦8数据库方案
  • Python序列化和反序列化
  • 如何用Docker部署ROS2
  • (C++)C++类和类的方法(基础教程)(与Python类的区别)
  • c++之基础B之sort排序(第三个参数没有)(第二课)
  • Fiddler中文教程 从入门到进阶的网络抓包与接口调试实战指南
  • Python Pandas.merge_asof函数解析与实战教程
  • VUE前端
  • [Agent开发平台] API网关 | 业务领域 | DTO格式 | 分页令牌
  • React 服务端渲染(SSR)详解
  • 使用 cron 配合 Docker 实现定时任务
  • 神经网络的并行计算与加速技术
  • 模型相关类代码回顾理解 | BatchNorm2d\fc.in_features\nn.Linear\torchsummary
  • Haproxy 七层代理深度解析