电磁波的传播

电磁波在介质表面的反射和折射

均匀线性介质中的电磁波

波动方程

2Eεμ2Et2=0

特解

E=E0ei(kzωt)B=B0ei(kzωt)

性质

  • 相速度,等相面的运动速度

    krωt=const.


    v=dr0dt=ωk

  • 波速
    将特解代入波动方程,得 ω k 的关系

    k2=εμω2

    代入相速度方程,得

    v=ωk=1εμ=ε0μ0εμc

    εε0>1 μμ01 ,所以介质中的波速始终小于光速。

  • B E 的关系

    B=kωE=εμek×E=1vek×E

介质表面的折射和反射

这里写图片描述

  • z=0 为界面,上方是介电常数和磁导率分别为 ε,μ 的均匀线性介质,下方的介电常数和磁导率是 ε,μ
  • 因为一般的电磁波可以分解成平面波,所以只考虑平面波的情形。
  • 平面斜入射的平面波的方程
    E=E0ei(krωt)B=1vek×E

    其中 k=(kx,0,kz)
  • 折射波的方程
    E=E0ei(krωt)B=μεek×E
  • 反射波的方程
    E′′=E′′0ei(k′′rω′′t)B′′=μεe′′k×E′′
  • 根据 E 在界面切向的连续性,得
    E0xei(kxxωt)+E′′0xei(k′′xx+k′′yyω′′t)=E0xei(kxx+kyyωt)

    这个式子在任意 (x,y,t) 都成立,所以有
    kx=k′′x=kxky=k′′y=0ω=ω′′=ω

    所以说反射和折射不改变入射波的频率。
  • 反射波的 k
    根据

    k=μεω

    入射波和反射波的环境和频率都相同,所以
    k′′=k

    所以
    k2x+k2z=k′′2x+k′′2z

    因为 k′′x=kz ,所以
    |k′′z|=|kz|,k′′z=kz

  • 折射波的 k
    折射波和入射波在不同的介质中,波速不同,波长也不同。

    λλ=kk=μεμε

  • 定义折射率
    n=μεμ0ε0

    则有
    λλ=nn
  • 折射定律
    由于
    sini=kxk,sini=kxk

    所以
    sinisini=kk=nn


    nsini=nsini

反射波和折射波的强度

三个波的相因子已经相等,所以还需要确定反射波和折射波的振幅。根据边界条件
E 连续,得

n×(E0+E′′0E0)=0

D 连续,得
n(D0+D′′0D0)=0

H 连续,得
n×(H0+H′′0E0)=0

B 连续,得
n(B0+B′′0B0)=0

在线性介质中,
n×(E0+E′′0E0)=0n(εE0+εE′′0εE0)=0n×(1μk×E0+1μk′′×E′′01μk×E0)=0n(k×E0+k′′×E′′0k×E0)=0

考虑两种特殊情形

  1. E 垂直于入射面
    利用两个切向连续的方程解得

    E0E0=2μncosiμncosi+μncosiE′′0E0=μncosiμncosiμncosi+μncosi

  2. E 平行于入射面
    利用两个法向连续的方程解得

    E0E0=2μncosiμncosi+μncosiE′′0E0=μncosiμncosiμncosi+μncosi

一般可以认为 μ=μ=1 ,所以有垂直入射时

E0E0=2ncosincosi+ncosiE′′0E0=ncosincosincosi+ncosi=sin(ii)sin(i+i)

平行入射时
E0E0=2ncosincosi+ncosiE′′0E0=ncosincosincosi+ncosi=tan(ii)tan(i+i)

特殊情形

在此基础上讨论几个特殊情形

半波损失

上述的电场振幅都是复数,但是它们的比值都是实数。如果比值为负,说明相位发生突变,称为半波损失
可以看出,折射波不会发生半波损失,只有垂直偏振的波( \? )从光疏介质向光密介质中入射时才会发生半波损失。

平行偏振的反射波的消失

平行偏振的反射波

E′′0E0=tan(ii)tan(i+i)

所以当 i+i=0 的时候, E′′0=0 ,这时反射波不存在。这个特殊角称为 布儒斯特角,记作 iB 。结合折射定律,得
taniB=nn

全反射

当电磁波从光密介质向光疏介质入射时,折射角大于入射角。由于折射角不大于 π2 ,所以入射角存在一个临界值。下面考虑入射角超过这个临界值的情形。

当入射角大于临界值时,折射定律仍然成立,不过此时无法定义折射角,折射定律的形式为

nsini=nkxk=nkxk(kx=kx)

此时 kx>k ,因为 k2=k2x+k2z ,所以 kz 是纯虚数, kz=iτ ,所以折射波可以表示成
E=E0eτzei(kxxωt)

