n阶常系数非齐次线性微分方程:

$$ y ^ { ( n ) } + p _ { 1 } y ^ { ( n - 1 ) } + p _ { 2 } y ^ { ( n - 2 ) } + .... + p _ { n - 1 } y ^ { \prime } + p _ { n } y = f ( x ) \neq 0 $$

记号:

$$\frac { d } { d x } = D , \frac { d ^ { 2 } } { d x ^ { 2 } } = D ^ { 2 } , ...., \frac { d ^ { n } } { d x ^ { n } } = D ^ { n }$$

则:

$$ y ^ { \prime } = \frac { d y } { d x } = D y , y ^ { \prime \prime } = \frac { d ^ { 2 } y } { d x ^ { 2 } } = D ^ { 2 } y , ...., y ^ { ( n ) } = \frac { d ^ { n } y } { d x ^ { n } } = D ^ { n } y$$

则n阶常系数非齐次线性微分方程可转化为:

$$\left( D ^ { n } + p _ { 1 } D ^ { n - 1 } + ...+ p _ { n _ { - 1 } } D + p _ { n } \right) y = f ( x )$$

记:

$$F ( D ) = D ^ { n } + p _ { 1 } D ^ { n - 1 } + ... + p _ { n - 1 } D + p _ { n }.$$

于是

$F ( D ) y = f ( x )$,从而得特解为:$y ^ { * } = \frac { f ( x ) } { F ( D ) } = \frac { 1 } { F ( D ) } f ( x )$

算子多项式性质

$$\begin{array} { l } { \text { (1)F(D) } \left( \alpha f _ { 1 } ( x ) + \beta f _ { 2 } ( x ) \right) = \alpha F ( D ) f _ { 1 } ( x ) + \beta F ( D ) f _ { 2 } ( x ) } \\ { \text { (2) } \left[ F _ { 1 } ( D ) + F _ { 2 } ( D ) \right] f ( x ) = F _ { 1 } ( D ) f ( x ) + F _ { 2 } ( D ) f ( x ) } \\ { ( 3 ) F ( D ) = F _ { 1 } ( D ) F _ { 2 } ( D ) } \\ { F ( D ) f ( x ) = F _ { 1 } ( D ) \left[ F _ { 2 } ( D ) f ( x ) \right] = F _ { 2 } ( D ) \left[ F _ { 1 } ( D ) f ( x ) \right] } \\ { \text { (4) } \left( \sum\limits _ { k = 0 } ^ { n } F _ { n - k } D ^ { k } \right) [ u ( x ) v ( x ) ] = \sum\limits _ { k = 0 } ^ { n } F _ { n - k } \left[ \sum\limits _ { m = 0 } ^ { k } C _ { k } ^ { m } \left( D ^ { m } u ( x ) \right) \left( D ^ { k - m } v ( x ) \right) \right] } \\ { ( 5 ) F \left( D _ { x } \right) f ( a x + b ) = \left. \left[ F \left( a D _ { u } \right) f ( u ) \right] \right| _ { u = a x + b } } \end{array}$$

算子多项式运算公式

$$\begin{array} { l } { \text { (1) } F ( D ) e ^ { k x } = e ^ { k ( x ) } F ( k ) } \\ { \text { (2) } F \left( D ^ { 2 } \right) \sin a x = \sin a x F \left( - a ^ { 2 } \right) } \\ { F \left( D ^ { 2 } \right) \cos a x = \cos a x F \left( - a ^ { 2 } \right) } \\ { \text { (3) } F ( D ) e ^ { k x } v ( x ) = e ^ { k x } F ( D + k ) v ( x ) } \\ { \text { (4) } F ( D ) x v ( x ) = x F ( D ) v ( x ) + F ^ { \prime } ( D ) v ( x ) } \end{array}$$

逆算子多项式性质

$$\begin{array} { l } { \text { (1) } \frac { 1 } { F ( D ) } F ( D ) f ( x ) = f ( x ) } \\ { \text { (2) } \frac { 1 } { F ( D ) } [ \alpha f ( x ) + \beta g ( x ) ] = \alpha \frac { 1 } { F ( D ) } f ( x ) + \beta \frac { 1 } { F ( D ) } g ( x ) } \\ { \text { (3) } F ( D ) = F _ { 1 } ( D ) F _ { 2 } ( D ) } \\ { \frac { 1 } { F ( D ) } f ( x ) = \frac { 1 } { F _ { 1 } ( D ) } \left[ \frac { 1 } { F _ { 2 } ( D ) } f ( x ) \right] = \frac { 1 } { F _ { 2 } ( D ) } \left[ \frac { 1 } { F _ { 1 } ( D ) } f ( x ) \right] } \end{array}$$

