JP2880151B1

JP2880151B1 - データ群の特性境界識別方法、電波伝搬推定方法および装置

Info

Publication number: JP2880151B1
Application number: JP10046410A
Authority: JP
Inventors: 裕也増井; 信也関澤
Original assignee: 株式会社ワイ・アール・ピー移動通信基盤技術研究所
Priority date: 1998-02-13
Filing date: 1998-02-13
Publication date: 1999-04-05
Anticipated expiration: 2018-02-13
Also published as: JPH11234219A

Abstract

【要約】【課題】解析者の主観に依存することなく特性の変化
する境界を推定する特性境界識別方法を提供する。【解決手段】横軸は送信局から受信局までの距離ｄ、
縦軸は伝搬損失値ｐである。ブレークポイント位置ｄ_i
を決め、送信局に対してブレークポイント以前となる領
域Ｒ_jと、ブレークポイント以遠となる領域Ｒ_j+1とに
分ける。各領域ごとに関数近似を行い、平均二乗誤差ｅ
_j（ｄ_i）を計算する処理を、ブレークポイント位置ｄ
_iを変化させて行う。評価関数Ｅ（ｄ_i）を、各ブレー
クポイント位置ｄ_iごとに計算する。評価関数Ｅ
（ｄ_i）は、平均二乗誤差ｅ_j（ｄ_i）に重み付け値ｗ
_jを付けた加算であり、この重み付け値ｗ_jにより正規
化される。評価関数Ｅ（ｄ_i）に基づいて適切なブレー
クポイント位置を識別判定する。

Description

【発明の詳細な説明】

【０００１】

【発明の属する技術分野】本発明は、特性の変化する境
界を有するデータ群からその特性の変化する境界を識別
する方法に関するものである。特に、移動無線通信シス
テムの回線設計において必要となる電波伝搬特性に関し
て実測デ−タに基づいた推定に用いる電波伝搬推定方法
および装置に関するものである。

【０００２】

【従来の技術】移動無線通信システム等におけるセル置
局設計等の回線設計においては、実測デ−タに基づいて
電波伝搬特性を推定することが必要である。小田恭弘、
田中哲「マイクロ波帯ストリ−トセル伝搬損失特性」、
電子情報通信学会技術報告Ａ・Ｐ９４−１０５（１９９
５−１）、ｐｐ．９−１４に述べられているように、一
般に地上波の見通し電波伝搬損失特性は、送信局から近
い領域では自由空間伝搬損に従い２乗則で減衰し、所定
の距離（ブレークポイント）よりも遠い領域ではおおよ
そ４乗則で減衰することが知られている。したがって、
このブレークポイントが特性の変化する境界となり、こ
のブレークポイントを識別判定することにより有効セル
範囲や他局干渉波影響の推定が行えるのでセル置局設計
等の無線回線設計に有効である。

【０００３】このブレークポイントを判定するには、実
測データに基づき、送信局から近い領域では２乗則減衰
に従う減衰係数で１本目の回帰線を引き、送信局から遠
い領域では２本目の回帰線として４乗則の近似線を最も
フィットする位置に引いて、その交点になった位置をブ
レークポイントとしていた。しかし、２本の回帰線を引
くときに、その引き方が解析者の主観に依存してしまう
ため、高精度にブレークポイントの位置を推定すること
が困難であった。

【０００４】

【発明が解決しようとする課題】本発明は、上述した問
題点を解決するためになされたもので、数値解析手法に
より評価関数を計算して特性境界を判定することによ
り、データの解析者の主観に依存することなく、特性の
変化する境界を推定することができるデータ群の特性境
界識別方法を提供することを目的とするものである。ま
た、この方法および装置を電波伝搬特性の推定に適用し
た電波伝搬推定方法および装置を提供することを目的と
するものである。

【０００５】

【課題を解決するための手段】請求項１に記載の発明に
おいては、特性の変化する境界を有するデータの群か
ら、前記境界を識別するための方法であって、（ａ）前
記特性の変化する境界を仮定するステップ、（ｂ）仮定
された前記境界により分割された各領域の前記データで
あって、前記境界から数えてほぼ同数のデータとなる範
囲のデータに対して、それぞれ関数近似処理を施して適
切な近似線を決定し、そのときの近似誤差を求めるステ
ップ、（ｃ）仮定された前記境界を変化させて、前記ス
テップ（ｂ）を順次繰り返すステップ、（ｄ）仮定され
た前記境界の各々に対して、前記各領域ごとの前記近似
誤差について重み付け総和を計算して評価値とするステ
ップ、（ｅ）仮定された前記境界の各々に対応する前記
評価値の大きさに基づいて前記特性の変化する境界の値
を識別するステップを有するものである。

【０００６】また、請求項２に記載の発明においては、
特性の変化する境界を有するデータの群から、前記境界
を識別するための方法であって、（ａ）前記特性の変化
する境界を仮定するステップ、（ｂ）仮定された前記境
界により分割された各領域の前記データであって、前記
境界から数えてほぼ同数のデータとなる範囲のデータに
対して、それぞれ関数近似処理を施して適切な近似線を
決定し、そのときの近似誤差を求めるステップ、（ｃ）
前記各領域ごとの前記近似誤差について重み付け総和を
計算して評価値とするステップ、（ｄ）仮定された前記
境界を変化させて、前記ステップ（ｂ）〜（ｃ）を順次
繰り返すステップ、（ｅ）仮定された前記境界の各々に
対応する前記評価値の大きさに基づいて前記特性の変化
する境界の値を識別するステップを有するものである。

【０００７】したがって、請求項１または請求項２に記
載の発明によれば、データの解析者の主観に依存するこ
となく、特性の変化する境界を推定することができる。
関数近似処理に用いるデータとして、特性の変化する境
界から数えてほぼ同数のデータとなる範囲のデータとし
たことから、隣接する領域での関数近似処理精度にばら
つきが生じる可能性が少なくなるとともに、特性の変化
する境界を推定する上で影響が少ないデータから関数近
似に使用しないようにすることから関数近似処理の信頼
性が増すことになる。

【０００８】請求項３に記載の発明においては、特性の
変化する境界を有するデータの群から、前記境界を識別
するための方法であって、（ａ）前記特性の変化する境
界を仮定するステップ、（ｂ）仮定された前記境界によ
り分割された各領域の前記データに対して、それぞれ関
数近似処理を施して適切な近似線を決定し、そのときの
近似誤差を求めるステップ、（ｃ）仮定された前記境界
を変化させて、前記ステップ（ｂ）を順次繰り返すステ
ップ、（ｄ）前記ステップ（ｂ）において求めた各領域
ごとの前記近似誤差について、仮定された前記境界を変
化させた範囲全体にわたって前記近似誤差の総和を求
め、前記総和の逆数を各領域ごとに与えられる重み値と
して設定するステップ、（ｅ）仮定された前記境界の各
々に対して、前記各領域ごとの前記近似誤差について前
記重み値により重み付けられた総和を計算して評価値と
するステップ、（ｆ）仮定された前記境界の各々に対応
する前記評価値の大きさに基づいて前記特性の変化する
境界の値を識別するステップを有するものである。した
がって、データの解析者の主観に依存することなく、特
性の変化する境界を推定することができる。評価関数を
正規化するための重み値を演算で求めることができるよ
うにしたことから、重み値を試行錯誤で決定する必要が
なくなる。

