常系数线性齐次微分方程

                     

贡献者: JierPeter

预备知识 指数函数(复数),线性微分方程的一般理论

   本节讨论的是齐次方程的解法,是常系数线性方程中最简单的情况。非齐次方程的几种解法请参见拉普拉斯变换与常系数线性微分方程文章。

  

未完成:引用所有常系数非齐次方程解法的文章。

1. 实数轴上的复值函数

   首先介绍一个很有用的概念,复值函数。

   复值函数和复变函数是不一样的。复变函数中的 “变” 指变量,因此复变函数是复数域到复数域的映射;但复值函数的自变量还是实数,它只是实数域到复数域的映射而已。

   由于任何复数都可以写成 $a+ \mathrm{i} b$ 的形式,其中 $a, b\in\mathbb{R}$,因此一个复值函数可以看成是两个实值函数的组合。

   设 $f(x), g(x)$ 是某给定区间上的连续实值函数,那么复值函数 $z(x)=f(x)+ \mathrm{i} g(x)$ 的极限被定义为 $\lim\limits_{x\to x_0}z(x)=\lim\limits_{x\to x_0}f(x)+ \mathrm{i} \lim\limits_{x\to x_0}g(x)$。进一步,$z(x)$ 的导函数是 $z'(x)=f'(x)+ \mathrm{i} g'(x)$。

   复数的指数由欧拉公式定义:对于实数 $a, b$,有

\begin{equation} \mathrm{e} ^{a+ \mathrm{i} b}= \mathrm{e} ^a(\cos b+ \mathrm{i} \sin b)~. \end{equation}

   如果 $K$ 是一个常数复数,那么 $ \mathrm{e} ^{Kt}$ 就是一个关于实变量 $t$ 的复值函数。这个函数继承了实值函数 $ \mathrm{e} ^{at}$ 的很多优良性质。

定理 1 

   设 $K, K_i$ 是复常数,$z(t), z_i(t)$ 是区间 $[a, b]$ 上可导的复值函数,那么我们有:

  • $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }(K_1z_1(t)+K_2z_2(t))=K_1z_1'(t)+K_2z_2'(t)$,
  • $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }(z_1(t)z_2(t))=z'_1(t)z_2(t)+z_1(t)z'_2(t)$,
  • $ \mathrm{e} ^{(K_1+K_2)t}= \mathrm{e} ^{K_1t} \mathrm{e} ^{K_2t}$,
  • $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \mathrm{e} ^{Kt}=K \mathrm{e} ^{Kt}$.

   另外,如果记 $\overline{K}$ 是复数 $K$ 的共轭1,那么我们还有

\begin{equation} \mathrm{e} ^{\overline{K}}=\overline{ \mathrm{e} ^K}~. \end{equation}

   复值函数也可以作常系数常微分方程的解。考虑到常微分方程的全体解构成一个线性空间,复值函数可以看成是把这个解空间拓展为一个复数域上的线性空间。如果 $\varphi(x)$、$\phi(x)$ 是某个常微分方程的实值解,那么 $ \mathrm{i} \varphi(x)$、$\varphi(x)+ \mathrm{i} \phi(x)$ 也都是同一个方程的解。

2. 常系数齐次线性微分方程

   由字面意义,可得定义:

定义 1 

   形如

\begin{equation} \frac{\mathrm{d}^n}{ \,\mathrm{d}{t} ^n}x(t)+a_1\frac{\mathrm{d}^{n-1}}{ \,\mathrm{d}{t} ^{n-1}}x(t)+\cdots+a_{n-1}\frac{\mathrm{d}}{ \,\mathrm{d}{t} }x(t)+a_nx(t)=0~ \end{equation}
的方程,称为常系数齐次线性微分方程。其中各 $a_i$ 都是常数。

   接下来我们分情况,详细讨论定义 1 中方程的解法。

无重根的特征方程

   求导运算中最方便的函数是什么?是 $ \mathrm{e} ^{Kt}$,因为其结果总是正比于原来的函数。假设存在复数 $\lambda$,使得 $x(t)= \mathrm{e} ^{\lambda t}$ 是式 3 的解,那么我们得到

