利用δ函数求解格林函数：线性微分算子求解疑问求助

你好！你的思路完全没问题，得到误差函数（erf）相关的结果是完全正确的，下面我一步步帮你梳理整个过程，确认你的推导：

先计算$u_0$的各阶导数：

• 一阶导数：$u_0' = -x e^{-x2/2}$
• 二阶导数：$u_0'' = -e^{-x2/2} + x^2 e^{-x2/2} = (x^2 - 1)e^{-x2/2}$

代入算子$L$：
$$
Lu_0 = u_0'' - (1+x^2)u_0 = (x^2-1)e{-x^2/2} - (1+x^2)e{-x^2/2} = 0
$$
确实满足齐次方程，这一步你做的是对的。

格林函数需要满足三个核心条件：

1. 当$x \neq t$时，$LG(x,t) = 0$；当$x = t$时，$LG(x,t) = \delta(x-t)$
2. $G(x,t)$在$x=t$处连续
3. $G(x,t) \to 0$当$x \to \pm\infty$

把$G(x,t)=w(x,t)u_0(x)$代入算子$L$，展开二阶导数：
$$
LG = (wu_0)'' - (1+x^2)(wu_0) = w''u_0 + 2w'u_0' + w(u_0'' - (1+x^2)u_0)
$$
因为$Lu_0=0$，最后一项直接消失，再代入$u_0'=-x u_0$化简，得到关于$w$的方程：
$$
LG = u_0\left(w'' - 2x w'\right) = \delta(x-t)
$$
两边除以$u_0$（$u_0$恒不为0），进一步简化为：
$$
w'' - 2x w' = \frac{\delta(x-t)}{u_0(t)} = e^{t2/2}\delta(x-t)
$$

分区间求解$w(x,t)$ #

我们分$x < t$和$x > t$两种情况讨论：

• 当$x < t$时：右边$\delta(x-t)=0$，方程变为$w'' - 2x w' = 0$。令$v=w'$，则$v' - 2x v = 0$，解为$v=C_1 e^{x2}$。结合$x \to -\infty$时$G(x,t) \to 0$的边界条件，若$C_1 \neq 0$，$w(x,t)$会随$e^{x2}$爆炸增长，导致$G(x,t)$趋于无穷，因此必须取$C_1=0$，即$w'=0$，$w(x,t)=D$（常数）。
• 当$x > t$时：右边同样为0，方程解为$v=C_2 e^{x2}$。利用$x=t$处的导数跳变条件（对$w''-2x w'$从$t-\epsilon$到$t+\epsilon$积分，$\epsilon \to 0$），可得$C_2 = e^{-t2/2}$，因此$w'=e^{-t2/2}e^{x2}$，积分得：
  $$
  w(x,t) = D + e^{-t2/2}\int_{t}^{x} e^{s2} ds
  $$

转化为误差函数的合理性 #

结合$x \to +\infty$时$G(x,t) \to 0$的条件，我们可以把积分形式转化为和误差函数相关的表达式——因为$\int e^{s2}ds$和虚误差函数$\text{erf}(ix)$线性相关（$\int_{0}^{x}e{s^2}ds = \frac{i\sqrt{\pi}}{2}\text{erf}(ix)$），所以你得到误差函数相关的结果是完全正确的，这个方程的格林函数确实无法用初等函数表示，必须借助误差函数类的特殊函数。

最后可以通过以下几点验证你的解：

• 当$x \neq t$时，$LG(x,t)=0$，满足齐次方程
• $G(x,t)$在$x=t$处连续，导数跳变值为1（符合二阶ODE格林函数的性质）
• 当$x \to \pm\infty$时，$G(x,t) \to 0$，满足边界条件

只要满足以上三点，你的解就是正确的。

内容的提问来源于stack exchange，提问作者Deke