矢量分析
矢量代数
标量:只有大小,没有方向的量
矢量:既有大小,又有方向的量
对于矢量A,我们常用eA表示与其同方向的单位矢量
矢量的加减运算遵循平行四边形定则,此处不做赘述
矢量的乘积运算有两种,点积和叉积
对于点积有
A⋅B=ABcosθ
对于叉积有
A×B=enABsinθ
点积之后的结果是一个标量,而叉积之后的结果是一个矢量,其方向遵循右手螺旋定则
矢量的叉积运算一般使用行列式进行计算
A×B=exAxBxeyAyByezAzBz
对于点积和叉积分别有标量三重积A⋅(B×C)和矢量三重积A×(B×C)
标量三重积有以下运算性质
A⋅(B×C)=B⋅(C×A)=C⋅(A×B)
矢量三重积有以下运算性质
A×(B×C)=B(A⋅C)−C(A⋅B)
三种常用的正交曲线坐标系
直角坐标系
比较熟悉,最常用的坐标系
三个坐标变量的变化范围是−∞<x,y,z<∞
位置矢量是
r=exx+eyy+ezz
其微分为
dr=exdx+eydy+ezdz
面元和体积元分别为
dSx=dydz,dSy=dxdz,dSz=dxdy
dV=dxdydz
圆柱坐标系
三个坐标变量的变化范围是0≤ρ<∞,0≤ϕ≤2π,−∞<z<∞
与直角坐标系的变换关系为
ρ=x2+y2,ϕ=arctan(y/x),z=z
或
x=ρcosϕ,y=ρsinϕ,z=z
位置矢量为
r=eρρ+ezz
其微分元矢量为
dr=eρdρ+eϕρdϕ+ezdz
面积元和体积元分别为
dSρ=ρdϕdz,dSϕ=dρdz,dSz=ρdρdϕ
dV=ρdρdϕdz
球坐标系
三个坐标变量的取值范围是
0≤r<∞,0≤θ≤π,0≤ϕ≤2π
与直角坐标系之间的变换关系
r=x2+y2+z2,θ=arccos(z/x2+y2+z2),ϕ=arctan(y/x)
或
x=rsinθcosϕ,y=rsinθsinϕ,z=rcosθ
其位置矢量是
r=err
微分元矢量为
dr=erdr+eθrdθ+eϕrsinθdϕ
面积元和体积元分别为
dSr=r2sinθdθdϕ,dSθ=rsinθdrdϕ,dSϕ=rdrdθ
dV=r2sinθdrdθdϕ
标量场的方向导数和梯度
标量场
若我们研究的物理量是一个标量,则由其所确定的场称为标量场
一个标量场可以用一个标量函数来表示,即
u=u(x,y,z,t)
或者由位置矢量来表示
u=u(r,t)
对于和时间无关的标量场,则可以表示为
u=u(r)
对于任意给定的常数c,方程u(r)=c就是等值面方程
等值面方程具有以下特点
- 常数c取不同的值,就可以得到不同的等值面,从而可形成等值面族
- 标量场的等值面族会充满场所在的整个空间
- 由于标量函数u(r)是单值的,一个点只能在一个等值面上,所以等值面之间互不相交
方向导数
设M0是标量场u(M)中的一点,从点M0出发引一条射线l,点M是射线l上的动点,到点M0的距离为Δl
当点M沿射线l趋近于M0时,比值Δlu(M)−u(M0)的极限称为标量场u(M)在点M0处沿l方向的方向导数,即
∂l∂u∣M0=Δl→0limΔlu(M)−u(M0)
方向导数的物理含义就是标量场u(M)在点M0处沿l方向的变化率
当∂l∂u>0时,标量场u(M)沿l方向是增加的;反之则是减小的
如果∂l∂u=0,则无变化
在直角坐标系当中,若射线l与x,y,z轴的夹角分别为α,β,γ,则有
dldx=cosα,dldy=cosβ,dldz=cosγ
式子当中的cosα,cosβ,cosγ为l的方向余弦,可得直角坐标系当中的方向导数的计算公式
∂l∂u=∂x∂ucosα+∂y∂ucosβ+∂z∂ucosγ
梯度
在直角坐标系当中,梯度的具体表达式如下
grad u=ex∂x∂u+ey∂y∂u+ez∂z∂u
此处我们定义∇为哈密顿算符,在直角坐标系当中有
∇=ex∂x∂+ey∂y∂+ez∂z∂
所以则有
grad u=∇u
在圆柱坐标系和球坐标系当中,梯度的计算公式为
∇u=eρ∂ρ∂u+eϕρ∂ϕ∂u+ez∂z∂u
∇u=er∂r∂u+eθr∂θ∂u+eϕrsinθ∂ϕ∂u
由以上分析可知,标量场的梯度有以下特点
- 通常称∇u为标量场u产生的梯度场
- 标量场u在给定定点M处沿任意方向el的方向导数等于该点的梯度∇u在方向el上的投影
- 标量场u在点M处的梯度垂直于过该点的等值面,且指向u(M)增加的方向
通量和散度
一个矢量场F可以被分解为三个标量场,即F(x,y,z)=exFx(x,y,z)+eyFy(x,y,z)+ezFz(x,y,z)
对于矢量场F(r),可以用一些有向曲线来形象的描述矢量在空间的分布,即矢量线
设M是矢量线上任意一点,位置矢量是r,则其微分矢量dr在点M处与矢量线相切,即
dr×F=0
则可推出矢量线方程
Fxdx=Fydy=Fzdz
在矢量场F当中,矢量穿过面元矢量的通量可定义为F⋅dS,将各个面元相加即可得到曲面S的通量
Φ=∫SF⋅dS
如果S是一个闭合曲面,则可表示为
Φ=∮SF⋅dS
在矢量场F中任意一点M作一个包围该点的任意闭合曲面S,当S以任意方式收缩到点M时,所限定的体积ΔV将趋于0
如果极限ΔV∮SF⋅dS存在,则将此极限称为矢量场F在点M处的散度,并记作div F,即
div F=ΔV→0limΔV∮SF⋅dS
⎩⎨⎧div F>0div F<0div F=0正通量源负通量源无通量源
散度定理
散度定理即高斯定理,它表明矢量场的散度在任意体积V上的体积分等于矢量场穿出限定该体积的闭合曲面S的通量,即
∫V∇⋅FdV=∮SF⋅dS
矢量场的环流和旋度
环流和环流面密度
矢量场F沿场中的一条有向闭合路径C的曲线积分
Γ=∮CF⋅dl
称为矢量场沿闭合路径C的环流
对于矢量场当中任意一点M,选取一个方向矢量en,并且以其为法向矢量作一面元矢量ΔS,边界为有向闭合路径C,若极限ΔS→0limΔS∮CF⋅dl存在,则称该极限为矢量场在点M处沿方向en的环流面密度,记作rotn F,即
rotn F=ΔS→0limΔS∮CF⋅dl
由此定义可知,环流面密度不仅和点M的位置有关,还和面元矢量的法向矢量有关
矢量场的旋度
矢量场F在点M处的旋度定义为一个矢量,以符号rot F表示,它在点M处沿方向en的分量等于矢量场在点M处沿方向en的环流面密度,即
en⋅rot F=rotn F
易知当en和rot F的方向相同时,en⋅rot F的值最大为∣rot F∣
因此,rot F的方向是使矢量场在M点处取得最大环流面密度的方向,其模值等于最大的环流面密度
直角坐标系,圆柱坐标系和球坐标系的旋度的计算公式分别为
rot F=ex(∂y∂Fz−∂z∂Fy)+ey(∂z∂Fx−∂x∂Fz)+ez(∂x∂Fy−∂y∂Fx)=∇×F
rot F=eρ (ρ1∂ϕ∂Fz−∂z∂Fϕ)+eϕ (∂z∂Fρ−∂ρ∂Fz)+ez ρ1[∂ρ∂(ρFϕ)−∂ϕ∂Fρ]
rot F=errsinθ1[∂θ∂(sinθFϕ)−∂ϕ∂Fθ]+eθr1[sinθ1∂ϕ∂Fr−∂r∂(rFϕ)]+eϕr1[∂r∂(rFθ)−∂θ∂Fr]
斯托克斯定理
矢量场F的旋度∇×F在曲面S上的面积分等于矢量场在限定曲面的闭合曲线C上的线积分,即
∫S∇×F⋅dS=∮CF⋅dl
无旋场的标量位
如果一个矢量场的旋度处处为0,即
∇×F≡0
则称该矢量场为无旋场,它是由散度源产生的
利用斯托克斯定理可知,无旋场沿闭合路径C的环流等于0,也就是说,无旋场内的曲线积分与路径无关,只和终点有关
标量场的梯度有一个重要的性质,即它的旋度恒等于0
∇×(∇u)≡0
也就是说,任何标量场的梯度场一定是无旋场
一个标量场的梯度可以由它的梯度完全确定
对于一个旋度处处为0的矢量场F,总可以把它表示为某一标量场的梯度,即如果∇×F≡0,则存在标量函数u,使得
F=−∇u
则函数u称为无旋场F的标量位函数,简称标量位,负号是为使其电磁场中电场强度E和标量电位φ的关系一致
设无旋场F的散度为
∇⋅F=g
其中g为已知的标量函数,可得
∇⋅(∇u)=−g
其中∇⋅(∇u)是对标量场u的二阶微分运算,称为拉普拉斯运算,记为∇⋅(∇u)=∇2u
∇2称为拉普拉斯算符,方程
∇2u=−g
称为标量泊松方程,当g=0时称为拉普拉斯方程
无散场的矢量位
如果一个矢量场F的散度处处为0,即
∇⋅F≡0
则称该矢量场为无散场
由散度定理可得,无散场通过任何闭合曲面S的通量等于0,即
∮SF⋅dS=∫v∇⋅FdV=0
矢量场的散度有一个重要的性质,就是旋度的散度恒等于0,即
∇⋅(∇×A)=0
也就是说,任何矢量场的旋度场一定是无散场
对于一个散度处处为0的矢量场F,总可以表示为某一矢量场的旋度,即如果∇⋅F≡0,则存在矢量函数A,使得
F=∇×A
函数A称为无散场F的矢量位函数,简称矢量位
定义矢量场的拉普拉斯运算为
∇2A=∇(∇⋅A)−∇×(∇×A)
方程
∇2A=−G
称为矢量泊松方程,在G=0的区域,则有
∇2A=0
称为矢量拉普拉斯方程
格林定理
格林第一恒等式定义如下
∫V(u∇2v+∇u⋅∇v)dV=∮Su∇v⋅endS
根据方向导数和梯度的关系,有
u∇v⋅en=u∂n∂v
将格林第一恒等式进行变换,可得
∫V(u∇2v−v∇2u)dV=∮S(u∂n∂v−v∂n∂u)dS
此方程称为格林第二恒等式
亥姆霍兹定理
在有限的区域V内,任一矢量场由它的散度,旋度和边界条件(场在空间上的分布)唯一地决定,且可表示为
F(r)=−∇u(r)+∇×A(r)
亥姆霍兹定理表明了
- 矢量场可以由一个标量函数的梯度和一个矢量函数的旋度之和来表示
- 一个矢量场可以表示为一个无散场和一个无旋场之和
- 如果区域V内矢量场的散度和旋度处处为0,则F由其在边界面S上的场分布完全确定
- 对于无界空间,只要矢量场满足∣F∣∝1/∣r−r′∣1+δ(δ>0),此时矢量场可以由其散度和旋度完全确定
只有在矢量场连续的区域内,∇⋅F和∇×F才有意义,因为它们都包含对于空间坐标的导数
电磁场的基本规律
电荷守恒定律
基本点电荷的电量为e=1.602×10−9C
电荷密度
- 电荷体密度
单位体积内所含电荷量的大小
- 电荷面密度
当电荷分布区域的厚度趋于0时,可以将这种电荷分布看作电荷面密度
在电荷体密度有限的表面上,其电荷面密度必然为0
在电荷面密度不为0的表面上,其电荷体密度的值则为无穷大
- 电荷线密度
当电荷连续分布于横截面积可以忽略的细线上时,电荷体密度和电荷面密度都趋于无穷大,此时需要用电荷线密度进行描述
- 点电荷
当带电体的尺寸远远小于观察点至带电体的距离时,带电体的形状和其中的电荷分布已经无关紧要了,可以抽象为一个带电荷的点,即点电荷
电流密度
- 电流密度矢量J(r,t)
在垂直于电荷运动的方向上取一个面积元enΔS,en为正电荷运动的方向,如果流过面积元的电流为Δi(r,t),则定义电流密度矢量为
J(r,t)=enΔS→0limΔSΔi(r,t)=endSdi(r,t)
- 面电流密度矢量JS(r,t)
若电荷在一个厚度h可以忽略的薄层是运动,所形成的电流称为面分布电流
定义面电流密度矢量为
JS(r,t)=endldi(r,t)
- 线电流
当电荷在一个横截面可以忽略的细线中运动时,可以把电流看作在一根无限细的线上流动,称为线电流
i=ρlv
其中ρl为线电流密度,v是电流速度
- 电流元
对于线电流i,沿电流流动的方向取一个线元矢量dl,将idl称为线电流的电流元
电荷守恒定律和电流连续性方程
电荷守恒定律:在一个与外界没有电荷交换的系统内,正负电荷的代数和在任何物理过程中始终保持不变
电流连续性:单位时间内从闭合面内流出的电荷量等于闭合面所限定的体积内的电荷减少量,即
∮SJ⋅dS=−dtd∫VρdV
此为电流连续性方程的积分形式,经过变换可得
∇⋅J+∂t∂ρ=0
此为电流连续性方程的微分形式
对于恒定电流而言,有∇⋅J=0,即恒定电流场是一个无散度的场,从任何闭合曲面穿过的恒定电流为0
真空中静电场的基本规则
库仑定律
在电场中放入一个试验电荷q0,设其受到的电场力为F,可定义该点的电场强度矢量E为
E=q0F
库仑定律的方程为
F12=e124πε0R122q1q2=4πε0R123q1q2R12
其中en是由q1指向q2的单位矢量,R12=r2−r1是由q1到q2的距离矢量,ε0=36π1×10−9F/m,称为真空的介电常数
静电场的电场力的叠加遵循平行四边形原则
从而可以推导出电场强度表达式为
E(r)=4πε0R3qR
其中R=r−r′是源点到场点的距离矢量
静电场的散度和旋度
静电场的高斯定理的微分形式为
∇⋅E=ε0ρ
积分形式为
∮SE⋅dS=ε01∫VρdV
静电场的旋度为
∇×E=0
所以静电场是无旋场
真空中恒定磁场的基本规律
安培力定律
若一个电量为q0的电荷以速度v在磁场中运动,它所受到的磁场力为
Fm=q0v×B
因此可得磁感应强度的大小为
B = \lim_{q_0\rightarrow 0}\frac{F_m|_\max}{q_0v}
可得电流元Idl在磁场中受到的磁场力为
dFm=Idl×B
对于电流元产生的磁场,毕奥-萨伐尔定律为
B(r)=4πμ0∮CR3Idl′×R
其中R=r−r′,R=∣r−r′∣
对于体电流密度为J(r′)的体分布电流,则磁感应强度为
B(r)=4πμ0∫VR3J(r′)×RdV′
同样,对于面电流密度为Js(r′)的面分布电流,磁感应强度为
B(r)=4πμ0∫SR3JS(r′)×RdS′
恒定磁场的散度和旋度
恒定磁场的散度为
∇⋅B(r)=0
可得恒定磁场为无散场
恒定磁场的旋度为
∇×B(r)=μ0J(r)
此为安培环路定理的微分形式,其积分形式为
∮CB(r)⋅dl=μ0I
媒质的电磁特性
电介质的极化特性
空间的电场强度E是自由电荷产生的电场强度E0与极化电荷产生的电场强度Ep的叠加,即
E=E0+Ep
极化强度矢量P定义如下
P=ΔV→0limΔV∑ipi
其中pi=qili表示体积ΔV中第i个分子的电偶极矩,极化强度矢量的单位是C/m2
若电介质内各点的P相同,则称为均匀极化,否则为非均匀极化
无论是均匀极化还是非均匀极化,电介质的表面都会出现面分布的极化电荷
但是对于均匀极化,电介质内部不会出现体分布的极化电荷,而非均匀极化会出现体分布的极化电荷
电介质表面上的极化电荷面密度为
ρSP=P⋅en
我们定义电位移矢量为
D=ε0E+P
其单位为C/m2
电介质中高斯定理的积分形式为
∮SD⋅dS=q
它表明电位移矢量穿过任一闭合面的通量等于该闭合面内总的自由电荷,而与极化电荷无关
电介质中高斯定理的微分形式为
∇⋅D=ρ
它表明电介质内任一点的电位移矢量的散度等于该点的自由电荷体密度,电位移矢量线从正的自由电荷出发,终于负的自由电荷
虽然电位移矢量的通量和散度都只与自由电荷有关,但不能认为电位移矢量只由自由电荷产生而与极化电荷无关
实际上电位移矢量是由自由电荷和极化电荷共同产生的电场强度矢量以及极化强度矢量所决定的
因此,电介质中静电场的基本方程的积分形式为
⎩⎨⎧∮SD⋅dS=q∮CE⋅dl=0
对于相同的电场强度E,不同电介质的极化情况是不相同的
大多数常见的电介质在电场的作用下发生极化时,其极化强度矢量与电场强度矢量成正比,表示为
P=ε0χeE
其中χe称为电介质的电极化率,是一个正实数
这类电介质的极化强度矢量P与电场强度矢量E的方向相同,电极化率和电场强度的方向无关,称为各向同性电介质
电极化率和电场强度大小无关的电介质称为线性介质,电极化率的值随电场强度的大小的不同而变化的电介质称为非线性介质
在电场强度不是很强的情况下,多数介质都是线性的,但是在强电场的作用下就可能呈现出非线性特性
对于各向同性介质,有
D=ε0E+χeε0E=εE
其中ε=ε0εr称为电介质的介电常数,εr=1+χe为电介质的相对介电常数
上述公式表示了各向同性电介质的本构关系,此关系表明在各向同性电介质中,电位移矢量D和电场强度E的方向相同
若电介质是线性的,则介电常数·与电场强度的大小无关,电位移矢量与电场强度的大小成正比
若电介质是均匀的,则介电常数处处相等,若是非均匀的,则介电常数是空间坐标的函数
对于各向异性电介质,电位移矢量和电场强度的方向通常是不相同的
磁介质的磁化特性
电子在自己的轨道上以恒定速度绕原子核运动,形成一个环形电流,它相当于一个磁偶极子,将其磁矩称为自旋磁矩
通常可以忽略原子的自旋,将磁介质的每个分子(或原子)等效于一个环形电流,称为分子电流/束缚电流,分子电流的磁矩称为分子磁矩,用pm表示,定义为
pm=iΔS
其中i为分子电流,ΔS=enΔS为分子电流所围的面积元矢量,其方向en与i流动的方向成右手螺旋关系
不同磁介质的分子磁矩不同,可能为0,可能不为0
磁介质可以分为顺磁体(分子磁矩不为0),抗磁体(分子磁矩为0)和铁磁体(磁畴)三种
描述磁介质磁化程度的物理量称为磁化强度矢量,用M表示,定义为单位体积内分子磁矩的矢量和
M=ΔV→0limΔV∑ipmi
其中pmi表示体积Δv内第i个分子的磁矩,M是一个矢量函数,单位是A/m,若磁介质当中各点M相同,称之为均匀磁化,否则称为非均匀磁化
磁介质被磁化后,其内部和表面可能出现宏观电流分布,称为磁化电流,可以表示为磁化电流密度的积分,即
IM=∫SJM⋅dS
JM=∇×M
磁化面电流密度矢量为
JSM=M×en
其中en为磁介质表面的法向单位矢量
磁介质被磁化后,磁感应强度等于传导电流I在真空中产生的磁感应强度和磁化电流在真空中产生的磁感应强度的叠加,即
B=B0+BM
安培环路定理的积分形式为
∮CH⋅dl=I
它表明磁场强度沿任意闭合路径的环量等于与该闭合路径交链的传导电流
其中H为磁场强度,如下所示
H=μ01B−M
由斯托克斯定理可得安培环路定理的微分形式
∇×H=J
由此可见,在磁介质中任一点的磁场强度等于该点的传导电流密度
磁化强度矢量M与磁场强度矢量H和磁介质的物理特性有关
实验表明,对于线性和各向同性磁介质,磁化强度和磁场强度成正比,表示为
M=χmH
其中χm称为磁介质的磁化率,不同的磁介质有不同的磁化率
由上式可得
B=(1+χm)μ0H=μrμ0H=μH
此式称为各向同性磁介质的本构关系,其中μ=μrμ0称为磁介质的磁导率,μr称为磁介质的相对磁导率
真空中χm=0,μr=1无磁化效应,M=0,B=μ0H
μr>1,χm>0的磁介质是顺磁体,μr<1,χm<0的磁介质是抗磁体
顺磁体和抗磁体的磁化效应都很弱,对于铁磁体,B和H的关系是非线性的,μr值可以很大
导电媒质的传导特性
导电媒质内部有许多能自由移动的带电粒子,它们在外电场的作用下可以定向移动形成电流
在线性和各向同性的导电媒质内,任意一点的电流密度矢量J与该点的电场强度E成正比,表示为
J=σE
其中σ称为媒质的电导率,上式是线性各向同性导电媒质的本构关系,也是欧姆定律的微分形式
电导率的值和媒质有关,且随温度变化
焦耳定律的微分形式为
pL=dVdPL=J⋅E
其中pL为单位体积的损耗功率,J=ρv,其中ρ为电荷体密度,v为电荷的平均速度
焦耳定律的积分形式为
PL=∫VpLdV=∫VJ⋅EdV
对于线性和各向同性的导体,焦耳定理的微分和积分形式可分别表示为
pL=σE⋅E=σE2
PL=∫VσE2dV
电磁感应定律和位移电流
法拉第电磁感应定律
实验表明,感应电动势与穿过回路所围面积的磁通量的时间变化率成正比,即
Ein=−dtdΦ=−dtd∫SB⋅dS
感应电流产生的磁通总是对原磁通的变化起阻碍作用
在空间中的总电场由感应电场和库伦电场的叠加
E=Ein+EC
由于有∮CE⋅dl=0,则有
∮CE⋅dl=−dtd∫SB⋅dS
这就是推广之后的法拉第电磁感应定律的积分形式
如果回路是静止的,则穿过回路的磁场会随时间变化,此时有
∮CE⋅dl=−∫S∂t∂B⋅dS
这就是静止回路位于时变磁场当中的法拉第电磁感应定律的积分形式
而可由斯托克斯定理推导出其微分形式
∇×E=−∂t∂B
法拉第电磁感应定律积分形式的一般形式为
∮CE⋅dl=−∫S∂t∂B⋅dS+∮C(v×B)⋅dl
其微分形式为
∇×E=−∂t∂B+∇×(v×B)
位移电流
之前的∇×H=J不适用于时变电磁场,需要进行修正
时变条件下的电流连续性方程为
Jd=∂t∂D
其中Jd称为位移电流密度
修正后的公式如下
∇×H=J+∂t∂D
时变电磁场当中安培环路定理的积分形式如下
∮CH⋅dl=∫S(J+∂t∂D)⋅dS
需要说明的是位移电流没有通常电流的概念,由D=ε0E+P可知
∂t∂D=ε0∂t∂E+∂t∂P
其中,∂t∂P是由于介质分子的极化电荷随时间变化而引起的,称为极化电流,而ε0∂t∂E则仅表示电场随时间的变化,并不对应任何电荷的运动
麦克斯韦方程组
麦克斯韦方程组的积分和微分形式
麦克斯韦方程组的积分形式为
∮CH⋅dl=∫SJ⋅dS+∫S∂t∂D⋅dS(1)
∮CE⋅dl=−∫S∂t∂B⋅dS(2)
∮SB⋅dS=0(3)
∮SD⋅dS=∫VρdV(4)
式(1)的含义是磁场强度沿任意闭合曲线的环量等于穿过以该闭合曲线为周界的任意曲面的传导电流和位移电流之和
式(2)的含义是电场强度沿任意闭合曲线的环量等于穿过以该闭合曲线为周界的任一曲面的磁通量变化率的负值
式(3)的含义是穿过任意闭合曲面的磁感应强度的通量恒等于0
式(4)的含义是穿过任意闭合曲面的电位移的通量等于该闭合面所包围的自由电荷的代数和
麦克斯韦方程组的微分形式为
∇×H=J+∂t∂D(5)
∇×E=−∂t∂B(6)
∇⋅B=0(7)
∇⋅D=ρ(8)
式(5)表明传导电流和时变电场要产生磁场,都是磁场的涡旋源
式(6)表明时变磁场要产生电场,是电场的涡旋源
式(7)表明磁场是无散场,磁感应线是闭合曲线
式(8)表明电荷要产生电场,是电场的散度源
由麦克斯韦方程组可以推导出恒定磁场方程和静电场方程如下
∇×H=J
∇⋅B=0
∇×E=0
∇⋅D=ρ
媒质的本构关系
媒质的本构关系也称为电磁场的辅助方程,在线性各向同性的导电媒质当中,本构关系为
D=εE
B=μH
J=σE
其中J表示在电场作用下在导电媒质中形成的传导电流密度
若媒质是均匀的且参数不随时间变化,则可以得到以下用E,H表示的麦克斯韦方程组
∇×H=Je+σE+ε∂t∂E
∇×E=−μ∂t∂H
∇⋅H=0
∇⋅E=ερ
其中Je为外加的电流密度
准静态电磁场
在时变的情况下,如果电磁场随时间变化非常缓慢时,若可以忽略∂t∂B或∂t∂D时,则这种电磁场称为准静态电磁场
它们的特点是电场和磁场都随时间变化,但是却具有一些静态场的性质
根据忽略的∂t∂B或∂t∂D的不同,可以分为电准静态电磁场和磁准静态电磁场
- 电准静态电磁场
此时麦克斯韦方程组可以近似为
∇×H=J+∂t∂D
∇×E≈0
∇⋅B=0
∇⋅D=ρ
可以利用静电场的公式去求解E和D
- 磁准静态电磁场
此时麦克斯韦方程组可以近似为
∇×H≈J
∇×E=−∂t∂B
∇⋅B=0
∇⋅D=ρ
可以利用恒定磁场的公式去求解B和H
电磁场的边界条件
边界条件的一般形式
设两种媒质的分界面为光滑曲面,分界面上任一点P处的法向量en由媒质2指向媒质1
一般表达式为
en×(H1−H2)=JS
en×(E1−E2)=0
en⋅(B1−B2)=0
en⋅(D1−D2)=ρS
其中JS为分界面上的传导面电流密度,ρS为分界面上的自由电荷面密度
若两种媒介的电导率均为有限值,则JS=limΔh→0JΔh=0,即分界面上不存在面电流分布,可得
en×(H1−H2)=0
当分界面上没有自由电荷面密度时,有
en⋅(D1−D2)=0
理想导体表面上的边界条件
设媒质2为理想导体,则在时变情况下,E2,D2,B2,H2均为0,则有
en×H1=JS
en×E1=0
en⋅B1=0
en⋅D1=ρS
理想介质分界面上的边界条件
设媒质1和媒质2为两种不同的理想介质,由于理想介质当中不存在自由电子,则有
en×(H1−H2)=0
en×(E1−E2)=0
en⋅(B1−B2)=0
en⋅(D1−D2)=0
静态电磁场及其边值问题的解
静电场分析
静电场的基本方程和边界条件
静电场的基本方程的积分形式为
⎩⎨⎧∮SD⋅dS=∫VρdV∮CE⋅dl=0
微分形式为
⎩⎨⎧∇⋅D=ρ∇×E=0
在线性,各向同性的电介质中,本构关系为
D=εE
静电场的边界条件的一般形式为
en×(E1−E2)=0或E1t=E2t
en⋅(D1−D2)=ρS或D1n−D2n=ρS
电位函数
- 电位和电位差
电场强度矢量E可以表示为标量函数φ的梯度,即
E(r)=−∇φ(r)
其中的φ(r)称为静电场的电位函数,简称为电位
电位函数不是唯一的,可以任意加上一个常数,电场强度E保持不变
可以得到电位的计算公式:在线性,各向同性的无界均匀媒质中,点电荷q的产生的电场的电位函数为
φ(r)=4πεRq+c
其中c为任意常数
点P,Q之间的电位差的物理意义是:把一个单位正电荷从点P沿任意路径移动到点Q的过程中,电场力所做的功
若选定Q点为电位参考点,则P点的电位为
φ(P)=∫PQE⋅dl
若场源电荷分布在有限区域,通常选定无限远处为电位参考点,即
φ(P)=∫P∞E⋅dl
- 静电位的微分方程和边界条件
静电位满足标量泊松方程
∇2φ(r)=−ερ(r)
在无自由电荷分布的区域中,即ρ=0,则电位函数满足拉普拉斯方程
∇2φ(r)=0
在静电场当中,由于导体的内部的电场为0,所以导体为等位体,导体表面为等位面
因此,在导体表面上的电位边界条件为
⎩⎨⎧φ=常数ε∂n∂φ=−ρS
导体系统的电容
电容器的电容定义为电荷量和电压之比
C=Uq
单位为法拉(F)
一个孤立导体可以看作是与无穷远处的另一个导体构成的电容器,其电容大小和电荷量以及电位有关,为
C=φq
在计算电容器的电容时,常常采用以下的计算步骤
- 假定两导体上分别带电荷+q和−q
- 根据假定的电荷求出电场强度E
- 由∫12E⋅dl求电位差U
- 求出比值C=Uq
设空间中存在N+1个导体,各导体所带的电荷量分别为q0,q1,q2,⋯,qN且所有导体的电荷量之和为0
则称该多导体为静电独立系统
若空间中的介质是线性的,则各导体的电位与各导体上的电荷之间也是线性关系
如果把电荷量为q0的导体作为电位参考点,则其余N个导体的电位φ1,φ2,⋯,φN可表示为
⎩⎨⎧φ1=α11q1+α12q2+⋯+α1NqNφ2=α21q1+α22q2+⋯+α2NqN⋯φN=αN1q1+αN2q2+⋯+αNNqN
其中αij称为电位系数,αii为自电位系数,αij为互电位系数
在计算αij的时候,可以令qj=0,而其余导体的电荷都为0,即
αij=qjφi(q1=⋯=qN=0)
所有的电位系数都具有对称性,即αij=αji
电容系数和部分电容暂且不写
静电场的能量
电场能量的公式为
We=21CU2=21qU
对于线性和各向同性介质,有D=εE,则有
We=21∫VεE2dV
静电力
Fi=−∂gi∂We(q为常数,dgi为虚位移)
导电媒质中的恒定电场分析
恒定电场的基本方程和边界条件
- 恒定电场的基本方程
由麦克斯韦方程组可以得出恒定电场满足的方程为
∇⋅D=ρ
∇×E=0
电场强度沿任一闭合路径的线积分恒为0,即
∮CE⋅dl=0
由于恒定电场不随时间变化,所以电荷的分布也不随时间变化,即∂t∂ρ=0
根据电流连续性方程可得
∮SJ⋅dS=0
相应的微分形式为
∇⋅J=0
在线性,各向同性的导电媒质当中,本构关系为
D=εE
J=σE
恒定电场也可以用电位梯度来表示
E=−∇φ
在均匀导电媒质当中(σ=常数)当中,电位函数满足拉普拉斯方程,即
∇2φ=0
导电媒质当中的电荷体密度为
ρ=J⋅∇(σε)
由此可见,在恒定电场当中,只有非均匀的导电媒质内部才可能存在电荷分布,而均匀导电媒质内的电荷密度处处为0,电荷只能分布在媒质的分界面上
- 恒定电场的边界条件
恒定电场的电位函数的边界条件是
φ1=φ2
σ1∂n∂φ1=σ2∂n∂φ2
恒定磁场分析
恒定磁场的基本方程和边界条件
由麦克斯韦方程组可以得出恒定磁场的基本方程的积分形式为
⎩⎨⎧∮CH⋅dl=∫SJ⋅dS∮SB⋅dS=0
微分形式为
⎩⎨⎧∇×H=J∇⋅B=0
这表明恒定磁场是无源有旋场,恒定电流是产生恒定磁场的涡旋源,磁力线是与源电流相交链的闭合曲线
线性,各向同性的磁介质的本构关系为
B=μH
在分界面上的B是连续的
矢量磁位
电磁场的位函数是根据电磁场性质引入的辅助函数,其数学定义与电磁场性质相关
B=∇×A(r)
恒定磁场中采用的规范条件为库伦规范
∇⋅A(r)=0
规范条件是人为引入的限定条件
矢量位的泊松方程为
∇2A=−μ0J
在无源区,J=0,则有矢量位的拉普拉斯方程
∇2A=0
矢量磁位可以通过矢量分解来求解
Ax=4πμ0∫VRJxdV
Ay=4πμ0∫VRJydV
Az=4πμ0∫VRJzdV
电感
磁通的为通过面积S的磁感应强度B
Φ=∫SB⋅dS=∮CA⋅dl
磁链是指与所有电流回路铰链的总磁通
单匝细线圈的磁链Ψ=Φ
多匝细线圈密绕的磁链Ψ=nΦ
粗导线回路的磁链Ψ=Ψo+Ψi
其中Ψo为外磁链,Ψi为内磁链
电流回路的电感是穿过某电流回路的磁链与回路中电流强度之比
L=IΦ
若某回路C载流为I,其产生的磁场穿过回路自身C所形成的自感磁链为,则定义回路C的自感系数为
L=IΨ(单位为H)
自感磁链是回路载流I在回路内自已产生的
回路自感仅与回路自身的几何形状、尺寸和媒质磁导率
有关,与回路中载流无关
若回路为N匝线圈密绕,则L = N^2L_0\quad\text{L_0为单匝线圈电感}
若回路导线直径较粗,则回路自感为内自感和外自感之和
L=Lo+Li
其中Lo为回路外自感,Li为回路内自感
回路C1与回路C2交链的磁通量为Ψ21,则回路C1对C2的互感系数为
M21=I1Ψ21
同理,回路C2对C1的互感系数为
M12=I2Ψ12
M12=I2Ψ12=4πμ∮C1∮C2R21dl2⋅dl1
互感的大小只与回路几何性质、相对位置和周围介质有关
M21=M12=M
恒定磁场的能量
若电流为体电流分布,则其在空间中产生的磁能为
Wm=21∫VJ⋅AdV
其中A为体电流J在dV处产生的磁位
V为整个空间
- 上式只适用于恒定磁场
- 被积函数21J⋅A不代表能量密度
对于N个回路系统,i回路自感为Lii,i回路与j回路间互感为Lij,i回路电流为Ii ,则磁回路系统的磁场能量为
Wm=21i=1∑Nj=1∑NLijIiIj
磁场的能量密度为
wm=21B⋅H
对应的磁场能量为
Wm=21∫VB⋅HdV
积分区域为电场所在的整个空间
时变电磁场
波动方程
电场强度矢量和磁场强度矢量满足的波动方程为
∇2E−με∂t2∂2E=μ∂t∂J+ε1∇ρ
∇2H−με∂t2∂2H=−∇×J
在无源区,ρ=0,J=0,则波动方程为
∇2E−με∂t2∂2E=0
∇2H−με∂t2∂2H=0
电磁场的位函数
矢量位和标量位
因为
∇⋅B=0
则可以将磁场B表示为一个矢量函数A的旋度,即
B=∇×A
其中的A则称为电磁场的矢量位
同样的,可以用一个标量函数φ的梯度来表示电场强度矢量E
通常规定矢量位A的散度为
∇⋅A=−με∂t∂φ
此式称为洛伦兹条件
达朗贝尔方程
在洛伦兹条件下,矢量位和标量位满足的微分方程称为达朗贝尔方程
∇2A−με∂t2∂2A=−μJ
∇2φ−με∂t2∂2φ=−ε1ρ
电磁能量守恒定律
电磁能量密度:单位体积中电磁场的能量,为电场能量和磁场能量之和
w=we+wm=21(εE2+μH2)
体积V内总能量为
W=∫V(21E⋅D+21H⋅B)dV
坡印廷定理物理意义:单位时间内流入体积V内的电磁能量等于体积V内增加的电磁能量与体积V内损耗的电磁能量之和
坡印廷定理微分形式:
−∇⋅(E×H)=E⋅J+∂t∂(21E⋅D+21H⋅B)
积分形式:
−∮S(E×H)⋅dS(流入体积V的电磁功率)=dtd∫V(21E⋅D+21H⋅B)dV(体积V内增加的电磁功率)+∫VE⋅JdV(体积V内损耗的电磁功率)
坡印廷矢量(电磁能流密度矢量)
坡印廷矢量的物理含义:通过垂直于能量传输方向单位面积的电磁功率(功率流密度)
S(t)=E(t)×H(t)
上式与时间有关,故也称瞬时坡印廷矢量
平均坡印廷矢量:将瞬时形式坡印廷矢量在一个周期内取平均
Sav=T1∫0TS(t)dt=T1∫0TE(t)×H(t)dt
瞬时坡印廷矢量反映某时刻的电磁能量流动情况
平均坡印廷矢量反映一个时间周期内的电磁能量传递情况
平均坡印廷矢量与时间无关
时谐电磁场
时谐电磁场是场量按正弦规律变化的电磁场,也叫正弦电磁场,表达式为
E(r,t)=E0cos[ωt+ϕ(r)]
H(r,t)=H0cos[ωt+ϕ(r)]
时谐电磁场可以通过复数形式来分离电磁场的空间变化和时间变化
A(r,t)=A0(r)cos[ωt+ϕ(r)]=Re{A0(r)ejϕ(r)ejωt}=Re[A(r)ejωt]
场量的复数形式转换为实数形式的方法
E=E0ejφ⟶×ejωtE0ej(ωt+φ)⟶取实部E0cos(ωt+φ)
在时谐电磁场当中,对于时间的求导运算为
∂t∂⟷jω
∂t2∂2⟷−ω2
可得麦克斯韦方程组的复数形式
⎩⎨⎧∇×H(r)=J(r)+jωD(r)∇×E(r)=−jωB(r)∇⋅D(r)=ρ(r)∇⋅B(r)=0
补充:
∇⋅J(r)=−jωρ(r)
亥姆霍兹方程
在无源区,有
⎩⎨⎧∇2E+k2E=0∇2H+k2H=0
其中k=ωμε
时谐电磁场的能量和坡印廷矢量
在复数形势下,瞬时坡印廷矢量的表达式为
S(r,t)=E(r,t)×H(r,t)=Re[E(r)ejωt]×Re[H(r)ejωt]
平均坡印廷矢量的表达式为
Sav(r)=21Re[E(r)×H∗(r)]
均匀平面波在无界媒质中的传播
理想介质中的均匀平面波
电磁波场量的求解
通过波动方程求解会过于复杂,我们通常会采用适当建系的方法来简化运算
- 平面波的等相位面为x-y平面
- 平面波的传播方向为z轴方向
- 电场矢量方向沿坐标轴方向,如x向
这样就可以让E只随着z坐标变化而变化
方程
∇2E+k2E=0
的通解为
E=Em+cos(ωt−kz)+Em−cos(ωt+kz)
对于
E=exEm−ejkz
其相伴的磁场为
H=η1(−ez)×E
其中
η=εμ
理想介质中均匀平面波的电场与磁场相互正交,且相位相同
对于z向波的平面波函数,有以下结论
e−jkz表示从+z方向传播的均匀平面波函数
ejkz表示从−z方向传播的均匀平面波函数
无界理想媒质中均匀平面波的传播特性
波长是空间相位差为2π的两个波阵面的间距,即
kλ=2π⟶λ=k2π=fμε1
相位常数:表示波传播单位距离的相位变化
k=λ2π
k的大小等于空间距离2π内所包含的波长数目,因此也称为波数
波矢量k:大小等于相位常数,方向为电磁波传播方向
相速:电磁波的等相位面在空间中的移动速度
vp=dtdz=kω=με1
vp=fλ
由公式易知,相速仅与媒质特性有关
真空中电磁波相速为光速
媒质本征阻抗(波阻抗):电场幅度和磁场幅度比
η=εμ
真空(空气)的本征阻抗为120π≈377Ω
对于沿ek方向传播的均匀平面电磁波,有
H=η1ek×E
E=ηH×ek
E,H,k三者相互垂直,且满足右手螺旋关系
电磁波的能量密度:w=εE2=μH2
电磁波的能流密度为
S(t)=E(t)×H(t)=ekη1E02cos2(ωt±k⋅r)
Sav=21Re[E×H∗]=2η1E02k
电场能量密度等于磁场能量密度
电场与磁场同相,电场振幅是磁场的η倍
均匀平面波在无界导电媒质中的传播
导电媒质中存在欧姆损耗,电介质被极化后存在极化损耗
通常用损耗角正切表征介质损耗特性,对导电媒质有
tanδσ=ωεσ=∣Jd∣∣J∣
介质损耗角正切等于导电媒质中传导电流与位移电流之比
⎩⎨⎧tanδσ≪1tanδσ≈1tanδσ≫1良绝缘体普通导电媒质良导体
媒质导电性的强弱与频率有关
由于在导电媒质中本征阻抗为复数,磁场相位滞后于电场相位j21arctanωεσ
在损耗媒质中,相速的公式为
vp=β(ω)ω
色散效应:相速随频率改变而改变的现象
色散波:具有色散效应的波
导电媒质中的电磁波为色散波
在导电媒质中,电场能量密度为
we=2εExm2e−2αzcos2(ωt−βz)
磁场能量密度为
wm=2μ∣ηc∣2Exm2e−2αzcos2(ωt−βz)
电磁波的平均能流密度为
Sav=2∣ηc∣Em2e−2αzcosφk(Em为电场振幅)
导电媒质中均匀平面波的磁场能量大于电场能量
在波的传播过程中,电场与磁场的振幅呈指数衰减
在导电媒质中,相速不仅与媒质参数有关,而且与频率有关,为色散波
磁场能量大于电场能量
在良导体中, 磁场相位滞后电场相位4π
在良导体中存在趋肤效应:高频电磁波只能存在于良导体的表面层内
趋肤深度δ:电磁波穿入良导体中,当波的幅度下降为表面处振幅的e1时,波在良导体中传播的距离
对于良导体:δ=2πλ
在良绝缘体中,有:
- 均匀平面波的衰减小
- 相位常数近似等于理想媒质中的相位常数
- 电场和磁场之间相位差比较小
电磁波的极化
基本概念
电磁波的极化:表征空间定点处电磁波电场矢量的时变规律
波的极化描述方法:在电磁波传播空间定点处,电场强度矢量的终端端点随时间变化的轨迹形状
基本形势:线极化,圆极化,椭圆极化
电场矢量的矢端方程
电磁波的极化由电场矢量的矢端方程揭示的轨迹形状来判断
以沿+z方向传播的均匀平面波为例,有
∣E(0,t)∣=Exm2cos2(ωt+ϕx)+Eym2cos2(ωt+ϕy)
电场矢端方程为
Exm2Ex2+Eym2Ey2−ExmEym2ExEycosΔϕ=sin2Δϕ
其中Δϕ=ϕy−ϕx
- 矢端的时间变化规律, 决定于各分量幅度和初相的大小
- 任意极化均可以由线极化得到
电磁波极化方式的判断
电磁波的极化状态取决于Ex和Ey的振幅Exm,Eym和相位差Δφ=φy−φx
对于沿+z方向传播的均匀平面波:
- 线极化 Δφ=0(1,3象限),±π(2,4象限)
- 圆极化 Δφ=±π/2,Exm=Eym,取+则为左旋圆极化,取−则为右旋圆极化
- 椭圆极化 其他情况 0<Δφ<π,左旋;−π<Δφ<0,右旋
任何一个线极化波、 圆极化波或椭圆极化波可分解成两个线极化波的叠加
任何一个线极化波都可以表示成旋向相反、振幅相等的两圆极化波的叠加
任何一个椭圆极化波也可以表示成旋向相反、振幅不等的两圆极化波的叠加