如何显式验证一维波动方程的推迟格林函数?
嘿,这个问题得用广义函数(分布)的微分规则来处理,因为Heaviside阶跃函数不是光滑函数,普通的微分在间断点处是没有定义的,咱们一步步来拆解:
首先明确目标:我们要验证把推迟格林函数 $G_{ret} = \frac{1}{2}\Theta(t - |x|)$ 代入波动方程左边后,结果等于 $\delta(t)\delta(x)$,也就是要证明:
$$\frac{\partial^2 G_{ret}}{\partial t^2} - \frac{\partial^2 G_{ret}}{\partial x^2} = \delta(t)\delta(x)$$
其中 $\Theta$ 是Heaviside阶跃函数,$\delta$ 是Dirac delta函数。
步骤1:计算对t的二阶偏导
先求一阶t偏导:阶跃函数 $\Theta(z)$ 的广义导数就是 $\delta(z)$,这里令 $z = t - |x|$,直接套用链式法则可得:
$$\frac{\partial G_{ret}}{\partial t} = \frac{1}{2}\delta(t - |x|)$$
再对t求二阶偏导,就是对上面的结果再做一次t的广义微分:
$$\frac{\partial^2 G_{ret}}{\partial t^2} = \frac{1}{2}\delta'(t - |x|)$$
这里的 $\delta'(z)$ 是delta函数的一阶分布导数。
步骤2:计算对x的二阶偏导
先求一阶x偏导,注意 $|x|$ 的导数是符号函数 $\text{sgn}(x) = \frac{x}{|x|}$(x≠0时),用链式法则:
$$\frac{\partial G_{ret}}{\partial x} = \frac{1}{2}\delta(t - |x|) \cdot \left(-\frac{d|x|}{dx}\right) = -\frac{1}{2}\text{sgn}(x)\delta(t - |x|)$$
接下来求二阶x偏导,这里要用到乘积的分布微分法则:对于两个广义函数 $f$ 和 $g$,有 $\frac{\partial(fg)}{\partial x} = \frac{\partial f}{\partial x}g + f\frac{\partial g}{\partial x}$。我们令 $f = -\frac{1}{2}\text{sgn}(x)$,$g = \delta(t - |x|)$,分两项计算:
第一项:$\frac{\partial f}{\partial x}g$
符号函数的广义导数是 $2\delta(x)$(因为 $\text{sgn}(x)$ 在x=0处的跃变是2,对应导数为 $2\delta(x)$),所以:
$$\frac{\partial f}{\partial x} = -\frac{1}{2} \cdot 2\delta(x) = -\delta(x)$$
代入后,只有当x=0时 $\delta(x)$ 非零,此时 $|x|=0$,$\delta(t - |x|) = \delta(t)$,因此这一项简化为:
$$-\delta(x)\delta(t)$$第二项:$f\frac{\partial g}{\partial x}$
先计算 $\frac{\partial g}{\partial x}$,用链式法则:
$$\frac{\partial \delta(t - |x|)}{\partial x} = \delta'(t - |x|) \cdot \left(-\frac{d|x|}{dx}\right) = -\text{sgn}(x)\delta'(t - |x|)$$
代入后,因为 $\text{sgn}^2(x) = 1$(x≠0时成立,x=0处不影响结果),所以这一项简化为:
$$\frac{1}{2}\delta'(t - |x|)$$
把两项加起来,就是对x的二阶偏导:
$$\frac{\partial^2 G_{ret}}{\partial x^2} = -\delta(t)\delta(x) + \frac{1}{2}\delta'(t - |x|)$$
步骤3:代入波动方程左边验证
现在把两个二阶偏导代入方程左边:
$$
\begin{align*}
&\frac{\partial^2 G_{ret}}{\partial t^2} - \frac{\partial^2 G_{ret}}{\partial x^2} \
=&\frac{1}{2}\delta'(t - |x|) - \left(-\delta(t)\delta(x) + \frac{1}{2}\delta'(t - |x|)\right) \
=&\frac{1}{2}\delta'(t - |x|) + \delta(t)\delta(x) - \frac{1}{2}\delta'(t - |x|) \
=&\delta(t)\delta(x)
\end{align*}
$$
你看,delta函数的导数项刚好抵消,剩下的就是我们预期的 $\delta(t)\delta(x)$,这就验证了这个推迟格林函数满足波动方程的格林函数条件!
如果觉得直接求分布导数有点抽象,也可以用测试函数法验证:找一个光滑的紧支集测试函数 $\phi(t,x)$,通过分部积分把格林函数的导数转移到测试函数上,最终也能得到积分结果等于 $\phi(0,0)$,这是广义函数理论里更严谨的验证方式,思路和上面一致。
备注:内容来源于stack exchange,提问作者Brain Stroke Patient

