1 基本概念
格林函数法是一类将线性微分方程问题转化为积分表示的分析方法。它把作用于未知函数的线性算子与外部源项、边界条件联系起来，用一个核函数来描述系统对单位激励的响应。借助这种表示，原本难以直接求解的非齐次问题，往往可以改写为形式更清晰的积分方程或显式解。
1.1 格林函数的定义
格林函数通常指线性微分算子在特定边界条件下，对点源输入所产生的响应函数。若把源项看作由许多微小点源叠加而成，那么未知解就可以表示为这些点源响应的加权积分。不同问题对应的格林函数并不相同，通常由算子形式、定义域形状以及边界条件共同决定。
1.2 线性算子与响应核
在格林函数框架中，微分方程常被写成线性算子作用于未知函数等于源项的形式。对应的格林函数可视为该算子的逆作用在核函数层面的表达，因此也常被称为响应核。它体现了输入与输出之间的线性映射关系，使问题具有更强的结构性。
1.3 齐次问题与非齐次问题
齐次问题是指方程右端源项为零的情形，非齐次问题则包含外力、热源、电荷密度等驱动项。格林函数法通常先研究齐次方程及其解空间，再将非齐次部分通过积分叠加引入。这样可以把整体解分解为边界贡献与源项贡献两部分。
1.4 边界条件的作用
边界条件决定格林函数的具体形式，也决定解是否唯一。对于同一个微分算子，若边界条件不同，得到的格林函数往往也不同。边界条件不仅约束解在区域边界上的行为，还会影响积分表示中的边界项，从而对最终结果产生决定性作用。
2 理论基础
格林函数法建立在线性微分方程、算子理论和叠加原理之上。其核心并不只是“求一个特殊函数”，而是利用线性结构把解的生成机制表达出来。理解这一方法，需要先把握基本解、逆算子以及线性组合在解构造中的意义。
2.1 线性微分方程理论
线性微分方程的一个重要特点是解的线性叠加仍为解。正因如此，点源响应可以累加成一般源项对应的解。对于常微分方程和偏微分方程，这种思想都具有普适性，只是具体的表示形式和构造技术有所不同。
2.1.1 常微分方程中的基本解
在常微分方程中，基本解通常对应于单位脉冲输入下的响应。它反映系统对局部激励的瞬时反应，并可通过积分把一般外力作用下的解表示出来。对于高阶线性方程，基本解与初值或边值条件结合后，能给出完整的解析表示。
2.1.2 偏微分方程中的基本解
在偏微分方程中，基本解常与空间维数和算子类型密切相关。例如，椭圆型、抛物型和双曲型方程的基本解结构各不相同。它们描述的是在无边界或理想化自由空间中，单位点源在时空中的传播或扩散方式。
2.2 算子方法
算子方法把微分方程看作抽象空间中的线性映射问题。未知函数是算子的作用对象，格林函数则体现为该映射的核表示。采用这一视角后，许多原本分散的技巧可以统一起来，便于分析可解性与构造解的步骤。
2.2.1 线性算子的逆
若线性算子在给定条件下可逆，则形式上可将未知解写成逆算子作用于源项。格林函数可以看作这种逆的具体积分表达。虽然在实际问题中逆算子未必能直接写出，但通过构造核函数，仍可实现等效的求解效果。
2.2.2 广义逆与积分算子表示
在某些情形下，算子可能不具备严格意义上的逆，或者逆只在特定函数空间中存在。这时常借助广义逆、伪逆或积分算子表示来处理。格林函数法正是将这种抽象逆关系具体化为积分核，从而为非标准边值问题提供可操作的表达。
2.3 叠加原理
叠加原理是格林函数法得以成立的关键。由于方程线性，多个点源的响应可以直接相加，最终形成对连续源分布的积分。该原理使得复杂载荷、分布式热源或电荷密度等问题可以分解为一系列简单响应的组合。
3 格林函数的构造
格林函数的构造通常依赖于方程类型、空间维数以及边界几何。实际处理中，常先求自由空间中的基本解，再根据边界条件进行修正。不同构造路线适用于不同问题，其中一维问题较易处理，多维和复杂区域则需要更多技巧。
3.1 一维情形
一维问题是理解格林函数法的起点。由于变量较少、边界结构清晰，常可以通过分段求解、连续性条件和跃迁条件直接构造格林函数。许多基础教材都以此为例展示方法的完整流程。
3.1.1 常系数方程的构造
对于常系数线性微分方程，格林函数可由特征根结构和匹配条件求得。通常先在源点两侧分别求齐次解，再利用函数连续性及导数的跳跃条件确定系数。这样的构造方式直观且便于写出显式公式。
3.1.2 初边值问题的构造
在初边值问题中，格林函数不仅要满足方程本身，还要满足初始条件和空间边界条件。构造时常将时间与空间变量分离，或通过变换将问题转化到更易处理的形式。最终得到的表示式通常同时包含初值项与边界输入项。
3.2 多维情形
多维问题的格林函数一般比一维情形复杂得多，因为区域形状和边界曲面会显著影响响应分布。构造时常结合对称性、坐标变换和特殊函数工具。若区域简单，可得到较明确的解析表达；若区域复杂，则通常转向近似或数值构造。
3.2.1 自由空间格林函数
自由空间格林函数描述的是无边界约束下点源的基本响应。它往往是多维问题构造的出发点，尤其适合处理无限域或局部近似问题。该函数具有较强的普适性，是许多边界修正方法的基础。
3.2.2 有界区域中的格林函数
在有界区域中，格林函数必须同时满足方程和边界条件，因此通常比自由空间形式多出修正项。其构造常依赖区域的几何特征，例如球域、矩形域或半空间。对于规则区域，可以借助解析技术得到较完整的表达；对于一般区域，则多采用展开或数值方式。
3.3 基本解与边界修正
实际边值问题中，基本解往往只能满足方程而不能自动满足边界条件，因此需要进行边界修正。修正的本质是向基本解中加入适当的齐次解，使其在边界上符合要求。这样得到的整体函数才是真正可用于积分表示的格林函数。
3.3.1 镜像法
镜像法是一种经典的构造技巧，常用于半空间或具有对称边界的区域。其思路是引入“虚拟源”或镜像源，使边界上的条件自动满足。该方法形象直观，常见于静电学和热传导问题中。
3.3.2 特征函数展开法
当区域边界较规则时，可以利用特征值问题的正交特征函数展开格林函数。通过将核函数写成级数形式，能够逐项满足边界条件与微分方程要求。该方法在矩形域、圆域以及某些定边界问题中尤其有效。
4 典型应用
格林函数法不仅是理论工具，也广泛用于具体方程的求解。它能够统一处理源项、边界项和初值项，在物理建模中尤其方便。许多经典方程都可通过格林函数写成积分解形式，从而便于分析和计算。
4.1 常微分方程求解
在常微分方程中，格林函数法常用于求解线性边值问题和受迫振动问题。相比直接积分或待定系数法，它在处理复杂边界条件时更具系统性。尤其当右端项形式较一般时，格林函数能够提供统一的解表示。
4.1.1 二阶边值问题
二阶线性边值问题是格林函数法最常见的应用之一。通过构造满足左右边界条件的核函数，可将未知函数表示为源项的积分。该表示不仅便于求解，还能直接反映解对不同位置载荷的敏感程度。
4.1.2 初值问题的积分表示
对于初值问题，格林函数可将系统响应写成初始状态与外力作用的卷积或积分形式。这样一来，初始位移、初始速度等信息都能被统一纳入表达式。该形式在动力系统和控制理论中也具有重要价值。
4.2 偏微分方程求解
偏微分方程中的格林函数法用途更广，常见于扩散、势场和波动问题。通过合适的核函数，可以把空间和时间变量的影响分离或组合起来，使复杂方程转化为更易处理的积分关系。
4.2.1 泊松方程
泊松方程描述静态场中源项与势函数之间的关系。格林函数可将势函数表示为源密度与核函数的积分，从而直接反映空间分布对场量的影响。它在电势、重力势和稳态热问题中都很常见。
4.2.2 拉普拉斯方程
拉普拉斯方程对应无源情形，常用于描述平衡态问题。格林函数在这里主要用于由边界数据重建区域内部的函数值。由于方程本身没有体源项，边界上的信息在解的确定中占据核心地位。
4.2.3 热方程
热方程中的格林函数反映热量随时间扩散的过程。其核函数通常具有平滑化特征，表示局部温度扰动会逐渐扩散到更大区域。利用该方法可以清晰描述初始温度分布与外部热源对系统演化的影响。
4.2.4 波动方程
波动方程对应传播型过程，格林函数体现信号或扰动的传播速度与传播路径。与热方程相比，波动问题中的响应更强调因果性和有限传播速度。借助格林函数，可分析脉冲激励在空间中的传播与反射。
4.3 物理场问题
在物理建模中，格林函数常被用来描述场量对源分布的响应。它在静电、热学和振动分析中都非常常见，能够把物理直觉与数学表示连接起来。通过核函数，可以直接看到“源在何处，场如何变化”。
4.3.1 静电势问题
静电势问题中，电荷分布与电势之间通常满足泊松方程或拉普拉斯方程。格林函数提供了从电荷密度到电势的积分关系，因此可用于研究导体附近的势分布或边界诱导效应。该方法在理论电磁学中具有基础地位。
4.3.2 热传导问题
热传导问题常借助格林函数研究温度场随时间和空间的变化。无论是稳态还是非稳态过程，只要方程保持线性，都可以用对应核函数来表示热源的影响。它在材料分析、工程传热和环境模拟中都十分常用。
4.3.3 受迫振动问题
受迫振动系统中，外力作用会引起结构或介质的响应。格林函数法可以把位移、速度或应力表示为外力历史的积分结果。该方式有助于分析共振、瞬态响应以及不同激励形式下的系统行为。
5 方法扩展
格林函数法并不局限于直接构造核函数，还与若干经典恒等式和算子理论密切相关。通过这些扩展工具，可以进一步推导积分表达、处理更一般边界条件，并建立时域与频域之间的联系。
5.1 格林第二公式
格林第二公式是连接区域内部函数值与边界积分的重要恒等式。它常用于推导格林函数的积分表示，并在证明唯一性时发挥作用。该公式把微分算子作用下的体积分与边界项联系起来，是边值问题分析中的核心工具。
5.2 格林第三公式
格林第三公式是在第二公式基础上进一步得到的表示公式，常用于将区域内函数表示为边界数据与源项的组合。对于满足适当光滑条件的函数，它能直接给出点值表达式。由此，区域内部的解可由边界信息和格林函数共同确定。
5.3 伴随算子与对偶形式
在更一般的线性问题中，伴随算子提供了从对偶角度理解格林函数的途径。通过引入伴随边值问题，可以处理非自伴算子或非对称边界条件下的情况。对偶形式也有助于建立更稳定的积分表示和理论估计。
5.4 时域与频域表示
格林函数既可以在时域中描述瞬态响应，也可以在频域中表征稳态或谐波响应。通过傅里叶变换或拉普拉斯变换，常可在两个表示之间互相转换。频域表示往往更便于代数处理，而时域表示则更直观地反映传播与演化过程。
6 数值方法中的应用
当解析构造过于困难时，格林函数仍可作为数值方法的基础。许多数值算法并非完全放弃解析结构，而是先把问题改写成积分方程，再对积分核进行离散化。这样既保留了边界信息，又有利于减少高维体网格的负担。
6.1 边界元方法
边界元方法以边界积分方程为核心，通常直接建立在格林函数或基本解之上。与体网格方法相比，它主要离散边界而非整个区域，因此在无限域和外部问题中很有优势。该方法特别适合势场和弹性问题。
6.2 积分方程离散化
将格林函数表示的积分方程离散化后，可转化为线性代数方程组求解。离散过程通常涉及数值积分、基函数选取和核函数近似。只要离散方案设计合理，原问题的边值约束和源项效应都能较好保留。
6.3 快速求解与近似算法
由于格林函数法常导致密集矩阵或高成本积分，快速算法便成为重要补充。常见思路包括核函数降阶、分块近似和多层加速技术。它们的目标是在尽量保持精度的同时，降低存储和计算开销。
7 优点与局限
格林函数法兼具理论美感与实际效用，但它并非对所有问题都同样方便。其优势主要体现在结构清晰和表示统一，而局限则来自构造难度、几何复杂性以及高维计算成本。理解这些特点，有助于在合适场景中使用该方法。
7.1 解析表示的优势
格林函数法最大的优点是能够给出明确的解析或半解析表示。这样的表达式便于分析解的性质，例如连续性、光滑性、衰减行为以及对源项的依赖关系。它也常作为检验数值结果正确性的基准。
7.2 对几何与边界条件的依赖
该方法对区域几何和边界条件较为敏感。边界越复杂，格林函数越难显式构造，所需修正也越繁琐。对于非规则区域、混合边界条件或强非均匀介质，方法往往需要与近似技术结合使用。
7.3 高维问题中的计算复杂度
在高维问题中，格林函数相关积分通常更难计算，存储量和运算量也会迅速增长。即便存在解析核，实际数值实现仍可能面临效率瓶颈。因此，在大规模工程计算中，常需借助降维、压缩或快速近似策略。
8 相关概念
格林函数法与若干基础概念密切相连。这些概念从不同角度解释了“响应”“展开”和“积分表示”的含义，也构成了理解该方法的知识网络。它们在数学物理与应用数学中都十分常见。
8.1 基本解
基本解是微分算子在点源作用下的特解，常被视为构造格林函数的起点。它通常不满足边界条件，但在无边界或局部分析中具有基础意义。许多格林函数都可以看作基本解加上边界修正项。
8.2 冲激响应
冲激响应描述系统对瞬时单位输入的输出，是工程和物理中非常重要的概念。格林函数在很多情形下正是冲激响应的数学表达。通过它可以把复杂输入拆解为大量短暂激励的叠加。
8.3 特征函数展开
特征函数展开是将解写成一组正交函数线性组合的方法。它常与格林函数法互补，尤其适用于边界规则的问题。二者结合时，往往可以得到级数形式的核函数或解表达式。
8.4 卷积表示
卷积表示是线性系统响应的一种标准形式，常见于时域分析和信号处理。格林函数法得到的积分表达在结构上与卷积十分相近，尤其在平移不变系统中几乎可以直接对应。它为理解“输入如何生成输出”提供了统一语言。