2
0说话人1: 哈喽，欢迎来到本期节目。先问你一个问题啊，如果一个物理系统受到的外力特别复杂，一会儿大一会儿小，形状还千奇百怪，咱们怎么求它的响应？
说话人2: 这还用问，硬算呗，把复杂的力代进方程里，跟它死磕。
说话人1: 磕得动吗？比如电荷随便分布、外力随便变化的那种。
说话人2: 呃……那估计得磕到牙都没了。
说话人1: 所以数学家们准备了一把万能钥匙，专门对付这种线性非齐次微分方程，叫格林函数。今天聊的这些内容，是李坚毅博士整理的数理话题，他把这套方法从历史到推导再到实例串成了一条特别清楚的线。
说话人2: 格林函数，听起来像某个人名。
说话人1: 答对了，它真的来自一个叫格林的人，乔治·格林。这位的经历说出来你可能不信，他正经上学只上了大概一年，小小年纪就辍学回家，在父亲的面包坊和风车磨坊里干活。
说话人2: 磨坊主？那他后来怎么成的数学家？
说话人1: 靠借书自学。他一边磨麦子一边啃拉普拉斯、拉格朗日、泊松这些人的著作，研究工作基本是在磨坊顶楼完成的。一八二八年，他自己掏钱印了一篇论文，据说只有五十一个人订阅，在他活着的时候几乎没人搭理。
说话人2: 这也太惨了，明珠暗投啊。
说话人1: 更悬的是他一八四一年去世以后，这篇论文差点彻底失传。直到四年后，二十一岁的开尔文勋爵把它重新挖出来，介绍给欧洲大陆的数学家，大家才发现这东西厉害。后来黎曼他们正式把里面的点源响应函数命名为格林函数。
说话人2: 等会儿，点源响应，这四个字我好像有点感觉，是不是说先搞清楚一个点产生的效果？
说话人1: 一点没错。要理解它，咱们得先从方程长什么样说起。一般的线性非齐次微分方程，可以念成这样：一个叫L的算子作用在未知函数y上，等于右边的源项f。你不用管符号，关键在这个L是线性的。
说话人2: 线性我熟，不就是直线那种感觉嘛。
说话人1: 差不多，精确说就两条规矩。第一，两个函数先加起来再让L作用，等于让L分别作用完再相加；第二，先放大倍数再作用，等于作用完再放大同样的倍数。合起来一句话：先做线性组合再作用，和先作用完再组合，结果完全一样。
说话人2: 听起来像句废话。
说话人1: 这句废话就是整座大厦的地基。它意味着算子可以跟求和、跟积分随便交换顺序。静电学里的泊松方程就是典型，负的拉普拉斯算子作用在电势上，等于电荷密度除以真空介电常数。因为线性，一堆电荷一起放产生的电势，跟每个电荷单独放再叠加，完全相同。
说话人2: 那怎么把一个随便形状的源给拆开呢？
说话人1: 问得太好了，这就要请出第二个主角，狄拉克δ函数。
说话人2: 希腊字母那个德尔塔？
说话人1: 对。它不是普通函数，是个广义函数，你可以把它想成一根又高又细的针：除了针尖那一个点之外，其他地方全是零，可它在全空间积分起来恰好等于一。它最实用的本事叫筛选性质，一个连续的源函数乘上δ函数再在全空间积分，会把源点那一位置的值原封不动地筛出来。
说话人2: 物理上怎么理解这个筛选？
说话人1: 这么说吧，任意一团电荷分布，都可以看成无穷多个点电荷连续加权叠在一起。位于某个位置的那个小点源，它的强度权重就是当地电荷密度乘上体积元。李博士在整理这部分内容的时候特意强调，δ函数的全部意义都由积分等式规定，千万别拿它当普通函数去代数字。
说话人2: 明白，它活着的意义就是待在积分号里。
说话人1: 总结到位。好，现在万事俱备。针对算子L，格林函数就定义成：在某个点放一个单位点源，也就是用δ函数描述的那种源，方程的解。记法上G括号x分号ξ，意思是ξ这个位置施加单位源，在x这个位置观测到的响应。
说话人2: 一个在ξ，一个在x，分得还挺清楚。
说话人1: 接下来就是最漂亮的一步。原方程右边的任意源f，先用δ函数的筛选性质写成积分，也就是无穷多个不同位置、不同强度点源的加权和；然后把点源换成它对应的格林函数，再利用L是线性的，把算子穿过积分号提到外面去。
说话人2: 等等，我好像看到门开了。
说话人1: 门已经开了。对比原方程，直接读出一个特解：未知函数等于格林函数和源函数在全空间做卷积积分。翻成人话就是，总响应等于每个点源响应按强度加权再累加。
说话人2: 哇，也就是说只要把单位点源这一个最标准的情况解明白，天底下所有源就都通关了？
说话人1: 对，这就是降维打击。解一个源，通吃一类方程。咱们马上看个实打实的例子，静电场。
说话人2: 来，我最爱听实算。
说话人1: 泊松方程的格林函数满足负拉普拉斯作用在G上等于δ。在三维无界空间、无穷远处电势趋于零的条件下，这个解就是四π乘距离分之一，也就是说单位点电荷的电势跟距离成反比，库仑定律的样子立刻出来了。
说话人2: 可是凭什么说这个表达式就代表单位点源呢？总不能它说是就是吧。
说话人1: 问得专业，咱们当场验证。以点源为球心取一个半径r的球面，球面面积是四πr平方；电势沿径向求负梯度得到电场，电场大小正好是四πr平方分之一。电场穿过整个闭合球面的总通量，就是场强乘以面积，四πr平方分之一乘四πr平方，等于一。
说话人2: 一！
说话人1: 对，不多不少正好一，这就是单位点源。要是半径换成别的数，r平方照样约掉，通量永远是一，这正是δ函数积分归一的样子。把格林函数代回通解公式，电势就是电荷密度除以四π真空介电常数乘距离，再对全空间积分，和库仑电势的叠加式一模一样。
说话人2: 格林当年研究的就是这个吧？
说话人1: 正是，静电位势边值问题就是这套方法的物理原型。不过这里藏着一个坑，格林函数不是唯一的。
说话人2: 啊？刚建立的信仰又晃了一下。
说话人1: 别慌，道理很简单。方程只规定了L作用G等于δ，你再给G加上任意一个齐次方程的解，也就是让L作用完等于零的函数，等式照样成立，因为δ加零还是δ。
说话人2: 所以到底哪个解算数？
说话人1: 由边界条件和初始条件拍板。格林本人在那篇一八二八年的论文里，就花了很大篇幅讨论导体边界怎么锁定位势。如果边界全是齐次约束，比如无穷远电势为零、振子初始位移和初始速度都是零，那卷积积分的结果会自动继承这套零条件，拿来就能用。
说话人2: 那边界上非要取非零值怎么办？
说话人1: 三条路：要么专门构造一个适配这套边界的格林函数；要么用齐次解去把边界上的偏差抵消掉；要么碰上对称几何体，用镜像法这种巧劲，在边界另一侧对称位置虚放一个源，边界条件自动满足。李博士在整理边界这部分时打了个特别形象的比方，说微分方程给的是一族解，边界条件干的事就是从这一大家子里点名点出唯一的那一个。
说话人2: 这个比喻好，方程负责生孩子，边界负责上户口。
说话人1: 哈哈，就是这个意思。前面看的都是静态问题，咱们最后把格林函数搬到时间里，看一个受迫简谐振子，感受一下它在动力学里怎么发威。
说话人2: 弹簧挂小球那个？
说话人1: 对，质量乘以加速度加上质量乘固有角频率平方再乘位移，等于外力。算子L就是质量乘二阶导加质量乘频率平方那一坨。现在设想在τ时刻，给系统来一个单位瞬时冲量。
说话人2: 咚的一下，敲一锤子？
说话人1: 没错，而且规定零初始条件，锤之前安安静静待着。这个冲量对应的格林函数叫推迟格林函数，它在τ时刻之前恒等于零，因果律嘛，锤子还没落，系统不可能有反应；τ时刻之后，就是一段振幅固定的简谐振动。
说话人2: 能给个具体数字感受一下吗？
说话人1: 太能了。假设质量就是一千克，固有角频率取一弧度每秒，对应周期大约六点二八秒，在第二秒末敲一个单位冲量。冲量等于动量变化，质量是一，所以小球瞬间获得一米每秒的速度增量，之后位移就是正弦函数sin括号t减二，t小于二的时候位移严格是零。
说话人2: 那要是外力不是一锤子，而是连续不断变化呢？
说话人1: 那就把连续外力想成无穷多次微小锤击，每一锤都激发一段同样的简谐振动，相位从落锤那一刻开始算，振幅跟当时的力成正比，最后把所有振动从过程起点积分到当前时刻。这个积分在力学里叫杜哈梅积分，本质就是时间版本的格林函数卷积。
说话人2: 静态空间里叠电势，动态时间里叠振动，内核居然是同一个动作。
说话人1: 对，这就是格林函数思想最迷人的地方，空间和时间、静电和振动，表面八竿子打不着，数学骨架完全一致。二十世纪施温格、费曼、朝永振一郎他们把格林函数搬进量子场论，粒子的传播直接用传播子描述，追根溯源还是磨坊顶楼上那套点源响应。
说话人2: 从一个磨面粉的磨坊，一路讲到量子电动力学，这条路也太浪漫了。
说话人1: 说到这儿，李坚毅博士在内容的最后留了一段感悟，我把原话念给你听：所谓复杂，不过是无数个简单的连续叠加；当你找到那个最基础的单元响应，整个世界给你的回声，就都藏在一次积分里了。
说话人2: 这句话值得品。先认清一个点，再理解整个场。
说话人1: 最后帮大家把今天的钥匙串收个口。第一，格林函数只对线性算子成立，非线性方程没法直接套；第二，定义式只给微分约束，必须配上边界或初始条件才有确定形式；第三，齐次边界直接卷积，非齐次边界记得用齐次解修正；第四，δ函数是广义分布，一切运算都以积分等式为准。
说话人2: 记住这四条，以后再碰上张牙舞爪的源，心里就不慌了。
说话人1: 对，把它切成点，逐个回应，再加权求和——复杂问题就是这么被温柔地拆开的。感谢收听，我们下期再见。
说话人2: 下期见。
