美文网首页简友广场我的大学如切如磋,如琢如磨
二阶变系数线性常微分方程的求解

二阶变系数线性常微分方程的求解

作者: Raow1 | 来源:发表于2021-02-25 20:44 被阅读0次

之前总结过一次常见的常微分方程的解法,不过里面并没有二阶变系数线性常微分方程的求解方法,现在补充如下。

对于f(x)y''+p(x)y'+q(x)y=g(x)形式的微分方程,主要求解步骤为:猜根;刘维尔公式;常数变易。

1. 猜根

此处的猜根是指对方程对应的齐次形式,且忽略初值条件。

  1. f(x)+p(x)+q(x)=0,则存在一特解y=e^x
  2. p(x)=f(x)+q(x),则存在一特解y=e^{-x}
  3. p(x)=-xq(x),则存在一特解y=x
  4. q(x)=0,则存在一特解y=c。(明显也是可降阶的微分方程)

2. 刘维尔公式

得到一个特解后,使用刘维尔公式y_2=y_1\int \frac{1}{y_1^2}e^{-\int P(x)\mathrm dx} \mathrm dx,或者另一形式的刘维尔公式w(x)=Ce^{-\int P(x)\mathrm dx}(以上P(x)=\frac{p(x)}{f(x)}),即可求得另一特解。
于是便得到了对应齐次方程的通解。

个人觉得第二种形式的刘维尔公式算起来更方便。

3. 常数变易

假如通过以上步骤得到齐次方程的通解为Y=C_1y_1(x)+C_2y_2(x),常数变易法令非齐次的通解为y=v_1(x)y_1(x)+v_2(x)y_2(x)。所以有
\begin{cases} & y_1v_1 '+y_2v_2' = 0 \\ & y_1'v_1'+y_2' v_2' = \frac{g(x)}{f(x)} \end{cases}
简写做,
\begin{pmatrix} y_1 & y_2 \\ y_1' & y_2' \end{pmatrix} \begin{pmatrix} v_1' \\ v_2' \end{pmatrix} = \begin{pmatrix} 0 \\ \frac{g(x)}{f(x)} \end{pmatrix}
于是便可以解得v_1',v_2',积分得到v_1,v_2后,代入即可得到非齐次方程的通解。

相关文章

  • 二阶变系数线性常微分方程的求解

    之前总结过一次常见的常微分方程的解法,不过里面并没有二阶变系数线性常微分方程的求解方法,现在补充如下。 对于形式的...

  • 常微分方程

    囊括了一阶线性微分方程,二阶常系数非齐次,二阶变系数,欧拉方程等主要题型。 2016-1-5. 求解柯西问题 显然...

  • 几道微分方程题

    13. 求欧拉方程: 令, 原方程化为 即,。此方程为二阶常系数齐次线性常微分方程。 其特征方程为, 所以, 所以...

  • 第23课 微分方程和exp(At)

    怎么求解一阶方程? 怎么求解一阶导数? 怎么求解常系数线性方程? 将它们转换成线性代数问题的思路是常系数线性方程解...

  • Sturm-Liouville问题定义

      考虑这样一个二阶线性常微分方程  其中被称为谱系数,,系数是区间上的实值函数,且满足自伴边界条件为 其中矩阵是...

  • 考研数学:微分方程考点总结归纳,附各类型微分方程解题对照表

    考研数学微分方程这部分的框架图如图 数一主要侧重于一阶微分方程及二阶常系数线性微分方程的求解,其他类型的方程考试中...

  • 使用Python+SymPy 求解线性方程组

    最近几天在复习线性代数,想用Python求解线性方程组,最开始想用SciPy,但是SciPy求解方程组好像要求系数...

  • 使用scipy库的root和fsolve函数求解方程

    求解非线性方程 使用scipy库的 root, fsolve 函数求解非线性方程。 求解传热方程(热辐射+热传导)...

  • 10-Python 科学计算_scipy 篇

    课程概要:1、非线性方程组2、数值积分3、常微分方程组 1、非线性方程组 2、数值积分(定积分) 3、常微分方程组

  • 常/偏微分方程求解

    常微分方程 1、四阶经典步长Runge-Kutta法求解高阶常微分方程 2、用ode15s函数、ode23s函数和...

网友评论

    本文标题:二阶变系数线性常微分方程的求解

    本文链接:https://www.haomeiwen.com/subject/mbnufltx.html