JP2005333451A - Antenna characteristic evaluation system - Google Patents

Antenna characteristic evaluation system Download PDF

Info

Publication number
JP2005333451A
JP2005333451A JP2004150381A JP2004150381A JP2005333451A JP 2005333451 A JP2005333451 A JP 2005333451A JP 2004150381 A JP2004150381 A JP 2004150381A JP 2004150381 A JP2004150381 A JP 2004150381A JP 2005333451 A JP2005333451 A JP 2005333451A
Authority
JP
Japan
Prior art keywords
antenna
directivity
wave propagation
radio wave
evaluation system
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.)
Pending
Application number
JP2004150381A
Other languages
Japanese (ja)
Inventor
Kenji Takagi
健治 高木
Takeya Nagano
健也 長野
Tomoaki Nishikido
友昭 西木戸
Yutaka Saito
裕 斎藤
Yoshio Koyanagi
芳雄 小柳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Panasonic Holdings Corp
Original Assignee
Matsushita Electric Industrial Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Matsushita Electric Industrial Co Ltd filed Critical Matsushita Electric Industrial Co Ltd
Priority to JP2004150381A priority Critical patent/JP2005333451A/en
Publication of JP2005333451A publication Critical patent/JP2005333451A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Monitoring And Testing Of Transmission In General (AREA)
  • Radio Transmission System (AREA)

Abstract

<P>PROBLEM TO BE SOLVED: To provide an antenna characteristic evaluation system capable of performing quick and highly accurate antenna characteristic evaluation while considering complicated antenna radiation characteristics and radio wave propagation characteristics including effect on human bodies. <P>SOLUTION: When an antenna directivity analysis part 11 in the antenna characteristic evaluation system 1 analyzes the directivity of an antenna by the FDTD method, the MOMENT method or the like which is an electro-magnetic field calculation method, a storage device 13 stores antenna complex directivity 12 such as amplitude and a phase obtained by the analytical results of the antenna directivity analysis part 11. A geometrical optical radio wave propagation characteristic analysis part 14 acquires the antenna complex directivity 12 stored in the storage device 13 and uses the directivity 12 as a directivity coefficient on a transmitting/receiving point. At the point, the geometrical optical radio wave propagation characteristic analysis part 14 excludes beams of directions of which directivity coefficients are small from calculation. <P>COPYRIGHT: (C)2006,JPO&NCIPI

Description

本発明は、無線局から送信された電波の伝搬特性を計算し、計算結果のデータに基づいてアンテナ特性の評価を行うアンテナ特性評価システムに関し、特に、無線局から送信された複数の光線を探索する幾何光学的解析手法を用いてアンテナ特性の評価を行うアンテナ特性評価システムに関する。   The present invention relates to an antenna characteristic evaluation system that calculates propagation characteristics of radio waves transmitted from a radio station and evaluates antenna characteristics based on calculation result data, and in particular, searches for a plurality of light beams transmitted from a radio station. The present invention relates to an antenna characteristic evaluation system for evaluating antenna characteristics using a geometric optical analysis technique.

近年広く利用されている携帯電話機は、さまざまな使用状態や電波伝搬環境に適応して指向性や偏波を制御できるようになっている。このような制御に対応できる携帯電話機用のアンテナの開発を行うには、人体への影響やアンテナの微細構造を考慮して計算されたアンテナの全立体角複素指向性と多重波伝搬パラメータとを統合して、実際の多重波伝搬環境における実効的なアンテナ性能を評価できる系(つまり、手段)が必要となる。   Mobile phones that have been widely used in recent years can control directivity and polarization in accordance with various usage conditions and radio wave propagation environments. In order to develop an antenna for a mobile phone that can handle such control, the total solid angle complex directivity of the antenna calculated in consideration of the influence on the human body and the fine structure of the antenna and the multi-wave propagation parameter are calculated. A system (that is, means) that can integrate and evaluate effective antenna performance in an actual multiwave propagation environment is required.

このようなアンテナ性能の評価系としては従来、ランダムフィールド法などの実験的手法が用いられているが、この方法では大規模な実験装置や多大な労力が必要となるなどの課題がある。そこで、このような課題に対応するためには計算機によるアンテナ特性評価方法が広く利用されている。このようなアンテナ特性評価方法は、一般的には、無線局から送信された電磁波の電界強度分布を求めるための手段として、電磁界解析手法あるいは幾何光学的計算手法を用いたものがある。例えば、電磁波の電界強度分布を求める手段として、下記の非特許文献1に記載されているFDTD(Finite Difference Time Domain:有限差分時間領域)法を用いた室内伝搬解析法や、複数の擬似光線の軌跡を探索する幾何光学的計算手法のレイラウンチング法などがある。レイラウンチング法においては、下記の特許文献1に開示されているような、送信点から送信される全光線を計算する前処理と、受信点に到達する光線を探索して受信強度を求める後処理との2つの処理に分けて計算処理を行い、電磁波の電界強度分布を計算する方法が知られている。あるいは、下記の特許文献2に開示されているような、送信点から送信される光線を任意の空間で分割して、複数のコンピュータで電磁波の電界強度分布を計算する方法などが知られている。
特開2001−136131号公報 特開2002−232348号公報 IEICE TRANS COMMUN VOL.E86-B, NO.6 JUNE 2003
Conventionally, an experimental method such as a random field method has been used as such an antenna performance evaluation system, but this method has problems such as requiring a large-scale experimental apparatus and a great deal of labor. Therefore, in order to cope with such a problem, an antenna characteristic evaluation method using a computer is widely used. Such an antenna characteristic evaluation method generally uses an electromagnetic field analysis method or a geometric optical calculation method as a means for obtaining the electric field strength distribution of an electromagnetic wave transmitted from a radio station. For example, as means for obtaining the electric field strength distribution of electromagnetic waves, an indoor propagation analysis method using the FDTD (Finite Difference Time Domain) method described in Non-Patent Document 1 below, There is a ray launching method, which is a geometric optical calculation method for searching a locus. In the ray launching method, as disclosed in Patent Document 1 below, pre-processing for calculating all rays transmitted from a transmission point, and post-processing for searching for rays reaching the reception point to obtain reception intensity There is known a method of calculating the electric field strength distribution of electromagnetic waves by performing calculation processing in two processes. Alternatively, as disclosed in Patent Document 2 below, a method of dividing a light beam transmitted from a transmission point in an arbitrary space and calculating a field intensity distribution of electromagnetic waves by a plurality of computers is known. .
JP 2001-136131 A JP 2002-232348 A IEICE TRANS COMMUN VOL.E86-B, NO.6 JUNE 2003

しかしながら、FDTD法を用いた電波伝搬解析方法では、一般的に、電波の周波数が高い、又は電界強度分布の解析領域が大きいと演算量が多くなり、屋外伝搬モデルなどの大規模解析領域において電波の伝搬解析を行うことが困難となる。また、幾何光学的計算手法を用いた電波伝搬特性の解析では、送受信点におけるアンテナは人体への影響を含んだアンテナ特性についての考慮がなされておらず、携帯電話機用アンテナの開発で重要な人体への影響を含んだ電波伝搬特性の解析を行うことが不可能である。さらに、電磁波の電界強度分布を計算するために複数の計算機を用いた計算高速化の手法では、設備投資に多大なコストがかかる等の課題がある。   However, in the radio wave propagation analysis method using the FDTD method, in general, when the frequency of radio waves is high or the analysis area of the electric field strength distribution is large, the amount of calculation increases. It becomes difficult to perform propagation analysis. In addition, in the analysis of radio wave propagation characteristics using geometric optical calculation methods, the antenna at the transmission / reception point does not consider the antenna characteristics including the effects on the human body, which is important for the development of mobile phone antennas. It is impossible to analyze radio wave propagation characteristics including effects on Furthermore, the calculation speed-up method using a plurality of computers for calculating the electric field strength distribution of the electromagnetic wave has a problem that the capital investment is very expensive.

本発明は、上記の事情に鑑みてなされたものであり、計算機によって人体などへの影響を含んだアンテナ放射特性と多重波伝搬環境とを統合して高速かつ高精度にアンテナ特性の評価を行うことが可能なアンテナ特性評価システムを提供することを目的とする。   The present invention has been made in view of the above circumstances, and evaluates antenna characteristics at high speed and with high accuracy by integrating the antenna radiation characteristics including effects on the human body and the multi-wave propagation environment by a computer. An object of the present invention is to provide an antenna characteristic evaluation system.

上記の課題を解決するために、本発明のアンテナ特性評価システムは、無線局から送信された電波の伝搬特性を計算してアンテナ特性の評価を行うアンテナ特性評価システムであって、第1アンテナ及び第2アンテナの離散的角度方位に対する複素指向性係数列を算出するアンテナ指向性解析手段と、送信点より複数の光線を放射し、その光線が受信点に到達するまでの軌跡を追跡することで所定の電波伝搬モデルにおける送受信点間の電波の伝搬特性を算出する幾何光学的電波伝搬特性解析手段とを備え、幾何光学的電波伝搬特性解析手段が、第1アンテナ及び第2アンテナの複素指向性係数列を送受信点におけるアンテナの指向性係数として用いることにより、アンテナ指向性解析手段と幾何光学的電波伝搬特性解析手段とを統合するように構成したことを特徴とする。   In order to solve the above problems, an antenna characteristic evaluation system of the present invention is an antenna characteristic evaluation system that calculates a propagation characteristic of a radio wave transmitted from a radio station and evaluates the antenna characteristic. By means of antenna directivity analysis means for calculating a complex directivity coefficient sequence for the discrete angle orientation of the second antenna, and by radiating a plurality of light rays from the transmission point and tracking the trajectory until the light rays reach the reception point Geometric optical wave propagation characteristic analyzing means for calculating propagation characteristics of radio waves between transmission and reception points in a predetermined radio wave propagation model, and the geometric optical wave propagation characteristic analyzing means is a complex directivity of the first antenna and the second antenna. By using the coefficient sequence as the antenna directivity coefficient at the transmission / reception point, the antenna directivity analysis means and the geometric optical wave propagation characteristic analysis means are integrated. Characterized by being configured to.

