二阶常系数非齐次线性微分方程：类型二

最新推荐文章于 2024-08-01 08:17:37 发布

Richard888888888

最新推荐文章于 2024-08-01 08:17:37 发布

阅读量449

点赞数 9

文章标签：机器学习线性代数人工智能

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/scdn1172/article/details/134658719

版权

$f(x)=p_{n}(x)e^{\alpha x}cos\beta x$ 或 $f(x)=p_{n}(x)e^{\alpha x}sin\beta x$

$\frac{d^{2}}{dx^{2}}+p\frac{dy}{dx}+qy=p_{n}(x)e^{\alpha x}cos\beta x$ (1)

$\frac{d^{2}}{dx^{2}}+p\frac{dy}{dx}+qy=p_{n}(x)e^{\alpha x}sin\beta x$ (2)

欧拉公式

$e^{i\beta x}=cos\beta x+isin\beta x$

$p_{n}(x)e^{(\alpha +i\beta) x}=p_{n}(x)cos\beta x+p_{n}(x)isin\beta x$

$\frac{d^{2}}{dx^{2}}+p\frac{dy}{dx}+qy=p_{n}(x)e^{(\alpha+i\beta ) x}$ (3)

设(3)有特解 $\widetilde{y}=\widetilde{y_{1}}+i\widetilde{y_{2}}$

$\widetilde{y_{1}}$ 是(1)的特解

$\widetilde{y_{2}}$ 是（2）的特解

例：

$\frac{d^{2}}{dx^{2}}-y=sinx$

解 $r^{2}-1=0 ,r=\pm 1$

所以其通解为

$Y=C_{1}e^{X}+C_{2}e^{-x}$

$n=0,\alpha =0,\beta =1$

$\frac{d^{2}}{dx^{2}}-y=e^{ix}$

这是属于上一篇提到的第一种情况

其特解

$\widetilde{y}=ae^{ix}$

$\frac{d\widetilde{y}}{dx}=iae^{ix}$

$\frac{d\widetilde{^{2}y}}{dx^{2}}=-ae^{ix}$

即 $-ae^{ix}-ae^{ix}=e^{ix}$

简化

-2a=1

a=-1/2

$\widetilde{y}=-\frac{1}{2}e^{ix}=-\frac{1}{2}(cosx+isinx)$

原方程特解为

$\widetilde{y}=-\frac{1}{2}sinx$

原方程通解为

$Y=Y+\widetilde{y}$

解法二

将 $P_{n}(x)e^{\alpha x}cos\beta x$ 或 $P_{n}(x)e^{\alpha x}sin\beta x$

改写为 $f(x)=[P^{1}_{n}(x)cos\beta x+P^{2}_{l}sin\beta x]e^{\alpha x}$

设 $\alpha +i\beta$ 是特征方程的k重根（k=0或1）

则方程特解形如：

$\widetilde{y}=x^{k}[R^{1}_{m(x)}cos\beta x+R^{2}_{m}(x)sin\beta x]e^{\alpha x}$

其中 m= max{n，l} ,指多项式的幂。

用解法二重求例题

$\frac{d^{2}}{dx^{2}}-y=sinx$

解

解 $r^{2}-1=0 ,r=\pm 1$

所以其通解为

$Y=C_{1}e^{X}+C_{2}e^{-x}$

$m=0,\alpha =0,\beta =1$

i不是特征方程的根，k=0

设 $\widetilde{y}=acosx+bsinx$ 为特解

$\frac{d\widetilde{y}}{dx}=-asinx+bcosx$

$\frac{d^{2}\widetilde{y}}{dx^{2}}=-acosx-bsinx$

所以代入原式

-acosx-bsinx-acosx-bsinx=sinx

-2acosx-(2b+1)sinx=0

因为 $\frac{sinx}{cosx}\neq$ 常数

所以线性无关。

所以 -2a=-(2b+1)=0

a=0,b=-1/2

例

$\frac{d^{2}}{dx^{2}}-y=e^{x}cos2x+sinx$

可以利用叠加原理，分别求特解。

1、先求齐次线性微分方程的通解

$\frac{d^{2}}{dx^{2}}-y=0$ 通解为 $Y=C_{1}e^{x}+C_{2}e^{-x}$

将右边的自由项，拆解成两部分。

$\frac{d^{2}}{dx^{2}}-y=e^{x}cos2x$ 的特解 $\widetilde{y_{1}}=-\frac{1}{8}e^{x}(cos2x-\frac{1}{8}sin2x)$

$\frac{d^{2}}{dx^{2}}-y=sinx$ 的特解 $\widetilde{y_{1}}=-\frac{1}{2} sinx$

所以所求方程为 $y=Y+\widetilde{y_{1}}+\widetilde{y_{2}}$

高阶方程的求解

例

$y^{'''}+3y^{''}+3y^{'}+y=e^{x}$

$r^{3}+3r^{2}+3r+1=0$

$(r+1)^{3}=0$ , $r_{1}=r_{2}=r_{3}=-1$

$y^{'''}+3y^{''}+3y^{'}+y=0$ 的通解 $Y=(C_{1}+C_{2}+C_{3}x^{2})e^{-x}$

$\alpha =1$ 不是特征方程的根。

可设特解为 $\widetilde{y}=ae^{-x}$

$\widetilde{y^{'}}=\widetilde{y^{''}}=\widetilde{y^{'''}}=ae^{-x}$ 代入原方程

$a+3a+3a+a=1$

a=1/8

Richard888888888

关注

9
点赞
踩
8

收藏

觉得还不错? 一键收藏
0
评论
二阶常系数非齐次线性微分方程：类型二

其中 m= max{n，l} ,指多项式的幂。是特征方程的k重根（k=0或1）所以 -2a=-(2b+1)=0。这是属于上一篇提到的第一种情况。可以利用叠加原理，分别求特解。1、先求齐次线性微分方程的通解。将右边的自由项，拆解成两部分。i不是特征方程的根，k=0。
复制链接

扫一扫

Richard888888888 CSDN认证博客专家 CSDN认证企业博客

码龄8年

75: 原创

8万+: 周排名

3万+: 总排名

15万+: 访问

: 等级

1406: 积分

53: 粉丝

100: 获赞

13: 评论

259: 收藏

私信

关注

热门文章

分类专栏

最新评论

二阶变系数线性微分方程
urkawaii: 这里并没有具体的写出通解Y的求解方式，而且举例非常不恰当，下面这个知乎链接作为通解的求解方式https://zhuanlan.zhihu.com/p/352960264 结合着看可以消灭99%的问题，如果猜根确实很特殊，只能通过该文的级数法试着去求解了。
电路中的叠加定理
Richard888888888: 电流源开路，就是最基本的并联电路了
电路中的叠加定理
m0_56344318: 你好，这个图c是怎么看电阻串并联关系的？
用keil5，STM32FO13,调用数组时，初始化为0，则会出现乱码现象
SunriseF: 一样的问题，初始化不赋值和赋值为0都会造成这个问题，最后只有赋值1了
互感确定同名端的方法
句霜: 非常好的教程，不过前三图片要是放到对应文字前阅读观感可能更好一点

最新文章

目录

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。