☰
介质界面极化电磁波传播特性:反射透射系数推导与工程应用
2026/10/7 22:17:57 网站建设 项目流程

简介:这份PDF文献《介质界面极化电磁波的传播特性研究》面向电磁场与微波技术、通信工程及雷达技术方向的学习者与研究人员,聚焦电磁波在介质分界面处的极化传播规律这一基础而关键的问题。资源包内仅含1个PDF文件,约249KB,为期刊论文全文,便于快速查阅与引用。文章从电磁场边界条件出发,系统推导垂直极化波与平行极化波的反射系数、透射系数表达式,并通过数值计算与模拟,讨论介质电容率对反射、折射系数的影响,分析斜滑投射下的全反射现象、平行极化波的“无反射”特性以及垂直极化反射与透射曲线的平行关系。内容还延伸至圆极化天线接收、极化域反隐身、大气气溶胶散射消偏等应用场景,兼具理论推导与仿真验证。目前已有94人学习,适合作为通信技术、雷达技术领域的技术开发参考与专业指导资料,也可为极化控制信息传播的后续研究提供理论依据。

1. 介质界面极化电磁波:从边界条件到工程直觉

搞天线、做雷达或者调射频链路的人,迟早会撞上一个绕不开的问题:电磁波打到两种介质的交界面上,到底有多少能量反射回来,多少透过去,相位翻不翻?很多人第一次接触这个,是在仿真软件里拖一个平面波源,设个入射角,看 S 参数曲线——但曲线为什么长那样,心里没底。渭南师范学院雷前召等人写的这篇《介质界面极化电磁波的传播特性研究》,做的就是把这层窗户纸捅破:从电磁场边界条件出发,推出垂直极化波和平行极化波在介质界面的反射系数与透射系数解析式,再用数值计算把光密到光疏、光疏到光密两种场景下的曲线画出来。它不教你怎么调天线,但它告诉你天线为什么这么调。适合谁看?做雷达极化处理、设计圆极化天线、搞遥感偏振分析的从业者,以及电磁场与电磁波课程学到界面反射这一章想找数值直觉的学生。全文的核心结论就几条,但每一条背后都有明确的物理图像和可复现的计算路径。

2. 两种极化波的反射透射系数:公式怎么来的、参数怎么代

2.1 垂直极化与平行极化的边界条件推导

任意极化的平面电磁波,总可以拆成两个正交分量的叠加:电场方向垂直于入射面的叫垂直极化波,电场方向平行于入射面的叫平行极化波。这个拆分不是数学游戏,而是因为这两种极化在界面上的行为完全不同——一个会出现无反射点,另一个不会。

推导的起点是电磁场边界条件。对于常见的非磁性介质,磁导率近似等于真空磁导率,界面两侧的电场切向分量连续、磁场切向分量连续。把入射波、反射波、透射波的场表达式代入,消去振幅,就能得到反射系数和透射系数。

垂直极化波的反射系数和透射系数分别为:

R⊥ = (cosθᵢ - √(ε₂/ε₁ - sin²θᵢ)) / (cosθᵢ + √(ε₂/ε₁ - sin²θᵢ))

T⊥ = 2cosθᵢ / (cosθᵢ + √(ε₂/ε₁ - sin²θᵢ))

平行极化波的对应表达式为:

R∥ = ((ε₂/ε₁)cosθᵢ - √(ε₂/ε₁ - sin²θᵢ)) / ((ε₂/ε₁)cosθᵢ + √(ε₂/ε₁ - sin²θᵢ))

T∥ = 2√(ε₂/ε₁)cosθᵢ / ((ε₂/ε₁)cosθᵢ + √(ε₂/ε₁ - sin²θᵢ))

这里 θᵢ 是入射角,ε₁ 是入射侧介质的电容率,ε₂ 是透射侧介质的电容率。注意,这四个公式里没有出现频率——在非色散介质假设下,反射和透射系数只与入射角和电容率比值有关,与频率无关。这一点在工程上很重要:如果你用不同频段测同一个界面,反射系数曲线应该重合,不重合就说明介质有色散或者测量有问题。