また、本発明のアンテナ特性評価システムは、さらに、送信点におけるアンテナの複素指向性係数列の中で、離散的角度方位に対する振幅が所定の閾値を越えるか否かを判定する振幅判定手段を備え、振幅判定手段の判定結果に応じて、送信点から光線を放射するか否かを決定するように構成したことを特徴とする。   In addition, the antenna characteristic evaluation system of the present invention further includes amplitude determination means for determining whether or not the amplitude with respect to the discrete angle azimuth exceeds a predetermined threshold in the complex directivity coefficient sequence of the antenna at the transmission point. According to the determination result of the amplitude determining means, it is configured to determine whether or not to emit a light beam from the transmission point.

また、本発明のアンテナ特性評価システムにおいては、アンテナ指向性解析手段は有限差分時間領域法を用いた構成としたことを特徴とする。また、本発明のアンテナ特性評価システムは、さらに、第1アンテナ及び第2アンテナの離散的角度方位に対する複素指向性係数列を測定するアンテナ指向性測定手段を備え、幾何光学的電波伝搬特性解析手段が、アンテナ指向性測定手段によって得られた複素指向性係数列を送受信点におけるアンテナの指向性係数として用いるように構成したことを特徴とする。   In the antenna characteristic evaluation system of the present invention, the antenna directivity analysis means is configured to use a finite difference time domain method. The antenna characteristic evaluation system according to the present invention further includes an antenna directivity measuring means for measuring a complex directivity coefficient sequence with respect to the discrete angle azimuths of the first antenna and the second antenna, and a geometric optical wave propagation characteristic analyzing means. However, the complex directivity coefficient sequence obtained by the antenna directivity measuring means is used as the antenna directivity coefficient at the transmission / reception point.

また、本発明のアンテナ特性評価システムは、幾何光学的電波伝搬特性解析手段が保有する複数の電波伝搬モデルを記憶する記憶手段と、第1アンテナ及び第2アンテナの複素指向性係数列を用いて、あらかじめ指定した複数の電波伝搬モデルに対して幾何光学的電波伝搬特性解析を実行する実行手段とを備える構成としたことを特徴とする。   Further, the antenna characteristic evaluation system of the present invention uses a storage means for storing a plurality of radio wave propagation models held by the geometric optical radio wave propagation characteristic analysis means, and a complex directivity coefficient sequence of the first antenna and the second antenna. The present invention is characterized by comprising an execution means for executing geometric optical wave propagation characteristic analysis on a plurality of radio wave propagation models specified in advance.

また、本発明のアンテナ特性評価システムにおいては、幾何光学的電波伝搬特性解析手段は、記憶手段にあらかじめ記憶された電波伝搬モデルに対して幾何光学的電波伝搬特性の解析計算を実行し、アンテナ指向性解析手段又はアンテナ指向性測定手段によって得られたアンテナの複素指向性係数列を用いて位相合成の計算を行うように構成したことを特徴とする。   Further, in the antenna characteristic evaluation system of the present invention, the geometric optical radio wave propagation characteristic analyzing unit executes an analysis calculation of the geometric optical radio wave propagation characteristic with respect to the radio wave propagation model stored in the storage unit in advance, and The phase composition calculation is performed using the complex directivity coefficient sequence of the antenna obtained by the characteristic analysis means or the antenna directivity measurement means.

本発明のアンテナ特性評価システムによれば、アンテナ指向性解析手段と幾何光学的電波伝搬特性解析手段とを統合することにより、実際の多重波伝搬環境における高精度なアンテナ特性評価を比較的短時間で行うことができる。また、幾何光学的電波伝搬特性解析手段において送信点から送信される光線の本数を軽減することができるので、アンテナ特性の計算速度が一段と向上する。さらに、人体への影響を含んだアンテナ特性と実際の多重波伝搬環境における電波の伝搬特性とを考慮したアンテナ特性評価を行うことが可能となる。   According to the antenna characteristic evaluation system of the present invention, highly accurate antenna characteristic evaluation in an actual multiwave propagation environment can be performed in a relatively short time by integrating the antenna directivity analysis means and the geometric optical wave propagation characteristic analysis means. Can be done. In addition, since the geometric optical wave propagation characteristic analyzing means can reduce the number of light rays transmitted from the transmission point, the calculation speed of the antenna characteristic is further improved. Furthermore, it is possible to perform antenna characteristic evaluation in consideration of antenna characteristics including effects on the human body and radio wave propagation characteristics in an actual multiwave propagation environment.

また、本発明のアンテナ特性評価システムによれば、アンテナ指向性解析手段に有限差分時間領域法を用いているので、人体への影響を含んだアンテナ指向性解析を比較的高精度に計算機で行うことが可能となり、アンテナ指向性解析と幾何光学的電波伝搬特性解析を同時に考慮したアンテナ特性評価を比較的短時間に行うことができる。さらに、アンテナ指向性測定手段によって得られた複素指向性係数列も送受信点におけるアンテナの指向性係数として用いることが可能であるため、実際に測定したアンテナ特性を用いた電波伝搬計算が行うことができ、実際の多重波伝搬環境におけるアンテナ特性評価を行うことができる。   Further, according to the antenna characteristic evaluation system of the present invention, the antenna directivity analysis means uses the finite difference time domain method, so that the antenna directivity analysis including the influence on the human body is performed with a computer with relatively high accuracy. Therefore, it is possible to perform antenna characteristic evaluation in consideration of antenna directivity analysis and geometric optical wave propagation characteristic analysis in a relatively short time. Furthermore, since the complex directivity coefficient sequence obtained by the antenna directivity measurement means can also be used as the antenna directivity coefficient at the transmission / reception point, radio wave propagation calculation using actually measured antenna characteristics can be performed. It is possible to evaluate antenna characteristics in an actual multiwave propagation environment.

また、本発明のアンテナ特性評価システムによれば、あらかじめ指定した複数の電波伝搬モデルに対して幾何光学的電波伝搬特性の解析を行って振幅や位相差などのアンテナ特性を計算することにより、市街地モデルや郊外モデルなどのあらゆる伝搬モデルにおけるアンテナ特性を評価することができる。さらに、本発明のアンテナ特性評価システムによれば、幾何光学的電波伝搬特性の解析を行うときにアンテナの複素指向性係数列を用いて位相合成計算を行うので、オペレータが幾何光学的電波伝搬解析計算を改めて行う必要がなくなる。したがって、幾何光学的電波伝搬解析計算で得られた伝搬パラメータと複素指向性係数列を用いて、比較的短時間で実際の電波伝搬環境におけるアンテナ特性の評価を行うことができる。   In addition, according to the antenna characteristic evaluation system of the present invention, by analyzing geometric optical wave propagation characteristics for a plurality of radio wave propagation models specified in advance and calculating antenna characteristics such as amplitude and phase difference, It is possible to evaluate antenna characteristics in all propagation models such as models and suburban models. Furthermore, according to the antenna characteristic evaluation system of the present invention, the phase synthesis calculation is performed using the complex directivity coefficient sequence of the antenna when the geometric optical radio wave propagation characteristic is analyzed, so that the operator can perform the geometric optical radio wave propagation analysis. There is no need to perform calculations again. Therefore, the antenna characteristics in the actual radio wave propagation environment can be evaluated in a relatively short time by using the propagation parameter and the complex directivity coefficient sequence obtained by the geometric optical wave propagation analysis calculation.

本発明にかかるアンテナ特性評価システムの骨子は、アンテナ指向性解析部により得られたアンテナ複素指向性を幾何光学的電波伝搬特性解析部の送受信地点における指向性係数として用いたことにある。これにより、実際の多重波伝搬環境における高精度なアンテナ特性評価を比較的短時間で行うことができる。したがって、計算機による高精度なアンテナ特性評価システムとして有効に利用することが可能となる。   The gist of the antenna characteristic evaluation system according to the present invention is that the antenna complex directivity obtained by the antenna directivity analysis unit is used as the directivity coefficient at the transmission / reception point of the geometric optical wave propagation characteristic analysis unit. Thereby, highly accurate antenna characteristic evaluation in an actual multiwave propagation environment can be performed in a relatively short time. Therefore, it can be effectively used as a highly accurate antenna characteristic evaluation system by a computer.

以下、図面を参照しながら、本発明におけるアンテナ特性評価システムの実施の形態の幾つかを詳細に説明する。なお、以下の説明では、各図面において共通する構成要素は同一符号を付し、かつ重複する説明は省略する。   Hereinafter, some embodiments of an antenna characteristic evaluation system according to the present invention will be described in detail with reference to the drawings. In the following description, components common to the drawings are given the same reference numerals, and redundant descriptions are omitted.