【０００９】請求項４に記載の発明においては、送受信
局間の伝搬距離に応じて受信電界強度をプロットした受
信電界強度データの群に基づいて前記受信電界強度の減
衰係数が急激に変化するブレークポイントの位置を推定
する電波伝搬推定方法であって、（ａ）前記ブレークポ
イントの位置を仮定するステップ、（ｂ）仮定された前
記ブレークポイントの位置により分割された各領域の前
記受信電界強度データに対して、それぞれ関数近似処理
を施して適切な近似線を決定し、そのときの近似誤差を
求めるステップ、（ｃ）仮定された前記ブレークポイン
トの位置を変化させて、前記ステップ（ｂ）を順次繰り
返すステップ、（ｄ）仮定された前記ブレークポイント
の位置の各々に対して、前記各領域ごとの前記近似誤差
について重み付け総和を計算して評価値とするステッ
プ、（ｅ）仮定された前記ブレークポイントの位置の各
々に対応する前記評価値の大きさに基づいて前記ブレー
クポイントの位置を識別するステップを有するものであ
る。

【００１０】また、請求項５に記載の発明においては、
送受信局間の伝搬距離に応じて受信電界強度をプロット
した受信電界強度データの群に基づいて前記受信電界強
度の減衰係数が急激に変化するブレークポイントの位置
を推定する電波伝搬推定方法であって、（ａ）前記ブレ
ークポイントの位置を仮定するステップ、（ｂ）仮定さ
れた前記ブレークポイントの位置により分割された各領
域の前記受信電界強度データに対して、それぞれ関数近
似処理を施して適切な近似線を決定し、そのときの近似
誤差を求めるステップ、（ｃ）前記各領域ごとの前記近
似誤差について重み付け総和を計算して評価値とするス
テップ、（ｄ）仮定された前記ブレークポイントの位置
を変化させて、前記ステップ（ｂ）〜（ｃ）を順次繰り
返すステップ、（ｅ）仮定された前記ブレークポイント
の位置の各々に対応する前記評価値の大きさに基づいて
前記ブレークポイントの位置を識別するステップを有す
るものである。

【００１１】したがって、請求項４または請求項５に記
載の発明によれば、データの解析者の主観に依存するこ
となく、ブレークポイントの位置を推定することができ
る。

【００１２】請求項６に記載の発明においては、送受信
局間の伝搬距離に応じて受信電界強度をプロットした受
信電界強度データの群に基づいて前記受信電界強度の減
衰係数が急激に変化するブレークポイントの位置を推定
する電波伝搬推定方法であって、（ａ）前記ブレークポ
イントの位置を仮定するステップ、（ｂ）仮定された前
記ブレークポイントの位置により分割された各領域の前
記受信電界強度データであって、前記ブレークポイント
から数えてほぼ同数の受信電界強度データとなる範囲の
受信電界強度データに対して、それぞれ関数近似処理を
施して適切な近似線を決定し、そのときの近似誤差を求
めるステップ、（ｃ）仮定された前記ブレークポイント
の位置を変化させて、前記ステップ（ｂ）を順次繰り返
すステップ、（ｄ）仮定された前記ブレークポイントの
位置の各々に対して、前記各領域ごとの前記近似誤差に
ついて重み付け総和を計算して評価値とするステップ、
（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップを有するものである。

【００１３】また、請求項７に記載の発明においては、
送受信局間の伝搬距離に応じて受信電界強度をプロット
した受信電界強度データの群に基づいて前記受信電界強
度の減衰係数が急激に変化するブレークポイントの位置
を推定する電波伝搬推定方法であって、（ａ）前記ブレ
ークポイントの位置を仮定するステップ、（ｂ）仮定さ
れた前記ブレークポイントの位置により分割された各領
域の前記受信電界強度データであって、前記ブレークポ
イントから数えてほぼ同数の受信電界強度データとなる
範囲の受信電界強度データに対して、それぞれ関数近似
処理を施して適切な近似線を決定し、そのときの近似誤
差を求めるステップ、（ｃ）前記各領域ごとの前記近似
誤差について重み付け総和を計算して評価値とするステ
ップ、（ｄ）仮定された前記ブレークポイントの位置を
変化させて、前記ステップ（ｂ）〜（ｃ）を順次繰り返
すステップ、（ｅ）仮定された前記ブレークポイントの
位置の各々に対応する前記評価値の大きさに基づいて前
記ブレークポイントの位置を識別するステップを有する
ものである。

【００１４】したがって、請求項６または請求項７に記
載の発明によれば、データの解析者の主観に依存するこ
となく、ブレークポイントの位置を推定することができ
る。関数近似処理に用いるデータとして、ブレークポイ
ントから数えてほぼ同数のデータとなる範囲のデータと
したことから、隣接する領域での関数近似処理精度にば
らつきが生じる可能性が少なくなるとともに、ブレーク
ポイントの位置を推定する上で影響が少ないデータから
関数近似に使用しないようにすることから関数近似処理
の信頼性が増すことになる。

【００１５】請求項８に記載の発明においては、送受信
局間の伝搬距離に応じて受信電界強度をプロットした受
信電界強度データの群に基づいて前記受信電界強度の減
衰係数が急激に変化するブレークポイントの位置を推定
する電波伝搬推定方法であって、（ａ）前記ブレークポ
イントの位置を仮定するステップ、（ｂ）仮定された前
記ブレークポイントにより分割された各領域の前記受信
電界強度データに対して、それぞれ関数近似処理を施し
て適切な近似線を決定し、そのときの近似誤差を求める
ステップ、（ｃ）仮定された前記ブレークポイントの位
置を変化させて、前記ステップ（ｂ）を順次繰り返すス
テップ、（ｄ）前記ステップ（ｂ）において求めた各領
域ごとの前記近似誤差について、仮定された前記ブレー
クポイントの位置を変化させた範囲全体にわたって前記
近似誤差の総和を求め、前記総和の逆数を各領域ごとに
与えられる重み値として設定するステップ、（ｅ）仮定
された前記ブレークポイントの位置の各々に対して、前
記各領域ごとの前記近似誤差について前記重み値により
重み付けられた総和を計算して評価値とするステップ、
（ｆ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップを有するものである。
したがって、データの解析者の主観に依存することな
く、ブレークポイントの位置を推定することができる。
評価関数を正規化するための重み値を演算で求めること
ができるようにしたことから、重み値を試行錯誤で決定
する必要がなくなる。

