离散傅里叶变换中边界条件的施加方法及特定稳态边值问题的求解疑问

你观察得非常准确——直接对带Dirichlet边界条件的有限差分方程做离散傅里叶变换（DFT）确实会遇到除以0的问题，核心原因是DFT天然假设信号是周期性的，而Dirichlet边界条件直接打破了这种周期性，导致频域表达式在直流分量（对应(k=0)）处出现奇异性。下面我会一步步讲解怎么正确结合边界条件来用傅里叶类变换求解你的问题：

一、先明确问题的核心矛盾 #

你的离散方程组其实分两部分：

• 内部点（(i=1,2,...,N-2)）：(\frac{u_{i+1} + u_{i-1} - 2u_i}{h^2} = 2)
• 边界点：(u_0=0)，(u_{N-1}=1)

DFT会默认 (u_N = u_0)、(u_{-1}=u_{N-1})，这和你的Dirichlet边界完全冲突，所以直接对整个序列做DFT会引入错误的边界假设，进而出现除以0的情况。

二、方法1：改用适配Dirichlet边界的离散正弦变换（DST） #

既然DFT不适合非周期边界，我们可以用第二类离散正弦变换（DST-II），它的基函数天然满足 (f(0)=f(L)=0)，完美匹配你的边界条件，步骤如下：

1. 对离散方程两边同时乘以 (h^2)，得到：
   [u_{i+1} + u_{i-1} - 2u_i = 2h^2]
2. 对等式两边做DST-II变换，利用DST的差分性质：
   二阶差分的DST变换为：(\mathcal{DST}{u_{i+1} + u_{i-1} - 2u_i} = -4\sin^2\left(\frac{k\pi}{2N}\right)\hat{u}_k)（其中(k=1,2,...,N-1)是DST的频点，(N)是离散点总数）
3. 计算右边常数项的DST变换：
   (\mathcal{DST}{2h^2} = 2h^2 \cdot \frac{2}{N} \sum_{i=0}^{N-1} \sin\left(\frac{k\pi i}{N-1}\right))（这个求和可以直接计算出解析表达式）
4. 解出频域分量(\hat{u}_k)：
   [\hat{u}_k = \frac{2h^2 \cdot \mathcal{DST}{1}}{-4\sin^2\left(\frac{k\pi}{2N}\right)}]
5. 最后做逆DST-II变换，就能得到满足Dirichlet边界条件的离散解，而且不会出现奇异性。

三、方法2：修正DFT的周期性假设，手动引入边界条件 #

如果你一定要用DFT，需要手动修正周期性假设带来的错误，步骤如下：

1. 拆分解为特解+齐次解：
   原微分方程的特解可以直接写出：(u_p(x) = x^2)（因为(u_p''=2)，且刚好满足(u_p(0)=0)、(u_p(1)=1)的边界条件，所以这个特解就是问题的精确解）。对于一般情况，如果特解不满足边界条件，设(u = u_p + u_h)，其中(u_h)满足齐次方程(u_{h,xx}=0)，且边界条件为(u_h(0)=0 - u_p(0))、(u_h(1)=1 - u_p(1))。
2. 修正DFT的边界误差：
   原差分算子在DFT下的表达式是：
   [\mathcal{F}{\text{二阶差分}} = -\frac{4\sin^2(kh/2)}{h2}]
   当(k=0)（直流分量）时，这个值为0，对应齐次方程的解，我们可以单独处理直流分量：结合边界条件解出直流分量的修正值，再对非直流分量做逆DFT，最后合并得到完整解。

四、验证你的问题的精确解 #

其实你的问题的精确解很容易直接求出：对(u_{xx}=2)积分两次得(u(x)=x^2 + Ax + B)，代入边界条件(u(0)=0)得(B=0)，(u(1)=1)得(1 + A = 1)即(A=0)，所以精确解是(u(x)=x^{2)，离散后就是(u_i=(ih)}2)（其中(h=1/(N-1))），你可以用这个结果验证上面两种方法的正确性。

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