(57)【要約】【課題】異なる複数の周波数帯域のリアルタイム監視を行なう。【解決手段】（１）１回目からｎ回目までのサンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得するステップと、（２）１回目からｎ回目までのデータから取得された周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１倍の周波数までの周波数スペクトラムを生成するステップと、（３）上記（２）のステップにより得られた１回目からｎ回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成するステップと、を備える。【選択図】図９
2020106478.tif 000002
【特許請求の範囲】
【請求項１】
アンダーサンプリング受信方式を用いたリアルタイム周波数モニタにおいて、サンプリング周波数をｎ−１（ｎは２以上の整数）回変えてサンプリングしたｎ個の時間軸データから、受信周波数スペクトラムを再生する際に、
（１）１回目からｎ回目までのサンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得するステップと、
（２）１回目からｎ回目までのデータから取得された周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までの周波数スペクトラムを生成するステップと、
（３）前記（２）のステップにより得られた１回目からｎ回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成するステップと、
を備えることを特徴とする周波数スペクトラム再生方法。
【請求項２】
請求項１に記載の周波数スペクトラム再生方法において、
前記サンプリング周波数をｎ−１回変えてサンプリングする動作を繰り返し行なう際に、ｋ（ｋはｎ以上の整数）回目の繰り返し時の再生を、
（１ａ）ｋ回目の前記サンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得するステップと、
（２ａ）ｋ回目のデータから得られた周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１倍の周波数までの周波数スペクトラムを生成するステップと、
（３ａ）前記（２ａ）のステップにより得られたｋ回目からｋ−ｎ＋１回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成するステップと、
を備えることを特徴とする周波数スペクトラム再生方法。
【請求項３】
アンダーサンプリング受信方式を用いたリアルタイム周波数モニタにおいて、
（１）１回目からｎ（ｎは２以上の整数）回目までのサンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得する手段と、
（２）１回目からｎ回目までのデータから取得された周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までの周波数スペクトラムを生成する手段と、
（３）前記（２）の手段により得られた１回目からｎ回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成する手段と、
を備え、
サンプリング周波数をｎ−１回変えてサンプリングしたｎ個の時間軸データから、受信周波数スペクトラムを再生することを特徴とする周波数スペクトラム再生装置。
【請求項４】
請求項３に記載の周波数スペクトラム再生装置において、
前記サンプリング周波数をｎ−１回変えてサンプリングする動作を繰り返し行なう際に、
（１ａ）ｋ（ｋはｎ以上の整数）回目の前記サンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得する手段と、
（２ａ）ｋ回目のデータから得られた周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１倍の周波数までの周波数スペクトラムを生成する手段と、
（３ａ）前記（２）の手段により得られたｋ回目からｋ−ｎ＋１回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成する手段と、
を備え、
ｋ回目の繰り返し時の再生を行なうことを特徴とする周波数スペクトラム再生装置。
【請求項５】
低雑音増幅器と無線周波数帯のアンチエリアスフィルタとの後段に並列に接続され、ｎ（ｎは２以上の整数）個の異なるサンプリング周波数でアンダーサンプリングする受信回路と、
１個目からｎ個目までの前記受信回路で取得してサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までのｎ個の周波数スペクトラムを取得する取得部と、
１個目からｎ個目までの前記受信回路で取得したｎ個の周波数スペクトラムをそれぞれナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までのｎ個の周波数スペクトラムを生成する第１生成部と、
前記第１生成部により得られた１個目からｎ個目までの前記受信回路の各周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成する第２生成部と、
を備えることを特徴とする受信機。
【請求項６】
請求項５に記載の受信機において、
前記第２生成部によって生成された前記１つの周波数スペクトラムを表示装置に表示させることを特徴とする受信機。
【請求項７】
（１）低雑音増幅器と無線周波数帯のアンチエリアスフィルタとの後段に並列にｎ（ｎは２以上の整数）個接続された受信回路において、各々異なるサンプリング周波数でアンダーサンプリングするステップと、
（２）１個目からｎ個目までの前記受信回路で取得してサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までのｎ個の周波数スペクトラムを取得するステップと、
（３）１個目からｎ個目までの前記受信回路で取得したｎ個の周波数スペクトラムをそれぞれナイキスト周波数毎に折り返してナイキスト周波数の２＋１ｍ（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までのｎ個の周波数スペクトラムを生成するステップと、
（４）前記（３）のステップにより得られた１個目からｎ個目までの前記受信回路の各周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成するステップと、
を備えることを特徴とする受信方法。
【請求項８】
請求項７に記載の受信方法において、
前記（４）のステップによって生成された前記１つの周波数スペクトラムを表示装置に表示させることを特徴とする受信方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本明細書に記載する技術は、周波数スペクトラム再生方法、周波数スペクトラム再生装置、受信機及び受信方法に関する。
【背景技術】
【０００２】
近年、ワイヤレスInternet of Things（ＩｏＴ）の普及により、異なるワイヤレスＩｏＴシステム間の干渉を避けるための、ワイヤレスＩｏＴ周波数帯域リアルタイムスペクトルモニタリング技術が開発されている。ワイヤレスＩｏＴ周波数帯域は、例えば、９２０ＭＨｚ帯域，２．４ＧＨｚ帯域及び５〜６ＧＨｚ帯域である。
【０００３】
異なるワイヤレスＩｏＴシステム間の干渉を避けるために、クロックスイッチング方式を用いたワイヤレスＩｏＴリアルタイムスペクトルモニタ用ダイレクトRadio Frequency（ＲＦ）アンダーサンプリング受信機が提唱される。これにより、クロック周波数の１／２（別言すれば、「ナイキスト限界」）を上回る監視周波数帯域幅で、１ミリ秒単位未満のリアルタイム性能が実現される。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２０１８−１７４４１５号公報
【非特許文献】
【０００５】
【非特許文献１】板谷聡子・丸橋建一・長谷川  淳・長谷川晃朗・雨海明博・尾関  敦・小林  宰・中島健智・山田亮太・辻  聡・江連裕一郎・伊藤  睦・冨田尚孝・児島史秀著、電子情報通信学会技術研究報告vol. 116, no. 479「製造現場における多種無線通信〜製造現場における無線通信の課題〜」RCS2016-303, pp. 85-90、一般社団法人電子情報通信学会、２０１７年３月
【非特許文献２】古市朋之・本良瑞樹・亀田卓・末松憲治著、「無線IoTリアルタイム周波数センサ用  ダイレクトRFアンダーサンプリング受信機のクロック周波数の検討」、2018年  電子情報通信学会総合大会、C-2-96、２０１８年３月
【発明の概要】
【発明が解決しようとする課題】
【０００６】
工場等の屋内空間では、異なるＩｏＴシステム間の干渉や、製造装置で発生する電磁ノイズとＩｏＴシステムとの間の干渉が問題となる場合がある。異なるワイヤレスＩｏＴシステム間の干渉を避けるためには、ＩｏＴ周波数帯域の信号及びノイズをミリ秒単位で監視するリアルタイム周波数検知技術が用いられる。
【０００７】
しかしながら、１２ＧＨｚを上回る高速アナログ／デジタル変換器（ＡＤＣ）が必要になるため、これらの帯域をナイキストサンプリングで一度にリアルタイム監視することは容易でない。
【０００８】
１つの側面では、本明細書に記載する技術は、異なる複数の周波数帯域のリアルタイム監視を行なうことを目的とする。
【課題を解決するための手段】
【０００９】
１つの側面において、周波数スペクトラム再生方法は、アンダーサンプリング受信方式を用いたリアルタイム周波数モニタにおいて、サンプリング周波数をｎ−１（ｎは２以上の整数）回変えてサンプリングしたｎ個の時間軸データから、受信周波数スペクトラムを再生する際に、（１）１回目からｎ回目までのサンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得するステップと、（２）１回目からｎ回目までのデータから取得された周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までの周波数スペクトラムを生成するステップと、（３）前記（２）のステップにより得られた１回目からｎ回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成するステップと、を備える。
【発明の効果】
【００１０】
１つの側面として、異なる複数の周波数帯域のリアルタイム監視を行なうことができる。
【図面の簡単な説明】
【００１１】
【図１】実施形態としての受信機の構成例を示すブロック図である。
【図２】図１に示した受信機におけるアンダーサンプリング処理を説明する図である。
【図３】図１に示した受信機におけるサンプリングクロックの切り替え前のアンダーサンプリング受信動作を説明する図である。
【図４】図１に示した受信機におけるサンプリングクロックの切り替え後のアンダーサンプリング受信動作を説明する図である。
【図５】図１に示した受信機における測定周波数帯域とサンプリング周波数との関係を例示するテーブルである。
【図６】図１に示した受信機におけるアンダーサンプリング受信時の観測周波数の計算結果を例示するグラフである。
【図７】図１に示した受信機におけるクロック切り替えを用いた周波数同定法を説明する図である。
【図８】図１に示した受信機におけるクロック周波数を例示するテーブルである。
【図９】図１に示した受信機におけるスペクトル再生処理を説明する図である。
【図１０】図１に示した受信機における観測スペクトルの第１の展開例を示す図である。
【図１１】図１に示した受信機における観測スペクトルの第２の展開例を示す図である。
【図１２】図１に示した受信機における観測スペクトルの第３の展開例及び復元されたスペクトルの一例を示す図である。
【図１３】図１に示した受信機におけるスペクトルのモニタ結果を例示する図である。
【図１４】変形例としての受信機の構成例を示すブロック図である。
【発明を実施するための形態】
【００１２】
以下、図面を参照して実施の形態を説明する。ただし、以下に示す実施形態はあくまでも例示に過ぎず、実施形態で明示しない種々の変形例や技術の適用を排除する意図はない。すなわち、本実施形態を、その趣旨を逸脱しない範囲で種々変形して実施することができる。
【００１３】
また、各図は、図中に示す構成要素のみを備えるという趣旨ではなく、他の構成要素を含むことができる。以下、図中において、同一の符号を付した部分は特に断らない限り、同一若しくは同様の部分を示す。
【００１４】
〔Ａ〕実施形態
図１は、実施形態としての受信機１００の構成例を示すブロック図である。
【００１５】
図１に示す受信機１００は、ダイレクトＲＦアンダーサンプリング受信方式を採用したリアルタイムスペクトルモニタとして機能する。
【００１６】
従来のダイレクトＲＦアンダーサンプリング受信機では、アンダーサンプリング後の第１ナイキスト領域での折り返しによる受信周波数のオーバーラップを防ぐため、受信帯域幅の合計はクロック周波数の１／２（別言すれば、「ナイキスト限界」）未満でなければならなかった。
【００１７】
しかしながら、図１に示す受信機１００では、９００ＭＨｚ〜６ＧＨｚの周波数帯域をカバーし、数百ＭＨｚのクロック周波数のみを必要とする。
【００１８】
図示する受信機１００は、例示的に、アンテナ１，Low Noise Amplifier（ＬＮＡ；低雑音増幅器）２，複数（図示する例では３つ）のBand-Pass Filter（ＢＰＦ）３１〜３３，Variable Gain Amplifier（ＶＧＡ）４，Received Signal Strength Indicator（ＲＳＳＩ）５，クロック（ＣＬＫ）６，Sample-And-Hold（Ｓ／Ｈ）回路７，Analog-to-Digital Converter（ＡＤＣ）８，Digital Signal Processor（ＤＳＰ）９及びモニタ１０を備える。
【００１９】
アンテナ１は、高周波信号を受信する。アンテナ１が受信する高周波信号は、マイクロ波帯の周波数を有する。マイクロ波帯の高周波信号としては、例えば、Ultra High Frequency（ＵＨＦ；極超短波）帯、センチメートル波、ミリ波帯、又は、サブミリ波帯の周波数を有するものが挙げられる。アンテナ１は、受信された高周波信号を、ＬＮＡ２へ出力する。
【００２０】
ＬＮＡ２は、アンテナ１から高周波信号が入力される。ＬＮＡ２は、入力された高周波信号を増幅して、増幅された高周波信号を、ＢＰＦ３１〜３３へ出力する。ここで、ＬＮＡ２は、マルチバンド高周波信号に含まれる入力信号を増幅することができる、広帯域の増幅器である。
【００２１】
ＢＰＦ３１〜３３は、ＬＮＡ２から増幅された高周波信号が入力される。ＢＰＦ３１〜３３は、２以上（図示する例では３つ）相互に並列接続されている。ＢＰＦ３１〜３３は、入力された高周波信号のうちの、所定の周波数帯域である通過帯域に含まれない成分を、当該通過帯域に含まれる成分よりも大きく減衰させた信号を出力するものである。
【００２２】
ＢＰＦ３１〜３３は、アンチエリアスフィルタの一例であり、それぞれマルチバンド高周波信号における異なった通過帯域を有する。ＢＰＦ３１〜３３が有する通過帯域の最高周波数は、アンダーサンプリング次数ｍ（ｍはアンダーサンプリング次数で１以上の整数）の周波数帯域（２ｍ−１）ｆｓ〜（２ｍ＋１）ｆｓに含まれている。ＢＰＦ３１〜３３は、入力された高周波信号を各ＢＰＦ３１〜３３によってろ波して、各ＢＰＦ３１〜３３の通過帯域に含まれるろ波された信号成分が合わせられた信号を、通過信号としてＶＧＡ４へ出力する。
【００２３】
図示する例において、ＢＰＦ３１は９２０ＭＨｚ帯域の信号を通過させ、ＢＰＦ３２は２．４ＧＨｚ帯域の信号を通過させ、ＢＰＦ３３は５ＧＨｚ帯域の信号を通過させる。
【００２４】
ＶＧＡ４は、ＢＰＦ３１〜３３からの通過信号を増幅して出力する。
【００２５】
ＲＳＳＩ５は、ＶＧＡ４が出力した信号の信号強度を測定して、ＶＧＡ４における通過信号の増幅のためにＶＧＡ４にフィードバックする。
【００２６】
Ｓ／Ｈ回路７には、ＢＰＦ３１〜３３を通過してＶＧＡ４によって増幅された高周波信号（別言すれば、「通過信号」）が入力されるとともに、ＣＬＫ６からクロック信号が入力される。Ｓ／Ｈ回路７は、ＣＬＫ６から入力されたクロック信号に同期したサンプリング周波数にて、ＢＰＦ３１〜３３から入力された通過信号に対するサンプリングを行ない、標本化信号を得る。
【００２７】
本実施形態では、クロック周波数ｆＣＬＫは、入力信号のナイキスト周波数よりも低い。従って、Ｓ／Ｈ回路７は、高周波信号に対するアンダーサンプリングを行なうことになる。また、本実施形態では、受信した高周波信号について、予め周波数変換を施したもの対してアンダーサンプリングを行うものではなく、Ｓ／Ｈ回路７が、直接アンダーサンプリングを行なうものである。Ｓ／Ｈ回路７は、サンプリング後の標本化信号をＡＤＣ８へ出力する。上述の通り、Ｓ／Ｈ回路７は、増幅及びろ波を受けた高周波信号として、アンダーサンプリング次数の異なる複数の帯域を有するマルチバンド高周波信号が入力されて、このマルチバンド高周波信号に対するアンダーサンプリングを行なうものである。
【００２８】
ＣＬＫ６は、クロック周波数ｆＣＬＫを有するクロック信号を生成する発振器である。ＣＬＫ６は、可変クロック信号源として機能する。ＣＬＫ６は、高周波信号に含まれる入力信号の帯域幅の２倍よりも低いクロック周波数ｆＣＬＫを有するクロック信号を生成し、生成したクロック信号をＳ／Ｈ回路７及びＡＤＣ８へ出力する。
【００２９】
ＡＤＣ８は、Ｓ／Ｈ回路７から標本化信号が入力される。ＡＤＣ８は、入力された標本化信号を、アナログ信号からデジタル信号に変換する。ＡＤＣ８は、変換後の標本化信号をＤＳＰ９に出力する。
【００３０】
ＤＳＰ９は、周波数スペクトラム再生装置の一例であり、Fast Fourier Transform（ＦＦＴ）回路９１及び受信信号周波数推定部９２を備える。
【００３１】
ＦＦＴ回路９１は、ＡＤＣ８から入力された標本化信号を高速フーリエ変換する。
【００３２】
受信信号周波数推定部９２は、高速フーリエ変換された標本化信号に基づき、受信機１００によって受信された信号の周波数を推定する。
【００３３】
モニタ１０は、表示装置の一例であり、液晶ディスプレイやOrganic Light-Emitting Diode（ＯＬＲＤ）ディスプレイ，Cathode Ray Tube（ＣＲＴ），電子ペーパーディスプレイ等であり、ＤＳＰ９による受信信号周波数の推定結果を出力する。モニタ１０は、受信信号の周波数スペクトラムを表示する。
【００３４】
図２は、図１に示した受信機１００におけるアンダーサンプリング処理を説明する図である。
【００３５】
図２に示す例では、受信周波数帯域が、ＺＯＮＥ１〜４で示されたｆＳ／２の帯域幅にそれぞれ分けられる。ナイキストサンプリングが実施されるＺＯＮＥ１に対して、ＺＯＮＥ２〜４が折り返されることにより、アンダーサンプリングが実現される。
【００３６】
これにより、低いサンプリングクロックで広帯域信号を受信可能であり、広帯域とリアルタイム性とが両立できる。ただし、どのナイキストゾーンから折り返されてきたのかが不明であり、受信周波数を特定できない。
【００３７】
図３は、図１に示した受信機１００におけるサンプリングクロックの切り替え前のアンダーサンプリング受信動作を説明する図である。
【００３８】
なお、本実施形態では、受信機１００と異なる送信装置により、それぞれ周波数ｆＲＦ１＝９２０ＭＨｚ、ｆＲＦ２＝２．４ＧＨｚ、及びｆＲＦ３＝５．８ＧＨｚの入力信号が含まれる高周波信号が送信される。そして、受信機１００は、この高周波信号を受信する。
【００３９】
まず、受信機１００のアンテナ１は、高周波信号を受信し、受信された高周波信号をＬＮＡ２へ出力する。
【００４０】
ＬＮＡ２は、アンテナ１から入力された高周波信号を増幅し、増幅された高周波信号をＢＰＦ３１〜３３へ出力する。
【００４１】
ＢＰＦ３１〜３３は、ＬＮＡ２から入力された高周波信号のうちの、通過帯域ｆ１，ｆ２，ｆ３に含まれる高周波信号を、ＶＧＡ４を介してＳ／Ｈ回路７へ出力する。なお、本実施形態では、ＢＰＦ３１の通過帯域ｆ１は、９１５ＭＨｚ〜９３０ＭＨｚである。また、ＢＰＦ３２の通過帯域ｆ２は、２．４ＧＨｚ〜２．５ＧＨｚである。また、ＢＰＦ３３の通過帯域ｆ３は、５．０ＧＨｚ〜５．８７５ＧＨｚである。これにより、ＢＰＦ３１〜３３の通過帯域外の信号が減衰されて、周波数ｆＲＦ１，ｆＲＦ２及びｆＲＦ３の信号を含むマルチバンド高周波信号がＳ／Ｈ回路７へ出力される。
【００４２】
ＣＬＫ６は、クロック周波数ｆＣＬＫを有するクロック信号を生成し、生成したクロック信号をＳ／Ｈ回路７へ出力する。図３に示す例では、クロック周波数ｆＣＬＫは１ＧＨｚである。
【００４３】
Ｓ／Ｈ回路７は、ＣＬＫ６から入力された１ＧＨｚのクロック信号に同期して、１ＧＨｚのサンプリング周波数ｆｓにて、ＢＰＦ３１〜３３から入力された高周波信号に対するアンダーサンプリングを行なう。Ｓ／Ｈ回路７は、サンプリング後の標本化信号をＡＤＣ８へ出力する。このとき、サンプリング周波数ｆｓの１／２（＝５００ＭＨｚ）よりも高い周波数を有するｎ次の周波数帯域の信号は、エリアシングによって、５００ＭＨｚ以下のベースバンドの領域に折り返される。
【００４４】
図３では、横軸に周波数を、縦軸に信号の振幅（別言すれば、「信号強度」）を示している。本実施形態では、図３に表すように、周波数ｆＲＦ１の信号は、エリアシングによって、周波数ｆＩＦ１＝８０ＭＨｚの信号に変換される。また、周波数ｆＲＦ２の信号は、周波数ｆＩＦ２＝１００ＭＨｚの信号に変換される。また、周波数ｆＲＦ３の信号は、周波数ｆＩＦ３＝２００ＭＨｚの信号に変換される。すなわち、Ｓ／Ｈ回路７によるアンダーサンプリングによって、周波数ｆＩＦ１，ｆＩＦ２及びｆＩＦ３の信号を含む標本化信号が得られる。
【００４５】
ＡＤＣ８は、Ｓ／Ｈ回路７から入力された標本化信号を、アナログ信号からデジタル信号に変換し、変換後の標本化信号をＤＳＰ９へ出力する。
【００４６】
ＤＳＰ９は、デジタル信号に変換された標本化信号から受信信号周波数を推定してモニタ１０に出力する。
【００４７】
図４は、図１に示した受信機１００におけるサンプリングクロックの切り替え後のアンダーサンプリング受信動作を説明する図である。
【００４８】
図４に示す例では、クロック周波数ｆＣＬＫは７００ＭＨｚである。図３及び図４に示すように、異なる周波数でサンプリングを行なうと、スペクトルの観測値も変化する。スペクトル観測値の変化から周波数を一意に識別することが可能である。
【００４９】
ある周波数ｆＲＦのＲＦ信号をサンプリング周波数ｆｓ１でアンダーサンプリングした際、第１ナイキストゾーンで観測される周波数ｆＭＥＳ１とｆＲＦとの関係は次の数式（１）で表わされる。
【００５０】
【数１】
2020106478.tif 000003
【００５１】
ｎ１はアンダーサンプリング次数（ｎ１＝１，２，・・・）である。
【００５２】
ｆＲＦを特定するにはサンプリングクロックをｆｓ２に切り替えて第１ナイキストゾーンで観測される周波数ｆＭＥＳ２が候補の各ｆＲＦに対して異なる値となればよい。
【００５３】
【数２】
2020106478.tif 000004
【００５４】
ｎ２は数式（２）を満たす整数（ｎ２＝０，１，・・・）である。
【００５５】
ｆＭＥＳ２が重複するｆｓ２の条件は次の数式（３）〜（５）で与えられ、これらの数式を満足しないようにｆｓ２を設定すればｆＲＦを特定できる。
【００５６】
【数３】
2020106478.tif 000005
【００５７】
【数４】
2020106478.tif 000006
【００５８】
【数５】
2020106478.tif 000007
【００５９】
ｎ１，１，ｎ１，２はｎ１を区別した表記であり、ｎ２，１，ｎ２，２はｎ２を区別した値である。
【００６０】
数式（３）を満足しないｆｓ２の条件は、以下の数式（６）及び（７）を満足するときである。ｆｓ１及びｆｓ２は、それぞれ１回目及び２回目のサンプリング周波数である。Ｍ及びＬは整数であり、ＭとＬとは互いに素である。Ｎは１回目のサンプリングにおける最大アンダーサンプリング次数であり、入力信号の最高周波数とサンプリング周波数とによって決定される。
【００６１】
数式（４）及び（５）にはｆＭＥＳ１が含まれており、その値を得るにはＦＦＴによる処理が必要である。ＦＦＴの処理時間を考慮すると、ミリ秒オーダーでのクロックの切り替えができずリアルタイム性が確保できない。そこで、ｆＭＥＳ１は未知として複数回（別言すれば、「２回以上」）のクロック切り替えが実行される。サンプリングクロックをｆｓ２からｆｓ３に切り替えたとき、数式（４）及び（５）を満たすｆＭＥＳ１がｆｓ２とｆｓ３とで異なれば、どちらか一方のサンプリングクロックで重複しても、もう一方のサンプリングクロックで重複を回避できるため、ｆＲＦが特定できる。
【００６２】
【数６】
2020106478.tif 000008
【００６３】
【数７】
2020106478.tif 000009
【００６４】
図５は、図１に示した受信機１００における測定周波数帯域とサンプリング周波数との関係を例示するテーブルである。図６は、図１に示した受信機１００におけるアンダーサンプリング受信時の観測周波数の計算結果を例示するグラフである。
【００６５】
図６は、図５に示す条件で、ｆＭＥＳ１に対する第２回及び第３回サンプリングの際に観測されるｆＭＥＳ２及びｆＭＥＳ３の計算結果である。ただし、第２回及び第３回のサンプリング周波数は数式（６）及び（７）を満たす値としている。
【００６６】
符号Ａ１はｆＲＦ＝９１５〜９３０ＭＨｚの周波数帯域を示し、符号Ａ２はｆＲＦ＝２４００〜２５００ＭＨｚの周波数帯域を示し、符号Ａ３はｆＲＦ＝５０００〜５２５０ＭＨｚの周波数帯域を示す。また、符号Ａ４はｆＲＦ＝５２５０〜５５００ＭＨｚの周波数帯域を示し、符号Ａ５はｆＲＦ＝５５００〜５７５０ＭＨｚの周波数帯域を示し、符号Ａ６はｆＲＦ＝５７５０〜５８７５ＭＨｚの周波数帯域を示す。
【００６７】
ｆＲＦは無線周波数を表わす。ｆＭＥＳ１，ｆＭＥＳ２及びｆＭＥＳ３は、ｆＲＦの信号をｆｓ１，ｆｓ２及びｆｓ３でそれぞれサンプリングした場合に第１ナイキストゾーンで観測される周波数を表わす。
【００６８】
図６において、１回目のサンプリング時に符号Ｂ１で示すｆＭＥＳ１で観測された信号は、２回目にサンプリングしても、符号Ｂ２で示すｆＭＥＳ２で観測されてしまい、元のＲＦ周波数を特定することができないことがある。このような特異点となる条件下においては３回目のサンプリングが必要となり、これにより符号Ｂ３及びＢ４の２つのｆＭＥＳ３に分離することができ、よって元のＲＦ周波数を特定することができる。今回の図５に示す条件においては図６から明らかなように、３回の異なるサンプリング周波数での観測により受信帯域内の全てのＲＦ信号の周波数を特定することができる。
【００６９】
図７は、図１に示した受信機１００におけるクロック切り替えを用いた周波数同定法を説明する図である。
【００７０】
図７に示す通り、３つの異なる周波数のクロックをサンプリングするアンダーサンプリング受信により、受信周波数を識別することが可能である。この測定サイクル時間Ｔｄ
はミリ秒単位よりも短い。なお、Ｔｄは、サンプリング周波数ｆｓ１の時間Ｔｄ１と、サンプリング周波数ｆｓ２の時間Ｔｄ２と、サンプリング周波数ｆｓ３の時間Ｔｄ３との和である。Ｔｄが１ｍｓよりも短ければ、リアルタイム性は満たされる。例えば、サンプリング時間が０．１ｍｓとなれば、分解能帯域幅（ＲＢＷ）は１０ｋＨｚとなる。クロック切り替え時間を無視すれば、サンプリング時間（別言すれば、「測定サイクル時間」）Ｔｄの合計は０．３ｍｓとなる。
【００７１】
本実施形態では、パルスパターン発生器（ＰＰＧ：例えば、M8020A Keysight）が高速スイッチングソースとして使用される。これにより、１６ＧＨｚを分周して得たクロック周波数に切り替えることが可能となる。各クロック周波数は、５００ＭＨｚ以下で数式（６）及び（７）を満たす様に設計される。
【００７２】
図８は、図１に示した受信機１００におけるクロック周波数を例示するテーブルである。
【００７３】
ＰＰＧクロック周波数の切り替え時間は、サンプリング時間に比べれば無視することが可能である。測定サイクル時間Ｔｄは２２．９３２マイクロ秒となる。よって、リアルタイム性能は十分に満たされる。
【００７４】
よって、本実施形態の受信機１００におけるＤＳＰ９によれば、ミリ秒単位以内のリアルタイム性能及び広帯域の両方を実現させて、受信信号周波数をモニタできる。
【００７５】
図９は、図１に示した受信機１００におけるスペクトル再生処理を説明する図である。
【００７６】
図９に示すように、第１回〜第３回の各サンプリングの観測スペクトルが高次ナイキストゾーンに展開され、展開スペクトルが生成される（符号Ｃ１参照）。そして、元のスペクトルはどの回のサンプリングにおいても展開スペクトルが一致するため、各回の展開スペクトルを最小値演算することでアンダーサンプリング前の元のスペクトルが復元（別言すれば、「再生」）できる（符号Ｃ２参照）。
【００７７】
すなわち、図１に示したＦＦＴ回路９１は、１回目からｎ（ｎは２以上の整数）回目までのサンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得する手段として機能する。
【００７８】
図１に示した受信信号周波数推定部９２は、１回目からｎ回目までのデータから取得された周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までの周波数スペクトラムを生成する手段として機能する。
【００７９】
また、受信信号周波数推定部９２は、得られた１回目からｎ回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成する手段としても機能する。
【００８０】
このように、受信機１００のＤＳＰ９は、サンプリング周波数をｎ−１回変えてサンプリングしたｎ個の時間軸データから、受信周波数スペクトラムを再生する。
【００８１】
これにより、異なる複数の周波数帯域のリアルタイム監視を行なうことができる。
【００８２】
図１０〜図１２は、図１に示した受信機１００における観測スペクトルの展開例及び復元されたスペクトルの一例を示す図である。
【００８３】
図１０に示すように、第１回サンプリングにおいて、クロック周波数ｆＣＬＫ１＝５００ＭＨｚとして、観測スペクトルが高次に展開されて展開スペクトルが生成される（符号Ｄ１参照）。
【００８４】
図１１に示すように、第２回サンプリングにおいて、クロック周波数ｆＣＬＫ２＝４１０．２５６ＭＨｚとして、観測スペクトルが高次に展開されて展開スペクトルが生成される（符号Ｄ２参照）。
【００８５】
図１２に示すように、第３回サンプリングにおいて、クロック周波数ｆＣＬＫ２＝３２６．５３０ＭＨｚとして、観測スペクトルが高次に展開されて展開スペクトルが生成される（符号Ｄ３参照）。
【００８６】
そして、図１２に示すように、図１０の符号Ｄ１，図１１の符号Ｄ２及び図１２の符号Ｄ３で示した展開スペクトルの最小値演算結果として、復元されたスペクトルが生成される（符号Ｄ４参照）。
【００８７】
ここで、第４回サンプリングの展開スペクトルを取得した場合には、第２回〜第４回サンプリングの展開スペクトルにより、復元したスペクトルを生成してよい。また、第５回サンプリングの展開スペクトルを取得した場合には、第３回〜第５回サンプリングの展開スペクトルにより、復元したスペクトルを生成してよい。更に、第６回サンプリングの展開スペクトルを取得した場合には、第４回〜第６回サンプリングの展開スペクトルにより、復元したスペクトルを生成してよい。このように、受信機１００のＤＳＰ９は、サンプリング周波数をｎ−１回変えてサンプリングする動作を繰り返し行なう。
【００８８】
すなわち、図１に示したＦＦＴ回路９１は、ｋ（ｋはｎ以上の整数）回目の前記サンプリング周波数でサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までの周波数スペクトラムを取得する手段として機能する。
【００８９】
図１に示した受信信号周波数推定部９２は、ｋ回目のデータから得られた周波数スペクトラムをナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１倍の周波数までの周波数スペクトラムを生成する手段として機能する。
【００９０】
また、受信信号周波数推定部９２は、得られたｋ回目からｋ−ｎ＋１回目までの周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成する手段としても機能する。
【００９１】
これにより、サンプリングを実施する毎に、復元スペクトルのモニタリングを実施できる。
【００９２】
図１３は、図１に示した受信機１００におけるスペクトルのモニタ結果を例示する図である。
【００９３】
図１３には、ＲＦＩＤ（９２０ＭＨｚ帯），無線ＬＡＮ及びBluetooth（登録商標）（２．４ＧＨｚ帯）デバイスをアンテナ１の近くに設置した場合のワイヤレスＩｏＴ周波数帯のスペクトログラムが示される。
【００９４】
符号Ｅ１に示すスペクトログラムにおいては、中心周波数９１９ＭＨｚのＲＦＩＤ信号が確認できる。
【００９５】
符号Ｅ２に示すスペクトログラムにおいては、無線ＬＡＮ及びBluetooth信号が確認できる。符号Ｅ３には、符号Ｅ２に示したスペクトログラムの時間軸をズームした２．４ＧＨｚ帯のスペクトログラムが示される。これにより、２２．９３２マイクロ秒の時間分解能でのスペクトルモニタリングが実現可能である。よって、ミリ秒単位のスペクトルモニタリングが実現される。
【００９６】
符号Ｅ４に示すスペクトログラムは、５ＧＨｚ帯域の全ての信号がノイズフロアレベル未満であることを示す。
【００９７】
〔Ｂ〕変形例
図１４は、変形例としての受信機２００の構成例を示すブロック図である。
【００９８】
図１４に示す受信機２００は、図１に示した受信機１００と同様に、アンテナ１，ＬＮＡ２，複数（図示する例では３つ）のＢＰＦ３１〜３３，ＶＧＡ４，ＲＳＳＩ５，ＤＳＰ９及びモニタ１０を備える。また、受信機２００は、ｎ組（ｎは２以上の整数；図示する例では３組）のＣＬＫ６，Ｓ／Ｈ回路７及びＡＤＣ８を備える。
【００９９】
受信機２００は、複数組のＣＬＫ６，Ｓ／Ｈ回路７及びＡＤＣ８を備えることにより、図９等に示した複数回のサンプリングを並行して実施する。図１４に示す受信機２００は、３組のＣＬＫ６，Ｓ／Ｈ回路７及びＡＤＣ８を備えているため、３回のサンプリングを並行して実施する。
【０１００】
すなわち、複数組のＣＬＫ６，Ｓ／Ｈ回路７及びＡＤＣ８は、ＬＮＡ２とＢＰＦ３１〜３３との後段に並列に接続され、ｎ（ｎは２以上の整数）個の異なるサンプリング周波数でアンダーサンプリングする受信回路として機能する。
【０１０１】
ＦＦＴ回路９１は、１個目からｎ個目までの受信回路で取得してサンプリングした時間軸データを高速フーリエ変換して、０からナイキスト周波数までのｎ個の周波数スペクトラムを取得する取得部として機能する。
【０１０２】
受信信号周波数推定部９２は、１個目からｎ個目までの受信回路で取得したｎ個の周波数スペクトラムをそれぞれナイキスト周波数毎に折り返してナイキスト周波数の２ｍ＋１（ｍはアンダーサンプリング次数で１以上の整数）倍の周波数までのｎ個の周波数スペクトラムを生成する第１生成部として機能する。
【０１０３】
また、受信信号周波数推定部９２は、得られた１個目からｎ個目までの受信回路の各周波数スペクトラムの周波数毎のデータそれぞれについて最小値演算を行ない、１つの周波数スペクトラムを生成する第２生成部としても機能する。
【０１０４】
これにより、複数のサンプリング周波数でのアンダーサンプリングを並行して実施できるため、モニタリング処理の時間を短縮できる。
【０１０５】
〔Ｃ〕その他
開示の技術は上述した各実施形態に限定されるものではなく、各実施形態の趣旨を逸脱しない範囲で種々変形して実施することができる。各実施形態の各構成及び各処理は、必要に応じて取捨選択することができ、あるいは適宜組み合わせてもよい。
【０１０６】
図１及び図１４に示した例においては、受信機１００及び２００が３つのＢＰＦ３１〜３３を備えることとしたが、これに限定されるものではない。受信機１００及び２００が備えるＢＰＦは２つ又は４つ以上であってもよい。これにより、種々の周波数帯域の信号をモニタリングできる。
【０１０７】
図９等に示した例においては、１〜３回のサンプリングによって元のスペクトルを復元することとしたが、これに限定されるものではない。４回以上のサンプリングによって、１つのスペクトルを復元してもよい。これにより、スペクトルの復元精度やノイズ耐性を向上できる。
【符号の説明】
【０１０８】
１００，２００  ：受信機
１              ：アンテナ
２              ：ＬＮＡ
３１，３２，３３：ＢＰＦ
４              ：ＶＧＡ
５              ：ＲＳＳＩ
６              ：ＣＬＫ
７              ：Ｓ／Ｈ回路
８              ：ＡＤＣ
９              ：ＤＳＰ
９１            ：ＦＦＴ回路
９２            ：受信信号周波数推定部
１０            ：モニタ
【図１】
2020106478.tif 000010
【図２】
2020106478.tif 000011
【図３】
2020106478.tif 000012
【図４】
2020106478.tif 000013
【図５】
2020106478.tif 000014
【図６】
2020106478.tif 000015
【図７】
2020106478.tif 000016
【図８】
2020106478.tif 000017
【図９】
2020106478.tif 000018
【図１０】
2020106478.tif 000019
【図１１】
2020106478.tif 000020
【図１２】
2020106478.tif 000021
【図１３】
2020106478.tif 000022
【図１４】
2020106478.tif 000023