【００１６】請求項９に記載の発明においては、送受信
局間の伝搬距離に応じて受信電界強度をプロットした受
信電界強度データの群に基づいて前記受信電界強度の減
衰係数が急激に変化するブレークポイントの位置を推定
する電波伝搬推定方法であって、（ａ）前記ブレークポ
イントの位置を仮定するステップ、（ｂ）仮定された前
記ブレークポイントの位置を境界として送信局に近い側
の領域では減衰係数を自由空間伝搬の減衰係数値に固定
し、遠い側の領域では減衰係数の値も関数近似処理での
フィッティング対象のパラメータとし、仮定された前記
ブレークポイントの位置により分割された各領域の前記
受信電界強度データに対して、それぞれ関数近似処理を
施して適切な近似線を決定し、そのときの近似誤差を求
めるステップ、（ｃ）仮定された前記ブレークポイント
の位置を変化させて、前記ステップ（ｂ）を順次繰り返
すステップ、（ｄ）仮定された前記ブレークポイントの
位置の各々に対して、前記各領域ごとの前記近似誤差に
ついて重み付け総和を計算して評価値とするステップ、
（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップを有するものである。

【００１７】また、請求項１０に記載の発明において
は、送受信局間の伝搬距離に応じて受信電界強度をプロ
ットした受信電界強度データの群に基づいて前記受信電
界強度の減衰係数が急激に変化するブレークポイントの
位置を推定する電波伝搬推定方法であって、（ａ）前記
ブレークポイントの位置を仮定するステップ、（ｂ）仮
定された前記ブレークポイントの位置を境界として送信
局に近い側の領域では減衰係数を自由空間伝搬の減衰係
数値に固定し、遠い側の領域では減衰係数の値も関数近
似処理でのフィッティング対象のパラメータとし、仮定
された前記ブレークポイントの位置により分割された各
領域の前記受信電界強度データに対して、それぞれ関数
近似処理を施して適切な近似線を決定し、そのときの近
似誤差を求めるステップ、（ｃ）前記各領域ごとの前記
近似誤差について重み付け総和を計算して評価値とする
ステップ、（ｄ）仮定された前記ブレークポイントの位
置を変化させて、前記ステップ（ｂ）〜（ｃ）を順次繰
り返すステップ、（ｅ）仮定された前記ブレークポイン
トの位置の各々に対応する前記評価値の大きさに基づい
て前記ブレークポイントの位置を識別するステップを有
するものである。

【００１８】したがって、請求項９または請求項１０に
記載の発明によれば、データの解析者の主観に依存する
ことなく、ブレークポイントの位置を推定することがで
きる。ブレークポイント以遠の領域では減衰係数値が大
きくばらつくため、これを関数近似処理においてフィッ
ティングパラメータとすることにより関数近似精度が向
上する。

【００１９】請求項１１に記載の発明においては、送受
信局間の伝搬距離に応じて受信電界強度をプロットした
受信電界強度データの群に基づいて前記受信電界強度の
減衰係数が急激に変化するブレークポイントの位置を推
定する電波伝搬推定方法であって、（ａ）前記ブレーク
ポイントの位置を仮定するステップ、（ｂ）仮定された
前記ブレークポイントの位置により分割された各領域の
前記受信電界強度データに対して、それぞれ関数近似処
理を施して適切な近似線を決定し、そのときの近似誤差
を求めるステップ、（ｃ）仮定された前記ブレークポイ
ントの位置を変化させて、前記ステップ（ｂ）を順次繰
り返すステップ、（ｄ）仮定された前記ブレークポイン
トの位置の各々に対して、前記各領域ごとの前記近似誤
差について重み付け総和を計算して評価値とするステッ
プ、（ｅ）仮定された前記ブレークポイントの位置の各
々に対応する前記評価値の大きさに基づいて、前記評価
値の最小値を求め、該最小値よりも所定の微少値分だけ
大きな評価値となる２個所の位置を求め、前記２個所の
位置の間の中間の位置を前記ブレークポイントの位置と
して識別するステップを有するものである。

【００２０】また、請求項１２に記載の発明において
は、送受信局間の伝搬距離に応じて受信電界強度をプロ
ットした受信電界強度データの群に基づいて前記受信電
界強度の減衰係数が急激に変化するブレークポイントの
位置を推定する電波伝搬推定方法であって、（ａ）前記
ブレークポイントの位置を仮定するステップ、（ｂ）仮
定された前記ブレークポイントの位置により分割された
各領域の前記受信電界強度データに対して、それぞれ関
数近似処理を施して適切な近似線を決定し、そのときの
近似誤差を求めるステップ、（ｃ）前記各領域ごとの前
記近似誤差について重み付け総和を計算して評価値とす
るステップ、（ｄ）仮定された前記ブレークポイントの
位置を変化させて、前記ステップ（ｂ）〜（ｃ）を順次
繰り返すステップ、（ｅ）仮定された前記ブレークポイ
ントの位置の各々に対応する前記評価値の大きさに基づ
いて、前記評価値の最小値を求め、該最小値よりも所定
の微少値分だけ大きな評価値となる２個所の位置を求
め、前記２個所の位置の間の中間の位置を前記ブレーク
ポイントの位置として識別するステップを有するもので
ある。

【００２１】したがって、請求項１１または請求項１２
に記載の発明によれば、データの解析者の主観に依存す
ることなく、ブレークポイントの位置を推定することが
できる。評価関数上の２個所の位置に基づいてブレーク
ポイントの位置を推定するため、評価値に微小な変動が
あっても、この影響を受けにくい。特に、評価値の最小
値近傍が平坦である場合に効果が大きい。

【００２２】請求項１３に記載の発明においては、電波
伝搬推定装置において、送受信局間の伝搬距離に応じて
受信電界強度をプロットした受信電界強度データ群を入
力するための入力装置と、前記入力装置により入力され
た前記受信電界強度データ群を記憶する記録装置と、前
記記録装置に記憶された前記受信電界強度データ群か
ら、前記受信電界強度の減衰係数が急激に変化するブレ
ークポイントの位置を識別する電波伝搬推定演算を実行
する演算処理装置と、前記演算処理装置により実行され
た電波伝搬推定演算の結果を出力する出力装置とを有す
る電波伝搬推定装置であって、前記演算処理装置は、ブ
レークポイントの位置を仮定し、仮定された前記ブレー
クポイントの位置を変化させる処理と、仮定された前記
ブレークポイントの位置の各々に関し、仮定された前記
ブレークポイントの位置で分割された各領域の前記受信
電界強度データに対して、それぞれ関数近似処理を施し
て近似線を決定し、そのときの近似誤差を求める処理
と、仮定された前記ブレークポイントの位置の各々に関
し、前記各領域ごとの前記近似誤差について重み付け総
和を計算して評価値とする処理と、仮定された前記ブレ
ークポイントの位置の各々に関する前記評価値の大きさ
に基づいて前記ブレークポイントの位置を識別する処理
を実行するものである。したがって、データの解析者の
主観に依存することなく、ブレークポイントの位置を推
定することができる。