2.2 数值计算:用 Python 把四条曲线画出来

公式有了,但光看表达式很难建立直觉。我一般会直接用 Python 把曲线画出来,改参数看变化。下面这段代码复现了论文中光密到光疏(ε₂/ε₁ < 1)和光疏到光密(ε₂/ε₁ > 1)两种场景下的反射、透射系数曲线。

import numpy as np import matplotlib.pyplot as plt def reflection_transmission(theta_i, eps_ratio): """ 计算垂直极化和平行极化波的反射系数与透射系数 theta_i: 入射角,弧度 eps_ratio: eps2 / eps1 """ cos_i = np.cos(theta_i) sin_i = np.sin(theta_i) sqrt_term = np.sqrt(eps_ratio - sin_i**2 + 0j) # 加0j防止负数开方报错 # 垂直极化 R_perp = (cos_i - sqrt_term) / (cos_i + sqrt_term) T_perp = 2 * cos_i / (cos_i + sqrt_term) # 平行极化 R_para = (eps_ratio * cos_i - sqrt_term) / (eps_ratio * cos_i + sqrt_term) T_para = 2 * np.sqrt(eps_ratio) * cos_i / (eps_ratio * cos_i + sqrt_term) return R_perp, T_perp, R_para, T_para # 入射角从0到90度 theta = np.linspace(0, np.pi/2, 500) # 场景一:光密到光疏,eps2/eps1 = 0.58 R_perp1, T_perp1, R_para1, T_para1 = reflection_transmission(theta, 0.58) # 场景二:光疏到光密,eps2/eps1 = 1.45 R_perp2, T_perp2, R_para2, T_para2 = reflection_transmission(theta, 1.45) fig, axes = plt.subplots(2, 2, figsize=(12, 10)) # 光密到光疏 - 垂直极化 axes[0,0].plot(theta, np.abs(R_perp1), '--', label='|R⊥|') axes[0,0].plot(theta, np.abs(T_perp1), '-', label='|T⊥|') axes[0,0].set_title('光密到光疏 (ε₂/ε₁=0.58) 垂直极化') axes[0,0].set_xlabel('入射角 (rad)') axes[0,0].legend() axes[0,0].grid(True) # 光密到光疏 - 平行极化 axes[0,1].plot(theta, np.abs(R_para1), '--', label='|R∥|') axes[0,1].plot(theta, np.abs(T_para1), '-', label='|T∥|') axes[0,1].set_title('光密到光疏 (ε₂/ε₁=0.58) 平行极化') axes[0,1].set_xlabel('入射角 (rad)') axes[0,1].legend() axes[0,1].grid(True) # 光疏到光密 - 垂直极化 axes[1,0].plot(theta, np.abs(R_perp2), '--', label='|R⊥|') axes[1,0].plot(theta, np.abs(T_perp2), '-', label='|T⊥|') axes[1,0].set_title('光疏到光密 (ε₂/ε₁=1.45) 垂直极化') axes[1,0].set_xlabel('入射角 (rad)') axes[1,0].legend() axes[1,0].grid(True) # 光疏到光密 - 平行极化 axes[1,1].plot(theta, np.abs(R_para2), '--', label='|R∥|') axes[1,1].plot(theta, np.abs(T_para2), '-', label='|T∥|') axes[1,1].set_title('光疏到光密 (ε₂/ε₁=1.45) 平行极化') axes[1,1].set_xlabel('入射角 (rad)') axes[1,1].legend() axes[1,1].grid(True) plt.tight_layout() plt.show()

这段代码的关键参数就两个:eps_ratio控制介质对比度,theta控制扫描范围。sqrt_term加了0j是为了处理全反射时根号内出现负数的情况——这时候反射系数模值会变成 1,透射系数归零。如果你把eps_ratio改成 0.68 或 1.75,就能复现论文图 1 和图 2 中不同电容率比值下的曲线变化。