<実施の形態1>
本発明におけるアンテナ特性評価システムの実施の形態1について図1から図5を用いて説明する。図1は、本発明の実施の形態1のアンテナ特性評価システムにおけるアンテナ指向性解析部と幾何光学的電波伝搬特性解析部の統合を示すブロック図である。図1において、アンテナ特性評価システム1は、アンテナ複素指向性12の解析を行うアンテナ指向性解析部(アンテナ指向性解析手段)11と、アンテナ指向性解析部11が解析したアンテナ複素指向性12を記憶する記憶装置(記憶手段)13と、記憶装置13に記憶されているアンテナ複素指向性12を取得して指向性係数として用いる幾何光学的電波伝搬特性解析部(幾何光学的電波伝搬特性解析手段)14とを備えた構成となっている。
<Embodiment 1>
Embodiment 1 of an antenna characteristic evaluation system according to the present invention will be described with reference to FIGS. FIG. 1 is a block diagram showing integration of an antenna directivity analysis unit and a geometric optical wave propagation characteristic analysis unit in the antenna characteristic evaluation system according to Embodiment 1 of the present invention. In FIG. 1, an antenna characteristic evaluation system 1 includes an antenna directivity analysis unit (antenna directivity analysis means) 11 that analyzes an antenna complex directivity 12, and an antenna complex directivity 12 analyzed by the antenna directivity analysis unit 11. A storage device (storage means) 13 for storing, and a geometric optical radio wave propagation characteristic analysis unit (geometric optical wave propagation characteristic analysis means) that acquires the antenna complex directivity 12 stored in the storage device 13 and uses it as a directivity coefficient ) 14.

次に、図1に示すアンテナ特性評価システム1の動作の概要を説明する。アンテナ特性評価システム1は、アンテナ指向性解析部11が、電磁界計算法であるFDTD法やMOMENT法などによってアンテナ指向性の解析を行うと、記憶装置13が、アンテナ指向性解析部11の解析結果によって得られた振幅や位相などのアンテナ複素指向性12を記憶する。さらに、幾何光学的電波伝搬特性解析部14が、記憶装置13に記憶されているアンテナ複素指向性12を取得して送受信点における指向性係数として用いる。このとき、幾何光学的電波伝搬特性解析部14は指向性係数が小さい方向の光線を計算から除外する。これによって、短時間かつ簡単な計算処理によってアンテナ特性評価を行うことができる。   Next, an outline of the operation of the antenna characteristic evaluation system 1 shown in FIG. 1 will be described. In the antenna characteristic evaluation system 1, when the antenna directivity analysis unit 11 analyzes the antenna directivity by the FDTD method or the MOMENT method that are electromagnetic field calculation methods, the storage device 13 analyzes the antenna directivity analysis unit 11. The antenna complex directivity 12 such as amplitude and phase obtained as a result is stored. Further, the geometric optical wave propagation characteristic analysis unit 14 acquires the antenna complex directivity 12 stored in the storage device 13 and uses it as the directivity coefficient at the transmission / reception point. At this time, the geometric optical wave propagation characteristic analysis unit 14 excludes light rays in a direction having a small directivity coefficient from the calculation. As a result, the antenna characteristics can be evaluated by a simple calculation process in a short time.

次に、フローチャートを用いて、図1に示すアンテナ特性評価システム1の動作を詳細説明する。図2は、図1に示す幾何光学的電波伝搬特性解析部が行う電波伝搬特性解析の前処理の流れを示すフローチャートである。また、図3は、図1に示す幾何光学的電波伝搬特性解析部が行う電波伝搬特性解析の後処理の流れを示すフローチャートである。つまり、図2及び図3は、図1のアンテナ特性評価システム1における幾何光学的電波伝搬特性解析部14が行う処理のフローチャートを示しており、図2に示すような送信点から送信される全光線を計算する処理を前処理と定義し、図3に示すような受信点に到達する光線を探索して受信強度を求める処理を後処理と定義することにする。   Next, the operation of the antenna characteristic evaluation system 1 shown in FIG. 1 will be described in detail using a flowchart. FIG. 2 is a flowchart showing a pre-processing flow of radio wave propagation characteristic analysis performed by the geometric optical wave propagation characteristic analysis unit shown in FIG. FIG. 3 is a flowchart showing a post-processing flow of radio wave propagation characteristic analysis performed by the geometric optical wave propagation characteristic analysis unit shown in FIG. That is, FIG. 2 and FIG. 3 show a flowchart of the processing performed by the geometric optical wave propagation characteristic analysis unit 14 in the antenna characteristic evaluation system 1 of FIG. 1, and all the data transmitted from the transmission point as shown in FIG. The processing for calculating the light rays is defined as preprocessing, and the processing for searching the light rays reaching the reception point as shown in FIG.

また、図4は、図2のフローチャートで示される前処理によって記憶装置に記憶されるデータ構成を示す図である。図4に示すように、電波伝搬特性解析の前処理によって、記憶装置13には、全光線の方程式41、全光線の進行角度42、全光線の偏波方程式43、全光線の反射面番号44、及び全光線の壁面番号45が記憶される。   FIG. 4 is a diagram showing a data configuration stored in the storage device by the preprocessing shown in the flowchart of FIG. As shown in FIG. 4, the pre-processing of the radio wave propagation characteristic analysis causes the storage device 13 to store the total light equation 41, the total light traveling angle 42, the total light polarization equation 43, and the total light reflection surface number 44. , And the wall number 45 of all rays.

また、図5は、本発明の第1の実施の形態における電波伝搬特性解析の計算を行うためのコンピュータシステムの構成図である。つまり、電波の電磁界計算及び図2や図3に示すような電波伝搬特性の解析処理を行うためのハードウェア構成は、例えば、図5に示すような一般的なコンピュータシステムを用いて実現することができる。   FIG. 5 is a configuration diagram of a computer system for performing calculation of radio wave propagation characteristic analysis in the first embodiment of the present invention. That is, the hardware configuration for performing the electromagnetic field calculation of radio waves and the analysis processing of the radio wave propagation characteristics as shown in FIG. 2 and FIG. 3 is realized using, for example, a general computer system as shown in FIG. be able to.

図5に示すように、このコンピュータシステムは、電波伝搬特性の演算処理を行うパーソナルコンピュータ(以下、パソコンと云う)などの制御装置51と、制御装置51の情報を取得して表示を行うCRTなどの表示装置52と、制御装置51の情報を取得して印刷を行うプリンタなどの印刷装置53と、データの入力を行うキーボードやマウスなどの入力装置54と、電波伝搬特性の解析結果を記憶するハードディスクなどの記憶装置13とを備えた構成となっている。なお、制御装置51は、請求項2に記載されている、送信点におけるアンテナの複素指向性係数列の中で離散的角度方位に対する振幅が所定の閾値を越えるか否かを判定する振幅判定手段と、請求項5に記載されている、アンテナの複素指向性係数列を用いてあらかじめ指定した複数の電波伝搬モデルに対して幾何光学的電波伝搬特性解析を実行する実行手段とを備えている。また、図5のコンピュータシステムは、脱着可能なCDやフロッピーディスク(登録商標)などの記憶媒体55を使用することができる。さらに、記憶装置13及び記憶媒体55は、電磁界計算及び図2や図3に示す電波伝搬特性の解析処理を実行するためプログラムやパラメータなどを記憶することができる。   As shown in FIG. 5, the computer system includes a control device 51 such as a personal computer (hereinafter referred to as a personal computer) that performs calculation processing of radio wave propagation characteristics, a CRT that acquires information from the control device 51 and displays the information. Display device 52, printing device 53 such as a printer that obtains and prints information from control device 51, input device 54 such as a keyboard or mouse for inputting data, and analysis results of radio wave propagation characteristics are stored. The configuration includes a storage device 13 such as a hard disk. Note that the control device 51 determines whether or not the amplitude with respect to the discrete angle azimuth exceeds a predetermined threshold in the complex directivity coefficient sequence of the antenna at the transmission point according to claim 2. And execution means for executing geometric optical wave propagation characteristic analysis on a plurality of wave propagation models specified in advance using the complex directivity coefficient sequence of the antenna. Further, the computer system of FIG. 5 can use a storage medium 55 such as a removable CD or a floppy disk (registered trademark). Furthermore, the storage device 13 and the storage medium 55 can store programs, parameters, and the like for executing electromagnetic field calculations and radio wave propagation characteristic analysis processing shown in FIGS.

以下、図2及び図3のフローチャートを用いて、本発明の実施の形態1における電波伝搬特性解析の計算手順を説明する。なお、必要に応じて図1から図5の各図面を参照しながらフローチャートを説明する。   Hereinafter, the calculation procedure of the radio wave propagation characteristic analysis in the first embodiment of the present invention will be described using the flowcharts of FIGS. 2 and 3. The flowchart will be described with reference to FIGS. 1 to 5 as necessary.