【００２３】

【発明の実施の形態】本発明のデータ群の特性境界識別
方法の実施の一形態を電波伝搬推定方法および装置を一
例として説明する。図１は、本発明の実施の一形態の電
波伝搬推定方法および装置の説明図である。図１（ａ）
は電波伝搬推定装置のシステム構成図、図１（ｂ）は電
波伝搬推定方法の処理手順の概要を示すフローチャート
である。図中、１は入力装置、２は演算処理装置、３は
記録装置、４は出力装置である。

【００２４】図１（ａ）において、送受信局間の伝搬距
離に応じて受信電界強度をプロットした受信電界強度デ
ータの群を入力装置１から演算処理装置２に入力して、
適宜記録装置３に記憶させる。そして、演算処理装置２
により電波伝搬推定演算を実行し、その結果を記録装置
３に記録するとともに出力装置４へ出力する。図１
（ｂ）を用いて、演算処理装置２における電波伝搬推定
の処理手順の概要を説明する。Ｓ１１において、測定に
より得られた受信電界強度データを入力する。この入力
データは送信局を基準とした送受信局間の距離と、各距
離での受信局の位置における伝搬損失値またはその短区
間中央値とする。なお、伝搬損失値の代わりに受信電力
値を用いても本質的に差異はない。

【００２５】Ｓ１２において、評価のために仮定する特
性の変化する境界、すなわち、評価するブレークポイン
ト位置ｄ_i（ｉ＝１〜Ｎ）を決め、送信局に対してブレ
ークポイント以前となる領域Ｒ_j（ｊ＝１）と、ブレー
クポイント以遠となる領域Ｒ_j（ｊ＝２）とに分ける。
そして、各領域ごとに対数関数を用いて関数近似を行
い、その時の平均二乗誤差ｅ_j（ｄ_i）（ｊ＝１，
２），（ｉ＝１〜Ｎ）を計算する処理を、評価するブレ
ークポイント位置ｄ_iを変化させて行う。

【００２６】次に、Ｓ１３において、Ｓ１２において得
られた平均二乗誤差ｅ_j（ｄ_i）を用いて評価関数Ｅ
（ｄ_i）を、各ブレークポイント位置ｄ_iごとに計算す
る。この評価関数Ｅ（ｄ_i）は、平均二乗誤差ｅ_j（ｄ
_i）に重み付け値ｗ_jを付けた加算であり、次の（１）
式で表され、この重み付け値ｗ_jを付けることにより正
規化される。

【００２７】

【数１】なお、この式（１）は、特性境界を１個所のブレークポ
イント位置として領域を２分割した場合の式である。一
般に、（Ｍ−１）個の特性境界を設けて領域をＭ分割す
る場合には、次の（２）式で表される。

【数２】Ｓ１４においては、評価関数Ｅ（ｄ_i）に基づいて適切
なブレークポイント位置ｄ_cを識別判定し、Ｓ１５にお
いて演算結果を出力する。

【００２８】図２は、図１に示した電波伝搬推定方法お
よび装置におけるブレークポイントの識別判定の説明図
である。図中、横軸は送信局から受信局までの距離ｄ、
縦軸は伝搬損失値ｐである。この伝搬損失値ｐとして
は、瞬時値データ系列に対して適切な区間で移動平均処
理を施して得られる中央値データ（短区間中央値）を用
いる。ｘ軸は対数、ｙ軸はデシベルで表示され、それゆ
え伝搬損失値ｐを距離ｄの対数関数を用いて近似したと
き、近似線は直線になる。

【００２９】ブレークポイント位置ｄ_iに対して、左側
の領域Ｒ_j（Ｍ−１個のブレークポイントの具体例で
は、ｊ＝１）は、送信局に近い側の領域であり、関数近
似処理を施すことにより最適に近似された近似線Ｐ_jお
よびそのときの平均二乗誤差ｅ_j（ｄ_i）が求められ
る。一方、左側の領域Ｒ_j+1に対しても同様の操作を適
用して、最適に近似された近似線Ｐ_j+1およびそのとき
の平均二乗誤差ｅ_j+1（ｄ_i）が求められる。

【００３０】図３は、図１（ｂ）のＳ１２における関数
近似処理手順を説明するためのフローチャートである。
Ｓ２１において、変数ｉ（ｉ＝１〜Ｎ）を１に初期化
し、ブレークポイント位置として仮定する候補の初期値
ｄ₁と、仮定するブレークポイント位置を変化させて調
査する回数Ｎを決定する。Ｓ２２においては、仮定する
ブレークポイント位置ｄ_iよりも送信局側の領域Ｒ_jの
データ群と、仮定するブレークポイント位置ｄ_iよりも
遠い領域Ｒ_j+1のデータ群の各々のデータ群について、
個別に対数関数を用いて関数近似処理を行い、平均二乗
誤差ｅ_j（ｄ_i）を計算する。

【００３１】Ｓ２３においては変数ｉを１加算する。引
き続くＳ２４において調査数Ｎ以下であれば、Ｓ２２に
処理を戻し、新たに仮定したブレークポイント位置ｄ_i
で同様に関数近似処理を行い、Ｓ２４においてｉがＮよ
り大きいと判断されるまで繰り返す。ブレークポイント
位置ｄ_iの候補の決め方としては、初期値ｄ₁を送信局
近傍の、ある距離Ｌ₀の位置に設定し、次の候補ｄ₂は
送信局から微少距離Δだけさらに離れたＬ₀＋Δとし、
そしてｉ番目の候補ｄ_iはＬ₀＋（ｉ−１）・Δと設定
するやり方が考えられる。

【００３２】図４は、図１（ｂ）のＳ１３において評価
関数Ｅ（ｄ_i）を求める式の説明図である。図中、３
１，３２は乗算器、３３は加算器である。左上の線図
は、評価するブレークポイントｄ_iよりも送信局側の領
域Ｒ₁を関数近似する式を領域Ｒ_iのデータ（データの
位置ｄをｄ_kとすると、ｄ_k≦ｄ_iのデータ）について
用いたときの平均二乗誤差ｅ₁（ｄ_i）の変化を表す図
である。ブレークポイントｄ₁からｄ_NまでＮ回実施し
た結果を示している。評価するブレークポイント位置ｄ
_iが正しいブレークポイント位置よりも送信局側にある
場合には、平均二乗誤差ｅ₁（ｄ_i）が比較的小さくな
る。そして、評価するブレークポイント位置ｄ_iが正し
いブレークポイント位置よりも遠くなって自由空間伝搬
の仮定から外れたデータまで含めて関数近似処理した場
合には、誤差が増大して行く傾向が現われる。