跑完这段代码,你会看到几个很直观的现象:光密到光疏时,入射角超过临界角后透射系数直接掉到零,反射系数升到 1;平行极化波的反射系数曲线会穿过零点,那个角度就是布儒斯特角;垂直极化波的反射系数始终不穿零。这些不是数学巧合,是边界条件决定的物理必然。

2.3 电容率比值怎么影响曲线形态

论文里明确给出了两组对比:光密到光疏时 ε₂/ε₁ 取 0.58 和 0.68,光疏到光密时取 1.45 和 1.75。对比这两组曲线,能看出一个规律:透射介质的电容率增大时,反射系数和透射系数一般都会变小。

这个结论听起来有点反直觉——透射侧电容率增大,意味着介质更“光密”,按理说反射应该更强才对。但注意,这里说的是电容率比值从 0.58 变到 0.68(仍然小于 1),或者从 1.45 变到 1.75(仍然大于 1),变化方向是比值向 1 靠拢。比值越接近 1,两种介质越相似,界面越“透明”,反射自然越弱。如果比值远离 1,比如从 0.58 变到 0.2,反射就会急剧增强。

工程上选介质的时候,如果你希望减少界面反射,要么让两侧电容率尽量接近,要么利用平行极化在布儒斯特角附近的无反射特性。前者是材料选择问题,后者是入射角控制问题。

3. 全反射、无反射与斜滑投射:三个必须分清的边界场景

3.1 临界角与全反射:光密到光疏的硬边界

当波从光密介质进入光疏介质,即 ε₁ > ε₂ 时,存在一个临界角 θc。入射角超过这个角度,透射波消失,入射波全被反射。临界角满足 sinθc = √(ε₂/ε₁)。

这个现象在光纤通信里是基础——纤芯折射率略高于包层,光在纤芯内多次全反射向前传播。但要注意,全反射并不意味着透射侧完全没有场。实际上在透射侧存在倏逝波,场强随深度指数衰减。论文里没有展开倏逝波的分析,但如果你做近场探测或者光纤耦合,这个衰减长度直接决定耦合效率。

全反射发生时,反射波与入射波相位相反,反射系数趋于 -1。论文里特别提到斜滑投射(入射角趋近 π/2)时,无论什么极化、什么介质组合,反射系数都趋于 -1,透射系数趋于 0。这解释了为什么地面雷达看低空目标时会有盲区——目标贴近地面,电磁波以极小的掠射角打到地面,几乎全被反射,目标回波被地面反射淹没。

3.2 布儒斯特角与平行极化无反射

平行极化波有一个垂直极化波没有的特性:存在一个入射角使反射系数为零。这个角度叫布儒斯特角,满足:

tanθ_B = √(ε₂/ε₁)

当入射角等于布儒斯特角时,平行极化波的反射波消失,全部能量透射过去。垂直极化波则没有这个现象——它的反射系数只在 ε₁ = ε₂ 时才为零,但那已经是没有界面了。

这个特性的工程价值很直接:如果你要拍玻璃后面的物体,用偏振镜把平行极化分量滤掉,玻璃表面的反射光就没了。论文里也提到,实际中常利用这个特性获取偏振光——自然光以布儒斯特角打到界面上,反射光中只剩下垂直极化分量,就是一个天然的偏振器。

但这里有个容易翻车的地方:布儒斯特角只对平行极化成立。如果你用圆极化天线接收,圆极化可以分解为垂直和平行两个分量,平行分量在布儒斯特角不反射,垂直分量照常反射,接收到的信号极化状态就变了。做极化匹配的时候,这个角度必须算进去。

3.3 垂直极化反射系数与透射系数的平行关系

论文里有一个容易被忽略但很实用的结论:垂直极化的反射系数和透射系数曲线是两条平行曲线,同一入射角处的透射系数与反射系数之差为 1。

这个“差为 1”不是模值之差,而是代数之差。因为垂直极化波的反射系数和透射系数满足 T⊥ - R⊥ = 1。你可以从公式直接验证:T⊥ - R⊥ = (2cosθᵢ - cosθᵢ + √(...)) / (cosθᵢ + √(...)) = 1。

