简介《介质界面极化电磁波的传播特性研究》是一篇学术论文围绕电磁波在介质界面的极化行为系统分析了垂直极化波与平行极化波的反射、透射系数及其随介质电容率的变化规律。资料源自期刊论文适合通信、雷达及电磁场专业的高年级学生与工程技术人员参考可用于理解全反射、无反射等现象的理论推导与数值模拟。压缩包内为1个PDF文件容量249KB内容结构完整含边界条件推导、反射透射系数公式及模拟曲线便于按需查阅或打印学习。目前已有94人学习该资源。通过阅读可掌握介质电容率对反射与透射系数的影响规律了解斜滑投射及光密/光疏介质条件下的传播特性差异为极化控制信息传播技术或相关反隐身、偏振分析应用提供理论支持。1. 介质界面极化电磁波为什么先读这一篇再碰极化匹配聊极化匹配之前先说一个反直觉的结论地面雷达之所以存在低空盲区很大一部分原因就藏在介质界面极化电磁波的反射规律里。斜滑投射时不论垂直极化还是平行极化反射系数都趋近 -1反射波与入射波相位相反直达波和地面反射波在天线处直接相消。这是 2014 年《内蒙古师范大学学报自然科学汉文版》一篇论文末尾给出的解释而整篇论文做的就是把介质界面上这两种极化波的反射系数、透射系数用电磁场边界条件推出来再给出一组可复现的数值模拟结果。对做电磁场课程设计、天线极化匹配、雷达链路仿真的人来说这份 PDF 最有价值的地方在于公式给得直接参数只有电容率比值 ε₂/ε₁场景又覆盖了全反射、无反射、斜滑投射这些最容易概念混淆的边界情况。你可以把它当作一份能落地到代码里的理论依据而不是一篇看过就忘的论文。下面我按“公式怎么来、结论怎么用、代码怎么写、坑在哪”的顺序把这四页纸拆开讲透。2. 从边界条件到反射透射系数四个公式决定你的曲线走向2.1 垂直极化与平行极化的定义先分清看的是哪个电场分量任意极化的平面电磁波都可以分解成两个正交分量电场方向垂直于入射面的叫垂直极化波电场方向在入射面内的叫平行极化波。这里的“入射面”由入射波矢和界面法线共同确定不是一个固定的水平面或垂直面。很多做天线仿真的同事第一次看这篇论文时都在这翻车——把垂直极化理解成“电场垂直于地面”但论文里的垂直极化是相对入射面定义的两者只在特殊场景下重合。在非磁性介质假设下μ₁≈μ₂≈μ₀界面两侧的磁导率差异被消掉最后决定反射、透射行为的只剩电容率的比值 ε₂/ε₁。电容率 ε 在中文教材里常叫介电常数英文习惯叫 permittivity实参就一个数相对介电常数。论文后面所有的曲线都是拿这个比值当自变量扫出来的。2.2 两个系数的具体形式把 ε₂/ε₁ 当作唯一变量设入射角为 θix ε₂/ε₁四种系数分别是垂直极化R⊥ (cosθi − √(x − sin²θi)) / (cosθi √(x − sin²θi))T⊥ 2cosθi / (cosθi √(x − sin²θi))平行极化R∥ (x·cosθi − √(x − sin²θi)) / (x·cosθi √(x − sin²θi))T∥ 2√x·cosθi / (x·cosθi √(x − sin²θi))这四个式子看着长结构其实很整齐。分子是两侧“法向阻抗”的差分母是两侧“法向阻抗”的和差别只在平行极化多乘了一个 x。根号里的 x − sin²θi 是从折射定律推出来的透射波在界面法线方向的波矢分量这决定了后面临界角、全反射一系列现象都藏在根号什么时候变负号。提示x 是透射介质电容率除以入射介质电容率千万别写反。x 1 表示波从光疏进入光密x 1 表示波从光密进入光疏。2.3 电容率比值的影响方向先记住 x 变大系数变小论文的摘要里有一句关键结论透射媒质电容率增大时反射系数和透射系数一般会变小反之变大。这个“一般”很值得琢磨。从物理直觉看x 越接近 1两种介质的电性质差异越小界面越“透明”反射自然弱。论文用两组数值验证了这一点从光密到光疏时取 x 0.58 和 0.68从光疏到光密时取 x 1.45 和 1.75每组都是较大的 x 对应较低的反射系数曲线。但“一般会变小”不是“所有角度都变小”。在临界角附近曲线会明显抬升反射和透射系数在临近全反射时一起增大。所以想直接用这句结论去预估某个特定角度的系数大概率会偏差。正确用法是先看 x 落在哪个区间再判断是否靠近临界角最后才用曲线取值。3. 三类边界现象临界角、布儒斯特角与斜滑投射的工程边界3.1 从光密到光疏临界角内反射增强临界角外全反射当 x 1也就是入射介质更密时根号里的 x − sin²θi 会存在一个过零点。令 sin²θc x这个角度就是临界角。入射角小于 θc 时反射系数随角度增大而增大入射角大于 θc 时根号变成负数透射波不再是向第二介质传播的行波R⊥ 和 R∥ 的模值都变成 1入射波全部被反射。论文图 1 模拟的就是这个区间。x 从 0.58 增大到 0.68 时两条曲线整体下移原因就是前面说的比值更接近 1界面差异变小。但注意一个容易误读的地方论文里说“入射角连续增大到临界角时反射系数、透射系数都增大”这里的透射系数并不代表有能量透射过去它在临界角处只是一个数学上的边界值。我们用这个公式做仿真时入射角一到 θc透射系数就是 2反射系数是 1两者之差仍然是 1但物理上已经没有行波透射了。3.2 平行极化的“无反射”与布儒斯特角垂直极化没有这个待遇令 R∥ 0可以得到一个非常干净的解析条件x·cosθi √(x − sin²θi)两边平方整理后得到 sin²θi x / (1 x)等价形式是 tanθi √x。这个角度叫布儒斯特角也叫起偏角。有意思的是平行极化出现无反射的条件里没有要求 x 大于还是小于 1。论文的模拟覆盖了光密到光疏和光疏到光密两个方向结论一致无论哪种情况平行极化波都存在无反射角度而垂直极化波除非 x 1两种介质完全相同否则永远没有反射系数过零的机会。这正是论文结论里“平行极化有无反射垂直极化没有”的由来。这个特性在工程里非常实用。一个无固定极化方向的波以布儒斯特角入射时反射波里只剩垂直极化分量反射过程本身就相当于一个偏振滤波器。后面会专门展开用法。3.3 斜滑投射的近极限行为R -1 的工程后果当入射角 θi 逼近 90°也就是波几乎贴着界面掠过时cosθi 趋于 0。代入任何一组公式都能看出R⊥ 和 R∥ 都趋于 -1T⊥ 和 T∥ 都趋于 0。这意味着所有极化状态的入射波都会被完全反射且反射波与入射波的电场相位差正好是 180°。论文里用这个结论解释了为什么“尽量倾斜观察物体表面时总是比较明亮”也解释了地面雷达的低空盲区低空目标反射的信号经地面反射后在雷达接收天线处与直达波干涉相消本质上就是掠射反射系数接近 -1 造成的。这种极限行为不依赖 x 的具体取值所以它对任何平坦介质界面都成立。做多径链路仿真时如果忽略这个 -1 号算出来的接收电平可能比真实值高出好几个 dB。4. 避坑与常见问题复现这篇论文前必须确认的五个细节4.1 现象画出来的曲线和论文完全相反临界角找不到了原因把 x ε₂/ε₁ 写成了 ε₁/ε₂。这会让光密到光疏的场景变成光疏到光密临界角自然消失布儒斯特角的位置也会整体偏移。这个错非常隐蔽因为公式本身不自检图形趋势看起来依然“合理”只是怎么都对不上论文的图 1 和图 2。解决写代码前先做一个基线校验。取 x 0.58θi 0 时R⊥ (1 − √0.58) / (1 √0.58) ≈ -0.121T⊥ ≈ 0.879。如果算出来符号不对第一件事查 x 的分子分母别急着调绘图参数。4.2 现象反射系数画出来是负数领导说“反射系数怎么能为负”原因反射系数为负只表示反射波电场与参考方向相差 180°是相位信息不是能量倒流。从光疏到光密时垂直极化的 R⊥ 在 0°~90° 范围内始终小于 0这是正常结果。有人为了“好看”给曲线取绝对值取完之后反射和透射曲线就不再是平行的差 1 关系后面所有核对项全部失效。解决直接画原始值。你只需要在图上加一条 0 参考线并在注释里写明“负值对应相位翻转”比取 abs 科学得多。4.3 现象θi 超过临界角后裸奔np.sqrt 算出 nan图形中断原因x 1 时sin²θi 可以比 x 大根号里是负数。数值上这叫全反射区严格来说 R 变成复数不能直接当实数画但 numpy 遇到负数开方只会返回 nan不给你任何报错曲线就出现一块空白。解决把画图角度范围限制在 [0, θc − 0.001]或者用 np.where 把 radicand 的负值替换成 np.nan。后者更安全函数在超界区域照样返回 nan画图时自动跳过不会污染正常区域。4.4 现象垂直入射时 R⊥ 和 R∥ 符号相反怀疑公式抄错了原因这不是错。垂直极化和平行极化对电场参考方向的约定不同垂直入射时两个反射系数本身就该差一个符号。论文给的公式里x 4 时 R⊥ ≈ -0.333R∥ ≈ 0.333这是电磁理论中 s 波和 p 波在界面反射的正常差异。很多第一次推导的人在这里困惑半天翻别的教材才发现是定义问题。解决做验证时只核对同一极化的极限值不要跨极化比较符号。想确认公式对不对用 θi 0 处的解析值、布儒斯特角位置、临界角位置三个锚点对照即可。4.5 现象平行极化也画“两条平行曲线”差 1 对不上原因T⊥ − R⊥ 1 是垂直极化的代数恒等式因为两个公式分母相同分子相减后恰好消掉。平行极化的分子是 x·cosθi 和 2√x·cosθi相减不会得到常数除非 x 1。论文里明确说的是“垂直极化的反射系数和透射系数曲线是两条平行曲线”没有说平行极化也是。解决凡是要用“差 1”来检查曲线形态的地方只针对垂直极化。平行极化用布儒斯特角过零和斜滑投射趋近 -1 这两个特性来核对别硬套。5. 用 Python 复现论文图 1 和图 2公式落地到曲线的完整步骤5.1 先搭好系数函数注意根号内的符号先准备一个矢量化的系数函数。输入入射角度和 x ε₂/ε₁返回四个系数。根号里的负值交给 numpy 的 where 处理避免直接 nan 报错。import numpy as np import matplotlib.pyplot as plt def fresnel_coefficients(theta_deg, x): # x eps2 / eps1入射介质为 eps1透射介质为 eps2 theta np.radians(theta_deg) c np.cos(theta) s2 np.sin(theta) ** 2 # 从光密到光疏时theta 超过临界角根号内为负 # 此时进入全反射区系数应视为复数这里用 nan 占位 radicand x - s2 m np.sqrt(np.where(radicand 0.0, radicand, np.nan)) # 垂直极化s 波 Rs (c - m) / (c m) Ts 2.0 * c / (c m) # 平行极化p 波 Rp (x * c - m) / (x * c m) Tp 2.0 * np.sqrt(x) * c / (x * c m) return Rs, Ts, Rp, Tp这个函数把论文公式 (1)(4) 直接翻译成了 NumPy 运算。theta_deg 是入射角单位用度是为了后面画图时横轴可读x 是电容率比值。返回顺序固定为 Rs, Ts, Rp, Tp对应垂直极化的反射、透射和平行极化的反射、透射。5.2 复现光密到光疏的两组曲线论文图 1 的模拟参数是 x 0.58 和 0.68。临界角由 sinθc √x 决定画图时让角度停在临界角前一点点避免根号进入负数区间。# 光密到光疏入射介质电容率大于透射介质 x1 0.58 theta_c1 np.degrees(np.arcsin(np.sqrt(x1))) # 临界角约 49.8 度 theta1 np.linspace(0, theta_c1 - 0.01, 400) Rs1, Ts1, Rp1, Tp1 fresnel_coefficients(theta1, x1) # 光密到光疏第二组x 更接近 1差异变小 x2 0.68 theta_c2 np.degrees(np.arcsin(np.sqrt(x2))) theta2 np.linspace(0, theta_c2 - 0.01, 400) Rs2, Ts2, Rp2, Tp2 fresnel_coefficients(theta2, x2) fig, axes plt.subplots(1, 2, figsize(11, 4.2)) # 左图垂直极化 axes[0].plot(theta1, Ts1, -, labelT_s, x0.58, lw2) axes[0].plot(theta1, Rs1, --, labelR_s, x0.58, lw2) axes[0].plot(theta2, Ts2, -, labelT_s, x0.68, lw2) axes[0].plot(theta2, Rs2, --, labelR_s, x0.68, lw2) axes[0].set_title(vertical polarization (s wave)) axes[0].set_xlabel(incident angle / deg) axes[0].set_ylabel(coefficient) axes[0].legend(fontsize9) # 右图平行极化 axes[1].plot(theta1, Tp1, -, labelT_p, x0.58, lw2) axes[1].plot(theta1, Rp1, --, labelR_p, x0.58, lw2) axes[1].plot(theta2, Tp2, -, labelT_p, x0.68, lw2) axes[1].plot(theta2, Rp2, --, labelR_p, x0.68, lw2) axes[1].axhline(0, colork, lw0.8) axes[1].set_title(parallel polarization (p wave)) axes[1].set_xlabel(incident angle / deg) axes[1].legend(fontsize9) plt.tight_layout() plt.show()这段代码里横轴没有连接到临界角因为论文只在临界角内给出了实数值。左图右图都叠加了两条 x 的曲线可以直接对比“比值增大时系数变小”的结论。注意右图里 R_p 曲线穿过 0 点那就是布儒斯特角的位置左图 R_s 始终不为 0这是两图的本质差异。5.3 换到光疏到光密再跑一遍论文图 2 用的 x 1.45 和 1.75。这个区间没有临界角入射角可以一直扫到接近 90°正好观察斜滑投射的极限行为。# 光疏到光密入射介质电容率小于透射介质 x3 1.45 x4 1.75 theta3 np.linspace(0, 89.5, 500) # 避开 90 度避免除零 Rs3, Ts3, Rp3, Tp3 fresnel_coefficients(theta3, x3) Rs4, Ts4, Rp4, Tp4 fresnel_coefficients(theta3, x4) fig, axes plt.subplots(1, 2, figsize(11, 4.2)) axes[0].plot(theta3, Ts3, -, labelT_s, x1.45, lw2) axes[0].plot(theta3, Rs3, --, labelR_s, x1.45, lw2) axes[0].plot(theta3, Ts4, -, labelT_s, x1.75, lw2) axes[0].plot(theta3, Rs4, --, labelR_s, x1.75, lw2) axes[0].set_title(vertical polarization (s wave)) axes[0].set_xlabel(incident angle / deg) axes[0].set_ylabel(coefficient) axes[0].legend(fontsize9) axes[1].plot(theta3, Tp3, -, labelT_p, x1.45, lw2) axes[1].plot(theta3, Rp3, --, labelR_p, x1.45, lw2) axes[1].plot(theta3, Tp4, -, labelT_p, x1.75, lw2) axes[1].plot(theta3, Rp4, --, labelR_p, x1.75, lw2) axes[1].axhline(0, colork, lw0.8) axes[1].set_title(parallel polarization (p wave)) axes[1].set_xlabel(incident angle / deg) axes[1].legend(fontsize9) plt.tight_layout() plt.show()跑完这组你会发现一个典型特征垂直极化的反射系数整个区间都是负的透射系数在 0 到 2 之间两条曲线严格平行同一角度差值恒为 1。平行极化这边R_p 先正后负穿过零点的角度就是布儒斯特角角度随 x 增大而增大。到 θi 90° 的极限位置四种系数全部收拢到 R -1、T 0。5.4 对着论文图核对四条判断画完图之后别急着收工拿下面这张核对表逐项过一遍能挡住八成“公式抄错但图看着挺顺”的情况。核对项预期结果对应代码位置垂直极化 θi 0 处反射系数(1 − √x) / (1 √x)负值Rs[0]垂直极化全区间T_s − R_s 恒等于 1Ts − Rs平行极化布儒斯特角tanθB √xR_p 过零Rp 与 0 交点斜滑投射 θi → 90°R_s R_p → -1T_s T_p → 0曲线尾端第一次复现时把这四项跑通基本可以确定系数函数没有大问题。看到垂直极化的平行差 1 曲线其实就是论文里的原话两条平行曲线。用它来检查公式比自己对着屏幕怀疑人生快得多。6. 把界面极化用在实操里偏振获取、极化匹配与低空盲区判断论文最后一段提到以布儒斯特角入射时反射波里只剩垂直极化分量。这个结论可以直接搬到光学实验和微波测试里。想得到一个相对纯净的线极化信号手边又没有偏振片时让波先经过一个介质界面调整入射角到 arctan(√x)反射波就是近乎纯净的垂直极化波。角度算不准时看反射波功率最小对应的入射角最省事这个习惯我在实验室里用过很多回。极化匹配也值得多说一句。圆极化天线能够接收任意取向的线极化波但代价是固有 3 dB 损耗。论文里强调了接收天线与入射波极化状态匹配的问题背后实质正是界面反射改变了波的极化构成。比如一个线极化波以布儒斯特角打到水面或玻璃面反射波极化方向会被“过滤”一次你再拿原来方向的线极化天线去收效率就对不上了。算链路预算时把界面看作一个极化变换器比事后打补偿系数更接近物理本质。低空盲区那个例子我建议做雷达仿真的朋友把斜滑投射写成一条固定规则地面反射路径的反射系数在低掠射角时直接取 -1不要用自由空间传播模型。这个近似可以解释很多实测里“目标明明在视距内接收机却收不到信号”的现象。换个场景如果把天线架高让反射点后移等效反射角变大反射系数不再接近 -1盲区深度就会明显下降。论文给的是现象工程上调整天线高度和极化方式就是解法。从那以后我拿到任何一份界面反射相关的论文或代码第一件事就是先在纸面上把 θi 0 和 θi π/2 这两个极限位置算一遍锚定曲线形态再谈数值模拟。这个小习惯帮我避掉了后面好几次把公式写反的低级失误希望帮到你。本文还有配套的精品资源点击获取 从边界条件到布儒斯特角：介质界面极化电磁波与Python仿真
简介《介质界面极化电磁波的传播特性研究》是一篇学术论文围绕电磁波在介质界面的极化行为系统分析了垂直极化波与平行极化波的反射、透射系数及其随介质电容率的变化规律。资料源自期刊论文适合通信、雷达及电磁场专业的高年级学生与工程技术人员参考可用于理解全反射、无反射等现象的理论推导与数值模拟。压缩包内为1个PDF文件容量249KB内容结构完整含边界条件推导、反射透射系数公式及模拟曲线便于按需查阅或打印学习。目前已有94人学习该资源。通过阅读可掌握介质电容率对反射与透射系数的影响规律了解斜滑投射及光密/光疏介质条件下的传播特性差异为极化控制信息传播技术或相关反隐身、偏振分析应用提供理论支持。1. 介质界面极化电磁波为什么先读这一篇再碰极化匹配聊极化匹配之前先说一个反直觉的结论地面雷达之所以存在低空盲区很大一部分原因就藏在介质界面极化电磁波的反射规律里。斜滑投射时不论垂直极化还是平行极化反射系数都趋近 -1反射波与入射波相位相反直达波和地面反射波在天线处直接相消。这是 2014 年《内蒙古师范大学学报自然科学汉文版》一篇论文末尾给出的解释而整篇论文做的就是把介质界面上这两种极化波的反射系数、透射系数用电磁场边界条件推出来再给出一组可复现的数值模拟结果。对做电磁场课程设计、天线极化匹配、雷达链路仿真的人来说这份 PDF 最有价值的地方在于公式给得直接参数只有电容率比值 ε₂/ε₁场景又覆盖了全反射、无反射、斜滑投射这些最容易概念混淆的边界情况。你可以把它当作一份能落地到代码里的理论依据而不是一篇看过就忘的论文。下面我按“公式怎么来、结论怎么用、代码怎么写、坑在哪”的顺序把这四页纸拆开讲透。2. 从边界条件到反射透射系数四个公式决定你的曲线走向2.1 垂直极化与平行极化的定义先分清看的是哪个电场分量任意极化的平面电磁波都可以分解成两个正交分量电场方向垂直于入射面的叫垂直极化波电场方向在入射面内的叫平行极化波。这里的“入射面”由入射波矢和界面法线共同确定不是一个固定的水平面或垂直面。很多做天线仿真的同事第一次看这篇论文时都在这翻车——把垂直极化理解成“电场垂直于地面”但论文里的垂直极化是相对入射面定义的两者只在特殊场景下重合。在非磁性介质假设下μ₁≈μ₂≈μ₀界面两侧的磁导率差异被消掉最后决定反射、透射行为的只剩电容率的比值 ε₂/ε₁。电容率 ε 在中文教材里常叫介电常数英文习惯叫 permittivity实参就一个数相对介电常数。论文后面所有的曲线都是拿这个比值当自变量扫出来的。2.2 两个系数的具体形式把 ε₂/ε₁ 当作唯一变量设入射角为 θix ε₂/ε₁四种系数分别是垂直极化R⊥ (cosθi − √(x − sin²θi)) / (cosθi √(x − sin²θi))T⊥ 2cosθi / (cosθi √(x − sin²θi))平行极化R∥ (x·cosθi − √(x − sin²θi)) / (x·cosθi √(x − sin²θi))T∥ 2√x·cosθi / (x·cosθi √(x − sin²θi))这四个式子看着长结构其实很整齐。分子是两侧“法向阻抗”的差分母是两侧“法向阻抗”的和差别只在平行极化多乘了一个 x。根号里的 x − sin²θi 是从折射定律推出来的透射波在界面法线方向的波矢分量这决定了后面临界角、全反射一系列现象都藏在根号什么时候变负号。提示x 是透射介质电容率除以入射介质电容率千万别写反。x 1 表示波从光疏进入光密x 1 表示波从光密进入光疏。2.3 电容率比值的影响方向先记住 x 变大系数变小论文的摘要里有一句关键结论透射媒质电容率增大时反射系数和透射系数一般会变小反之变大。这个“一般”很值得琢磨。从物理直觉看x 越接近 1两种介质的电性质差异越小界面越“透明”反射自然弱。论文用两组数值验证了这一点从光密到光疏时取 x 0.58 和 0.68从光疏到光密时取 x 1.45 和 1.75每组都是较大的 x 对应较低的反射系数曲线。但“一般会变小”不是“所有角度都变小”。在临界角附近曲线会明显抬升反射和透射系数在临近全反射时一起增大。所以想直接用这句结论去预估某个特定角度的系数大概率会偏差。正确用法是先看 x 落在哪个区间再判断是否靠近临界角最后才用曲线取值。3. 三类边界现象临界角、布儒斯特角与斜滑投射的工程边界3.1 从光密到光疏临界角内反射增强临界角外全反射当 x 1也就是入射介质更密时根号里的 x − sin²θi 会存在一个过零点。令 sin²θc x这个角度就是临界角。入射角小于 θc 时反射系数随角度增大而增大入射角大于 θc 时根号变成负数透射波不再是向第二介质传播的行波R⊥ 和 R∥ 的模值都变成 1入射波全部被反射。论文图 1 模拟的就是这个区间。x 从 0.58 增大到 0.68 时两条曲线整体下移原因就是前面说的比值更接近 1界面差异变小。但注意一个容易误读的地方论文里说“入射角连续增大到临界角时反射系数、透射系数都增大”这里的透射系数并不代表有能量透射过去它在临界角处只是一个数学上的边界值。我们用这个公式做仿真时入射角一到 θc透射系数就是 2反射系数是 1两者之差仍然是 1但物理上已经没有行波透射了。3.2 平行极化的“无反射”与布儒斯特角垂直极化没有这个待遇令 R∥ 0可以得到一个非常干净的解析条件x·cosθi √(x − sin²θi)两边平方整理后得到 sin²θi x / (1 x)等价形式是 tanθi √x。这个角度叫布儒斯特角也叫起偏角。有意思的是平行极化出现无反射的条件里没有要求 x 大于还是小于 1。论文的模拟覆盖了光密到光疏和光疏到光密两个方向结论一致无论哪种情况平行极化波都存在无反射角度而垂直极化波除非 x 1两种介质完全相同否则永远没有反射系数过零的机会。这正是论文结论里“平行极化有无反射垂直极化没有”的由来。这个特性在工程里非常实用。一个无固定极化方向的波以布儒斯特角入射时反射波里只剩垂直极化分量反射过程本身就相当于一个偏振滤波器。后面会专门展开用法。3.3 斜滑投射的近极限行为R -1 的工程后果当入射角 θi 逼近 90°也就是波几乎贴着界面掠过时cosθi 趋于 0。代入任何一组公式都能看出R⊥ 和 R∥ 都趋于 -1T⊥ 和 T∥ 都趋于 0。这意味着所有极化状态的入射波都会被完全反射且反射波与入射波的电场相位差正好是 180°。论文里用这个结论解释了为什么“尽量倾斜观察物体表面时总是比较明亮”也解释了地面雷达的低空盲区低空目标反射的信号经地面反射后在雷达接收天线处与直达波干涉相消本质上就是掠射反射系数接近 -1 造成的。这种极限行为不依赖 x 的具体取值所以它对任何平坦介质界面都成立。做多径链路仿真时如果忽略这个 -1 号算出来的接收电平可能比真实值高出好几个 dB。4. 避坑与常见问题复现这篇论文前必须确认的五个细节4.1 现象画出来的曲线和论文完全相反临界角找不到了原因把 x ε₂/ε₁ 写成了 ε₁/ε₂。这会让光密到光疏的场景变成光疏到光密临界角自然消失布儒斯特角的位置也会整体偏移。这个错非常隐蔽因为公式本身不自检图形趋势看起来依然“合理”只是怎么都对不上论文的图 1 和图 2。解决写代码前先做一个基线校验。取 x 0.58θi 0 时R⊥ (1 − √0.58) / (1 √0.58) ≈ -0.121T⊥ ≈ 0.879。如果算出来符号不对第一件事查 x 的分子分母别急着调绘图参数。4.2 现象反射系数画出来是负数领导说“反射系数怎么能为负”原因反射系数为负只表示反射波电场与参考方向相差 180°是相位信息不是能量倒流。从光疏到光密时垂直极化的 R⊥ 在 0°~90° 范围内始终小于 0这是正常结果。有人为了“好看”给曲线取绝对值取完之后反射和透射曲线就不再是平行的差 1 关系后面所有核对项全部失效。解决直接画原始值。你只需要在图上加一条 0 参考线并在注释里写明“负值对应相位翻转”比取 abs 科学得多。4.3 现象θi 超过临界角后裸奔np.sqrt 算出 nan图形中断原因x 1 时sin²θi 可以比 x 大根号里是负数。数值上这叫全反射区严格来说 R 变成复数不能直接当实数画但 numpy 遇到负数开方只会返回 nan不给你任何报错曲线就出现一块空白。解决把画图角度范围限制在 [0, θc − 0.001]或者用 np.where 把 radicand 的负值替换成 np.nan。后者更安全函数在超界区域照样返回 nan画图时自动跳过不会污染正常区域。4.4 现象垂直入射时 R⊥ 和 R∥ 符号相反怀疑公式抄错了原因这不是错。垂直极化和平行极化对电场参考方向的约定不同垂直入射时两个反射系数本身就该差一个符号。论文给的公式里x 4 时 R⊥ ≈ -0.333R∥ ≈ 0.333这是电磁理论中 s 波和 p 波在界面反射的正常差异。很多第一次推导的人在这里困惑半天翻别的教材才发现是定义问题。解决做验证时只核对同一极化的极限值不要跨极化比较符号。想确认公式对不对用 θi 0 处的解析值、布儒斯特角位置、临界角位置三个锚点对照即可。4.5 现象平行极化也画“两条平行曲线”差 1 对不上原因T⊥ − R⊥ 1 是垂直极化的代数恒等式因为两个公式分母相同分子相减后恰好消掉。平行极化的分子是 x·cosθi 和 2√x·cosθi相减不会得到常数除非 x 1。论文里明确说的是“垂直极化的反射系数和透射系数曲线是两条平行曲线”没有说平行极化也是。解决凡是要用“差 1”来检查曲线形态的地方只针对垂直极化。平行极化用布儒斯特角过零和斜滑投射趋近 -1 这两个特性来核对别硬套。5. 用 Python 复现论文图 1 和图 2公式落地到曲线的完整步骤5.1 先搭好系数函数注意根号内的符号先准备一个矢量化的系数函数。输入入射角度和 x ε₂/ε₁返回四个系数。根号里的负值交给 numpy 的 where 处理避免直接 nan 报错。import numpy as np import matplotlib.pyplot as plt def fresnel_coefficients(theta_deg, x): # x eps2 / eps1入射介质为 eps1透射介质为 eps2 theta np.radians(theta_deg) c np.cos(theta) s2 np.sin(theta) ** 2 # 从光密到光疏时theta 超过临界角根号内为负 # 此时进入全反射区系数应视为复数这里用 nan 占位 radicand x - s2 m np.sqrt(np.where(radicand 0.0, radicand, np.nan)) # 垂直极化s 波 Rs (c - m) / (c m) Ts 2.0 * c / (c m) # 平行极化p 波 Rp (x * c - m) / (x * c m) Tp 2.0 * np.sqrt(x) * c / (x * c m) return Rs, Ts, Rp, Tp这个函数把论文公式 (1)(4) 直接翻译成了 NumPy 运算。theta_deg 是入射角单位用度是为了后面画图时横轴可读x 是电容率比值。返回顺序固定为 Rs, Ts, Rp, Tp对应垂直极化的反射、透射和平行极化的反射、透射。5.2 复现光密到光疏的两组曲线论文图 1 的模拟参数是 x 0.58 和 0.68。临界角由 sinθc √x 决定画图时让角度停在临界角前一点点避免根号进入负数区间。# 光密到光疏入射介质电容率大于透射介质 x1 0.58 theta_c1 np.degrees(np.arcsin(np.sqrt(x1))) # 临界角约 49.8 度 theta1 np.linspace(0, theta_c1 - 0.01, 400) Rs1, Ts1, Rp1, Tp1 fresnel_coefficients(theta1, x1) # 光密到光疏第二组x 更接近 1差异变小 x2 0.68 theta_c2 np.degrees(np.arcsin(np.sqrt(x2))) theta2 np.linspace(0, theta_c2 - 0.01, 400) Rs2, Ts2, Rp2, Tp2 fresnel_coefficients(theta2, x2) fig, axes plt.subplots(1, 2, figsize(11, 4.2)) # 左图垂直极化 axes[0].plot(theta1, Ts1, -, labelT_s, x0.58, lw2) axes[0].plot(theta1, Rs1, --, labelR_s, x0.58, lw2) axes[0].plot(theta2, Ts2, -, labelT_s, x0.68, lw2) axes[0].plot(theta2, Rs2, --, labelR_s, x0.68, lw2) axes[0].set_title(vertical polarization (s wave)) axes[0].set_xlabel(incident angle / deg) axes[0].set_ylabel(coefficient) axes[0].legend(fontsize9) # 右图平行极化 axes[1].plot(theta1, Tp1, -, labelT_p, x0.58, lw2) axes[1].plot(theta1, Rp1, --, labelR_p, x0.58, lw2) axes[1].plot(theta2, Tp2, -, labelT_p, x0.68, lw2) axes[1].plot(theta2, Rp2, --, labelR_p, x0.68, lw2) axes[1].axhline(0, colork, lw0.8) axes[1].set_title(parallel polarization (p wave)) axes[1].set_xlabel(incident angle / deg) axes[1].legend(fontsize9) plt.tight_layout() plt.show()这段代码里横轴没有连接到临界角因为论文只在临界角内给出了实数值。左图右图都叠加了两条 x 的曲线可以直接对比“比值增大时系数变小”的结论。注意右图里 R_p 曲线穿过 0 点那就是布儒斯特角的位置左图 R_s 始终不为 0这是两图的本质差异。5.3 换到光疏到光密再跑一遍论文图 2 用的 x 1.45 和 1.75。这个区间没有临界角入射角可以一直扫到接近 90°正好观察斜滑投射的极限行为。# 光疏到光密入射介质电容率小于透射介质 x3 1.45 x4 1.75 theta3 np.linspace(0, 89.5, 500) # 避开 90 度避免除零 Rs3, Ts3, Rp3, Tp3 fresnel_coefficients(theta3, x3) Rs4, Ts4, Rp4, Tp4 fresnel_coefficients(theta3, x4) fig, axes plt.subplots(1, 2, figsize(11, 4.2)) axes[0].plot(theta3, Ts3, -, labelT_s, x1.45, lw2) axes[0].plot(theta3, Rs3, --, labelR_s, x1.45, lw2) axes[0].plot(theta3, Ts4, -, labelT_s, x1.75, lw2) axes[0].plot(theta3, Rs4, --, labelR_s, x1.75, lw2) axes[0].set_title(vertical polarization (s wave)) axes[0].set_xlabel(incident angle / deg) axes[0].set_ylabel(coefficient) axes[0].legend(fontsize9) axes[1].plot(theta3, Tp3, -, labelT_p, x1.45, lw2) axes[1].plot(theta3, Rp3, --, labelR_p, x1.45, lw2) axes[1].plot(theta3, Tp4, -, labelT_p, x1.75, lw2) axes[1].plot(theta3, Rp4, --, labelR_p, x1.75, lw2) axes[1].axhline(0, colork, lw0.8) axes[1].set_title(parallel polarization (p wave)) axes[1].set_xlabel(incident angle / deg) axes[1].legend(fontsize9) plt.tight_layout() plt.show()跑完这组你会发现一个典型特征垂直极化的反射系数整个区间都是负的透射系数在 0 到 2 之间两条曲线严格平行同一角度差值恒为 1。平行极化这边R_p 先正后负穿过零点的角度就是布儒斯特角角度随 x 增大而增大。到 θi 90° 的极限位置四种系数全部收拢到 R -1、T 0。5.4 对着论文图核对四条判断画完图之后别急着收工拿下面这张核对表逐项过一遍能挡住八成“公式抄错但图看着挺顺”的情况。核对项预期结果对应代码位置垂直极化 θi 0 处反射系数(1 − √x) / (1 √x)负值Rs[0]垂直极化全区间T_s − R_s 恒等于 1Ts − Rs平行极化布儒斯特角tanθB √xR_p 过零Rp 与 0 交点斜滑投射 θi → 90°R_s R_p → -1T_s T_p → 0曲线尾端第一次复现时把这四项跑通基本可以确定系数函数没有大问题。看到垂直极化的平行差 1 曲线其实就是论文里的原话两条平行曲线。用它来检查公式比自己对着屏幕怀疑人生快得多。6. 把界面极化用在实操里偏振获取、极化匹配与低空盲区判断论文最后一段提到以布儒斯特角入射时反射波里只剩垂直极化分量。这个结论可以直接搬到光学实验和微波测试里。想得到一个相对纯净的线极化信号手边又没有偏振片时让波先经过一个介质界面调整入射角到 arctan(√x)反射波就是近乎纯净的垂直极化波。角度算不准时看反射波功率最小对应的入射角最省事这个习惯我在实验室里用过很多回。极化匹配也值得多说一句。圆极化天线能够接收任意取向的线极化波但代价是固有 3 dB 损耗。论文里强调了接收天线与入射波极化状态匹配的问题背后实质正是界面反射改变了波的极化构成。比如一个线极化波以布儒斯特角打到水面或玻璃面反射波极化方向会被“过滤”一次你再拿原来方向的线极化天线去收效率就对不上了。算链路预算时把界面看作一个极化变换器比事后打补偿系数更接近物理本质。低空盲区那个例子我建议做雷达仿真的朋友把斜滑投射写成一条固定规则地面反射路径的反射系数在低掠射角时直接取 -1不要用自由空间传播模型。这个近似可以解释很多实测里“目标明明在视距内接收机却收不到信号”的现象。换个场景如果把天线架高让反射点后移等效反射角变大反射系数不再接近 -1盲区深度就会明显下降。论文给的是现象工程上调整天线高度和极化方式就是解法。从那以后我拿到任何一份界面反射相关的论文或代码第一件事就是先在纸面上把 θi 0 和 θi π/2 这两个极限位置算一遍锚定曲线形态再谈数值模拟。这个小习惯帮我避掉了后面好几次把公式写反的低级失误希望帮到你。本文还有配套的精品资源点击获取