数学物理方法 · 数学物理方程篇·第五章完整讲义
学习总目录

数学物理方法 / 数学物理方程篇

数学物理方程篇·第五章学习讲义

基本解和解的积分表达式

δ抽样 → 单点响应 → 格林函数与镜像 → 热核 → 波动核 → 降维与传播 → 推迟势。按教材顺序展开概念、例题、证明主线与解题条件;附章末全部12个大题的解法与结果,以及自编检查题。

严镇军教材配套讲读教材325—374页36道自编检查题离线SVG公式与图解

这一章不是再塞进一堆互不相干的公式。它只做一件事:先回答“在一个点放一个单位扰动,会发生什么”,再用叠加求任意分布的扰动。 δ负责描述“单位点扰动”,基本解负责描述“响应”,积分负责把各点的响应加起来。

按教材第325—374页的顺序展开,包含§5.1、§5.2的四个小节、§5.3和§5.4的四个小节。为方便学习,本讲义再拆为16节。所有解释、图解与检查题均重新编写;第5.15节对应原书章末全部12个大题及其小题。第5.16节提供本书学习收尾用的训练与自检。

第一轮必须掌握: δ抽样与缩放;基本解和格林函数的定义;三种积分项分别接收什么数据;半空间镜像;热核;一维达朗贝尔、三维球面平均、二维圆盘积分;有限传播和推迟时间。第二轮再补: 挖小球的严格证明、三重傅里叶积分的每一行、一般分布的拓扑定义。它们有理论用途,但不应挡住你先把标准题做对。

本章是原书正文的最后一章。读完文件,只表示有了完整学习材料;真正的完成标准,是合上讲义仍能识别题型、写出核、代入数据并核查条件。

学习顺序:先跟着例题写关键步骤,再独立完成检查题,最后按习题路线回到原书。证明段落均说明用途,不要求第一轮逐行背诵。

数学物理方程篇·第五章 / 5.1

先理解点源再认识δ

教材§5.1,第325—328页

5.1 · 教材§5.1,第325—328页

01 整章地图:源、核、积分各做一件事

先想一个最朴素的问题:如果知道一滴墨水在水中怎样散开,能不能预测一整条墨线怎样散开?在线性模型中,可以。把墨线切成许多小段,每段近似一个点源,把它们的影响加起来;小段无限细,就得到积分。

数学上,把观察点写作 ,源点写作 。 表示“在 放一个单位扰动,过了 后,在 测得的响应”。初始分布为 时,小段贡献是 ,于是

这里不是把 当万能黑箱:它必须满足同一个方程以及相应的边界条件。全空间常系数问题的核常只依赖 ,才能简写成卷积。有限区域一般要保留两个空间变量。

图 P5-1 把分布拆成点源,再把每份响应相加
有边界时K通常依赖x和ξ两点;全空间常系数情形才常简写成x−ξ。线性:每个点的响应按源强加权初始分布 φ(ξ)切成小段,每段近似一点源强是 φ(ξ)dξ单位响应 K(t,x,ξ)源点ξ → 观察点x满足同一方程与边界把所有源点的贡献相加u(t,x) = ∫ K(t,x,ξ) φ(ξ) dξ观察点x保持不动;ξ在积分区域里走遍所有源点。

有边界时K通常依赖x和ξ两点;全空间常系数情形才常简写成x−ξ。

做到什么程度: 看到积分要能指出谁是观察点、谁是积分变量、谁是已知数据、谁是核。暂时不用背任何三维积分公式。

5.1 · 教材§5.1,第325—328页

02 为什么“零点以外为零”还不够定义δ

普通函数如果只在一个点非零,通常积分仍是零。单位点质量却要求总质量是1。因此,不能把“,其余为零”当作普通函数的精确定义。真正可用于计算的定义是它对光滑测试函数的作用:

这是抽样性质:把 在源点的值取出来。严格理论中写成 ,测试函数通常取光滑且具有紧支集的函数。你现在只需把δ看成“放进积分后执行抽样的对象”。

用宽度 、高度 的矩形近似:

它的面积始终为1;当 时,。这是积分意义下趋向δ,不是给δ规定一个普通函数的点值。高斯窄峰也能做同一件事。

图 P5-2 δ近似越窄越高,面积始终为1
让ε趋零,积分对连续测试函数取出f(0),不是逐点得到一个无限高的普通函数。矩形近似:宽 2ε,高 1/(2ε),面积 1−ε+εε=1;高度=0.5;总面积=1。δ由积分作用定义;图中每条矩形仍是普通近似函数。

让ε趋零,积分对连续测试函数取出f(0),不是逐点得到一个无限高的普通函数。

必须避免: 把δ的“高度”代入常规四则运算、计算 、随意约掉δ。一般分布之间的乘法并没有自动定义。

5.1 · 教材§5.1,第325—328页

03 点质量、点热量、瞬时冲量怎么写

总质量 集中在 ,线密度写成 。因为对任意包含该点的区间积分得到 。集中热量 放进均匀细杆,若横截面积 、密度 、比热 ,温度初值满足

所以

δ前的系数不能只凭“单位点源”四个字猜。先写守恒量,再确定系数。教材把某些热量取为 ,并取单位截面积,才使初温恰为δ。

力在 给物体一个冲量 ,则写 。把 跨过 积分,可得 。δ描述的是有限的总冲量,不能理解成“无限大力持续普通时间”。

做题用途: 把题目中的“集中在一点”“瞬时作用”“总量为某数”翻译成正确的初值或源项。空间δ与时间δ的积分变量不同,不要混用。

5.1 · 教材§5.1,第325—328页

04 抽样前先检查根在哪里

因为源点2在积分区间内部。若改成 ,结果是0。如果根恰落在端点,必须交代端点约定;对称正则化常产生半权重,但这不是所有单边定义都自动遵循的规则。本讲义没有声明约定时,抽样公式均取根严格位于区间内部。

不要漏掉积分测度的单位: 的量纲是长度的倒数,使 无量纲;三维δ的量纲是长度的负三次方。

5.1 · 教材§5.1,第325—328页 · 自编检查题

练习1 用总热量确定δ前的系数

均匀无限杆的 ,在 瞬时输入热量18。写出温度初值,并计算 。

展开解析

由 ,得 。第二个积分抽样得到 。先确定总量,再抽样,顺序不能倒过来。

5.1 · 教材§5.1,第325—328页 · 自编检查题

练习2 δ是否真的等于一个无限高的函数

矩形近似的宽度减为原来的一半,高度和面积怎样变?这是否给出了δ在0点的普通函数值?

展开解析

高度变为两倍,面积仍为1。它给出的是对测试函数积分的极限,不能给出一个普通实数或“无限实数”作为δ的点值。记住面积和抽样作用即可。

数学物理方程篇·第五章 / 5.2

δ的运算逐条推出来

教材§5.1,第328—331页

5.2 · 教材§5.1,第328—331页

05 平移、偶性、卷积单位元

源点从0移到 ,写成 。把 换成 不改变作用,因为δ关于自己的中心对称。与光滑函数相乘时有

证明只需两边乘任意测试函数 后积分,均得到 。这叫分布相等,并不是比较它们每个点的数值。

卷积则给出

第二个式子请自己写成积分验证:。用途: 初值为δ时,卷积公式直接留下核本身;证明初值恢复时,δ把 还原成 。

5.2 · 教材§5.1,第328—331页

06 缩放公式:为什么一定有绝对值

对 ,

换元 后,积分测度提供 。若 ,积分上下限也反向;总量不能变成负数,所以这里必须是绝对值。

更一般地, 在各零点 附近光滑且 ,则

例如

流程: 解根 → 检查哪些根在积分域 → 每根除以导数绝对值 → 代入被抽样函数。 在0处是重根,不能硬套该公式。

图 P5-3 复合δ先找简单根,再除以导数绝对值
完整流程:解根,筛选区间,除以导数绝对值,再抽样。示例:δ(x²−4),两个简单根是 −2 与 +2−2+2g′(−2)=−4g′(+2)=+4每根权重 = 1/|g′| = 1/4根不在积分域内就不贡献;重根不能套简单根公式。