【００３３】一方、左下の線図は、評価するブレークポ
イント位置ｄ_iより遠方の領域Ｒ₂を関数近似処理する
式を領域Ｒ_i+1のデータ（データの位置ｄをｄ_kとする
と、ｄ_k＞ｄ_iのデータ）について用いたときの平均二
乗誤差ｅ₂（ｄ_i）の変化を表す図である。ブレークポ
イント位置ｄ₁からｄ_NまでＮ回実施した結果を示して
いる。この場合、評価するブレークポイント位置ｄ_iが
送信局に近いときには、自由空間伝搬に従うデータも含
めて関数近似処理を行うため、平均二乗誤差ｅ
₂（ｄ_i）が大きくなる。逆に、評価するブレークポイ
ント位置ｄ_iが送信局から遠くなり自由空間伝搬に従う
データが除かれて行くと誤差が小さくなる。これらの平
均二乗誤差ｅ₁（ｄ_i）とｅ₂（ｄ_i）とを、乗算器３
１，３３において、式（１）に従い、重み付け値ｗ₁，
ｗ₂を乗算した後、加算器３２において重み付け加算処
理して、右中央の線図に示す評価関数Ｅ（ｄ_i）とす
る。なお、図１（ｂ）のＳ１３における評価関数Ｅ（ｄ
_i）の計算は、重み付け値ｗ₁，ｗ₂が予め決められて
いる場合に、図１（ｂ）のＳ１２と同時に行うこともで
きる。すなわち、図３のステップＳ２２とＳ２３の間に
おいて、ｄ_iごとに行うこともできる。

【００３４】しかし、上述した電波伝搬推定方法および
装置には、次のような改善の余地がある。（１）図１（ｂ）のＳ１２において、関数近似処理に用
いる総データ個数としては、領域Ｒ_jと領域Ｒ_j+1では
必ずしも一致しない。関数近似処理に用いるデータ量が
特性の変化する境界であるブレークポイント以前と以遠
とで異なっていると、関数近似処理精度にばらつきが生
じる可能性がある。（２）図１（ｂ）のＳ１２において、関数近似処理に用
いる減衰係数値として一律に固定値を設定すると、実測
データの減衰係数値は、ブレークポイント以遠の領域Ｒ
_j+1においては、大きくばらつくため、必ずしも最適な
近似にはならない。（３）図１（ｂ）のＳ１３において、評価関数Ｅ
（ｄ_i）を求める際の重み付け値ｗ_jを試行錯誤で決定
すると、そのチューニングに手間がかかる。（４）図１（ｂ）のＳ１４において、ブレークポイント
位置の識別処理では、評価関数Ｅ（ｄ_i）が最小値とな
る位置をブレークポイントと判定していたため、評価値
Ｅ（ｄ_i）の微小な変動により誤差を生じやすい。

【００３５】最初に、上述した（１）の関数近似処理精
度の改善について説明する。図２において、領域Ｒ_jに
おけるデータ数ｎ_jと領域Ｒ_j+1におけるデータ数ｎ
_j+1とを比較し、少ないデータ数の方の領域では全デー
タを用いて関数近似処理を施し、他方の領域ではブレー
クポイント位置ｄ_iから始めて上述した少ない総データ
数と等しいデータ数あるいは、ほぼ同数となるまでのデ
ータだけを用いて関数近似処理を行うこととする。

【００３６】この図では、領域Ｒ_jおよび領域Ｒ_j+1に
おいて楕円枠で囲んだ範囲のデータ数ｎが、ブレークポ
イント位置ｄ_iを挟んで等しくなり、この範囲内のデー
タに対して関数近似処理を施す。隣接した領域Ｒ_j，領
域Ｒ_j+1で等しいデータ個数ｎを用いて各々における関
数近似処理を行うことにより、各領域Ｒ_j，領域Ｒ_j+1
での誤差が同等に得られ、それらを総合した評価関数Ｅ
（ｄ_i）が適切に求められる。ブレークポイント位置ｄ
_iの近傍のデータは、ブレークポイントを推定する上で
影響が大きい。また、送信局から遠ざかった位置のデー
タは、受信電界強度が低下し、ノイズ信号成分が多く含
まれることになるので、上述したように、ブレークポイ
ント位置ｄ_iから取り込むデータ範囲を決めて行くこと
により推定の信頼性が増すことになる。

【００３７】なお、少ないデータ個数の方の領域でも、
必ずしも全てのデータを用いることはなく、隣接する各
領域Ｒ_j，領域Ｒ_j+1でブレークポイントｄ_iの境界か
ら数えて同数のデータ数となるまでの範囲のデータを用
いて関数近似処理を施せばよい。

【００３８】次に、上述した（２）の近似関数の設定の
改善について説明する。送信局に近い側の領域では障害
物が無く、ほぼ自由空間伝搬とみなせるので減衰係数を
自由空間伝搬の減衰係数である２乗に固定して関数近似
処理を行い、逆側の領域では一般に交通トラフィック等
に応じて見通し率が劣化するので減衰係数の値が大きく
ばらつくため、減衰係数の値自体も関数近似処理のフィ
ッティングパラメータとする。すなわち、図１（ｂ）の
ステップ（ａ）において、ブレークポイント位置ｄ_iを
境界として送信局に近い側の領域Ｒ_j（ｊ＝１）では減
衰係数を自由空間伝搬の減衰係数値に固定して関数近似
処理を行い、遠い側の領域Ｒ_j+1（ｊ＝１）では減衰係
数の値も関数近似処理でのフィッティング対象のパラメ
ータとする。送信局に近い側の領域Ｒ_j（ｊ＝１）での
近似式は、デシベル表現で次の（３）式となる。

【００３９】

【数３】また、遠い側の領域Ｒ_j+1（ｊ＝１）での近似式は、同
じくデシベル表現で次の（４）式となる。ここでｃ₁、
ｃ₂、ｃ₃は、関数近似のパラメータである。

【数４】

【００４０】ただし、ブレークポイント位置ｄ_iにおい
て、両近似式が交差することが望ましいから、両近似式
の値がｄ＝ｄ_iにおいて等しくなるように関数近似する
ことが望ましい。

【００４１】式（１）に、式（３），式（４）の近似式
を結びつける式は、次の（５）式である。各領域Ｒ₁，
Ｒ₂にけるデータ数は既に等しくｎ個としている。