最初に、図2に示す電波伝搬特性解析の前処理の流れを説明する。まず、パラメータ入力の手順において、オペレータは、光線の最大反射回数Rmax、閾値、周波数、壁面材質、放射光線数、及び送信位置などの各種のパラメータを入力装置54から入力する(ステップST1)。また、壁面方程式導出の手順において、制御装置51は、入力されたパラメータに対応する壁面の形状に基づいて壁面方程式を導き出し、この壁面方程式を記憶装置13に記憶させる(ステップST2)。さらに、制御装置51は、複素指向性の各々の方位における振幅導出の手順において、ステップST1で入力されたパラメータの放射光線数Nにより求めたアンテナ複素指向性12の各々の角度ごとに電波の振幅を求める(ステップST3)。   First, the flow of preprocessing for radio wave propagation characteristic analysis shown in FIG. 2 will be described. First, in the parameter input procedure, the operator inputs various parameters such as the maximum number of reflections Rmax of light rays, a threshold value, a frequency, a wall surface material, the number of radiation rays, and a transmission position from the input device 54 (step ST1). In the procedure for deriving the wall surface equation, the control device 51 derives the wall surface equation based on the shape of the wall surface corresponding to the input parameter, and stores this wall surface equation in the storage device 13 (step ST2). Further, the control device 51 determines the amplitude of the radio wave for each angle of the antenna complex directivity 12 obtained from the number of radiated rays N of the parameter input in step ST1 in the procedure for deriving the amplitude in each direction of the complex directivity. Is obtained (step ST3).

次に、制御装置51は、電波の振幅が所定の閾値以上であるか否かを判定する(ステップST4)。ここで、制御装置51は、電波の振幅が所定の閾値以上ならば(ステップST4でYESの場合)、放射光線方程式を導き出す(ステップST5)。一方、電波の振幅が所定の閾値未満ならば(ステップST4でNOの場合)、放射光線数Nの光線番号nを一つ加算してn+1とし(ステップST6)、ステップST3に戻ってアンテナ複素指向性12の各々の方位における振幅の導出を継続する。このようにして、入力された各方位のアンテナ複素指向性12を用いて放射する光線を減少させることにより、電波伝搬特性解析の計算時間を短縮することができる。   Next, the control device 51 determines whether or not the amplitude of the radio wave is greater than or equal to a predetermined threshold (step ST4). Here, if the amplitude of the radio wave is equal to or greater than a predetermined threshold value (in the case of YES in step ST4), the control device 51 derives a radiation ray equation (step ST5). On the other hand, if the amplitude of the radio wave is less than the predetermined threshold value (in the case of NO in step ST4), one ray number n of the number of radiated rays N is added to be n + 1 (step ST6), and the process returns to step ST3 and the antenna complex pointing Continue to derive the amplitude in each orientation of gender 12. In this way, the calculation time of the radio wave propagation characteristic analysis can be shortened by reducing the rays emitted using the antenna complex directivity 12 of each azimuth input.

なお、ステップST5における放射光線方程式の導出処理では、送信点を中心としてステップST1のパラメータ入力での放射光線数Nにより求めた所定の離散角度ごとに全方位に放射する光線の方程式を求めたが、光線番号nが光線数Nより小さければ、ステップST7の手順のように光線番号nを1つ加算して(つまり、光線番号をn+1にして)、ステップST3の振幅導出手順に戻り、前述の各ステップを繰り返して放射光線方程式の導出処理を行う。   In the radiation beam equation derivation process in step ST5, an equation of light rays radiated in all directions is obtained for each predetermined discrete angle obtained from the number N of radiation rays at the parameter input in step ST1 with the transmission point as the center. If the ray number n is smaller than the ray number N, the ray number n is incremented by 1 as in the procedure of step ST7 (that is, the ray number is set to n + 1), and the procedure returns to the amplitude derivation procedure of step ST3. Each step is repeated to derive the radiation ray equation.

次に、ステップST5における放射光線方程式の導出ステップで求めた放射光線方程式に対して、振幅や位相などのアンテナ複素指向性12を乗算する(ステップST8)。そして、反射光線方程式導出の手順において、ある特定の光線の反射壁面を決定し、その反射壁面によって所望の反射光線方程式を導き出す(ステップST9)。このような反射光線方程式の導出処理を行う度に反射回数Rを1ずつ加算して行き(つまり、反射光線方程式の導出処理を行うごとに反射回数をR+1として行き)、最大反射回数Rmaxとなるまで反射光線方程式の導出処理を繰り返す(ステップST10)。又は、全放射光線Nαにおける反射光線方程式を求めるまで光線番号nを1ずつ加算して行き(つまり、反射光線方程式の導出処理を行うごとに光線番号をn+1として行き)、最大放射光線数Nmaxとなるまで導出処理を繰り返す(ステップST11)。   Next, the radiation ray equation obtained in the step of deriving the radiation ray equation in step ST5 is multiplied by the antenna complex directivity 12 such as amplitude and phase (step ST8). Then, in the procedure for deriving the reflected ray equation, a reflection wall surface of a specific ray is determined, and a desired reflected ray equation is derived from the reflection wall surface (step ST9). Each time such a reflected ray equation derivation process is performed, the number of reflections R is incremented by 1 (that is, each time the reflected ray equation is derived, the number of reflections is R + 1), and the maximum number of reflections Rmax is obtained. The process of deriving the reflected ray equation is repeated until step ST10. Alternatively, the ray number n is incremented by 1 until the reflected ray equation for the total emitted ray Nα is obtained (that is, the ray number is set to n + 1 every time the reflected ray equation is derived), and the maximum number of emitted rays Nmax The derivation process is repeated until it becomes (step ST11).

図2に示すような電波伝搬特性解析の前処理によって、記憶装置13には、ステップST1のパラメータ入力の手順において入力されたデータと、図4に示すような全光線の方程式41、全光線の進行角度42、全光線の偏波方程式43、全光線の反射面番号44、及び全光線の壁面番号45が記憶される。なお、全放射光線数Nαは、ステップST1で送信点においてパラメータ入力された放射光線本数Nから、ステップST4で振幅が閾値より大きいときに除外された光線の本数を差し引いたものである。   By the pre-processing of the radio wave propagation characteristic analysis as shown in FIG. 2, the data input in the parameter input procedure of step ST1 and the total ray equation 41 as shown in FIG. The traveling angle 42, the polarization equation 43 of all rays, the reflection surface number 44 of all rays, and the wall number 45 of all rays are stored. Note that the total number of emitted rays Nα is obtained by subtracting the number of rays excluded when the amplitude is larger than the threshold value in step ST4 from the number N of emitted rays input as parameters at the transmission point in step ST1.

次に、図3に示す電波伝搬特性解析の後処理の流れを説明する。まず、パラメータ入力の手順において、オペレータは、図2に示す前処理においてアンテナ解析により得られたアンテナ複素指向性12と受信位置及び電波の周波数をパラメータとして入力装置54から入力する(ステップST21)。そして、制御装置51は、受信球導出の手順において、記憶装置13から読み込んだ光線番号aの受信地点までの距離dを導き出し、送信角度間隔をφとしたとき、後処理における受信球の半径rを次の式(1)によって求める(ステップST22)。
r=21/2d・sin(φ/2) (1)
Next, the post-processing flow of radio wave propagation characteristic analysis shown in FIG. 3 will be described. First, in the parameter input procedure, the operator inputs the antenna complex directivity 12 obtained by the antenna analysis in the preprocessing shown in FIG. 2, the reception position, and the frequency of the radio wave from the input device 54 as parameters (step ST21). Then, the control device 51 derives the distance d to the reception point of the ray number a read from the storage device 13 in the procedure for deriving the reception sphere, and when the transmission angle interval is φ, the radius r of the reception sphere in post-processing. Is obtained by the following equation (1) (step ST22).
r = 2 1/2 d · sin (φ / 2) (1)

次に、制御装置51は、反射光線方程式と受信点を中心とする半径rの球との交点を求めることにより到達光線の探索を行い、光線が受信点に到達したか否かを判定する(ステップST23)。ここで、光線が受信点に到達した場合は(ステップST23でYESの場合)、受信点に到達した光線に光線到来角度に対応したアンテナ複素指向性12を乗算する(ステップST24)。一方、光線番号aの光線が受信点に到達していないときは(ステップST23でNOの場合)、光線番号aに1を加算してa+1にして(ステップST25)、ステップST22の受信球導出のステップに戻り、次の光線番号(a+1)の受信球導出処理に移行する。   Next, the control device 51 searches for the reaching ray by obtaining the intersection of the reflected ray equation and a sphere having a radius r centered on the receiving point, and determines whether or not the ray has reached the receiving point ( Step ST23). Here, when the light beam reaches the reception point (YES in step ST23), the light beam reaching the reception point is multiplied by the antenna complex directivity 12 corresponding to the light beam arrival angle (step ST24). On the other hand, when the ray with the ray number a has not reached the reception point (in the case of NO at step ST23), 1 is added to the ray number a to become a + 1 (step ST25), and the received sphere is derived at step ST22. Returning to the step, the process proceeds to the reception sphere derivation process for the next ray number (a + 1).

また、ステップST24で、受信点に到達した光線に光線到来角度に対応したアンテナ複素指向性12を乗算した場合は、受信点に到達する経路は光線ごとに異なるため、それぞれの光線には位相差や振幅差が生じる。したがって、これらの経路が異なる光線の振幅について位相差を考慮して加算し位相合成を行う(ステップST26)。そして、全ての到来角度の光線について位相合成を行うまでステップST26からステップST25の帰還ループを辿り、ステップST22からステップST26までの処理を繰り返す。   Further, in step ST24, when the light beam that has reached the reception point is multiplied by the antenna complex directivity 12 corresponding to the light beam arrival angle, the path that reaches the reception point differs for each light beam. And amplitude difference occurs. Therefore, the phase synthesis is performed by adding the amplitudes of the light beams having different paths in consideration of the phase difference (step ST26). Then, the process from step ST22 to step ST26 is repeated by following the feedback loop from step ST26 to step ST25 until the phase synthesis is performed for the rays of all the arrival angles.