完整流程:解根,筛选区间,除以导数绝对值,再抽样。

5.2 · 教材§5.1,第328—331页

07 δ的导数:负号从分部积分来

导数的定义是

更一般地

它保留了分部积分规则。考试最常见错误是把导数也当成普通抽样,漏掉负号。例如 。

两个重要恒等式:

第二个式子对测试函数 的作用为 。更一般地,

注意区分 和 :前者是δ导函数在负自变量处的拉回,等于 ;后者还要乘链式法则的 ,所以等于 。

5.2 · 教材§5.1,第328—331页

08 傅里叶表示不是一个普通收敛积分

沿用上一章约定:

那么 ,形式上写作

这个积分不是逐点收敛的普通积分。可以先截断为 ,得到 ,再在测试函数积分意义下取极限。另一种方式是乘上 ,得到归一化高斯,再令 。

为何要讲: 后面把 代入变换后的方程,才有最简单的初值;三维波核中的δ也由这种形式出现。第一轮不要求你证明分布空间中的收敛定理,但必须知道不能随意交换这些振荡积分。

图 P5-4 频域高斯截断对应空间中的单位窄峰
橙色示较小ε;越接近频域常数1,空间越接近δ。正则化:先保留高斯频率窗,再取极限频域乘子exp(−εk²);ε逐渐减小空间近似exp[−x²/(4ε)] / √(4πε)频率窗变宽空间峰变窄,总量1未加窗的指数积分须按分布解释,不能当普通收敛积分。

橙色示较小ε;越接近频域常数1,空间越接近δ。

5.2 · 教材§5.1,第328—331页 · 自编检查题

练习3 同一题里同时检查根和雅可比

计算 ,以及 。

展开解析

根为 ,导数绝对值均为4。第一个积分为 ;第二个只保留正根,等于 。不能先把负根忽略,再去看区间。

5.2 · 教材§5.1,第328—331页 · 自编检查题

练习4 δ导数与乘积

计算 ,并把 化成以2为中心的δ及其导数。

展开解析

积分为 。乘积恒等式给出 。再对测试函数作用可检查两边均为 。

数学物理方程篇·第五章 / 5.3

多维δ和有限区间的点源展开

教材§5.1末,第330—334页

5.3 · 教材§5.1末,第330—334页

09 多维δ先和体积元一起看

二维点源写成 ,三维再乘 。简记 ,其中 是空间点。抽样为

光滑可逆换元 后,在源点 附近有

这和体积元 正好抵消。极坐标且源点 时,可在角坐标的局部图中写

写 也等价,因为δ把它取在 。角变量按周期理解,不能把同一点重复计数。原点没有唯一角坐标,此式不能直接拿 代入;旋转对称的二维原点δ可用径向表示并结合 理解。

图 P5-5 极坐标面积元与δ中的雅可比互相抵消
小扇区的面积由半径r产生雅可比;点源分布中相应出现倒数。点源强度不能随着坐标系统而改变dA = r dr dθδ中的倒数因子δ二维 → δ径向δ角向/r₀在源点 r=r₀ 抽样乘 dA 后因子抵消本图取 r₀>0;极点没有唯一角坐标。δ不是普通密度数值;它必须连同积分测度一起换元。

小扇区的面积由半径r产生雅可比;点源分布中相应出现倒数。

5.3 · 教材§5.1末,第330—334页

10 正交展开中,δ把系数积分变成函数值

设边界条件对应的一组实特征函数为 ,范数平方为 。在分布意义下,

你已经会算普通函数的傅里叶系数:。现在令 ,抽样后只剩 。这就是基本解与分离变量法的接头处。

两端温度固定为0用正弦;两端绝热用余弦,并保留常数模态。对区间 ,绝热δ展开为

这里“等于”不能理解为在 逐点求和得到一个普通数。

5.3 · 教材§5.1末,第330—334页

11 教材的中点加热杆,完整走一遍

杆长 ,两端绝热,初温 。空间特征函数为 。常数项系数 ;非零模态系数为 。奇数模态消失,偶数 给出 。所以

这是温度方程 的解。逐项积分,余弦项为0,故总积分始终为1; 后只有 留下。你可以同时检查总热量和最终均温。

必须掌握: 常数项为什么不衰减,δ怎样产生特征系数。无需把这个特定中点公式当成新公式背。

图 P5-6 边界决定正弦或余弦,绝热必须保留零模态
示意n=1模态。右图端点斜率为0,不要求端点函数值为0。同一个点源,不同边界筛选不同的空间基零值边界:正弦Xₙ=sin(nπx/L)没有常数模态零导数边界:余弦Xₙ=cos(nπx/L),含n=0常数模态对热方程不衰减0L0Lδ的系数 = Xₙ(ξ)/∫Xₙ² dx;不要忘记除以范数。波动的常数模态:单位初速度产生t,而不是指数衰减。

示意n=1模态。右图端点斜率为0,不要求端点函数值为0。

5.3 · 教材§5.1末,第330—334页

12 教材的正弦振荡积分其实是在找延迟

教材例5在 时计算含 的二重积分。核心恒等式是

取 ;在 且 时,第二个δ没有根。余下 ,根为 ,雅可比为 。因此教材前置系数为 的积分,结果就是 。

意义: 不需要硬算振荡积分;它把输入延迟了传播时间 。若 ,根在积分域外,响应为0。交换积分应通过正则化解释,不能把不收敛的普通积分直接用绝对收敛定理处理。

5.3 · 教材§5.1末,第330—334页 · 自编检查题

练习5 一维区间δ的正弦系数

两端固定温度的 杆,初温为 。写出温度的特征级数。

展开解析

正弦范数平方为 ,系数为 。故 。没有常数模态;热可由两端流失,不能要求总热量一直守恒。

5.3 · 教材§5.1末,第330—334页 · 自编检查题

练习6 多维δ的换元

设 。把原点二维δ写成 中的形式,并检验其总量。

展开解析

雅可比绝对值是6,所以 。乘体积元 后总积分为1。单独变δ、不变面积元,是错误来源。

数学物理方程篇·第五章 / 5.4

静态基本解是算子的点源响应

教材§5.2.1,第334—336页

5.4 · 教材§5.2.1,第334—336页

13 先写定义,再分清两个符号约定

设 是常系数线性空间微分算子。它的基本解 满足

于是 的一个特解为 ,因为 。它只保证方程,并不自动满足指定边界;还可能加任意齐次解。因此“求到基本解”不等于“解完任意边值问题”。

本章格林函数沿用教材另一个约定:,配上零狄利克雷边界。请在纸上固定下面这张表:

对象 定义 三维奇异部分
拉普拉斯算子的基本解
本讲义的狄利克雷格林函数 ,边界为0

若别的教材把 当算子或让 ,相应公式会换符号。考试先写定义,比孤立背一个正负号可靠。

图 P5-7 本章的基本解和格林函数采用不同源符号
其他教材可以选择不同约定;只要定义和表示公式成套一致。先看定义右端的正负号,再写核算子Δ的基本解 ΓΔΓ = +δ三维 Γ = −1/(4πr)用来解 Δu=F:Γ*F本章格林函数 GΔG = −δ,边界为0三维奇异部分 +1/(4πr)用来解 Δu=−f:源项 +∫Gf同一个拉普拉斯,右端换号,核也换号。二维对应:Γ=ln r/(2π);G的奇异部分=−ln r/(2π)。

其他教材可以选择不同约定;只要定义和表示公式成套一致。

5.4 · 教材§5.2.1,第334—336页

14 三维为什么是负的1除以4πr

源点以外应满足 。径向函数满足 ,解为 。取无穷远衰减的形式 ,再对包含原点的小球积分:

因此 。也就是

计算 的拉普拉斯只能得到0,不会看见原点的δ。必须补上小球通量,才能定出源强和符号。这里的小球外法向朝径向外;若在“挖掉小球”的区域上用格林公式,小球内边界的区域外法向反而朝向球心。

