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

二阶变系数线性微分方程

article 2024/9/21 14:38:19

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/news/156876.html

相关文章：

  • 光伏设计方案中最重要的是什么？
  • 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）——数据类型转换
  • 记一次ThreadPoolTaskExecutor的坑
  • 2023年道路运输企业主要负责人证模拟考试题库及道路运输企业主要负责人理论考试试题
  • IRS辅助的隐蔽通信 （IRS aided covert communication）
  • csapp-linklab之第3阶段“输出学号”实验报告（强弱符号）
  • qt 安装
  • [C/C++]数据结构 堆排序(详细图解)
  • C++ 基础篇
  • 预约按摩小程序有哪些功能特点？
  • autojs-ui悬浮按钮模板
  • 【C语言】存储类型说明符——auto、static、extern、register