EEM-RTM
理論説明書
株式会社 EEM
目次
1 概要...................................................................................3
2 レイローンチング法.....................................................................4
3 伝達公式と送受信アンテナ...............................................................6
3.1 伝達公式..........................................................................6
3.2 アンテナ指向性....................................................................6
3.3 アンテナ利得......................................................................8
3.4 送信アンテナの偏波................................................................8
4 Fresnel反射係数.......................................................................9
5 多層壁面の反射・透過係数..............................................................11
6 直接波の判定法........................................................................13
7 鏡影点の計算法........................................................................14
8 回折波の計算法........................................................................15
9 遅延特性..............................................................................19
参考文献.................................................................................21
2
1 概要
電磁波シミュレータEEM-RTMは波長に比べて十分に広い領域の電波伝搬特性を光学的手法(レイトレー
シング法)によって数値解析するものである。
解析物体は任意形状の面の集まりとし、波長に比べ十分大きいものとする。その物性値(誘電率、導電
率、厚さ)を個別に指定することができる。
送信アンテナと受信アンテナの指向性を個別に指定することができる。
壁の反射・透過特性には、光学的反射・透過モデルを使用する。
伝搬経路の計算方法にはレイローンチング法を用いる。
すなわち、送信点から全方向に多数のレイを放射し、各レイの多重反射の軌跡を追跡し、受信点の近
くを通過する経路を求める。
壁の透過を考慮する場合は、壁に当たったレイは反射波と透過波の二つに枝分かれする。
レイローンチング法は波の直進性を仮定しているので、回折波については別途、伝搬経路を探索する。
このとき、反射+回折する経路も考慮する。なお、回折特性には楔モデルの回折理論を用いる。
以上の方法で、直接波、反射波、回折波、透過波をすべて計算することができる。
3
2 レイローンチング法
レイローンチング法とは、図2.1のように送信点から全方向に多数のレイ(光線)を放射し、各レイが
面で反射される経路を追跡し、指定された受信点の近傍を通過したとき、受信されたものとみなす方
法である。
図2.1 レイローンチング法
反射点r0を通り単位方向ベクトルが u の直線と、受信点rとの距離dは次式で計算される。
d=∣u× r−r0∣ (2.1)
この値が受信円の半径より小さいとき受信したとみなす。ただし、隣接するレイの間隔は伝搬距離に
比例して大きくなるので、受信円の半径を伝搬距離に比例させる。
伝搬経路を漏れなく求めるに、受信円の半径を少し大きめにとり、後で重複分を削除する。
すなわち、伝搬経路の反射面の番号を保存しておき、異なる伝搬経路の反射面の番号がすべて一致す
ることはないことから重複分を削除する。
その後、受信円の広がりによる誤差をなくすために伝搬経路を厳密に評価する。
図2.2は 2回反射の例であり、送信点Aから出た波が面1と面2で2回反射して受信点Bに到達するも
のとする。点Aの面1に関する鏡影をA'、点Bの面2に関する鏡影をB'とする。
正しい伝搬経路が存在するためには以下の条件がすべて成り立たなければならない。
(1)線分A'B'が面2を横切る。その交点を反射点 Qとする。
(2)線分A'Qが面1を横切る。その交点を反射点 Pとする。
(3)線分APを直接波が伝搬する。
(4)線分 PQを直接波が伝搬する。
(5)線分 QBを直接波が伝搬する。
なお、直接波が伝搬する条件、および、鏡影点の計算方法については後述する。
4
面1
面2
受信点
送信点
受信されるレイ
図2.2 2回反射の概念図
放射レイの緯度方向の分割数を Nとすると、第一放射レイの総数Mは次式で評価できる。
M=∑i=1
N {2N sin iN }
≈2N 2∫0
1
sin x dx N≫1
= 4N 2
(2.2)
壁による透過波を考える場合は、図2.3のように壁に当たるごとに反射波と透過波の2つに分岐する。
図2.3 反射+透過
5
面2の鏡影B'
反射点 Q
送信点A 面2
反射点 P
面1
受信点B面1の鏡影A'
壁
壁
壁
送信点
3 伝達公式と送受信アンテナ
3.1 伝達公式
Friisの伝達公式を、複数の送信アンテナがあり、さらにそれぞれ複数の伝搬経路がある場合に拡張す
る。伝搬経路の位相差を考慮しない電力和は
P r=∑n [∑i G t ,n ,iGr , n ,i∣Rn , i∣
2{ 4 l n , i
exp − l n , i}2]P t ,n (3.1)
となり、伝搬経路の位相差を考慮した電力和は
P r=∣∑n [∑i E t , n , i E r ,n ,i Rn ,i 4 ln , i exp {− jk l n , i}]exp jnP t , n∣
2
(3.2)
となる。ここで、総和 nは送信アンテナについてとり、総和iはそれぞれの送信アンテナから出るすべ
ての伝搬経路についてとる。また、Pr:受信電力、Pt,n:送信アンテナnの送信電力、φn: 送信アンテ
ナnの位相、λ:波長、k=2π/λ:波数、α:空間の減衰定数、ln,i:伝搬経路長である。Rn,iは反射・
透過・回折係数の積であり複素数である。さらに、Gt,n,i Gr,n,i Et,n,i Er,n,iは送信アンテナnか
らの伝搬経路 iの送信方向および受信方向の電力絶対利得または電界絶対利得である。アンテナ指向
性の位相は一般に観測方向で異なるので、Et,n,iとEr,n,iは複素数になる。電力絶対利得は電界絶対
利得の2乗である。
G t ,n ,i=∣E t ,n ,i∣2
G r ,n ,i=∣E r , n ,i∣2 (3.3)
また、無指向性アンテナの受信電力Prと受信電界 Eには下記の関係がある。[6]
E2
120=4 P r
2 (3.4)
3.2 アンテナ指向性
送信アンテナまたは受信アンテナに以下のような指向性を設定することができる。
以下では送信アンテナと受信アンテナの扱いは同じなので、電界絶対利得Et,Erを単にEと表す。
観測方向を , とすると観測方向の単位ベクトルは次式になる。
r=x sincosy sinsinz cos (3.5)
6
図3.1 極座標
(1)無指向性
E ,=1 (3.6)
(2)ダイポール
Z方向を向いたダイポールアンテナであり、3dBビーム幅を wとする。
E (θ ,ϕ )={K cos{π2
θ−π/2w } (|θ−π/2|≤w )
0 (|θ−π/2|>w )(3.7)
ここで、Kは電界利得であり、計算方法は後述する。
(3)ビーム
ビームの中心方向を 0,0 、θ方向とφ方向の3dBビーム幅を w ,w とする。
E ,={K cos{2
−0
w }cos{2
−0
w } ∣−0∣≤w , ∣−0∣≤w0 others
(3.8)
ここで、Kは電界利得であり、計算方法は後述する。
式(3.8)はパラメータによっては非物理的な特性になるので、指向性を図で確認することが必要である。
(4)ファイル指定
全方向の指向性(θ成分とφ成分の振幅と位相)をファイルで指定する。さらに Z軸の周りに回転する
ことができる。
7
Z
X
Yφ
θr
r
(1)無指向性 (2)ダイポール (3)ビーム
図3.2 アンテナ指向性
3.3 アンテナ利得
アンテナの最大放射方向(ボアサイト)を max ,max とし、
Emax≡E max ,max (3.9)
とおく。このとき電界利得は次式で計算される。
K={ 4 Emax2
∫0
2
d ∫0
sin d E2 ,}12
(3.10)
この Kを用いると3.2(2)(3)は絶対利得になる。3.2(4)では最大値で正規化した後 Kをかけて絶対利得
にする。
3.4 送信アンテナの偏波
送信アンテナの偏波方向が垂直偏波・水平偏波・円偏波のとき直接波の電界は次式のようになる。
E ,={ E , 垂直偏波E , 水平偏波
E ,± j
2円偏波
(3.11)
3.2(1)(2)(3)では偏波を上のいずれかに指定する。
ただし、ダイポールアンテナなどを除いては、一般的に全方向で同じ偏波特性を持つことはないので、
式(3.11)は近似的なものである。
より正確に偏波特性を表すには3.2(4)を用いて全方向指向性をファイルで指定することが必要である。
なお、式(3.1)(3.2)の G r , n , i , E r , n , i は受信アンテナの偏波不整合損を含んでいる。
8
4 Fresnel反射係数
均質な半空間の反射係数はFresnelの公式によって計算される。[1]
図4.1のように座標系をとる。入射波を E i とすると、反射波 E r は次式で計算される。
E r=V Ei⋅v i vrH Ei⋅hi hr (4.1)
V=−sin2− cos
−sin2 cosTM波、 P波 (4.2)
H=cos−−sin2
cos−sin2TE波、 S波 (4.3)
=r−j
0
=2 f , 0=8.854×10−12 [F /m] (4.4)
cos=∣n⋅r i∣ (4.5)
h i=hr=r i×n
∣r i×n∣(4.6)
v i=hi×ri (4.7)
r r=ri−2 n⋅ri n (4.8)
vr=rr×hr (4.9)
図 4.1 平面による反射
9
θθ
反射波入射波
反射面 εr,σ
n
hi
vi
ri hrvr
rr
式(4.2)(4.3)は特殊な場合として、以下の性質をもっている
V =0= H =0=−−1
10
(4.10)
V = /2=1, H = /2=−1(4.11)
V ∞=H ∞=−1 PEC (4.12)
図 4.2に ϵr=2,10,80 , σ=0 のときの V ,H の入射角特性を示す。
σ=0のとき、 H は常に負であり、 V は B のとき負、 B のとき正である。ここで、
B はBrewster角と呼ばれ、次式で計算される。
B=arctan r (4.13)
図 4.2 反射係数の入射角特性(σ=0)
10
0 10 20 30 40 50 60 70 80 90
-1.0
-0.8
-0.6
-0.4
-0.2
0.0
0.2
0.4
0.6
0.8
1.0
θ[度 ]
ΓV
ΓH
εr=10
εr=2
εr=80
5 多層壁面の反射・透過係数
図5.1に示す多層壁面の反射・透過係数は特性行列法(F行列法)によって計算される。[2][3][4]
層の数を Nとし、各層の複素誘電率、複素透磁率、厚さを n , n , hn とする。両端の領域は真空とす
る。
各境界での入射角は Snellの法則より次式で計算される。
√ϵnμnsinθn=√ϵ0μ0sinθ (n=1,...N) (5.1)
第 n層の特性行列は次式で与えられる。
M n=[ cosδn jαn sinδnjαn
−1sinδn cosδn ] (n=1,...N) (5.2)
αn={√μnϵn cosθn (TM波 , P波)
√μnϵn
1cosθn
(TE波 ,S波) (n=0,...,N) (5.3)
n=n nhncos n n=1,⋯, N (5.4)
これより全体の特性行列が次式で得られる。
M=M 1M 2⋯M N=[m11 m12
m21 m22] (5.5)
[ E i
H i]=M [ E t
H t ] (5.6)
以上から、多層壁全体の電界反射・透過係数(r,t)は次式で計算される。
r=(m11α0+m12)−α0(m21α0+m22)(m11α0+m12)+α0(m21α0+m22) (5.7)
t=2α0
(m11α0+m12)+α0(m21α0+m22) (5.8)
11
図5.1 多層壁面モデル
12
θ
z0
θ2
θN
z1
z2
zN-1
zN
θ1
θ
z
・・・
Ei,Hiε0,μ0
Et,Htε0,μ0
h1M1ε1μ1
h2M2ε2μ2
hNMNεN
μN
6 直接波の判定法
2点間を直接波が伝搬するか否かは、2点を結ぶ線分上に物体が存在するか否かで判定する。
平面は3角形の集合で表現されるので、ここでは3角形を考える。送信点と受信点の位置を
r a= xa , ya , z a , rb= xb , yb , zb とし、3角形の3頂点の位置を r i=x i , y i , z i i=1,2,3 とすると、受信点に直接波が届かない条件は以下のようになる。
伝搬経路と3角形との交点を2通りで表現すると、
r a1rb−r a=r12r 2−r13r3−r1 (6.1)
となり、これを成分に分けると、
[ xa−xb x2−x1 x3−x1y a− yb y2− y1 y3− y1z a−zb z 2−z1 z3−z1
][1
2
3]=[ xa−x1
ya− y1za−z1
] (6.2)
となる。これを解いて、
011, 02, 03, 231 (6.3)
のとき受信点に直接波が届かない。
図6.1 伝搬経路と3角形
13
送信点 ra
r3
r2
r1
受信点 rb
7 鏡影点の計算法
4角形による鏡影点は以下のようにして計算する。
図7.1のように4角形の4頂点の位置ベクトルをr1~r4とし、点 Pの鏡影点を P'、線分 PP'と4角形と
の交点を Qとする。4角形の法線ベクトルを n とすると、点 Qと点 P'は
rQ=r P n(7.1)
r P '=r P2 n(7.2)
n= r3−r1× r4−r2∣ r3−r1× r4−r2∣ (7.3)
とあらわされ、条件
rQ−r1⋅n=0 (7.4)
から未知数αが求まり、式(7.2)より鏡影点は次式になる。
r P '=r p−2 {r P−r1⋅n } n (7.5)
式(7.5)は n の符号を変えても変わらない。すなわち面の表裏を知る必要がない。また、頂点1は2,
3,4に変えても同じである。
図7.1 4角形による鏡影点
14
P
3
4Q
2
1
P'
n
8 回折波の計算法
直接波が到達しない点には回折波が到達する。回折波の大きさは楔モデルの回折理論[5][6]で計算す
る。
図8.1のように、送信点と受信点の間に物体が存在し、直接波がさえぎられ、回折波が到達するもの
とする。物体を楔モデルとみなすと、回折波の大きさは、楔の稜線から送信点と受信点と結ぶ直線ま
での最短距離で決まる。
図8.1 回折波
上で述べた距離は以下のようにして計算される。
送信点および受信点の位置ベクトルr1,r2、楔の稜線の両端の位置ベクトルをr3,r4とする。
送信点と受信点を結ぶ線分上の点を P、楔の稜線上の点を Qとすると、
r P=r1u r2−r1 (8.1)
rQ=r3v r 4−r3 (8.2)
となり、PQ間の距離は
PQ 2=∣{r1u r2−r1}−{r 3v r 4−r 3}∣2
=au2buvcv2duev f(8.3)
a=∣r 2−r1∣2
b=−2 r2−r 1⋅r 4−r 3c=∣r 4−r3∣
2
d=−2 r2−r 1⋅r3−r1e=−2 r4−r3⋅ r1−r3f=∣r 3−r 1∣
2
(8.4)
となる。u,vに関する2次形式(8.3)を変形すると、
PQ2=aubvd2a
2
4 ac−b2
4a v 2ae−bd4a c−b2
2
bde−cd 2−a e2
4a c−b2 f
(8.5)
となる。式(8.4)より、
15
回折波
受信点 r2
送信点 r1
r3
r4
Q
P
l1
l2
a0, 4a c−b20 (8.6)
であるから、式(8.5)の平方根は
u=be−2c d
4a c−b2 (8.7)
v=bd−2a e
4a c−b2 (8.8)
のとき、最小値
h= bde−cd 2−ae2
4a c−b2 f 12
(8.9)
をとる。
なお、上で求めた点 Qを1回回折波が通るには、送信点と点 Q、および点 Qと受信点の間に物体が存在
してはならない。
以上の判定をすべての物体の辺について行う。
回折波の大きさは以下の式で計算される。[4][5]
Ed=E0F ± 2 1l 1 1
l 2 h (8.10)
F z ≡ 1
2∣I z ∣
(8.11)
I z ≡∫z
∞
e− j
2u2
du −∞ z∞(8.12)
ここで、式(8.10)の E0は物体が存在しないときの直接波の大きさ、 l 1 は送信点と楔の距離、 l 2は楔と受信点の距離であり、符号は送受信間が見通し外のとき正、見通し内のとき負である。
式(8.11)は z=−∞ ,0,∞ のとき F z =1,0.5,0 になる。
式(8.12)の計算には以下の近似式[7]を用いる。
I z =e− j
2z 2
{g z − j f z } 0z∞ (8.13)
I −z =1− j −I z (8.14)
f z = 10.926z
21.792z3.104z2 z 0z∞ , ∣ z ∣2×10−3 (8.15)
g z= 1
24.142z3.492z26.670z3 z 0z∞ , ∣ z ∣2×10−3 (8.16)
2回回折波については、図8.2のように同じ物体の2つの稜線(上部または側面)を経由するものだけを
考え、図8.3のように異なる物体の稜線を経由するものを考えない。後者は計算時間が物体の数の2乗
に比例するためである。
16
図8.2 計算される2回回折波
図8.3 計算されない2回回折波
2回回折波の大きさについては式(8.10)に対応する式が以下のようになる。
Ed=E0F 2 1l 1 1
l2 h1F 2 1l 2 1
l3 h2 (8.17)
図 8.4 2回回折波断面図
17
送信点受信点
送信点 受信点
物体
物体
(1)上部経由 (2)側面経由
物体1
物体2送信点
受信点
送信点受信点
回折点回折点
l1
l2
l3h1
h2
回折と反射を含む場合は、図8.5のように[反射+回折]と[回折+反射]を計算する。
この場合、計算時間は物体の数の2乗に比例する。
図8.5 回折+反射の計算モデル
18
(1)反射 +回折 (2)回折 +反射
9 遅延特性
ある受信点で横軸を時間、縦軸を電力としてプロットしたものを遅延プロファイルと呼ぶ。
光学的反射モデルでは遅延プロファイルは図9.1のように離散的なものになる。
図9.1 遅延プロファイル
受信電力の時間変化を p(t)とすると、電力和P、平均遅延TD、遅延スプレッドSが次式で定義される。
[6][8]
P=∫ pt dt=∑i
pt i (9.1)
T D=1P∫ (t−t 0) p(t )dt= 1
P∑i
(t i−t 0) p(t i) (9.2)
S=[ 1P∫ {t−t 0−T D}2 pt dt ]12
=[ 1P∫ t−t 02 p t dt−T D
2 ]12
=[ 1P∑i
ti−t 02 p t i−T D
2 ]12
(9.3)
ここで、t0は第 1波(見通し内では直接波)が到来する時刻である。
仮に遅延プロファイルが連続量であり次式のような指数関数であると仮定すると
p (t )=1τ e−t−t0
τ (t>t 0) (9.4)
公式
∫0
∞
xne−xdx=n! (n≥0) (9.5)
から
19
p(t)
tt0
t1
t2
t3
t4
P=∫0
∞
p(t )dt=1
∫0
∞
(t−t 0) p(t )dt=τ
∫0
∞
(t−t0)2 p(t )dt=2 τ2
となり、
T D=S=τとなる。すなわち、遅延プロファイルが指数関数のときは平均遅延と遅延スプレッドは一致し、受信
電力が1/eになる時間である。
受信電圧和(式(3.2)の||内、複素数)V(x)を受信点を移動して求めたものから次式により空間相関係数
が計算される。
d =∣∫V x V * xd dx∣
∫∣V x∣2dx(*:複素共役) (9.6)
V x=∫V x , t dt=∑i
V x , t i (9.7)
20
参考文献
[1]J.D.Kraus, "Electromagnetics, 4th edition", McGraw-Hill, 1992.
[2]M.ボルン、E.ウォルフ「光学の原理 I」東海大学出版、1985
[3]飯島泰蔵監修「電磁界の近代解析法」電子通信学会、1979
[4]安達三郎、米山務「電波伝送工学」コロナ社、1981
[5]W.L.Stutzman and G.A.Thiele, "Antenna Theory and Design, 2nd edition", John Wiley &
Sons, 1998.
[6]進士昌明編著「無線通信の電波伝搬」電子情報通信学会、1992
[7]M.Abramowitz and I.Stegun ed., "Handbook of Mathematical Functions", Dover, 1965.
[8]S.R.Saunders, "Antennas and propagation for wireless communication systems", John Wiley
& Sons, 1999.
21