图 P5-8 点源强度由小球总通量定出
本图画的是完整小球自身的外法向,通量归一化决定系数与符号。三维:球面面积 4πε² × 径向导数点源ε对 Γ=C/r∂ᵣΓ=−C/ε²总通量=−4πC若源是+δ,C=−1/(4π)只算源点以外的ΔΓ=0,无法确定点源强度。挖孔区域的内边界外法向朝球心,与本图方向相反。

本图画的是完整小球自身的外法向,通量归一化决定系数与符号。

5.4 · 教材§5.2.1,第334—336页

15 傅里叶求法:把微分算子换成一个多项式

对全空间方程做空间傅里叶变换,,,所以形式上

逆变换、用球坐标积分,就得到上一单元的结果。这段教材计算要传达的不是三重积分有多长,而是:常系数微分算子在频域中变成乘法,基本解由其符号的倒数产生。

遇到分母的零点时,逆变换的解释可能需要额外条件;不能对所有 都不加说明地写“除完就唯一”。例如亥姆霍兹方程还涉及辐射条件。第一轮掌握频域方程、符号与归一化,球坐标积分的机械步骤可放第二轮。

5.4 · 教材§5.2.1,第334—336页

16 用一个径向源检验符号

求球 内 、 的正则径向解, 为常数。由 ,原点正则排除 项,得到

当 ,内部为正,符合 为正且源积分 为正的直觉。若公式算出整个球内负值,先检查你到底在解 还是 。

5.4 · 教材§5.2.1,第334—336页 · 自编检查题

练习7 基本解与齐次解

若 ,光滑函数 满足 ,那么 还是基本解吗?这能说明边界条件可省略吗?

展开解析

还是,因为 。这恰好说明基本解本身未必唯一,边界或无穷远条件仍要用来确定具体解。不能把“是基本解”误当成“是该题唯一解”。

5.4 · 教材§5.2.1,第334—336页 · 自编检查题

练习8 通量检查三维系数

候选核 满足 。用小球通量求 。

展开解析

通量为 ,应等于 ,所以 。若目标改成 ,则 。这一步能在没有查表时重建符号。

数学物理方程篇·第五章 / 5.5

格林函数把边界和源一起装进解

教材§5.2.2,第336—341页

5.5 · 教材§5.2.2,第336—341页

17 格林函数比全空间基本解多做了什么

区域为 ,边界为 。对固定源点 ,定义

把它拆成 。第一项负责源点的奇异性; 在整个区域内调和,负责把边界值抵消成0。 的边界值恰为 。

因此,求格林函数是在做一件非常具体的事:保留区域内唯一的单位点源,用一个没有内部新源的修正满足边界。 镜像法就是用区域外的虚拟点源构造这个调和修正。虚拟点源若落进区域,会给方程添上不该有的δ。

5.5 · 教材§5.2.2,第336—341页

18 第二格林公式:只保留本章要用的一行

对足够光滑的 ,散度定理给出

这里 是区域外法向,。令 ,考虑 、。左边成为 ;右边因 ,只留下 。移项得到

这就是本章静态题的主公式。 第一项接收内部源 ;第二项接收给定的边界值 。法向导数对边界积分变量 求,不是随手对观察点 求。

图 P5-9 内部源与边界值从不同入口组成解
格林函数本身边界为0;边界数据通过它的法向导数进入解。Δu=−f,边界u=g:两份数据分别进入内部源体积分:+∫G fG满足 ΔG=−δ零边界部分由源驱动边界值边界积分:−∫g ∂ₙG对边界积分变量求外法向内部无源的调和响应叠加得到所求u回代方程与每一段边界如果题目是Δu=+F,体积分应改为−∫GF。

格林函数本身边界为0;边界数据通过它的法向导数进入解。

严格证明中, 在源点不光滑,不能直接把它塞进只适用于处处光滑函数的定理。可先挖掉半径 的小球,在剩余区域应用公式,再令 ;也可使用分布理论。为什么讲这一点: 解释“δ从哪里进公式”和内边界符号;第一轮不要求逐行证明极限估计。

5.5 · 教材§5.2.2,第336—341页

19 对称性是静态自伴问题的互易性

同一个狄利克雷拉普拉斯问题满足 。证明主线是对两个源点的格林函数应用第二格林公式:两者都在边界为0,边界项消失;内部两个δ各抽样一次,得到交换点位后的相等。

用处一: 核的表达式虽然长,交换观察点和源点后应不变,可用于验算。用处二: 在表示公式中调换变量时不会换成另一个核。

不要把这个结论扩展为“所有有时间的核都在时间上对称”。因果热核、推迟波核只允许过去影响未来。空间算子若非自伴,也通常要引入伴随核,不能照搬这里的互易性。

5.5 · 教材§5.2.2,第336—341页

20 做格林函数题的四项检查

  1. 内部方程:除了真实源点外是否调和;源点强度是否为 。
  2. 边界:把每一段边界代入, 是否为0。
  3. 位置:所有虚拟源是否都在区域外。
  4. 解的条件:外部或无界区域是否还满足要求的衰减、增长或有界性条件。

若边界给的是法向导数,不能直接使用上面的零狄利克雷格林函数公式。在有界纯诺伊曼问题中还存在常数零模态,并需满足通量与源的相容条件。这里强调它,是为了防止你一见“格林函数”就不看边界类型。

5.5 · 教材§5.2.2,第336—341页 · 自编检查题

练习9 从表示公式核对源的符号

已知 、。若题目给的是 、,写解的表示式。

展开解析

令上一公式的 ,得 。可以对体积分施加拉普拉斯检验:。不要只改边界项或只凭记忆判断。

5.5 · 教材§5.2.2,第336—341页 · 自编检查题

练习10 外法向朝哪边

区域为上半空间 ,在平面 上 等于 还是 ?若挖去一个小球,小球边界的区域外法向朝哪里?

展开解析

上半空间的外面在下方,故 。挖孔后,小球内部属于区域外面,所以孔壁的区域外法向朝球心,与小球自身的径向外法向相反。

数学物理方程篇·第五章 / 5.6

三维镜像法从平面做到球

教材§5.2.3,第341—346页

5.6 · 教材§5.2.3,第341—346页

21 半空间:异号镜像用来消去边界值

取 ,源点 ,,镜像点 。定义 、,则

在 上两距离相同,所以边界值为0。镜像点在下半空间,所以它在研究区域内没有新增δ。真实源前的正号满足 。

若要零法向导数,常用同号镜像,让法向导数相消。但那是另一种边界条件,不能把“镜像总是异号”当口诀。这里的全空间点势还应结合无穷远行为,选取符合所求问题的核。

图 P5-10 零值边界用异号像,零导数边界用同号像
左图取垂直截面。两条虚线等长,异号势在平面上抵消。在边界上等距:用符号选择要抵消的量真实源 +1镜像 −1z=0零值边界等距 → 两个势相减G=(1/r₁−1/r₂)/(4π)虚拟源必须在区域外零导数边界则用同号像上半空间的平面外法向朝下:∂ₙ=−∂z。

左图取垂直截面。两条虚线等长,异号势在平面上抵消。

5.6 · 教材§5.2.3,第341—346页

22 半空间泊松核:边界积分中的负号变成正核

对观察点 ,,边界点记作 。由外法向 及镜像核求导,得到

因此 、 的解为该正核与 的平面积分。先试 :核的积分应为1。极坐标检验为

做到什么程度: 会从镜像写出 ,会沿正确法向求导,能用常数边界检查归一化。半空间泊松核不是二维上半平面的 ,维数改变了分母指数。

5.6 · 教材§5.2.3,第341—346页

23 球内镜像的位置不是普通反射

球半径 ,源点 距球心 ,。反演点和权重分别为

格林函数为

在 上,,所以两项抵消。距离反演和电量权重缺一不可。 球心源点 不能直接代入发散的反演坐标,要取极限,得到 。

图 P5-11 球面镜像同时改变位置和权重
图取ρ<R,反演点在球外。源在球心时取极限,不能直接除以ρ=0。反演位置:ρρ*=R²;镜像权重:−R/ρOξξ*边界上的距离比恒为R/ρG = [1/|x−ξ| − (R/ρ)/|x−ξ*|] / (4π)球面反演不是把源点向外平移同一段距离。