次に、電界強度導出の手順において位相合成された光線の電界強度を求め、この電界強度を記憶装置13に記憶させる(ステップST27)。このようにして、ある受信地点において位相合成された光線の電界強度を求めたら、受信地点mを追尾しながら(つまり、受信地点をm+1にして)、オペレータが設定したP個の受信地点を移動し(ステップS28)、ステップST22の受信球の導出処理に戻って前述と同様のステップで処理を繰り返す。そして、結果出力の手順において、各々の受信地点における受信電界強度を結果出力として記憶装置13に記憶させる(ステップST29)。   Next, the electric field strength of the light beam subjected to phase synthesis in the procedure for deriving the electric field strength is obtained, and this electric field strength is stored in the storage device 13 (step ST27). In this way, when the electric field strength of the phase-combined light beam at a certain reception point is obtained, the reception point m is tracked (that is, the reception point is set to m + 1), and the P reception points set by the operator are moved. (Step S28), the process returns to the reception sphere derivation process of Step ST22, and the process is repeated in the same steps as described above. In the result output procedure, the received electric field intensity at each reception point is stored in the storage device 13 as a result output (step ST29).

以上説明したように、本発明の実施の形態1のアンテナ特性評価システムによれば、アンテナ指向性解析部11により得られたアンテナ複素指向性12を送信点及び受信点において指向性係数として使用し、送信点において振幅の小さい方位の光線を放射しない方法を用いることにより、人体などへの影響を含んだ複雑な特性を持つアンテナの実フィールドにおける高精度な電波伝搬特性の解析を行うことができる。さらに、本発明の実施の形態1のアンテナ特性評価システムによれば、送信点において用いたアンテナ複素指向性12により放射光線数を減少させることが可能となる。例えば、八木アンテナ(登録商標)などの単一指向性のアンテナ特性を用いることにより、光線の放射数を40%程度まで減少させることができる。   As described above, according to the antenna characteristic evaluation system of Embodiment 1 of the present invention, the antenna complex directivity 12 obtained by the antenna directivity analysis unit 11 is used as the directivity coefficient at the transmission point and the reception point. By using a method that does not radiate light beams with a small amplitude at the transmission point, it is possible to analyze radio wave propagation characteristics with high accuracy in the real field of antennas that have complex characteristics including effects on the human body. . Furthermore, according to the antenna characteristic evaluation system of Embodiment 1 of the present invention, the number of radiation rays can be reduced by the antenna complex directivity 12 used at the transmission point. For example, by using unidirectional antenna characteristics such as Yagi Antenna (registered trademark), the number of rays emitted can be reduced to about 40%.

このようにして、アンテナ複素指向性を求めるアンテナ指向性解析部11が電磁界計算法であるFDTD法を用いることにより、人体などへの影響を含んだ複雑な特性を持つアンテナ複素指向性12をより正確に計算することができる。なお、本発明におけるアンテナ特性評価システム1は、FDTD法やMOMENT法などのアンテナ指向性解析部11と幾何光学的電波伝搬特性解析部14とを同一のパソコンシステムに備える構成としたが、本発明はこのようなシステムに限定されるものではない。例えば、アンテナ指向性解析部11と幾何光学的電波伝搬特性解析部14がそれぞれ別のパソコンシステムにあって、アンテナ複素指向性12を記憶装置13や脱着可能な記憶媒体55を媒介として受け渡す構成としても本発明のアンテナ特性評価システムを実現することができる。   In this way, the antenna directivity analysis unit 11 for obtaining the antenna complex directivity uses the FDTD method, which is an electromagnetic field calculation method, so that the antenna complex directivity 12 having complicated characteristics including the influence on the human body and the like can be obtained. It can be calculated more accurately. The antenna characteristic evaluation system 1 according to the present invention is configured to include the antenna directivity analysis unit 11 such as the FDTD method and the MOMENT method and the geometric optical wave propagation characteristic analysis unit 14 in the same personal computer system. Is not limited to such a system. For example, the antenna directivity analysis unit 11 and the geometric optical wave propagation characteristic analysis unit 14 are provided in different personal computer systems, and the antenna complex directivity 12 is transferred via the storage device 13 or the removable storage medium 55 as a medium. However, the antenna characteristic evaluation system of the present invention can be realized.

<実施の形態2>
次に、図6から図10を用いて、本発明の実施の形態2におけるアンテナ特性評価システムについて説明する。図6は、本発明の実施の形態2における電波伝搬特性解析の計算を行うためのコンピュータシステムの構成図である。つまり、この図は、本発明の実施の形態2におけるアンテナ特性解析装置と計算機の構成を示している。図6に示す実施の形態2のアンテナ特性評価システムにおけるコンピュータシステムの構成は、図5に示す実施の形態1の構成に対して、ネットワークアナライザなどのアンテナ特性測定装置61と被測定アンテナ62及び受信アンテナ(又は送信アンテナ)63が追加させたものであり、一般的なパソコンシステムによって構成することができる。上記以外の構成要素である制御装置51、表示装置52、印刷装置53、入力装置54、記憶装置13、及び記憶媒体55は図5に示す実施の形態1と同様であるのでその説明は省略する。なお、記憶装置13及び記憶媒体55は、電磁界計算及び電波伝搬解析処理を実行するためプログラムやパラメータなどを記憶することができると共に、アンテナ特性測定装置61で得られた振幅や位相などのアンテナ複素指向性12を記憶することができる。
<Embodiment 2>
Next, the antenna characteristic evaluation system according to Embodiment 2 of the present invention will be described with reference to FIGS. FIG. 6 is a configuration diagram of a computer system for performing a radio wave propagation characteristic analysis calculation in the second embodiment of the present invention. That is, this figure shows the configuration of the antenna characteristic analysis apparatus and the computer according to Embodiment 2 of the present invention. The configuration of the computer system in the antenna characteristic evaluation system of the second embodiment shown in FIG. 6 is different from that of the first embodiment shown in FIG. 5 in the antenna characteristic measuring device 61 such as a network analyzer, the antenna to be measured 62 and the reception. An antenna (or transmission antenna) 63 is added, and can be configured by a general personal computer system. The control device 51, the display device 52, the printing device 53, the input device 54, the storage device 13, and the storage medium 55, which are constituent elements other than those described above, are the same as those in the first embodiment shown in FIG. . The storage device 13 and the storage medium 55 can store programs, parameters, and the like for executing electromagnetic field calculation and radio wave propagation analysis processing, and antennas such as amplitude and phase obtained by the antenna characteristic measurement device 61. The complex directivity 12 can be stored.

また、図7は、図6の記憶装置が行う電波伝搬特性の解析処理を機能的に示したブロック図である。つまり、この図は、本発明の実施の形態2においてアンテナ特性解析装置からの位相や振幅などの情報を電波伝搬特性の解析処理のために受け渡す状態を示す機能ブロック図である。図7に示す電波伝搬特性の解析処理の機能ブロックは、アンテナ複素指向性12及び複数の電波伝搬モデル71を記憶する記憶装置13と、被測定アンテナ62及び受信アンテナ(又は送信アンテナ)63を備えるアンテナ特性測定装置(アンテナ指向性測定手段)61と、電波伝搬特性解析の前処理30及び後処理31を実行する幾何光学的電波伝搬特性解析部14とを備えた構成となっている。   FIG. 7 is a block diagram functionally showing the radio wave propagation characteristic analysis process performed by the storage device of FIG. That is, this figure is a functional block diagram showing a state in which information such as phase and amplitude from the antenna characteristic analysis apparatus is delivered for radio wave propagation characteristic analysis processing in Embodiment 2 of the present invention. The functional block of the radio wave propagation characteristic analysis processing shown in FIG. 7 includes a storage device 13 that stores the antenna complex directivity 12 and a plurality of radio wave propagation models 71, a measured antenna 62, and a reception antenna (or transmission antenna) 63. The antenna characteristic measuring device (antenna directivity measuring means) 61 and a geometric optical wave propagation characteristic analyzing unit 14 that executes pre-processing 30 and post-processing 31 of radio wave propagation characteristic analysis are provided.

次に、図6及び図7を用いて実施の形態2におけるアンテナ特性評価システムの動作を説明する。まず、制御装置51は、ネットワークアナライザなどのアンテナ特性測定装置61を用いて電波の振幅や位相などのアンテナ複素指向性12を測定し、その測定情報(つまり、アンテナ複素指向性12)を記憶装置13に記憶させる。なお、記憶装置13には、あらかじめ、実施の形態1で定義した幾何光学的解析法の前処理30や後処理31を行うプログラムや複数の電波伝搬モデル71が記憶されている。そして、記憶装置13は、アンテナ特性測定装置61より取り込んで自己に記憶されたアンテナ複素指向性12を入力パラメータとして、幾何光学的電波伝搬特性解析部14の前処理30及び後処理31に渡す。つまり、記憶装置13が媒体となって、アンテナ特性測定装置61が測定したアンテナ複素指向性12を、電波伝搬特性解析を行うための前処理30及び後処理31に渡す。その後、オペレータは、例えば市街地モデルや郊外モデルなどのあらかじめ記憶された複数の電波伝搬モデル71の中から所望する1つの電波伝搬モデル71を選択する。   Next, the operation of the antenna characteristic evaluation system according to Embodiment 2 will be described with reference to FIGS. First, the control device 51 measures the antenna complex directivity 12 such as the amplitude and phase of radio waves using an antenna characteristic measuring device 61 such as a network analyzer, and stores the measurement information (that is, the antenna complex directivity 12) as a storage device. 13 is stored. The storage device 13 stores in advance a program for performing pre-processing 30 and post-processing 31 of the geometric optical analysis method defined in Embodiment 1, and a plurality of radio wave propagation models 71. Then, the storage device 13 passes the antenna complex directivity 12 acquired from the antenna characteristic measurement device 61 and stored therein as an input parameter to the preprocessing 30 and the postprocessing 31 of the geometric optical wave propagation characteristic analysis unit 14. That is, the storage device 13 serves as a medium, and the antenna complex directivity 12 measured by the antenna characteristic measuring device 61 is passed to the preprocessing 30 and the postprocessing 31 for performing the radio wave propagation characteristic analysis. Thereafter, the operator selects one desired radio wave propagation model 71 from a plurality of radio wave propagation models 71 stored in advance, such as an urban area model and a suburban model.

