電気通信大学

アンテナ・伝播研究会(H27.12.10)

電波伝搬：物理現象とモデリング～先達の足跡を辿る～
唐沢 好男
電気通信大学1

電気通信大学

発表の内容
 電波伝搬研究概観（５０年（以上）を振り返って） 先達の足跡を辿る（源流から潮流へ）・フェージング確率分布：
仲上 稔氏
・移動伝搬モデル：
奥村 善久氏・移動通信の三位一体研究思想： 池上 文夫氏
・降雨減衰モデル：
森田 和夫氏 まとめ

2

電気通信大学

電波伝搬の研究： 理論・実測・モデリング

半世紀前は、「伝搬はシステムを制す」の時代
3電気通信大学

我が国の電波伝搬研究 (URSI-F発表件数)推移
伝見
搬通
外URSI-F 300回 (1987.12: 小口知宏)（唐沢感覚）

基礎
伝見
搬通
内通衛
信星
屋移
内動通通
信信
セリンモシーント
グ 1955

50年前

6065

30年前

7075

8085MIMOOFDMなどシステム指向2000

154電気通信大学

我が国の先達がその構築に主導的役割を担い、かつ、今日の無線通信の基礎となっている伝搬研究とその結晶である伝搬モデルを取り上げる。（不易流行）
移動伝搬理論を支える確率分布：仲上ｎ分布（今日の仲上・ライス分布）と仲上ｍ分布を生み出した仲上稔氏の研究（1940年代）
移動通信伝搬の概念構築と回線設計手法を打ち立てた
奥村善久氏の研究開発（1960年代）
広帯域移動通信時代を予測してアンテナ(A)・伝搬(P)・システム(S)の融合研究（A・P・S三位一体研究）を提唱した池上文夫氏の思想（1980年代）
無線通信の重要課題、降雨減衰のモデル化に取り組んだ
森田和夫氏の研究開発（1970年代）
5

電気通信大学

フェージングの確率分布の研究
仲上稔氏（国際電気通信株式会社）の研究（１９４０年代）
短波のフェージングの研究が移動通信のマルチパスフェージングの基礎的確率分布として花開いている
仲上が生み出した３つの確率分布
ｎ分布 （仲上･ライス分布）
ｍ分布ｑ分布 （仲上・ホイト分布）6電気通信大学

仲上ｎ分布
歴史【短波のフェージングの研究】
仲上稔, “位相ならびに振幅が不規則に変化する多数の振動の合成振幅に関する研究（其の二）,”電気通信学会誌, no. 202,1940.01.
【信号と雑音の研究】
S. O. Rice, “Mathematical analysis of random noise,” Bell
System Tech. Jour. (BSTJ), vol. 23, no. 3, 1994, and vol. 24, no.
1, pp. 46–156, 1945.
ライス分布 → ライス・仲上分布 → 仲上･ライス分布7電気通信大学

仲上ｎ分布
仲上・ライス分布の現代の理解xrf r (r ,  )  f x ( x) f y ( y )
y
r

x
y


f x  N (a,  2 )ay a 2  2ar cos   r 2 exp  
2222
rf r (r )  20

f r (r ,  )d0axΔr

 a  r   ar  Iexp
22  02 2  r2

r

2f y  N (0,  2 )f ( x, y)  f x ( x) f y ( y)この導出はライス流

8電気通信大学

仲上ｎ分布

誕生

仲上は振幅が異なる多数の正弦波の合成振幅rの確率分布を考えた
ハンケル変換形の特性関数（フーリエ・ベッセル変換）を利用


F (k )   J 0 (rk ) f (r )dr0

f (r )  r  k J 0 (rk ) F (k )dk0

振幅r1, r2, 位相差  の正弦波の合成波の包絡線（振幅）はr  r12  r22  2r1r2 cos 
ハンケル変換形の特性関数

F (k )  J 0 (kr1 ) J 0 (kr2 )9電気通信大学

仲上ｎ分布
振幅が異なる多数の正弦波の合成波のハンケル変換形特性関数

  2k 2  ( N  1)
F (k )   J 0 (kri )  exp  2 i 1

NN

22
    ri  / 2
 i1 振幅aの強い正弦波が加わった信号強度のハンケル変換型特性関数の逆変換（確率密度関数）は

  2k 2  dkf (r )  r  kJ0 (rk )J 0 (ak ) exp  02  a 2  r 2   ar r
I 2 exp  2  0
2  2    

同時代の洋の東西で、異なる目的、異なる手法で、後の時代の移動伝搬の基本式：仲上・ライス分布が生まれた

10電気通信大学

仲上ｍ分布

そのルーツ

短波のイベント毎の相対頻度特性

相対頻度特性：特殊方眼紙

0 .0 1

レイリー分布（ｍ＝１）

x = 1 + 2 – exp (2c / M)
M(M  20 log e)

0 .1

y = ln p' (for c  0)
= – ln p' (for c > 0)p' 1.0

レイリー分布0 .1

yx
0 .0 1-4 0-3 0-2 0-1 00

10c ( dB )ｍ分布の誕生の直線になる

どのイベントでもy = m x の直線になる
11電気通信大学

仲上m分布

現代の理解

m2m
 m 2
2 m 1m分布
f r   mr
exp   r （振幅の分布）Ω Γ m Ω r zm  0.5z  r2
ガンマ分布
（電力、SNRの分布）mm
 m m 1g z   mz exp   z Ω Γ m Ω m分布は、・数学モデルとしてガンマ分布から変換されたものではなく「発見的手法」によって生まれた・仲上･ライス分布を一定の範囲で近似する・他の関数と組み合わせた変換をするときの解析性に優れている
12

電気通信大学

確率分布の相互関係