图取ρ<R,反演点在球外。源在球心时取极限,不能直接除以ρ=0。

5.6 · 教材§5.2.3,第341—346页

24 球内泊松公式与球心平均值

对 、边界点 ,

在球心,分子为 、距离为 ,核为 。所以球心的调和值等于球面边界值的面积平均。这是前面球谐展开得到的结果,在这里从格林函数重新出现。

例如边界给 ,无需做面积分: 本身调和并满足边界,故是解。也可由唯一性验证。先识别简单调和函数,再决定是否要积分。 题目要求推导泊松公式时,才需要完整求法向导数。

图 P5-12 球心泊松核均匀,球心值等于边界平均
图为球的截面与赤道示意,平均使用球面面积元,不是截面圆周长度。球心没有方向偏好:每块等面积边界权重相同观察点在球心泊松核 = 1/(4πR²)积分边界值再除以面积u(0)=球面平均例如平均 ξᵢ² = R²/3观察点离开球心后,各边界点的权重不再相同。

图为球的截面与赤道示意,平均使用球面面积元,不是截面圆周长度。

5.6 · 教材§5.2.3,第341—346页 · 自编检查题

练习11 球镜像两件事都要算

球半径2,源点 。求反演点、镜像权重,并验证边界点 处抵消。

展开解析

反演点 ,权重 。真实距离1,镜像距离2,所以 。镜像点位于球外;如果误取同权重,则边界不会为0。

5.6 · 教材§5.2.3,第341—346页 · 自编检查题

练习12 球心值不需要做长积分

单位球面边界为 ,求调和解在球心的值。

展开解析

球心值等于球面平均。对称性使三个平方的平均相同,且球面上 ,故各平均为 。答案 。注意边界多项式 本身不调和,不能直接当整个球内的解。

数学物理方程篇·第五章 / 5.7

二维对数核与矩形谱展开

教材§5.2.4,第346—352页

5.7 · 教材§5.2.4,第346—352页

25 二维基本解用圆周通量定系数

二维径向调和函数为 。要求 ,小圆周通量为 ,故

三维是 ,二维是 ,不能只换积分符号而保留原核。二维物理图像常对应三维中沿一个方向均匀延伸的线源。对数中的长度严格说应除以参考长度;改变参考长度只增加常数,不改变拉普拉斯和点源强度。

图 P5-13 二维对数核和三维倒距离核由不同通量定出
基本解分别为ln r/(2π)与−1/(4πr);本章格林函数的奇异部分取反号。通量定系数:周长与面积的维数不同二维圆周径向调和:A+B ln r周长:2πr导数:B/r单位通量给 B=1/(2π)三维球面径向调和:A+B/r面积:4πr²导数:−B/r²单位通量给 B=−1/(4π)这里统一取 ΔΓ=+δ。改变维数不能只改积分符号,核与测度必须一起改变。

基本解分别为ln r/(2π)与−1/(4πr);本章格林函数的奇异部分取反号。

5.7 · 教材§5.2.4,第346—352页

26 上半平面:两个对数合并成距离比

区域 ,真实源 、镜像源 。本章约定的格林函数为

真实源附近为 ,所以 。边界两距离相同,。法向导数给出上半平面泊松公式:

例如 且 ,要求有界并向高处衰减,解为 。用微分和边界检查即可,无需每次重新计算有理函数积分。

5.7 · 教材§5.2.4,第346—352页

27 圆盘:二维反演还要减一个常数

半径 的圆盘,源点距圆心 ,反演点同样是 ,但二维格林函数为

在边界上距离比为 ,额外的 才使对数变为0。与三维的乘权重不同,这里权重进入对数,等价于修正常数。圆心源为 。

圆盘泊松公式为

若 ,解为 。这与第二、三章用分离变量求出的结果一致。不必把不同方法的同一个解背成两份知识。

图 P5-14 圆盘反演的距离比还要乘ρ除以R
图为平面圆盘。省略ρ/R会留下非零常数边界,虽然内部点源强度仍正确。二维圆盘:要让边界上的对数恰好为ln1ρR²/ρG = (1/2π) ln[ ρ|x−ξ*| / (R|x−ξ|) ]边界距离比为R/ρ乘ρ/R后为1,边界值为0二维通过对数中的常数修正;三维是像源强度修正。

图为平面圆盘。省略ρ/R会留下非零常数边界,虽然内部点源强度仍正确。

5.7 · 教材§5.2.4,第346—352页

28 矩形没有简单镜像时,回到特征函数

设 ,四边为零边值。取

范数平方为 ,且 。因此

为什么除以特征值: 拉普拉斯作用在每个模态上会乘 ,正好还原 的展开。系数漏掉范数,或把特征值乘到分子,是最常见的两种错法。

若 ,零边值,立刻得到 。格林函数的双重级数,是把所有这种单模态响应一次装进去。

图 P5-15 矩形格林函数是每个模态除以特征值
对G再施加Δ,−λ抵消分母,得到−δ的特征展开。δ展开 → 每个空间模态单独响应矩形模态Φₘₙ=sin(mπx/A)sin(nπy/B),范数平方AB/4−ΔΦₘₙ=λₘₙΦₘₙ,λₘₙ>0G的每个系数δ系数:Φₘₙ(ξ)/范数平方格林系数:再除以λₘₙ边界由模态自动满足;范数和特征值是两个不同因子。

对G再施加Δ,−λ抵消分母,得到−δ的特征展开。

5.7 · 教材§5.2.4,第346—352页 · 自编检查题

练习13 对数核的符号

候选 在上半平面边界为0。它是否满足本章的 ?

展开解析

不满足。真实源的对数项系数为正,所以产生 。边界对并不保证内部源符号对。应把距离比倒过来,即整体乘 。

5.7 · 教材§5.2.4,第346—352页 · 自编检查题

练习14 用一个模态解矩形泊松方程

单位正方形中 ,四边为0。求解。

展开解析

该模态的特征值为 ,所以 。代入拉普拉斯得到右端,四边都为0。无需把全部双重级数都展开。

数学物理方程篇·第五章 / 5.8

一阶时间方程与热核的通用公式

教材§5.3,第352—357页

5.8 · 教材§5.3,第352—357页

29 一阶时间只需要一个初值

对全空间常系数问题 ,初值为 。定义 、,则

卷积只对空间做,后面还单独对源发生的时间 积分。第一项让初始状态演化 ;第二项让 时刻新产生的一小份源演化剩余的 。

写开一次,避免用星号掩盖变量:

它的积分上限不能变成无穷大,否则把未来源也算进现在。对常系数算子可以利用平移不变性;有边界时要换成满足该边界的 。

图 P5-16 τ时刻产生的源只演化t减τ的时间
箭头表示每份源被推进到同一观察时刻,不表示额外的空间传播轨迹。历史叠加:新生源只经历剩余的时间0现在tτ1τ2τ3τ4源 f(τ,ξ)到现在的响应 U(t−τ)*f(τ)τ越接近t,演化时间t−τ越短。积分只取0≤τ≤t;未来的源不参与现在的解。

箭头表示每份源被推进到同一观察时刻,不表示额外的空间传播轨迹。

5.8 · 教材§5.3,第352—357页

30 为什么源项的系数恰好是1

令 ,用带变上限积分的求导法则:

关键是 ,所以第一项抽样还原 。。初值项同理满足方程,且 。把两项相加,就是教材定理的证明主线。

需要学到哪里: 会独立求这一行,知道为什么是 ,知道最上端产生的是 。证明中交换微分和积分所需的正则性、增长条件可放第二轮;但遇到不收敛的数据,不能默认公式仍给经典解。

5.8 · 教材§5.3,第352—357页

31 从频域推导任意维热核

取 、。变换后

于是 。各坐标的高斯逆变换相乘,得到

一维分母为 ;二维为 ;三维为 。不要把三维指数的幂写成一维的。 一维教材的 与上述写法完全相同。