このような本発明の実施の形態2のアンテナ特性評価システムによれば、ネットワークアナライザなどのアンテナ特性測定装置61より得られた振幅や位相などのアンテナ複素指向性12を、記憶装置13(例えば、パソコンなど)を媒体として幾何光学的電波伝搬特性解析部14の前処理30や後処理31の送受信地点に用いている。これにより、オペレータが電波伝搬モデル71を選択することにより、電波伝搬モデル71の詳細なパラメータをその都度入力することなく、比較的短時間で実際の多重波環境におけるアンテナ特性を計算することができる。   According to the antenna characteristic evaluation system of the second embodiment of the present invention, the antenna complex directivity 12 such as the amplitude and phase obtained from the antenna characteristic measuring device 61 such as a network analyzer is stored in the storage device 13 (for example, It is used as a transmission / reception point for the pre-processing 30 and post-processing 31 of the geometric optical radio wave propagation characteristic analyzing unit 14 using a computer or the like as a medium. Thus, when the operator selects the radio wave propagation model 71, the antenna characteristics in the actual multiwave environment can be calculated in a relatively short time without inputting detailed parameters of the radio wave propagation model 71 each time. .

なお、前処理30において送信点におけるアンテナ複素指向性12を所定のものに限定してあらかじめ前処理30と後処理31を実行しておき、その後、アンテナ特性測定装置61を用いて得られた受信点におけるアンテナ複素指向性12を到達光線に乗算して位相合成するように構成すれば、さらに短時間で実際の多重波環境におけるアンテナ特性を計算することができる。   In the pre-processing 30, the antenna complex directivity 12 at the transmission point is limited to a predetermined value, and the pre-processing 30 and the post-processing 31 are executed in advance, and then the reception obtained using the antenna characteristic measuring device 61 is performed. If the antenna complex directivity 12 at the point is multiplied by the arrival light and phase-combined, the antenna characteristics in the actual multiwave environment can be calculated in a shorter time.

次に、フローチャートを用いて本発明の実施の形態2におけるアンテナ特性評価システムの動作を説明する。図8は、本発明の実施の形態2における電波伝搬特性解析の後処理Aの流れを示すフローチャートである。また、図9は、図8に示す電波伝搬特性解析の後処理Aによって記憶装置に記憶されるデータ構成を示す図である。図9に示すように、記憶装置13には、データ80として、到達光線の方程式81、到達光線の進行角度82、到達光線の偏波方程式83、及び到達光線の伝搬距離84が記憶される。さらに、図10は、本発明の実施の形態2における電波伝搬特性解析の後処理Bの流れを示すフローチャートである。したがって、図8、図9及び図10を用いて実施の形態2における電波伝搬特性解析の後処理の流れを説明する。   Next, the operation of the antenna characteristic evaluation system according to Embodiment 2 of the present invention will be described using a flowchart. FIG. 8 is a flowchart showing a flow of post-processing A of radio wave propagation characteristic analysis according to Embodiment 2 of the present invention. FIG. 9 is a diagram showing a data configuration stored in the storage device by post-processing A of the radio wave propagation characteristic analysis shown in FIG. As shown in FIG. 9, the storage device 13 stores, as data 80, a reaching light equation 81, a reaching light traveling angle 82, a reaching light polarization equation 83, and a reaching light propagation distance 84. Furthermore, FIG. 10 is a flowchart showing a flow of post-processing B of radio wave propagation characteristic analysis in Embodiment 2 of the present invention. Therefore, the flow of post-processing of radio wave propagation characteristic analysis in the second embodiment will be described with reference to FIG. 8, FIG. 9, and FIG.

最初に、図8の電波伝搬特性解析の後処理Aについて説明する。まず、パラメータ入力の手順において、オペレータは、受信位置及び電波の周波数をパラメータとして入力装置54から入力する(ステップST31)。そして、制御装置51は、受信球導出の手順において、記憶装置13から読み込んだ光線番号aの受信地点までの距離dを導き出し、送信角度間隔をφとしたとき、後処理Aにおける受信球の半径rを前述の式(1)によって求める(ステップST32)。   First, post-processing A of radio wave propagation characteristic analysis in FIG. 8 will be described. First, in the parameter input procedure, the operator inputs the reception position and radio wave frequency as parameters from the input device 54 (step ST31). Then, the control device 51 derives the distance d to the reception point of the ray number a read from the storage device 13 in the procedure for deriving the reception sphere, and when the transmission angle interval is φ, the radius of the reception sphere in the post-processing A r is obtained by the above-described equation (1) (step ST32).

次に、制御装置51は、反射光線方程式と受信点を中心とする半径rの球との交点を求めることにより到達光線の探索を行い、電波が受信点に到達したか否かを判定する(ステップST33)。ここで、電波が受信点に到達した場合は(ステップST33でYESの場合)、到達光線データ出力の手順に移行する(ステップST34)。一方、光線番号aの光線が受信点に到達していないときは(ステップST33でNOの場合)、光線番号aに1を加算してa+1にして(ステップST35)、ステップST32の受信球導出のステップに戻り、次の光線番号(a+1)の受信球導出処理に移行する。   Next, the control device 51 searches for a reaching ray by obtaining an intersection of the reflected ray equation and a sphere having a radius r centered on the receiving point, and determines whether or not the radio wave has reached the receiving point ( Step ST33). Here, when the radio wave arrives at the reception point (YES in step ST33), the procedure proceeds to the procedure for outputting the reaching ray data (step ST34). On the other hand, when the ray with the ray number a has not reached the reception point (NO in step ST33), 1 is added to the ray number a to become a + 1 (step ST35), and the reception sphere derivation in step ST32 is performed. Returning to the step, the process proceeds to the reception sphere derivation process for the next ray number (a + 1).

また、ステップST34の到達光線データ出力の手順に移行した場合は、図9に示すように、到達光線の方程式81、到達光線の進行角度82、到達光線の偏波方程式83、及び到達光線の伝搬距離84からなるデータ80を記憶装置13に記憶させる(ステップST34)。そして、放射光線数Nαが光線番号aより小さくなったか否か(つまり、Nα<aか否か)を判定し(ステップST36)、Nα<aでなければ(ステップST36でNOの場合)、光線番号aに1を加算してa+1にして(ステップST35)、前述のステップST32の受信球導出の手順に戻り、Nα<aとなるまでステップST32からステップST36の処理を繰り返す。   Further, when the procedure proceeds to the procedure of the reaching ray data output in step ST34, as shown in FIG. 9, the reaching ray equation 81, the reaching ray traveling angle 82, the reaching ray polarization equation 83, and the propagation of the reaching ray. Data 80 consisting of the distance 84 is stored in the storage device 13 (step ST34). Then, it is determined whether or not the number of emitted light rays Nα is smaller than the light ray number a (that is, whether or not Nα <a) (step ST36). If Nα <a (NO in step ST36), the light ray 1 is added to the number a to make it a + 1 (step ST35), the process returns to the reception sphere derivation procedure of step ST32 described above, and the processing from step ST32 to step ST36 is repeated until Nα <a.

一方、放射光線数Nαが光線番号aより小さい場合、つまり、Nα<aの場合は(ステップST36でYESの場合)、予め設定されたP個の受信地点を移動して前述の各ステップと同様の処理を行い、これから移動する受信地点の受信点数Pが受信地点mを下回ったか否か(つまり、P<mか否か)を判定し(ステップST37)、P<mでなければ(ステップST37でNOの場合)、受信地点mを追尾しながら(つまり、受信地点をm+1にして)、オペレータが設定したP個の受信地点を移動し(ステップST38)、ステップST32の受信球の導出処理に戻って前述と同様のステップで処理を繰り返す。このようにして、これから移動する受信地点の受信点数Pが受信地点mを下回ったら、つまり、P<mとなったら(ステップST37でYESの場合)、処理を終了する。   On the other hand, if the number of radiated rays Nα is smaller than the ray number a, that is, if Nα <a (in the case of YES in step ST36), P preset reception points are moved and the same as the above-described steps. It is determined whether or not the number P of reception points to be moved from now on is less than the reception point m (that is, whether P <m) (step ST37). If P <m is not satisfied (step ST37). In the case of NO), while tracking the reception point m (that is, the reception point is set to m + 1), the P reception points set by the operator are moved (step ST38), and the reception sphere derivation process in step ST32 is performed. Return and repeat the process in the same steps as described above. In this way, when the number P of reception points to be moved from now on falls below the reception point m, that is, when P <m (in the case of YES in step ST37), the processing is terminated.