这个关系在实验测量里很有用。如果你只能测反射系数,透射系数可以直接算出来,不用另做一套测量。反过来也一样。但注意,这个关系只对垂直极化成立,平行极化不满足这么简洁的线性关系。

4. 避坑与排查:极化电磁波计算中的五个常见翻车点

4.1 把电容率比值和折射率比值搞混

现象:用公式算出来的反射系数和仿真对不上,曲线形状对但数值差一截。

原因:电容率 ε 和折射率 n 的关系是 n = √(εμ)。在非磁性介质中 μ ≈ μ₀,所以 n ≈ √ε。但很多人直接把折射率比值代入公式,忘了电容率比值是折射率比值的平方。比如 n₂/n₁ = 1.5,那 ε₂/ε₁ = 2.25,不是 1.5。

解决:代公式前先确认你手里的是电容率还是折射率。如果是折射率,先平方再代入。论文里所有 ε₂/ε₁ 都是电容率比值,不是折射率比值。

4.2 全反射时根号内负数导致程序报错

现象:入射角超过临界角后,Python 报RuntimeWarning: invalid value encountered in sqrt,曲线断掉。

原因:np.sqrt()对负数返回 NaN。全反射时 ε₂/ε₁ - sin²θᵢ < 0,根号内是负数。

解决:把sqrt_term的计算改成np.sqrt(eps_ratio - sin_i**2 + 0j),加一个复数单位,让 numpy 返回复数结果。反射系数取模值就是 1,透射系数取模值就是 0,物理上对应全反射。

4.3 布儒斯特角只对平行极化成立,别套到垂直极化上

现象:用垂直极化波找布儒斯特角,发现反射系数怎么都不为零,以为公式推错了。

原因:布儒斯特角是平行极化波特有的。垂直极化波的反射系数表达式里,分子是 cosθᵢ - √(...),这个式子只有在 cosθᵢ = √(...) 时才为零,而那个条件对应的是 ε₁ = ε₂,不是布儒斯特角。

解决:找布儒斯特角只用平行极化公式。垂直极化波没有无反射角,这是物理事实,不是计算错误。

4.4 斜滑投射时忽略相位反转

现象:做低空目标探测仿真时,地面反射波和目标回波叠加后信号完全抵消,以为是目标隐身了。

原因:斜滑投射时反射系数趋于 -1,反射波与入射波相位相反。如果地面反射波和目标回波路径差恰好对应半波长,两者相位相反、幅度相近,就会相消。

解决:做低空探测仿真时,地面反射系数必须取 -1 而不是 +1。如果忽略这个相位反转,仿真结果会完全错误。论文里明确写了“反射波与入射波相位相反”,这不是可选项,是必须项。

4.5 把透射系数和折射系数混为一谈

现象:论文摘要里出现了“折射系数”这个词,有人以为透射系数就是折射系数,直接混用。

原因:论文摘要中“反射系数和折射系数一般会变小”这里的“折射系数”实际指的就是透射系数,是同一篇文章里的用词不统一。但在其他文献里,“折射系数”可能指折射率,含义完全不同。

解决:看文献时先确认术语定义。本文中“折射系数”等同于“透射系数”,但换一篇文献就不一定了。自己写代码或报告时,统一用“透射系数”避免歧义。

5. 从公式到工程:极化匹配与雷达低空盲区的验证方法

5.1 用反射系数曲线做天线极化匹配检查

天线极化匹配是射频链路里的一个隐形杀手。发射天线是垂直极化,接收天线是水平极化,理论上接收功率为零——但实际中总有一些交叉极化分量,不会真的为零,但损耗可能高达 20 dB 以上。更隐蔽的情况是:发射和接收都是圆极化,但传播路径上有一个斜入射的介质界面,平行极化分量在布儒斯特角附近被削弱,到达接收端时极化状态已经偏了。

