JP3605229B2 - Digital phase modulation method display method - Google Patents
Digital phase modulation method display method Download PDFInfo
- Publication number
- JP3605229B2 JP3605229B2 JP13281796A JP13281796A JP3605229B2 JP 3605229 B2 JP3605229 B2 JP 3605229B2 JP 13281796 A JP13281796 A JP 13281796A JP 13281796 A JP13281796 A JP 13281796A JP 3605229 B2 JP3605229 B2 JP 3605229B2
- Authority
- JP
- Japan
- Prior art keywords
- data
- phase modulation
- digital phase
- domain data
- time domain
- 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
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Description
【0001】
【発明の属する技術分野】
本発明は、未知の無線機からの受信信号のデジタル位相変調方式を判定するためのデジタル位相変調方式判定表示方法に関する。
【0002】
【従来の技術及び発明が解決しようとする課題】
無線機のデジタル位相変調方式として、PSK、QPSK、π/4シフトQPSK、16QAM、256QAM等の種々のデジタル位相変調方式が存在する。このようなデジタル位相変調方式は、受信信号の搬送波周波数及び位相が既知であれば、その精確な周波数及び位相で受信信号を同期検波し、その検波したデータをグラフ表示することにより容易に確認することができる。すなわち、各デジタル位相変調方式に応じて、位相の変化点に対応した特有のシンボル・パターンが決まっているからである。このように、従来の信号解析システムでは、デジタル位相変調方式が既知の入力信号を解析することにより、精確な搬送波周波数及び位相に基づいて適正な同期検波を行うことができる。
【0003】
しかし、未知の無線機からの入力信号の場合には、デジタル位相変調方式が不明であり、その信号の搬送波周波数及び位相も未知なので、その入力信号を適正に同期検波することができない。入力信号から搬送波信号を抽出するためには、デジタル位相変調方式が既知でなければならず、未知のデジタル位相変調方式の入力信号に対しては、考えられる全ての方式に夫々対応した多数の回路を用意し、適正な同期検波ができるまで切替使用する如き方法しかないことになるが、このような試行錯誤の解析方法は極めて煩雑でコストも嵩む欠点がある。したがって、未知の無線機からの入力信号を精確に同期検波してグラフ表示することは困難であり、未知の無線機からの受信信号からその未知の無線機のデジタル位相変調方式を判定することも極めて困難であった。
【0004】
したがって、本発明の目的は、未知の無線機からの入力信号のデジタル位相変調方式を精確かつ簡単に判定することが可能なデジタル位相変調方式判定表示方法を提供することである。
【0005】
【課題を解決する為の手段】
本発明は、未知の無線機からの入力信号のデジタル位相変調方式を自動的に判定するデジタル位相変調方式判定表示方法を提供する。前提として、受信信号の同相成分(I成分)及び直交成分(Q成分)を直交座標平面又は極座標平面上に表示すると、デジタル位相変調方式に応じてシンボル・パターンが異なることを利用するものである。よって、未知の無線機からの入力信号からデジタル位相変調方式を判定するための表示を行うため、以下のような手順を採用した。先ず、未知のデジタル位相変調信号の周波数領域データを生成し、この周波数領域データから上記未知のデジタル位相変調信号の見積搬送波周波数を求める。次に、この見積搬送波周波数に基づいて未知のデジタル位相変調信号の時間領域データを検波演算し、時間領域の直交座標データを求め、この直交座標データの各点において何れか一方の座標が極値となるシンボル点を求め、これらのシンボル点のみを複数の分割領域で構成された座標平面上にプロットし、上記複数の分割領域の中で上記シンボル点がプロットされなかった分割領域の数を計数するという一連の動作を見積搬送波周波数を変化させながら順次繰り返し実行する。この結果、シンボル点がプロットされなかった分割領域の数が最大となるときの見積搬送波周波数が真の搬送波周波数として求められる。この搬送波周波数に基づいて、未知のデジタル位相変調信号の時間領域データを検波演算し、この結果得られたデータを表示スクリーン上に表示する。
【0006】
本発明によれば、時間領域の直交座標データを求める検波演算は、比較的計算負荷が重いので、計算負荷を低減する方法も採用できる。すなわち、全てのシンボル点を直交座標データから求めるのではなく、直交座標データを極座標データに座標変換し、前に求めたシンボル点に対応するデータを極座標データから抽出する。これにより、直交座標データからシンボル点を直接求める計算負荷が格段に軽減し、処理の高速化が図られる。
【0007】
【発明の実施の形態】
図4は、本発明のデジタル位相変調方式判定表示方法を使用するのに好適な解析システムの構成の一例を示すブロック図である。アンテナ10は、未知又は既知の無線機の出力信号の電波を受信するものである。このアンテナ10により受信された信号は、受信器及び周波数変換器として機能するダウンコンバータ12に供給され、中間周波数信号に変換される。この中間周波数信号は、デジタル・スペクトラム・アナライザ14に供給され、デジタル・データ列(時間領域)に変換され、更にFFT処理によりデジタル周波数成分データに変換され、メモリ(図示せず)に記憶される。その後、メモリ内のデジタル周波数成分データは、コンピュータ16に送られる。このコンピュータ16は、表示スクリーン18を有し、測定結果、解析結果、システムの設定状態等のシステム全体の制御に必要な情報の表示をする。ダウンコンバータ12としては、例えば、ソニー・テクトロニクス社製の2784型スペクトラム・アナライザ等が好適であり、デジタル・スペクトラム・アナライザ14としては、ソニー・テクトロニクス社製の3055型リアルタイム・スペクトラム・アナライザ等が好適である。しかし、その他の適当な測定機器を使用できることは勿論である。ただし、デジタル・スペクトラム・アナライザ14としては、解析精度の向上のために時間的に切れ目なくデータを解析することが望ましいので、リアルタイム解析機能を有することが好適である。なお、デジタル・スペクトラム・アナライザにおけるリアルタイム解析機能については、当業者には周知であり、本願発明と直接関係するものではないので説明は省略する。
【0008】
データ解析は、コンピュータ16がソフトウェアに従って実行する。図1及び図2は、本発明のデジタル位相変調方式判定表示方法に従ってコンピュータ16が実行する処理の流れ図の一例である。先ず、アンテナ10及びダウン・コンバータ12を介してデジタル・スペクトラム・アナライザ14に入力された未知の無線機からの信号の周波数領域データを生成する(ステップ100)。具体的には、入力信号をデジタイズして時間領域データを生成し、FFT(高速フーリエ変換)処理を実行して周波数領域データを生成する。次に、生成した周波数領域データに基づき、見積搬送波周波数を計算する(ステップ102)。この搬送波周波数の見積もり計算は、周波数領域データの振幅の全体の累算値の2分の1の値に対応する周波数値を採用する。このような見積もり計算で求めた見積搬送波周波数の値は、実際の入力信号の搬送波周波数に一致するとは限らないが、後述するような処理により、徐々に真の値に近似させていくことができる。
【0009】
次に、見積搬送波周波数に基づいて入力信号の時間領域データを検波演算し、時間領域の直交座標データを求める(ステップ104)。この場合に使用する入力信号の時間領域データは、デジタル・スペクトラム・アナライザ14で使用したデータをそのまま用いても良いし、ステップ100で生成した周波数領域データに対して逆FFT演算処理を実行して求めても良い。また、検波演算とは、ハードウエアによる検波と類似の処理をコンピュータ16の演算により行うことを意味する。すなわち、時間領域の直交座標データは、I(同相)軸とQ(直交)軸の直交座標平面上で表されるものである。この直交座標データは、入力信号のデジタル位相変調方式が既知であり、搬送波周波数及び位相も既知であれば、精確に同期検波が可能で、I−Q座標平面上にデータをプロットすると、各デジタル位相変調方式に対応したシンボル・パターンが表示されるものであるが、本発明においては、この時点ではデジタル位相変調方式も搬送波周波数も未知であり、近似的な検波演算を行っているに過ぎない。具体的には、検波演算は、以下のような数式1及び数式2に基づいて実行される。
【0010】
【数1】
a=cosωx
b=sinωx
ω=2πFT
F:周波数領域幅の中心周波数と見積搬送波周波数との差
T:サンプリング周期
x:データ変数
【数2】
I(x)=a・i(x)+b・q(x)
Q(x)=a・q(x)−b・i(x)
i(x):時間領域の同相成分データ
q(x):時間領域の直交成分データ
I(x):直交座標データのI軸成分データ
Q(x):直交座標データのQ軸成分データ
上記の検波演算を実行することにより、時間領域の直交座標データI(x)及びQ(x)が得られる。ここで得られた直交座標データは、精確な同期検波により得られたものではなく、近似的な見積もり計算の結果であることに留意されたい。
【0011】
次に、この直交座標データの各点において、何れか一方の座標が極値となるシンボル点を求める(ステップ106)。シンボル点は、位相変調における位相の変化時点に対応すると考えられる点であり、位相の変化時点において、同相成分座標I(x)又は直交成分データQ(x)の何れか一方が極値(極大値又は極小値)となると考えられるからである。
【0012】
次に、ステップ106で求めたシンボル点を複数の分割領域で構成された座標平面上にプロットする(ステップ108)。この場合の座標平面は、直交座標平面でも極座標平面でも構わない。どのような座標平面でも、複数の小さな分割領域で区分されていれば良い。また、各分割領域の形状も特に限定する必要もないので、採用する座標平面に都合の良い分割領域を採用できる。また、ここで言う「プロット」という処理は、実際に座標平面をスクリーン上に表示して行う必要は全くない。座標平面上の座標データを記録するだけで十分である。つまり、その座標平面上の複数の領域の各々にシンボル点が存在するか否かをコンピュータ16が判断できれば十分であり、観測者が判断する必要はない。
【0013】
次に、座標平面を構成する複数の分割領域の中で、シンボル点が1つもプロットされなかった分割領域の数を計数する(ステップ110)。その後、見積搬送波周波数を所定値だけ変化させる(ステップ112)。以上のステップ104から112までの処理を見積搬送波周波数の検索範囲に亘って順次繰り返し実行する。検索範囲の全ての処理が終了すると(ステップ114)、シンボル点がプロットされなかった分割領域の数が最大となる見積搬送波周波数を求め、これを真の搬送波周波数と見なし(ステップ116)、その周波数に基づいて入力信号の時間領域データを再度検波演算する(118)。この結果得られたデータを表示スクリーン上の座標平面(I−Q座標平面)にグラフ表示する(ステップ120)。この結果、入力信号のデジタル位相変調方式に対応したシンボル・パターンが表示されるので、通常の知識を有する当業者は、一目瞭然に未知の入力信号のデジタル位相変調方式を判断することができる。
【0014】
見積搬送波周波数の検索範囲は、最初に求めた見積搬送波周波数を中心とした所定の周波数範囲であり、例えば、周波数領域データ幅が200kHzの場合で6.25kHzである。そして、検索処理の際の見積搬送波周波数の検索ステップは、例えば、約4Hzである。これらの値は一例であって、他の検索範囲及び検索ステップで実行できることは勿論である。
【0015】
上述の処理において、ステップ104で実行する検波演算の処理は、三角関数の計算を含んでいるので、比較的時間がかかる。そこで、全ての検索範囲内の検索ステップ毎に数式1及び数式2の計算を実行するのは必ずしも適切ではない。勿論、理論的には全てステップ104の検波演算を行えば良いのであるが、精度を殆ど犠牲にすることなく、もっと高速に計算を行う方法を採用しても良い。
【0016】
すなわち、検索範囲の全体に亘り、検索ステップ毎に全て検波演算を実行するのではなく、検索範囲を例えば複数等分し、各検索分割領域毎に1回ずつ検波演算を実行し、その結果得られたシンボル点のデータを利用して、各検索分割領域内の残りの検索ステップの計算を簡単化することができる。一例として、最初に計算した見積搬送波周波数を中心として20の検索分割領域を設定し、各検索分割領域毎に80回の検索ステップを実行する場合を説明する。この場合、ステップ104の検波演算を実行するのは、全体で20回(各検索分割領域毎に1回)だけであり、各検索分割領域毎に残りの79回は、近似計算を行う。この場合の処理手順を図3の流れ図に示す。
【0017】
図3の流れ図は、図1のステップ108〜110までの処理に替わるものであり、計算負荷を軽減するための処理である。なお、図2の処理は、20の分割領域の各々における処理を示すものであり、検索範囲の全体では、図3の処理を検索範囲を変化させつつ20回繰り返すことになる。
【0018】
図2において、図1のステップ104で求めた直交座標データを極座標データ(r,φ)に座標変換する(ステップ200)。次に、極座標データから直交座標データのシンボル点に対応するデータ点を抽出する(ステップ202)。この処理は、単純に直交座標データのシンボル点と序列が同じ位置にあるデータ点を極座標データから抽出するだけなので、従来のように、面倒な検波演算及びシンボル点の探索処理を省略した極めて簡単な計算である。次に、極座標データから抽出したシンボル点に対応するデータ点を複数の分割領域で構成された座標平面上にプロットする(ステップ204)。次に、シンボル点及びシンボル点に対応するデータ点がプロットされなかった分割領域の数を計数する(ステップ206)。更に、極座標データ(r,φ)の偏角成分φのみを所定角度だけ変化させた極座標データ(r,φ+Δφ)を生成する(ステップ208)。この処理は、図1の場合には面倒な検波演算を実行して直交座標データを生成するステップ104の処理に対応するものであり、極座標データ(r,φ)の動径成分rを変化させず、偏角成分φのみを微小変化させるだけで、検索精度を落とさずに計算処理の負荷を大幅に軽減している。
【0019】
以上のステップ200〜208の処理は、判断ステップ210により、1つの検索分割領域の処理が終了するまで繰り返し実行される。すなわち、この例では、1つの検索分割領域の処理が80回繰り返されることになる。検索分割領域の処理は20回繰り返されるので、全体では80×20=1600回の処理が実行される。この処理の中で、図1の検波演算(ステップ104)が実行されるのは20回のみであり、図3の処理が79×20=1580回繰り返される。この結果、計算負荷が格段に軽減し、処理の高速化が実現できる。なお、上述の計算処理の実行回数は任意であり、この事例に限定されないことは勿論である。
【0020】
図5〜図8は、本発明により得られたコンピュータ16の表示スクリーン18上における表示結果のハードコピー出力例を示している。「I−Q」と表示された座標平面表示において、濃度の高い(濃い)点の集合部分がシンボル点又はシンボル点に対応するデータ点を表している。シンボル点以外の線状の表示は、シンボル点以外のデータに基づく表示であり、位相の変化の様子を表しているが、デジタル位相変調方式の判定のためにはシンボル点又はシンボル点に対応するデータ点の表示があれば十分である。図5では、「QPSK」のシンボル・パターン表示が明瞭に示され、図6では、「16QAM」のシンボル・パターンが明瞭に表示され、図7では「64QAM」のシンボル・パターン表示が明瞭に示され、図8では、「1/πシフトQPSK」のシンボル・パターン表示が明瞭に示されている。以上の表示から、当業者は一目瞭然に未知の入力信号のデジタル位相変調方式を判定することが可能である。
【0021】
以上、本発明の好適実施例を説明したが、本発明は、上述の実施例のみに限定されるものではなく、本発明の要旨から逸脱することなく、種々の変形及び修正を加え得ることは当業者には明らかである。
【図面の簡単な説明】
【図1】本発明の方法を示す流れ図の一部分である。
【図2】本発明の方法を示す図1の部分以外の流れ図である。
【図3】本発明の方法の他の実施の形態の一例を示す流れ図である。
【図4】本発明を適用するのに好適なシステムの構成を示すブロック図である。
【図5】本発明に基づくコンピュータ・ディスプレイ表示の一例を示す中間調画像である。
【図6】本発明に基づくコンピュータ・ディスプレイ表示の他の例を示す中間調画像である。
【図7】本発明に基づくコンピュータ・ディスプレイ表示の他の例を示す中間調画像である。
【図8】本発明に基づくコンピュータ・ディスプレイ表示の更に別の例を示す中間調画像である。
【符号の説明】
10 アンテナ
12 ダウンコンバータ
14 デジタル・スペクトラム・アナライザ
16 コンピュータ
18 表示スクリーン[0001]
TECHNICAL FIELD OF THE INVENTION
The present invention relates to a digital phase modulation method determination display method for determining a digital phase modulation method of a signal received from an unknown wireless device.
[0002]
Problems to be solved by the prior art and the invention
Various digital phase modulation methods such as PSK, QPSK, π / 4 shift QPSK, 16 QAM, and 256 QAM are available as digital phase modulation methods for wireless devices. In such a digital phase modulation method, if the carrier frequency and phase of the received signal are known, the received signal is synchronously detected at the accurate frequency and phase, and the detected data is easily confirmed by displaying the data graphically. be able to. That is, a unique symbol pattern corresponding to a phase change point is determined according to each digital phase modulation method. As described above, in the conventional signal analysis system, an appropriate synchronous detection can be performed based on a precise carrier frequency and phase by analyzing an input signal whose digital phase modulation method is known.
[0003]
However, in the case of an input signal from an unknown radio, the digital phase modulation method is unknown, and the carrier frequency and phase of the signal are also unknown, so that the input signal cannot be properly synchronously detected. In order to extract a carrier signal from an input signal, the digital phase modulation method must be known, and for an unknown digital phase modulation input signal, a large number of circuits corresponding to all conceivable methods respectively. There is only a method of switching and using until a proper synchronous detection can be performed, but such a trial-and-error analysis method has a drawback that it is extremely complicated and increases the cost. Therefore, it is difficult to accurately detect an input signal from an unknown wireless device by synchronous detection and display the graph, and it is also possible to determine the digital phase modulation method of the unknown wireless device from a received signal from the unknown wireless device. It was extremely difficult.
[0004]
SUMMARY OF THE INVENTION Accordingly, an object of the present invention is to provide a digital phase modulation method determination and display method capable of accurately and easily determining a digital phase modulation method of an input signal from an unknown wireless device.
[0005]
[Means for solving the problem]
The present invention provides a digital phase modulation method determination display method for automatically determining a digital phase modulation method of an input signal from an unknown wireless device. As a premise, when the in-phase component (I component) and the quadrature component (Q component) of the received signal are displayed on a quadrature coordinate plane or a polar coordinate plane, the fact that the symbol pattern differs depending on the digital phase modulation method is used. . Therefore, in order to perform display for determining a digital phase modulation method from an input signal from an unknown wireless device, the following procedure is employed. First, frequency domain data of an unknown digital phase modulation signal is generated, and an estimated carrier frequency of the unknown digital phase modulation signal is obtained from the frequency domain data. Next, the time domain data of the unknown digital phase modulation signal is detected and calculated based on the estimated carrier frequency to obtain rectangular coordinate data in the time domain. At each point of the rectangular coordinate data, one of the coordinates is an extreme value. Are determined, and only these symbol points are plotted on a coordinate plane composed of a plurality of divided areas, and the number of divided areas in which the symbol points are not plotted among the plurality of divided areas is counted. Are sequentially repeated while changing the estimated carrier frequency. As a result, the estimated carrier frequency when the number of divided regions in which the symbol points are not plotted becomes the maximum is obtained as the true carrier frequency. Based on this carrier frequency, time domain data of the unknown digital phase modulation signal is detected and calculated, and the resulting data is displayed on a display screen.
[0006]
According to the present invention, since the detection calculation for obtaining the rectangular coordinate data in the time domain has a relatively heavy calculation load, a method of reducing the calculation load can also be adopted. That is, not all the symbol points are obtained from the rectangular coordinate data, but the rectangular coordinate data is converted into polar coordinate data, and data corresponding to the previously obtained symbol point is extracted from the polar coordinate data. As a result, the calculation load for directly finding the symbol points from the orthogonal coordinate data is significantly reduced, and the processing speed is increased.
[0007]
BEST MODE FOR CARRYING OUT THE INVENTION
FIG. 4 is a block diagram showing an example of the configuration of an analysis system suitable for using the digital phase modulation system determination display method of the present invention. The
[0008]
The data analysis is executed by the
[0009]
Next, detection processing is performed on the time domain data of the input signal based on the estimated carrier frequency to obtain rectangular coordinate data in the time domain (step 104). As the time domain data of the input signal used in this case, the data used in the digital spectrum analyzer 14 may be used as it is, or the frequency domain data generated in
[0010]
(Equation 1)
a = cosωx
b = sin ωx
ω = 2πFT
F: Difference between the center frequency of the frequency domain width and the estimated carrier frequency T: Sampling period x: Data variable
I (x) = a · i (x) + b · q (x)
Q (x) = a · q (x) −bi · (x)
i (x): In-phase component data in time domain q (x): Quadrature component data in time domain I (x): I-axis component data of quadrature coordinate data Q (x): Q-axis component data of quadrature coordinate data By performing the detection operation, the rectangular coordinate data I (x) and Q (x) in the time domain are obtained. Note that the orthogonal coordinate data obtained here is not obtained by accurate synchronous detection, but is the result of an approximate estimation calculation.
[0011]
Next, for each point of the orthogonal coordinate data, a symbol point at which any one of the coordinates has an extreme value is determined (step 106). The symbol point is a point considered to correspond to the phase change point in the phase modulation. At the phase change point, one of the in-phase component coordinate I (x) and the quadrature component data Q (x) is an extreme value (maximum value). Value or local minimum value).
[0012]
Next, the symbol points obtained in
[0013]
Next, the number of divided regions in which no symbol point is plotted among a plurality of divided regions constituting the coordinate plane is counted (step 110). Thereafter, the estimated carrier frequency is changed by a predetermined value (step 112). The
[0014]
The search range of the estimated carrier frequency is a predetermined frequency range centered on the first estimated carrier frequency. For example, when the frequency domain data width is 200 kHz, the search range is 6.25 kHz. The search step of the estimated carrier frequency in the search process is, for example, about 4 Hz. These values are merely examples, and can be executed in other search ranges and search steps.
[0015]
In the above-described processing, the processing of the detection calculation performed in
[0016]
In other words, the detection operation is not performed for each search step over the entire search range, but the search range is divided into a plurality of equal parts, for example, and the detection calculation is performed once for each search divided region, and the result is obtained. Using the data of the obtained symbol points, the calculation of the remaining search steps in each search divided area can be simplified. As an example, a case will be described in which 20 search division regions are set around the estimated carrier wave frequency calculated first, and 80 search steps are executed for each search division region. In this case, the detection calculation in
[0017]
The flowchart of FIG. 3 replaces the processing of
[0018]
In FIG. 2, the orthogonal coordinate data obtained in
[0019]
The processing of the
[0020]
5 to 8 show examples of hard copy output of the display result on the
[0021]
As described above, the preferred embodiments of the present invention have been described. However, the present invention is not limited to only the above-described embodiments, and various changes and modifications can be made without departing from the gist of the present invention. It will be clear to those skilled in the art.
[Brief description of the drawings]
FIG. 1 is a part of a flow chart illustrating the method of the present invention.
FIG. 2 is a flowchart other than the part of FIG. 1 illustrating the method of the present invention.
FIG. 3 is a flowchart showing an example of another embodiment of the method of the present invention.
FIG. 4 is a block diagram showing a configuration of a system suitable for applying the present invention.
FIG. 5 is a halftone image showing an example of a computer display display according to the present invention.
FIG. 6 is a halftone image showing another example of a computer display display according to the present invention.
FIG. 7 is a halftone image showing another example of a computer display display according to the present invention.
FIG. 8 is a halftone image showing yet another example of a computer display display according to the present invention.
[Explanation of symbols]
Claims (8)
(b)該周波数領域データから上記未知のデジタル位相変調信号の見積搬送波周波数を求め、
(c)該見積搬送波周波数に基づいて未知のデジタル位相変調信号の時間領域データを検波演算し、時間領域の直交座標データを求め、
(d)該直交座標データの各点において何れか一方の座標が極値となるシンボル点を求め、
(e)これらのシンボル点のみを複数の分割領域で構成された座標平面上にプロットし、
(f)上記複数の分割領域の中で上記シンボル点がプロットされなかった分割領域の数を計数し、
(g)上記見積搬送波周波数を変化させ、
(h)上記ステップ(c)〜(g)を繰り返し、上記シンボル点がプロットされなかった分割領域の数が最大となるときの上記見積搬送波周波数に基づいて、上記未知のデジタル位相変調信号の時間領域データを検波演算し、
(i)この結果得られたデータを表示スクリーン上に表示することを特徴とするデジタル位相変調方式判定表示方法。(A) generating frequency domain data of an unknown digital phase modulation signal,
(B) determining an estimated carrier frequency of the unknown digital phase modulation signal from the frequency domain data;
(C) detecting and calculating time domain data of the unknown digital phase modulation signal based on the estimated carrier frequency to obtain rectangular coordinate data in the time domain;
(D) finding a symbol point at which one of the coordinates is an extreme value at each point of the rectangular coordinate data;
(E) plotting only these symbol points on a coordinate plane composed of a plurality of divided areas,
(F) counting the number of divided areas in which the symbol points are not plotted among the plurality of divided areas;
(G) changing the estimated carrier frequency,
(H) repeating the above steps (c) to (g) to determine the time of the unknown digital phase modulation signal based on the estimated carrier frequency when the number of divided areas where the symbol points are not plotted is maximized. Detect and calculate the area data,
(I) A digital phase modulation method determination display method, wherein the data obtained as a result is displayed on a display screen.
該時間領域データをFFT演算することにより上記周波数領域データを生成することを特徴とする請求項1記載のデジタル位相変調方式判定表示方法。In the step (a), time domain data of the unknown digital phase modulation signal is generated,
2. The method according to claim 1, wherein the frequency domain data is generated by performing an FFT operation on the time domain data.
(b)該周波数領域データから上記未知のデジタル位相変調信号の見積搬送波周波数を求め、
(c)該見積搬送波周波数に基づいて未知のデジタル位相変調信号の時間領域データを検波演算し、時間領域の直交座標データを求め、
(d)該直交座標データの各点において何れか一方の座標が極値となるシンボル点を求め、
(e)上記直交座標データを極座標データに座標変換し、
(f)該極座標データから上記シンボル点に対応するデータ点を抽出し、
(g)これらの抽出したデータ点のみを複数の分割領域で構成された極座標平面上にプロットし、
(h)上記複数の分割領域の中でデータ点がプロットされなかった分割領域の数を計数し、
(i)上記見積搬送波周波数を変化させ、
(j)上記ステップ(c)〜(i)を繰り返し、上記データ点がプロットされなかった分割領域の数が最大となるときの上記見積搬送波周波数に基づいて、上記未知のデジタル位相変調信号の時間領域データを検波演算し、
(k)この結果得られたデータを表示スクリーン上に表示することを特徴とするデジタル位相変調方式判定表示方法。(A) generating frequency domain data of an unknown digital phase modulation signal,
(B) determining an estimated carrier frequency of the unknown digital phase modulation signal from the frequency domain data;
(C) detecting and calculating time domain data of the unknown digital phase modulation signal based on the estimated carrier frequency to obtain rectangular coordinate data in the time domain;
(D) finding a symbol point at which one of the coordinates is an extreme value at each point of the rectangular coordinate data;
(E) converting the rectangular coordinate data into polar coordinate data,
(F) extracting a data point corresponding to the symbol point from the polar coordinate data;
(G) plotting only these extracted data points on a polar coordinate plane composed of a plurality of divided areas,
(H) counting the number of divided areas in which data points are not plotted among the plurality of divided areas;
(I) changing the estimated carrier frequency,
(J) repeating the above steps (c) to (i) to determine the time of the unknown digital phase modulation signal based on the estimated carrier frequency when the number of divided areas where the data points are not plotted is at a maximum. Detect and calculate the area data,
(K) A digital phase modulation method determination display method, wherein the data obtained as a result is displayed on a display screen.
該時間領域データをFFT演算することにより上記周波数領域データを生成することを特徴とする請求項5記載のデジタル位相変調方式判定表示方法。In the step (a), time domain data of the unknown digital phase modulation signal is generated,
6. The method according to claim 5, wherein the frequency domain data is generated by performing an FFT operation on the time domain data.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP13281796A JP3605229B2 (en) | 1996-04-30 | 1996-04-30 | Digital phase modulation method display method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP13281796A JP3605229B2 (en) | 1996-04-30 | 1996-04-30 | Digital phase modulation method display method |
Publications (2)
Publication Number | Publication Date |
---|---|
JPH09298567A JPH09298567A (en) | 1997-11-18 |
JP3605229B2 true JP3605229B2 (en) | 2004-12-22 |
Family
ID=15090269
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP13281796A Expired - Fee Related JP3605229B2 (en) | 1996-04-30 | 1996-04-30 | Digital phase modulation method display method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP3605229B2 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3934607B2 (en) * | 2001-07-06 | 2007-06-20 | 株式会社アドバンテスト | MEASUREMENT DATA DISPLAY DEVICE, METHOD, PROGRAM, AND RECORDING MEDIUM CONTAINING THE PROGRAM |
US8615382B2 (en) * | 2011-01-27 | 2013-12-24 | Tektronix, Inc. | Test and measurement instrument with common presentation of time domain data |
-
1996
- 1996-04-30 JP JP13281796A patent/JP3605229B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JPH09298567A (en) | 1997-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2086704C (en) | Method for evaluating similarity of signals having a carrier frequency offset | |
EP2048509B1 (en) | Modulation signature trigger | |
US20090154589A1 (en) | Systems and methods for signal modulation and demodulation using phase | |
JP2004007692A (en) | Frequency offset estimator | |
JP3271504B2 (en) | Frequency estimation circuit and AFC circuit using the same | |
USRE49560E1 (en) | Method and a measuring device for measuring broadband measurement signals | |
CN107305223B (en) | Improved phase difference frequency estimation method | |
CN109067680B (en) | Carrier frequency offset estimation method and device of baseband signal | |
JP2022021603A (en) | Receiving device and receiving method, and mobile terminal test apparatus provided with the receiving device | |
JP3605229B2 (en) | Digital phase modulation method display method | |
JP2020010195A (en) | Frequency estimation device and tracking receiver | |
US7697634B2 (en) | Interpolation of complex signals | |
JP3950242B2 (en) | Offset QPSK modulation analysis method | |
JP5155381B2 (en) | Signal generator, signal generation system, and signal generation method | |
US9537690B1 (en) | Method and apparatus for extraction of baseband waveform from amplitude modulated signal via time domain sampling | |
JP3819592B2 (en) | 64QAM, 256QAM modulation analysis method | |
JP6586183B2 (en) | Frame synchronization apparatus, measurement apparatus including the same, frame synchronization method, and measurement method | |
JPH07245634A (en) | Method and system for description of characteristic of phase-modulated signal | |
JPH11234354A (en) | Measurement method for modulation accuracy of dpsk modulator | |
JP2003244263A (en) | Signal processor | |
JP4512835B2 (en) | Constellation display method and apparatus | |
JP4105079B2 (en) | Phase difference calculation circuit | |
JP2833877B2 (en) | Modulation error measurement device | |
JP3360583B2 (en) | Phase modulation signal detection device | |
CN113162652A (en) | Method and module for detecting frequency of spread spectrum signal |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20040914 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20041001 |
|
R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20071008 Year of fee payment: 3 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20081008 Year of fee payment: 4 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091008 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20091008 Year of fee payment: 5 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20101008 Year of fee payment: 6 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20111008 Year of fee payment: 7 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121008 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20121008 Year of fee payment: 8 |
|
FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20131008 Year of fee payment: 9 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
LAPS | Cancellation because of no payment of annual fees |