热传导方程的分离变量法与傅里叶级数展开

热传导方程，或扩散方程，是描述热量、物质等扩散过程的基本偏微分方程。分离变量法是求解这类线性方程的一种经典方法，其核心思想是将多变量问题转化为多个单变量问题，而傅里叶级数展开则是确定最终解中未知系数的关键工具。

1. 明确问题与方程形式

设定一根长度为 $L$ 的均匀细杆，其温度分布 $u(x,t)$ 随位置 $x$ 和时间 $t$ 变化。杆的两端保持恒温（例如 0 度），且已知初始时刻的温度分布。该问题可用一维热传导方程及定解条件描述。

写出控制方程（热传导方程）：
$$ \frac{\partial u}{\partial t} = k \frac{\partial^2 u}{\partial x^2}, \quad 0 < x < L, \quad t > 0 $$
其中 $k > 0$ 是热扩散系数，表征材料的导热能力。

列出定解条件：

• 边界条件（两端温度固定为0）：
  $$ u(0, t) = 0, \quad u(L, t) = 0, \quad t > 0 $$
• 初始条件（初始温度分布为已知函数 $f(x)$）：
  $$ u(x, 0) = f(x), \quad 0 \leq x \leq L $$

2. 应用分离变量法

假设解具有变量分离的形式，即 $u(x,t) = X(x) \cdot T(t)$。其中 $X(x)$ 仅是 $x$ 的函数，$T(t)$ 仅是 $t$ 的函数。

代入偏微分方程：
$$ X(x) \frac{dT(t)}{dt} = k \frac{d^2 X(x)}{dx^2} T(t) $$

分离变量，将方程两边同时除以 $k \cdot X(x) \cdot T(t)$：
$$ \frac{1}{k} \frac{1}{T(t)} \frac{dT(t)}{dt} = \frac{1}{X(x)} \frac{d^2 X(x)}{dx^2} $$
现在，方程左边仅是 $t$ 的函数，右边仅是 $x$ 的函数。要使等式对任意 $x$ 和 $t$ 成立，两边必须等于同一个常数。设这个常数为 $-\lambda$（负号是为了后续求解方便而引入）。

3. 求解特征值问题

分离变量后，我们得到两个常微分方程：

1. 关于空间部分 $X(x)$ 的方程：
   $$ \frac{d^2 X}{dx^2} + \lambda X = 0 $$
2. 关于时间部分 $T(t)$ 的方程：
   $$ \frac{dT}{dt} + k\lambda T = 0 $$

将边界条件 $u(0,t)=0$ 和 $u(L,t)=0$ 应用于空间部分，得到：
$$ X(0) = 0, \quad X(L) = 0 $$

现在，求解带边界条件的 $X(x)$ 方程。这是一个施图姆-刘维尔问题，其解的形式取决于常数 $\lambda$ 的符号。

• 分析可知，当 $\lambda > 0$ 时，存在非零解满足零边界条件。
• 令 $\lambda = \left(\frac{n\pi}{L}\right)^2$，其中 $n$ 是正整数。
• 得到一族非零解（特征函数）：
  $$ X_n(x) = \sin\left(\frac{n\pi x}{L}\right), \quad n=1, 2, 3, ... $$

将每个 $\lambda_n = \left(\frac{n\pi}{L}\right)^2$ 代入时间方程，其解为：
$$ T_n(t) = e^{-k\lambda_n t} = e^{-k\left(\frac{n\pi}{L}\right)^2 t} $$

4. 构造级数解与叠加原理

根据线性齐次方程的叠加原理，原方程满足边界条件的通解可以写成所有特解的线性组合：
$$ u(x, t) = \sum_{n=1}^{\infty} B_n \sin\left(\frac{n\pi x}{L}\right) e^{-k\left(\frac{n\pi}{L}\right)^2 t} $$
这里的 $B_n$ 是待定的常数系数。

5. 利用傅里叶级数确定系数

应用初始条件 $u(x,0) = f(x)$：
$$ u(x, 0) = \sum_{n=1}^{\infty} B_n \sin\left(\frac{n\pi x}{L}\right) = f(x), \quad 0 \leq x \leq L $$

这表明，函数 $f(x)$ 在区间 $[0, L]$ 上需要展开成正弦傅里叶级数。利用正弦函数的正交性，即：
$$ \int_0^L \sin\left(\frac{m\pi x}{L}\right) \sin\left(\frac{n\pi x}{L}\right) dx = \begin{cases} 0, & m \neq n \\ \frac{L}{2}, & m = n \end{cases} $$

两边同时乘以 $\sin\left(\frac{m\pi x}{L}\right)$，并从 $0$ 到 $L$ 积分：
$$ \int_0^L f(x) \sin\left(\frac{m\pi x}{L}\right) dx = \sum_{n=1}^{\infty} B_n \int_0^L \sin\left(\frac{n\pi x}{L}\right) \sin\left(\frac{m\pi x}{L}\right) dx = B_m \cdot \frac{L}{2} $$

解出系数 $B_n$：
$$ B_n = \frac{2}{L} \int_0^L f(x) \sin\left(\frac{n\pi x}{L}\right) dx $$

最终，问题的完整解由级数表达式和系数公式共同给出。