JP4869187B2 - 伝搬経路計算装置、伝搬経路計算方法および記録媒体 - Google Patents
伝搬経路計算装置、伝搬経路計算方法および記録媒体 Download PDFInfo
- Publication number
- JP4869187B2 JP4869187B2 JP2007235627A JP2007235627A JP4869187B2 JP 4869187 B2 JP4869187 B2 JP 4869187B2 JP 2007235627 A JP2007235627 A JP 2007235627A JP 2007235627 A JP2007235627 A JP 2007235627A JP 4869187 B2 JP4869187 B2 JP 4869187B2
- Authority
- JP
- Japan
- Prior art keywords
- propagation path
- zenith angle
- calculation
- refractive index
- electron density
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
本発明は、電離層中を伝搬する電波の伝搬経路を計算するための伝搬経路計算装置、伝搬経路計算方法およびその伝搬経路計算方法を実現するためのプログラムを記録した記録媒体に関する。
従来、HF帯域(10〜30MHz)を用いた通信が行われている。このHF帯域を用いた通信においては、送信元で射出された電波は、地球の上空に存在する電離層で反射され、電離層反射波として見通し距離外の遠方に伝搬する。電波は、一般に、周波数が低いほど電離層で反射しやすくなり、また、電離層を通過または反射するときに減衰する。
ところで、HF帯域を用いた通信は電離層の状態に影響されるので、遠方の相手に電波が届くかどうかを確認するために、電波の出射角度、周波数、時間帯および自己の位置などを変数として、電波の伝搬経路を計算により求めることが行われている。この伝搬経路の計算においては、電波の伝搬経路を辿るための計算を行うポイント(以下、「計算ポイント」という)が順次設定される。ある計算ポイントから次の計算ポイントまでの距離または電波の伝搬時間を「ステップ」と呼ぶ。
なお、関連する技術として、非特許文献1は、電離層における電波の伝搬について開示しており、また、非特許文献2は、電波の伝搬経路を計算する技術を開示している。
"International Reference Ionosphere 2000" Radio Science, Volume 36, Number 2, Page261-275, March/April 2001 "Effects of on Whistle-Mode Ray Tracing" RADIO SCIENCE, Vol 1(New Series), No. 3, Page269-283, March 1966
"International Reference Ionosphere 2000" Radio Science, Volume 36, Number 2, Page261-275, March/April 2001 "Effects of on Whistle-Mode Ray Tracing" RADIO SCIENCE, Vol 1(New Series), No. 3, Page269-283, March 1966
しかしながら、上述した従来の電波の伝搬経路を計算する技術では、電波の伝搬経路のうちの電離層の部分、特に反射点近辺では、ステップを細かくしないと、電離層における反射状態を再現できなくなる。
すなわち、ステップが荒いと、次の計算ポイントが、本来反射しなければならない場所より高い位置に到達してしまい、屈折率が複素数になるなど通常の状態ではなくなる。その一方、反射点近辺の経路を正確に再現するために、伝搬経路の全体にわたってステップを細かくすると、膨大な計算時間が必要になる。
本発明の課題は、電波の反射点近辺における伝搬経路を正確に計算することができ、しかも、計算時間を短くできる伝搬経路計算装置、伝搬経路計算方法およびその伝搬経路計算方法を実現するためのプログラムを記録した記録媒体を提供することにある。
上記課題を解決するために、請求項1記載の発明に係る伝搬経路計算装置は、電波の伝搬経路の現在の計算ポイントにおける電子密度を求める電子密度計算部と、電子密度計算部で求められた電子密度に基づき屈折率を求める屈折率計算部と、屈折率計算部で求められた屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求める天頂角計算部と、天頂角計算部で求められた天頂角に応じて、前記現在の計算ポイントから次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めるステップ計算部を備えたことを特徴とする。
また、請求項2記載の発明に係る伝搬経路計算方法は、電波の伝搬経路の現在の計算ポイントにおける電子密度を求め、該求められた電子密度に基づき屈折率を求め、該求められた屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求め、該求められた天頂角に応じて、前記現在の計算ポイントから次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めることを特徴とする。
また、請求項3記載の発明に係るコンピュータ読み取り可能な記録媒体は、電波の伝搬経路の現在の計算ポイントにおける電子密度を求め、該求められた電子密度に基づき屈折率を求め、該求められた屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求め、該求められた天頂角に応じて、前記現在の計算ポイントから次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めるためのプログラムを記憶したことを特徴とする。
請求項1〜請求項3記載の発明によれば、電波の伝搬経路の現在の計算ポイントにおける電子密度を求め、電子密度に基づき屈折率を求め、屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求め、さらに、求めた天頂角に応じて、次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めるので、天頂角が大きくなる電離層の反射点近辺では、ステップが小さくなって細かい計算が可能になり、伝搬経路を正確に計算することができる。その一方、天頂角が小さくなる電離層の反射点近辺以外では、ステップが大きくなるので、計算時間を小さくすることができる。
以下、本発明の実施の形態を、図面を参照しながら詳細に説明する。
図1は、本発明の実施例1に係る伝搬経路計算装置の構成を示す図である。この伝搬経路計算装置は、例えばパーソナルコンピュータのソフトウェア処理によって実現された電子密度計算部11、屈折率計算部12、天頂角計算部13およびステップ計算部14を備えている。
電子密度計算部11は、電波の伝搬経路として設定される各計算ポイントにおける電子密度を計算する。図2(a)は、電離層を形成するD層、E層、F1層およびF2層における電子密度の変化を示す図である。この電子密度計算部11は、電子密度アルゴリズムに従って計算により電子密度を求めるように構成することもできるし、図2(a)に示すような電子密度と高さHとの関係を示すデータをデータベースとして予め格納しておくことにより構成することもできる。
図2(b)は、図2(a)に示すような電子密度を有する電離層を電波が伝搬する様子を示す図である。所定の仰角el及び周波数で射出された電波は、電離層の電子密度に応じた屈折率で経路を変更しながらF2層の下部で反射点に到達し、この反射点で反射された電波は、電離層の電子密度に応じた屈折率で経路を変更しながら伝播し、地表に到達する。この電子密度計算部11で計算された電子密度は、屈折率・伝播方向計算部12に送られる。
屈折率・伝播方向計算部12は、電子密度計算部11から送られてくる電子密度に基づき各計算ポイントにおける屈折率を計算し、屈折率に応じた伝播方向ベクトルを計算する。この屈折率・伝播方向計算部12で計算された屈折率・伝播方向ベクトルは、天頂角計算部13に送られる。
天頂角計算部13は、屈折率・伝播方向計算部12から送られてくる伝播方向ベクトルに基づき、各計算ポイントにおける電波の伝搬経路の天頂角、つまりステップの伝搬方向ベクトルと鉛直方向との成す角を計算する。この天頂角計算部13で計算された天頂角は、ステップ計算部14に送られる。
ステップ計算部14は、天頂角計算部13から送られてくる天頂角に応じて、次の計算ポイントまでのステップ(距離または伝搬時間)を計算する。この際、天頂角が大きくなるに連れて小さくなるように変化するステップが求められる。
具体的には、ステップ計算部14は、(1)式に示すフェルミ分布関数を利用してステップstepを求める。この実施例1に係る伝搬経路計算装置は、フェルミ分布関数をステップの計算に利用したこと、および、フェルミ準位(α)、温度に相当する部分(β)の設定方法に特徴がある。
ここで、
τ0:初期ステップ値
z:ステップの伝搬方向ベクトルと鉛直方向との成す角(天頂角)
α:初期ステップ値を半分にする天頂角(0<α<π/2)
β:stepのα近辺における変化量を決める値
βは下記の(2)式によって求められる。
τ0:初期ステップ値
z:ステップの伝搬方向ベクトルと鉛直方向との成す角(天頂角)
α:初期ステップ値を半分にする天頂角(0<α<π/2)
β:stepのα近辺における変化量を決める値
βは下記の(2)式によって求められる。
ここで、
el:送信(受信)仰角(0≦el≦π/2)
γ:定数(0<γ<π/2)
図3は、(1)式により求められるステップstepが天頂角zの大きさに応じて変化する様子を示す図である。天頂角zが大きくなるとステップstepが小さくなる様子が示されている。
el:送信(受信)仰角(0≦el≦π/2)
γ:定数(0<γ<π/2)
図3は、(1)式により求められるステップstepが天頂角zの大きさに応じて変化する様子を示す図である。天頂角zが大きくなるとステップstepが小さくなる様子が示されている。
図4は、各計算ポイントにおける伝播方向ベクトルと天頂角との関係を示す図であり、仰角elで射出された電波の高さが大きくなると屈折率が大きくなって天頂角zが大きくなる様子(z1<z2<z3<z4・・・)が示されている。
次に、上記のように構成される本発明の実施例1に係る伝搬経路計算装置の動作を、伝搬経路計算処理を中心に、図5に示すフローチャートを参照しながら説明する。なお、図5のフローチャートに示す伝搬経路計算処理は、計算ポイントが電離層に存在する場合の動作を中心に示している。
まず、初期ステップ値τ0が設定される(ステップS11)。この初期ステップ値τ0は、ステップ計算部14に予め格納しておくように構成することもできるし、図示しない入力装置からステップ計算部14に設定するように構成することもできる。
次いで、電子密度が求められる(ステップS13)。電子密度が基準値以上であるかどうかにより電離層に入ったかどうかが調べられる(ステップS12)。すなわち、順次にステップを計算した結果、計算ポイントの高さが電離層の高さに到達したかどうかが調べられる。ステップS12において、電離層内に計算ポイントがないことが判断されると、処理は終了する。この場合、図示は省略するが、電波の伝搬経路が電離層以外である場合の処理、つまり荒いステップでの伝搬経路の計算が行われる。
一方、ステップS12において、電離層内に計算ポイントがあることが判断されると、次いで、屈折率伝播方向ベクトルが求められる(ステップS14)。
すなわち、屈折率・伝播方向計算部12は、電子密度計算部11から送られてくる電子密度に基づき、現在の計算ポイントにおける屈折率・伝播方向ベクトルを計算し、天頂角計算部13に送る。
次いで、伝播方向ベクトルから天頂角が求められる(ステップS15)。すなわち、天頂角計算部13は、屈折率・伝播方向計算部12から送られてくる屈折率および現在の計算ポイントに対する電波の伝播方向ベクトルに基づき、現在の計算ポイントにおける電波の伝搬経路の天頂角を計算し、ステップ計算部14に送る。
次いで、ステップ計算が行われる(ステップS16)。すなわち、ステップ計算部14は、天頂角計算部13から送られてくる天頂角に応じて、現在の計算ポイントから次の計算ポイントまでのステップ(距離または伝搬時間)を計算する。次いで、次の計算ポイントが設定される(ステップS17)。すなわち、ステップ計算部14は、ステップS16で計算されたステップに基づき、次の計算ポイントを設定する。
次いで、電子密度が求められ(ステップS13)、電離層から出たかどうかが調べられる(ステップS12)。ステップS12において、計算ポイントが電離層内にないことが判断されると、ステップS14に戻り、上述した処理が繰り返される。これにより、図4に示すように、ステップstepが順次に計算され、電波の伝搬経路が求められる。
一方、ステップS12において、計算ポイントが電離層にないと判断されると、処理は終了する。この場合、図示は省略するが、電波の伝搬経路が電離層以外である場合の処理が行われる。
以上説明したように、本発明の実施例1に係る伝搬経路計算装置によれば、天頂角が大きくなる電離層の反射点近辺では、ステップが小さくなって細かい計算が可能になり、伝搬経路を正確に計算することができる。その一方、天頂角が小さくなる電離層の反射点近辺以外では、ステップが大きくなるので計算時間を小さくすることができる。
なお、図5のフローチャートで示す手順は、これをコンピュータプログラムで実現してハードディスク等の記録媒体に格納しておき、必要に応じてコンピュータ上で実行するように構成することができる。
本発明は、例えばHF帯域で通信を行う通信装置に適用可能である。
11 電子密度計算部
12 屈折率・伝播方向計算部
13 天頂角計算部
14 ステップ計算部
12 屈折率・伝播方向計算部
13 天頂角計算部
14 ステップ計算部
Claims (3)
- 電波の伝搬経路の現在の計算ポイントにおける電子密度を求める電子密度計算部と、
前記電子密度計算部で求められた電子密度に基づき屈折率を求める屈折率計算部と、
前記屈折率計算部で求められた屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求める天頂角計算部と、
前記天頂角計算部で求められた天頂角に応じて、前記現在の計算ポイントから次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めるステップ計算部と、
を備えたことを特徴とする伝搬経路計算装置。 - 電波の伝搬経路の現在の計算ポイントにおける電子密度を求め、
該求められた電子密度に基づき屈折率を求め、
該求められた屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求め、
該求められた天頂角に応じて、前記現在の計算ポイントから次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めることを特徴とする伝搬経路計算方法。 - 電波の伝搬経路の現在の計算ポイントにおける電子密度を求め、
該求められた電子密度に基づき屈折率を求め、
該求められた屈折率に基づき現在の計算ポイントにおける電波の伝搬経路の天頂角を求め、
該求められた天頂角に応じて、前記現在の計算ポイントから次の計算ポイントまでの距離または伝搬時間を表し且つ天頂角が大きくなるに連れて小さくなるように変化するステップを求めるためのプログラムを記憶したことを特徴とするコンピュータ読み取り可能な記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007235627A JP4869187B2 (ja) | 2007-09-11 | 2007-09-11 | 伝搬経路計算装置、伝搬経路計算方法および記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2007235627A JP4869187B2 (ja) | 2007-09-11 | 2007-09-11 | 伝搬経路計算装置、伝搬経路計算方法および記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2009069982A JP2009069982A (ja) | 2009-04-02 |
JP4869187B2 true JP4869187B2 (ja) | 2012-02-08 |
Family
ID=40606196
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2007235627A Expired - Fee Related JP4869187B2 (ja) | 2007-09-11 | 2007-09-11 | 伝搬経路計算装置、伝搬経路計算方法および記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP4869187B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104132999B (zh) * | 2014-07-31 | 2017-01-04 | 中国人民解放军后勤工程学院 | 基于前期到达波的圆柱形储液罐底板声源定位方法 |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6587403B2 (ja) * | 2014-03-25 | 2019-10-09 | 三菱電機株式会社 | 位置標定方法および位置標定装置 |
EP3692652A4 (en) * | 2017-10-04 | 2021-06-30 | Skywave Networks LLC | SETTING OF GEARS BASED ON DIRECT IONOSPHERIC DETECTION |
WO2020003513A1 (ja) * | 2018-06-29 | 2020-01-02 | 三菱電機株式会社 | レーダ装置 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH04253282A (ja) * | 1991-01-29 | 1992-09-09 | Fujitsu Ltd | 蜃気楼画像生成方法 |
JP2723813B2 (ja) * | 1995-02-16 | 1998-03-09 | 防衛庁技術研究本部長 | 収束帯内目標の水平距離算出方式 |
JP3993557B2 (ja) * | 2003-11-10 | 2007-10-17 | 三菱重工業株式会社 | 電磁界分布シミュレーション方法およびその装置 |
-
2007
- 2007-09-11 JP JP2007235627A patent/JP4869187B2/ja not_active Expired - Fee Related
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104132999B (zh) * | 2014-07-31 | 2017-01-04 | 中国人民解放军后勤工程学院 | 基于前期到达波的圆柱形储液罐底板声源定位方法 |
Also Published As
Publication number | Publication date |
---|---|
JP2009069982A (ja) | 2009-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
GB2532542A (en) | Risk quantification for policy deployment | |
Alvarez et al. | A discontinuous finite element formulation for Helmholtz equation | |
JP2023509572A (ja) | 車両を検出するための方法、装置、電子機器、記憶媒体およびコンピュータプログラム | |
JP4530898B2 (ja) | 電波伝搬の推定プログラム、電波伝搬の推定方法、この方法を実行する装置 | |
JP4869187B2 (ja) | 伝搬経路計算装置、伝搬経路計算方法および記録媒体 | |
TW201610730A (zh) | 點雲網格簡化系統及方法 | |
JP2008124938A (ja) | チャネル特性解析装置及び方法 | |
US20180314978A1 (en) | Learning apparatus and method for learning a model corresponding to a function changing in time series | |
JP5790646B2 (ja) | パープレキシティ算出装置 | |
CN111684541A (zh) | 信息处理装置、程序、工艺处理执行装置及信息处理系统 | |
WO2009073334A2 (en) | Method and apparatus for enhancing the accuracy and speed of a ray launching simulation tool | |
Zhao | Improvement on the numerical dispersion of 2-D ADI-FDTD with artificial anisotropy | |
JP2015162188A (ja) | データ解析装置及び方法 | |
KR101492197B1 (ko) | 순차적 근사 최적 설계 방법 및 장치 | |
JP2011223486A (ja) | 電界強度推定装置、方法及びプログラム | |
US20160148246A1 (en) | Automated System for Safe Policy Improvement | |
US20090018807A1 (en) | Hybrid Method for Enforcing Curvature Related Boundary Conditions in Solving One-Phase Fluid Flow Over a Deformable Domain | |
CN114049530A (zh) | 混合精度神经网络量化方法、装置及设备 | |
US9444778B2 (en) | Method and apparatus for propagating a message in a social network | |
US20090167756A1 (en) | Method and apparatus for computation of wireless signal diffraction in a three-dimensional space | |
KR100961875B1 (ko) | 수직 비행경로 자동생성 방법 | |
KR100329518B1 (ko) | 무선망 설계시스템에서의 전파해석방법 | |
JP2009128025A (ja) | 測定装置、測定方法、ならびに、プログラム | |
JP5150547B2 (ja) | 置局設計支援装置及びコンピュータプログラム | |
US10651890B2 (en) | Device and method for wireless communication and communication terminal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20091207 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20111025 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20111115 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20141125 Year of fee payment: 3 |
|
LAPS | Cancellation because of no payment of annual fees |