电磁场课后习题验证指南:从麦克斯韦方程组到Python复核

发布时间:2026/9/19 14:20:52
电磁场课后习题验证指南:从麦克斯韦方程组到Python复核 简介这是一份面向电磁场理论课程学习者与备考者的课后习题答案资料内容围绕导行电磁波展开系统梳理纵向场法、TE/TM/TEM波判定、矩形波导模式与截止波长、同轴传输线场分布、圆柱介质波导与光纤、谐振腔等核心知识点。每章配套小结便于快速回顾公式与结论题目解答则能帮助验证理解、突破计算难点。压缩包内共1个PDF文件大小1.74MB便于下载后按章节查阅。目前已有832人学习下载适合电磁场、微波工程及通信技术相关专业学生自学与考前复习使用。1. 电磁场课后习题答案只是起点可验证才算学会电磁场课后习题和普通电路题最大的区别在于你很难从结果直接看出对错。算出一个电场强度是 3.2 kV/m 还是 3.2 V/m推导过程里每一步却都可能“看起来是对的”。这也解释了为什么“答案 每章小总结”这类资料比单纯给数值更有用——对答案时真正应该对照的是每章把知识压成的检查框架。下文从一线工程师视角把电磁场习题的常见题型、数值复核方法和参数陷阱拆开讲覆盖正在修课的学生、转做射频天线的软件工程师以及需要回炉电磁场理论的硬件开发者。目标只有一个拿到任何一份习题答案都能独立判断它是否成立。2. 电磁场习题的题型地图从麦克斯韦方程组到三类必考计算2.1 麦克斯韦方程组是习题的唯一主线麦克斯韦方程组是电磁场课后习题的唯一主线。静电场、恒定磁场、时变电磁场、传输线、波导和辐射全部是这一组方程在不同前置条件下的特例。教材章节顺序几乎都按这条主线展开所以“每章小总结”的实质就是重新强调当前章节对应的那组简化前提。静电场只要求 ∇×E 0 和 ∇·D ρ电场无旋意味着标量电位 φ 存在问题就转化为解拉普拉斯方程或泊松方程到了时变场∇×E -∂B/∂t 让电场不再无旋电位方法必须升级为矢量磁位与标量电位的组合。理解这个递进关系比背公式重要因为习题的难点从来不是代入数值而是判断“当前场景允许用哪一层简化”。典型错误是把恒定电流场直接套用静电场的高斯定理忽略了 ∇·J 0 的电流连续性对场分布的限制结果在导体分界面得出错误的电流密度。写章节总结时我会要求自己先用一句话回答“本章在什么条件下成立”再列公式没有这句前提的总结检索价值会大打折扣。2.2 边界条件四条关系放倒大多数分界面题边界条件由麦克斯韦方程组的积分形式导出但课后习题里它们以“分界面两侧场量关系”出现。四条边界条件、适用场景和常见失分点可以整理成一张表边界条件表达式适用场景常见失分点切向电场连续n̂ × (E₁ - E₂) 0任意两种介质分界面忘记导体内 E 0法向电位移跳变n̂ · (D₁ - D₂) ρs分界面存在自由面电荷默认 D 法向连续切向磁场跳变n̂ × (H₁ - H₂) Js分界面存在面电流面电流方向判反法向磁感应连续n̂ · (B₁ - B₂) 0所有磁介质分界面误写为 H 法向连续做题时的判断顺序我一般固定为先看两侧介质类型再看分界面上有没有自由面电荷或面电流最后决定哪条条件能用、哪条要带源项。很多答案对不上就是把“无源分界面”和“有源分界面”混用了。比如导体-介质分界面上导体内电场为零介质一侧的切向电场也必须为零而法向电场由面电荷密度决定数值可以很大这个结论在传输线内导体和波导壁上反复出现。用一个很小的 Python 片段验证无源分界面上的两组条件import numpy as np # 介质1 和 介质2 中的电场界面法向为 x 方向 E1 np.array([10.0, 3.0, 0.0]) # V/m E2 np.array([ 2.5, 3.0, 0.0]) eps1, eps2 1.0, 4.0 # 相对介电常数 # 切向分量y、z必须相等E 切向连续 assert abs(E1[1] - E2[1]) 1e-9 assert abs(E1[2] - E2[2]) 1e-9 # 法向分量必须满足 D1n D2n无自由面电荷 D1n eps1 * E1[0] D2n eps2 * E2[0] print(fD1n{D1n:.2f}, D2n{D2n:.2f}, 残差{abs(D1n-D2n):.2e})代码逻辑是把题目给的两组场量代进边界条件做“对账”前两个断言检查切向电场连续后两行检查法向电位移连续。实际使用时把题目算出来的 E1、E2 替换进去残差接近 1e-9 就说明这一步自洽不满足答案十有八九是在分界面处理上出了问题。参数上注意 eps 取相对介电常数即可因为真空介电常数 ε0 在两侧相同时会被约掉所以代码里不出现。2.3 带每章小总结的答案应该倒着读拿到一份带每章小总结的电磁场课后习题答案不要急着翻题目我一般先做三个动作。第一通读小结把每个公式的限定条件圈出来第二从本章挑 3 到 5 道代表性习题独立做一遍卡住超过十五分钟再参考答案第三核对时只看推导链条直到最后一步才看数值。这套顺序针对的是习题答案最常见的错误形态中间步骤的符号或方向错了最终数值只差一个系数很容易被当成粗心而漏过。提示小结里“注意当……时……不成立”这类句子才是复习时的核心资产。公式可以现查条件判断只能靠平时积累。3. 用 Python 复核电磁场课后习题答案数值验证集与数量级巡检3.1 最小依赖numpy matplotlib 覆盖九成静场题复核电磁场习题答案不需要重型仿真软件常见做法是只装 numpy 和 matplotlib把答案给的场分量代进去做数值检查。微波频段的复杂结构再考虑 CST 或 HFSS但课后题几乎都可以用解析或半解析方法验证。最基本的检查是旋度和散度静电场必须无旋恒定磁场必须无散这两条性质不受坐标系选择影响非常适合程序化验证。import numpy as np # 给定二维静电场 E (2x, -2y)验证无旋条件 # 无旋要求 curl(E) dEy/dx - dEx/dy 0 x, y np.meshgrid(np.linspace(-1, 1, 201), np.linspace(-1, 1, 201)) Ex 2.0 * x Ey -2.0 * y curl_z np.gradient(Ey, x[1, 0] - x[0, 0], axis1) \ - np.gradient(Ex, y[1, 0] - y[0, 0], axis0) print(max|curl(E)| , np.abs(curl_z).max())逻辑说明np.gradient默认用二阶中心差分近似偏导axis1表示沿 x 方向axis0表示沿 y 方向网格步长必须显式传入否则默认步长为 1会把导数整体放大。输出接近 1e-10 说明场无旋电位函数存在数值大于 1e-3 则说明该场不满足静电场条件答案需要重推。把Ex、Ey换成任意题目给出的场分量这套骨架可以直接复用。验证磁场无散同理改算np.gradient(Bx) np.gradient(By)看散度是否为零。3.2 传输线反射系数的复平面自洽校验import cmath Z0 50.0 # 特征阻抗单位 ohm ZL 75.0 25.0j # 负载阻抗单位 ohm Gamma (ZL - Z0) / (ZL Z0) rho abs(Gamma) vswr (1 rho) / (1 - rho) print(fGamma {Gamma:.4f}) print(f|Gamma| {rho:.4f}, VSWR {vswr:.4f}) assert rho 1.0 1e-12, 无源负载的反射系数模值不可能大于 1参数说明Z0 由传输线几何与填充介质决定同轴线 Z0 (60/√εr)·ln(b/a)其中 b 是外导体半径、a 是内导体半径εr 是填充介质的相对介电常数ZL 是终端负载题目没给复数形式时实部也会走同一套公式。复核分两步走先看模值无源负载下 |Γ| 必须小于等于 1否则公式里 ZL - Z0 或 ZL Z0 的符号、共轭处理有误再看三个锚点ZL Z0 时 Γ 0、VSWR 1ZL 开路时 Γ → 1ZL 短路时 Γ → -1。这三个锚点能一次排除反射系数公式的大部分笔误比把 VSWR 算到小数点后四位更可靠。3.3 数量级巡检表三秒排除离谱答案电磁场习题答案往往差在数量级而不是数值本身。我复核时先对数量级再对符号最后才核对精确值。一份常用的数量级巡检表检查对象合理区间说明空气介质电场强度1 V/m ~ 1 MV/m超过 3 MV/m 达到空气击穿自由空间波阻抗376.73 Ω看到 377 直接对应 η0同轴/微带特征阻抗10 ~ 150 Ω教科书习题极少超出波导截止频率1 ~ 100 GHz与波导尺寸直接相关天线远场距离r ≥ 2D²/λ小于该值必须用近场方法这张表不用背放在手边即可。很多答案看起来算到了小数点后三位实际因为 ε0 和 μ0 的幂次弄反结果差了 10¹² 倍。数量级这一关过掉之前不值得去核对具体数值。3.4 复杂化简用 SymPy 兜底遇到分式化简、三角恒等式和带复数模值的推导手算容易滑档。常见做法是让符号计算顶一遍import sympy as sp c, t, d sp.symbols(c t d, positiveTrue) expr sp.sqrt((c*t - d)**2 1) / (1 c*t*d) print(sp.simplify(expr))说明把习题推导里要化简的表达式抄进expr用sp.simplify或sp.expand验证两个中间表达式是否恒等。positiveTrue这个假设对开根号和绝对值化简很关键不加它 SymPy 会保留大量分支结果看起来“不乱但用不上”。符号计算适合验证代数步骤不适合替你做物理判断边界条件和适用场景仍然要自己确认。4. 四类高频电磁场计算题的拆解模板参数、量纲与边界条件4.1 静电场高斯定理与泊松方程的选择判据静电场题的第一判断是体系对称性。无限长线电荷、无限大平面电荷、球对称分布直接用 ∮E·dS Q/ε0 会比解偏微分方程快一个量级两个偏心导体球、非均匀介质中的电荷对称性被破坏必须退到泊松方程 ∇²φ -ρ/ε 加边界条件数值求解。import numpy as np N 200 # 离散网格数 L 1.0 # 求解区间长度m V0, V1 0.0, 10.0 # 左右边界电位V rho 1e-9 # 体电荷密度C/m^3 eps0 8.854e-12 # 真空介电常数 h L / (N - 1) A np.zeros((N, N)) b np.zeros(N) for i in range(1, N - 1): A[i, i - 1] 1 / h**2 A[i, i] -2 / h**2 A[i, i 1] 1 / h**2 b[i] -rho / eps0 A[0, 0] A[-1, -1] 1 b[0], b[-1] V0, V1 V np.linalg.solve(A, b) E -np.gradient(V, h) print(fV_max {V.max():.3f} V, E_max {np.abs(E).max():.3f} V/m)这块代码对应平行板电容器间含均匀体电荷的典型题。参数说明N取 200 时二阶中心差分在光滑问题上误差在 1e-4 量级再增大网格对答案精度提升有限边界行单独赋单位值保证系数矩阵非奇异。np.gradient(V, h)的第二个参数是网格步长漏掉它电场就会整体差一个因子1/h。验证解的正确性把V代回差分方程看残差再把E做散度应还原为rho/eps0。4.2 传输线特征阻抗公式的适用边界传输线课后题的高频参数是特征阻抗 Z0但同一个符号在不同几何下公式完全不同。同轴线内是 TEM 模Z0 (60/√εr)·ln(b/a)微带线上方空气、下方介质要用有效介电常数 εeff 修正场的非均匀性。微带线的常用解析公式Pozar 教材形式import math def microstrip_z0(w, h, er): u w / h # 有效介电常数数值介于 1 和 er 之间 eps_eff (er 1) / 2 (er - 1) / 2 / math.sqrt(1 12 / u) if u 1: z0 60 / math.sqrt(eps_eff) * math.log(8 / u u / 4) else: z0 120 * math.pi / (math.sqrt(eps_eff) * (u 1.393 0.667 * math.log(u 1.444))) return z0, eps_eff # FR4 常用参数er4.4, h1.6mm, w3.0mm z0, ee microstrip_z0(3.0, 1.6, 4.4) print(fZ0 {z0:.1f} ohm, eps_eff {ee:.3f})参数陷阱有三个。第一u w/h决定走哪个分支u 小于等于 1 和大于 1 的公式不能互换混用会让对数项出现符号错误第二εeff 永远在 1 和 εr 之间算出超出这个区间的值说明公式抄错第三铜箔厚度对窄线w/h 较小时才明显常规 FR4 板可以忽略但高频习题里如果明确给了厚度参数要先用等效线宽修正再做计算。算出结果后和上一章 3.2 节的反射系数校验连起来用整个传输线大题的自洽性就闭环了。4.3 矩形波导截止频率与模式编号的对应矩形波导的截止频率公式是 fc_mn (1/(2√με))·√((m/a)² (n/b)²)。空气填充时 √με 1/c主模 TE10 的截止波长是 2a。课后题最常见的失分点是 a、b 对调波导约定 a ≥ bTE10 是唯一能单模工作的模式工作频率必须落在 fc(TE10) 和 fc(TE20)或 fc(TE01)之间。模式场结构截止频率典型用途TE10宽边方向半个驻波c/(2a)矩形波导主模TE20宽边方向一个整驻波c/a高次模抑制TE01窄边方向半个驻波c/(2b)低损耗长距离传输TM11纵向电场和磁场并存(c/2)·√(1/a² 1/b²)天线馈源解题时先辨识模式下标再代公式m 对应宽边 a 方向的半驻波数n 对应窄边 b 方向。题里出现“TE10 模式单模传输”本质上是在给你宽边 a 和频率区间两个约束反推波导尺寸时直接用不等式联立。4.4 辐射与天线远场距离这道隐形前置辐射类的课后习题答案很多错在适用范围拿了远场公式代入的却是近场距离。工程上远场条件取 r ≥ 2D²/λD 是天线最大尺寸。小于这个距离球面波不能近似成平面波弗里斯公式 Pr Pt·Gt·Gr·(λ/(4πr))² 不成立。验证代码如下import math f 2.4e9 # 频率Hz wl 3e8 / f # 波长m D 0.3 # 天线最大尺寸m r 10.0 # 接收距离m r_far 2 * D**2 / wl print(f远场距离 {r_far:.2f} m, 当前距离 {r:.2f} m) assert r r_far, 当前距离不满足远场条件弗里斯公式不适用 # 电尺寸检查D 远小于波长时应按电小天线处理 print(fD/lambda {D / wl:.2f})参数说明f 单位取 Hz波长用光速除以频率得到D 和 r 保持同一单位。第一段输出直接对比当前距离和远场距离断言失败意味着答案必须改用近场或数值方法。第二行电尺寸检查是补充分支当 D/λ 小于约 0.1 时天线接近电小天线增益模型要换这也是雷达方程和链路预算题里常见的隐藏扣分点。5. 把每章小总结做成可检索 PDF 的四个结构技巧做这类资料时我最后才写每章小总结不是因为它不重要而是想清楚全章的题型和易错点之后总结才写得准。下面四个结构技巧是从检索和复习两个用途倒推出来的。5.1 公式索引表优先于正文段落章节总结最常见的失败形态是写成叙述体流水账比如“本章介绍了静电场的基本方程”这类文字在检索时几乎没有命中价值。我的做法是只留三类内容公式、限定条件、题型判断。开头放一张三列索引表公式、物理量说明、适用条件。复习时按条件反查公式比从正文里翻快一个量级。这张表也是 PDF 全文搜索命中率最高的区域值得把最常用的四条边界条件和三大电位方程放进去。5.2 条件限定句单独成块麦克斯韦方程组只有四个方程但工程场景的简化版本非常多。“当频率远低于……”“当介质无耗……”“当结构尺寸远小于波长……”这些限定句才是分界线。把每章出现的限定条件抽出来用 Markdown 的引用块或 PDF 的提示框单独呈现检索时一眼可见。这种写法的附加好处是逼自己回答“这个公式什么时候不能用”比把公式抄十遍都管用。5.3 题型-方法映射表收尾每章末尾加两列映射表左边题目特征右边解法入口。例如“球对称电荷分布 → 高斯定理直接积分”“平行板间含介质分层 → 边界条件逐层递推”“矩形波导求单模带宽 → 截止频率不等式联立”。这张表等于把整章的解题流程压成一张决策树也是后续复习时定位原文的锚点。题目特征描述要写题干里真正出现的话术而不是自己总结的抽象名词这样检索命中率才高。5.4 用 Pandoc 生成带书签的 PDFpandoc chapter2.md \ --pdf-enginexelatex \ -V CJKmainfontNoto Serif CJK SC \ --toc --toc-depth2 \ -o chapter2.pdf参数说明--toc生成可展开的书签目录--toc-depth2控制书签到二级标题中文字体参数避免 xelatex 渲染中文报错。生成后建议用pdfinfo chapter2.pdf检查元数据再用pdftk chapter2.pdf dump_data确认书签完整。章号和节号直接写进标题文本这样即使脱离 PDF 书签按文件名和书签文本也能定位到具体的电磁场课后习题章节。本文还有配套的精品资源点击获取