图 P5-17 热核降峰展宽,总量守恒
灰线为参考时刻;实线为所选时刻。降低的峰被扩大的宽度补偿。一维热核,a=1:总积分始终为1t=0.2;中心高度=0.6308;典型宽度随√t增长。任何t>0时全空间都有非零响应;没有有限传播边界。

灰线为参考时刻;实线为所选时刻。降低的峰被扩大的宽度补偿。

5.8 · 教材§5.3,第352—357页

32 热核必须有的三个性质

第一,,且 。第二,固定 时, 有 ;但对初值卷积则 ,不能把前者的逐点极限拿进积分得到0。第三,峰高降低,典型宽度按 增长。

一维归一化核的二阶矩为 。因此总热量守恒不表示峰值守恒。增加 会加快扩散。

对任意 ,热核在全部空间为正,所以经典热方程没有有限传播速度。这是模型性质;不能给热核人为加上波动方程的 截断。

5.8 · 教材§5.3,第352—357页 · 自编检查题

练习15 时间脉冲应该先传播多久

零初值热方程的源是 ,。写 时的解。

展开解析

当 时为0;当 时为 。源发生前不能有响应;发生后传播的时间是 。恰在冲击时刻以分布方式解释,不需要规定峰的普通点值。

5.8 · 教材§5.3,第352—357页 · 自编检查题

练习16 热核宽度与总量

把 增为两倍,固定时刻一维热核的中心高度、典型宽度和总积分各怎样变化?

展开解析

高度减半,典型宽度加倍,总积分仍为1。因为归一化系数含 ,指数的空间尺度含 。只改指数、不改前系数,会破坏总量。

数学物理方程篇·第五章 / 5.9

热核题从积分走到可检查的结果

教材§5.3应用与本讲义补充

5.9 · 教材§5.3应用与本讲义补充

33 高斯初温:先配方,再读物理量

求 、,。卷积为两个高斯的积分。将指数配方,或使用高斯积分 ,得到

先检查 恢复初温;再检查总积分始终是 。不要误以为初始高斯总量必为1,它有没有归一化取决于前面的系数。

做题要点: 系数、指数中的时间尺度要配套。若空间整体平移到 ,把 换成 即可。

5.9 · 教材§5.3应用与本讲义补充

34 正弦初值、常数源和反应项

若初值 ,不用做卷积积分:它是拉普拉斯的特征函数,解为 。零初值且源为常数 ,全空间均匀解为 ;这也由 得到。

方程改成 ,令 ,就把反应项去掉。点初值的核为

它的总积分为 ,不再总量守恒。若再加源 ,源积分中乘的是 ,因为每份源只经历余下的演化时间。

5.9 · 教材§5.3应用与本讲义补充

35 有边界时,热核也可以做镜像

半轴 、齐次边界,令源点 。固定零温边界用

绝热边界用

前者在 两值抵消;后者在 两导数抵消。它们必须与半轴上的数据做 ,不能再拿全轴积分重复计算。

对 ,空间总积分保持1;对 ,小于1并随时间减少,因为边界会吸收热。有限区间的核也可按正弦/余弦特征级数构造。全空间热核、镜像热核、特征级数是同一个思路在不同边界下的实现。

图 P5-18 热核的镜像相减或相加满足不同端点条件
零温边界允许热流失;绝热边界保持半轴总量。半轴x>0:同一条热核的正负镜像零温:两核相减Kᴰ=H(x−ξ)−H(x+ξ)在x=0值为0绝热:两核相加Kᴺ=H(x−ξ)+H(x+ξ)在x=0导数为000真实源ξ=1,t=0.4,a=1;积分空间仅取半轴。

零温边界允许热流失;绝热边界保持半轴总量。

5.9 · 教材§5.3应用与本讲义补充 · 自编检查题

练习17 高斯初值检验

取 、初温 。写出解和中心温度,并求中心温度降为初始一半的时刻。

展开解析

解为 。中心温度 。令其为 ,得到 ,所以 。

5.9 · 教材§5.3应用与本讲义补充 · 自编检查题

练习18 反应项的时间不能放错

零初值方程 ,源在空间均匀,。求解并说明为什么结果不含 。

展开解析

每一时刻的空间积分 ,故 。均匀数据使扩散项为0,也可直接解 验证。 极限为 。

数学物理方程篇·第五章 / 5.10

二阶时间方程先分清两个初值

教材§5.4.1,第357—362页

5.10 · 教材§5.4.1,第357—362页

36 波动基本解约定的是单位初速度

对 ,初位移 、初速度 。教材定义的 满足

它与热核的初始条件不同。通用公式为

请按初值记:初位移前面多一个时间导数,初速度直接乘核,外力再积时间。 原因不是约定强迫你死记,而是 、:第一项在0时恢复 ;第二项本身初值为0,时间导数恢复 。

源项的一阶时间导数没有端点贡献,因为 ;二阶求导才由 产生 。这和热方程在一阶求导就出现 的机制一致,但时间阶数不同。

图 P5-19 位移、速度、外力的三个入口
位移项的时间导数作用在整个核卷积上,不只作用在一个前系数。先定义 W(0)=0,Wₜ(0)=δ,再按数据组装初位移 φ∂t[W(t)*φ]:在0时恢复φ初速度 ψW(t)*ψ:在0时的时间导数恢复ψ外力 f∫₀ᵗ W(t−τ)*f(τ)dτ:二次求导产生f热方程只需一个初值;波动方程必须核对两个初值。

位移项的时间导数作用在整个核卷积上,不只作用在一个前系数。

5.10 · 教材§5.4.1,第357—362页

37 一维核是一个不断变宽的平台

取 。频域解满足 、、,故

其中 按极限解释为 。对 逆变换得

不要把这个平台误当作初位移的传播形状。它对应的是初速度δ,初位移0。它的空间总积分是 ,因为给单位初速度后总位移的初始增长率为1。

5.10 · 教材§5.4.1,第357—362页

38 把平台核代回去,就是达朗贝尔公式

将 与 卷积,得到 。初位移项多一个时间导数,用变限积分法则:

因此

再加源贡献

题目处理: 先分成位移、速度、外力三份,再分别处理。不要在总公式中一口气代入而忘了某个初值。若源是常数 ,外力贡献应为 ,这是检验系数最简单的办法。

5.10 · 教材§5.4.1,第357—362页

39 一维有限区间的零模态不能照搬成正弦

对两端零导数的区间 ,初位移0、初速度δ,特征展开给

常数空间模态的时间方程为 ,其单位初速度响应是 ,所以出现 。若换成零值边界,改用正弦且没有零模态。任何候选答案都必须回到端点求导或代值检查。 本书该章习题3(2)的答案印刷存在这一类不匹配,详见5.15节。

5.10 · 教材§5.4.1,第357—362页 · 自编检查题

练习19 位移与速度到底放哪里

全轴波动方程 ,,。求解。若交换成 ,结果怎样?

展开解析

第一题为 。第二题为 ,假设 。在 分别求值和求导可判断:不能把两个时间因子交换。 时分别为1和 。

5.10 · 教材§5.4.1,第357—362页 · 自编检查题

练习20 外力积分的系数

全轴波动方程初位移、初速度都为0,源 为常数。用双重积分算出解。

展开解析

内区间长度为 ,故源项为 。对时间求两次导数为 ,空间二阶导数为0,初值也正确。

数学物理方程篇·第五章 / 5.11

三维波核与球面平均

教材§5.4.1例1,第359—362页

5.11 · 教材§5.4.1例1,第359—362页

40 三维核为什么集中在球面

三维频域仍有 。逆变换时利用球对称性,角积分将三重积分降成径向正弦积分;再用δ的傅里叶表示,得到

这不是普通函数,不应该去问球面上的有限高度。它说明点初速度的响应只在半径 的球面上。δ对半径抽样,必须结合球坐标体积元 ,不能把 单独看作源点处可直接代入的数。

一个归一化检查:

这与一维速度核的总积分相同。若你的结果为1,往往把速度核和热核混淆;若得到 ,往往漏掉一个 。