\begin{equation} \mathrm{e} ^{\lambda t} \left(\lambda^n+a_1\lambda^{n-1}+\cdots+a_{n-1}\lambda+a_n \right) =0~. \end{equation}

   由于 $ \mathrm{e} ^{\lambda t}$ 处处不为 $0$,式 4 就可以化为

\begin{equation} \lambda^n+a_1\lambda^{n-1}+\cdots+a_{n-1}\lambda+a_n=0~. \end{equation}

   式 5 就被称为式 3 特征方程(characteristic equation)辅助方程(auxiliary equation)

   如果从式 5 中解出 $\lambda$2,那么 $x(t)= \mathrm{e} ^{\lambda t}$ 就是式 3 的一个解。根据预备知识中的定理 2 ,如果式 5 有 $n$ 个互不相同的解 $\lambda_i$,那么 $\{ \mathrm{e} ^{\lambda_i t}\}$ 就构成了式 3 的一组基解。这样,我们就把解常系数齐次线性微分方程化为解代数方程了。

   事实上,由于我们常把高阶求导写成低阶求导的乘积形式,比如 $\frac{\mathrm{d}^3}{ \,\mathrm{d}{t} ^3}=\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }\cdot\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }\cdot\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }$,,因此式 5 形式上相当于将 $\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }$ 替换为 $\lambda$ 并消去 $x(t)$ 后的结果。之所以我们可以将 $\lambda$ 视为一个数字去解,也是因为 $ \mathrm{e} ^{\lambda t}$ 在微分算子作用下的特殊性。

例 1 

   考虑方程

\begin{equation} \frac{\mathrm{d}^2}{ \,\mathrm{d}{t} ^2}x(t)-\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }x(t)-2x(t)=0~. \end{equation}
其特征方程为
\begin{equation} (\lambda-2)(\lambda+1)=0~, \end{equation}

   因此其通解为

\begin{equation} x(t)=A \mathrm{e} ^{2t}+B \mathrm{e} ^{-t}~. \end{equation}
其中 $A, B$ 是常数。

有重根的情况

   有的时候,特征方程存在重根,导致无法解出 $n$ 个不同的 $\lambda$,这时候该怎么办呢?我们先观察以下例子。

例 2 

   考虑方程

\begin{equation} \frac{\mathrm{d}^2}{ \,\mathrm{d}{t} ^2}x(t)-2\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }x(t)+x(t)=0~. \end{equation}

   其特征方程为

\begin{equation} (\lambda-1)^2=0~, \end{equation}
只能解出一个二重根 $\lambda=1$。按照之前的讨论,我们只能给出一个基解 $ \mathrm{e} ^{t}$,还有一个解怎么办呢?

   尝试把 $t \mathrm{e} ^{t}$ 代进式 9 试试,我们发现这也是一个解。这样,我们就有了两个线性无关解,也就解出了式 9

   为什么例 2 中给出的 $t \mathrm{e} ^{t}$ 也构成一个解呢?

   我们首先考虑 $k$ 重根里最简单的情况,即这个根为 $\lambda_0=0$。这意味着特征方程里有一个因子 $\lambda^k$,进而意味着式 3 里有一个因子 $\frac{\mathrm{d}^k}{ \,\mathrm{d}{t} ^k}= \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \right) ^k$,故原方程的形式应为

\begin{equation} \left(\frac{\mathrm{d}^n}{ \,\mathrm{d}{t} ^n}+a_1\frac{\mathrm{d}^{n-1}}{ \,\mathrm{d}{t} ^{n-1}}+\cdots+a_{n-k}\frac{\mathrm{d}^k}{ \,\mathrm{d}{t} ^k} \right) x(t)=0~. \end{equation}
或更直观地,
\begin{equation} \begin{aligned} & \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \right) ^k \left(\frac{\mathrm{d}^{n-k}}{ \,\mathrm{d}{t} ^{n-k}}+a_1\frac{\mathrm{d}^{n-k-1}}{ \,\mathrm{d}{t} ^{n-k-1}}+\cdots+a_{n-k} \right) x(t)\\=&0\\ =& \left(\frac{\mathrm{d}^{n-k}}{ \,\mathrm{d}{t} ^{n-k}}+a_1\frac{\mathrm{d}^{n-k-1}}{ \,\mathrm{d}{t} ^{n-k-1}}+\cdots+a_{n-k} \right) \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \right) ^kx(t)~, \end{aligned} \end{equation}
相当于我们有两个方程
\begin{equation} \left\{\begin{aligned} & \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \right) ^kx(t)=0\\ & \left(\frac{\mathrm{d}^{n-k}}{ \,\mathrm{d}{t} ^{n-k}}+a_1\frac{\mathrm{d}^{n-k-1}}{ \,\mathrm{d}{t} ^{n-k-1}}+\cdots+a_{n-k} \right) x(t)=0~. \end{aligned}\right. \end{equation}

   事实上,如果将特征方程拆分成 $(\lambda-\lambda_1)^{k_1}(\lambda-\lambda_2)^{k_2}\cdots$ 的形式,原方程实际上就是在解每一个 $ \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_i \right) ^{k_i}x(t)=0$。因此,我们只需要考虑形如 $ \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) ^kx(t)=0$ 的方程即可。

   在最简单情形下,我们要考虑的方程是