逆算子多项式运算公式

$$\begin{array} { l } { \text { (1) } \frac { 1 } { F ( D ) } e ^ { \lambda x } = \frac { e ^ { \lambda x } } { F ( \lambda ) } ( F ( \lambda ) \neq 0 ) } \\ { \text { (2) } \frac { 1 } { F \left( D ^ { 2 } \right) } \sin \omega x = \frac { \sin \omega x } { F \left( - \omega ^ { 2 } \right) } \left( F \left( - \omega ^ { 2 } \right) \neq 0 \right) } \\   { \text { (3) }\frac { 1 } { F \left( D ^ { 2 } \right) } \cos \omega x = \frac { \cos \omega x } { F \left( - \omega ^ { 2 } \right) } \left( F \left( - \omega ^ { 2 } \right) \neq 0 \right) } \end{array}$$

算子解法及例题

类型一

$f ( x ) = f _ { k } ( x )$

解法:当$F ( 0 ) \neq 0$时$y ^ { * } = \frac { 1 } { F ( D ) } f _ { k } ( x ) = Q _ { k } ( D ) f _ { k } ( x )$

当$F ( 0 ) = 0$时设$\boldsymbol { F } ( \boldsymbol { D } ) = \boldsymbol { D } ^ { r } \overline { \boldsymbol { F } } ( \boldsymbol { D } )$其中$\overline { \boldsymbol { F } } ( 0 ) \neq 0$则

$y ^ { * } = \frac { 1 } { F ( D ) } f _ { k } ( x ) = \frac { 1 } { D ^ { r } } \left[ \frac { 1 } { \overline { F } ( D ) } f _ { k } ( x ) \right] = \frac { 1 } { D ^ { r } } \left[ Q _ { k } ( D ) f _ { k } ( x ) \right]$

例1 $$y ^ { \prime \prime } + 5 y ^ { \prime } + 6 y = x ^ { 2 }$$

解:写成$( D + 2 ) ( D + 3 ) y = x ^ { 2 }$

故对应齐次方程$( D + 2 ) ( D + 3 ) y = 0$通解为$$y _ { 1 } ( x ) = C _ { 1 } e ^ { - 2 x } + C _ { 2 } e ^ { - 3 x }$$

求$y ^ { * } ( x )$

$$( D + 2 ) ( D + 3 ) y ^ { * } = x ^ { 2 }$$

$$\begin{aligned} y ^ { * } ( x ) & = \frac { 1 } { ( D + 2 ) ( D + 3 ) } x ^ { 2 } \\ & = \left( \frac { 1 } { D + 2 } - \frac { 1 } { D + 3 } \right) x ^ { 2 } \end{aligned}$$

$$\begin{array} { l } { = \frac { 1 } { D + 2 } x ^ { 2 } - \frac { 1 } { D + 3 } x ^ { 2 } } \\ { = \frac { 1 } { 2 } \frac { 1 } { 1 + \frac { D } { 2 } } x ^ { 2 } - \frac { 1 } { 3 } \frac { 1 } { 1 + \frac { D } { 3 } } x ^ { 2 } } \\ { = \frac { 1 } { 2 } \left( 1 - \frac { D } { 2 } + \frac { D ^ { 2 } } { 4 } - \cdots \right) x ^ { 2 } } \\ { \quad - \frac { 1 } { 3 } \left( 1 - \frac { D } { 3 } + \frac { D ^ { 2 } } { 9 } - \cdots \right) x ^ { 2 } } \end{array}$$

$$\begin{array} { l } { = \frac { 1 } { 2 } \left( x ^ { 2 } - \frac { 1 } { 2 } \left( x ^ { 2 } \right) ^ { \prime } + \frac { 1 } { 4 } \left( x ^ { 2 } \right) ^ { \prime \prime } \right) } \\ { \quad - \frac { 1 } { 3 } \left( x ^ { 2 } - \frac { 1 } { 3 } \left( x ^ { 2 } \right) ^ { \prime } + \frac { 1 } { 9 } \left( x ^ { 2 } \right) ^ { \prime \prime } \right) } \\ { = \frac { 1 } { 2 } \left( x ^ { 2 } - x + \frac { 1 } { 2 } \right) - \frac { 1 } { 3 } \left( x ^ { 2 } - \frac { 2 x } { 3 } + \frac { 2 } { 9 } \right) } \\ { = \frac { 1 } { 6 } x ^ { 2 } - \frac { 5 } { 18 } x + \frac { 19 } { 108 } } \end{array}$$