图 P5-20 速度核在一二维内部有值,三维集中在球面
三维示意为球壳截面,不是有限厚度、有限高度的实体薄层。单位初速度的自由波核;支撑半径为atat球面δ壳δ(r−at)/(4πar)三个维数的空间总积分都等于t。一二维核是可积函数;三维核须按分布与球面测度理解。

三维示意为球壳截面,不是有限厚度、有限高度的实体薄层。

5.11 · 教材§5.4.1例1,第359—362页

41 从δ球壳积分到球面平均

固定观察点 ,令球面 ,定义

球坐标换元后

所以三维自由波动解为

这称为三维泊松公式或基尔霍夫公式的球面平均形式。时间导数作用在整个 上,包括球面半径随时间改变的影响。不能只对前面的 求导。

5.11 · 教材§5.4.1例1,第359—362页

42 用球平均做多项式题

把球面点写成 ,。球面平均满足 、。因此

若 、,则

检验:,。在三维里系数是3;二维同题系数是2,不能跨维数套答案。

若 是调和函数并且球平均公式适用,则 ,零初速度时解保持为 。这是快速做题的捷径。

5.11 · 教材§5.4.1例1,第359—362页 · 自编检查题

练习21 三维核的单位初速度检查

若把 错写为 ,它的空间积分是多少?相应初速度总量是多少?

展开解析

积分为 ,在0时的导数为 ,所以初速度总量为 ,不是1。应补上分母中的 。这个检查比直接处理原点δ容易。

5.11 · 教材§5.4.1例1,第359—362页 · 自编检查题

练习22 初速度为二次多项式

三维全空间波动方程 。求解。

展开解析

用速度项 。于是 ,;。注意不能再加一次时间导数,否则初值会错。

数学物理方程篇·第五章 / 5.12

降维法与二维波动

教材§5.4.2,第362—365页

5.12 · 教材§5.4.2,第362—365页

43 降维不是把三维公式里的z删掉

若三维初值与 无关,则解也可取与 无关,满足二维波动方程。但三维球面投影到二维平面时,球面的上、下两半都投到一个圆盘上,且面积元发生变化。

球半径 ,投影半径 ,一半球面的面积元为

两半球面提供因子2。把它带进三维速度项,就得二维核

两个关键点: 积分域是圆盘,不是圆周;边缘有可积奇异性,不是一个有限的边缘值。分母的根号来自投影面积元,因此不能凭外形把三维球壳直接降成二维圆环。

图 P5-21 两个半球面投到圆盘,边缘产生根号权重
二维根号核的几何来源是球面投影面积元,且两半球贡献相等。球面投影到圆盘:上、下两面都要计算面积元变化每半球:dS=at dA/√(a²t²−r²)两半球提供因子2圆盘边缘权重变大降维要积掉横向变量,不能只删除一个坐标。

二维根号核的几何来源是球面投影面积元,且两半球贡献相等。

5.12 · 教材§5.4.2,第362—365页

44 二维完整公式:初位移仍然有时间导数

令 ,则

求时间导数时,积分区域也在动,边缘核又奇异,不能直接把一个普通端点值代入变限求导公式。安全做法是先积分,或令 ,把区域固定成单位圆盘,再求导。

例如二维速度为常数1,极坐标积分给

所以零初位移时 。它同时验证了二维核的总积分为 。

5.12 · 教材§5.4.2,第362—365页

45 教材圆盘初位移例题:为什么会出现负的尾部

初速度0,初位移在单位圆盘内为1、外为0。只求中心 ,对 ,先算位移核的积分再求导:

当 ,积分为 ,故 。当 ,积分为 ,故

后一个量为负,且 时趋于0。波动方程不保非负,所以这不是符号错误。 附近的奇异来自初值在圆周上的跳跃,不能声称此处是普通光滑解;平滑初值会改变这个尖锐行为。

本例要学的: 分段由圆盘与初值支集的交叠关系决定;先求积分再求时间导数。无需背这个特定单位圆盘的答案。

5.12 · 教材§5.4.2,第362—365页 · 自编检查题

练习23 二维核的总积分

用极坐标验证 。

展开解析

角积分 消掉前系数,剩 。可令 直接积分。边界奇异虽发散,整体可积。

5.12 · 教材§5.4.2,第362—365页 · 自编检查题

练习24 二维二次初位移

二维波动方程 、。不做圆盘积分,构造解并验证。

展开解析

设 。则 ,,所以 。解为 。三维对应 ,由维数决定。

数学物理方程篇·第五章 / 5.13

传播范围与有无后效

教材带星号§5.4.3,第365—370页

5.13 · 教材带星号§5.4.3,第365—370页

46 影响区域和依赖区域是两个观察方向

影响区域: 固定一个源点,问它在时刻 能影响哪些观察点。依赖区域: 固定一个观察点和时刻,问它依赖哪些初始点。对平移不变的自由波,二者使用同一距离条件,但看图的中心不同。

一维、二维的速度核分别支撑在区间 和圆盘 ;三维速度核支撑在球面 。三者都有有限传播速度 ,因为核在 为0。有限传播不等于只有波前有响应。

5.13 · 教材带星号§5.4.3,第365—370页

47 三维紧支集信号何时开始、何时结束

假设初位移和初速度都集中在有界区域 ,且没有外力、边界反射。固定观察点 ,令到 的最近、最远距离为 。

当 ,依赖球面还没碰到支集,响应为0;当 ,球面已完全越过支集,响应重新为0。中间时段有可能有响应,但也可能由于数据抵消而为0,不能把“可能非零”说成“必非零”。

这就是三维自由波的强惠更斯性质:点源响应在波前,紧支集信号经过后没有持续的内部尾部。这里说的是标准常系数、全空间、无外力的三维波方程;加边界、阻尼、势项后不能原封不动推广。

5.13 · 教材带星号§5.4.3,第365—370页

48 二维为何有尾巴,一维又要看哪种初值

二维依赖域是整个圆盘。圆盘覆盖了初值支集以后,仍然包含它,所以一般继续有衰减响应;有明确前沿,却一般没有明确结束时刻。这叫后效。并不是说每一组特殊初值都必定产生非零尾巴。

一维达朗贝尔公式中,初位移项只取两个端点,紧支集初位移会分成左右移动的波;初速度项则在整个区间积分。例如紧支集初速度的总积分非零,固定观察点在区间完全包住支集后会留下常数位移。这说明一维也不能笼统套三维“波后恢复零”的结论。

图 P5-22 三维紧支集自由波有结束窗,二维一般留尾部
只讨论全空间标准自由波、紧支集初值;几何相交表示可能非零,数据抵消仍可能使解为零。三维:依赖球面穿过支集后,信号结束观察点初值支集距离时间窗dₘᵢₙ/a:开始可能到达dₘₐₓ/a:最后可能到达当前外球面已越过小支集,所有球面数据为0。

只讨论全空间标准自由波、紧支集初值;几何相交表示可能非零,数据抵消仍可能使解为零。

为什么这部分要学: 能不算积分就判断“还没到”“已经过去”“可能有尾巴”,直接解决传播时间和影响域题。教材虽标星,图像和结论仍值得第一轮掌握;完整几何证明可第二轮再看。

5.13 · 教材带星号§5.4.3,第365—370页 · 自编检查题

练习25 三维信号时间窗

源初值支集为半径2的球,球心到观察点距离5,波速 。在没有外力和边界时,三维自由波在哪个时间窗内可能非零?

展开解析

最近距离3、最远距离7,所以可能非零的时间窗为 。窗外必为0;窗内是否真的非零取决于数据。前后临界时刻若有不光滑初值,需要按对应解的意义解释。

5.13 · 教材带星号§5.4.3,第365—370页 · 自编检查题

练习26 同一个支集换成二维

把上题改成二维圆盘,波速仍为3。是否能断言 后解恒为0?

展开解析

不能。 时依赖圆盘与支集不相交,必为0; 后圆盘覆盖支集,但积分仍在支集上进行,一般存在尾部。只凭距离无法断言尾部系数非零或决定其符号。

数学物理方程篇·第五章 / 5.14

推迟势是空间距离换成时间延迟

教材§5.4.4,第371—372页