沿界面传播,且幅度随着透射深度衰减。

导电介质中的电磁波

电磁波在导体中会产生电流,产生焦耳热,导致波的耗散和衰减。

E=0×E=BtB=0×B=μ0j+εμEt=μ0σE+εμEt

由此得到的波动方程多了一个阻尼项

2EμσEtεμ2Et2=0

方程的特解仍然可以写成

E=E0ei(krωt)

代回方程,可以得到 k ω 的关系。
k2=k2x+k2y+k2z=εμω2+iμσω

相比于介质中的关系,这是一个复数关系,因此 ω k 中至少有一个分量是复数。

扰动电场产生的电磁波

如果 t=0 时导体内有一个扰动电场 E0(r) ,将电场用傅里叶变换展开,得

E0(r)=eikrE0k(r)dk

其中,每个 E0k(r) 是一个单频的平面波, k 为实数。它们都满足上面的推导,所以每个分量的 ω 都是复数。设 ω=ω0iω1 ,那么

Ek=E0kei(krωt)=E0keω1tei(krω0t)

这是一个振幅随时间衰减的波。

导体中的透射波

若电磁波从导体外面射入导体(假设垂直入射),那么和上节的结论相同

ω=ω=ω′′

所以 k 是复数。
不妨设透射波的 k=(0,0,kz) kz=kz0+iτ ,那么电波的方程
E=E0eτzei(kz0zωt)

波幅不随时间衰减,而随距离(透射深度)衰减,所以这仍然是一个稳态的波。 kz0 描述透射波的传播, τ 描述衰减的深度。


k2=k2x+k2y+k2z=εμω2+iμσω

两边的实部和虚部分别相等,得

k2z0τ2=εμω2kz0τ=12μσω

上式说明总有 kz0>τ ,分两种情况讨论

不良导体


τkz0

即衰减很弱,那么
σ2εω

σ 很小,这种情况对应着不良导体。此时解得
kz0=εμω=μεμ0ε0kτ=12μεσ

和不导电介质中基本相同。

良导体

如果 τkz0 ,那么

k2z0τ2kz0τ

所以
σ2εω

σ 很大,对应良导体的情形。此时解得
kz0=τ=12μσωεμω=μεμ0ε0k

波长比在不导电介质中短很多,且振幅衰减很快,穿透深度 1τ 和波长相当。

下面考虑振幅。设入射波是沿 x 方向的线偏振波,由E H 的连续性,并忽略 μ μ0 的差别,得到

E0+E′′0=E0k(E0E′′0)=kE0

解得
E′′0E0=kkk+kE0E0=2kk+k

对于良导体, k=k0(1+i) , k0k=12μσωμ0ε0ω=σ2ε0ω
所以 kk=σε0ωeiπ41
所以
E′′0E0=1E0E0=2ε0ωσ1

这说明入射波几乎不能进入良导体,几乎完全被反射了。

理想导体

σ 的导体称为理想导体,理想导体内没有电磁场。在导体的边界条件中

nE=Σεn×E=0nB=0n×B=μK

E B 都是指导体外面的电磁场。
所以导体外面的电场垂直于表面;磁场平行于表面;导体能通过调节表面电荷密度和面电流密度保证内部电磁场为0.

微波在波导管内的传播

高频电磁波的传播不能直接用导线,否则能量损失严重。微波常常用波导管来传输。现在讨论一个矩形波导管模型。
这里写图片描述

建模

设一个截面为矩形的波导管, x,y 方向的长度分别为 a,b ,在波导管中沿 z 轴方向传播的电磁波满足的方程

2E1c22Et2=0

由于理想导体的特性,边界条件

Ey|x=0,a=0Ex|y=0,b=0Ez|x=0,a;y=0,b=0

求解

方程的特解为

E=E0(x,y)ei(kzzωt)

其中
E0x=[αx1cos(kxx)+αx2sin(kxx)][βx1cos(kyy)+βx2sin(kyy)]E0y=[αy1cos(kxx)+αy2sin(kxx)][βy1cos(kyy)+βy2sin(kyy)]E0z=[αz1cos(kxx)+αz2sin(kxx)][βz1cos(kyy)+βz2sin(kyy)]

代入边界条件,得
E0x=[A1cos(kxx)+A2sin(kxx)]sin(kyy)E0y=sin(kxx)[B1cos(kyy)+B2sin(kyy)]E0z=C1sin(kxx)sin(kyy)kx=mπaky=nπak2x+k2x+k2x=ω2c2

