電磁波工学
柴田幸司
第8回
電磁波の伝搬特性Ⅱ
本章の目的
・産業や通信に用いられる電磁波は宇宙的な規模での振る舞いを考えると、その周波 数によって空間を伝搬する性質などが異なる。 ・よって、特に電離層での振る舞いを例に、その違いについて理解する。 電離層伝搬に関連する周波数 MF(中波) ・・・ 300kHz~ 3MHz HF(短波SW) ・・・ 3MHz~30MHz上に行くほど電子密度が高い
電離層が等価的 な誘電体となる 超長波 3kHz~30kHz(波長:10~100km) 電離層の影響をあまり受けず、直接波とし て伝搬(水中も伝搬) LF 長波 30kHz~300kHz(波長1~10km) 船舶・航空機用ビーコン、標準電波(JJY) 昼間はD層で ほとんど吸収 短波(3MHz~30MHz) 中波 (300kHz~3MHz)~電磁波の電離層への垂直入射~
電離層(fp=√N)電磁波が電離層に垂直入射した場合、入射
させる周波数 f が電子サイクロトロンの周波
数 f
Hと衝突周波数νに比べて充分大きい時、
平面波に対する電離層の屈折率は
電子密度の 非常に濃い場所 (プラズマ状態) 2 2 2 0 21
1
f
f
m
N
e
n
p
で与えられる。ここで、諸定数は以下の通りである。
f : 送信する電磁波の周波数
[Hz]
f
p:その電子密度での電離層の臨界
(プラズマ)周波数[Hz]
N : 電離層の電子密度 [個/m
3]
e : 電子の電荷量 1.6×10
-19[C]
m : 電子の質量 9.109×10
-31[kg]
(SI単位系による)
ε
o: 真空の誘電率
8.854・10
-12短波帯の電波伝搬
あとで導出 電荷が沢山あるので誘電体 と考えることが出来る 1より小さい
2 2 29
1
81
1
f
N
f
N
n
f
2
であるから、先のプラズマの屈折率nにω=2πfを代入してみると
2 2 0 2 2 0 2 2 0 24
1
)
2
(
1
1
f
m
N
e
f
m
N
e
m
N
e
n
2 2 0 24
1
f
N
m
e
となる。さらに であるから与式に e2 o m 4 2 80.403 を代入して となる。そこで, を得る。つまり、fp 以上の周波数の 電磁波は電離層で反射されず突き ぬけることになる。nの変形によるfおよびf
pの分離
定数 2 21
f
f
n
p pf
N
9
とおけば e =1.6・10-19 m=9.109・10-31 ε0 =8.854・10-12N
9
をあらたに 電離層の臨界周波数として 1より小さい垂直入射では屈折率がn=0を満たす点で電磁波は全反射するから、与式にn=0を
代入すれば
2 2 9 1 0 f N となる。よって、この両辺を2乗すれば
2 29
1
0
f
N
となり
f
f
p
9
N
を得る。
これが電離層で反射する限界(臨界)周波数[Hz]と電子密度[個/m
3]との関係で
ある。
一例として、8MHzの電波を上空に送信して垂直入射させた電磁波が電離層にて
反射する為の電子密度は
]
m
個
[
10
901
7
9
10
8
9
3 11 2 6 2/
.
f
N
以上となる必要がある。
送信電磁波の周波数fと電子密度との関係
(電子密度のNを用いた臨界周波数の導出)
プラズマの中では n<1だがn=0だと 全反射 これを満足するためには
2 2 9 1 f N となる必要なので 2
2 9 N f をさらに変形して f 9 N の時全反射なので、f
N
p
9
であるから p f f なる関係、すなわち臨界周波数よりも送信する周波数の方が小さいと全反射 √の中が+になるためには最大電子密度N
maxと臨界周波数
N
max 今、上空に行くほど高くなる電子密度を持つ 電離層について、一番上部での電子密度を Nmaxとする。 ここで、宇宙に飛び出さず反射する限界の周 波数をあらためて臨界周波数と呼び max9 N
f
c
で表される。つまり、周波数を高くするにつれて 徐々に高い電離層で反射し、fc 以上では反射 せずに電離層を通過することになる。f
1f
2f
3=f
Cf
4f
1< f
2< f
3< f
4 上にいくほど電子密度が高いから電子 イオン(電子が抜けた分子) 中性の原子、分子
プラズマ状態
x
t
x
0
T
p : p
プラズマ角周波数
電離層
0 2 2
m
Ne
p
N: 粒子密度 e: 粒子の電荷量 m: 粒子の質量+
‐
電離層発生のメカニズム
N2 O2上空の分子
・・・(1) p p T 1
コンデンサ・モデル
誘電率静電容量
d
S
C
電離層を誘電体と仮定
プラズマ状態
での誘電率
2 2 01
p 電子密度高度
大 小 大 小d
C
s
誘電体
e
e:
誘電率
e
N
・・・(2) 0
0
e
0 0 真空の誘電率よりも小さい よって
2 2 2 2 2 2 01
2
2
1
1
f
f
f
f
p p p r
なので 2 21
f
f
n
p となる。 比誘電率は1より小さいであるが、J=0としてあらたに電磁波が電離層を伝搬中の変位電流Ji を c d i
J
J
J
t
D
ここに、Jd ・・・自由空間での電界の変位電流 Jc ・・・電子の運動による対流電流(プラス ゙マ中に発生) と定義する。そして電磁波の変位電流 Jdは伝搬中の電磁波の電界の振幅をE とし、時間因子も合わせて表示すれば t j dj
Ee
J
0 -m : 質量 e : 電荷量 v : 速度 と表現される。 電界の時間変化が 変位電流を生み出す t j
Ee
e
dt
dv
m
一方、JC =N・evを求めることを考えると、電界Eによっ て電子は加速されるが、その際に粒子として次の運 動方程式を満足する必要がある。 電界の印加による加速分 初速度を零とすれば、上記微分 方程式を変形してみると t jEe
m
e
dt
dv
であるから、vを求めるためこの両辺を tで積分すれば t j t je
j
E
m
e
dt
e
E
m
e
v
1
となるから、これを整理すれば運動方程式からの電離層の屈折率の導出
電離層内の屈折率は内部に発生する変位電流を整理すれば求まる。まず、アンペア・ マクスウェルの法則はt
H
J
D
,D
E
0
電界の時間変化により発生する磁界の中心には変位電流が流 れていると仮定 積分した結果t j
Ee
m
j
e
v
t j t j c d iEe
m
j
Ne
Ee
j
J
J
J
0 2 となるので、電離層中を電磁波が伝搬時の 変位電流Ji は結局 t jEe
m
Ne
j
0 2 2 01
等価的な比誘電率(1より小さい) (これの√をとれば屈折率になる) を得る。 よって電子が速度vで移動した時の電子の運動による対流電流Jcは電離層内の電子密度 Nを用いて t j cEe
m
j
Ne
Nev
J
2
となる。10 102 103 104 105 106 107 (N/cm3) 電子密度 50 70 90 130 C D E F1 F2 50 100 500 (Km) 高 度
夜
日中
H 地表 t0 t1m
t
t
H
1 03
10
82
p
電離層における高度分布と測定法
i1 r2 1
2
1 sini1 2 sinr2 1 2 1 2 sin sinr i i1 r2 1
2
i1 r2 1
2
1 2
1 2
電離層(プラズマ媒質中)におけるスネルの法則
~光の屈折~ ・・・(3) i1 r2 0
0
大気 電離層 一般的な媒質 電離層でではε r <1なのでi
0
i r sin sin 0 1 sin 0 i 1 sinr 90 i 90 r 0 sin i i
0
90 r 0 sin r 90 r 1 sini r 0 i sin i 0
射入射時の電離層における屈折および反射条件
電離層(地表)に平行に進む条件 i r sin sin 0 0
r電離層の反射
電離層の反射の条件 ・・・(5) 下部は誘電率が低いので 海面での反射は起こらない1 2 i1 r2 1
2
i2 3
r3 2 2 1 1 sini sinr 3 3 2 2 sini sinr 2 2 r i 3 3 2 2 sinr sinr 3 3 1 1 sini sinr 1 3 1 3 sin sinr i 2 3 i r
0
i r sin sin 0 電離層内における多重屈折
・・・(4) となり、多重層におけるスネルの法則は の条件では となるから、これを整理すれば となる。電離層(fp) Q P R i i 電離層を利用した長距離通信では、送受信間の距離と反射により可能となる周波数の把握が重要 である。図に示す様に電磁波が電離層により反射する場合、垂直入射(i=0)時には反射する臨界周波 数(その周波数以下なら反射する)はf=fpであるが、射入射する場合にはスネルの法則からfとfpには 以下の関係が成り立つ。
)
cos(
1
)
sec(
i
f
i
f
f
p
p f : 送信する電磁波の周波数(Hz) fp :電離層の臨界(プラズマ)周波数(Hz) N : 電離層の電子密度(個/S) i : 入射角電離層への斜入射時の反射条件
・・・(1) 以上のことから、(1)式を満足するように入射角iが大きい場 合には反射による通信が可能な周波数 fは高くなることが 分かる。0
i(rad)
i(deg)
0
) cos( 1 i1
4 45
1.414
2 90
∞
3 60
2.0
)
cos(
1
i
の性質 この関係を正割の法則と呼び、fなる周波数の電波が電離層 へ斜入射した時に反射するための最小角度 i を定める。つま り、その角度以上では電磁波は反射する。 入射角が大きいほど同じ周波数でも 反射しやすい(スネルより) f=fp×∞より、無限大の周波数も反射する 実際には f=fp (垂直入射)(1)式の導出
まず、点Pから発射した電磁波がQに対して角度i0 で入射した場合、入射角i0 と透過角i1とはスネルの 法則より次式が成り立つ)
sin(
)
sin(
0 1 1 0i
n
i
n
ここで、noおよびn1は自由空間および電離 層の屈折率として no=1 , 2 2 11
f
f
n
p で与えられる。 電離層(fp) Q P R i0 i0n
0n
1 i1 ・・・(2)ここでi0 =90°の時、i1 =90°すなわち
sin(
)
1
1
i
の時に電磁波は全反射するので 2 2 1 1 1 0)
sin(
)
1
sin(
f
f
n
i
n
i
p 1
01
22f
f
i
sin
p (2)式はと変形できる。そこで、両辺の2乗をとれば 2 2 2 0 2
1
1
)
(
sin
f
f
f
f
i
p p となり、これを整理すれば)
(
cos
)
(
sin
1
2 0 2 0 2i
i
f
f
p
となる。さらにf, fp >0である から両辺のルートはf
cos(
i
0)
f
p
となり、これをfについて整理すれば)
sec(
)
cos(
i
0f
i
0f
f
p
p
を得る。 池澤俊冶郎, “電波工学概説”, 森北出版, PP.117-120.0 sin i 2 2
1
sin
pi
)
1
(
2 2 0
p反射の条件
1
20
2
p p
0
7MHz 10.5MHz 14MHz 21MHzMHz
7
pf
i
i
p 2 2 2 2cos
sin
1
i
f
f
p pcos
1
・・・(5) ・・・(2)反射周波数と入射角との関係のまとめ
コセカント法則 ・・・(6)H p
f
i D: 距離 R H i cos R 2 2 ) 2 (D H R H
R
i
f
f
p p
cos
1
2 2)
2
(
1
D
H
H
p
2)
2
(
1
H
D
p
2)
2
(
1
H
D
f
f
m
p
MUF (Maximum Usable Freq. )
・・・(6)
・・・(7) 距離dにて通信が可能な最大周波数を
f
m とすればf=0.8fp f=2fp f=4fp 参考資料:電波受験
i
f
f
pcos
・・・(6)’i
(N/cm3) 電子密度 50 70 90 130 C D E F1 F2 50 100 500 (Km) 高 度
夜
日中
D層減衰
F
1,F
2 i1 r2 1
2
例題1
無線中継所において、周波数15MHzを用い11時に、また周波数8MHzを用いて24時 に通信をしたい。11時および24時における電離層の臨界プラズマ周波数をそれぞれ 13MHzおよび6MHzとすれば、その時刻に受信可能かを調べよ。 但し、送受信間の距離を600km、電離層の高さを400kmとする。また、大地および電 離層を平面とし、電離層の電子密度および送受信電界強度の変化は考慮しないものと する。 正割の法則より
i
cos
f
i
sec
f
f
p
p1
であるから、
0
.
87
15
13
cos
f
f
i
p を得る。ここで図における臨界→入射角 i での飛躍距離dを求めると、電離層の高さHは 400kmであるから
800
2
400
2
/
tan
d
d
H
d
i
となる。なお、ここでのdは通信できる最小距離であり、d=800・tan(i)として求められる。なお、 計算時にはtan(i)をcos(i)で表現する必要がある。 まず、11時にて与えられた反射する最小入射角度をiとすれば、15MHzではこれによりdは
448
87
.
0
)
87
.
0
(
1
800
tan
800
2
i
d
kmとなり、この値は *1 *1
2 21
1
cos
tan
より1
1
2 2
cos
tan
となるので、両辺のルートをとった 電波を反射できる最小角度iにおける通信可能な最短距離を意味する。1
2 2
cos
sin
の両辺に1/cos2θをかけると
2 2 21
1
cos
cos
sin
となり、これを整理すれば (a) 通信可 (b) 通信不可 H 電離層(fp ) 突き抜け現象 i i1 i0 (b) (a) dH d/2 i
H
/
D
i
tan
2
H d/2 i 伝搬長なる関係を用いれば
cos
cos
cos
cos
cos
tan
2 2 2 21
1
1
1
を得る。これより、11時ではd=448km<送受信間距離600kmであり、dよりも大きい為に入射 角i0をiより大きくすることができ電波が受信可能となる。 次に、24時では同様に
0
.
75
15
6
cos
f
f
i
p となるので 飛躍距離dは
704
75
.
0
)
75
.
0
(
1
800
)
cos(
)
cos(
1
tan
800
2
i
i
i
d
kmを得る。 以上のことから、24時ではd=704km>送受信間距離600kmより、dより距離が狭い場合に は入射角度入射角i1はiより小さくなり、電波は電離層で反射されず通信不可能となる。
cos
cos
cos
cos
tan
2 2 21
1
f[MHz]
送信周波数
f
p[MHz]
垂直入射時の臨
界周波数
通信
15
13
○
8
6
×
例題2
電離層(E層)の臨界周波数fp が4MHzの時、400kmの距離をE層の反射で伝搬する最高使用 周波数(MUF)を求めよ。但し、E層の高さを100kmとする。
h
r
i
cos
f
i
sin
f
f
p
p
610
4
1
MUF
E層 i d=400km h=100km r r T Q R 図において 223 2 2 2 h d r kmとなるので 6 3 3 610
92
.
8
10
100
10
223
10
4
Hz
8
.
92
MHz となる。 注: 最低使用周波数 (LUF)はD層の電波吸 収を考慮して決められる。
r
h
i
cos
より対流圏伝搬
地球の球面を考慮した多層分割モデルを右 図に示す。第0層と1層との間にスネルの法則 を適用すれば、入射角inと波数knとの関係に おいて、次式が得られる。)
sin(
)
sin(
0 1 1 0i
k
i
k
ここで、各層の屈折率と波数には以下の 関係がある。 0k
n
k
i
i
次に、各層は充分薄いものとし、地球を球形 とみなせば、地球の中心Oから地表までの距 離R0 と第1層までの距離R1 で作られる三角形 において、正弦の法則から次式が求まる。)
sin(
)
sin(
0 1 0 0i
R
i
R
O 地球の中心 0 1 2 R0 R1 R2 R3 k0, n0 k1, n1 k2, n2 対地球面 i0 i 0 ’ i1 i2 i1 ’ i2 ’ ・・・(1) ・・・(2) ・・・(3) この3角形 において(1)および(3)式を用いて、第0層から第1層への入射角i0’を消去することを考えると