通解为:$$\begin{array} { l } { y ( x ) = y _ { 1 } ( x ) + y ^ { * } ( x ) } \\ { = C _ { 1 } e ^ { - 2 x } + C _ { 2 } e ^ { - 3 } + \frac { 1 } { 6 } x ^ { 2 } - \frac { 5 } { 18 } x + \frac { 19 } { 108 } } \end{array}$$

类型二:$f ( x ) = e ^ { 2 x } f _ { k } ( x )$

解法:$y ^ { * } = \frac { 1 } { F ( D ) } \left[ e ^ { \lambda x } f _ { k } ( x ) \right] = e ^ { \lambda x } \frac { 1 } { F ( \lambda + D ) } f _ { k } ( x )$

例2:$$\left( D ^ { 2 } - 6 D + 10 \right) y = e ^ { 3 x }$$

$$y ^ { * } = \frac { 1 } { P ( D ) } e ^ { 3 x } = \frac { 1 } { D ^ { 2 } - 6 D + 10 } e ^ { 3 x } = \frac { 1 } { 3 ^ { 2 } - 6 \cdot 3 + 10 } e ^ { 3 x } = e ^ { 3 x }$$

通解为:

$$y = e ^ { 3 x } \left( C _ { 1 } \cos x + C _ { 2 } \sin x + 1 \right)$$

类型三:$$f ( x ) = f _ { k } ( x ) e ^ { \lambda x } \cos \omega x,f ( x ) = f _ { k } ( x ) e ^ { \lambda x } \sin \omega x$$

解法:

考虑方程

$\boldsymbol { F } ( \boldsymbol { D } ) \boldsymbol { y } = \boldsymbol { e } ^ { ( \lambda + i \omega ) x } \boldsymbol { f } _ { k } ( \boldsymbol { x } )$特解的实部和虚部

例3:求

$\left( \boldsymbol { D } ^ { 2 } - 2 \boldsymbol { D } + 5 \right) \boldsymbol { y } = \boldsymbol { x } e ^ { \boldsymbol { x } } \sin 2 \boldsymbol { x }$

解:考虑$\left( D ^ { 2 } - 2 D + 5 \right) y = e ^ { ( 1 + 2 i ) x } x$特解

$$y ^ { * } = \frac { 1 } { \left( D ^ { 2 } - 2 D + 5 \right) } e ^ { ( 1 + 2 i ) x } x$$

$$\begin{array} { l } { = e ^ { ( 1 + 2 i ) x } \frac { 1 } { ( D + 1 + 2 i ) ^ { 2 } - 2 ( D + 1 + 2 i ) + 5 } x } \\ { = e ^ { ( 1 + 2 i ) x } \left( - \frac { x ^ { 2 } } { 8 } i + \frac { x } { 16 } \right) } \\ { = e ^ { x } \left[ \left( \frac { x ^ { 2 } } { 8 } \sin 2 x + \frac { x } { 16 } \cos 2 x \right) + i \left( - \frac { x ^ { 2 } } { 8 } \cos 2 x + \frac { x } { 16 } \sin 2 x \right) \right] } \end{array}$$

$$y _ { 1 } ^ { * } = e ^ { x } \left( - \frac { x ^ { 2 } } { 8 } \cos 2 x + \frac { x } { 16 } \sin 2 x \right)$$

转载于:https://www.cnblogs.com/Keyon-16/p/10298163.html

