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

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

时间:2024-10-22 18:31:42浏览次数:7  
标签:rx cfrac beta 齐次 7.7 微分方程 alpha 高等数学 mathrm

在二阶齐次线性微分方程

\[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)\) 为二阶变系数齐次线性微分方程

要找微分方程 \((2)\) 的通解,可以先求出它的两个解 \(y_1, y_2\) ,如果它们之比不为常数,即 \(y_1\) 与 \(y_2\) 线性无关,那么 \(y = C_1 y_1 + C_2 y_2\) 就是方程 \((2)\) 的通解。

当 \(r\) 为常数时,指数函数 \(y = \mathrm{e}^{rx}\) 和它的各阶导数都只相差一个常数因子。由于指数函数有这个特点,因此用 \(y = \mathrm{e}^{rx}\) 来尝试,看能否选取适当的常数 \(r\) ,使 \(y = \mathrm{e}^{rx}\) 满足方程 \((2)\) 。

将 \(y = \mathrm{e}^{rx}\) 求导,得到

\[y' = \mathrm{e}^{rx}, \quad y'' = r^2 \mathrm{e}^{rx} \]

把 \(y, y'\) 代入方程 \((2)\) ,得

\[(r^2 + pr + q) \mathrm{e}^{rx} = 0 \]

由于 \(\mathrm{e}^{rx} \neq 0\) ,所以

\[r^2 + pr + q = 0 \tag{3} \]

由此可见,只要 \(r\) 满足代数方程 \((3)\) ,函数 \(y = \mathrm{e}^{rx}\) 就是微分方程 \((2)\) 的解。代数方程 \((3)\) 叫做微分方程 \((2)\) 的特征方程

特征方程 \((3)\) 的两个根 \(r_1, r_2\) 可以用公式

\[r_{1, 2} = \cfrac{-p \pm \sqrt{p^2 - 4q}}{2} \]

求出。它们有三种不同的情形:
(1)当 \(p^2 - 4q > 0\) 时,\(r_1, r_2\) 是两个不相等的实根

\[r_1 = \cfrac{-p + \sqrt{p^2 - 4q}}{2}, \quad r_2 = \cfrac{-p - \sqrt{p^2 - 4q}}{2} \]

(2)当 \(p^2 - 4q = 0\) 时,\(r_1, r_2\) 是两个相等的实根

\[r_1 = r_2 = - \cfrac{p}{2} \]

(3)当 \(p^2 - 4q < 0\) 时,\(r_1, r_2\) 是一对共轭复根

\[r_1 = \alpha + \beta \mathrm{i}, \quad r_2 = \alpha - \beta \mathrm{i} \]

其中

\[\alpha = - \cfrac{p}{2}, \quad \beta = \cfrac{\sqrt{4q - p^2}}{2} \]

相应地,微分方程 \((2)\) 的通解也有三种不同的情_形。
(1)特征方程有两个不相等的实根:\(r_1 \neq r_2\)
由上面的讨论可知,\(y_1 = \mathrm{e}^{r_1x}, y_2 = \mathrm{e}^{r_2x}\) 是微分方程 \((2)\) 的两个解,并且 \(\cfrac{y_2}{y_1} = \cfrac{\mathrm{e}^{r_2x}}{\mathrm{e}^{r_1x}} = \mathrm{e}^{(r_2 - r_1)x}\) 不是常数,因此微分方程 \((2)\) 的通解为

\[y = C_1 \mathrm{e}^{r_1x} + C_2 \mathrm{e}^{r_2x} \]

(2)特征方程有两个相等的实根:\(r_1 = r_2\)
这时,只得到微分方程 \((2)\) 的一个解

\[y_1 = \mathrm{e}^{r_1x} . \]

为了得出微分方程 \((2)\) 的通解,还需求出另一个解 \(y_2\) ,并要求 \(\cfrac{y_2}{y_1}\) 不是常数。设 \(\cfrac{y_2}{y_1} = u(x)\),即 \(y_2 = \mathrm{e}^{r_1x} u(x)\) 。下面来求 \(u(x)\) 。将 \(y_2\) 求导,得

\[\begin{align*} y_2' &= \mathrm{e}^{r_1x}(u' + r_1 u) \\ y_2'' &= \mathrm{e}^{r_1x}(u'' + 2r_1u' + r_1^2 u') \end{align*} \]

将 \(y_2, y_2'\) 和 \(y_2''\) 代入微分方程 \((2)\) ,得

\[\mathrm{e}^{r_1x}[(u'' + 2r_1 u' + r_1^2 u) + p(u' + r_1 u) + qu] = 0 \]

约去 \(\mathrm{e}^{r_1x}\) ,合并同类项,得

\[u'' + (2r_1 + p)u' + (r_1^2 + pr_1 + q)u = 0 \]

由于 \(r_1\) 是特征方程 \((3)\) 的二重根。因此 \(r_1^2 + pr_1 + q = 0\) 且 \(2r_1 + p = 0\),于是得

\[u'' = 0 \]

因为这里只要得到一个不为常数的解,所以不妨取 \(u = x\) ,由此可得到微分方程 \((2)\) 的另一个解

\[y_2 = x \mathrm{e}^{r_1x} \]

从而微分方程 \((2)\) 的通解为

\[y = C_1 \mathrm{e}^{r_1x} + C_2 x \mathrm{e}^{r_1x} \]

\[y = (C_1 + C_2 x) \mathrm{e}^{r_1x} . \]

(3)特征方程有一对共轭复根:\(r_1 = \alpha + \beta \mathrm{i}, r_2 = \alpha - \beta \mathrm{i}(\beta \neq 0)\)
这时,\(y_1 = \mathrm{e}^{(\alpha + \beta \mathrm{i})x}, y_2 = \mathrm{e}^{(\alpha - \beta \mathrm{i})x}\) 是微分方程 \((2)\) 的两个解,但它们是复值函数形式。为了得到实值函数形式的解,先利用欧拉公式 \(\mathrm{e}^{i \theta} = \cos \theta + \mathrm{i} \sin \theta\) 把 \(y_1, y_2\) 改写为

\[\begin{align*} y_1 &= \mathrm{e}^{(\alpha + \beta \mathrm{i})x} = \mathrm{e}^{\alpha x} \cdot \mathrm{e}^{\beta x \mathrm{i}} = \mathrm{e}^{\alpha x} (\cos \beta x + \mathrm{i} \sin \beta x), \\ y_2 &= \mathrm{e}^{(\alpha - \beta \mathrm{i})x} = \mathrm{e}^{\alpha x} \cdot \mathrm{e}^{- \beta x \mathrm{i}} = \mathrm{e}^{\alpha x} (\cos \beta x - \mathrm{i} \sin \beta x) \end{align*} \]

由于复值函数 \(y_1\) 与 \(y_2\) 之间成共轭关系,因此,取它们的和除以2就得到它们的实部,取它们的差除以 \(2 \mathrm{i}\) 就得到它们的虚部。由于方程 \((2)\) 的解符合叠加原理,所以实值函数

\[\overline{y_1} = \cfrac{1}{2} (y_1 + y_2) = \mathrm{e}^{\alpha x} \cos \beta x, \quad \overline{y_2} = \cfrac{1}{2 \mathrm{i}} (y_1 - y_2) = \mathrm{e}^{\alpha x} \sin \beta x \]

还是微分方程 \((2)\) 的解,且 \(\cfrac{\overline{y_1}}{\overline{y_2}} = \cfrac{\mathrm{e}^{\alpha x} \cos \beta x}{\mathrm{e}^{\alpha x} \sin \beta x} = \cot \beta x\) 不是常数,所以微分方程 \((2)\) 的通解为

\[y = \mathrm{e}^{\alpha x} (C_1 \cos \beta x + C_2 \sin \beta x). \]

综上所述,求二阶常系数齐次线性微分方程

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

的通解的步骤如下:

第一步 写出微分方程 \((2)\) 的特征方程

\[r^2 + pr + q = 0 \tag{3} \]

第二步 求出特征方程 \((3)\) 的两个根 \(r_1, r_2\) .

第三步 根据特征方程 \((3)\) 的两个根的不同情形,按照下列表格写出微分方程 \((2)\) 的通解:

特征方程 \(r^2 + pr + q = 0\) 的两个根 \(r_1, r_2\) 微分方程 \(y'' + py' + qy = 0\) 的通解
两个不相等的实根 \(r_1, r_2\) \(y = C_1 \mathrm{e}^{r_1 x} + C_2 \mathrm{e}^{r_2 x}\)
两个相等的实根 \(r_1 = r_2\) \(y = (C_1 + C_2 x)\mathrm{e}^{r_1 x}\)
一对共轭复根 \(r_{1, 2} = \alpha \pm \beta \mathrm{i}\) \(y = \mathrm{e}^{\alpha x} (C_1 \cos \beta x + C_2 \sin \beta x)\)

二阶常系数齐次线性微分方程所用的方法以及方程的通解的形式,可以推广到 \(n\) 阶常系数齐次线性微分方程,简单地叙述如下:
\(n\) 阶常系数齐次线性微分方程的一般形式是

\[y^{(n)} + p_1 y^{(n - 1)} + p_2 y^{(n - 2)} + \cdots + p_{n - 1}y' + p_n y = 0 \tag{4} \]

其中 \(p_1, p_2, \cdots, p_{n - 1}, p_n\) 都是常数。

有时我们用记号 \(\mathrm{D}\) (叫做微分算子)表示对 \(x\) 求导的运算 \(\cfrac{\mathrm{d}}{\mathrm{d}x}\) ,把 \(\cfrac{\mathrm{d}y}{\mathrm{d}x}\) 记作 \(\mathrm{D}y\), 把 \(\cfrac{\mathrm{d}^n y}{\mathrm{d}x^n}\) 记作 \(\mathrm{D}^n y\),并把方程 \((4)\) 记作

\[(\mathrm{D}^n + p_1 \mathrm{D}^{n - 1} + \cdots + p_{n - 1} \mathrm{D} + p_n)y = 0 \tag{5} \]

\[L(\mathrm{D}) = \mathrm{D}^n + p_1 \mathrm{D}^{n - 1} + \cdots + p_{n - 1} \mathrm{D} + p_n \]

\(L(\mathrm{D})\) 叫做微分算子 \(\mathrm{D}\) 的 \(n\) 次多项式。于是方程 \((5)\) 可记作

\[L(\mathrm{D})y = 0 \]

如同讨论二阶常系数齐次线性微分方程那样,令 \(y = \mathrm{e}^{rx}\) 。由于 \(\mathrm{D}\mathrm{e}^{rx} = r \mathrm{e}^{rx}, \cdots, \mathrm{D}^n \mathrm{e}^{rx} = r^n \mathrm{e}^{rx}\) ,故 \(L(\mathrm{D}) \mathrm{e}^{rx} = L(r) \mathrm{e}^{rx}\)。因此把 \(y = \mathrm{e}^{rx}\) 代入方程 \((5)\),得

\[L(r) \mathrm{e}^{rx} = 0 \]

由此可见,如果选取 \(r\) 是 \(n\) 次代数方程 \(L(r) = 0\) 即

\[r^n + p_1 r^{n - 1} + p_2 r^{n - 2} + \cdots + p_{n - 1} r + p_n = 0 \tag{6} \]

的根,那么作出的函数 \(y = \mathrm{e}^{rx}\) 就是方程 \((5)\) 的一个解。
方程 \((6)\) 叫做方程 \((5)\) 的特征方程。

根据特征方程的根,可以写出其对应的微分方程的解如下:

特征方程的根 微分方程通解中的对应项
单实根 \(r\) 给出一项:\(C \mathrm{e}^{rx}\)
一对单复根 \(r_{1, 2} = \alpha \pm \beta \mathrm{i}\) 给出两项:\(\mathrm{e}^{\alpha x} (C_1 \cos \beta x + C_2 \sin \beta x)\)
\(k\) 重实根 \(r\) 给出 \(k\) 项:\(\mathrm{e}^{rx} (C_1 + C_2 x + \cdots + C_k x^{k - 1})\)
一对 \(k\) 重复根 \(r_{1, 2} = \alpha \pm \beta \mathrm{i}\) 给出 \(2k\) 项: \(\mathrm{e}^{\alpha x} [(C_1 + C_2 x + \cdots + C_k x^{k - 1}) \cos \beta x + (D_1 + D_2 x + \cdots + D_k x^{k - 1}) \sin \beta x]\)

从代数学知道,\(n\) 次代数方程有 \(n\) 个根(重根按重复数计算),而特征方程的每一个根都对应着通解中的一项,且每项哥含一个任意常数,这样就得到 \(n\) 阶常系数齐次线性微分方程的通解

\[y = C_1 y_1 + C_2 y_2 + \cdots + C_n y_n . \]

标签:rx,cfrac,beta,齐次,7.7,微分方程,alpha,高等数学,mathrm
From: https://www.cnblogs.com/mowenpan1995/p/18493488/gdsx7-7cxsqcxxwffc

相关文章

  • 高等数学 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)......
  • 高等数学 5.5 反常积分的审敛法 Γ函数
    文章目录无穷限反常积分的审敛法无界函数的反常积分审敛法三、Γ\GammaΓ函数无穷限反常积分的审敛法定理1设函数f(x)f(x)f(x)在区间[a,+∞)[a,+\infty)[a,+∞)上连续,且f(x)⩾0f(x)\geqslant0f(x)⩾0.若函数F(x)=∫axf(t)dtF(x)=\int_a^xf(t)\mathrm{d}t......