電磁波工学
柴田幸司
平面波の反射と透過
左図に示す様に、平面波が境界面に対して垂直に入射する 場合を考える。この時の入射波をE i , H iとすると、入射波は 境界において、透過波E t, H tとE r, H rとに分解される。この 時の透過量をT,反射量をΓとおくと、領域2における媒質の 誘電率に対して透過量Tおよび反射量Γは決定される。よっ て、その量を先の電磁波の境界条件から求めてみる。垂直入射の場合
そのためまず透過係数は、 入射波と透過波の比として i t E E T と定義する。 次に反射係数は、入射波と反射波の比として i r E E と定義する。 Boundary Plane Medium (ガラスなど) Region1 Region2 z y Ei , Hi Er, Hr Et, Ht ε0 , μ0 ε, μ そして、これら透過係数や反射係数を計算する為のEyとHx成分を持つ平面波の各電 磁界成分は、マクスウェルの方程式より変形される次式の波動方程式の解としてEyの みを用いて次の様に置くことが出来る。 Ey 屈折率と誘電率 との関係 0 0 v c n c・・・媒質中の光速度 0r , 0r v・・・真空中の光速度0
2 2 0 0 2 2
t
z
E
E
t k z
j yi iE
e
E
0
t k z
j yr rE
e
E
0
t kz
j yt tE
e
E
t k z
j yi xi iH
Y
E
e
H
0 0
t k z
j yr xr rH
Y
E
e
H
0 0
t kz
j yt xt tH
YE
e
H
理由
x yH
E
0 0
スカラ表現されたマクスウェル方程式の解より電界と磁界には なる関係が あるので、H
xE
y 0 0
・・・
(1)
・・・
(2)
・・・
(3)
・・・
(4)
・・・
(5)
・・・
(6)
となる。ここで、自由 空間インピーダンスは o oY
Z
1
0 0
0 0 0
k
ここで、k0および kは自由空間および媒質中の波数であり、周波数と誘電率によって 決定される波長と関連する定数としてマクスウェルの方程式から波動方程式を導出す れば、その一般解の構成要素として,
k
を得る。 現時点では未知数であり、境界条件を考慮することにより定数として決定される。 Ey と直交している
sin j cos ej であり 波動として振舞うY
Z
1
また である。 Γ、Tを求める方針 1. (1)~(6)式に接線成分の連続条件を適用して 2. 各領域の電磁界の和が等しい性質より、2つの 方程式を組み立てる。 3. そして、Γを求める時には透過電界であるEt を消去し Tを求める時には反射電界であるErを消去して関係式を導出する。 よって与式は y xY
E
H
0 となり、(4)および(6)式は直ちに求まる。 なお、磁界の反射波は境界で更に向きが反転して(5)式のごとく( - )が無くなる。 ここで、誘電体境界面において、電磁界には
H
H
K
n
1
2
0
E
1E
2n
2 1n
E
E
n
・・・ (7) ・・・ (8) ・・・ (9) なる関係が存在する。境界条件について
異なる媒質の境界(不連続部)では、マクスウェルの方程式はそのまま適用
出来ない。 →境界面に境界条件を適用
媒質内における 電磁波の伝搬 Boundary Plane Medium (ガラスなど) Region1 Region2 Region1 Region2 ε1 , μ1 , σ1 ε2 , μ2 , σ2 n t τ ⊿ C S 小さな面Sとこれを囲む閉路C(⊿C=0)につ いて拡張されたマクスウェルアンペア(1.4) 式を適用
H
H
K
n
1
2
・・・
(1.16)
K ・・・ 境界面を面電流が流れて いる時の電流密度 [A/m](1.16)式
・・・
境界面に面電流が流れている場合、その分の不連続が生ずる。
・・・
(1.4)
変位電流 線積分
S Cdt
dS
d
d
s
J
D
n
H
電流Jにより 有限の値同様に小さな面Sとこれを囲む閉路C(
⊿C=0)
についてファラデー(1.7)式を適用
n
E
1
E
2
0
・・・
(1.17)
→
電界の接線成分は連続
領域1および2共に完全導体でなければ、面電流=0となる。
0
H
1H
2n
0
E
1E
2n
・・・
(1.20)
不連続部における
一般的な境界条件
(1.20)式を変形すれば
n
H
1
n
H
2(1.20)式の意味
x
y
z
n
Hx Hyとなる。これは、領域1および2の単位法線ベクトル
nに垂
直な磁界の2成分(H
xおよびH
y)が、境界面では等しいと
いうことを意味している。
S Cdt
dS
d
dt
d
d
s
B
n
E
・・・
(1.7)
誘電体の境界条件
境界面 媒質2 (ガラスなど) 領域1 領域2 z y 媒質1 (空気)入射波
透過波
媒質の誘電率に関連して透過
H
H
K
n
1
2
K ・・・ 境界面を面電流が流れて いる時の密度 [A/m]境界において、磁界はKの不連続が生じ
る。但し、完全金属以外ではK=0となり
0
E
1E
2n
となる。すなわち、電界の接線成分は連
続(境界での電界の接線成分は等しい)
2 1n
E
E
n
2 1n
H
H
n
が成り立つ。また、電界は
アンペアマクスウェルの積分形より
ファラデーの積分形より
n
z
y
x
E
yE
x 1E
n
とは
n…単位法線ベクトル
(境界面と垂直なベクトル)
1E
n
とは
nに垂直なE
yおよびE
x成分の
電界のこと
2 1n
E
E
n
について
つまり、境界において
2 1 y yE
E
2 1 x xE
E
領域名
が成り立つ
領域1
領域2
すなわち、(8)式は境界において電界の接線成分が等しいことを表しており 2 1 x x
E
E
,
となる。次に、(9)式のKは境界が完全導体では無い時には零であり
0
H
1H
2n
より
n
H
1
n
H
2 となるから 結局、磁界の接線成分も等しくなり,
を得る。 2 1 y yE
E
2 1 x xH
H
H
1y
H
y2接線成分の連続条件の導入
よって、領域1の電磁界の和と領域2の電磁界の和が等しいことより t r iE
E
E
,
H
i
H
r
H
t も成り立つ。 この条件より、Eyに関して
j
t kz
yt z k t j yr z k t j yie
E
e
E
e
E
0
0
・・・ (12) が成り立つ。・・・
(10)
・・・
(11)
また、Hxに関してはH
xi
H
xr
H
xt より
j
t kz
yt z k t j yr z k t j yie
Y
E
e
Y
E
e
E
Y
0 0 0 0 を変形して
j
t kz
yt z k t j yr z k t j yie
E
e
Y
E
e
E
Y
0
0
0 ・・・ (13) となる。(12)と(13)式について z=0, t=0の場合にはe0=1となるので、eの項を省略すれば改めて yt yr yi
E
E
E
・・・
(14)
yt yr yiE
YE
E
Y
0(
)
・・・
(15)
反射係数の導出
これらの式より、まず反射係数ΓはEyiとEyrとの比であるから透過電界であるEytを消 去することを考える。すなわち、(14)式を(15)式に代入すれば
)
(
)
(
0E
yiE
yrY
E
yiE
yrY
となるので、これを整理すれば yr yr yi yiY
E
Y
E
Y
E
E
Y
0
0
→(
)
(
)
0 0Y
E
Y
Y
Y
E
yi
yr
と変形される。これより、反射係数は
0 0 0 0 0 0 0 0 0 0Y
Y
Y
Y
E
E
E
E
yi yr i r を得る。・・・
(17)
・・・
(16)
と書き直せる。i t
E
E
T
透過係数Tは入射と透過の比であるから、(14)および(15)式からEyrを消去すること を考える。まず、(15)式を yi yt yrYE
Y
E
E
Y
0
0
→Y
0E
yr
Y
0E
yi
YE
yt と変形する。一方、(14)式は yi yt yrE
E
E
・・・ (18) となるので これを(18)式に代入すれば
yt yi yi ytE
Y
E
Y
E
E
Y
0
0
を得る。 これを展開すれば yt yi yi ytY
E
Y
E
YE
E
Y
0
0
0
yi yi yt ytY
E
Y
E
Y
E
E
Y
0
0
0
yi ytY
Y
Y
E
E
0
2
0
より となる。 ここで、透過係数は と定義されているのでy方向の電界について透過係数の導出
となり、さらに展開して0 0 0 0 0 0 0
2
2
・・・
(19)
であり、この式にYおよびY0 を代入して与式は を得る。また、たとえばガラスなどの誘電体を仮定してμ=μ0とすると 0 0 0 0 0 0 0 0
0 0 0 0 0 0 02
2
T
となり、これが2媒質による平面波の反射係数及び透過係数の公式である。Y
Y
Y
E
E
E
E
T
yi yt i t
0 02
Boundary Plane Medium (ガラスなど) Region1 Region2 z y
E
i, H
iE
r, H
rE
t, H
tε
0ε
そこで、具体的に領域2の媒質の比誘電率が与え られた時の反射係数について考える。もし、媒質 2の比誘電率がεr =1の場合には0
1
1
1
1
1
1
0 0 0 0
r r r r
となるから次式を得る。1
1
1
2
1
2
2
0 0 0
r rT
また、媒質2の比誘電率がεr =2の場合には171
.
0
414
.
2
414
.
0
1
414
.
1
414
.
1
1
1
2
2
1
829
.
0
414
.
2
2
414
.
1
1
2
2
1
2
T
となる。
r
0
Γ
、Tの計算例
等価回路との関係
Boundary Plane Medium (ガラスなど) Region1 Region2 z yE
i, H
iE
r, H
rE
t, H
tε
0ε
交流電圧源 電源抵抗z=0
Z
Z
0 なお、この様なモデルは等価回路に置き換 えることが出来て、ΓおよびTはそれぞれ 0 0 0 0 02
1
1
1
2
2
Z
Z
Z
Z
Z
Z
Y
Y
Y
T
0 0 0 0 0 01
1
1
1
Z
Z
Z
Z
Z
Z
Z
Z
Y
Y
Y
Y
となる。つまり、インピーダンス不連続 問題に帰着する。 領域1 領域2同様の手法により、厚みを持った媒質による透過波も計算が可能である。 Ei , Hi Et, Ht z=0 z=d z=dでの境界条件から領域2と3との電磁界の関係式を組み立て、領域1と2 に関する関係式も含めて4元連立方程式からEtおよびHtを導出 → ガラス窓などによる光の透過係数の計算が可能 Region1 Region2 Region3
Boundary Plane Medium (ガラスなど) Region1 Region2 z y 左図に示す様に、自由空間ε0中の電界および磁界が εなる誘電媒質に対して入射面と水平に平面波がθiな る角度で入射する場合を考える。 この時の入射波はyz面に対して時間項は省略してy方 向およびz方向成分の合成として
E
i, H
iE
r, H
r
j k ysin i zcos i i i oˆ
cos
ˆ
sin
e
E
y
z
1 iE
ysin i zcos i k j oˆ
e
Z
E
1 1x
Η
i と表現出来る。ここで、E0 は電界強度の初期値である。 ・・・ (2a1) ・・・ (2a2)電磁波の媒質への斜め入射
水平偏波の場合
θi θrE
t, H
t θt Ezi Eyi Hxiε
0ε
z方向成分 y方向成分電界が入射断面であるyz面に平行して振動
E0 i cos ˆzE0 i sin ˆyE0 つまり、水平偏波の場合には入射面であるyz面に対し て電界Eは平行でありEy およびEz が存在する。一方、 磁界はフレミングの左手の法則よりHx のみが存在する。 さらに、境界面からの反射波をとおく。この式は反射量の表現式であるから、zの符号のみを逆にすることを留意する。ここで k1およびZ1は領域1における波数および特性インピーダンスであり、空気(真空)中において 0 0 0 1
c
k
0 0 1
Z
,
e
sin
ˆ
cos
ˆ
E
j k ysin r zcos r r r o
y
z
1 rE
,
xˆ
( sin cos ) 1 1 y r z r k j oe
Z
E
rH
・・・
(2a3)
・・・
(2a4)
となる。また、θrは反射角でありΓは境界における反射係数である。 さらに、領域2への透過波を,
角周波数を真空中の光の速度で割った値 とおけば、T は透過係数である。ここで、k2 およびZ2 は領域2における波数および特性イン ピーダンスとして・・・
(2a5)
,
)
sin
ˆ
cos
ˆ
(
j k2 (ysin t zcos t) t t oT
e
E
y
z
E
t.
xˆ
( sin cos ) 2 2 y t z t k j oT
e
Z
E
tH
・・・
(2a6)
Z方向のみ向きが変化する0 1 2
k
1 0 2
Z
で与えられる。なお、現時点においてはΓおよびT 、さらにθr , θtは未知数である。 以上の様に定義した電磁界について、z=0点における接線成分の連続性より 0 ) 2 ( 0 ) 1 (
y z z yE
E
・・・ (2a7) 0 ) 2 ( 0 ) 1 (
x z z xH
H
・・・ (2a8) を適用する。その際に入射、反射および透過波の関係を考慮すれば、垂直入射の時と同様 0 0 0
yr z
yt z z yiE
E
E
および 0 0 0
xr z
xt z z xiH
H
H
なる関係を得る。 さらに、z=0における境界面に平行な成分について上記の式に(2a1)~(2a6)式を代入すれ ば、zに関する関数はすべて消え以下の2式を得る。,
cos
cos
cos
1 sin i 1 sin r j k2 ysin tt y k j r y k j i
e
e
T
e
.
1
sin 2 sin 1 sin 1 2 1 1 y i j k y r j k y t k je
Z
T
e
Z
e
Z
・・・ (2a10) ・・・ (2a9)(2a9)および(2a10)式の両辺はy座標の関数であり、 この時点では未知数であるΓおよびTは定まらな い(つまり、Γ、Tは一定とならない)。 ここで、先程の境界条件が成り立つということは、 図に示す様にz=0における境界面に沿ってy= y1 ,y2 ,y3 ・・・yn において同じ状態にあることである。 つまり、この2つの式がすべてのy座標について 成り立つ為には指数関数の変数同士が相等しい 必要があり、その為に次式を満足する必要がある。
,
sin
sin
sin
1 2 1y
ik
y
rk
y
tk
・・・
(2a11)
このyに関する式の解は、一般にスネルの反射および屈折に関する法則として、次式で表さ れることが知られている。,
r i
・・・
(2a12)
.
sin
sin
2 1 ik
tk
・・・
(2a13)
.
sin
sin
2 1k
k
i t
周波数は関係無 い(光でも電波 でも同じ) 上記の式は、(2.9)および(2.10)式の位相項がいかなるyに対しても境界の両面にて同じ比 で変化することを表しており、しばしばこれを位相の整合条件と呼ぶ。 ここで、θt に着目すると よりi t
sin
k
k
sin
2 1
より、透過波の角度は
i tsin
k
k
sin
2 1 1 となる。 式(2a9)および(2a10)に(2a12)および(2a13)の条件を適用し、両式からΓまたはTだけに ついて整理すれば、以下の2媒質における反射および透過係数を得る。 t i i j k y t y k j i y k j ie
e
T
e
1 sin 1 sin 2 sincos
cos
cos
まず、(2a9)式に i r
を代入して
i j k y t t y k j i ie
T
e
1 sin 2 sincos
cos
cos
1
cos
cos
cos
sin sin 2 1
i i t y k j y k jT
e
e
t i
・・(2a11a) (y=0の場合について) となるので、これを 整理すれば y=0の場合には 0 1 0 e e となるから となる。一方、(2a10)式にy=0を代入して両辺に-を掛けると 2 1 11
Z
T
Z
Z
・・(2a11b) となるので、これをTについて整理すれば
1
1 2Z
Z
T
・・(2a11c) となる。また、Γにて整理すれば2 1 1
1
Z
T
Z
Z
,
2 1 11
Z
T
Z
Z
・・(2a11d) となる。そこで、これらの関係よりまず反射係数Γを求めることを考える。すなわち、 Tを消去するために(2a11c)式を(2a11a)式に代入して
t i iZ
Z
cos
cos
cos
1
1 2
i i t tZ
Z
Z
Z
cos
cos
cos
cos
1 2 1 2
t i i tZ
Z
Z
Z
cos
cos
cos
cos
1 2 1 2
i t t iZ
Z
Z
Z
cos
cos
cos
cos
1 2 1 2
t i i t t i i tZ
Z
Z
Z
Z
Z
Z
Z
cos
cos
cos
cos
cos
cos
cos
cos
2 1 1 2 1 2 1 2
となる。次に、透過係数TはΓを消去するために(2a11d)式を(2a11a)式に代入して i i t
Z
T
Z
Z
T
cos
cos
1
cos
2 1 1
t i iT
iZ
Z
T
cos
cos
cos
cos
2 1 i i t
T
Z
Z
T
cos
2
cos
cos
2 1
i i tZ
Z
T
cos
cos
2
cos
2 1
i t i i t iZ
Z
Z
Z
Z
T
cos
cos
cos
2
cos
cos
cos
2
1 2 2 2 1
を得る。 からi
t
i
t
Z
Z
Z
Z
cos
cos
cos
cos
1
2
1
2
・・・ (2a14)i
t
i
Z
Z
Z
T
cos
cos
cos
2
1
2
2
・・・ (2a15) となる。ここで、電磁界が垂直入射する場合には θi =θr = θt=0であるから 1 2 1 2Z
Z
Z
Z
1
2
2
2
Z
Z
Z
T
となり、先の垂直入射における 結果と一致する。 以上をまとめれば この偏波ではブリュスター角( θb )という特別な角度が存在し、反射係数はΓ=0な ることが知られており、この条件について考えてみる。(2a14)式にてΓ=0となる為には 分子が零になる必要があるので,
cos
cos
1 2 tZ
bZ
を得る。さらに入射角をブリュスター角と してあらたにθi =θb と置き換えれば0
cos
cos
1 2
t
Z
i
Z
・・・*1 b tZ
Z
cos
cos
2 1 となる。 なる関係を得る。よって*1式をcosθtについて整理すると ・・・*4 一方、三角関数の 公式である t t
2sin
1
cos
となる。ここで、(2a13)式のスネルの法則は1
sin
cos
2
t
2
t
なる関係をcosθt について整理すると ・・・*2 よって、この式をθb のみの関係式に変形することを考える。そのために以下の手順 を実行する。 1.三角関数の公式よりcosθt をsinθt で表現し 2.スネルの法則よりこれをsinθi で表現し 3.θi =θb としてθb で表現し 4.このcosθt を*4式に代入する,
sin
1
cos
2 2 2 2 1 b tk
k
t ik
sin
sin
k
1
2
であるから、これをsinθt について整理すれば t ik
k
sin
sin
2 1
を得るので、θi =θb なる関係より※2式にこの式を代入すれば ・・・*3 を得る。これより、*3に*4式を代入してcosθtを消去すれば となる。よって、,
sin
1
cos
2 2 2 2 1 2 1 b bk
k
Z
Z
0
sin
1
cos
2 2 2 2 1 2 1
b bk
k
Z
Z
を満たす様な b
がブリュスター角となる。 |Γ|1.0
0.2
θ
90
0
b
垂直偏波の場合
垂直偏波の場合には電界のベクトルは入射面のyz断面に対して垂直となる。まず、入射する 電磁界は,
)
sin
ˆ
cos
ˆ
(
( sin cos ) 1 1 y i z i k j i i oe
Z
E
y
z
H
i,
ˆ
j k1 (ysin i zcos i) oe
E
x
E
i ・・・ (2b1) とおくことが出来る。ここで、k1およびZ1は領域1における波数、特性インピーダンスであり 0 0 1
k
0 0 1
Z
となる。これより、反射および透過係数は ・・・ (2b2),
)
sin
ˆ
cos
ˆ
(
( sin cos ) 1 1 y r z r k j r r oe
Z
E
y
z
H
r ・・・ (2b3),
ˆ
j k1 (ysin r zcos r ) oe
E
x
E
r ・・・ (2b4) ・・・ (2b5),
ˆ
j k2 (ysin t zcos t ) oT
e
E
x
E
t,
)
sin
ˆ
cos
ˆ
(
( sin cos ) 2 2 y t z t k j t t oe
Z
T
E
y
z
H
t ・・・ (2b6) とおく。ここで、ΓおよびT は反射および透過係数であり、領域2について 0 2
k
0 2Z
となる。なお、現時点においてはΓおよびT 、さらにはθr , θtは未知数である。 以上の様に定義した電磁界について、Z=0点における電磁界の接線成分の連続性 (境界条件)より 0 ) 2 ( 0 ) 1 (
x z z xE
E
・・・ (2b7) 0 ) 2 ( 0 ) 1 (
y z z yH
H
・・・ (2b8) を適用すれば、先ほどと同様に以下の2式を得る。,
sin sin sin 1 2 1 y i j k y r j k y t k je
T
e
e
r i j k y r y k j ie
Z
e
Z
sin 1 sin 1 2 1cos
cos
1
・・・ (2b9) ・・・ (2b10).
cos
sin 2 2 y t k j te
Z
T
水平偏波の時と同様に、ExおよびHyがZ=0点においてすべてのyに対して連続
である為には、位相の整合条件より、以下のスネルの法則が成立する。
,
sin
sin
sin
1
2
1
i
k
r
k
t
k
・・・
(2b11)
式(2b9)および(2b10)に(2a11)式の条件を適用することにより、それぞれ、
ΓおよびTについて整理すれば以下の2媒質における反射および透過係数を
得る。
t
i
t
i
Z
Z
Z
Z
cos
cos
cos
cos
1
2
1
2
・・・
(2a12)
t
i
i
Z
Z
Z
T
cos
cos
cos
2
1
2
2
・・・
(2a13)
,
cos
cos
12 i
Z
tZ
となる。また、スネルの法則より
sin
sin
sin
,
2 1 1 i