A1,A2,,B1,B2,C1 都是待定的复常数。代入横波性条件 E=0 ,得
kx(A1sinkxk+A2coskxx)sinkyy+kysinkxx(B1sinkyy+B2coskyy)+ikzC1sin(kxx)sin(kyy)=0

对任意的点都成立,可以得
A2=0B2=0kxA1+kyB1ikzC1=0

所以
E0x=A1cos(kxx)sin(kyy)E0y=B1sin(kxx)cos(kyy)E0z=C1sin(kxx)sin(kyy)kx=mπaky=nπak2x+k2y+k2z=ω2c2

磁场

磁场 B 的形式也是

B=B0(x,y)ei(kzzωt)

根据 Bt=×E

得到
B0x=1ω(\?)(B1kz+iC1ky)sin(kxx)cos(kyy)B0y=1ω(A1kz+iC1kx)cos(kxx)sin(kyy)B0z=iω(B1kxA1ky)cos(kxx)cos(kyy)

可以看出,当给定 ω,m,n 之后, kz 也确定了,还需要 A1,B1,C1 中四个实参量,就可以得到确定的解。

横电/磁波型

可以引入两种基本波型作为基。满足

C1=0
的波作为一种基本波型,满足
B1kz=A1ky
的波作为另一种基本波型。前一种电矢量垂直于传播方向称为 横电型波 TEmn ,后一种磁矢量垂直于传播方向,称为 横磁型波 TMmn

频率和波导管尺寸的关系

由于

k2x+k2y+k2z=(mπa)2+(nπb)2+k2z=ω2c2

的限制,在波导管内要传播频率 ω 的波, m,n 都有上限,具体个数和波导管的尺寸有关。而对一根固定的波导管来说,因为 m,n 不能同时为 0 ,所以ω有下限,即
ωmin=cπl,l=max{a,b}

波长
λmax=2πcωmin=2l

所以波导管能传送的最大波长是 2l


本篇主要参考俞允强《电动力学简明教程》

  • 7
    点赞
  • 22
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
电磁波传播的MATLAB程序可以用于模拟和分析电磁波在不同介质传播,包括空气、水、地球等。下面是一个简单的电磁波传播的MATLAB程序示例: ```MATLAB % 定义常数 c = 3e8; % 光速 mu = pi*4e-7; % 磁导率 epsilon = 8.854e-12; % 介电常数 % 定义模拟参数 f = 1e9; % 频率 lambda = c/f; % 长 dx = lambda/10; % 空间间隔 dt = dx/c/2; % 时间间隔 nt = 200; % 模拟时间步数 nx = 200; % 空间步数 % 初始化场 E = zeros(nx, nt); % 电场 H = zeros(nx, nt); % 磁场 % 定义介质 sigma = 0.01; % 电导率 epsilon_r = 4; % 相对介电常数 % 计算介质参数 epsilon_eff = epsilon_r * epsilon; mu_eff = mu; % 初始化介质 sigma_eff = sigma * ones(nx, 1); epsilon_eff = epsilon_eff * ones(nx, 1); mu_eff = mu_eff * ones(nx, 1); % 计算系数 ca = (1 - sigma_eff*dt./(2*epsilon_eff))./(1 + sigma_eff*dt./(2*epsilon_eff)); cb = (dt./(epsilon_eff*dx))./(1 + sigma_eff*dt./(2*epsilon_eff)); da = (1 - sigma_eff*dt./(2*mu_eff))./(1 + sigma_eff*dt./(2*mu_eff)); db = (dt./(mu_eff*dx))./(1 + sigma_eff*dt./(2*mu_eff)); % 模拟电磁波传播 for i = 2:nt % 更新电场 E(1, i) = sin(2*pi*f*dt*i); % 定义源 E(2:end-1, i) = ca(2:end-1).*E(2:end-1, i-1) + cb(2:end-1).*(H(3:end, i-1) - H(2:end-1, i-1)); E(end, i) = E(end-1, i); % 更新磁场 H(1:end-1, i) = da(1:end-1).*H(1:end-1, i-1) - db(1:end-1).*(E(2:end, i) - E(1:end-1, i)); H(end, i) = H(end-1, i); end % 绘制结果 figure; subplot(2, 1, 1); imagesc(E); xlabel('时间步数'); ylabel('空间步数'); title('电场分布'); colorbar; subplot(2, 1, 2); imagesc(H); xlabel('时间步数'); ylabel('空间步数'); title('磁场分布'); colorbar; ``` 这个简单的程序模拟了电磁波在一个具有一定电导率和相对介电常数的介质传播,并且绘制了电场和磁场的分布情况。你可以根据自己的需要修改模拟参数和介质参数,以及增加其他功能。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值