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

线性微分方程与非线性微分方程

方程一

d x d t = x \frac{dx}{dt} = x dtdx=x

这是一个一阶线性常微分方程,可以直接分离变量求解。

将变量分离:

d x x = d t \frac{dx}{x} = dt xdx=dt

两边积分:

∫ 1 x d x = ∫ 1 d t ⇒ ln ⁡ ∣ x ∣ = t + C \int \frac{1}{x} \, dx = \int 1 \, dt \Rightarrow \ln |x| = t + C x1dx=1dtlnx=t+C

两边取指数:

∣ x ∣ = e t + C = e C ⋅ e t |x| = e^{t + C} = e^C \cdot e^t x=et+C=eCet
A = ± e C A = \pm e^C A=±eC,即为常数,得通解:

x ( t ) = A e t x(t) = A e^t x(t)=Aet

方程二

d 2 x d t 2 = cos ⁡ x \frac{d^2x}{dt^2} = \cos x dt2d2x=cosx

这个是一个二阶常微分方程,并且由于右侧是 cos ⁡ x \cos x cosx,它是非线性的,因此不能通过常规线性微分方程方法求得解析解。但我们可以采用能量守恒法(乘上 d x d t \frac{dx}{dt} dtdx 积分)将其化为一阶微分方程。

步骤一:将其降阶

设:

v = d x d t ⇒ d 2 x d t 2 = d v d t = d v d x ⋅ d x d t = v ⋅ d v d x v = \frac{dx}{dt} \Rightarrow \frac{d^2x}{dt^2} = \frac{dv}{dt} = \frac{dv}{dx} \cdot \frac{dx}{dt} = v \cdot \frac{dv}{dx} v=dtdxdt2d2x=dtdv=dxdvdtdx=vdxdv

代入原方程:

v ⋅ d v d x = cos ⁡ x v \cdot \frac{dv}{dx} = \cos x vdxdv=cosx

两边同时积分:

∫ v d v = ∫ cos ⁡ x d x ⇒ v 2 2 = sin ⁡ x + C \int v \, dv = \int \cos x \, dx \Rightarrow \frac{v^2}{2} = \sin x + C vdv=cosxdx2v2=sinx+C

即:

( d x d t ) 2 = 2 sin ⁡ x + C 1 ⇒ d x d t = ± 2 sin ⁡ x + C 1 \left(\frac{dx}{dt}\right)^2 = 2 \sin x + C_1 \Rightarrow \frac{dx}{dt} = \pm \sqrt{2\sin x + C_1} (dtdx)2=2sinx+C1dtdx=±2sinx+C1

步骤二:形式上积分得到 t

此时我们得到了:

d x 2 sin ⁡ x + C 1 = ± d t \frac{dx}{\sqrt{2 \sin x + C_1}} = \pm dt 2sinx+C1 dx=±dt

积分两边:

∫ d x 2 sin ⁡ x + C 1 = ± t + C 2 \int \frac{dx}{\sqrt{2 \sin x + C_1}} = \pm t + C_2 2sinx+C1 dx=±t+C2

这个积分通常不能用初等函数表示,它涉及椭圆积分,所以解析解一般只能写成隐式表达式

t = ∫ d x 2 sin ⁡ x + C 1 + C t = \int \frac{dx}{\sqrt{2 \sin x + C_1}} + C t=2sinx+C1 dx+C

相关文章:

  • Windows查看和修改IP,IP互相ping通
  • 智能机器人在物流行业的应用:效率提升与未来展望
  • 【QT】编写第一个 QT 程序 对象树 Qt 编程事项 内存泄露问题
  • 【SystemC初认识】SystemC是什么?有哪些主要组件?如何简单使用?
  • Java高阶程序员学习计划(详细到天,需有一定Java基础)
  • TimeDistill:通过跨架构蒸馏的MLP高效长期时间序列预测
  • LeRobot 项目部署运行逻辑(四)——control_robot.py
  • 使用 Spring Data Redis 实现 Redis 数据存储详解
  • L35.【LeetCode题解】转置矩阵(C语言)
  • 11.Spring Boot 3.1.5 中使用 SpringDoc OpenAPI(替代 Swagger)生成 API 文档
  • 2025.4.28 Vue.js 学习笔记
  • Rancher 2.6.3企业级容器管理平台部署实践
  • 百家号等新媒体私信入口是否可以聚合到企业微信的客服,如何实现
  • E. Unpleasant Strings【Educational Codeforces Round 178 (Rated for Div. 2)】
  • SpringAI实现AI应用-搭建知识库
  • 核心技能:ArcGIS洪水灾害普查、风险评估及淹没制图
  • 【数学建模国奖速成系列】优秀论文绘图复现代码(三)
  • X²+1素数素数
  • 《Python实战进阶》 No46:CPython的GIL与多线程优化
  • 直播美颜SDK是什么?跨平台美颜SDK开发与接入全解析
  • 遍体鳞伤就是击不倒,这是国米老男孩最后的倔强
  • “五一”假期首日国铁郑州局迎大客流,预计发送旅客逾95万人次
  • 医学统计专家童新元逝世,终年61岁
  • 特朗普加征关税冲击波:美国零售、汽车、航空、科技企业纷纷预警业绩波动
  • 笔墨如何“构城”?上海美院城市山水晋京展出
  • 解放日报:中国大模型企业的发展机遇已经到来