我一般会这样做:先确定传播路径上主要界面的入射角和电容率比值,用第 2 章的代码算出该角度下垂直和平行分量的反射、透射系数,然后看两个分量的幅度比和相位差。如果幅度比偏离 1 超过 3 dB,或者相位差偏离 90 度超过 20 度,圆极化就退化了,接收端需要调整极化匹配策略。

下面这段代码可以快速算出一个给定入射角和电容率比值下的极化退化情况:

import numpy as np def polarization_degradation(theta_deg, eps_ratio): """ 计算经过介质界面后圆极化的退化程度 theta_deg: 入射角,度 eps_ratio: eps2/eps1 返回:幅度比(dB),相位差(度) """ theta = np.radians(theta_deg) cos_i = np.cos(theta) sin_i = np.sin(theta) sqrt_term = np.sqrt(eps_ratio - sin_i**2 + 0j) # 垂直极化透射系数 T_perp = 2 * cos_i / (cos_i + sqrt_term) # 平行极化透射系数 T_para = 2 * np.sqrt(eps_ratio) * cos_i / (eps_ratio * cos_i + sqrt_term) amp_ratio_dB = 20 * np.log10(np.abs(T_para) / np.abs(T_perp)) phase_diff_deg = np.degrees(np.angle(T_para) - np.angle(T_perp)) return amp_ratio_dB, phase_diff_deg # 示例:入射角60度,光疏到光密 eps2/eps1=2.25 amp, phase = polarization_degradation(60, 2.25) print(f"幅度比: {amp:.2f} dB, 相位差: {phase:.2f} 度")

跑一下这个例子,你会看到幅度比和相位差都不是零,圆极化经过界面后变成了椭圆极化。如果这个界面在天线近场或者传播主路径上,链路预算里的极化损耗必须把这部分算进去,不能只算天线本身的交叉极化。

5.2 地面雷达低空盲区的量化估算

论文在结束语里提到,斜滑投射时反射系数趋于 -1,可以用来解释地面雷达为什么不能发现低空目标。这个解释是定性的,但工程上需要定量估算盲区范围。

思路是这样的:雷达波打到地面,反射波与直射波在目标处叠加。如果目标高度为 h,雷达高度为 h_r,目标距离为 R,地面反射路径与直射路径的路径差约为 ΔR ≈ 2h·h_r/R。当 ΔR = λ/2 时,反射波与直射波相位相反,叠加后场强抵消,目标处场强出现极小值。对应的目标高度就是盲区高度。

用这个公式可以快速估算:假设雷达高度 10 m,目标距离 10 km,波长 0.1 m(S 波段),盲区高度 h = λR/(4h_r) = 0.1×10000/(4×10) = 25 m。也就是说,25 m 以下的目标,直射波和地面反射波在目标处相消,雷达很难发现。

这个估算假设地面反射系数为 -1,实际地面不是理想导体,反射系数模值小于 1,盲区不会完全消掉,但会显著衰减。论文里没有给出这个估算公式,但这是从斜滑投射反射系数趋于 -1 这个结论直接推出来的工程用法。

5.3 一个我踩过的坑:忽略透射系数的相位

早期做多层介质透射仿真时,我只关注透射系数的模值,觉得相位不重要。结果做多层匹配的时候,每层厚度按模值优化,实际插损比仿真大了好几个 dB。后来才发现,透射系数的相位直接决定多层结构的干涉效果,模值匹配了但相位不对,层与层之间的反射波不能相消,插损自然下不来。

从那以后我每次做多层介质仿真,都强制把透射系数的相位一起输出,和模值放在同一张图里看。模值决定单层损耗,相位决定多层干涉,两个都不能丢。论文里的公式给的是复数表达式,取模只是其中一步,相位信息一直在里面,别自己把它扔掉。

希望这些从公式到代码再到踩坑的记录,能帮你在遇到介质界面极化问题时少走几步弯路。

本文还有配套的精品资源,点击获取

需要专业的网站建设服务?

联系我们获取免费的网站建设咨询和方案报价,让我们帮助您实现业务目标

立即咨询