【数５】

【００４２】次に、上述した（３）の重み付け値ｗ_jの
決定の改善について説明する。本発明の実施の一形態で
は、重み付け値ｗ_jとして、評価する全てのブレークポ
イントｄ_iに対する各平均二乗誤差ｅ₁（ｄ_i），ｅ₂
（ｄ_i）の総和の逆数を適用する。重みの計算式を、次
の（６）式に示す。

【数６】この重み付け値の設定により、式（５）の右辺の１項目
と２項目とが正規化されて、評価関数が適切に与えられ
る。重み付け値ｗ_jを試行錯誤で決める場合に比べてそ
のパラメータチューニングが省力化される。

【００４３】図５は、本発明の実施の一形態における、
重み付け総和の計算方法の処理手順を示すフローチャー
トである。図１（ｂ）に示したＳ１２，Ｓ１３の処理ス
テップをこの図５のフローチャートに置き換える。Ｓ４
１においては、ブレークポイント位置を示す変数ｉおよ
び、ブレークポイント位置ｄ_iで分割されたＭ個（Ｍ＝
２）の隣接領域Ｒ_j（ｊ＝１，２）での平均二乗誤差ｅ
_j（ｄ_i）の総和を求めるパラメータＳ_j（ｊ＝１，
２）を初期化する。なお、（Ｍ−１）個の特性境界があ
る一般的な場合、Ｓ_j（ｊ＝１，・・・，Ｍ）は、領域
Ｒ_j（ｊ＝１，・・・，Ｍ）で平均二乗誤差ｅ
_j（ｄ_i）の総和を求めるパラメータとなる。

【００４４】Ｓ４２においては、領域Ｒ_jを示す変数ｊ
を１とする。Ｓ４３においては、仮定するブレークポイ
ント位置ｄ_iを境界として分割される領域Ｒ_jのデータ
群を設定し、Ｓ４４においては、これらのデータを関数
近似して平均二乗誤差ｅ_jを計算する。この計算の際、
上述した（１）の関数近似処理精度の改善および上述し
た（２）の近似関数の設定の改善を行うことができる。
Ｓ４５においては、Ｓ４３で計算した平均二乗誤差ｅ_j
（ｄ_i）を、総和を求めるパラメータＳ_jに加算する。
Ｓ４６においては、ｊに１を加算し、Ｓ４７においては
ｊがＭ（ブレークポイント位置の例では、Ｍ＝２）を超
えないときにはＳ４２に処理を戻し、Ｓ４３〜Ｓ４７の
処理を全領域Ｍについて終了するまで繰り返す。ｊがＭ
を超えたときにはＳ４８に処理を進め、Ｓ４８において
は、ｉを１ずつ加算し、Ｓ４９においては、ｉがＮを超
えないときにはＳ４２に処理を戻し、Ｎを超えたときに
はＳ５０に処理を進める。

【００４５】したがって、評価するブレークポイント位
置ｄ_iがＮ個所になるまで平均二乗誤差ｅ_j（ｄ_i）の
総和を求めるパラメータＳ_j（ｊ＝１〜Ｍ）を調べる。
Ｓ５０においては、式（６）に従い、各領域Ｒ_j（ｊ＝
１〜Ｍ）で得られた平均二乗誤差ｅ_j（ｄ_i）の総和を
求めるパラメータＳ_jの逆数を評価関数の各項の重み付
け値ｗ_jに設定する。Ｓ５１においては、評価関数Ｅ
（ｄ_i）を、式（２）、より具体的には、式（５）に従
って計算する。このようにして、評価関数の重み付け値
の設定は、誤差を正規化する係数とすることで一意に決
定することができる。

【００４６】次に、上述した（４）のブレークポイント
位置の識別処理について説明する。図６は、本発明の実
施の形態におけるブレークポイント位置の判定方法を説
明する線図である。図中、横軸は仮定したブレークポイ
ント位置ｄ_iであり、縦軸は評価関数Ｅ（ｄ_i）の値で
ある。ｄ_mの位置で評価関数Ｅ（ｄ_m）の値が最小値Ｅ
_minになったとする。実測データに基づいてこの評価関
数を計算してみると、このｄ_mの位置の近傍が平坦にな
る場合がある。そのため、ブレークポイント位置を識別
判定するときに、評価関数計算の際の微小な誤差の影響
を受けることになり、必ずしも上述したｄ_mの位置が正
しいブレークポイント位置とはならない。

【００４７】そこで、最小値Ｅ_minよりも適切な微少値
ΔＥだけ大きな評価値を決め、この評価値Ｅ_min＋ΔＥ
となる位置をｄ_mの前後で1個所ずつ求め、その２個所
ｄ_a，ｄ_bの中心位置ｄ_c＝（ｄ_a＋ｄ_b）／２を正し
いブレークポイントの位置であると判定する。２個所を
基準にとって判定処理を行うため、より高精度な識別が
可能となる。この微少値ΔＥとしては、例えば、評価関
数Ｅ（ｄ_i）の最大値の約５％程度とした。

【００４８】図７は、図６を参照して説明したブレーク
ポイント位置の判定方法の処理を説明するためのフロー
チャートである。Ｓ６１においては、評価関数Ｅ
（ｄ_i）の最小値Ｅ_minを探索する。次に、Ｓ６２にお
いては、適切な微小値ΔＥを決定し、Ｓ６３において
は、評価値がＥ_min＋ΔＥとなる位置ｄ_a及びｄ_bを調
べる。Ｓ６４においては、ｄ_aとｄ_b の中心位置（ｄ
_a＋ｄ_b）／２を計算して、これをブレークポイント位
置ｄ_cとする。

【００４９】上述した説明では、送受信局間の伝搬距離
に応じて受信電界強度を測定値に基づいてプロットした
受信電界強度データに基づいて、前記受信電界強度の減
衰係数が急激に変化する１個所のブレークポイントの位
置を識別する電波伝搬推定方法および装置を例にして説
明した。しかし、（Ｍー１）個所の特性の変化する境界
を有する一般的な場合に、この境界および関数近似式の
推定についても同様な処理が適用可能である。

【００５０】

【発明の効果】上述した説明から明らかなように、本発
明によれば、ブレークポイント位置などの特性境界の推
定を人の主観に頼らず高精度かつ容易に行うことができ
るという効果がある。さらに、関数近似精度のばらつき
の減少、関数近似精度の向上、評価関数の微小な誤差の
影響の軽減などを図ることができる。

【図面の簡単な説明】

【図１】本発明の実施の一形態の電波伝搬推定方法およ
び装置の説明図である。図１（ａ）は電波伝搬推定装置
のシステム構成図、図１（ｂ）は電波伝搬推定方法の処
理手順の概要を示すフローチャートである。