次に、図10の電波伝搬特性解析の後処理Bについて説明する。まず、オペレータは、アンテナ特性測定装置61で得られたアンテナ複素指向性12を入力装置54から入力する(ステップST41)。そして、制御装置51は、図8の後処理Aで記憶装置13に記憶されたデータ80(つまり、到達光線の方程式81、到達光線の進行角度82、到達光線の偏波方程式83、及び到達光線の伝搬距離84)の読み込みを行う(ステップST42)。さらに、制御装置51は、これらのデータ80に対して到来する光線に対応した角度のアンテナ複素指向性12を乗算する(ステップST43)。   Next, post-processing B of radio wave propagation characteristic analysis in FIG. 10 will be described. First, the operator inputs the antenna complex directivity 12 obtained by the antenna characteristic measuring device 61 from the input device 54 (step ST41). Then, the control device 51 stores the data 80 (that is, the reaching ray equation 81, the reaching ray traveling angle 82, the reaching ray polarization equation 83, and the reaching ray) stored in the storage device 13 in the post-processing A of FIG. Is read (step ST42). Further, the control device 51 multiplies these data 80 by the antenna complex directivity 12 having an angle corresponding to the incoming light beam (step ST43).

このようにして、データ80を読み込むべき到来光線数Qが所定の到来光線番号qを下回ったか否か(つまり、Q<qか否か)を判定し(ステップST44)、Q<qでなければ(ステップST44でNOの場合)、到来光線番号qに1を加算してq+1とし(ステップST45)、ステップST42に戻ってQ<qとなるまでステップST42からステップST44の処理を繰り返す。   In this way, it is determined whether or not the number of incoming rays Q from which data 80 is to be read has fallen below a predetermined incoming ray number q (that is, whether Q <q) (step ST44). (In the case of NO in step ST44), 1 is added to the incoming ray number q to obtain q + 1 (step ST45), and the process returns to step ST42 and repeats the processing from step ST42 to step ST44 until Q <q.

このようにして前述の処理をQ回繰り返し、データ80を読み込むべき到来光線数Qが所定の到来光線番号qを下回ったら、つまり、Q<qとなったら(ステップST44でYESの場合)、位相合成の手順においてQ本の到来光線を位相合成し(ステップST46)、さらに、電界強度導出の手順において受信点における電界強度を算出する(ステップST47)。   In this way, the above-described processing is repeated Q times, and if the number of incoming rays Q from which data 80 is to be read falls below a predetermined incoming ray number q, that is, if Q <q (in the case of YES in step ST44), the phase The Q incoming rays are phase-combined in the combining procedure (step ST46), and the electric field strength at the receiving point is calculated in the electric field strength derivation procedure (step ST47).

次に、受信地点の受信点数Pが受信地点mを下回ったか否か(つまり、P<mか否か)を判定し(ステップST48)、P<mでなければ(ステップST48でNOの場合)、受信地点mを追尾しながら(つまり、受信地点をm+1にして)、オペレータが設定したP個の受信地点を移動し(ステップST49)、ステップST42の到来光線のデータ読み込み処理に戻って前述と同様のステップで処理を繰り返す。以上の手順によって、これから移動する受信地点の受信点数Pが受信地点mを下回ったら、つまり、P<mとなったら(ステップST48でYESの場合)、処理を終了する。言い換えれば、あらかじめ設定されたP個の受信地点を移動して同様の処理を行う。   Next, it is determined whether or not the number of reception points P at the reception point is lower than the reception point m (that is, whether P <m) (step ST48). If P <m is not satisfied (NO in step ST48). Then, while tracking the reception point m (that is, the reception point is set to m + 1), the P reception points set by the operator are moved (step ST49), and the process returns to the incoming ray data reading process of step ST42. The process is repeated in the same steps. If the number P of reception points to be moved from now on is less than the reception point m by the above procedure, that is, if P <m (YES in step ST48), the process is terminated. In other words, the same processing is performed by moving P reception points set in advance.

以上説明したように、図10に示す後処理Bのみを計算することにより、あらかじめ定められた実際の多重波環境におけるアンテナ特性を比較的短時間で求めることができ、アンテナ特性の比較検討を容易に行うことができる。また、後処理Bで必要なデータ80は受信点に到達する光線のみであるため、アンテナ特性を評価する際に計算機システムの記憶装置13に記憶されるデータ量は数十メガバイト程度となる。このため、アンテナ特性の計算量も低減されるので、ノートパソコンなどの比較的計算性能が低い計算機システムにおいても、アンテナ特性の計算を行うことが可能となる。なお、以上述べた説明では、アンテナ特性測定装置61としてネットワークアナライザで測定されたアンテナ複素指向性を用いたが、これに限ることなく、電磁界計算によって得られたアンテナ複素指向性を用いてもよい。   As described above, by calculating only the post-processing B shown in FIG. 10, it is possible to obtain the antenna characteristics in a predetermined actual multiwave environment in a relatively short time, and it is easy to compare and examine the antenna characteristics. Can be done. Further, since the data 80 necessary for the post-processing B is only the light beam reaching the reception point, the amount of data stored in the storage device 13 of the computer system when evaluating the antenna characteristics is about several tens of megabytes. For this reason, since the amount of calculation of antenna characteristics is also reduced, it is possible to calculate antenna characteristics even in a computer system with relatively low calculation performance such as a notebook personal computer. In the above description, the antenna complex directivity measured by the network analyzer is used as the antenna characteristic measuring device 61. However, the present invention is not limited to this, and the antenna complex directivity obtained by the electromagnetic field calculation may be used. Good.

本発明にかかるアンテナ特性評価システムは、アンテナ指向性解析部により得られたアンテナ複素指向性を幾何光学的電波伝搬特性解析部の送受信地点における指向性係数として用いることにより、実際の多重波伝搬環境における高精度なアンテナ特性評価を行うことができる。したがって、計算機による高精度なアンテナ特性評価システムとして有効に利用することができる。   The antenna characteristic evaluation system according to the present invention uses an antenna complex directivity obtained by an antenna directivity analysis unit as a directivity coefficient at a transmission / reception point of a geometric optical radio wave propagation characteristic analysis unit, thereby realizing an actual multiwave propagation environment. The antenna characteristics can be evaluated with high accuracy. Therefore, it can be effectively used as a highly accurate antenna characteristic evaluation system by a computer.

本発明の実施の形態1のアンテナ特性評価システムにおけるアンテナ指向性解析部と幾何光学的電波伝搬特性解析部の統合を示すブロック図The block diagram which shows integration of the antenna directivity analysis part and geometric optical wave propagation characteristic analysis part in the antenna characteristic evaluation system of Embodiment 1 of this invention 図1に示す幾何光学的電波伝搬特性解析部が行う電波伝搬特性解析の前処理の流れを示すフローチャートThe flowchart which shows the flow of the pre-processing of the radio wave propagation characteristic analysis which the geometric optical wave propagation characteristic analysis part shown in FIG. 1 performs 図1に示す幾何光学的電波伝搬特性解析部が行う電波伝搬特性解析の後処理の流れを示すフローチャートThe flowchart which shows the flow of the post-process of the radio wave propagation characteristic analysis which the geometric optical wave propagation characteristic analysis part shown in FIG. 1 performs 図2のフローチャートで示される前処理によって記憶装置に記憶されるデータ構成を示す図The figure which shows the data structure memorize | stored in a memory | storage device by the pre-process shown by the flowchart of FIG. 本発明の実施の形態1における電波伝搬特性解析の計算を行うためのコンピュータシステムの構成図Configuration diagram of a computer system for performing calculation of radio wave propagation characteristic analysis in Embodiment 1 of the present invention 本発明の実施の形態2における電波伝搬特性解析の計算を行うためのコンピュータシステムの構成図Configuration diagram of a computer system for performing calculation of radio wave propagation characteristic analysis in Embodiment 2 of the present invention 図6の記憶装置が行う電波伝搬特性の解析処理を機能的に示したブロック図Block diagram functionally showing the radio wave propagation characteristic analysis processing performed by the storage device of FIG. 本発明の実施の形態2における電波伝搬特性解析の後処理Aの流れを示すフローチャートThe flowchart which shows the flow of the post-process A of the electromagnetic wave propagation characteristic analysis in Embodiment 2 of this invention 図8に示す電波伝搬特性解析の後処理Aによって記憶装置に記憶されるデータ構成を示す図The figure which shows the data structure memorize | stored in a memory | storage device by the post-process A of the radio wave propagation characteristic analysis shown in FIG. 本発明の実施の形態2における電波伝搬特性解析の後処理Bの流れを示すフローチャートThe flowchart which shows the flow of the post-process B of the radio wave propagation characteristic analysis in Embodiment 2 of this invention

符号の説明Explanation of symbols

