JP2007336078A - 変調波識別装置、変調波識別方法、及び受信装置 - Google Patents
変調波識別装置、変調波識別方法、及び受信装置 Download PDFInfo
- Publication number
- JP2007336078A JP2007336078A JP2006163666A JP2006163666A JP2007336078A JP 2007336078 A JP2007336078 A JP 2007336078A JP 2006163666 A JP2006163666 A JP 2006163666A JP 2006163666 A JP2006163666 A JP 2006163666A JP 2007336078 A JP2007336078 A JP 2007336078A
- Authority
- JP
- Japan
- Prior art keywords
- probability density
- density distribution
- distribution
- modulated wave
- calculated
- 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
Links
Images
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
- Circuits Of Receivers In General (AREA)
Abstract
【課題】変調波がマルチキャリアであるか否かを簡易に識別する。
【解決手段】変調波から抽出されたベースバンド信号が入力されると(ステップS21)、このベースバンド信号の時間軸波形をサンプリングし、このサンプリング値から4次のキュムラントを計算する(ステップS22)。また、ベースバンド信号の振幅分布を示す確率密度分布を求め、この確率密度分布と、基準とする確率密度分布との間のクルバック・ライブラーの距離を計算する(ステップS23)。そして、4次のキュムラントの計算値により前記サンプリング値がガウス分布に近いと判定され(ステップS24)、かつクルバック・ライブラーの距離の計算値に基づいて上記求められた確率密度分布が基準とする確率密度分布と類似していると判定される(ステップS25)場合に、変調波をマルチキャリアと識別する(ステップS26)。
【選択図】 図2
【解決手段】変調波から抽出されたベースバンド信号が入力されると(ステップS21)、このベースバンド信号の時間軸波形をサンプリングし、このサンプリング値から4次のキュムラントを計算する(ステップS22)。また、ベースバンド信号の振幅分布を示す確率密度分布を求め、この確率密度分布と、基準とする確率密度分布との間のクルバック・ライブラーの距離を計算する(ステップS23)。そして、4次のキュムラントの計算値により前記サンプリング値がガウス分布に近いと判定され(ステップS24)、かつクルバック・ライブラーの距離の計算値に基づいて上記求められた確率密度分布が基準とする確率密度分布と類似していると判定される(ステップS25)場合に、変調波をマルチキャリアと識別する(ステップS26)。
【選択図】 図2
Description
この発明は、未知の変調波を受信し、この変調波の変調方式を自動的に識別する変調波識別装置、変調波識別方法、及び受信装置に関する。
通信の分野において、変調方式の識別は、目的とする信号の自動受信や自動監視等に幅広く応用されている。変調方式の識別処理には、例えば、周波数スペクトラムの形状から各変調方式の有する特徴を抽出し、この特徴を用いて変調方式を識別する手法が提案されている。これにより、通信諸元が未知である複数の変調方式の受信信号を自動的に識別し、その変調方式を特定できるため、単一装置にて、複数の変調方式の識別が可能となり、装置の小型化、及び高性能化を実現が可能となる(特許文献1を参照。)。
特許第3576880号公報
ところが、実際に伝搬されている変調波の周波数スペクトラムは、時間的に変動しており、スペクトラムの形状から変調波がマルチキャリアであるかシングルキャリアであるのかを識別するのは困難である。このため、マルチキャリアであるという仮定のもとに所定の周波数間隔でサブキャリアに分離してから分析する必要があり、変調波が目的とする信号であるか否かが分離後でないと判断できない。
この発明は上記事情に着目してなされたもので、その目的とするところは、変調波がマルチキャリアであるか否かを簡易に識別することができる変調波識別装置、変調波識別方法、及び受信装置を提供することにある。
上記目的を達成するためにこの発明に係わる変調波識別装置は、受信された変調波のベースバンド信号を入力して時間軸波形をサンプリングし、このサンプリング値から4次のキュムラント(Cumlant)を計算する第1の計算手段と、前記ベースバンド信号を入力して振幅分布を示す確率密度分布を求め、この確率密度分布と基準とする確率密度分布との間のクルバック・ライブラーの距離(Kullback-Leibler Divergence)を計算する第2の計算手段と、前記第1の計算手段により計算された4次のキュムラントの計算値と所定の第1の閾値とを比較し前記サンプリング値がガウス分布に近似されるか否かを判定する第1の判定手段と、前記第2の計算手段により計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定手段と、第1の判定手段において前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定手段において前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別手段とを具備することを特徴とする。
この発明の変調波識別方法は、変調波を受信する変調波識別装置で用いられ、前記受信された変調波からベースバンド信号を抽出する抽出ステップと、前記抽出されたベースバンド信号の時間軸波形をサンプリングし、このサンプリング値から4次のキュムラントを計算する第1の計算ステップと、前記抽出されたベースバンド信号の振幅分布を示す確率密度分布を求め、この確率密度分布と、基準とする確率密度分布との間のクルバック・ライブラーの距離を計算する第2の計算ステップと、前記第1の計算ステップにより計算された4次のキュムラントの計算値と所定の第1の閾値とを比較し前記サンプリング値がガウス分布に近似されるか否かを判定する第1の判定ステップと、前記第2の計算ステップにより計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定ステップと、第1の判定ステップにおいて前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定ステップにおいて前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別ステップとを具備することを特徴とする。
また、この発明の受信装置は、変調波を受信する受信手段と、前記受信された変調波からベースバンド信号を抽出する抽出手段と、前記抽出されたベースバンド信号の時間軸波形をサンプリングし、このサンプリング値から4次のキュムラントを計算する第1の計算手段と、前記抽出されたベースバンド信号の振幅分布を示す確率密度分布を求め、この確率密度分布と、基準とする確率密度分布との間のクルバック・ライブラーの距離を計算する第2の計算手段と、前記第1の計算手段により計算された4次のキュムラントの計算値と所定の第1の閾値とを比較し前記サンプリング値がガウス分布に近似されるか否かを判定する第1の判定手段と、
前記第2の計算手段により計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定手段と、第1の判定手段において前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定手段において前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別手段とを具備することを特徴とする。
前記第2の計算手段により計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定手段と、第1の判定手段において前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定手段において前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別手段とを具備することを特徴とする。
マルチキャリアの信号は、複数のキャリアに異なる内容の信号が含まれることが特徴である。このため、信号全体の確率密度分布は複素ガウス分布に近くなる。この性質を用いて、受信信号全体がガウス分布に近いことを検出しマルチキャリアと判断する。これにより、受信信号をサブキャリアに分離することなく、入力された変調波がマルチキャリアであるか否かを簡易な方法で判断することが可能となる。
したがってこの発明によれば、変調波がマルチキャリアであるか否かを簡易に識別することができる変調波識別装置、変調波識別方法、及び受信装置を提供することができる。
以下、図面を参照しながら本発明の実施の形態を詳細に説明する。
図1は、この発明に係わる変調波識別装置を備えた受信装置の一実施形態の構成を示すブロック図である。この受信装置100は、アンテナ1と、受信部2と、A/D変換部3と、FFT処理部4と、受信フィルタ5と、変調波識別部6と、分離部7と、復調部8とを備える。
図1は、この発明に係わる変調波識別装置を備えた受信装置の一実施形態の構成を示すブロック図である。この受信装置100は、アンテナ1と、受信部2と、A/D変換部3と、FFT処理部4と、受信フィルタ5と、変調波識別部6と、分離部7と、復調部8とを備える。
アンテナ1に到来した変調信号は、受信部2により受信信号として取り込まれ、A/D変換部3を介してデジタル信号に変換される。変換されたデジタル信号は、FFT処理部4により周波数軸上のスペクトラムに変換され、この周波数スペクトラムから目的とする信号の出現を検知する。検知においては目的の信号の持続時間や帯域幅を考慮して不要な信号を排除する。そして、検知された目的とする信号から中心周波数を特定し、中心周波数に基づいて受信フィルタ5より目的の信号のベースバンド信号を抽出する。
さて、抽出されたベースバンド信号は、変調波識別部6に入力される。図2は、変調波識別部6の処理手順とその内容を示すフローチャートである。
変調波識別部6は、ステップS21において受信フィルタ5からベースバンド信号が入力されると、まず、ステップS22において4次のキュムラントを計算する。図3は、4次のキュムラントの計算式を示す図である。ベースバンド信号の時間軸波形をサンプリングし、このサンプリングした振幅の値を入力値として4次のキュムラントCum4を計算する。なお、E1〜E4は、期待値の計算を示し、xは入力信号を示す変数である。4次のキュムラントの値は、信号のガウス性が高いほど0に近づく性質を有する。
変調波識別部6は、ステップS21において受信フィルタ5からベースバンド信号が入力されると、まず、ステップS22において4次のキュムラントを計算する。図3は、4次のキュムラントの計算式を示す図である。ベースバンド信号の時間軸波形をサンプリングし、このサンプリングした振幅の値を入力値として4次のキュムラントCum4を計算する。なお、E1〜E4は、期待値の計算を示し、xは入力信号を示す変数である。4次のキュムラントの値は、信号のガウス性が高いほど0に近づく性質を有する。
また、変調波識別部6は、ステップS23においてクルバック・ライブラーの距離を計算する。図4は、クルバック・ライブラーの距離を説明する図である。クルバック・ライブラーの距離Lは、2つの確率分布の距離の尺度であり、類似度を測ることができる。ここで、p(x)は基準となる確率密度関数、q(x)は入力信号の振幅分布を示す確率密度関数である。ここでは、基準とする確率密度関数には、複素ガウス分布を振幅分布に変換した場合の分布であるレイリー分布を使用する。つまり、Lの値が0に近づくほど入力信号の振幅分布が基準分布であるレイリー分布に近いと言える。
ステップS24において、変調波識別部6は、ステップS22で求められた4次のキュムラントの値と予め定められた閾値とを比較し、サンプリング値がガウス分布に近似されるか否かを判定する。また、変調波識別部6は、ステップS25において、ステップS23で求められたクルバック・ライブラーの距離の計算値と予め定められた閾値とを比較し、入力信号の振幅分布がレイリー分布に類似しているか否かを判定する。
変調波識別部6は、ステップS26において、ステップS22においてサンプリング値がガウス分布に近似されると判定され、かつステップS23において入力信号の振幅分布がレイリー分布に類似していると判定された場合に、受信信号をマルチキャリアと判定する。なお、各閾値はマルチキャリア以外(シングルキャリア)の場合の値の分布を求めた後、マルチキャリアと境界線が引ける条件を用いて決定する。
図5に、4次のキュムラント及びクルバック・ライブラーの距離の計算結果の例を示す。図5の(a)及び(b)は、シングルキャリアの場合の計算例である。図5の(c)は、マルチキャリア(キャリア数8)の場合の計算例である。同図を比較しても明らかなように、4次のキュムラント及びクルバック・ライブラーの距離の計算値から、受信信号がマルチキャリアであるか否かを簡単に判別することができる。これにより、受信信号がシングルキャリアと分かればサブキャリアに分離する処理を省略することができる。
変調波識別部6において、変調方式がマルチキャリアであると判定された場合は、受信信号は分離部7においてサブキャリアにそれぞれ分離される。分離された各サブキャリアは復調部8によりそれぞれ復調され、受信データとして出力される。一方、変調波識別部6においてシングルキャリアと判定された場合には、受信信号は復調部8に入力され、復調処理されのち受信データとして出力される。
図6に受信信号の周波数スペクトルの一例を示す。この図からわかるように、実際の受信信号の周波数スペクトラムは、スペクトラムの形状から受信信号がマルチキャリアであるかシングルキャリアであるのかを識別するのは困難である。このため、受信信号をマルチキャリアであるという仮定のもとに所定の周波数間隔でサブキャリアに分離してから分析する必要があり、受信信号が所望の信号であるか否かが分離後でないと分からなかった。
これに対し、上記実施形態では、マルチキャリアの信号は、複数のキャリアに異なる信号が含まれるため、信号全体の確率密度分布が複素ガウス分布に近くなるという特徴に注目し、受信信号全体がガウス分布に近いことを検出してマルチキャリアと判断する。これにより、受信信号をサブキャリアに分離することなく、入力された変調波がマルチキャリアであるか否かを簡易な方法で判断することが可能となる。これにより、受信された変調波が目的とする信号であるか否かを簡単に識別することが可能となる。
なお、この発明は上記実施形態に限定されるものではない。例えば、上記実施形態において、4次のキュムラントによる判定とクルバック・ライブラーの距離による判定のいずれかに重み付けをしてマルチキャリアであるか否かを判定するようにしてもよい。また、目的とする信号の経験値から各閾値を変更することもできる。その他にも、受信装置100の構成及び動作手順とその内容についてもこの発明の要旨を逸脱しない範囲で種々に変形して実施できる。
要するにこの発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。
100…受信装置、1…アンテナ、2…受信部、3…A/D変換部、4…FFT処理部、5…受信フィルタ、6…変調波識別部、7…分離部、8…復調部。
Claims (6)
- 受信された変調波のベースバンド信号を入力して時間軸波形をサンプリングし、このサンプリング値から4次のキュムラント(Cumlant)を計算する第1の計算手段と、
前記ベースバンド信号を入力して振幅分布を示す確率密度分布を求め、この確率密度分布と基準とする確率密度分布との間のクルバック・ライブラーの距離(Kullback-Leibler Divergence)を計算する第2の計算手段と、
前記第1の計算手段により計算された4次のキュムラントの計算値と所定の第1の閾値とを比較し前記サンプリング値がガウス分布に近似されるか否かを判定する第1の判定手段と、
前記第2の計算手段により計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定手段と、
第1の判定手段において前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定手段において前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別手段と
を具備することを特徴とする変調波識別装置。 - 前記第2の計算手段は、前記基準とする確率密度分布としてレイリー分布を用いることを特徴とする請求項1記載の変調波識別装置。
- 変調波を受信する変調波識別装置で用いられ、
前記受信された変調波からベースバンド信号を抽出する抽出ステップと、
前記抽出されたベースバンド信号の時間軸波形をサンプリングし、このサンプリング値から4次のキュムラントを計算する第1の計算ステップと、
前記抽出されたベースバンド信号の振幅分布を示す確率密度分布を求め、この確率密度分布と、基準とする確率密度分布との間のクルバック・ライブラーの距離を計算する第2の計算ステップと、
前記第1の計算ステップにより計算された4次のキュムラントの計算値と所定の第1の閾値とを比較し前記サンプリング値がガウス分布に近似されるか否かを判定する第1の判定ステップと、
前記第2の計算ステップにより計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定ステップと、
第1の判定ステップにおいて前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定ステップにおいて前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別ステップと
を具備することを特徴とする変調波識別方法。 - 前記第2の計算ステップは、前記基準とする確率密度分布としてレイリー分布を用いることを特徴とする請求項3記載の変調波識別方法。
- 変調波を受信する受信手段と、
前記受信された変調波からベースバンド信号を抽出する抽出手段と、
前記抽出されたベースバンド信号の時間軸波形をサンプリングし、このサンプリング値から4次のキュムラントを計算する第1の計算手段と、
前記抽出されたベースバンド信号の振幅分布を示す確率密度分布を求め、この確率密度分布と、基準とする確率密度分布との間のクルバック・ライブラーの距離を計算する第2の計算手段と、
前記第1の計算手段により計算された4次のキュムラントの計算値と所定の第1の閾値とを比較し前記サンプリング値がガウス分布に近似されるか否かを判定する第1の判定手段と、
前記第2の計算手段により計算されたクルバック・ライブラーの距離の計算値と所定の第2の閾値とを比較し前記求められた確率密度分布が前記基準とする確率密度分布に類似しているか否かを判定する第2の判定手段と、
第1の判定手段において前記サンプリング値がガウス分布に近似されると判定され、かつ第2の判定手段において前記求められた確率密度分布が前記基準とする確率密度分布に類似していると判定される場合に、前記変調波がマルチキャリアであると識別する識別手段と
を具備することを特徴とする受信装置。 - 前記第2の計算手段は、前記基準とする確率密度分布としてレイリー分布を用いることを特徴とする請求項5記載の受信装置。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006163666A JP2007336078A (ja) | 2006-06-13 | 2006-06-13 | 変調波識別装置、変調波識別方法、及び受信装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2006163666A JP2007336078A (ja) | 2006-06-13 | 2006-06-13 | 変調波識別装置、変調波識別方法、及び受信装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2007336078A true JP2007336078A (ja) | 2007-12-27 |
Family
ID=38935148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2006163666A Pending JP2007336078A (ja) | 2006-06-13 | 2006-06-13 | 変調波識別装置、変調波識別方法、及び受信装置 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2007336078A (ja) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010521841A (ja) * | 2007-03-16 | 2010-06-24 | トムソン ライセンシング | マルチパス環境において変調を分類する装置及び方法 |
JP2013002878A (ja) * | 2011-06-14 | 2013-01-07 | Toshiba Corp | 受信装置及び背景雑音推定方法 |
JP2013187844A (ja) * | 2012-03-09 | 2013-09-19 | Toshiba Corp | 信号識別装置 |
US8693380B2 (en) | 2007-02-27 | 2014-04-08 | Thomson Licensing | Method and apparatus for power management in WLAN |
US20180063733A1 (en) * | 2016-08-26 | 2018-03-01 | Allen-Vanguard Corporation | System and Method for Cellular Network Identification |
US10027422B2 (en) | 2016-05-20 | 2018-07-17 | Fujitsu Limited | Optical signal receiving device and method of receiving an optical signal |
-
2006
- 2006-06-13 JP JP2006163666A patent/JP2007336078A/ja active Pending
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8693380B2 (en) | 2007-02-27 | 2014-04-08 | Thomson Licensing | Method and apparatus for power management in WLAN |
JP2010521841A (ja) * | 2007-03-16 | 2010-06-24 | トムソン ライセンシング | マルチパス環境において変調を分類する装置及び方法 |
US8385473B2 (en) | 2007-03-16 | 2013-02-26 | Thomson Licensing | Apparatus and method for classifying modulations in multipath environments |
JP2013002878A (ja) * | 2011-06-14 | 2013-01-07 | Toshiba Corp | 受信装置及び背景雑音推定方法 |
JP2013187844A (ja) * | 2012-03-09 | 2013-09-19 | Toshiba Corp | 信号識別装置 |
US9071469B2 (en) | 2012-03-09 | 2015-06-30 | Kabushiki Kaisha Toshiba | Signal identification device |
US10027422B2 (en) | 2016-05-20 | 2018-07-17 | Fujitsu Limited | Optical signal receiving device and method of receiving an optical signal |
US20180063733A1 (en) * | 2016-08-26 | 2018-03-01 | Allen-Vanguard Corporation | System and Method for Cellular Network Identification |
GB2556370A (en) * | 2016-08-26 | 2018-05-30 | Allen Vanguard Corp | System and method for cellular network identification |
GB2556370B (en) * | 2016-08-26 | 2019-12-11 | Allen Vanguard Corp | System and method for cellular network identification |
US10548037B2 (en) * | 2016-08-26 | 2020-01-28 | Allen-Vanguard Corporation | System and method for cellular network identification |
AU2017219069B2 (en) * | 2016-08-26 | 2021-12-23 | Allen-Vanguard Corporation | System and method for cellular network identification |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Yang et al. | Cyclostationary feature detection based spectrum sensing algorithm under complicated electromagnetic environment in cognitive radio networks | |
CN108632185B (zh) | 一种船舶vdes系统的ais体制解调方法及解调系统 | |
US8619909B2 (en) | Signal detector using matched filter for training signal detection | |
US9584347B2 (en) | Methods and systems for rapid detection of digital radio signals | |
JP2007336078A (ja) | 変調波識別装置、変調波識別方法、及び受信装置 | |
Hassan et al. | Automatic modulation recognition using wavelet transform and neural network | |
CN102780534B (zh) | 一种通用的电磁信号自动搜索方法 | |
CN105207965A (zh) | 一种vhf/uhf频段的调制自动识别方法 | |
CN112003803B (zh) | 一种vhf、uhf频段航空电台信号的侦收设备 | |
US9083569B2 (en) | Apparatus and method for detecting broadcasting signal | |
JP2016119644A (ja) | 無線受信装置、電子機器、及び無線受信方法 | |
CN115347978B (zh) | 一种识别ais帧头数据的方法、装置及存储介质 | |
CN103236863B (zh) | 采用了硬限幅器的fm数字调制解调电路的静噪处理方法 | |
CN108880720B (zh) | 收音设备自动搜台方法、装置及收音设备 | |
CN108055090B (zh) | 一种ais台站的时频功率稳定性检测方法 | |
CN108243131B (zh) | 解调方法、装置、频谱检测仪及计算机可读取存储介质 | |
JP2007189420A (ja) | 受信処理装置 | |
US11523366B2 (en) | Terminal number estimating system, terminal specifying system, terminal number estimating device, terminal number specifying device, and processing method | |
CN112865882A (zh) | 一种抗频偏的dmr对讲机信号快速识别方法 | |
CN110235377B (zh) | 解调器 | |
JP2008085894A (ja) | シンボル速度検出装置およびプログラム | |
KR101303256B1 (ko) | 모르스 신호의 실시간 탐지 해독 장치 및 방법 | |
KR102005616B1 (ko) | Ofdm 기반의 dab 시스템 및 그것을 이용한 모드 및 프레임 동기 방법 | |
CN113098817B (zh) | 一种基于Goertzel算法的标签信号解调方法与系统 | |
Li et al. | Automatic digital modulation recognition based on euclidean distance in hyperspace |