1 依赖

阅读本文你需要:

  • 学过一点微分方程
  • 脑子

2 写在前面

最近在阅读 Stein 的分析四部曲,Real Analysis 读到 Hilbertian Space 有点读不下去了,破案了原来是我代数没好好学。于是又开始从 Fourier Analysis 读起。这篇文章就取自于 Fourier Analysis 的第一章。Stein 本人是建议,如果想要快速开始傅里叶分析的学习,可以直接跳过第一章,开始第二章,不过我还是读了下去,然后欣喜地发现——第一章非常精彩。它解决了我长久以来的一个疑惑:傅里叶是怎么想出傅里叶级数的?

历史上,傅里叶是研究热传导方程后开始了傅里叶级数的研究。(实际上在他之前已经有几位数学家,例如欧拉、拉格朗日,提出过类似的命题)这在这本 Fourier Analysis 中也有涉及。不过我在这里仅介ban绍yun从波动方程,“变出”傅里叶级数的方法。

3 波动方程

波动方程是形如:

\frac{\partial^2 U}{\partial x^2}=c^2\frac{\partial^2 U}{\partial t^2}.

的一类偏微分方程。不失一般性,我们可以作变量代换让上面的方程变为下面的形式,以简化问题:

\frac{\partial^2 u}{\partial x^2}=\frac{\partial^2 u}{\partial t^2}.

其中 $0\le x\le \pi$ 代表位置,$0\le t$ 代表时间。那么,如何解决这样的问题呢?有两种简单的思路。

3.1. 转化为行波

作变量代换:

\eta=x+t, \xi = x-t.

由链式法则立刻有:

\begin{aligned}
u(x,t)&=v(\eta,\xi)\\&=u((\eta+\xi)/2, (\eta-\xi)/2).\\
v_\eta&=u_x/2+u_y/2.\\
v_{\eta\xi}&=u_{xx}/4-u_{xy}/4+u_{yx}/4-u_{yy}/4\\
&=(u_{xx}-u_{yy})/4=0.
\end{aligned}

于是:

\int\int{v_{\eta\xi}}\mathrm{d}\eta\mathrm{d}\xi=\int\int{0}\mathrm{d}\eta\mathrm{d}\xi=F(\eta)+G(\xi).

将 $x, t$ 代回:

u(x,t)=F(x+t)+G(x-t).

注意到这里的 $F, G$ 是任意的一阶可微函数,并且根据物理意义可以自然地进行奇延拓。这个式子的意义就是,波动方程的解是两个波速相同(为1),方向相反的行波的叠加。

3.2. 转化为驻波

我们接下来考虑这种形式的方程解:

u(x,t)=\psi(x)\phi(t).

这被称为变量分离法。将上式代入我们的波动方程:

\psi''\phi=\phi''\psi.

假设 $\psi, \phi > 0$,则可得:

\frac{\psi''}{\psi}=\frac{\phi''}{\phi}.

由于该式左边为关于 $x$ 的函数,右边为关于 $t$ 的函数,若要该式恒成立:

\exists \lambda, s.t. \frac{\psi''}{\psi}=\frac{\phi''}{\phi}=\lambda.

我们获得一组方程:

\left\{  
             \begin{array}{}  
             \psi''-\lambda\psi=0,\\  
             \phi''-\lambda\phi=0.&  \\  
             \end{array}  
\right.  

由于波函数起源于实际问题,实际上我们有隐含的条件:

u(0,t)=u(\pi,t).

故而上述方程组中我们要求 $\lambda<0$,且 $\sqrt{-\lambda}\in\mathbb{Z}$,

\begin{aligned}
&\text{令} -m^2 \text{ 表示 }\lambda, \text{ 最终的解为:}\\
&\left\{  
             \begin{array}{}  
             \psi(x)=\hat{C}\sin{mx}, &m\in\mathbb{Z}.\\  
             \phi(t)=C_1\sin{mt}+C_2\cos{mt},&  \\  
             \end{array}
\right.\\
\end{aligned}

注意到上面的解对于每一个 $m$ 均为该方程的通解,我们知道微分方程的解可以表示为所有通解的线性组合,于是我们猜测,方程的解有如下的级数形式(先不考虑敛散性):

u(x,t)=\sum_{m=-\infty}^{+\infty}(A_m\sin{mt}+B_m\cos{mt})\sin{mx}.

令 $t=0$,我们猜测

\begin{aligned}
u(x,t)&=\sum_{m=-\infty}^{+\infty}\hat{A_m}\sin{mx}\\
&=F(x)+G(x)=f(x).
\end{aligned}

根据 $f(x)$ 的任意性,和波动方程的物理意义,我们猜测,上述方程可以对任意 $[-\pi, \pi]$ 上的奇函数满足。我们进一步猜测,任意的奇函数可以表示为正弦级数。类似地,猜测对任意偶函数,存在这样的余弦级数展开:

\begin{aligned}
g(x)&=\sum_{m=-\infty}^{+\infty}\tilde{A_m}\cos{mx}.\\
\end{aligned}

于是,由于任意函数都可以表示成奇函数与偶函数的叠加,我们最终猜测,任意函数有如下形式的三角级数展开:

h(x)=\sum_{m=-\infty}^{+\infty}(A_m\cos{mx}+B_m\sin{mx})

如果对上面这个结果感兴趣,不妨去翻一翻 Stein 书第一章最后的热传导方程部分。这里就不详细介绍。

结语

确实,如读者所见,这篇文章(包括 Stein 本人的原文),并没有涉及到硬核的数学证明,整篇文章充满了各种猜测,详细的证明被作者留到了第二章介绍。但我依然认为第一章介绍的这些,是具有启发性、趣味性的、无比精彩的内容。

Offensive77

A little learning is a dangerous thing.

2 Comments

发表评论

您的电子邮箱地址不会被公开。 必填项已用*标注