微分算子求解微分方程相关推荐

  1. AI攻破高数核心,1秒内精确求解微分方程、不定积分,性能远超Matlab

    栗子 鱼羊 发自 海边边  量子位 报道 | 公众号 QbitAI 大家都知道,AI (神经网络) 连加减法这样的简单算术都做不好: 可现在,AI已经懂得微积分,把魔爪伸向你最爱的高数了. 它不光会求 ...

  2. MATLAB从入门到精通-欧拉法与梯形法求解微分方程(含MATLAB源码)

    前言 以下是我为大家准备的几个精品专栏,喜欢的小伙伴可自行订阅,你的支持就是我不断更新的动力哟! MATLAB-30天带你从入门到精通 MATLAB深入理解高级教程(附源码) tableau可视化数据 ...

  3. matlab求微分方程同届,Matlab学习——求解微分方程(组)

    介绍: 1.在 Matlab 中,用大写字母 D 表示导数,Dy 表示 y 关于自变量的一阶导数,D2y 表示 y 关于自变量的二阶导数,依此类推.函数 dsolve 用来解决常微分方程(组)的求解问 ...

  4. AI攻破高数核心,1秒内精确求解微分方程、不定积分

    大家都知道,AI (神经网络) 连加减法这样的简单算术都做不好: 可现在,AI已经懂得微积分,把魔爪伸向你最爱的高数了. 它不光会求不定积分: 还能解常微分方程: 一阶二阶都可以. 这是Faceboo ...

  5. 栅格单元值的选取方法_计算机求解微分方程的六大数值计算方法

    今天介绍关于计算机求解微分方程的六大数值计算的方法. 1.有限元法 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解 ...

  6. python求解微分方程_python能解微分方程吗

    Sympy Sympy是一个数学符号库(sym代表了symbol,符号),包括了积分,微分方程等各种数学运算方法,为python提供了强大的数学运算支持.对于图像来说,虽然都是做离散的计算,操作最多的 ...

  7. 龙格库塔法和欧拉法求解微分方程的比较

    文章目录 计算机如何理解连续系统的动态特性? 欧拉法求解微分方程 龙格库塔法求解微分方程 MATLAB代码编写和仿真效果 计算机如何理解连续系统的动态特性? 一般连续系统的动态特性可以由一个微分方程, ...

  8. matlab求解微分方程6,牛津大学出版社数学应用软件作业6 用Matlab求解微分方程(组)的解析解和数值解:经典爱情语录大全...

    漳州理工职业学院-酒会礼仪 注:上机作业文件夹以自己的班级姓名学号命名,文件夹包括如下上机报告和Matlab程序. 上机报告模板如下: 佛山科学技术学院 上 机 报 告 课程名称 数学应用软件 上机项 ...

  9. 龙格库塔(runge-kutta,RK)法求解微分方程

    求解微分方程的意思就是,已知导函数,求原函数. 先声明一点,欧拉法.中值法.龙哥库塔法求解微分方程,得出的结果不是表达式,而是一系列离散点. 一.欧拉法递推 问题描述:已知y'=f(x,y),求y(x ...

  10. matlab输入数据作方程,用MATLAB函数编写并求解微分方程

    实验八:用MATLAB 函数编写并求解微分方程 一.实验原理 为了对连续系统进行方针,首先需要建立其数学模型,然后利用计算机求这些数学模型,从而得出数学模型的数值解.由于连续系统是通过微分方程老建模的 ...

最新文章

  1. XamarinEssentials教程首选项Preferences判断项目是否存在
  2. .NET中DLL“没有可放置在工具箱的组件”—FreeTextBox
  3. linux通过html显示图片,如何在Linux上将HTML页面转化成Png图片?
  4. 尼尔机器人技能快捷键_《尼尔机械纪元》连招操作技巧
  5. 你与数据科学家只差这26条python技巧
  6. 猎豹浏览器怎么设置双击关闭网页 双击关闭网页设置方法
  7. Linux学习笔记-动态库的生成
  8. C/C++获取指定网口的IP地址
  9. 【CTR】ESMM:多任务联合学习
  10. Matlab 绘制箱线图
  11. Atitit  循环(loop), 递归(recursion), 遍历(traversal), 迭代(iterate).
  12. 模拟生成随机四位GB2312的汉字验证码(python3.X)
  13. access做mysql前端,ASP+ACCESS留言板制作详细教程
  14. 8086、80286、80386
  15. 电磁仿真软件CST2020 下载链接及安装破·解教程笔记
  16. 最强大脑记忆曲线(10)—— 信息统计
  17. 定时任务框架Quartz
  18. Python数据类型及操作
  19. 8-1 用QPainter绘图(Painting with QPainter)
  20. Python pywin32 学习笔记(一)

热门文章

  1. 利用SMS轻松实现资产管理,SMS2003系列之四
  2. 一文教会你使用Markdown来制作一份简洁明了的求职简历
  3. 苹果Mac视频特效合成神器:Blackmagic Fusion Studio
  4. 区块链之比特币的私钥,公钥和地址是什么?
  5. shell脚本循环嵌套
  6. ARKit何以从同类技术中胜出?
  7. 虚拟化存储解决方案需要确定的三件事
  8. 《CCNA学习指南:Cisco网络设备互连(ICND1)(第4版)》——第1章网络功能
  9. openstack instance bootmgr is missing 问题 修复
  10. bottle框架学习(四)之模版进阶使用