IEEE Xplore “Nakagami distribution”
1430件のヒット → 世界的なモデル

13移動伝搬・移動通信研究のパイオニア
奥村善久氏（当時：電電公社通研）

電気通信大学
移動伝搬の原点（1962 ~ 1965）：
奥村カーブ

奥村カーブ
東京タワー設置の送信アンテナ

移動伝搬実験測定コース

写真類出典：
奥村幸彦：
通ソマガジン No. 29,
2014より受信電界強度測定車

14電気通信大学

移動伝搬基本構造モデル【奥村モデル】
瞬時値
長区間中央値

積算的確率過程（対数正規分布系）

短区間中央値

加算的確率過程（正規分布系）

複雑な移動伝搬環
境を物理現象の持つ性質を採り入れた形で組み立てており、見事

瞬時変動（レイリーフェージング）

短区間
中央値
中央値
（対数正規シャドーイング）
長区間
中央値

距離変動
（奥村カーブ）距離 (km)15

電気通信大学

移動伝搬研究（奥村カーブと奥村モデル）を含めた『世界初のセルラー電話ネットワーク、システム及び標準規格に対する先駆的貢献』により、2013年、全米工学アカデミーより、工学のノーベル賞と言われる「チャールズ・スターク・ドレイパー賞 （Charles Stark DraperPrize）」受賞。日本人研究者として初めての快挙。

写真類出典：
奥村幸彦：
通ソマガジン No. 29,
2014より16

電気通信大学 ）移動通信の研究哲学： A・P・S 三位一体の研究（1980～
池上文夫氏（当時：京都大学教授）

システムに様々な影響を与える電波伝搬要因

17電気通信大学

移動通信の研究哲学： A・P・S 三位一体の研究

PPASAPS

A

三位一体

送受信装置
(a) 第Ⅰ世代
伝搬特性をあるがままに受け入れ、装置は伝搬を、伝搬は装置を与えられたものとする受動的立場

S(b) 第Ⅱ世代
伝搬とアンテナは一つの伝送路を構成し、両者は独立なものではなく、互いに関係して良好な伝送路の実現を目指すとする能動的立場

A：アンテナ
P：伝搬
S：システム(c) 第Ⅲ世代
伝送特性の劣化の原因は伝送路のみにあるのではなく、それに反応して動くシステムにもあり、A・P・S三位一体研究が必須
18電気通信大学

ディジタル移動通信 A・P・S 三位一体研究の取り組み

19電気通信大学

降雨減衰推定法： 森田和夫氏（当時：電電公社通研）の研究（1960年代後半～1970年代）

・周波数が10GH程度以上で減衰が無視できなくなり高くなるほど致命的問題に・地上回線・衛星回線に共通の問題
・日本は、中緯度地域の中では降雨が多く、この問題が特に深刻に
20

電気通信大学

降雨減衰推定の一般的アプローチ（地上・衛星回線共通）

R( p)等価通路長 Le [km]A( p)  k R( p) Le ( p)n

時間率p%での降雨減衰[dB]

降雨強度の確率分布や等価通路長を定式化
降雨強度R[mm/h]での減衰係数[dB/km] するのが研究課題
21電気通信大学

降雨減衰推定：森田・樋口モデル（1975）［樋口伊佐夫氏（統計数理研究所）］
特徴： 物理現象を数式化した理論モデル

降雨強度分布： ガンマ分布（パラメータ値n, b が地域依存）（1940～1950、1971～1973年代データ利用）
降雨の空間相関特性：  (d )  exp(  d ) （パラメータ値:  ）区間雨量（通路上の各点の雨量のn乗の積分値）＝電波の減衰量がガンマ分布近似で表されることを理論解析で示した
22

電気通信大学

降雨減衰推定：森田・樋口モデル（1975）b enf ( A) An 1 exp  b e A (0  A[dB])Γ (n e )een e  n D 2 / M , be  b D /( M )4 3
1 6 M  2  D  2D   D   exp    D  2     n, b : 降雨強度のガンマ分布パラメータ, : 降雨減衰係数
: 降雨強度空間相関パラメータ、 D: 区間距離 (km)
筆者は複雑な自然現象を理論モデルで組み上げた森田の着眼点に敬意を表している。
実際の計算では、空間相関の形やパラメータ値αの地域依存性のため無線通信の国際標準化機関であるCCIR（現ITU-R）での世界標準には至らなかったが、長く我が国の無線回線（特にアクセス回線）降雨減衰推定法の計算法として実用に供されるとともに、その後の降雨減衰推定法高度化の礎となってきた。
23

電気通信大学

降雨減衰推定法の進展（高精度化）
森田・樋口モデル（1975）： 降雨強度・降雨減衰： ガンマ分布近似

細矢等改良モデル（1984）： 森田・樋口モデルのパラメータの最適化
20GHz帯回線への適用

細矢の降雨強度分布モデル（1988）： 降雨強度をM分布で近似
AMeDAS時間降水量データからの降雨強度推定

小野・唐沢モデル（2008）： 降雨強度・降雨減衰： M分布近似
年変化特性を考慮した安全係数概念の導入（我が国のマイクロ波帯の固定無線システムの回線設計基準に採用）

24電気通信大学

むすびに代えて電波伝搬の研究は、対流圏・電離圏、固定通信・衛星通信・移動通信の多岐にわたり、また、通信・放送・リモートセンシング・電磁干渉など異なる視点からの取り組みも多く、限られた紙面や時間で、その歴史を網羅するのは講演者の力が及ばない
そこで、現在の潮流の水源地を開拓した四人のパイオニア（1940, 1960, 1970, 1980年代）の研究に焦点を当てた。
ここで取り上げることができなかった新規性・有効性の高い電波伝搬研究も数多くある。

25