5.14 · 教材§5.4.4,第371—372页

49 从时间积分推到推迟时间

三维零初值、源为 的波动方程有

固定源点 ,δ在时间上的根为 。根在 内要求 ;对 的导数绝对值为 ,再产生一个 。于是

这里的两个 来自不同地方:一个原本就在三维速度核中;另一个来自对时间抽样的雅可比。把它们都交代清楚,分母的 就不容易背错。

图 P5-23 只把过去光锥上的源送到现在的观察点
图用距离代替全部空间方向,是过去波锥的半个径向截面。源在负时间不存在时,还要限制r<at。观察点的过去锥:τ=t−r/a时间距离r现在:r=0,时刻tτ=t−r/aτ=t−r/aτ=t−r/a距离越远,取的发出时刻越早;还要乘1/(4πa²r)。

图用距离代替全部空间方向,是过去波锥的半个径向截面。源在负时间不存在时,还要限制r<at。

5.14 · 教材§5.4.4,第371—372页

50 只取发出时刻的源,不能取到达时刻的源

距观察点 的源需要传播 的时间,所以现在 看到的是过去 的源。若源只在 开始作用,未到达的点要排除。也可将 对负时间延拓为0,然后在全空间积分。

例如点源 ,观察点 ,令 。当 时

当 时为0。源点本身需要分布或奇异解解释,不能在上式直接令 。若 为有限时长信号,则到达观察点的是同样时长、平移了 的信号,幅度含 。

5.14 · 教材§5.4.4,第371—372页

51 其他维数的外力不能照搬三维推迟势

二维零初值源响应为

二维核在圆盘内都有值,没有只在单一传播时间上抽样的δ,所以一般不能缩成三维那样的单个推迟时刻。要保留历史积分。

核对全章: 源项永远是“单位源响应 × 实际源强,再把允许的时空点加起来”。静态只积空间,动态还要处理时间;边界、初始位移、初速度分别有自己的入口。

5.14 · 教材§5.4.4,第371—372页 · 自编检查题

练习27 点源到达时间

波速2,点源与观察点相距6,源强在发出时刻 为 ,其余为0。观察点何时可能收到信号,幅度是多少?

展开解析

传播时间为3,所以接收区间为 。幅度为 。区间端点按信号阶跃的约定解释。不能把发出区间 当成接收区间。

5.14 · 教材§5.4.4,第371—372页 · 自编检查题

练习28 近处和远处看见的是不同历史

同一源强 从0时开始。两个观察点距源分别为 与 。它们在同一时刻看见的相位和幅度关系是什么?

展开解析

都已到达以后,两者分别取 与 。远点的幅度因子减半,并多延迟 ,相位差为 。不能只比较 而忽略推迟时间。

数学物理方程篇·第五章 / 5.15

原书章末12个大题逐项解题

教材习题,第372—374页

5.15 · 教材习题,第372—374页

52 原书题1:五个δ恒等式都用测试函数证明

以下默认测试函数光滑且有紧支集,所有恒等式均在分布意义下成立。

(1) :与 配对得 。

(2) :两边与 配对都等于 。

(3) ,:换元 ,两边对 的作用都为 。负的 会反向积分区间,结果仍含绝对值。

(4) :由 ,而 。

(5) :左边对 的作用是 。

评分点: 不能只写“由δ在0以外为0可知”,那不足以证明分布恒等式。写出任意测试函数的积分才是完整的一步。

5.15 · 教材习题,第372—374页

53 原书题2:坐标变换中的δ

设 在对应点附近光滑可逆,雅可比 。待证

证明时将右边乘 ,对 积分,δ把它取在 ,得到 。这正是左边与 的配对,所以成立。教材写 时,δ会自动将其取在源点。

极坐标 ,故 时为 。角变量按周期或单一坐标图处理,避免 与 重复。源在极点时用径向分布处理,不能除以 。

5.15 · 教材习题,第372—374页

54 原书题3:有限杆点初值与有限弦点初速度

(1) 热方程,零值边界。 取 ,范数平方 ,初值δ的系数为 ,所以

每项在两端为0, 的级数初值按分布意义恢复δ。对 指数衰减使它成为平滑温度。

(2) 波动方程,零导数边界。 空间应为余弦,且包含常数模态。令 ,得到

在 为0;时间导数的初值是余弦δ展开。对 求导后每项带 ,两端均为0。空间积分等于 ,也与单位总初速度一致。

原书答案勘误说明: 扫描本印刷第398页对应答案的非零模态含 ,与第373页题面 不匹配。这里按题面给出余弦解;不能同时保留 又把空间基写成正弦。端点求导检验比机械对答案可靠。

5.15 · 教材习题,第372—374页

55 原书题4:各向异性与双调和基本解

本题统一取 ,不是 。

(1) ,。令 ,算子变成二维拉普拉斯,而 ,故可取

它与 只差常数。必须出现 ,因为δ也随坐标缩放。

(2) 二维 。对 ,,再用 ,可取

第一步没有隐藏的点δ: 的小圆通量趋于0。第二步才由 产生点源。

(3) 三维 。因 (分布意义下此步无点δ),而 ,所以

检查符号:再施加一次拉普拉斯,。

5.15 · 教材习题,第372—374页

56 原书题5:亥姆霍兹基本解还需说明选择

对 ,取 。空间傅里叶变换形式上给 。分母在球面 为0,需要指定分布解释或无穷远条件,不能声称其普通逆积分绝对收敛。

不另加辐射条件时,可取实基本解

验证: 时径向方程 ; 的主奇异项为 ,通量为1,零阶项的小球积分趋于0,所以源为 。

若采用时间因子 并要求向外辐射,可选 。它与实核的差是齐次解。本题必会: 频域分母、近源归一化、为何基本解不是不加条件就唯一。辐射条件的完整谱论推导可暂缓。

5.15 · 教材习题,第372—374页

57 原书题6:三种三维区域的镜像格林函数

本题取 、零值边界。

(1) 四分之一空间 。 源 ,另设 、、。则

反射一次换号,反射两次恢复正号。在任一坐标平面上成对等距抵消,三个镜像均在区域外。

(2) 上半球 。 先用球内格林函数 处理球面,再令 处理平面:

在球面上两项各自为0,在 上两项相等。展开后有真实源、反演源、平面反射源、反演后反射源四项;只有真实源在上半球内部。球面反演权重仍为 ,不能漏掉。

(3) 层状空间 。 真实源 ,。平移像点 与 ,则

按每个括号成对求和。 处通过 等距抵消, 处通过 抵消。一次反射无法同时满足两面,故需要无限像列。无穷远方向取符合衰减要求的解。

5.15 · 教材习题,第372—374页

58 原书题7:二维象限与半圆盘

本题仍取 。二维点源势改用 ,而不是三维的 。

(1) 第一象限。 真实源为 ,镜像为 、、,相应距离为 。反射符号与三维相同:

在 上 ,在 上 ;两处均为0。

(2) 单位上半圆盘。 用单位圆盘格林函数 ,源 、、,平面反射为 。取

若 、,写开为

两个圆盘核的常数修正相消。圆弧上各自为0,直径上相等相减。解法思想: 先满足曲边,再用反射补平边;不必重新猜四个源的权重。

5.15 · 教材习题,第372—374页

59 原书题8:带反应的热方程

方程 ,基本解的初值为δ。令 ,则 ,。因此

检验总积分为 ,并满足 、。不要以热方程的守恒性否定这里的增长或衰减;多出的 本来就在产生或消耗总量。

5.15 · 教材习题,第372—374页

60 原书题9:输运与阻尼波的基本解

(1) ,初值 。 沿特征线 ,全导数为 。因此

相应初值核为 。必须核对漂移方向:当 ,形状 向右移。

(2) ,初位移0、初速度 。令 ,逐项求导后得到

频域速度核为 。其逆变换为

所以

教材提示中的一般参数式需要成套处理缩放:逆变换 对应 。在本题 ,于是得到上式。第一轮无需证明整张贝塞尔积分变换表,但必须会指数替换、参数对应与初值检查。

