当前位置： 首页 > article >正文

二阶变系数线性微分方程

article 2024/11/13 16:06:31

1、变量替换法

欧拉方程

a_{0}x^{n}\frac{d^{n}y}{d^{n}x}+a_{1}x^{n-1}\frac{d^{n-1}y}{d^{n-1}x}+.......a_{n-1}x\frac{dy}{dx}+a_{n}y=f(x)

a_{0},a_{1},.....a_{n}是常数，f(x)是已知的函数。

二阶欧拉方程

a_{0}x^{2}\frac{d^{2}y}{d^{2}x}+a_{1}x\frac{dy}{dx}+a_{3}y=f(x)  (1)

当x>0时，令x=e^{t},则t=lnx

\frac{dy}{dx}=\frac{dy}{dt}\frac{dt}{dx}=\frac{dy}{dt}\frac{1}{x}

\frac{d^{2}y}{dx^{2}}=\frac{d^{2}y}{dt^{2}}\frac{dt}{dx}=\frac{dy}{dt}\frac{1}{x}+\frac{dy}{dt}\frac{-1}{x^{2}}=(\frac{dy}{dx}-\frac{dy}{dt})\frac{1}{x^{2}}

代入（1）中，

a_{0}\frac{d^{2}y}{dt^{2}}-a_{0}\frac{dy}{dt}+a_{1}\frac{dy}{dt}+a_{2}y=f(e^{t})

a_{0}\frac{d^{2}y}{dt^{2}}+(a_{1}-a_{0})\frac{dy}{dt}+a_{2}y=f(e^{t}).这样就把欧拉方程，化成了二阶常系数非齐次微分方程

当x<0时，令x=-e^{t}, t=ln(-x)=ln|x|

a_{0}\frac{d^{2}y}{dt^{2}}-a_{0}\frac{dy}{dt}+a_{1}\frac{dy}{dt}+a_{2}y=f(-e^{t})

 例题

x^{2}\frac{d^{2}y}{dt^{2}}+x\frac{dy}{dt}=6lnx-\frac{1}{x}

解:令x=e^{t},则t=lnx

代入上面的推导得

\frac{d^{2}y}{dt^{2}}=6t-e^{-t}

\frac{dy}{dt}=3t^{2}+e^{-t}+C_{1}

y=t^{3}-e^{-t}+C_{1}t+C_{2}

所以

y=(lnx)^{3}-\frac{1}{x}+C_{1}lnx+C_{2}

2、降阶法

\frac{d^{2}y}{dt^{2}}+p(x)\frac{dy}{dt}+q(x)y=0 （1）

齐次线性微分方程都是有解的

设(1)有一个已经的非零解y_{1}，

令y=y_{1}u,其中u=u(x)是一个待定函数。

则y^{'}=y^{'}_{1}u+y_{1}u^{'}

y^{''}=y^{''}_{1}u+2y^{'}_{1}u^{'}+y_{1}u^{''}

代入（1）

y_{1}u^{''}+(2y^{'}_{1}+py_{1})u^{'}+(y^{''}_{1}+py^{'}_{1}+qy_{1})u=0

因为y_{1}是解，代入(1)中，公式恒成立，所以

y^{''}_{1}+py^{'}_{1}+qy_{1}=0 成立

所以

y_{1}u^{''}+(2y^{'}_{1}+py_{1})u^{'}=0

转换成一个以u为函数，x自变量的二阶微分方程。

阶数没有阶，再次引入新的变量

令z=u^{'}

y_{1}\frac{dz}{dx}=-(2y^{'}_{1}+py_{1})z

转换成一个以z为函数，x自变量的一阶可分离变量方程。

\frac{dz}{z}=-2\frac{dy_{1}}{y_{1}}-pdx

两边求积分

ln|z|=-2ln|y_{1}|-\int p(x)dx+ln|C_{2}|

ln|z|=ln|\frac{1}{y_{1}^{2}}|-lne^{-\int p(x)dx}+ln|C_{2}|

z=\frac{C_{2}}{y^{2}_{1}}e^{-\int p(x)dx}  z=0 也是解，即C_{2}=0的情形

所以u=C_{1}+C_{2}\int \frac{1}{y^{2}_{1}}e^{-\int p(x)dx}

所以(1)的通解为y=y_{1}[C_{1}+C_{2}\int \frac{1}{y^{2}_{1}}e^{-\int p(x)dx}]

刘维尔公式


http://www.kler.cn/a/156876.html

相关文章：

  • 常用中间件介绍
  • 后端：Aop 面向切面编程
  • C语言第十一周课——函数的调用
  • sqoop import将Oracle数据加载至hive，数据量变少，只能导入一个mapper的数据量
  • FMC 扩展子卡6 路 422,8 组 LVDS,8 路 GPIO
  • Spring资源加载模块，原来XML就这，活该被注解踩在脚下 手写Spring第六篇了
  • 光伏设计方案中最重要的是什么？
  • 2、Linux_远程操作
  • 深入了解Java Duration类，对时间的精细操作
  • 龙迅#LT6911GX是一款高性能HDMI2.1至MIPI或LVDS芯片，支持图像处理，DSC压缩和嵌入式LPDDR4 旋转功能！
  • 【springboot原理篇】Bean的加载方式，面试必看
  • sprintf VS snprintf 函数
  • Realtek蓝牙Android10.0移植结束后的基本测试和常见问题分析
  • Linux halt命令教程：如何安全地关闭你的系统（附详细实例和注意事项）
  • 浅谈基于能耗评价指标的医院智能配电能效管理分析
  • 新版IDEA中，module模块无法被识别，类全部变成咖啡杯无法被识
  • Leetcode 538 把二叉搜索树转换为累加树
  • 管理Android12系统的WLAN热点
  • OpenAI发布一周年，那些声称超过它的模型都怎么样了？
  • 如何知道B站各分区直播数据趋势？
  • MySQL进阶_EXPLAIN重点字段解析
  • 语音芯片的BUSY状态指示功能特征：提升用户体验与系统稳定性的关键
  • JAVA Spring boot Process finished with exit code 0
  • golang channel执行原理与代码分析
  • 基于Langchain的txt文本向量库搭建与检索
  • 菜鸟学习日记（python）——数据类型转换