【図２】図１に示した電波伝搬推定方法および装置にお
けるブレークポイントの識別判定の説明図である。

【図３】図１（ｂ）のＳ１２における関数近似処理手順
を説明するためのフローチャートである。

【図４】図１（ｂ）のＳ１３において評価関数Ｅ
（ｄ_i）を求める式の説明図である。

【図５】本発明の実施の一形態における、重み付け総和
の計算方法の処理手順を示すフローチャートである。

【図６】本発明の実施の形態における、ブレークポイン
ト位置の判定方法を説明する線図である。

【図７】図６を参照して説明したブレークポイント位置
の判定方法の処理を説明するためのフローチャートであ
る。

【符号の説明】

１入力装置、２演算処理装置、３記録装置、４
出力装置、３１，３２乗算器、３３加算器

───────────────────────────────────────────────────── フロントページの続き (56)参考文献特開平８−163012（ＪＰ，Ａ) 特開平10−200429（ＪＰ，Ａ) (58)調査した分野(Int.Cl.⁶，ＤＢ名) H04B 17/00

Claims

(57)【特許請求の範囲】

【請求項１】特性の変化する境界を有するデータの群
から、前記境界を識別するための方法であって、次の
（ａ）〜（ｅ）のステップを有するデータ群の特性境界
識別方法。（ａ）前記特性の変化する境界を仮定するステップ、（ｂ）仮定された前記境界により分割された各領域の前
記データであって、前記境界から数えてほぼ同数のデー
タとなる範囲のデータに対して、それぞれ関数近似処理
を施して適切な近似線を決定し、そのときの近似誤差を
求めるステップ、（ｃ）仮定された前記境界を変化させて、前記ステップ
（ｂ）を順次繰り返すステップ、（ｄ）仮定された前記境界の各々に対して、前記各領域
ごとの前記近似誤差について重み付け総和を計算して評
価値とするステップ、（ｅ）仮定された前記境界の各々に対応する前記評価値
の大きさに基づいて前記特性の変化する境界の値を識別
するステップ。
【請求項２】特性の変化する境界を有するデータの群
から、前記境界を識別するための方法であって、次の
（ａ）〜（ｅ）のステップを有するデータ群の特性境界
識別方法。（ａ）前記特性の変化する境界を仮定するステップ、（ｂ）仮定された前記境界により分割された各領域の前
記データであって、前記境界から数えてほぼ同数のデー
タとなる範囲のデータに対して、それぞれ関数近似処理
を施して適切な近似線を決定し、そのときの近似誤差を
求めるステップ、（ｃ）前記各領域ごとの前記近似誤差について重み付け
総和を計算して評価値とするステップ、（ｄ）仮定された前記境界を変化させて、前記ステップ
（ｂ）〜（ｃ）を順次繰り返すステップ、（ｅ）仮定された前記境界の各々に対応する前記評価値
の大きさに基づいて前記特性の変化する境界の値を識別
するステップ。
【請求項３】特性の変化する境界を有するデータの群
から、前記境界を識別するための方法であって、次の
（ａ）〜（ｆ）のステップを有するデータ群の特性境界
識別方法。（ａ）前記特性の変化する境界を仮定するステップ、（ｂ）仮定された前記境界により分割された各領域の前
記データに対して、それぞれ関数近似処理を施して適切
な近似線を決定し、そのときの近似誤差を求めるステッ
プ、（ｃ）仮定された前記境界を変化させて、前記ステップ
（ｂ）を順次繰り返すステップ、（ｄ）前記ステップ（ｂ）において求めた各領域ごとの
前記近似誤差について、仮定された前記境界を変化させ
た範囲全体にわたって前記近似誤差の総和を求め、前記
総和の逆数を各領域ごとに与えられる重み値として設定
するステップ、（ｅ）仮定された前記境界の各々に対して、前記各領域
ごとの前記近似誤差について前記重み値により重み付け
られた総和を計算して評価値とするステップ、（ｆ）仮定された前記境界の各々に対応する前記評価値
の大きさに基づいて前記特性の変化する境界の値を識別
するステップ。
【請求項４】送受信局間の伝搬距離に応じて受信電界
強度をプロットした受信電界強度データの群に基づいて
前記受信電界強度の減衰係数が急激に変化するブレーク
ポイントの位置を推定する電波伝搬推定方法であって、
次の（ａ）〜（ｅ）のステップを有する電波伝搬推定方
法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置により分
割された各領域の前記受信電界強度データに対して、そ
れぞれ関数近似処理を施して適切な近似線を決定し、そ
のときの近似誤差を求めるステップ、（ｃ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）を順次繰り返すステップ、（ｄ）仮定された前記ブレークポイントの位置の各々に
対して、前記各領域ごとの前記近似誤差について重み付
け総和を計算して評価値とするステップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項５】送受信局間の伝搬距離に応じて受信電界
強度をプロットした受信電界強度データの群に基づいて
前記受信電界強度の減衰係数が急激に変化するブレーク
ポイントの位置を推定する電波伝搬推定方法であって、
次の（ａ）〜（ｅ）のステップを有する電波伝搬推定方
法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置により分
割された各領域の前記受信電界強度データに対して、そ
れぞれ関数近似処理を施して適切な近似線を決定し、そ
のときの近似誤差を求めるステップ、（ｃ）前記各領域ごとの前記近似誤差について重み付け
総和を計算して評価値とするステップ、（ｄ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）〜（ｃ）を順次繰り返すステ
ップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項６】送受信局間の伝搬距離に応じて受信電界
強度をプロットした受信電界強度データの群に基づいて
前記受信電界強度の減衰係数が急激に変化するブレーク
ポイントの位置を推定する電波伝搬推定方法であって、
次の（ａ）〜（ｅ）のステップを有する電波伝搬推定方
法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置により分
割された各領域の前記受信電界強度データであって、前
記ブレークポイントから数えてほぼ同数の受信電界強度
データとなる範囲の受信電界強度データに対して、それ
ぞれ関数近似処理を施して適切な近似線を決定し、その
ときの近似誤差を求めるステップ、（ｃ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）を順次繰り返すステップ、（ｄ）仮定された前記ブレークポイントの位置の各々に
対して、前記各領域ごとの前記近似誤差について重み付
け総和を計算して評価値とするステップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項７】送受信局間の伝搬距離に応じて受信電界
強度をプロットした受信電界強度データの群に基づいて
前記受信電界強度の減衰係数が急激に変化するブレーク
ポイントの位置を推定する電波伝搬推定方法であって、
次の（ａ）〜（ｅ）のステップを有する電波伝搬推定方
法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置により分
割された各領域の前記受信電界強度データであって、前
記ブレークポイントから数えてほぼ同数の受信電界強度
データとなる範囲の受信電界強度データに対して、それ
ぞれ関数近似処理を施して適切な近似線を決定し、その
ときの近似誤差を求めるステップ、（ｃ）前記各領域ごとの前記近似誤差について重み付け
総和を計算して評価値とするステップ、（ｄ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）〜（ｃ）を順次繰り返すステ
ップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項８】送受信局間の伝搬距離に応じて受信電界
強度をプロットした受信電界強度データの群に基づいて
前記受信電界強度の減衰係数が急激に変化するブレーク
ポイントの位置を推定する電波伝搬推定方法であって、
次の（ａ）〜（ｆ）のステップを有する電波伝搬推定方
法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントにより分割され
た各領域の前記受信電界強度データに対して、それぞれ
関数近似処理を施して適切な近似線を決定し、そのとき
の近似誤差を求めるステップ、（ｃ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）を順次繰り返すステップ、（ｄ）前記ステップ（ｂ）において求めた各領域ごとの
前記近似誤差について、仮定された前記ブレークポイン
トの位置を変化させた範囲全体にわたって前記近似誤差
の総和を求め、前記総和の逆数を各領域ごとに与えられ
る重み値として設定するステップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対して、前記各領域ごとの前記近似誤差について前記重
み値により重み付けられた総和を計算して評価値とする
ステップ、（ｆ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項９】送受信局間の伝搬距離に応じて受信電界
強度をプロットした受信電界強度データの群に基づいて
前記受信電界強度の減衰係数が急激に変化するブレーク
ポイントの位置を推定する電波伝搬推定方法であって、
次の（ａ）〜（ｅ）のステップを有する電波伝搬推定方
法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置を境界と
して送信局に近い側の領域では減衰係数を自由空間伝搬
の減衰係数値に固定し、遠い側の領域では減衰係数の値
も関数近似処理でのフィッティング対象のパラメータと
し、仮定された前記ブレークポイントの位置により分割
された各領域の前記受信電界強度データに対して、それ
ぞれ関数近似処理を施して適切な近似線を決定し、その
ときの近似誤差を求めるステップ、（ｃ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）を順次繰り返すステップ、（ｄ）仮定された前記ブレークポイントの位置の各々に
対して、前記各領域ごとの前記近似誤差について重み付
け総和を計算して評価値とするステップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項１０】送受信局間の伝搬距離に応じて受信電
界強度をプロットした受信電界強度データの群に基づい
て前記受信電界強度の減衰係数が急激に変化するブレー
クポイントの位置を推定する電波伝搬推定方法であっ
て、次の（ａ）〜（ｅ）のステップを有する電波伝搬推
定方法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置を境界と
して送信局に近い側の領域では減衰係数を自由空間伝搬
の減衰係数値に固定し、遠い側の領域では減衰係数の値
も関数近似処理でのフィッティング対象のパラメータと
し、仮定された前記ブレークポイントの位置により分割され
た各領域の前記受信電界強度データに対して、それぞれ
関数近似処理を施して適切な近似線を決定し、そのとき
の近似誤差を求めるステップ、（ｃ）前記各領域ごとの前記近似誤差について重み付け
総和を計算して評価値とするステップ、（ｄ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）〜（ｃ）を順次繰り返すステ
ップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて前記ブレークポ
イントの位置を識別するステップ。
【請求項１１】送受信局間の伝搬距離に応じて受信電
界強度をプロットした受信電界強度データの群に基づい
て前記受信電界強度の減衰係数が急激に変化するブレー
クポイントの位置を推定する電波伝搬推定方法であっ
て、次の（ａ）〜（ｅ）のステップを有する電波伝搬推
定方法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置により分
割された各領域の前記受信電界強度データに対して、そ
れぞれ関数近似処理を施して適切な近似線を決定し、そ
のときの近似誤差を求めるステップ、（ｃ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）を順次繰り返すステップ、（ｄ）仮定された前記ブレークポイントの位置の各々に
対して、前記各領域ごとの前記近似誤差について重み付
け総和を計算して評価値とするステップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて、前記評価値の
最小値を求め、該最小値よりも所定の微少値分だけ大き
な評価値となる２個所の位置を求め、前記２個所の位置
の間の中間の位置を前記ブレークポイントの位置として
識別するステップ。
【請求項１２】送受信局間の伝搬距離に応じて受信電
界強度をプロットした受信電界強度データの群に基づい
て前記受信電界強度の減衰係数が急激に変化するブレー
クポイントの位置を推定する電波伝搬推定方法であっ
て、次の（ａ）〜（ｅ）のステップを有する電波伝搬推
定方法。（ａ）前記ブレークポイントの位置を仮定するステッ
プ、（ｂ）仮定された前記ブレークポイントの位置により分
割された各領域の前記受信電界強度データに対して、そ
れぞれ関数近似処理を施して適切な近似線を決定し、そ
のときの近似誤差を求めるステップ、（ｃ）前記各領域ごとの前記近似誤差について重み付け
総和を計算して評価値とするステップ、（ｄ）仮定された前記ブレークポイントの位置を変化さ
せて、前記ステップ（ｂ）〜（ｃ）を順次繰り返すステ
ップ、（ｅ）仮定された前記ブレークポイントの位置の各々に
対応する前記評価値の大きさに基づいて、前記評価値の
最小値を求め、該最小値よりも所定の微少値分だけ大き
な評価値となる２個所の位置を求め、前記２個所の位置
の間の中間の位置を前記ブレークポイントの位置として
識別するステップ。
【請求項１３】送受信局間の伝搬距離に応じて受信電
界強度をプロットした受信電界強度データ群を入力する
ための入力装置と、前記入力装置により入力された前記受信電界強度データ
群を記憶する記録装置と、前記記録装置に記憶された前記受信電界強度データ群か
ら、前記受信電界強度の減衰係数が急激に変化するブレ
ークポイントの位置を識別する電波伝搬推定演算を実行
する演算処理装置と、前記演算処理装置により実行された電波伝搬推定演算の
結果を出力する出力装置とを有する電波伝搬推定装置で
あって、前記演算処理装置は、ブレークポイントの位置を仮定し、仮定された前記ブレ
ークポイントの位置を変化させる処理と、仮定された前記ブレークポイントの位置の各々に関し、
仮定された前記ブレークポイントの位置で分割された各
領域の前記受信電界強度データに対して、それぞれ関数
近似処理を施して近似線を決定し、そのときの近似誤差
を求める処理と、仮定された前記ブレークポイントの位置の各々に関し、
前記各領域ごとの前記近似誤差について重み付け総和を
計算して評価値とする処理と、仮定された前記ブレークポイントの位置の各々に関する
前記評価値の大きさに基づいて前記ブレークポイントの
位置を識別する処理を実行するものである、ことを特徴
とする電波伝搬推定装置。