5.15 · 教材习题,第372—374页

61 原书题10:二维非齐次波动的积分式

初位移、初速度都为0,所以只保留源项。令 ,则

得分要件: 正确的 、正确的历史时间 、随 缩小的圆盘、正确的根号。用常数源 检验:空间积分为 ,时间积分为 。

5.15 · 教材习题,第372—374页

62 原书题11:再降一维得到达朗贝尔

把二维数据取为与 无关。对二维核中的横向变量积分,固定 ,有

当 成立,外部为0。这正是一维平台核。再对初位移核卷积求时间导数,就得到

这里降维做了真正的横向积分;只删变量会丢掉几何权重,得不到正确的一维系数。

5.15 · 教材习题,第372—374页

63 原书题12:四道直接使用公式的题

(1) 高斯初温。 对 、,第33单元的 给出

(2) 二维波,初位移 ,初速度0。 因 、,时间展开在二次项终止:

检验 ,两个初值也正确。虽然多项式不在 ,核在有限圆盘上的积分仍有意义,且可直接验证所得经典解。

(3) 二维波,源 ,初位移0、初速度 。 因 ,设 ,得到 ,故

(4) 三维波,源 ,初位移和初速度都为 。 同样 ,所以 ,得到

这组题真正考什么: 公式的三份数据是否放对位置,以及能否看出拉普拉斯使多项式降次。能用有限模态或多项式直接求解时,不要强迫自己做复杂积分。

数学物理方程篇·第五章 / 5.16

最后八题与全书学习收尾

本章综合训练与掌握标准

5.16 · 本章综合训练与掌握标准

64 把整章压成一张方法选择表

题目触发条件 先求什么 最后怎样组装
全空间 再加满足附加条件的齐次解
、给定边界值 、 源体积分减边界法向导数积分
一阶时间 初值卷积加历史源积分
二阶时间 位移项求时间导数,速度项直接卷积,再加源
平面、球、圆等简单边界 镜像或反演修正 先检查所有边界,再写数据积分
矩形等有限区域 特征函数和范数 静态除以特征值,动态乘相应时间响应
图 P5-24 先看区域与时间阶数,再决定解的组装方式
方法由题目数据决定;算完后检验条件,是比只看答案形状更可靠的习惯。先辨认数据和区域,再选表示式静态 + 有边界源:Δu=−f边界:u=gG的体积分与边界积分动态 + 全空间时间一阶:初值与外力时间二阶:位移、速度、外力选对维数的响应核简单边界平面反射、球面反演先核查像点位置和符号有限区域用特征基展开δ范数、零模态不能省最后固定做:方程、初值、边界、归一化或通量四项检查。

方法由题目数据决定;算完后检验条件,是比只看答案形状更可靠的习惯。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习29 δ乘积不要丢掉导数项

化简 。

展开解析

按 ,结果为 。用任意 检查作用为 。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习30 给定核后恢复一个静态解

单位球内 ,边界值为2,且解在球心正则。求解并说明两份数据各贡献多少。

展开解析

源的零边值响应为 ;边界常数的调和响应为2。总解 。拉普拉斯为 , 时为2。这两份分别对应格林源积分与泊松边界积分。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习31 半轴热核的总量

初温 ,。半轴零温边界与绝热边界在很长时间后的总热量分别怎样?说明理由。

展开解析

绝热核 在半轴积分始终为1。零温核相减,总积分为 ,随 趋于0。后者可将两个高斯积分换元得到,反映边界吸收热量。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习32 三维球面平均的混合项

三维自由波动方程初位移 ,初速度 。求解。

展开解析

为调和多项式,球面平均不变;常数初速度给出 。解为 。时间二阶导数和空间拉普拉斯都为0,初值正确。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习33 二维波动与恒定源

二维全空间 ,初位移 ,初速度0。求解。

展开解析

初位移部分为 ,源部分为 ,所以总解为 。直接代入 检查。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习34 推迟势中的δ缩放

说明 对 积分时为什么产生 。若 ,结果为何为0?

展开解析

内函数对 的导数是 ,根 。简单根公式给 。 时根为负,不在 内,所以没有抽样贡献。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习35 三维能否出现永久尾部

有人说:“三维点源响应只在球面,所以任何三维有外力问题的解最后必为0。”指出缺少的条件。

展开解析

持续外力会不断发出新信号,不能从单次脉冲的球面支撑推断总响应必终止。恢复零的时间窗结论还依赖无外力、紧支集初值、全空间标准自由波、无边界反射等条件;无界初值也可能一直提供球面数据。

5.16 · 本章综合训练与掌握标准 · 自编检查题

练习36 独立识别三类题的第一步

分别写第一步:A. 有限绝热杆中点加热;B. 半球内零边界点源;C. 二维全空间零初值、给定随时变化外力的波动。不要计算最终积分。

展开解析

A. 用余弦特征函数含常数模态,将δ投影,时间乘热衰减因子。B. 用球内格林函数减去源点关于平面的反射核,并核查球面和平面均为0。C. 写二维圆盘波核的时空卷积,空间半径为 ,再对 积分。三题分别由边界、几何、时间阶数决定入口。

5.16 · 本章综合训练与掌握标准

65 到什么程度算学会本章

基础过关: 不看答案做对检查题1—8,能说出δ是如何抽样的;面对基本解先写它满足 还是 。

标准习题过关: 能独立写出半空间与球内格林函数;能正确组装源项、边界项和初值项;能完成原书题3、8、9(1)、10、12,最后逐项验初值和边界。错因要记录为具体动作,例如“忘了余弦的常数模态”,不要只写“粗心”。

提高过关: 能推导二维降维核和三维推迟势;能处理复合镜像、双调和基本解与亥姆霍兹核;能交代正则性、支撑、无穷远条件和辐射选择。原书题4—7及9(2)适合放在第二轮集中做。

不要求第一轮背下所有核的每个积分推导。要求你能用一个短检验重建关键系数:静态核用小球通量,热核用总积分1,波动速度核用总积分 ,边界核用代值或求法向导数。

5.16 · 本章综合训练与掌握标准

66 全书整理完以后,下一步怎样学

这本书复变函数篇七章、数学物理方程篇五章的正文材料至此均已整理成讲义。学习可以分三次往回走:

  1. 方法骨架: 复积分与留数 → 傅里叶/拉普拉斯变换 → 方程、初值、边界分类。每一章先做最基本的例题,能说出为什么选该方法。
  2. 方程主线: 建模与定解条件 → 分离变量及正交投影 → 特殊函数 → 积分变换 → 本章基本解。重点检查边界如何筛选特征函数,以及初值如何决定系数。
  3. 混合练习: 把题目抄到另一页,遮住章标题,先判断区域、边界、时间阶数、源和初值,再求解。若离开章标题就选不出方法,优先练判断,不要继续背公式。

每次做题都固定写五行:题目数据 → 选法理由 → 核或特征基 → 系数/积分 → 验证。 做错后定位具体失败的那一行,回看对应单元,再用一题同类小题确认已经改正。

本讲义附的36道检查题都有解析,但应先遮住答案独立写。原书本章12个大题也给了路线、推导要点和结果;这不能代替你亲自完成一次计算。真正能让脑子清楚的,是把每个符号对应到“它在题目里负责什么”。

来源、编排与掌握要求

依据用户提供的严镇军《数学物理方法》扫描本数学物理方程篇·第五章整理;正文印刷页325—372上部,章末习题372—374。扫描PDF文件页序等于印刷页码加6。按原书节序覆盖各主题,以重新讲解、推导主线和例题说明组织,采用原创讲解,不逐字复述原文。自编检查题与原书习题分开编号,原书题号可在本章习题路线中查找;第五章另按原书12个大题逐项给出解法、结果与适用条件;自编检查题另列。

第一轮优先掌握能独立做题所需的概念、步骤和适用条件;背景、证明及拓展均说明用途和暂缓深究的边界。数学物理方程篇·第五章含36道自编检查题及完整解析。公式与SVG内嵌,可离线阅读;屏幕答案默认折叠,打印时自动展开。整理日期:2026年10月22日。