格林函数方法

物理问题引入：为什么微分方程可以变成积分卷积？

在电磁学中求解静电势时，泊松方程是一个空间偏微分方程：

\nabla^{2} ϕ (r) = - \frac{ρ (r)}{ε_{0}}

nabla^2 phi(r) = -rho(r)varepsilon_0

如果你直接面对复杂的连续带电体（如带电云层、带电导体），尝试直接求解偏微分方程是极其困难的．然而，大一普通物理就教给过我们一个积分公式：

ϕ (r) = \frac{1}{4 π ε_{0}} ∭ \frac{ρ (r^{'})}{| r - r^{'} |} d^{3} r^{'}

phi(r) = 14pivarepsilon_0iiint rho(r')|r - r'| d^3r'

** 为什么一个二阶偏微分方程能够被彻底逆转为一个求和积分？乔治·格林（George Green）在 1828 年创立的 格林函数法（Green's Function Method）** 揭示了这种转化背后的普适数学思想：先求单个点源的激发场，再通过连续积分线性叠加出全空间任意场源的响应！

1. 格林函数的基本定义与代数思想

设线性自共轭微分算子为 $L$ L，考虑带有复杂非齐次源项 $f (r)$ f(r) 的定解方程：

L u (r) = f (r) (r \in Ω)

L u(r) = f(r) quad (r in Omega)

且在边界 $\partial Ω$ dOmega 上满足齐次边界条件（如 ${u |}_{\partial Ω} = 0$ left.uright|_dOmega = 0）．

1.1 点源响应方程

定义 格林函数 $G (r, r^{'})$ G(r, r') 为位于 $r^{'}$ r' 处的理想单位点源 在观测点 $r$ r 所激发出的场响应：

L_{r} G (r, r^{'}) = δ (r - r^{'})

L_r G(r, r') = delta(r - r')

其中 $G (r, r^{'})$ G(r, r') 在边界上满足与原问题相同的齐次边界条件．

1.2 积分卷积直接给出通解

利用狄拉克 $δ$ delta 函数的筛选性质，任意连续物理源 $f (r)$ f(r) 均可表示为点源的连续叠加：

f (r) = \int_{Ω} f (r^{'}) δ (r - r^{'}) d r^{'}

f(r) = int_Omega f(r') delta(r - r') dr'

根据微分算子 $L$ L 的线性叠加性，总解直接由点源响应的积分卷积给出：

u (r) = \int_{Ω} G (r, r^{'}) f (r^{'}) d r^{'}

u(r) = int_Omega G(r, r') f(r') dr'

微分算子的逆运算

从线性代数的视角看：若将微分算子视作无限维连续矩阵 $L$ L，方程为 $L u = f$ L u = f．则格林函数 $G$ G 本质上就是微分算子的 逆矩阵（逆算子）： $u = L^{- 1} f$ u = L^-1 f！

2. 格林函数的对称互易性（Reciprocity）

由 Sturm–Liouville 自共轭性及格林第二恒等式，对于满足齐次自伴边界条件的体系，格林函数在观测点 $r$ r 与源点 $r^{'}$ r' 的对换下严格对称：

G (r, r^{'}) = G (r^{'}, r)

G(r, r') = G(r', r)

物理图像（声学与电磁学互易定理）：在 $r^{'}$ r' 处敲一下音叉，在 $r$ r 处测得的声压；严格等于把音叉移到 $r$ r 处敲击，在 $r^{'}$ r' 处测得的声压！

3. 本征函数展开法（谱表示）

若微分算符 $L$ L 对应的自共轭本征值问题已被解出：

L ψ_{n} (r) = λ_{n} ψ_{n} (r)

L psi_n(r) = lambda_n psi_n(r)

且本征函数系正交归一完备 $\int ψ_{n}^{*} (r) ψ_{m} (r) d r = δ_{n m}$ int psi_n^*(r)psi_m(r)dr = delta_nm，根据完备性， $δ$ delta 函数可展开为：

δ (r - r^{'}) = \sum_{n = 1}^{\infty} ψ_{n} (r) ψ_{n}^{*} (r^{'})

delta(r - r') = sum_n=1^inf psi_n(r) psi_n^*(r')

尝试将格林函数按此基底展开： $G (r, r^{'}) = \sum_{n = 1}^{\infty} c_{n} (r^{'}) ψ_{n} (r)$ G(r, r') = sum_n=1^inf c_n(r') psi_n(r)．代入点源方程：

\sum_{n = 1}^{\infty} c_{n} (r^{'}) λ_{n} ψ_{n} (r) = \sum_{n = 1}^{\infty} ψ_{n} (r) ψ_{n}^{*} (r^{'}) ⟹ c_{n} (r^{'}) = \frac{ψ_{n}^{*} (r^{'})}{λ_{n}}

sum_n=1^inf c_n(r') lambda_n psi_n(r) = sum_n=1^inf psi_n(r) psi_n^*(r') implies c_n(r') = psi_n^*(r')lambda_n

由此得到格林函数的 本征函数谱表示：

G (r, r^{'}) = \sum_{n = 1}^{\infty} \frac{ψ_{n} (r) ψ_{n}^{*} (r^{'})}{λ_{n}}

G(r, r') = sum_n=1^inf psi_n(r) psi_n^*(r')lambda_n

这一公式将微分方程的求逆与线性代数的谱分解（特征值倒数求和）完美融为一体！

4. 物理经典范例

4.1 静电场镜像法：半空间格林函数的构造

考虑在无穷大接地导体平板（ $z = 0$ z=0 面，电势 $ϕ = 0$ phi=0）上方 $z^{'} = d$ z'=d 处放置一点电荷 $q$ q．要求解满足导体面边界的半空间格林函数，直接利用 镜像点源构造：

在实源点 $r^{'} = (0, 0, d)$ r' = (0, 0, d) 放置单位正源；
在虚源点 $r^{″} = (0, 0, - d)$ r'' = (0, 0, -d) 放置单位镜像负源；

G (r, r^{'}) = \frac{1}{4 π \sqrt{x^{2} + y^{2} + (z - d)^{2}}} - \frac{1}{4 π \sqrt{x^{2} + y^{2} + (z + d)^{2}}}

G(r, r') = 14pix^2 + y^2 + (z - d)^2 - 14pix^2 + y^2 + (z + d)^2

当 $z = 0$ z=0 时，两项严格抵消， $G |_{导体面} = 0$ G|_导体面 = 0！普物中神乎其技的「镜像法」，本质上就是通过寻找几何对称虚源来强行满足格林函数的齐次边界条件！

4.2 电动力学时变波动方程与推迟势

对于真空中四维达朗贝尔波动方程 $◻ A = - μ_{0} J$ square A = -mu_0 J，由于光速 $c$ c 有限，响应必须服从因果律（响应时间 $t$ t 晚于激发时间 $t^{'}$ t'）．利用傅里叶变换与围道积分，解出波动方程的 推迟格林函数（Retarded Green's Function）：

G (r, t; r^{'}, t^{'}) = \frac{1}{4 π | r - r^{'} |} δ (t - t^{'} - \frac{| r - r^{'} |}{c})

G(r, t; r', t') = 14pi |r - r'| deltaleft( t - t' - |r - r'|c right)

将此格林函数与电流密度 $J (r^{'}, t^{'})$ J(r', t') 卷积，立刻导出发射电磁辐射的著名 推迟势（Retarded Potentials）：

A (r, t) = \frac{μ_{0}}{4 π} ∭ \frac{J (r^{'}, t - \frac{| r - r^{'} |}{c})}{| r - r^{'} |} d^{3} r^{'}

A(r, t) = mu_04pi iiint Jleft(r', t - |r-r'|cright)|r - r'| d^3r'

5. 学习衔接

下一大阶段：掌握柱坐标与球坐标下的特殊函数解法，进入特殊函数族导学；
量子力学接口：格林函数在量子力学微扰展开中对应 传播子（Propagator）与预解算子．

本页面最近更新：2026/10/1 10:49:01，更新历史
发现错误？想一起完善？在 GitHub 上编辑此页！
本页面贡献者：Leafuke, Physics Learning Wiki Team
本页面的全部内容在 CC BY-SA 4.0 和 SATA 协议之条款下提供，附加条款亦可能应用