1 アンテナ特性評価システム
11 アンテナ指向性解析部
12 アンテナ複素指向性
13 記憶装置
14 幾何光学的電波伝搬特性解析部
30 前処理
31 後処理
41 全光線の方程式
42 全光線の進行角度
43 全光線の偏波方程式
44 全光線の反射面番号
45 全光線の壁面番号
51 制御装置
52 表示装置
53 印刷装置
54 入力装置
55 記憶媒体
61 アンテナ特性測定装置
62 被測定アンテナ
63 受信アンテナ(又は送信アンテナ)
71 電波伝搬モデル
80 データ
81 到達光線の方程式
82 到達光線の進行角度
83 到達光線の偏波方程式
84 到達光線の伝搬距離
DESCRIPTION OF SYMBOLS 1 Antenna characteristic evaluation system 11 Antenna directivity analysis part 12 Antenna complex directivity 13 Storage device 14 Geometric optical wave propagation characteristic analysis part 30 Preprocessing 31 Postprocessing 41 Equation of all rays 42 Progression angle of all rays 43 Deviation of all rays Wave equation 44 Reflecting surface number of all rays 45 Wall number of all rays 51 Control device 52 Display device 53 Printing device 54 Input device 55 Storage medium 61 Antenna characteristic measuring device 62 Antenna to be measured 63 Receiving antenna (or transmitting antenna)
71 Radio propagation model 80 Data 81 Reaching ray equation 82 Traveling angle of reaching ray 83 Polarization equation of reaching ray 84 Propagating distance of reaching ray

Claims (6)

無線局から送信された電波の伝搬特性を計算してアンテナ特性の評価を行うアンテナ特性評価システムであって、
第1アンテナ及び第2アンテナの離散的角度方位に対する複素指向性係数列を算出するアンテナ指向性解析手段と、
送信点より複数の光線を放射し、その光線が受信点に到達するまでの軌跡を追跡することで所定の電波伝搬モデルにおける送受信点間の電波の伝搬特性を算出する幾何光学的電波伝搬特性解析手段とを備え、
前記幾何光学的電波伝搬特性解析手段が、前記第1アンテナ及び前記第2アンテナの前記複素指向性係数列を前記送受信点におけるアンテナの指向性係数として用いることにより、前記アンテナ指向性解析手段と前記幾何光学的電波伝搬特性解析手段とを統合することを特徴とするアンテナ特性評価システム。
An antenna characteristic evaluation system that evaluates antenna characteristics by calculating the propagation characteristics of radio waves transmitted from a radio station,
Antenna directivity analysis means for calculating a complex directivity coefficient sequence for discrete angular orientations of the first antenna and the second antenna;
Geometric optical wave propagation characteristics analysis that calculates the propagation characteristics of radio waves between transmission and reception points in a given radio wave propagation model by radiating multiple rays from the transmission point and tracing the trajectory until the rays reach the reception point Means and
The geometric optical radio wave propagation characteristic analyzing means uses the complex directivity coefficient sequence of the first antenna and the second antenna as a directivity coefficient of the antenna at the transmission / reception point, so that the antenna directivity analysis means and the antenna directivity analysis means An antenna characteristic evaluation system that integrates geometric optical wave propagation characteristic analysis means.
さらに、前記送信点におけるアンテナの前記複素指向性係数列の中で、離散的角度方位に対する振幅が所定の閾値を越えるか否かを判定する振幅判定手段を備え、
前記振幅判定手段の判定結果に応じて、前記送信点から光線を放射するか否かを決定することを特徴とする請求項1に記載のアンテナ特性評価システム。
Furthermore, an amplitude determination means for determining whether or not an amplitude with respect to a discrete angle azimuth exceeds a predetermined threshold in the complex directivity coefficient sequence of the antenna at the transmission point,
2. The antenna characteristic evaluation system according to claim 1, wherein whether or not to radiate a light beam from the transmission point is determined according to a determination result of the amplitude determination unit.
前記アンテナ指向性解析手段は有限差分時間領域法を用いることを特徴とする請求項1又は請求項2に記載のアンテナ特性評価システム。   The antenna characteristic evaluation system according to claim 1 or 2, wherein the antenna directivity analysis means uses a finite difference time domain method. さらに、前記第1アンテナ及び前記第2アンテナの離散的角度方位に対する前記複素指向性係数列を測定するアンテナ指向性測定手段を備え、
前記幾何光学的電波伝搬特性解析手段が、前記アンテナ指向性測定手段によって得られた前記複素指向性係数列を送受信点におけるアンテナの指向性係数として用いることを特徴とする請求項1又は請求項2に記載のアンテナ特性評価システム。
And antenna directivity measurement means for measuring the complex directivity coefficient sequence with respect to discrete angular orientations of the first antenna and the second antenna,
3. The geometric optical wave propagation characteristic analyzing means uses the complex directivity coefficient sequence obtained by the antenna directivity measuring means as an antenna directivity coefficient at a transmission / reception point. Antenna characteristic evaluation system described in 1.
前記幾何光学的電波伝搬特性解析手段が保有する複数の電波伝搬モデルを記憶する記憶手段と、
前記第1アンテナ及び前記第2アンテナの前記複素指向性係数列を用いて、あらかじめ指定した前記複数の電波伝搬モデルに対して幾何光学的電波伝搬特性解析を実行する実行手段と、
を備えることを特徴とする請求項1から請求項4のいずれかに記載のアンテナ特性評価システム。
Storage means for storing a plurality of radio wave propagation models possessed by the geometric optical radio wave propagation characteristic analysis means;
Executing means for performing geometric optical wave propagation characteristic analysis on the plurality of wave propagation models specified in advance using the complex directivity coefficient sequences of the first antenna and the second antenna;
The antenna characteristic evaluation system according to any one of claims 1 to 4, further comprising:
前記幾何光学的電波伝搬特性解析手段は、
前記記憶手段にあらかじめ記憶された前記電波伝搬モデルに対して幾何光学的電波伝搬特性の解析計算を実行し、
前記アンテナ指向性解析手段又は前記アンテナ指向性測定手段によって得られたアンテナの前記複素指向性係数列を用いて位相合成の計算を行うことを特徴とする請求項5に記載のアンテナ特性評価システム。
The geometric optical radio wave propagation characteristic analyzing means includes:
Execute geometric optical wave propagation characteristic analysis calculation for the radio wave propagation model stored in advance in the storage means,
6. The antenna characteristic evaluation system according to claim 5, wherein a phase synthesis calculation is performed using the complex directivity coefficient sequence of the antenna obtained by the antenna directivity analysis means or the antenna directivity measurement means.
JP2004150381A 2004-05-20 2004-05-20 Antenna characteristic evaluation system Pending JP2005333451A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2004150381A JP2005333451A (en) 2004-05-20 2004-05-20 Antenna characteristic evaluation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2004150381A JP2005333451A (en) 2004-05-20 2004-05-20 Antenna characteristic evaluation system

Publications (1)

Publication Number Publication Date
JP2005333451A true JP2005333451A (en) 2005-12-02

Family

ID=35487790

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2004150381A Pending JP2005333451A (en) 2004-05-20 2004-05-20 Antenna characteristic evaluation system

Country Status (1)

Country Link
JP (1) JP2005333451A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020024165A (en) * 2018-08-08 2020-02-13 株式会社Nttドコモ Radiated power estimation method

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2020024165A (en) * 2018-08-08 2020-02-13 株式会社Nttドコモ Radiated power estimation method
JP7016303B2 (en) 2018-08-08 2022-02-04 株式会社Nttドコモ Radiation power estimation method

Similar Documents

Publication Publication Date Title
CN112612024A (en) Microwave array rapid imaging method
JP4530898B2 (en) Radio wave propagation estimation program, radio wave propagation estimation method, and apparatus for executing the method
US11619701B2 (en) Satellite tracking system and method thereof
JP5155934B2 (en) Radio station location search device
Bahador et al. Three-stage Lamb-wave-based damage localization algorithm in plate-like structures for structural health monitoring applications
JP2007101376A (en) System and method for estimating radio wave arrival condition, and program
Ollander et al. Simulation-based analysis of multipath delay distributions in urban canyons
JP2003318811A (en) Device, method, program for estimating and calculating reception electric field intensity and recording medium
CN112069713B (en) Near-field scattering characteristic modeling method, electronic device and storage medium
Varela et al. Planar Wide Mesh Scanning using Multi-Probe Systems
Fu et al. A low SNR and fast passive location algorithm based on virtual time reversal
JP2005333451A (en) Antenna characteristic evaluation system
US20090167756A1 (en) Method and apparatus for computation of wireless signal diffraction in a three-dimensional space
Zuo et al. Position error detection and compensation for far-field radar cross-section measurement
KR20160077973A (en) Apparatus and method for simulating analysis about tactical air navigation system
KR100329518B1 (en) Prediction and analysis method of cell planning for wireless communication
JP2002107397A (en) Shielding material detecting method in characteristic estimation for radio wave propagation, and system and recording medium thereof
JP2001153906A (en) Instrument and method for measuring antenna
JP7315349B2 (en) Position estimation device and position estimation method
JP4372287B2 (en) Radio wave propagation analysis method and apparatus, and storage medium
Fabrizio et al. Single site geolocation method for a linear array
US10379198B2 (en) Determining positions of transducers for receiving and/or transmitting wave signals
JP7241993B2 (en) Radio wave source position estimation device, radio wave source position estimation system, and radio wave source position estimation method
JP2006135844A (en) Method and device for geometrical-optical radiowave propagation characteristics analysis, and recording medium
KR102359201B1 (en) METHOD AND APPARATUS FOR ANALYZING PROPAGATION CHARACTERISTIC BASED ON ANTENNA REFLECTOR Of GEOSTATIONARY SATELLITE AND SPACE PROPAGATION ENVIRONMENT