含狄拉克δ函数边界条件的偏微分方程求解困惑

我完全理解你查过类似问题但还是卡壳的感受——带狄拉克δ函数的边界条件确实容易绕晕，咱们先从细节入手理清楚：

首先纠正一个很容易犯的小笔误：你写的Y(y)表达式里变量用了x，应该是y才对！正确的形式应该是：
$$Y(y)=A\cosh\left(\frac{k\pi}{a}y\right) + B\sinh\left(\frac{k\pi}{a}y\right)$$
这个小错误可能直接导致你混淆了边界条件的逻辑，咱们接着往下说。

根据第一个边界条件$Y(0)=0$，代入后$\cosh(0)=1$，$\sinh(0)=0$，所以$A=0$，$Y(y)$简化为：
$$Y(y)=B\sinh\left(\frac{k\pi}{a}y\right)$$
这部分你是对的，但问题出在第二个边界条件的理解上：你把$U(x,b)=K\delta(x-x_0)$直接当成了$Y(b)=K\delta(x-x_0)$，这是核心误区！

因为咱们用分离变量法得到的是$U(x,y)=X(x)Y(y)$，所以$U(x,b)=X(x)\cdot Y(b)$——这里$Y(b)$是一个常数（因为b是给定的常数），而$\delta(x-x_0)$是关于x的分布函数，显然不可能让一个常数等于它。正确的逻辑是：这个边界条件是在给x方向的特征函数$X(x)$赋予系数，需要用特征函数的正交性来处理。

通常分离变量法中，x方向会得到一组正交的特征函数族，比如假设$X_k(x)$是满足x方向齐次边界条件的特征函数（比如常见的$\sin\left(\frac{k\pi x}{a}\right)$），我们需要把解展开成级数形式：
$$U(x,y)=\sum_{k=1}^\infty X_k(x)Y_k(y)$$

把这个级数代入边界条件$U(x,b)=K\delta(x-x_0)$，得到：
$$\sum_{k=1}^\infty X_k(x)Y_k(b)=K\delta(x-x_0)$$

接下来利用特征函数的正交性：对任意两个特征函数$X_k(x)$和$X_m(x)$，有$\int_0^a X_k(x)X_m(x)dx=N_k\delta_{km}$（其中$N_k$是归一化常数，$\delta_{km}$是克罗内克δ，仅当k=m时为1，否则为0）。我们给等式两边同时乘以$X_m(x)$，并在$[0,a]$上积分：

左边积分结果为：
$$\int_0^a \sum_{k=1}^\infty X_k(x)X_m(x)Y_k(b)dx = Y_m(b)\cdot N_m$$

右边积分利用狄拉克δ函数的筛选性质，结果为：
$$\int_0^a K\delta(x-x_0)X_m(x)dx=KX_m(x_0)$$

由此可以解出$Y_m(b)=\frac{KX_m(x_0)}{N_m}$，再结合$Y_k(y)=B_k\sinh\left(\frac{k\pi}{a}y\right)$，代入$y=b$可得$B_k=\frac{Y_k(b)}{\sinh\left(\frac{k\pi}{a}b\right)}=\frac{KX_k(x_0)}{N_k\sinh\left(\frac{k\pi}{a}b\right)}$。

最后把系数代入级数，就能得到满足所有条件的解了。

总结一下，你之前的主要问题是误将分离变量后的乘积形式边界条件拆成了Y(b)单独等于δ函数，忽略了X(x)的存在；另外变量笔误也加重了混淆，修正这两点后，利用特征函数正交性就能顺利处理δ函数的边界条件啦。

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