电磁学乱七八糟的符号(四)
@(study)[Maxe, markdown_study, LaTex_study]
author:何伟宝
这里重点是一般传输规律和矩形波导,chapter6 电磁波的传输
纵向场量法
说白了也就是从麦克斯韦方程式抽象出我们需要的波动方程,流程如下:

矢量波动方程
在无源自由空间场量中(由麦克斯韦方程式):
∇2E+k2E=0(1.1)
∇2H+k2H=0(1.2)
在波导中,设电磁波沿着z轴传输:
E(x,y,z)=E(x,y)eγr(1.3)
H(x,y,z)=H(x,y)eγr(1.4)
其中有:
行波因子γ
在上一章说过他也是一个传播常数,当γ为实数时,信号衰减.虚数时信号传播,且波数为其虚部
矢量分解
这里有意地把纵横量分开了:
E=(axEx+ayEy)+azEz
H=(axHx+ayHy)+azHz
顺便把拉普拉斯算符∇也分开:
∇t2=(∂x2∂2+∂y2∂2)+∂z2∂2=∇xy2+∂z2∂2
标量波动方程
将矢量分解的三个方程先带入(1.3)(1.4)再代入(1.1)(1.2),只截取纵向分量得:
∇xy2Ez+(k2+γ2)Ez=0
∇xy2Hz+(k2+γ2)Hz=0
再将上式改写成(1.3)(1.4)形式:
Ez(x,y,z)=Ez(x,y)E−γz
Hz(x,y,z)=Hz(x,y)H−γz
考虑麦克斯韦方程的旋度式:
∇×E=−jωμH
∇×H=jωεE
联立上四式可得六个标量方程:
∂y∂Ez+γEy=−jωμHx(标量1)
−γEx−∂x∂Ez=−jωμHy(标量2)
∂x∂Ey−∂x∂Ex=−jωμHz(标量3)
千万不要慌,由麦克斯韦美好的对称性可以知道,我们只要算一对叉乘就可以了,由对称性:
∂y∂Hz+γHy=jωεEx(标量4)
−γHx−∂x∂Hz=jωεEy(标量5)
∂x∂Hy−∂x∂Hx=jωεEz(标量6)
纵横关系式
联立以上六式可得(解这个会有点痛苦,但是这不重要)纵横关系式:
Ex=−kc21(γ∂x∂Ez+jωμ∂y∂Hz)(e.x)
Ey=−kc21(γ∂y∂Ez−jωμ∂x∂Hz)(e.y)
Hx=−kc21(γ∂x∂Hz−jωμ∂y∂Ez)(h.x)
Ex=−kc21(γ∂x∂Hz+jωμ∂x∂Ez)(h.y)
其中:
kc2=k2+γ2
如果不用书本的表示方法的话,你会发现一点公式的美学…
自此,纵向常量法就介绍完成了.这里的重点在于纵横关系式
各种导波的一般传输特性
概述
这一小节解决的问题是,某种电磁波要在波导中传输的存在可能性问题.重点有TEM,TE,TM波等.并且提供假设各种波存在的时候,怎么求解波动方程的思路.
TEM横电磁波
还是回到我们熟悉的波动方程,可以把上面的纵横关系式:
∇xy2Ez+kc2Ez=0(波动1)
∇xy2Hz+kc2Hz=0(波动2)
显然这一节的教材安排是不合理的…因为在TEM波中:
Ez=0,Hz=0
显然代入纵横关系式中,全军覆没…所以分析横电磁波的时候不能采用纵向常量法得到直接表达式
这时候我们可以代入得到纵横关系式前面一点的关系式中:
kc2=0或γ2+k2=0(2.1)
∇xy2E(x,y)=0∇xy2H(x,y)=0(tem)
那么我们就可以知道,代入纵横关系式会凉凉的原因是,(tem)他看上去就是一个静态场的方程,用麦克斯韦旋度式便变成0了.
这也启发我们,在求解TEM波动方程的时候,之需要先算出导波的横向分布函数,再乘以纵向传播因子e−γz就可以得到波动方程了.而且并不是每一种波导都会有TEM模.
存在条件
首先说明的一点是:空心波导只能传输TM或TE波,不能传输TEM波,因为在无外源的无限长导体空管中不可能存在静电场
书上P175,结合来看吧…(懒得打字)
TEM传播常数和相速
由(2.1)可知
γ=α+jβ=jk=jωεμ
解得
α=0,β=ωεμ
所以相速为:
v=βω=εμ1
可以看出TEM模导行波是与频率无关的非色散波
TEM的波阻抗
由(标量2)和(标量6)并代入TEM的定义式:
γEx=jωμHy
γHy=jωεEx
代入γ=jωεμ得(注意,求解不是联立.只要用其中一条代入就行了)
ZTEM=HyEx=εμ=η
可以看出,ZTEM和频率是没有关系的.
所以:TEM模在任何频率下都能传播非色散横电磁波
TE nor TM
在TM波中,Ez̸=0和Hz=0.所以只需要代入(波动1),同理:
在TE波中,Hz̸=0和Ez=0.所以只需要代入(波动2)
存在条件
可以看出,无论是哪一种,kc2都不会等于0,所以:
γ2+k2̸=0
被称为波导中TM波和TE波的存在条件。
传播常数和截止频率
由传播因子e−jγz可以知道,在e−γz→0时,传播截止.这个时候有γ→1
所以有:
γ=kc2−ωc2εμ=0
解得:
fc=2πεμkc
其中,fc被称为截止频率或临界频率(c to cut),所以反过来求γ得:
γ=⎩⎨⎧jk1−(ffc)2=jβf>fckc1−(ffc)2=αf<fc
可以看出:
当f<fc时,传播因子变成了e−αz,是一个衰减型凋落场
当f>fc时,传播因子变成了e−jβz,表示一个传播型色散行波
相速和波导波长
当f>fc时,因为是一个色散波,我们可以来讨论一下他的相速,由上面:
β=k1−(ffc)2
所以有,相速:
vp=βω=1−(ffc)2v>v
波导内波导行波的波长称为波导波长:
λg=β2π=k2π1−(ffc)21=1−(ffc)2λ>λ
表明导行波是与频率有关的色散行波
波阻抗
TM波
由纵横关系式,结合tm波的特征可得:
Ex=−kc2γ∙∂x∂Ez
Ey=−kc2γ∙∂y∂Ez
Hx=kc2jωε∙∂y∂Ez
Ey=−kc2jωε∙∂x∂Ez
所以定义TM波的波阻抗为:
ZTM=HyEx=Hx−Ey=jωεγ
消去γ得:
ZTM=⎩⎨⎧η1−(ffc)2=RTM,f>fc−jωεkc1−(ffc)2=−jXcTM,f<fc
TE波
按照TM波的套路,代入Ez=0得:
ZTM=⎩⎪⎨⎪⎧η1−(ffc)21=RTE,f>fcjkcμω1−(ffc)21=jXcTE,f<fc
互易性
由上面可以得出:
ZTM∙ZTE=η2=(ZTEM)2
可以看到TE和TM波的波阻抗具有互易性
矩形波导
这里也是要做纵横关系式求解的最后一步,代入边界条件
由前面就可以知道,矩形波导不能传播TEM波
首先假设矩形波导的数学模型:

