利用复傅里叶级数求解带周期性边界条件的热方程

嘿，我来帮你理顺这个思路！你现在的方向其实没大错，但在假设解的形式上绕了点弯路，咱们重新一步步来，把逻辑理清楚：

首先先明确我们要解的问题：

$k>0,$
$$\begin{cases}
u_t = ku_{xx}, -\pi \le x \le \pi, t>0\
u(x,0)=f(x)\
u(-\pi,t)=u(\pi,t), u_x(-\pi,t)=u_x(\pi,t), t>0
\end{cases}$$

你想到用复傅里叶级数是完全正确的，因为周期性边界条件刚好和复傅里叶的基函数 $e^{inx}$ 完美匹配——这些基函数本身就满足 $e^{in(-\pi)}=e{in\pi}$，它们的导数 $in e^{inx}$ 也满足周期性条件，天生适配这个问题的边界约束。

正确的解形式假设 #

你之前的问题在于，先去假设了 $u_{xx}(x,t)$ 的形式，这有点舍本逐末了。我们应该直接对整个解 $u(x,t)$ 做假设：既然初始条件 $f(x)$ 可以展开为复傅里叶级数，那对于任意时刻 $t$，$u(x,t)$ 也可以表示为关于 $x$ 的复傅里叶级数，只是系数会随时间 $t$ 变化，也就是：
$$u(x,t) = \sum_{n=-\infty}^{\infty} c_n(t) e^{inx}$$
这里的 $c_n(t)$ 是依赖于时间的系数，当 $t=0$ 时，$c_n(0)$ 就是你定义的 $f(x)$ 的复傅里叶系数：
$$c_n(0) = \frac{1}{2\pi}\int_{-\pi}^{\pi} f(x)e^{-inx}dx$$

代入热方程推导系数的微分方程 #

接下来把这个假设的 $u(x,t)$ 代入热方程 $u_t = ku_{xx}$，分别计算两边的表达式：

1. 计算左边 $u_t$：对时间 $t$ 逐项求导（热方程的解满足逐项求导的收敛条件，放心用）：
   $$u_t = \sum_{n=-\infty}^{\infty} c_n'(t) e^{inx}$$
2. 计算右边 $ku_{xx}$：对 $x$ 求二阶偏导，每一项 $e^{inx}$ 的二阶导数是 $(in)^2 e^{inx} = -n^2 e^{inx}$，所以：
   $$u_{xx} = \sum_{n=-\infty}^{\infty} c_n(t) (-n^2) e^{inx}$$
   两边乘 $k$ 就得到：
   $$ku_{xx} = \sum_{n=-\infty}^{\infty} c_n(t) (-k n^2) e^{inx}$$

现在热方程要求左右两边相等，而复傅里叶级数的系数是唯一的——也就是说，两边级数中对应同一个 $n$ 的系数必须相等，这样我们就得到了每个 $c_n(t)$ 满足的常微分方程：
$$c_n'(t) = -k n^2 c_n(t)$$

求解系数的微分方程 #

这是个非常简单的一阶线性常微分方程，通解就是指数形式：
$$c_n(t) = c_n(0) e^{-k n^2 t}$$
而我们已经知道 $c_n(0)$ 是 $f(x)$ 的复傅里叶系数，把它代回去就得到了 $c_n(t)$ 的具体形式。

得到最终解 #

最后把 $c_n(t)$ 代回一开始假设的 $u(x,t)$ 级数中，就得到了热方程的解：
$$u(x,t) = \sum_{n=-\infty}^{\infty} \left( \frac{1}{2\pi}\int_{-\pi}^{\pi} f(\xi)e^{-in\xi}d\xi \right) e^{-k n^2 t} e^{inx}$$
或者也可以写成更简洁的形式：
$$u(x,t) = \sum_{n=-\infty}^{\infty} c_n e^{-k n^2 t} e^{inx}$$
这里的 $c_n$ 就是你定义的 $f(x)$ 的复傅里叶系数。

回头看你的思路问题 #

你之前先去假设 $u_{xx}(x,t)$ 的形式，相当于跳过了对整个解的构造，直接去碰二阶导数，这就把自己绕住了。其实核心逻辑和分离变量法是一致的：分离变量法里我们找 $u(x,t)=X(x)T(t)$，而这里的复傅里叶级数，本质上是把所有满足边界条件的 $X(x)=e^{inx}$ 都叠加起来，对应的 $T(t)$ 就是我们解出来的 $e^{-k n^2 t}$，只是用级数的形式统一表示了所有可能的分离变量解的叠加。

备注：内容来源于stack exchange，提问作者JAEMTO