常微分二阶线性齐次微分方程的通解推导
对于形如 y′′+py′+qy=0y''+py'+qy=0y′′+py′+qy=0 的常微分二阶线性齐次微分方程,需要先求出两个线性无关的实数域特解 y1y_1y1 和 y2y_2y2 (即 y1y2≠C\frac{y_1}{y_2}\ne Cy2y1=C ),再将两个特解叠加得到通解 y=C1y1+C2y2y=C_1y_1+C_2y_2y=C1y1+C2y2。可以通过求解特征方程来求出微分方程的特解。
设解的形式为 y=erxy=e^{rx}y=erx ,将其带入微分方程可得 (r2+pr+q)erx=0(r^2+pr+q)e^{rx}=0(r2+pr+q)erx=0 。由于erxe^{rx}erx不为0, 可得 r2+pr+q=0r^2+pr+q=0r2+pr+q=0 ,该等式被称为微分方程的特征方程。特征方程是一个一元二次代数方程,其解可由求根公式得到,有三种情况:
(1) Δ=p2−4q>0\Delta = p^2-4q>0Δ=p2−4q>0,存在两个不相等的实根 r1≠r2r_1\ne r_2r1=r2
则两个特解为 y1=er1x,y2=er2xy_1=e^{r_1x}, y_2=e^{r_2x}y1=er1x,y2=er2x,通解为 y=C1er1x+C2er2xy=C_1e^{r_1x}+C_2e^{r_2x}y=C1er1x+C2er2x
(2)Δ=p2−4q=0\Delta = p^2-4q=0Δ=p2−4q=0,存在两个相等的实根 r=r1=r2r=r_1=r_2r=r1=r2
由根的判别公式可知,此时r=−p2r=-\frac{p}{2}r=−2p
此时只有一个特解 y1=erxy_1=e^{rx}y1=erx,因此需要构建出另一个线性无关的特解,设 y2=u(x)erxy_2=u(x)e^{rx}y2=u(x)erx,可知 y2y1=u(x)≠C\frac{y_2}{y_1}=u(x)\ne Cy1y2=u(x)=C,满足线性无关条件
将 y2=u(x)erxy_2=u(x)e^{rx}y2=u(x)erx 带入原微分方程 y′′+py′+qy=0y''+py'+qy=0y′′+py′+qy=0 可得
(u′′+2u′r+ur2+u′p+urp+uq)erx=0(u′′+2u′r+u′p)erx=0(因为 u(r2+rp+q)=0)\begin{aligned} (u''+2u'r+ur^2+u'p+urp+uq)e^{rx}&=0\\ (u''+2u'r+u'p)e^{rx}&=0\quad(\text{因为 }u(r^2+rp+q)=0) \end{aligned} (u′′+2u′r+ur2+u′p+urp+uq)erx(u′′+2u′r+u′p)erx=0=0(因为 u(r2+rp+q)=0)
将 r=−p2r=-\frac{p}{2}r=−2p 带入上式,可得 u′′erx=0⇒u′′=0u''e^{rx}=0\Rightarrow u''=0u′′erx=0⇒u′′=0
在知道二阶导数 u′′(x)u''(x)u′′(x) 等于0后,不妨设 u(x)=xu(x)=xu(x)=x (也可以设 u(x)=kx+bu(x)=kx+bu(x)=kx+b,在叠加特解后可以通过合并常数项得到相同通解),则得到特解 y2=xerxy_2=xe^{rx}y2=xerx,有通解 y=(C1+C2x)erxy=(C_1+C_2x)e^{rx}y=(C1+C2x)erx
(3)Δ=p2−4q<0\Delta = p^2-4q<0Δ=p2−4q<0,存在共轭复根 r1=a+bi,r2=a−bir_1=a+bi, r_2=a-bir1=a+bi,r2=a−bi
即有特解 y1=e(a+bi)x,y2=e(a−bi)xy_1=e^{(a+bi)x}, y_2=e^{(a-bi)x}y1=e(a+bi)x,y2=e(a−bi)x, 根据欧拉公式 ebi=cos(b)+isin(b)e^{bi}=cos(b)+isin(b)ebi=cos(b)+isin(b) ,可将特解化为 y1=eax(cosbx+isinbx),y2=eax(cosbx−isinbx)y_1=e^{ax}(cosbx+isinbx), y_2=e^{ax}(cosbx-isinbx)y1=eax(cosbx+isinbx),y2=eax(cosbx−isinbx), 则有通解 y=eax[D1(cosbx+isinbx)+D2(cosbx−isinbx)]y=e^{ax}[D_1(cosbx+isinbx)+D_2(cosbx-isinbx)]y=eax[D1(cosbx+isinbx)+D2(cosbx−isinbx)]。
以上通解含复数,需要将复数项消去。固定通解中的常数项,可以得到两个不含复数的特解:
令 D1=D2=12D_1=D_2=\frac{1}{2}D1=D2=21, 可得 u1=eaxcosbxu_1=e^{ax}cosbxu1=eaxcosbx
令 D1=12i,D2=−12iD_1=\frac{1}{2i}, D_2=-\frac{1}{2i}D1=2i1,D2=−2i1, 可得 u2=eaxsinbxu_2=e^{ax}sinbxu2=eaxsinbx
则通解为 y=C1u1+C2u2=eax(C1cosbx+C2sinbx)y=C_1u_1+C_2u_2=e^{ax}(C_1cosbx+C_2sinbx)y=C1u1+C2u2=eax(C1cosbx+C2sinbx)
参考资料
常系数二阶线性齐次微分方程的求解 - 知乎 (zhihu.com)
Differential Equations - Complex Roots (lamar.edu)
Differential Equations - Repeated Roots (lamar.edu)
常微分二阶线性齐次微分方程的通解推导相关推荐
- 二阶线性齐次微分方程的解与对应的二阶线性非齐次微分方程的解的关系
以下简称 二阶线性齐次微分方程 为 齐 二阶线性非齐次微分方程 为 非齐 重要的性质.定理(共6条): 证明1:若y1.y2.y3是非齐的解,a.b.c为常数且a+b+c=0,y=a y1+by2+c ...
- 三阶齐次线性方程求通解_非齐次线性方程通解求法------常数变易法.ppt
非齐次线性方程通解求法------常数变易法 高阶线性微分方程 小结.思考题 降阶法与常数变易法 线性微分方程的解的结构 概念的引入 解 受力分析 一.概念的引入 物体自由振动的微分方程 强迫振动的方 ...
- 高数_第5章常微分方程_二阶线性微分方程解的结构
二阶线性微分方程的一般形式是 y'' + p(x)y' + q(x)y = f(x) -------- (1) 这里所谓的线性是指未知函数 y 及其 一阶导y', ...
- (10/300)二阶线性微分方程的解的性质和结构(二)
先占个坑,明天来填....明天能有时间吗????但是我把写博客这个事情必须坚持下去,等闲下来的时候一定要留几篇存稿. +++++++++++++++++++2020年3月20日08点27分++++++ ...
- 三阶齐次线性方程求通解_已知一个三阶常系数线性齐次微分方程的特征根
[简答题]有人说:"电容器带电多电容就大,带电少电容就小,不带电则没有电容."这种说法对吗?为什么? [填空题]思维导图由 英国大脑基金会总裁,被誉为 的英国的 东尼 . 博赞发明 ...
- 二阶龙格库塔公式推导_二阶常系数齐次线性方程通解推导(涉及常数变易法和欧拉公式)...
欧拉恒等式 二阶微分方程明显比一阶难了很多,下面三图详细地对二阶常系数齐次线性方程的通解进行了推导. 有几下几点需要注意: 1.理解思路. 求二阶常系数齐次线性方程的解,一开始是靠猜的,因为以e为底数 ...
- 二阶齐次线性微分方程的通解可以表示成两个线性无关解的线性组合
tex源代码如下: 1 \documentclass[a4paper, 12pt]{article} 2 \usepackage{amsmath,amsfonts,bm} 3 \usepackage{ ...
- 二阶线性微分方程解的结构(齐次与非齐次)+ 常数变易法 | 高阶微分方程(二)
一.线性微分方程的解的结构 1.1 二阶齐次线性方程 y′′+P(x)y′+Q(x)y=0(1)y''+P(x)y'+Q(x)y=0 \tag{1} y′′+P(x)y′+Q(x)y=0(1) 定理1 ...
- (2-1)二阶线性常微分方程相关理论
本文主要涉及内容如下: 1. 二阶齐次线性常微分方程的常点与奇点 2. 二阶齐次线性常微分方程在常点邻域内的解 1. 二阶齐次线性常微分方程的常点与奇点 二阶齐次线性常微分方程的标准形式为: { d ...
- 二阶边值问题的数值解matlab,二阶线性微分方程边值问题的MATLAB求解
2012 年 3 月第 26 卷 第 1 期 阴 山 学 刊 YINSHAN ACADEMIC JOURNAL Mar. 2012 Vo1. 26 No. 1 二阶线性微分方程边值问题的 MATLAB ...
最新文章
- AI一分钟 | 妈呀!连地铁都开始无人驾驶了,飞机还远吗;北京无人驾驶新规出台,终于知道李彦宏该不该被罚了(12月19日)
- R语言ggplot2可视化:使用scale_shape_identity函数显示pch点形状、使用geom_text函数为pch形状添加标(plot characters)
- Linux 下 Redis 安装详解
- mybatis必知必会一
- 怎样用计算机添加文章标题,计算机论文题目怎么定-易指做帮写网
- 如何关闭父窗体?C#
- db2 命令选项解释
- lua笔记之userdata
- 记录一次nginx502/504问题解决过程
- 什么时候不需要回表?
- Python可视化库
- 【教程New】IDM安装使用教程--便捷下载网页视频的好方法
- Android申请权限(相机权限和读写权限)
- K13318 破译信息
- PHP - 垃圾回收机制收集
- 利用python爬取王者荣耀皮肤壁纸
- android源码编译出现No private recovery resources for TARGET_DEVICE解决方法
- 关于Adams驱动函数单位与符号d的问题
- Note For Linux By Jes(13)-认识与分析日志文件
- 唤醒时间过长 Android,关于android:功耗问题之过多唤醒源-wakeups