(57)【特許請求の範囲】
【請求項１】
識別対象の無線機から送信された無線信号を受信する信号受信部と、
前記信号受信部により受信された無線信号を短時間フーリエ変換し、前記無線信号の短時間フーリエ変換結果を示すフーリエ変換信号を出力するフーリエ変換部と、
前記フーリエ変換部から出力されたフーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形を算出し、前記エネルギーの時間波形から前記無線信号の立ち上がり時刻を検出する立ち上がり検出部と、
前記フーリエ変換部から出力されたフーリエ変換信号の中から、前記立ち上がり検出部により検出された立ち上がり時刻を含む時間帯のフーリエ変換信号を抽出し、抽出したフーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出したフーリエ変換信号の変動を特徴量として算出する特徴量算出部と、
前記特徴量算出部により算出された特徴量に基づいて、前記無線機を識別する無線機識別部と
を備えた無線機識別装置。
【請求項２】
前記信号受信部は、既知の無線機から送信された無線信号を受信し、
前記フーリエ変換部は、前記信号受信部により受信された既知の無線機の無線信号を短時間フーリエ変換し、既知の無線機の無線信号の短時間フーリエ変換結果を示すフーリエ変換信号を出力し、
前記立ち上がり検出部は、前記フーリエ変換部から出力された既知の無線機に係るフーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形を算出し、当該エネルギーの時間波形から既知の無線機の無線信号の立ち上がり時刻を検出し、
前記特徴量算出部は、前記フーリエ変換部から出力された既知の無線機に係るフーリエ変換信号の中から、前記立ち上がり検出部により検出された立ち上がり時刻を含む時間帯のフーリエ変換信号を抽出し、抽出した既知の無線機に係るフーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出した既知の無線機に係るフーリエ変換信号の変動を特徴量として算出し、算出した既知の無線機に係る特徴量をデータベース部に格納し、
前記無線機識別部は、前記特徴量算出部により算出された識別対象の無線機に係る特徴量と、前記データベース部により格納されている既知の無線機に係る特徴量とを比較し、特徴量の比較結果に基づいて、識別対象の無線機を識別することを特徴とする請求項１記載の無線機識別装置。
【請求項３】
特徴量算出用のパラメータとして、前記時間帯の長さを示すパラメータと、前記設定周波数成分を示すパラメータとを設定する第１のパラメータ設定部を備え、
前記特徴量算出部は、前記第１のパラメータ設定部により設定された特徴量算出用のパラメータを用いて、識別対象の無線機に係る特徴量及び既知の無線機に係る特徴量のそれぞれを算出することを特徴とする請求項２記載の無線機識別装置。
【請求項４】
前記第１のパラメータ設定部は、互いに異なる既知の無線機から送信された無線信号のそれぞれが前記信号受信部に受信された際に、前記特徴量算出部によりそれぞれ算出される特徴量が重ならないように、あるいは、前記特徴量算出部によりそれぞれ算出される特徴量の重なりが最も小さくなるように、前記特徴量算出用のパラメータを更新することを特徴とする請求項３記載の無線機識別装置。
【請求項５】
前記立ち上がり検出部は、前記スペクトログラムに含まれている複数の周波数における電力のうち、同一時刻の電力同士をそれぞれ合計することで、前記エネルギーの時間波形として、各々の時刻における電力の合計値を算出し、無線信号の立ち上がり時刻として、前記電力の合計値が閾値に到達した時刻を検出することを特徴とする請求項１記載の無線機識別装置。
【請求項６】
前記特徴量算出部は、
第１から第３の特徴量として、前記複素時間信号における振幅の分散値、歪度及び尖度をそれぞれ算出し、第４から第６の特徴量として、前記複素時間信号における位相の分散値、歪度及び尖度をそれぞれ算出し、第７から第９の特徴量として、前記複素時間信号における周波数の分散値、歪度及び尖度をそれぞれ算出することを特徴とする請求項１記載の無線機識別装置。
【請求項７】
前記第１から第９の特徴量の重みを示す重み付けパラメータを設定する第２のパラメータ設定部を備え、
前記特徴量算出部は、前記第２のパラメータ設定部により設定された重み付けパラメータに従って前記第１から第９の特徴量の重み付けを実施することを特徴とする請求項６記載の無線機識別装置。
【請求項８】
信号受信部が、識別対象の無線機から送信された無線信号を受信し、
フーリエ変換部が、前記信号受信部により受信された無線信号を短時間フーリエ変換し、前記無線信号の短時間フーリエ変換結果を示すフーリエ変換信号を出力し、
立ち上がり検出部が、前記フーリエ変換部から出力されたフーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形を算出し、前記エネルギーの時間波形から前記無線信号の立ち上がり時刻を検出し、
特徴量算出部が、前記フーリエ変換部から出力されたフーリエ変換信号の中から、前記立ち上がり検出部により検出された立ち上がり時刻を含む時間帯のフーリエ変換信号を抽出し、抽出したフーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出したフーリエ変換信号の変動を特徴量として算出し、
無線機識別部が、前記特徴量算出部により算出された特徴量に基づいて、前記無線機を識別する
無線機識別方法。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、無線機を識別する無線機識別装置及び無線機識別方法に関するものである。
【背景技術】
【０００２】
無線機から送信される無線信号の立ち上がりなどの過渡応答は、個体差があることが知られており、近年では、無線機から送信される無線信号を解析することで、無線機を識別する無線機識別装置が提案されている。
例えば、無線機識別装置は、無線機から送信される無線信号の立ち上がりの特徴量を算出し、算出した特徴量に基づいて、無線機を識別する。
無線信号の立ち上がりの特徴量としては、立ち上がり時の無線信号における振幅等の分散値、歪度及び尖度などがある。【０００３】
ただし、信号対雑音比 （Ｓｉｇｎａｌ−ｔｏ−Ｎｏｉｓｅ  Ｒａｔｉｏ：ＳＮＲ）が低い環境下では、立ち上がりの特徴量の算出精度が低下して、無線機の識別精度が低下する。このため、無線機から送信される無線信号を短時間フーリエ変換することで、ＳＮＲを改善してから、立ち上がりの特徴量を算出する無線機識別装置が以下の特許文献１に開示されている。
この無線機識別装置では、立ち上がりの特徴量として、立ち上がりの時定数を算出している。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特開２００２−２６８２６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
従来の無線機識別装置は、立ち上がりの時定数に基づく無線機の識別方法を採用している。立ち上がりの時定数に基づく無線機の識別方法では、無線機の種類については識別することが可能である。しかし、同一種類の無線機は、異なる無線機であっても、立ち上がりの時定数の差異が小さいため、無線機の個体を識別することができないという課題があった。【０００６】
この発明は上記のような課題を解決するためになされたもので、無線機の個体識別を行うことができる無線機識別装置及び無線機識別方法を得ることを目的とする。
【課題を解決するための手段】
【０００７】
この発明に係る無線機識別装置は、識別対象の無線機から送信された無線信号を受信する信号受信部と、信号受信部により受信された無線信号を短時間フーリエ変換し、無線信号の短時間フーリエ変換結果を示すフーリエ変換信号を出力するフーリエ変換部と、フーリエ変換部から出力されたフーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形を算出し、エネルギーの時間波形から無線信号の立ち上がり時刻を検出する立ち上がり検出部と、フーリエ変換部から出力されたフーリエ変換信号の中から、立ち上がり検出部により検出された立ち上がり時刻を含む時間帯のフーリエ変換信号を抽出し、抽出したフーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出したフーリエ変換信号の変動を特徴量として算出する特徴量算出部とを設け、無線機識別部が、特徴量算出部により算出された特徴量に基づいて、無線機を識別するようにしたものである。
【発明の効果】
【０００８】
この発明によれば、フーリエ変換部から出力されたフーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形を算出し、エネルギーの時間波形から無線信号の立ち上がり時刻を検出する立ち上がり検出部と、フーリエ変換部から出力されたフーリエ変換信号の中から、立ち上がり検出部により検出された立ち上がり時刻を含む時間帯のフーリエ変換信号を抽出し、抽出したフーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出したフーリエ変換信号の変動を特徴量として算出する特徴量算出部とを設け、無線機識別部が、特徴量算出部により算出された特徴量に基づいて、無線機を識別するように構成したので、無線機の個体を識別することができる効果がある。
【図面の簡単な説明】
【０００９】
【図１】この発明の実施の形態１による無線機識別装置を示す構成図である。
【図２】この発明の実施の形態１による無線機識別装置を示すハードウェア構成図である。
【図３】無線機識別装置における信号受信部１を除く構成要素がソフトウェア又はファームウェアなどで実現される場合のコンピュータのハードウェア構成図である。
【図４】この発明の実施の形態１による無線機識別方法における学習処理を示すフローチャートである。
【図５】フーリエ変換部２による学習用無線信号ｓｇ（ｎ）及び識別用無線信号ｓｄ（ｎ）のＳＴＦＴを示す説明図である。
【図６】立ち上がり検出部３による立ち上がり時刻Ｔ０の検出処理を示す説明図である。
【図７】特徴量算出部７による特徴量の算出処理を示す説明図である。
【図８】図８Ａは、特徴量算出用のパラメータが適正でないために、種類（Ａ）の無線機に係る特徴量ベクトル〇と、種類（Ｂ）の無線機に係る特徴量ベクトル◇と、種類（Ｃ）の無線機に係る特徴量ベクトル△とが一部重なっている例を示す説明図、図８Ｂは、特徴量算出用のパラメータが適正であるために、種類（Ａ）の無線機に係る特徴量ベクトル〇と、種類（Ｂ）の無線機に係る特徴量ベクトル◇と、種類（Ｃ）の無線機に係る特徴量ベクトル△とが重なっていない例を示す説明図である。
【図９】第１のデータベース部９に格納されているシリアルナンバーが付加されている学習用無線信号を示す説明図である。
【図１０】第２のデータベース部１０に格納される複数の既知の無線機に係る特徴量の一例を示す説明図である。
【図１１】この発明の実施の形態１による無線機識別方法における識別処理を示すフローチャートである。
【図１２】無線機識別部１１による無線機の識別処理を示す説明図である。
【発明を実施するための形態】
【００１０】
以下、この発明をより詳細に説明するために、この発明を実施するための形態について、添付の図面に従って説明する。
【００１１】
実施の形態１．
図１は、この発明の実施の形態１による無線機識別装置を示す構成図である。
図２は、この発明の実施の形態１による無線機識別装置を示すハードウェア構成図である。
図１及び図２において、信号受信部１は、例えば図２に示す信号受信回路２１で実現される。
信号受信部１は、既知の無線機から送信された無線信号を学習用無線信号として受信し、受信した学習用無線信号をアナログ信号からデジタル信号に変換することで、デジタルの学習用無線信号をフーリエ変換部２及び第１のデータベース部９に出力する。
また、信号受信部１は、識別対象の無線機から送信された無線信号を識別用無線信号として受信し、受信した識別用無線信号をアナログ信号からデジタル信号に変換することで、デジタルの識別用無線信号をフーリエ変換部２に出力する。
【００１２】
フーリエ変換部２は、例えば図２に示すフーリエ変換回路２２で実現される。
フーリエ変換部２は、信号受信部１から出力されたデジタルの学習用無線信号を短時間フーリエ変換（ＳＴＦＴ：Ｓｈｏｒｔ  Ｔｉｍｅ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）し、学習用無線信号の短時間フーリエ変換結果であるＳＴＦＴ結果を示すフーリエ変換信号（以下、学習用フーリエ変換信号と称する）を立ち上がり検出部３及び特徴量算出部７に出力する処理を実施する。
また、フーリエ変換部２は、信号受信部１から出力されたデジタルの識別用無線信号をＳＴＦＴし、識別用無線信号のＳＴＦＴ結果を示すフーリエ変換信号（以下、識別用フーリエ変換信号と称する）を立ち上がり検出部３及び特徴量算出部７に出力する処理を実施する。
ＳＴＦＴは、一定時間毎に異なる時刻の無線信号をそれぞれ切り出し、切り出した各々の無線信号をそれぞれ高速フーリエ変換（ＦＦＴ：Ｆａｓｔ  Ｆｏｕｒｉｅｒ  Ｔｒａｎｓｆｏｒｍ）する手法である。
【００１３】
立ち上がり検出部３は、例えば図２に示す立ち上がり検出回路２３で実現される。
立ち上がり検出部３は、フーリエ変換部２から出力された学習用フーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形（以下、学習用のエネルギー時間波形と称する）を算出し、学習用のエネルギー時間波形から学習用無線信号における立ち上がり時刻を検出する処理を実施する。
また、立ち上がり検出部３は、フーリエ変換部２から出力された識別用フーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形（以下、識別用のエネルギー時間波形と称する）を算出し、識別用のエネルギー時間波形から識別用無線信号における立ち上がり時刻を検出する処理を実施する。
【００１４】
パラメータ設定部４は、第１のパラメータ設定部５及び第２のパラメータ設定部６を備えており、例えば図２に示すパラメータ設定回路２４で実現される。
第１のパラメータ設定部５は、特徴量算出用のパラメータとして、立ち上がり検出部３により検出された立ち上がり時刻を含む時間帯の長さである時間範囲Δｔを示すパラメータと、設定周波数成分である設定周波数ビンを示すパラメータとを設定する処理を実施する。
第２のパラメータ設定部６は、特徴量算出部７により算出される第１から第９の特徴量の重みを示す重み付けパラメータを設定する処理を実施する。
また、第２のパラメータ設定部６は、フーリエ変換部２によりＳＴＦＴが行われる際に用いられる窓関数のサンプル数及び窓関数をシフトさせるサンプル数を設定する処理を実施する。
【００１５】
特徴量算出部７は、例えば図２に示す特徴量算出回路２５で実現される。
特徴量算出部７は、パラメータ設定部４により設定された特徴量算出用のパラメータから、立ち上がり時刻を含む時間帯及び設定周波数ビンを認識する。
特徴量算出部７は、フーリエ変換部２から出力された学習用フーリエ変換信号の中から、立ち上がり検出部３により検出された学習用無線信号における立ち上がり時刻を含む時間帯の学習用フーリエ変換信号を抽出する処理を実施する。
そして、特徴量算出部７は、抽出した学習用フーリエ変換信号に含まれている設定周波数ビンの複素時間信号から、抽出した学習用フーリエ変換信号の変動を特徴量として算出する処理を実施する。
【００１６】
また、特徴量算出部７は、フーリエ変換部２から出力された識別用フーリエ変換信号の中から、立ち上がり検出部３により算出された識別用無線信号における立ち上がり時刻を含む時間帯の識別用フーリエ変換信号を抽出する処理を実施する。
そして、特徴量算出部７は、抽出した識別用フーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出した識別用フーリエ変換信号の変動を特徴量として算出する処理を実施する。
【００１７】
データベース部８は、第１のデータベース部９及び第２のデータベース部１０を備えており、例えば図２に示すデータベース回路２６で実現される。
第１のデータベース部９は、信号受信部１から出力されたデジタルの学習用無線信号を格納する。
第２のデータベース部１０は、既知の無線機に係る特徴量として、特徴量算出部７により算出された学習用フーリエ変換信号の変動を格納する。
【００１８】
無線機識別部１１は、例えば図２に示す無線機識別回路２７で実現される。
無線機識別部１１は、特徴量算出部７により算出された識別対象の無線機に係る特徴量と、第２のデータベース部１０により格納されている既知の無線機に係る特徴量とを比較する処理を実施する。
無線機識別部１１は、識別対象の無線機に係る特徴量と既知の無線機に係る特徴量との比較結果に基づいて、識別対象の無線機を識別する処理を実施する。
【００１９】
図１では、無線機識別装置の構成要素である信号受信部１、フーリエ変換部２、立ち上がり検出部３、パラメータ設定部４、特徴量算出部７、データベース部８及び無線機識別部１１のそれぞれが、図２に示すような専用のハードウェアで実現されるものを想定している。即ち、信号受信回路２１、フーリエ変換回路２２、立ち上がり検出回路２３、パラメータ設定回路２４、特徴量算出回路２５、データベース回路２６及び無線機識別回路２７で実現されるものを想定している。
【００２０】
ここで、データベース回路２６は、例えば、ＲＡＭ（Ｒａｎｄｏｍ  Ａｃｃｅｓｓ  Ｍｅｍｏｒｙ）、ＲＯＭ（Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、フラッシュメモリ、ＥＰＲＯＭ（Ｅｒａｓａｂｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒｙ）、ＥＥＰＲＯＭ（Ｅｌｅｃｔｒｉｃａｌｌｙ  Ｅｒａｓａｂｌｅ  Ｐｒｏｇｒａｍｍａｂｌｅ  Ｒｅａｄ  Ｏｎｌｙ  Ｍｅｍｏｒy）などの不揮発性又は揮発性の半導体メモリ、磁気ディスク、フレキシブルディスク、光ディスク、コンパクトディスク、ミニディスク、あるいは、ＤＶＤ（Ｄｉｇｉｔａｌ  Ｖｅｒｓａｔｉｌｅ  Ｄｉｓｃ）が該当する。
また、信号受信回路２１、フーリエ変換回路２２、立ち上がり検出回路２３、パラメータ設定回路２４、特徴量算出回路２５及び無線機識別回路２７は、例えば、単一回路、複合回路、プログラム化したプロセッサ、並列プログラム化したプロセッサ、ＡＳＩＣ（Ａｐｐｌｉｃａｔｉｏｎ  Ｓｐｅｃｉｆｉｃ  Ｉｎｔｅｇｒａｔｅｄ  Ｃｉｒｃｕｉｔ）、ＦＰＧＡ（Ｆｉｅｌｄ−Ｐｒｏｇｒａｍｍａｂｌｅ  Ｇａｔｅ  Ａｒｒａｙ）、または、これらを組み合わせたものが該当する。
【００２１】
無線機識別装置における信号受信部１を除く構成要素は、専用のハードウェアで実現されるものに限るものではなく、無線機識別装置における信号受信部１を除く構成要素が、ソフトウェア、ファームウェア、または、ソフトウェアとファームウェアとの組み合わせで実現されるものであってもよい。
ソフトウェア又はファームウェアはプログラムとして、コンピュータのメモリに格納される。コンピュータは、プログラムを実行するハードウェアを意味し、例えば、ＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）、中央処理装置、処理装置、演算装置、マイクロプロセッサ、マイクロコンピュータ、プロセッサ、あるいは、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）が該当する。
【００２２】
図３は、無線機識別装置における信号受信部１を除く構成要素がソフトウェア又はファームウェアなどで実現される場合のコンピュータのハードウェア構成図である。
無線機識別装置における信号受信部１を除く構成要素がソフトウェア又はファームウェアなどで実現される場合、データベース部８をコンピュータのメモリ３１上に構成するとともに、フーリエ変換部２、立ち上がり検出部３、パラメータ設定部４、特徴量算出部７及び無線機識別部１１の処理手順をコンピュータに実行させるためのプログラムをメモリ３１に格納し、コンピュータのプロセッサ３２がメモリ３１に格納されているプログラムを実行するようにすればよい。
【００２３】
また、図２では、無線機識別装置の構成要素のそれぞれが専用のハードウェアで実現される例を示し、図３では、無線機識別装置がソフトウェアやファームウェアなどで実現される例を示しているが、無線機識別装置における一部の構成要素が専用のハードウェアで実現され、残りの構成要素がソフトウェアやファームウェアなどで実現されるものであってもよい。
【００２４】
次に動作について説明する。
＜学習時の処理＞
最初に、既知の無線機から送信された無線信号を学習する際の処理内容を説明する。
図４は、この発明の実施の形態１による無線機識別方法における学習処理を示すフローチャートである。
信号受信部１は、既知の無線機から送信された無線信号を学習用無線信号として受信し、受信した学習用無線信号をアナログ信号からデジタル信号に変換する（図４のステップＳＴ１）。
学習用無線信号には、既知の無線機を識別する個体情報であるシリアルナンバーが付加されているものとする。
信号受信部１は、デジタルの学習用無線信号ｓｇ（ｎ）をフーリエ変換部２及び第１のデータベース部９に出力する。ｓｇ（ｎ）におけるｎは、学習用無線信号のサンプル番号である。例えば、図９には、サンプル番号ｎが１〜４の学習用無線信号が示されている。
第１のデータベース部９には、シリアルナンバーが付加されている学習用無線信号ｓｇ（ｎ）が格納される。
図９は、第１のデータベース部９に格納されているシリアルナンバーが付加されている学習用無線信号を示す説明図である。
図９の例では、シリアルナンバーが同じ無線機であっても、受信時のＳＮＲ又はサンプリング周波数などの受信条件が異なっている複数の学習用無線信号が格納されている。
【００２５】
フーリエ変換部２は、以下の式（１）に示すように、信号受信部１から出力されたデジタルの学習用無線信号ｓｇ（ｎ）をＳＴＦＴし、学習用無線信号ｓｇ（ｎ）のＳＴＦＴ結果を示すフーリエ変換信号である学習用フーリエ変換信号Ｓｇ（ｍ，ｋ）を立ち上がり検出部３及び特徴量算出部７に出力する（図４のステップＳＴ２）。
0006567242.tif 000002
【００２６】
式（１）において、ｗ（ｎ−ｍ）は、窓関数であり、例えば、矩形窓又はハミング窓が考えられる。
ｍは、窓関数ｗ（ｎ−ｍ）をシフトさせるサンプル数、ｍ０は、窓関数ｗ（ｎ−ｍ）の適用区間の初めのサンプル番号、Ｍは、窓関数ｗ（ｎ−ｍ）のサンプル数、ｋは、周波数ビン番号である。
なお、窓関数ｗ（ｎ−ｍ）のサンプル数Ｍ及び窓関数ｗ（ｎ−ｍ）をシフトさせるサンプル数ｍのそれぞれは、第２のパラメータ設定部６によって事前に設定される。
【００２７】
ここで、図５は、フーリエ変換部２による学習用無線信号ｓｇ（ｎ）及び識別用無線信号ｓｄ（ｎ）のＳＴＦＴを示す説明図である。
学習用無線信号ｓｇ（ｎ）及び識別用無線信号ｓｄ（ｎ）は、時間と振幅との関係を示す２次元信号である。識別用無線信号ｓｄ（ｎ）については後述する。
ＳＴＦＴ結果を示す学習用フーリエ変換信号Ｓｇ（ｍ，ｋ）及び識別用フーリエ変換信号Ｓｄ（ｍ，ｋ）は、時間と周波数及び電力との関係を示す３次元信号である。識別用フーリエ変換信号Ｓｄ（ｍ，ｋ）については後述する。
フーリエ変換部２により学習用無線信号ｓｇ（ｎ）及び識別用無線信号ｓｄ（ｎ）のそれぞれがＳＴＦＴされることで、学習用無線信号ｓｇ（ｎ）及び識別用無線信号ｓｄ（ｎ）のそれぞれが、コヒーレント積分されるので、ＳＮＲが改善される効果が得られる。
【００２８】
立ち上がり検出部３は、学習用無線信号ｓｇ（ｎ）の立ち上がり時刻Ｔ０を検出する処理を実施する（図４のステップＳＴ３）。
図６は、立ち上がり検出部３による立ち上がり時刻Ｔ０の検出処理を示す説明図である。
以下、立ち上がり検出部３による立ち上がり時刻Ｔ０の検出処理を具体的に説明する。
【００２９】
まず、立ち上がり検出部３は、以下の式（２）に示すように、フーリエ変換部２から出力された学習用フーリエ変換信号Ｓｇ（ｍ，ｋ）の２乗値｜Ｓｇ（ｍ，ｋ）｜２をスペクトログラムＳＰｇ（ｍ，ｋ）として算出する。
0006567242.tif 000003
【００３０】
次に、立ち上がり検出部３は、以下の式（３）に示すように、スペクトログラムＳＰｇ（ｍ，ｋ）におけるエネルギーの時間波形である学習用のエネルギー時間波形Ｅｇ（ｍ）を算出する。
学習用のエネルギー時間波形Ｅｇ（ｍ）は、スペクトログラムＳＰｇ（ｍ，ｋ）に含まれている複数の周波数における電力のうち、同一時刻の電力同士がそれぞれ合計されることで算出される各々の時刻における電力の合計値である。
0006567242.tif 000004
式（３）において、ｋ１は、スペクトログラムＳＰｇ（ｍ，ｋ）に含まれている複数の周波数のうち、電力を合計する周波数の下限周波数を示す周波数ビン番号、ｋｎは、スペクトログラムＳＰｇ（ｍ，ｋ）に含まれている複数の周波数のうち、電力を合計する周波数の上限周波数を示す周波数ビン番号である。
ｋ＝ｋ１，・・・，ｋｎである。
【００３１】
次に、立ち上がり検出部３は、学習用のエネルギー時間波形Ｅｇ（ｍ）と閾値Ｅｔｈとを比較する。
閾値Ｅｔｈとしては、例えば、学習用のエネルギー時間波形Ｅｇ（ｍ）の最大値から設定値Ｅ０を減算した値が考えられる。
立ち上がり検出部３は、学習用のエネルギー時間波形Ｅｇ（ｍ）と閾値Ｅｔｈとの比較結果を参照することで、最初は、ノイズレベルのエネルギーを示していた学習用のエネルギー時間波形Ｅｇ（ｍ）が、閾値Ｅｔｈに到達したか否かを判定する。
立ち上がり検出部３は、学習用のエネルギー時間波形Ｅｇ（ｍ）が閾値Ｅｔｈに到達した時刻を、学習用無線信号ｓｇ（ｎ）の立ち上がり時刻Ｔ０として検出する。
この実施の形態１では、立ち上がり検出部３が、学習用無線信号ｓｇ（ｎ）の立ち上がり時刻Ｔ０を検出する際、スペクトログラムＳＰｇ（ｍ，ｋ）に含まれている複数の周波数における電力のうち、同一時刻の電力同士をそれぞれ合計するようにしている。これにより、立ち上がり時に電力が大きく変動する周波数が含まれている場合でも、当該電力の変動の影響が低減されるため、立ち上がり時刻Ｔ０の検出精度が向上する。
【００３２】
特徴量算出部７は、第１のパラメータ設定部５により設定された特徴量算出用のパラメータから、立ち上がり時刻Ｔ０を含む時間帯（ｔ１〜ｔ２）及び設定周波数ビンｋｂを認識する。
第１のパラメータ設定部５による特徴量算出用のパラメータの設定処理については後述する。
次に、特徴量算出部７は、図６に示すように、フーリエ変換部２から出力された学習用フーリエ変換信号Ｓｇ（ｍ，ｋ）の中から、立ち上がり検出部３により検出された学習用無線信号ｓｇ（ｎ）の立ち上がり時刻Ｔ０を含む時間帯（ｔ１〜ｔ２）の学習用フーリエ変換信号ｓｇ（ｎ）ｔ１〜ｔ２を抽出する。
【００３３】
次に、特徴量算出部７は、図７に示すように、抽出した学習用フーリエ変換信号ｓｇ（ｎ）ｔ１〜ｔ２に含まれている複素時間信号の中から、設定周波数ビンｋｂの複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２の切り出しを行う。
電子出願の都合上、明細書の文章中では、ｓの文字の上に“＾”の記号を付することができないので、ｓｇハット（ｎ）ｔ１〜ｔ２のように表記している。
図７は、特徴量算出部７による特徴量の算出処理を示す説明図である。
特徴量算出部７は、切り出した複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２から、抽出した学習用フーリエ変換信号ｓｇ（ｎ）ｔ１〜ｔ２の変動を特徴量として算出する（図４のステップＳＴ４）。
【００３４】
以下、特徴量算出部７による特徴量の算出処理を具体的に説明する。
特徴量算出部７は、切り出した複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２から、以下の式（４）に示すように、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２の瞬時振幅ａｇ（ｎ）を算出する。
0006567242.tif 000005
式（４）において、Ｉｇ（ｎ）は、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２の実部、Ｑｇ（ｎ）は、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２の虚部である。
【００３５】
また、特徴量算出部７は、切り出した複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２から、以下の式（５）に示すように、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２の瞬時位相φｇ（ｎ）を算出する。
0006567242.tif 000006
また、特徴量算出部７は、切り出した複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２から、以下の式（６）に示すように、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２の瞬時周波数ｆｇ（ｎ）を算出する。
0006567242.tif 000007
【００３６】
次に、特徴量算出部７は、以下の式（７）に示すように、第１の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時振幅ａｇ（ｎ）の分散値ｐｇ１を算出する。
0006567242.tif 000008
0006567242.tif 000009
特徴量算出部７は、以下の式（９）に示すように、第２の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時振幅ａｇ（ｎ）の歪度ｐｇ２を算出する。
0006567242.tif 000010
0006567242.tif 000011
特徴量算出部７は、以下の式（１１）に示すように、第３の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時振幅ａｇ（ｎ）の尖度ｐｇ３を算出する。
0006567242.tif 000012
【００３７】
次に、特徴量算出部７は、以下の式（１２）に示すように、第４の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時位相φｇ（ｎ）の分散値ｐｇ４を算出する。
0006567242.tif 000013
0006567242.tif 000014
特徴量算出部７は、以下の式（１４）に示すように、第５の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時位相φｇ（ｎ）の歪度ｐｇ５を算出する。
0006567242.tif 000015
0006567242.tif 000016
特徴量算出部７は、以下の式（１６）に示すように、第６の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時位相φｇ（ｎ）の尖度ｐｇ６を算出する。
0006567242.tif 000017
【００３８】
次に、特徴量算出部７は、以下の式（１７）に示すように、第７の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時周波数ｆｇ（ｎ）の分散値ｐｇ７を算出する。
0006567242.tif 000018
0006567242.tif 000019
特徴量算出部７は、以下の式（１９）に示すように、第８の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時周波数ｆｇ（ｎ）の歪度ｐｇ８を算出する。
0006567242.tif 000020
0006567242.tif 000021
特徴量算出部７は、以下の式（２１）に示すように、第９の特徴量として、複素時間信号ｓｇハット（ｎ）ｔ１〜ｔ２における瞬時周波数ｆｇ（ｎ）の尖度ｐｇ９を算出する。
0006567242.tif 000022
【００３９】
特徴量算出部７は、第２のパラメータ設定部６により設定された第１から第９の特徴量ｐｇｉ（ｉ＝１，２，・・・，９）の重みｗｉを示す重み付けパラメータを取得する。
第１から第９の特徴量ｐｇｉの中には、無線機の識別に有効な特徴量と、識別に有効でない特徴量とが存在している。このため、識別に有効な特徴量には大きな重みが設定され、識別に有効でない特徴量には小さな重みが設定される。
例えば、ｗｉ（ｉ＝１，２，・・・，９）は、０以上１以下の値であり、ｗ１＋ｗ２＋・・・＋ｗ９＝１である。
第２のパラメータ設定部６による重み付けパラメータの設定方法は、特に問わないが、例えば、第１から第９の特徴量ｐｇｉの分散値で、第１から第９の特徴量ｐｇｉを規格化し、第１から第９の特徴量ｐｇｉを規格化した値を重みとして用いる方法が考えられる。
【００４０】
特徴量算出部７は、以下の式（２２）に示すように、第１から第９の特徴量ｐｇｉと重みｗｉとをそれぞれ乗算することで、第１から第９の特徴量ｐｇｉの重み付けを行う。
ｐ’ｇｉ＝ｐｇｉ×ｗｉ                      （２２）
ｉ＝１，２，・・・，９
特徴量算出部７は、重み付け後の第１から第９の特徴量ｐ’ｇｉを並べたベクトルを、既知の無線機の特徴量ベクトルｐｇとして算出する。
特徴量算出部７は、算出した既知の無線機の特徴量ベクトルｐｇを、学習用無線信号に付加されているシリアルナンバーと一緒に、第２のデータベース部１０に格納する（図４のステップＳＴ５）。
また、特徴量算出部７は、既知の無線機の特徴量ベクトルｐｇを第２のデータベース部１０に格納した旨を示す情報を第１のパラメータ設定部５に出力する。
【００４１】
ここで、第１のパラメータ設定部５による特徴量算出用のパラメータの設定処理について説明する。
図８は、特徴量算出用のパラメータと複数の既知の無線機に係る特徴量ベクトルとの関係を示す説明図である。
図８の例では、〇は、種類（Ａ）の既知の無線機に係る特徴量ベクトルｐｇ、◇は、種類（Ｂ）の既知の無線機に係る特徴量ベクトルｐｇ、△は、種類（Ｃ）の既知の無線機に係る特徴量ベクトルｐｇである。
図８の例では、種類（Ａ）に属する無線機の個体の数が１０であるため、１０個の特徴量ベクトル〇が得られている。
また、種類（Ｂ）に属する無線機の個体の数が１０であるため、１０個の特徴量ベクトル◇が得られており、種類（Ｃ）に属する無線機の個体の数が１０であるため、１０個の特徴量ベクトル△が得られている。
【００４２】
図８Ａは、特徴量算出用のパラメータが適正でないために、種類（Ａ）の無線機に係る特徴量ベクトル〇と、種類（Ｂ）の無線機に係る特徴量ベクトル◇と、種類（Ｃ）の無線機に係る特徴量ベクトル△とが一部重なっている例を示している。
既知の無線機に係る特徴量ベクトル〇と特徴量ベクトル◇と特徴量ベクトル△とが重なっている場合、これらの特徴量ベクトル〇，◇，△と、識別対象の無線機に係る特徴量ベクトルとを比較しても、特徴量ベクトル〇，◇，△の中で、識別対象の無線機に係る特徴量ベクトルと最も類似している特徴量ベクトルを高精度に特定することができない。
したがって、特徴量ベクトル〇，◇，△の中で、識別対象の無線機に係る特徴量ベクトルと最も類似している特徴量ベクトルを高精度に特定する上で、既知の無線機に係る特徴量ベクトル〇と特徴量ベクトル◇と特徴量ベクトル△とが重なっていないことが望ましい。
図８Ｂは、特徴量算出用のパラメータが適正であるために、種類（Ａ）の無線機に係る特徴量ベクトル〇と、種類（Ｂ）の無線機に係る特徴量ベクトル◇と、種類（Ｃ）の無線機に係る特徴量ベクトル△とが重なっていない例を示している。
【００４３】
第１のパラメータ設定部５は、特徴量算出用のパラメータの初期設定として、時間帯（ｔ１〜ｔ２）の長さである時間範囲Δｔを示すパラメータと、設定周波数ビンｋｂを示すパラメータとを任意の値に設定する。
時間帯（ｔ１〜ｔ２）については、学習用のエネルギー時間波形Ｅｇ（ｍ）に基づいて設定する方法が考えられる。
例えば、学習用無線信号ｓｇ（ｎ）の立ち上がり時刻Ｔ０よりも前の時刻において、学習用のエネルギー時間波形Ｅｇ（ｍ）の最大値よりもエネルギーが３０ｄＢ下がっている時刻をｔ１とする。また、学習用無線信号ｓｇ（ｎ）の立ち上がり時刻Ｔ０よりも後の時刻において、学習用のエネルギー時間波形Ｅｇ（ｍ）の最大値よりもエネルギーが３０ｄＢ下がっている時刻をｔ２とする。
また、ｔ２については、学習用フーリエ変換信号ｓｇ（ｎ）の周波数方向の変動に基づいて設定する方法が考えられる。
例えば、学習用フーリエ変換信号ｓｇ（ｎ）の過渡応答時の周波数変動が、特定の周波数範囲内に収まった時刻をｔ２とする。
【００４４】
設定周波数ビンｋｂについては、学習用のエネルギー時間波形Ｅｇ（ｍ）において、平均電力が最大値をとる周波数ビンを採用する方法が考えられる。
また、設定周波数ビンｋｂについては、学習用のエネルギー時間波形Ｅｇ（ｍ）において、平均電力が閾値以上の周波数ビンを検出して、検出した１以上の周波数ビンに係る周波数の平均を求め、平均の周波数の周波数ビンを採用する方法が考えられる。
【００４５】
第１のパラメータ設定部５は、特徴量算出部７から既知の無線機の特徴量ベクトルｐｇを第２のデータベース部１０に格納した旨を示す情報を受けると、第２のデータベース部１０から、複数の既知の無線機の特徴量ベクトルｐｇを取得する。
【００４６】
第１のパラメータ設定部５は、複数の既知の無線機に係る特徴量ベクトルｐｇの間に重なりがあるか否かを判定し、例えば、図８Ｂに示すように、複数の既知の無線機に係る特徴量ベクトルｐｇの間に重なりがなければ、先に設定している特徴量算出用のパラメータを有効として、パラメータの設定処理を終了する。
第１のパラメータ設定部５は、例えば、図８Ａに示すように、複数の既知の無線機に係る特徴量ベクトルｐｇの間に重なりがあれば、特徴量算出用のパラメータとして、先に設定している時間範囲Δｔ又は設定周波数ビンｋｂのうち、少なくとも１つを更新する。
特徴量算出部７は、第１のパラメータ設定部５により更新された特徴量算出用のパラメータを用いて、再度、複数の既知の無線機の特徴量ベクトルｐｇをそれぞれ算出する。
第１のパラメータ設定部５は、特徴量算出部７により再度算出された複数の既知の無線機に係る特徴量ベクトルｐｇの間に重なりがあるか否かを判定する。
複数の既知の無線機に係る特徴量ベクトルｐｇの間に重なりがなくなるまで、第１のパラメータ設定部５による特徴量算出用のパラメータの更新処理と、特徴量算出部７による特徴量ベクトルｐｇの算出処理とが繰り返される。
なお、パラメータの更新処理と、特徴量ベクトルｐｇの算出処理とを繰り返しても、複数の既知の無線機に係る重なりがなくならない場合には、重なりが最も小さくなったパラメータを適正なパラメータとして、パラメータの設定処理を終了する。
【００４７】
第１のパラメータ設定部５により特徴量算出用のパラメータが更新されることで、特徴量算出用のパラメータが適正な値になると、図８Ｂに示すように、第２のデータベース部１０に格納される複数の既知の無線機に係る特徴量ベクトルｐｇの間の重なりがなくなる、もしくは、重なりが最も小さくなる。
特徴量算出用のパラメータが更新されたのち、特徴量算出部７により再度算出された複数の既知の無線機に係る特徴量ベクトルｐｇが第２のデータベース部１０に格納される際、不適正な特徴量算出用のパラメータに従って算出されている複数の既知の無線機に係る特徴量ベクトルｐｇは、第２のデータベース部１０から削除される。
図１０は、第２のデータベース部１０に格納される複数の既知の無線機に係る特徴量の一例を示す説明図である。
図１０では、既知の無線機の種類が、種類（Ａ）、種類（Ｂ）の例を示しており、また、種類（Ａ）に属する無線機が、シリアルナンバー（１）の個体、シリアルナンバー（２）の個体の例を示している。
また、図１０には、特徴量算出用のパラメータ及び重み付けパラメータの一例を示している。
【００４８】
この実施の形態１では、第１のパラメータ設定部５によって、特徴量算出用のパラメータが適正な値になるまで、特徴量算出用のパラメータが更新される例を示しているが、特徴量算出用のパラメータの更新タイミングは、これに限るものではない。
例えば、信号受信部１によって、既知の無線機から新たに送信された学習用無線信号が受信されることで、特徴量算出部７によって、新たに特徴量が算出されたタイミングで、第１のパラメータ設定部５が、特徴量算出用のパラメータを更新するようにしてもよい。
【００４９】
＜識別時の処理＞
次に、識別対象の無線機から送信された無線信号を用いて、無線機を識別する際の処理内容を説明する。
図１１は、この発明の実施の形態１による無線機識別方法における識別処理を示すフローチャートである。
信号受信部１は、識別対象の無線機から送信された無線信号を識別用無線信号として受信し、受信した識別用無線信号をアナログ信号からデジタル信号に変換する（図１１のステップＳＴ１１）。
信号受信部１は、デジタルの識別用無線信号ｓｄ（ｎ）をフーリエ変換部２に出力する。ｓｄ（ｎ）におけるｎは、識別用無線信号のサンプル番号である。
【００５０】
フーリエ変換部２は、以下の式（２３）に示すように、信号受信部１から出力されたデジタルの識別用無線信号ｓｄ（ｎ）をＳＴＦＴする（図１１のステップＳＴ１２）。
フーリエ変換部２は、識別用無線信号ｓｄ（ｎ）のＳＴＦＴ結果を示すフーリエ変換信号である識別用フーリエ変換信号Ｓｄ（ｍ，ｋ）を立ち上がり検出部３及び特徴量算出部７に出力する（図１１のステップＳＴ１２）。図５は、識別用無線信号ｓｄ（ｎ）のＳＴＦＴについても表している。
0006567242.tif 000023
立ち上がり検出部３は、識別用無線信号ｓｄ（ｎ）の立ち上がり時刻Ｔ０を検出する処理を実施する（図１１のステップＳＴ１３）。
以下、立ち上がり検出部３による立ち上がり時刻Ｔ０の検出処理を具体的に説明する。
【００５１】
まず、立ち上がり検出部３は、以下の式（２４）に示すように、フーリエ変換部２から出力された識別用フーリエ変換信号Ｓｄ（ｍ，ｋ）の２乗値｜Ｓｄ（ｍ，ｋ）｜２をスペクトログラムＳＰｄ（ｍ，ｋ）として算出する。
0006567242.tif 000024
【００５２】
次に、立ち上がり検出部３は、以下の式（２５）に示すように、スペクトログラムＳＰｄ（ｍ，ｋ）におけるエネルギーの時間波形である識別用のエネルギー時間波形Ｅｄ（ｍ）を算出する。
識別用のエネルギー時間波形Ｅｄ（ｍ）は、スペクトログラムＳＰｄ（ｍ，ｋ）に含まれている複数の周波数における電力のうち、同一時刻の電力同士がそれぞれ合計されることで算出される各々の時刻における電力の合計値である。
0006567242.tif 000025
式（２５）において、ｋ１は、スペクトログラムＳＰｄ（ｍ，ｋ）に含まれている複数の周波数のうち、電力を合計する周波数の下限周波数を示す周波数ビン番号、ｋｎは、スペクトログラムＳＰｄ（ｍ，ｋ）に含まれている複数の周波数のうち、電力を合計する周波数の上限周波数を示す周波数ビン番号である。
ｋ＝ｋ１，・・・，ｋｎである。
【００５３】
次に、立ち上がり検出部３は、識別用のエネルギー時間波形Ｅｄ（ｍ）と閾値Ｅｔｈとを比較する。
立ち上がり検出部３は、識別用のエネルギー時間波形Ｅｄ（ｍ）と閾値Ｅｔｈとの比較結果を参照することで、最初は、ノイズレベルのエネルギーを示していた識別用のエネルギー時間波形Ｅｄ（ｍ）が、閾値Ｅｔｈに到達したか否かを判定する。
立ち上がり検出部３は、識別用のエネルギー時間波形Ｅｄ（ｍ）が閾値Ｅｔｈに到達した時刻を、識別用無線信号ｓｄ（ｎ）の立ち上がり時刻Ｔ０として検出する。
この実施の形態１では、立ち上がり検出部３が、識別用無線信号ｓｄ（ｎ）の立ち上がり時刻Ｔ０を検出する際、スペクトログラムＳＰｄ（ｍ，ｋ）に含まれている複数の周波数における電力のうち、同一時刻の電力同士をそれぞれ合計するようにしている。これにより、立ち上がり時に電力が大きく変動する周波数が含まれている場合でも、当該電力の変動の影響が低減されるため、立ち上がり時刻Ｔ０の検出精度が向上する。
【００５４】
特徴量算出部７は、第１のパラメータ設定部５により設定された特徴量算出用のパラメータから、立ち上がり時刻Ｔ０を含む時間帯（ｔ１〜ｔ２）及び設定周波数ビンｋｂを認識する。
この特徴量算出用のパラメータは、第１のパラメータ設定部５によって、複数の既知の無線機に係る特徴量ベクトルｐｇの間に重なりがなくなるように、もしくは、重なりが最も小さくなるように設定されたパラメータである。
次に、特徴量算出部７は、図６に示すように、フーリエ変換部２から出力された識別用フーリエ変換信号Ｓｄ（ｍ，ｋ）の中から、立ち上がり検出部３により検出された識別用無線信号ｓｄ（ｎ）の立ち上がり時刻Ｔ０を含む時間帯（ｔ１〜ｔ２）の識別用フーリエ変換信号ｓｄ（ｎ）ｔ１〜ｔ２を抽出する。
【００５５】
次に、特徴量算出部７は、図７に示すように、抽出した識別用フーリエ変換信号ｓｄ（ｎ）ｔ１〜ｔ２に含まれている複素時間信号の中から、設定周波数ビンｋｂの複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２の切り出しを行う。
電子出願の都合上、明細書の文章中では、ｓの文字の上に“＾”の記号を付することができないので、ｓｄハット（ｎ）ｔ１〜ｔ２のように表記している。
特徴量算出部７は、切り出した複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２から、抽出した識別用フーリエ変換信号ｓｄ（ｎ）ｔ１〜ｔ２の変動を特徴量として算出する（図１１のステップＳＴ１４）。
【００５６】
以下、特徴量算出部７による特徴量の算出処理を具体的に説明する。
特徴量算出部７は、切り出した複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２から、以下の式（２６）に示すように、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２の瞬時振幅ａｄ（ｎ）を算出する。
0006567242.tif 000026
式（２６）において、Ｉｄ（ｎ）は、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２の実部、Ｑｄ（ｎ）は、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２の虚部である。
【００５７】
また、特徴量算出部７は、切り出した複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２から、以下の式（２７）に示すように、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２の瞬時位相φｄ（ｎ）を算出する。
0006567242.tif 000027
また、特徴量算出部７は、切り出した複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２から、以下の式（２８）に示すように、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２の瞬時周波数ｆｄ（ｎ）を算出する。
0006567242.tif 000028
【００５８】
次に、特徴量算出部７は、以下の式（２９）に示すように、第１の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時振幅ａｄ（ｎ）の分散値ｐｄ１を算出する。
0006567242.tif 000029
0006567242.tif 000030
特徴量算出部７は、以下の式（３１）に示すように、第２の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時振幅ａｄ（ｎ）の歪度ｐｄ２を算出する。
0006567242.tif 000031
0006567242.tif 000032
特徴量算出部７は、以下の式（３３）に示すように、第３の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時振幅ａｄ（ｎ）の尖度ｐｄ３を算出する。
0006567242.tif 000033
【００５９】
次に、特徴量算出部７は、以下の式（３４）に示すように、第４の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時位相φｄ（ｎ）の分散値ｐｄ４を算出する。
0006567242.tif 000034
0006567242.tif 000035
特徴量算出部７は、以下の式（３６）に示すように、第５の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時位相φｄ（ｎ）の歪度ｐｄ５を算出する。
0006567242.tif 000036
0006567242.tif 000037
特徴量算出部７は、以下の式（３８）に示すように、第６の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時位相φｄ（ｎ）の尖度ｐｄ６を算出する。
0006567242.tif 000038
【００６０】
次に、特徴量算出部７は、以下の式（３９）に示すように、第７の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時周波数ｆｄ（ｎ）の分散値ｐｄ７を算出する。
0006567242.tif 000039
0006567242.tif 000040
特徴量算出部７は、以下の式（４１）に示すように、第８の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時周波数ｆｄ（ｎ）の歪度ｐｄ８を算出する。
0006567242.tif 000041
0006567242.tif 000042
特徴量算出部７は、以下の式（４３）に示すように、第９の特徴量として、複素時間信号ｓｄハット（ｎ）ｔ１〜ｔ２における瞬時周波数ｆｄ（ｎ）の尖度ｐｄ９を算出する。
0006567242.tif 000043
【００６１】
特徴量算出部７は、第２のパラメータ設定部６により設定された第１から第９の特徴量ｐｄｉ（ｉ＝１，２，・・・，９）の重みｗｉを示す重み付けパラメータを取得する。
特徴量算出部７は、以下の式（４４）に示すように、第１から第９の特徴量ｐｄｉと重みｗｉとをそれぞれ乗算することで、第１から第９の特徴量ｐｄｉの重み付けを行う。
ｐ’ｄｉ＝ｐｄｉ×ｗｉ    （４４）
ｉ＝１，２，・・・，９
特徴量算出部７は、重み付け後の第１から第９の特徴量ｐ’ｄｉを並べたベクトルを、識別対象の無線機の特徴量ベクトルｐｄとして算出する。
特徴量算出部７は、算出した識別対象の無線機の特徴量ベクトルｐｄを無線機識別部１１に出力する。
【００６２】
無線機識別部１１は、特徴量算出部７により算出された識別対象の無線機に係る特徴量ベクトルｐｄと、第２のデータベース部１０により格納されている複数の既知の無線機に係る特徴量ベクトルｐｇとをそれぞれ比較する。
そして、無線機識別部１１は、識別対象の無線機に係る特徴量ベクトルｐｄと、複数の既知の無線機に係る特徴量ベクトルｐｇとの比較結果に基づいて、識別対象の無線機を識別する（図１１のステップＳＴ１５）。
【００６３】
以下、無線機識別部１１による無線機の識別処理を具体的に説明する。
図１２は、無線機識別部１１による無線機の識別処理を示す説明図である。
図１２では、説明の簡単化のため、特徴量算出部７により算出される第１から第９の特徴量のうち、第１から第３の特徴量だけを例示している。
【００６４】
無線機識別部１１は、以下の式（４５）に示すように、識別対象の無線機に係る特徴量ベクトルｐｄと、第２のデータベース部１０により格納されている複数の既知の無線機に係る特徴量ベクトルｐｇとの間のマハラノビス距離ｄｉｊ，ｈをそれぞれ算出する。
ここでは、説明の便宜上、既知の無線機がＨ個あり、Ｈ個の無線機に係る特徴量ベクトルｐｇが第２のデータベース部１０にそれぞれ格納されているものとする。
Ｈ個の無線機の中には、例えば、種類（Ａ）に属する無線機、種類（Ｂ）に属する無線機、あるいは、種類（Ｃ）に属する無線機が含まれている。
0006567242.tif 000044
式（４５）において、ｘｉ，ｈは、ｎ次元の特徴量空間における既知の無線機ｈ（ｈ＝１，・・・，Ｈ）に係る特徴量ベクトルｐｇの位置、ｘｊは、ｎ次元の特徴量空間における識別対象の無線機に係る特徴量ベクトルｐｄの位置である。
（ｘｉ，ｈ−ｘｊ）は、ｎ次元の特徴量空間における２点間のユークリッド距離、Ｓは、既知の無線機ｈから送信される学習用無線信号ｓｇ（ｎ）の特徴量空間における分布の分散共分散行列である。
マハラノビス距離ｄｉｊ，ｈは、学習用無線信号ｓｇ（ｎ）の特徴量空間における分布の分散による影響が正規化された空間での距離となる。
【００６５】
無線機識別部１１は、算出したＨ個のマハラノビス距離ｄｉｊ，ｈを相互に比較し、Ｈ個のマハラノビス距離ｄｉｊ，ｈの中で、最小のマハラノビス距離ｄｉｊ，ｈを特定する。
無線機識別部１１は、Ｈ個の無線機の中で、最小のマハラノビス距離ｄｉｊ，ｈに係る既知の無線機が、識別対象の無線機と最も類似していると判断する。
無線機識別部１１は、識別対象の無線機の識別結果として、最小のマハラノビス距離ｄｉｊ，ｈに係る既知の無線機のシリアルナンバーを出力する。これにより、無線機の種類だけでなく、無線機の個体識別が可能になる。
図１２の例では、３つの既知の無線機に係る特徴量ベクトルｐｇの位置が示されている。３つの既知の無線機の中で、特徴量ベクトルｐｇの位置が●である既知の無線機とのマハラノビス距離が最小であるため、識別対象の無線機は、特徴量ベクトルｐｇの位置が●の無線機であると識別される。
【００６６】
ここでは、無線機識別部１１が、マハラノビス距離ｄｉｊ，ｈを用いて、無線機を識別する例を示しているが、これに限るものではなく、例えば、最小二乗法による線形識別、サポートベクターマシン（ＳＶＭ：Ｓｕｐｐｏｒｔ  Ｖｅｃｔｏｒ  Ｍａｃｈｉｎｅ）、あるいは、ニューラルネットワークを用いて、無線機を識別するようにしてもよい。
【００６７】
以上で明らかなように、この実施の形態１によれば、フーリエ変換部２から出力されたフーリエ変換信号のスペクトログラムにおけるエネルギーの時間波形を算出し、エネルギーの時間波形から無線信号の立ち上がり時刻を検出する立ち上がり検出部３と、フーリエ変換部２から出力されたフーリエ変換信号の中から、立ち上がり検出部３により検出された立ち上がり時刻を含む時間帯のフーリエ変換信号を抽出し、抽出したフーリエ変換信号に含まれている設定周波数成分の複素時間信号から、抽出したフーリエ変換信号の変動を特徴量として算出する特徴量算出部７とを設け、無線機識別部１１が、特徴量算出部７により算出された特徴量に基づいて、無線機を識別するように構成したので、無線機の個体を識別することができる効果を奏する。
【００６８】
なお、本願発明はその発明の範囲内において、実施の形態の任意の構成要素の変形、もしくは実施の形態の任意の構成要素の省略が可能である。
【産業上の利用可能性】
【００６９】
この発明は、無線機を識別する無線機識別装置及び無線機識別方法に適している。
【符号の説明】
【００７０】
１  信号受信部、２  フーリエ変換部、３  立ち上がり検出部、４  パラメータ設定部、５  第１のパラメータ設定部、６  第２のパラメータ設定部、７  特徴量算出部、８  データベース部、９  第１のデータベース部、１０  第２のデータベース部、１１  無線機識別部、２１  信号受信回路、２２  フーリエ変換回路、２３  立ち上がり検出回路、２４  パラメータ設定回路、２５  特徴量算出回路、２６  データベース回路、２７  無線機識別回路、３１  メモリ、３２  プロセッサ。
【図１】
0006567242.tif 000045
【図２】
0006567242.tif 000046
【図３】
0006567242.tif 000047
【図４】
0006567242.tif 000048
【図５】
0006567242.tif 000049
【図６】
0006567242.tif 000050
【図７】
0006567242.tif 000051
【図８】
0006567242.tif 000052
【図９】
0006567242.tif 000053
【図１０】
0006567242.tif 000054
【図１１】
0006567242.tif 000055
【図１２】
0006567242.tif 000056