长a宽b壁导体
先上一张图辅助一下大家后面看边界条件的法向还是切向

TM(图的右边)
边界条件:
(∂x2∂2+∂y2∂2+kc2)Ez(x,y)=0
{Ez∣x=0=0,Ez∣x=a=0Ez∣y=0=0,Ez∣y=b=0
其中kc2=γ2+k2称为截止波数.
公式的意义是很明确的:
传播TM波的时候矩形波导的边界都没有电场强度
以下是我以为的原因(有异议可以评论,大家互相学习一下)
- 一个原因(一对边)在于,边界条件中,法向的电场强度连续,而理想导体内部没有电磁场
- 另一对边是因为,上一章说过的趋肤效应导致的,而由于是σ=∞所以就为0了
纵向解
由于我们想求的纵横关系式中,x和y是独立分开的.所以假设:
Ez(x,y)=X(x)Y(y)
代入波动方程并化成常微分方程得:
dx2d2X+kx2X=0
dx2d2Y+ky2Y=0
其中:kc2=kx2+ky2
显然特征方程的根是两个纯虚数,故设通解:
X(x)=Asinkxx+Bcoskxx
Y(y)=Csinkyy+Dcoskyy
分别代入边界条件可得(书上P176):
Ez(x,y)=E0sinamπxsinbnπy,m,n=1,2,3......
其中:E0=AC由激励源强度确定
大概的思路是先带入x=0和y=0那两条,算出B,D=0再代入剩下两条即可.
横向解
现在求出了Ez的表达式,显然,代入一般情况可得:
Ex=−kc2γ(amπ)E0cosamπxsinbnπy
Ey=−kc2γ(bnπ)E0sinamπxsinbnπy
Hx=kc2jωε(bnπ)E0sinamπxcosbnπy
Hy=kc2jωε(amπ)E0sinbnπxcosamπy
其中:
kc=γ2+k2=kx2+ky2=(amπ)2+(bnπ)2
由TE,TM的存在条件可以知道,当m=n=0时,方程无意义
TE(图的左边)
由于和TM是同一个套路,这里就直接给公式了:
边界条件
(∂x2∂2+∂y2∂2+kc2)Hz(x,y)=0
{Hz∣x=0=0,Hz∣x=a=0Hz∣y=0=0,Hz∣y=b=0
纵向解
Hz(x,y)=H0cosamπxcosbnπy,m,n=1,2,3......
横向解
Ex=kc2jωε(bnπ)H0cosamπxsinbnπy
Ey=−kc2jωε(amπ)H0sinamπxcosbnπy
Hx=kc2γ(amπ)H0cosbnπysinbmπx
Hy=kc2γ(bnπ)H0sinamπxsinbnπy
同理:m=n=0时,公式无意义
横场分布的物理特性
这里对应的是P178,下面列举出来只作复习回想用:
- 沿x,y的驻波性和z向的行波性
- 平面波的非均匀性
- 场的多模性
- 模式的兼并性
- 模式的阶次性
导波的纵场传输特性*
截止性(高通特性)
之前在一般传输特性就讲过这个问题,只是k可以由m和n给出,所以回代得:
kc=γ2+k2=kx2+ky2=(amπ)2+(bnπ)2
fc−2πεμkc=2εμ1(am)2+(bn)2
λc=kc2π=(am)2+(bn)22
色散性和滤波性
由上一个性质可以知道,在截取频率之前的波形都会因为传播常数的实部不为0而全部被去掉
所以当f>fc时(α=0):
β=ω2εμ−(amπ)2−(bnπ)2
λg=β2π=ω2εμ−(amπ)2−(bnπ)22π
vp=βω=ω2εμ−(amπ)2−(bnπ)2ω
阻抗双重性
这个由截止性就知道,低于截止频率的波阻抗呈阻性,高于的呈电抗性:
ZTM=jωεγ={ωε1ω2εμ−(amπ)2−(bnπ)2=RTM,f>fc−jωε1(amπ)2+(bnπ)2−ω2εμ=−jXcTM,f<fc
ZTE=γjωμ=⎩⎨⎧ωμ1ω2εμ−(amπ)2−(bnπ)21=RTM,f>fcjωμ(amπ)2+(bnπ)2−ω2εμ1=jXcTM,f<fc
主模TE10的传输特性
用主模传输的重点问题在于单模传输 单模传输 单模传输 单模传输
场分布
至于为什么TE10是主模的话,就不说了,你只要把 m,n的各个值代进去纵横关系式,就可以知道了
Ey=πωμaH0sinaπxcos(ωt−βz−2π)
Hx=πβaH0sinaπxcos(ωt−βz+2π)
Hz=H0cosaπxcos(ωt−βz)
…其他三个为0…
传输特性
根据前面说的那些,代入m=1,n=0得:
fc=2aεμ1
λc=2a
β=k1−(ffc)2=ω2εμ−(aπ)2
λg=β2π=k2π1−(ffc)21=ω2εμ−(aπ)22π
vp=βω=1−(ffc)2v=ω2εμ−(aπ)2ω
ZTE=η1−(ffc)21=ωμω2εμ−(aπ)21

结语
因为这里写了比较多的波动方程,所以会有点长!
如果你想请我吃个南五的话