\begin{equation} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \right) ^kx(t)=0~, \end{equation}
易验证它的一个基本解组是 $\{1, t, t^2, \cdots, t^{k-1}\}$。这是因为形如 $t^m$ 的多项式函数,能够在 $ \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} } \right) ^k$ 算子作用下变成 $0$,其中 $k>m$。

   那么,对于一般的情况,即

\begin{equation} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) ^kx(t)=0~, \end{equation}
我们要考虑的就是什么样的 $x(t)$ 能被 $ \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) ^k$ 变成 $0$。

   答案自然是形如 $t^m \mathrm{e} ^{\lambda_0 t}$ 的函数。

   首先观察到

\begin{equation} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) t^m \mathrm{e} ^{\lambda_0 t}=mx^{m-1} \mathrm{e} ^{\lambda_0 t}~, \end{equation}
因此每个 $ \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) $ 算子就像是在单独给 $t^m$ 降阶。

   于是

\begin{equation} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) ^mt^m \mathrm{e} ^{\lambda_0 t}=m! \mathrm{e} ^{\lambda_0 t}~, \end{equation}
再作用一次,就得到
\begin{equation} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_0 \right) ^{m+1}t^m \mathrm{e} ^{\lambda_0 t}=0~. \end{equation}
因此,对于任意 $m< k$,$t^m \mathrm{e} ^{\lambda_0 t}$ 都是式 15 的解。

   这样,我们就得到齐次线性微分方程的完整通解了:

定理 2 齐次线性微分方程的通解

   将齐次线性微分方程式 3 分解为如下形式:

\begin{equation} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_1 \right) ^{k_1} \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_2 \right) ^{k_2}\cdots \left(\frac{ \,\mathrm{d}{}} { \,\mathrm{d}{t} }-\lambda_r \right) ^{k_r}x(t)=0~. \end{equation}
换句话说,就是式 3 的特征方程有根 $\lambda_1, \lambda_2, \cdots, \lambda_r$,且 $\lambda_i$ 的重数为 $k_i$。

   那么式 3 或者说式 19 的基本解组为

\begin{equation} \{ \mathrm{e} ^{\lambda_i t}, t \mathrm{e} ^{\lambda_i t}, \cdots, t^{k_i-1} \mathrm{e} ^{\lambda_i t} \}~, \end{equation}
其中 $i\in\{1, 2, 3, \cdots, r\}$。


1. ^ 即如果存在实数 $a, b$ 使 $K=a+ \mathrm{i} b$,那么 $\overline{K}=a- \mathrm{i} b$。
2. ^ $\lambda$ 可以是复数,这就是为什么我们要先介绍复值函数的概念。


致读者: 小时百科一直以来坚持所有内容免费无广告,这导致我们处于严重的亏损状态。 长此以往很可能会最终导致我们不得不选择大量广告以及内容付费等。 因此,我们请求广大读者热心打赏 ,使网站得以健康发展。 如果看到这条信息的每位读者能慷慨打赏 20 元,我们一周就能脱离亏损, 并在接下来的一年里向所有读者继续免费提供优质内容。 但遗憾的是只有不到 1% 的读者愿意捐款, 他们的付出帮助了 99% 的读者免费获取知识, 我们在此表示感谢。

                     

友情链接: 超理论坛 | ©小时科技 保留一切权利