首页 > 其他分享 >高等数学 7.8常系数非齐次线性微分方程

高等数学 7.8常系数非齐次线性微分方程

时间:2024-10-24 11:03:22浏览次数:1  
标签:mathrm 多项式 cfrac 齐次 7.8 高等数学 omega 特征方程 lambda

目录

二阶常系数非齐次线性微分方程的一般形式是

\[y'' + py' + qy = f(x) \tag{1} \]

其中 \(p, q\) 是常数

由之前的内容可知,求二阶常系数非齐次线性微分方程的通解,归结为求对应的齐次方程

\[y'' + py' + qy = 0 \tag{2} \]

的通解和非齐次方程 \((1)\) 本身的一个特解。

二阶常系数齐次线性微分方程的通解的求法已在 7.7 节得到解决,本节只需讨论求二阶常系数非齐次线性微分方程的一个特解 \(y^*\) 的方法。

本节只介绍当方程 \((1)\) 中的 \(f(x)\) 取两种常见形式时求 \(y^*\) 的方法。这种方法的特点是不用积分就可求出 \(y^*\) ,它叫做 待定系数法 。\(f(x)\) 的两种形式是
(1)\(f(x) = \mathrm{e}^{\lambda x} P_m(x)\) ,其中 \(\lambda\) 是常数,\(P_m(x)\) 是 \(x\) 的一个 \(m\) 次多项式

\[P_m (x) = a_0 x^m + a_1 x^{m - 1} + \cdots + a_{m - 1} x + a_m ; \]

(2)\(f(x) = \mathrm{e}^{\lambda x}[P_l (x) \cos \omega x + Q_n (x) \sin \omega x]\) ,其中 \(\lambda, \omega\) 是常数,\(\omega \neq 0, P_l (x), Q_n (x)\) 分别是 \(x\) 的 \(l\) 次、\(n\) 次多项式,且仅有一个可为零。

一、\(f(x) = \mathrm{e}^{\lambda x} P_m(x)\) 型

方程 \((1)\) 的特解 \(y^*\) 是使 \((1)\) 称为恒等式的函数。怎样的函数能使 \((1)\) 成为恒等式呢?因为 \((1)\) 的右端 \(f(x)\) 是多项式 \(P_m (x)\) 与指数函数 \(\mathrm{e}^{\lambda x}\) 的乘积,而多项式与指数函数乘积的导数仍为多项式与指数函数的乘积,因此,我们推测 \(y^* = R(x) \mathrm{e}^{\lambda x}\) (其中 \(R(x)\) 为某个多项式)可能是方程 \((1)\) 的特解。把 \(y^*, y^{*'}\) 及 \(y^{*''}\) 代入方程 \((1)\) 然后考虑能否选取适当的多项式 \(R(x)\) ,使 \(y^* = R(x) \mathrm{e}^{\lambda x}\) 满足方程 \((1)\) 。为此,将

\[\begin{align*} y^* &= R(x) \mathrm{e}^{\lambda x} \\ y^{*'} &= \mathrm{e}^{\lambda x} [\lambda R(x) + R' (x)] \\ y^{*''} &= \mathrm{e}^{\lambda x}[\lambda^2 R(x) + 2 \lambda R'(x) + R''(x)] \end{align*} \]

代入方程 \((1)\) 并消去 \(\mathrm{e}^{\lambda x}\) ,得

\[R''(x) + (2 \lambda + p)R'(x) + (\lambda^2 + p \lambda + q) R(x) = P_m (x) \tag{3} \]

(i)如果 \(\lambda\) 不是 \((2)\) 的特征方程 \(r^2 + pr + q = 0\) 的根,即 \(r^2 + pr + q \neq 0\) ,由于 \(P_m (x)\) 是一个 \(m\) 次多项式,要使 \((3)\) 的两端相等,那么可令 \(R(x)\) 为另一个 \(m\) 次多项式:

\[R_m (x) = b_0 x^m + b_1 x^{m - 1} + \cdots + b_{m - 1} x + b_m , \]

代入 \((3)\) 式,比较等式两端 \(x\) 同次幂的系数,就得到以 \(b_0, b_1, \cdots, b_m\) 作为未知数的 \(m + 1\) 个方程的联立方程组。从而可以确定这些 \(b_i (i = 0, 1, \cdots, m)\) ,并得到所求的特解 \(y^* = R_m(x) \mathrm{e}^{\lambda x}\) 。

(ii)如果 \(\lambda\) 是特征方程 \(r^2 + pr + q = 0\) 的单根,即 \(r^2 + pr + q = 0\) 但 \(2 \lambda + p \neq 0\) ,要使 \((3)\) 的两端恒等,那么 \(R'(x)\) 必须是 \(m\) 次多项式。此时可令

\[R(x) = x R_m (x) , \]

并可用同样的方法来确定 \(R_m (x)\) 的系数 \(b_i (i = 0, 1, \cdots, m)\) 。

(iii)如果 \(\lambda\) 是特征方程 \(r^2 + pr + q = 0\) 的重根,即 \(r^2 + pr + q = 0\) 且 \(2 \lambda + p = 0\) ,要使 \((3)\) 的两端恒等,那么 \(R'' (x)\) 必须是 \(m\) 次多项式。此时可令

\[R(x) = x^2 R_m(x) , \]

并用同样的方法来确定 \(R_m (x)\) 中的系数。

综上所述,有以下结论:
如果 \(f(x) = \mathrm{e}^{\lambda x} P_m (x)\) ,那么二阶常系数非齐次线性微分方程 \((1)\) 具有形如

\[y^* = x^k R_m (x) \mathrm{e}^{\lambda x} \tag{4} \]

的特解,其中 \(R_m(x)\) 是与 \(P_m (x)\) 同次 (\(m\) 次)的多项式,而 \(k\) 按 \(\lambda\) 不是特征方程的根、是特征方程的单根或是特征方程的重根依次取为 \(0, 1\) 或 \(2\) 。

上述结论可推广到 \(n\) 阶常系数非齐次线性微分方程,但要注意 \((4)\) 式中的 \(k\) 是特征方程含根 \(\lambda\) 的重复次数(即若 \(\lambda\) 不是特征方程的根,则 \(k\) 取为 \(0\);若 \(\lambda\) 是特征方程的 \(s\) 重根,则 \(k\) 取为 \(s\))。

二、\(f(x) = \mathrm{e}^{\lambda x} [P_l (x) \cos \omega x + Q_n (x) \sin \omega x]\) 型

应用欧拉公式

\[\cos \theta = \cfrac{1}{2} (\mathrm{e}^{i \theta} + \mathrm{e}^{-i \theta}), \quad \sin \theta = \cfrac{1}{2 \mathrm{i}} (\mathrm{e}^{i \theta} - \mathrm{e}^{-i \theta}), \]

把 \(f(x)\) 表示成复变指数的形式,有

\[\begin{align*} f(x) &= \mathrm{e}^{\lambda x} \left( P_l \cos \omega x + Q_n \sin \omega x \right) \\ &= \mathrm{e}^{\lambda x} \left( P_l \cfrac{\mathrm{e}^{\omega x \mathrm{i}} + \mathrm{e}^{-\omega x \mathrm{i}}}{2} + Q_n \cfrac{\mathrm{e}^{\omega x \mathrm{i}} - \mathrm{e}^{-\omega x \mathrm{i}}}{2 \mathrm{i}} \right) \\ &= \left( \cfrac{P_l}{2} + \cfrac{Q_n}{2\mathrm{i}} \right) \mathrm{e}^{(\lambda + \omega \mathrm{i})x} + \left( \cfrac{P_l}{2} - \cfrac{Q_n}{2\mathrm{i}} \right) \mathrm{e}^{(\lambda - \omega \mathrm{i})x} \\ &= P(x) \mathrm{e}^{(\lambda + \omega \mathrm{i})x} + \overline{P} (x) \mathrm{e}^{(\lambda - \omega \mathrm{i})x} , \end{align*} \]

其中

\[P(x) = \cfrac{P_l}{2} + \cfrac{Q_n}{2 \mathrm{i}} = \cfrac{P_l}{2} - \cfrac{Q_n}{2} \mathrm{i}, \quad \overline{P}(x) = \cfrac{P_l}{2} - \cfrac{Q_n}{2 \mathrm{i}} = \cfrac{P_l}{2} + \cfrac{Q_n}{2} \mathrm{i} \]

是互成共轭的 \(m\) 次多项式(即它们对应的系数是共轭复数),而 \(m = \max{\{ l, n \}}\) 。

应用之前的结果,对于 \(f(x)\) 中的第一项 \(P(x) \mathrm{e}^{(\lambda + \omega \mathrm{i})x}\) ,可以求出一个 \(m\) 次多项式 \(R_m (x)\) ,使得 \(y_1^* = x^k R_m \mathrm{e}^{(\lambda + \omega \mathrm{i})x}\) 为方程

\[y'' + py' + qy = P(x) \mathrm{e}^{(\lambda + \omega \mathrm{i})x} \]

的特解,其中 \(k\) 按 \(\lambda + \omega \mathrm{i}\) 不是特征方程的根或是特征方程的单根依次取 \(0\) 或 \(1\) 。由于 \(f(x)\) 的第二项 \(\overline{P} (x) \mathrm{e}^{(\lambda - \omega \mathrm{i})x}\) 与第一项 \(P(x) \mathrm{e}^{(\lambda + \omega \mathrm{i})x}\) 成共轭,所以与 \(y_1^*\) 成共轭的函数 \(y_2^* = x^k \overline{R}_m \mathrm{e}^{(\lambda - \omega \mathrm{i})x}\) 必然是方程

\[y'' + py' + qy = \overline{P} (x) \mathrm{e}^{(\lambda - \omega \mathrm{i})x} \]

的特解,这里 \(\overline{R}_m\) 表示与 \(R_m\) 成共轭的 \(m\) 次多项式。于是方程 \((1)\) 具有形如

\[y^* = x^k R_m \mathrm{e}^{(\lambda + \omega \mathrm{i})x} + x^k \overline{R}_m \mathrm{e}^{(\lambda - \omega \mathrm{i})x} \]

的特解。上式可写为

\[\begin{align*} y^* &= x^k \mathrm{e}^{\lambda x} \left( R_m \mathrm{e}^{\omega x \mathrm{i}} + \overline{R}_m \mathrm{e}^{-\omega x \mathrm{i}} \right) \\ &= x^k \mathrm{e}^{\lambda x} [R_m (\cos \omega x + \mathrm{i} \sin \omega x) + \overline{R}_m (\cos \omega x - \mathrm{i} \sin \omega x)] \end{align*} \]

由于括号内的两项是互成共轭的,相加后即无虚部,所以可以写成实函数的形式

\[y_* = x^k \mathrm{e}^{\lambda x} [R_m^{(1)} (x) \cos \omega x + R_m^{(2)} \sin \omega x] \]

综上所述,我们有如下结论:
如果 \(f(x) = \mathrm{e}^{\lambda x}[P_l (x) \cos \omega x + Q_n (x) \sin \omega x]\) ,那么二阶常系数非齐次线性微分方程 \((1)\) 的特解可设为

\[y^* = x^k \mathrm{e}^{\lambda x} [R_m^{(1)} (x) \cos \omega x + R_m^{(2)} \sin \omega x] \tag{5} \]

其中 \(R^{(1)}_m (x), R^{(2)}_m (x)\) 是 \(m\) 次多项式,\(m = \max{\{ l, n \}}\),而 \(k\) 按 \(\lambda + \omega \mathrm{i}\) (或 \(\lambda - \omega \mathrm{i}\))不是特征方程的根或是特征方程的根依次取\(0\)或\(1\) 。

上述结论可以推广到 \(n\) 阶常系数非齐次线性微分方程,但要注意 \((5)\) 式中的 \(k\) 是特征方程中含根 \(\lambda + \omega \mathrm{i}\) (或 \(\lambda - \omega \mathrm{i}\))的重复次数。

标签:mathrm,多项式,cfrac,齐次,7.8,高等数学,omega,特征方程,lambda
From: https://www.cnblogs.com/mowenpan1995/p/18499178/gdsx7-8cxsfqcxxwffc

相关文章

  • 高等数学 7.7常系数齐次线性微分方程
    在二阶齐次线性微分方程\[y''+P(x)y'+Q(x)y=0\tag{1}\]中,如果\(y',y\)的系数\(P(x),Q(x)\)均为常数,即\((1)\)式成为\[y''+py'+qy=0\tag{2}\]其中\(p,q\)是常数,那么称\((2)\)为二阶常系数齐次线性微分方程。如果\(p,q\)不全为常数,就称\((1......
  • 高等数学 7.6高阶线性微分方程
    目录一、线性微分方程的解的结构*二、常数变易法方程\[\cfrac{\mathrm{d}^2y}{\mathrm{d}x^2}+P(x)\cfrac{\mathrm{d}y}{\mathrm{d}x}+Q(x)=f(x)\tag{1}\]叫做二阶线性微分方程。当方程右端\(f(x)\equiv0\)时,方程叫做齐次的;当\(f(x)\not\equiv0\)时,方程叫做非......
  • 高等数学 7.5可降阶的高阶微分方程
    目录一、\(y^{(n)}=f(x)\)型的微分方程二、\(y''=f(x,y')\)型的微分方程三、\(y''=f(y,y')\)型的微分方程一、\(y^{(n)}=f(x)\)型的微分方程微分方程\[y^{(n)}=f(x)\tag{1}\]的右端仅含有自变量\(x\)。容易看出,只要把\(y^{(n-1)}\)作为新的未知函数,那......
  • 高等数学 7.4一阶线性微分方程
    @目录一、线性方程*二、伯努利方程一、线性方程方程\[\cfrac{\mathrm{d}y}{\mathrm{d}x}+P(x)y=Q(x)\tag{1}\]叫做一阶线性微分方程,因为它对于未知函数\(y\)及其导数是一次方程。如果\(Q(x)\equiv0\),那么方程\((1)\)称为齐次的;如果\(Q(x)\not\equiv0\),那么方......
  • 高等数学 7.3 齐次方程
    目录一、齐次方程*二、可化为齐次的方程一、齐次方程如果一阶微分方程可化成\[\cfrac{\mathrm{d}y}{\mathrm{d}x}=\varphi\left(\cfrac{y}{x}\right)\tag{1}\]的形式,那么就称这方程为齐次方程。在齐次方程\[\cfrac{\mathrm{d}y}{\mathrm{d}x}=\varphi\left(\cfrac......
  • 高等数学 7.2 可分离变量的微分方程
    讨论一阶微分方程\[y'=f(x,y)\tag{1}\]的一些解法。一阶微分方程有时也写成如下的对称形式:\[P(x,y)\mathrm{d}x+Q(x,y)\mathrm{d}y=0\tag{2}\]在方程\((2)\)中,变量\(x\)与\(y\)对称,它既可以看作是以\(x\)为自变量\(y\)为因变量的方程\[\cfrac{\mathr......
  • 高等数学 7.1 微分方程的基本概念
    一般地,凡表示未知函数、未知函数的导数与自变量之间的关系的方程,叫做微分方程,有时也简称方程。微分方程中所出现的未知函数的最高阶导数的阶数,叫做微分方程的阶。一般地,\(n\)阶微分方程的形式是\[F(x,y,y',\cdots,y^{(n)})=0\tag{1}\]这里必须指出,在方程\((1)\)中,\(......
  • 高等数学 6.2 定积分在几何学上的应用
    目录一、平面图形的面积1.直角坐标情形2.极坐标情形二、体积1.旋转体体积2.平行截面面积为已知的立体的体积三、平面曲线的弧长一、平面图形的面积1.直角坐标情形我们已经知道,由曲线\(y=f(x)(f(x)\geqslant0)\)及直线\(x=a,x=b(a<b)\)与\(x\)轴所围成的曲边......
  • 高等数学 6.1 定积分的元素法
    在定积分的应用中,经常采用所谓的元素法。为了说明这种方法,先回顾一下曲边梯形的面积问题。设\(f(x)\)在区间\([a,b]\)上连续且\(f(x)\geqslant0\),求以曲线\(y=f(x)\)为曲边、底为\([a,b]\)的曲边梯形的面积\(A\)。把这个面积\(A\)表示为定积分\[A=\int_a^......
  • 高等数学 5.4反常积分
    文章目录一、无穷限的反常积分二、无界函数的反常积分一、无穷限的反常积分设函数f(x)f(x)......