(57)【要約】【課題】波形位置合わせにおいて、ノイズを含む波形に対して高精度に位置ずれ量を推定し、サンプリング間隔以下の位置ずれの量も検出できる、波形の位置合わせ方法を提供する。【解決手段】波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、位相合わせ手段は、入力波形信号及び参照波形信号それぞれを離散フーリエ変換して位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し位相限定相関関数を算出し、波形間の位置ずれ量を検出する。【選択図】図２ａ
2007240170.tif 000002
【特許請求の範囲】
【請求項１】
波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、
前記位相合わせ手段は、
外部からの、入力波形信号を離散フーリエ変換し算出するステップ及び参照波形信号を離散フーリエ変換し算出するステップとを有し、該入力波形信号及び参照波形信号をそれぞれ離散フーリエ変換した際に分離される位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し相関関数を検出することを特徴とする波形の位置合わせ方法。
【請求項２】
波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、
前記位相合わせ手段は、
外部からの、入力波形信号を離散フーリエ変換し算出するステップ、予め設定された参照波形信号を離散フーリエ変換し算出するステップとを有し、該入力波形信号及び参照波形信号をそれぞれ離散フーリエ変換した際に分離される位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し相関関数を検出し、前記相関関数に相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする波形の位置合わせ方法。
【請求項３】
前記位相合わせ手段は、フーリエ変換の性質から理論的に算出される相関ピークの形状を近似する相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする請求項２に記載の波形の位置合わせ方法。
【請求項４】
前記位相合わせ手段は、２次関数を用いた相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする請求項２に記載の波形の位置合わせ方法。
【請求項５】
前記位相合わせ手段は、ガウス関数を用いた相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする請求項２に記載の波形の位置合わせ方法。
【請求項６】
波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、
前記位相合わせ手段は、
外部からの、入力波形信号を離散フーリエ変換し算出するステップ及び参照波形信号を離散フーリエ変換し算出するステップあるいは外部からの、入力波形信号を離散フーリエ変換し算出するステップ及び予め設定された参照波形信号を離散フーリエ変換し算出するステップ、いづれかの入力波形信号及び参照波形信号をそれぞれ離散フーリエ変換した際に分離される位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し相関関数を検出し、前記相関関数に相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定し、前記位置ずれ量を用いて波形リサンプリングすることで該入力波形の位置ずれ補正した波形を算出することを特徴とする波形の位置合わせ方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ディジタルオシロスコープ等の電子計測機器から得られる波形データ解析処理に関し、特に位相限定相関もしくはそれに準ずる計算式により、サンプリング分解能を超える精度での波形データの位置合わせを可能とする、波形の位置合わせ方法に関するものである。
【背景技術】
【０００２】
近年、急速なインターネットや電子メールの活用、移動通信端末等の情報通信ネットワーク技術の発展にともない、高度情報化社会が実現されてきている。これは電子化された情報を情報通信ネットワークによって、自己の通信機器端末より相手の通信機器端末へ高速、且つ、大量に双方向にも通信できる技術である。この情報通信ネットワークによる利便性が高まり、産業の効率化、時間の短縮化などが図られ、社会構造の変化をも促すものとなっている。
【０００３】
一方、技術が進み、大多数の人達が簡便に利用できることと同時に別の課題を生じさせている。即ち、音声だけの通信システムの漏洩とは異なり、情報通信ネットワーク上では電子化された信号へ不正アクセスを行うことで認証情報等を引き出す犯罪が生じており、プライバシー保護、セキュリティ対策等充分なセキュリティ対策を立てる必要が生じている。そのため、秘密情報管理は重要な技術と位置付けられている。
【０００４】
従来、情報通信ネットワークにおいては、送受信を行う通信装置にセキュリティ対策を行うための手法が備えられていることが多く、一般的には送受信を行う装置には伝送信号に対してディジタル信号を用いた暗号化処理を施す場合が多い。この暗号化処理は送受信装置内、暗号装置である、LSI(集積回路)に読込まれた専用アプリケーションの実行により任意の伝送信号にその処理を施していた。しかし、送受信装置間に接続された情報通信ネットワーク上に何らかの手段を用いて送受信装置の暗号装置の信号動作等を物理的に観察することができるとなれば、サイドチャネル攻撃（side-channel attack）を受けることになり、送受装置内における秘密情報を盗聴・解読されることでセキュリティに対する安全性が確保できなくなる恐れがある。秘密情報を奪うサイドチャネル攻撃としては、システムが計算している処理時間を観察するタイミング攻撃(timing attack)、システムが計算している時の消費電力を統計処理することで秘密情報(秘密鍵やPIN等)を観察する電力解析攻撃(power analysis attack)、暗号装置から漏洩する電磁波を観察解析する電磁波解析攻撃等が明らかになっている。
【０００５】
今後、身の回りに存在する情報関連機器がネットワークを介して接続される、ユビキタス情報社会が現実のものとなれば、暗号装置の目的にした攻撃、または情報通信ネットワーク運用の脆弱性を狙った犯罪により重大な社会的脅威となる可能性が高く、特にサイドチャネル攻撃に対しての有効な暗号処理システムの開発が望まれている。
【０００６】
暗号処理システムの開発のためには、想定される暗号攻撃方法の潜在的な能力及び限界を詳細に評価する必要があり、暗号処理システムの安全性を評価・検証する技術の確立が不可欠であると言える。上術した信号動作の観察にはディジタル測定機器を用いることが必要であるため、ディジタル測定機器から取得されるデータを前提とした評価方式が必要となる。
【０００７】
従来、このディジタル測定機器から取得されるデータへの評価技術分野においては画像補正を行う先行技術が開示されている。例えば、特許文献１には、対象画像を高解像度化する際に用いられる参照画像との間の動きを推定する動き動作ステップと、対象画像標本化時に発生するエイリアス成分を動き量に基づいて推定し、エイリアス成分を用いて対象画像を高解像度化する高解像度ステップとを含み、対象画像とのこの対象画像に対して微小に位置ずれした参照画像とを用いてエイリアス予測を行うことによって、対象画像を高解像度化することが記載されている。
【０００８】
また、このディジタル測定機器から取得されるデータへの評価技術分野でのパターン検出する方法としての別の先行技術があり、例えば、特許文献２には、予め登録された登録画像を入力画像中から検索する画像処理におけるパターン検出方法等において、逆フーリエ変換等の処理を行うことなく、入力画像中から登録画像のパターン検出を行えることが記載されている。
【０００９】
【特許文献１】特開２００５−７２９３６号公報
【特許文献２】特開２００４−３２６７２８号公報
【非特許文献１】K. Takita、 T. Aoki、 Y. Sasaki、 T. Higuchi、 and K. Kobayashi、 “High-accuracy subpixel image registration based on phase-only correlation、” IEICE Transactions on Fundamentals of Electronics、 Communications and Computer Sciences、 Vol. E86-A、 No. 8、 pp. 1925-1934、 August 2003.
【発明の開示】
【発明が解決しようとする課題】
【００１０】
上述したように、情報関連機器がネットワークを介して他の情報関連機器に接続し通信を行う以上暗号装置は必要であり、サイドチャネル攻撃のうち、システムの計算処理時の消費電力波形を観察する電力解析攻撃に関しては、統計解析手法等の知識のある者がオシロスコープ、パーソナルコンピュータなどの比較的入手し易い電子機器を揃えることで実行を可能としてしまう問題があった。また、この電力解析攻撃を実行した場合、その実行痕跡が残らないことも問題であった。以上から電力解析攻撃への耐性を検証するため、暗号処理システムの安全性評価方法を確立しなければならないという問題があった。
【００１１】
本発明は、前記問題を解決するため、検証する波形のノイズによる変形の影響を受けない、サンプリング間隔以下でも波形取得の際の位置ずれの量を検出できる、位相限定相関法に基づく高精度な波形位置合わせ方法を提案することを目的とする。この方法は、離散フーリエ変換した波形より得られる位相成分に着目した位相限定相関関数が２つの信号波形間でデルタ関数に近いピーク特性を示すことを利用し、ノイズを含む波形に対して高精度に位置ずれ量を推定する方法である。さらには、検証する波形に合わせて相関ピークモデルを当てはめることで、サンプリング分解能を越える精度で位置ずれの量を推定する方法である。
【課題を解決するための手段】
【００１２】
前記目的を達成するため、本発明の請求項１に記載の波形の位置合わせ方法は、
波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、
前記位相合わせ手段は、
外部からの、入力波形信号を離散フーリエ変換し算出するステップ及び参照波形信号を離散フーリエ変換し算出するステップとを有し、該入力波形信号及び参照波形信号をそれぞれ離散フーリエ変換した際に分離される位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し相関関数を検出することを特徴とする。
【００１３】
本発明の請求項２に記載の波形の位置合わせ方法は、
波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、
前記位相合わせ手段は、
外部からの、入力波形信号を離散フーリエ変換し算出するステップ、予め設定された参照波形信号を離散フーリエ変換し算出するステップとを有し、該入力波形信号及び参照波形信号をそれぞれ離散フーリエ変換した際に分離される位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し相関関数を検出し、前記相関関数に相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする。
【００１４】
本発明の請求項３に記載の波形の位置合わせ方法は、請求項２に記載の波形の位置合わせ方法であって、
前記位相合わせ手段は、フーリエ変換の性質から理論的に算出される相関ピークの形状を近似する相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする。
【００１５】
本発明の請求項４に記載の波形の位置合わせ方法は、請求項２に記載の波形の位置合わせ方法であって、
前記位相合わせ手段は、２次関数を用いた相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする。
【００１６】
本発明の請求項５に記載の波形の位置合わせ方法は、請求項２に記載の波形の位置合わせ方法であって、
前記位相合わせ手段は、ガウス関数を用いた相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定することを特徴とする。
【００１７】
本発明の請求項６に記載の波形の位置合わせ方法は、
波形信号の位相による位置合わせ処理を行う位相合わせ手段を有する波形の位置合わせ方法において、
前記位相合わせ手段は、
外部からの、入力波形信号を離散フーリエ変換し算出するステップ及び参照波形信号を離散フーリエ変換し算出するステップあるいは外部からの、入力波形信号を離散フーリエ変換し算出するステップ及び予め設定された参照波形信号を離散フーリエ変換し算出するステップ、いづれかの入力波形信号及び参照波形信号をそれぞれ離散フーリエ変換した際に分離される位相信号を抽出し、得られた位相差信号を逆離散フーリエ変換し相関関数を検出し、前記相関関数に相関ピークモデルを当てはめることで波形信号間の位置ずれの量を推定し、前記位置ずれ量を用いて波形リサンプリングすることで該入力波形の位置ずれ補正した波形を算出することを特徴とする。
【発明の効果】
【００１８】
以上述べたように本発明によれば、波形信号と参照波形信号を離散フーリエ変換した波形より得られる位相差信号を逆離散フーリエ変換することから、サンプリング分解能を越える精度で波形信号間の位置ずれ量を補正することができる。また、電力解析や電磁波解析の統計処理に悪影響を与える波形を排除することが可能である。更には、暗号装置のソフトウェアに対する電力・電磁波解析に適用すれば、大幅な測定精度の向上を図ることができ、従来では解析できなかった低いサンプリングレート、あるいは同等のサンプリングレートではより少ない波形サンプル数でも測定ができる。
【００１９】
位相限定相関法による波形位置合わせ手法においては、さまざまな暗号アルゴリズムに対して、ソフトウェアだけでなくハードウェアによる実装であっても解析処理に高い効果がある。また、電力解析だけでなく電磁波解析等にも利用できる。さらにはトリガ信号やシステムクロックの検出が難しい場合であっても、高い精度での解析処理が期待できる。以上の特徴から、位相限定相関法による波形位置合わせ手法は、汎用性に優れた解析ソフトウェアの構築を可能とする。
【発明を実施するための最良の形態】
【００２０】
本発明の実施の形態に係る位相限定相関（ＰＯＣ：Phase-Only Correlation）法を用いた波形位置合わせ方法について図面を参照しながら説明する。本明細書における位置合わせ方法とは、主に統計処理による波形解析のために位相限定相関法を用いた波形の位置合わせ方法を実現したソフトウェアを指すものである。
【００２１】
［波形の位置合わせ方法］
本発明の実施の形態は、波形位置合わせにおいては外部からの入力波形信号と参照する波形信号を登録し、その信号を離散フーリエ変換し、振幅信号と位相信号に分離したうえで形状情報(変化成分)が含まれる位相信号のみを抽出する。次に得られた位相差信号を逆離散フーリエ変換し位相限定相関関数を算出するものである。なお、参照波形信号は、入力波形信号と同様に外部から与えてもよい。また、位相差信号は合成位相スペクトルとも呼ばれる。
【００２２】
本発明の実施の形態に係る、波形の位置合わせ方法について図面を参照して説明する。
図１は波形の位置合わせの方法の構成を示す一実施の形態のブロック図である。尚、(ａ)の波形位置合わせ装置１は入力波形信号、参照波形信号を外部より取り込み、演算による位置ずれ量の推定処理を行った上で、入力波形信号の位置ずれ補正を行った波形信号を検出するものである。また、(ｂ)の波形位置合わせ装置１は入力波形信号を外部より取り込み、且つ参照波形信号は波形位置合わせ装置１に内蔵される記憶部(特に図示せず)より任意に選択でき、演算による位置ずれ量の推定処理を行った上で、入力波形信号の位置ずれ補正を行った波形信号を検出するものである。
【００２３】
図１(ａ)について、詳細に説明する。波形位置合わせ装置１は位相限定相関関数算出手段２、波形位置解析手段３、波形位置補正手段４から構成される。
【００２４】
次に図１(ａ)での処理内容について説明する。波形位置合わせ装置１は外部より入力波形信号、参照波形信号を取り込み、位相限定相関関数算出手段２へ入力される。位相限定相関関数算出手段２では入力波形信号及び参照波形信号の取得領域を指定したうえで、それぞれ離散フーリエ変換(ＤＦＴ：Discrete Fourier Transform)して位相信号のみを抽出し、得られた位相差信号を逆離散フーリエ変換(ＩＤＦＴ：Inverse Discrete Fourier Transform)することで位相限定相関関数を算出する。波形位置解析手段３では、位相限定相関関数算出手段２により得られる位相限定相関関数に相関ピークモデルを当てはめることで、波形信号間の位置ずれの量（相関ピーク座標とも呼ばれる）を推定する。波形位置補正手段４では、波形位置解析手段３で得られる位置ずれ量を用いて、入力波形信号の位置ずれ補正を行った波形信号を算出することになる。
【００２５】
図１(ｂ)での処理内容について説明する。図１(ａ)とは波形位置合わせ装置１に組み込まれる参照波形信号情報が異なる。例えば、図１(ａ)の参照波形信号は外部から得られる波形信号であり、２つの波形信号を比較する目的であるから、参照波形信号は入力波形信号に近い値が入力されるものである。これに対して、図１(ｂ)の参照信号選択手段５には、予め基準となる波形信号が登録されており、入力波形信号に応じて任意の基準波形信号を参照波形信号として利用するものである。
【００２６】
次に、図２を参照して位相限定相関法を用いた波形位置合わせを行うための処理を示すフローを説明する。
【００２７】
まず、図２(ａ)のステップＳ１では、外部からの入力波形信号及び参照波形信号を検出する。ステップＳ２では入力波形信号及び参照波形信号の取得領域を指定する。ステップＳ３ではＳ２からの入力波形信号、基準となる参照波形信号の離散フーリエ変換処理を行う。ステップＳ４ではＳ３で処理された入力波形信号の位相信号と参照波形信号の位相信号のみを抽出する。次にステップＳ５では、入力波形信号と参照波形信号の位相差信号の逆離散フーリエ変換として、位相限定相関関数を算出する。このとき、入力される波形信号に応じて、（１）位相差信号に対してスペクトル重み付け関数を適用した位相限定相関関数（スペクトル重み付け型位相限定相関関数）、（２）位相差信号に対して逆離散フーリエ変換する周波数帯域を制限した位相限定相関関数（帯域制限位相限定相関関数）、もしくは（３）位相差信号に対してスペクトル重み付け関数を適用した帯域制限位相限定相関関数を算出してもよい。ステップＳ６では算出された相関関数のピークから波形信号間の位置ずれの量を推定する。ステップＳ７では、Ｓ６で得られた位置ずれの量を用いた入力波形信号の波形リサンプリングにより、位置ずれ補正した波形信号を算出するものである。
【００２８】
また、図２(ｂ)のステップＳ１では、外部からの入力波形信号を検出する。ステップＳ２では入力波形信号の取得領域を指定する。ステップＳ３では入力波形信号を離散フーリエ変換する。ステップＳ４ではＳ３で算出された入力波形信号の位相信号のみを抽出する。次にステップＳ５では基準となる参照波形信号を検出する。ステップＳ６では参照波形信号を離散フーリエ変換する。ステップＳ７ではＳ６で算出された参照波形信号の位相信号のみを抽出する。ステップＳ８では入力波形信号と参照波形信号の位相差信号の逆離散フーリエ変換として、位相限定相関関数を算出する。このとき、入力される波形信号に応じて、（１）位相差信号に対してスペクトル重み付け関数を適用した位相限定相関関数（スペクトル重み付け型位相限定相関関数）、（２）位相差信号に対して逆離散フーリエ変換する周波数帯域を制限した位相限定相関関数（帯域制限位相限定相関関数）、もしくは（３）位相差信号に対してスペクトル重み付け関数を適用した帯域制限位相限定相関関数を算出してもよい。ステップＳ９では算出された相関関数のピークから波形信号間の位置ずれの量を推定する。ステップＳ１０では、Ｓ９で得られた位置ずれの量を用いた入力波形信号の波形リサンプリングにより、位置ずれ補正した波形信号を算出するものである。
【００２９】
［位相限定相関法を用いた位置ずれ量検出の手法］
次に本発明の実施の形態に係る位相限定相関法を用いた位置ずれ量検出の手法について説明する。
Ｎ点の２つの信号波形をｆ（ｎ）、ｇ（ｎ）とする。
但し、定式化の便宜上、離散時間のインデックスをｎ＝−Ｍ，Λ，Ｍとし、信号波形の長さをＮ＝２Ｍ＋１とする。これらの信号波形の離散フーリエ変換をそれぞれＦ（ｋ）及びＧ（ｋ）として次の数式１で与える。
【００３０】
【数１】
2007240170.tif 000003
【数１ａ】
2007240170.tif 000004
【数１ｂ】
2007240170.tif 000005
【数１ｃ】
2007240170.tif 000006
【００３１】
但し、数式１ａである。ここで、ＡＦ（ｋ）及びＡＧ（ｋ）は、それぞれ信号波形ｆ（ｎ）及びｇ（ｎ）の振幅成分、数式１ｂ及び数式１ｃはそれぞれの信号の位相成分である。
一般性を失うことなく離散周波数のインデックスをｋ＝−Ｍ，Λ，Ｍとすることができる。このとき、合成位相スペクトルＲＦＧ（ｋ）は、
【００３２】
【数２】
2007240170.tif 000007
【数２ａ】
2007240170.tif 000008
【００３３】
ここで、数式２ａはＧ（ｋ）の複素共役であり、θＦＧ＝θＦ（ｋ）−θＧ（ｋ）である。このように、ＲＦＧ（ｋ）はＦ（ｋ）及びＧ（ｋ）の位相差を表す。ｆ（ｎ）とｇ（ｎ）の位相限定相関関数（以下ＰＯＣ関数と呼ぶ）ｒｆｇ（ｎ）はＲＦＧ（ｋ）のＩＤＦＴとして、次のように表される。
【００３４】
【数３】
2007240170.tif 000009
【００３５】
次に連続時間で定義された信号波形ｆｃ（ｔ）を考える。ここで、ｔは実数である。いま、δをｔに関する微小移動量を表す実数とすると、連続時間でｆｃ（ｔ）をδだけ微小移動した波形はｆｃ（ｔ−δ）と表現できる。これらの連続時間信号ｆｃ（ｔ）及びｆｃ（ｔ−δ）を標本化間隔Ｔで標本化した離散時間信号をそれぞれｆ（ｎ）及びｇ（ｎ）とし、次の数式４で定義する。
【００３６】
【数４】
2007240170.tif 000010
【００３７】
但し、ｎ＝−Ｍ，Λ，Ｍとする。以下では単純化するため、Ｔ＝１とする。
このとき、離散時間信号ｆ（ｎ）及びｇ（ｎ）に関するＰＯＣ関数を用いて、連続時間での微小移動量δを推定する問題を考える。
但し、δは、離散時間においてサンプリング間隔以下の移動量に対応するものと考える。ここで、ｆ（ｎ）及びｇ（ｎ）のＤＦＴＦ（ｋ）及びＧ（ｋ）の間には次の近似が成り立つ。
【００３８】
【数５】
2007240170.tif 000011
【００３９】
数式５が近似であるのは、連続時間信号と離散時間信号に対するフーリエ変換の性質の違いに起因する。このとき、ｆ（ｎ）及びｇ（ｎ）の合成位相スペクトルＲＦＧ（ｋ）及びＰＯＣ関数ｒｆｇ（ｎ）は、次のように表せる。
【００４０】
【数６】
2007240170.tif 000012
【００４１】
ここで、α＝＜１である。数式６は、波形間に微小移動量δがある場合のＰＯＣ関数の一般形を表している。相関ピークの座標は波形間の位置ずれを表し、相関ピークの高さαは波形間の類似度の指標となる。波形の変化によってαの値は変化する。このように微小移動した２つの波形のＰＯＣ関数は1サンプルの幅の急峻なピークをもつ。以上から、ＰＯＣ関数の相関ピークを検出することにより、波形間の移動量及び類似度を評価することができる。
【００４２】
前記のＰＯＣ関数の算出においては、入力される波形信号に応じて、（１）スペクトル重み付け型位相限定相関関数、（２）帯域制限位相限定相関関数、もしくは（３）スペクトル重み付け関数を適用した帯域制限位相限定相関関数を算出してもよい。ここで、帯域制限位相限定相関関数を用いた場合、波形間の位置ずれ量は、相関ピーク座標の値に帯域幅に対する波形信号の長さの比を乗じたものとなる。
【００４３】
［移動量（位置ずれ量）推定の高精度化手法］
次に図３を用いて、高精度な移動量（位置ずれ量）推定の手法について述べる。図３の黒丸（黒い点の座標）はＰＯＣ関数を計算して得られた相関ピーク近傍のデータを示す。２つの波形ｆ（ｎ）及びｇ（ｎ）のＰＯＣ関数を計算することによって得られるｒｆｇ（ｎ）はｎ＝−Ｍ，Λ，Ｍの離散点のみのデータである。このとき、相関ピークモデルを実データに当てはめる（フィッティングする）ことで、波形のサンプル間に存在する真のピーク座標を推定することができる。図３の例では、ピークモデルの当てはめによりδ＝０．５付近が波形信号間の位置ずれの量として検出される。相関ピークモデルとしては、式(９)（数式６参照）のようなフーリエ変換の性質から理論的に算出される相関ピークの形状を近似するモデルを利用する。スペクトル重み付け型位相限定相関関数を用いた場合、スペクトル重み付け関数と式(８)（数式６参照）で表される合成位相スペクトルＲＦＧ（ｋ）の積をＩＤＦＴして算出される相関ピークモデルを利用する。帯域制限位相限定相関関数を用いた場合、通常のＰＯＣ関数と同様の相関ピークモデルを利用する。スペクトル重み付け関数を適用した帯域制限位相限定相関関数を用いた場合、スペクトル重み付け関数と式(８)（数式６参照）で表される合成位相スペクトルＲＦＧ（ｋ）の積を帯域制限してＩＤＦＴすることにより算出される相関ピークモデルを利用する。ここで、相関ピークモデルには、式(９)（数式６参照）及び前記方法で算出される相関ピークモデル以外のモデル(２次関数やガウス関数など)を利用することも可能である。
【００４４】
上述した方法により検出された位置ずれ量を用いて波形の位置合わせを行う。波形位置合わせとは、２つの信号波形ｆ（ｎ）及びｇ（ｎ）の位置ずれ量をδとすると、ｇ（ｎ）をδ分だけずらしたｇ’（ｎ）を求めることである。例えば、周波数領域でｇ（ｎ）の位相を回転させることによりｇ’（ｎ）を以下のように算出できる。
【００４５】
位置ずれ量δのとき、ｇ’（ｎ）のＤＦＴＧ’（ｋ）は次のように近似できる。
【００４６】
【数７】
2007240170.tif 000013
【００４７】
このとき、ｇ’（ｎ）はＧ’（ｋ）のＩＤＦＴとして次の数式８で求まる。
【００４８】
【数８】
2007240170.tif 000014
【００４９】
波形位置合わせの方法には、前記以外にも、双三次補間法などの各種補間方法を利用できる。また、測定機器の内部で検出した位置ずれ量δを補正することも可能と考えられる。提案する手法による波形位置合わせの例を図４に示す。図４(ａ)は、マイクロプロセッサ(８ＭＨｚ動作)の消費電力を１ＧＳａ/ｓで取得した２つの波形である。ここで、取得した波形間には位置ずれが確認できる。位相限定相関法及びその高精度化手法を適用した結果、２つの波形の位置ずれ量はδ＝１．５５５５であった。図４(ｂ)に位置合わせ後の波形を示す。このように、同形ではない信号波形間の位置ずれに対しサンプリング分解能を越える精度で補正することができる。
【００５０】
図５には本発明における位置ずれ量検出と波形リサンプリングによる位置ずれ補正した信号波形を示す線図を示す。
【００５１】
［差分電力解析への応用］
次に本発明で提案する波形位置合わせ手法を用いた差分電力解析（ＤＰＡ：Differential Power Analysis）について述べる。暗号装置への電力解析とは、暗号処理によって変化する消費電力に着目して鍵などの秘密情報が読み取れるかどうかを判定することである。ＤＰＡとは、暗号処理による消費電力の変化が小さく、測定誤差やノイズなどの影響で測定波形から秘密情報を直接見分けることが困難な場合に、数千から数万パターンの消費電力波形を統計処理して、わずかな秘密情報を増幅させる解析法である。鍵のビットパターンの違いが公開鍵暗号ほど電力波形に現れない共通鍵暗号の実装に対して、単純電力解析(ＳＰＡ：Simple Power Analysis)よりも有効である。その解析能力の高さから暗号装置の安全性を評価する上で欠かすことのできない解析の一つであり、高精度な実施技術の確立が課題となっている。
【００５２】
DＰＡでは、秘密情報の抽出のため電力波形データへの統計的な処理を必要とする。そのため、入力データの異なる複数の電力波形を同じ処理が行われる正確なタイミングで取得しなければならない。しかし、実際の測定で完全に同一のタイミングを実現することは難しい。まず、暗号処理に同期した都合のよいトリガ信号が得られるとは限らない。ＰＬＬを実装したＬＳＩでは内部クロックが外部クロック等の制御信号と必ずしも同期している保証はない。また、トリガを得られたとしても、そのジッタ等により取り込み誤差が生じる。結果として、取得された電力波形データには常に位置ずれが含まれていることになる。その位置ずれは、測定機器のサンプリング間隔以下と微小な場合でもしばしば解析結果に大きな影響を与える。この電力波形間の位置ずれの問題を解決することができれば、波形取得のタイミング精度に影響されない強力なＤＰＡを実施できると考えられる。
【００５３】
以下では、暗号評価プラットホームであるＩＮＳＴＡＣ−８のプロセッサ上にソフトウェア実装したＤＥＳ（共通鍵暗号の一つ）に対して、本発明の波形位置合わせ手法を用いたＤＰＡを行った実験について記載する。本発明を用いたＤＰＡの概要を図６に示す。尚、従来のＤＰＡは、クロックやトリガ信号等によって計測時に波形の位置を合わせている。本発明を用いたＤＰＡでは、同様に複数の波形を取得後、取得した波形に対してさらに位相限定相関法で極めて高い精度の位置合わせを行うものとした。図６における位置合わせでは、まず、取得した波形の中から位置合わせの基準となる波形を参照波形として登録する。
【００５４】
次に、上述した位相限定相関関数及び高精度な位置ずれ量推定手法を用いて、参照波形と他の波形との位置ずれ量を検出する。得られた位置ずれ量を用いて波形リサンプリングを行い、位置ずれ補正した波形を算出する。この位置ずれを補正した波形の様子を図５に示す。補正後の処理は従来のＤＰＡと同様である。まず、秘密鍵を予測し、その秘密鍵で得られる中間データを計算する。次に、その中間データのうち１ビットが０か１かで消費電力波形を２グループに分けて平均を取る。秘密鍵の予想が合っており、数千から数万パターンの消費電力波形を正しく振り分けることができるのであれば、それを集めて平均を取った波形にピークが現れる。予想が誤っているとランダムに振り分けられた波形の平均を取ることになるため、そのようなピークは現れない。 
【００５５】
図７(ａ)、（ｂ）にそれぞれ従来のＤＰＡ結果と本発明を用いたＤＰＡ結果を示す。ここでは、２００ＭＳａ／ｓで取得した波形の場合を示す。図７(ａ)のように、従来のＤＰＡでは誤った鍵の推定時にもっとも高いピークが現れた。一方、本発明を用いたＤＰＡでは、正しい鍵推定の場合にのみ高いピークが現れた。ここで、図７(ａ)の方が（ｂ）と比べてピークとノイズの比が低下していることにも注意されたい。このように、従来では解析できない低いサンプリングレートでも、本発明を用いることにより、精度のよい解析が可能となることが確認された。
【００５６】
サンプリングレートが１００ＭＳａ／ｓ、２００ＭＳａ／ｓ、４００ＭＳａ／ｓ及び１ＧＳａ／ｓにおけるＤＰＡの失敗率と波形サンプル数の関係を図８(ａ)（ｂ）に示す。(ａ)は従来のＤＰＡ、（ｂ）は本発明を用いたＤＰＡの結果である。縦軸のエラービット数は、推定に用いた出力全３２ビットのうち鍵推定に誤ったビットの数を示す。３２の出力すべてで推定を誤った場合のエラービット数は３２となる。エラービット数は０に下がらなくても、部分的に鍵候補を絞ることができるため、電力解析に対する耐性を示す指標となる。横軸は統計処理に用いた波形サンプル数である。図８から、本発明を用いたＤＰＡがすべてのサンプリングレートで従来のＤＰＡよりも高い精度で鍵の推定に成功している。２００ＭＳａ／ｓとサンプリング分解能が低く従来法ではほとんど鍵が求まらないような場合でも、本発明によって５０％以上の鍵推定に成功した。また、４００ＭＳａ／ｓの場合、鍵推定の精度を同じとすると、本発明を用いたＤＰＡは従来の２／３程度と少ないサンプル数しか必要としないことも図８から判る。
【産業上の利用可能性】
【００５７】
本発明は、統計処理のような複数の波形を処理する波形解析であれば汎用的に活用することができるため、暗号処理システムの安全性評価のための利用のみならず、ディジタル測定機器に直接組み込むことや、ディジタル測定機器とオフラインで接続する波形データの加工処理ソフトに組み込むことで各種の統計処理による波形解析を高精度に行うことが可能となる。 
【図面の簡単な説明】
【００５８】
【図１】本発明の一実施の形態に係る位相限定相関法を用いた波形位置合わせを行うためのブロック図である。
【図２ａ】本発明の一実施の形態に係る位相限定相関法を用いた波形位置合わせを行うための処理を示すフローチャート図である。
【図２ｂ】本発明の一実施の形態に係る位相限定相関法を用いた波形位置合わせを行うための処理を示すフローチャート図である。
【図３】ＰＯＣ関数を計算して得られた相関ピーク近傍のデータと相関ピークモデルを当てはめた例を示す図である。
【図４】本発明における波形の位置合わせの例を示す図である。
【図５】本発明における位置ずれ量検出と波形のリサンプリングによる位置ずれ補正を示す線図である。
【図６】本発明における波形の位置合わせを用いたＤＰＡのブロック図である。
【図７】従来のＤＰＡと本発明における波形の位置合わせを用いたＤＰＡの結果を示す図である。
【図８】従来のＤＰＡと本発明における波形の位置合わせを用いたＤＰＡの失敗率と波形サンプル数の関係を示す図である。
【符号の説明】
【００５９】
１      波形位置合わせ装置
２      位相限定相関関数算出手段
３      波形位置解析手段
４      波形位置補正手段
５      参照信号選択手段
【図１】
2007240170.tif 000015
【図２ａ】
2007240170.tif 000016
【図２ｂ】
2007240170.tif 000017
【図３】
2007240170.tif 000018
【図４】
2007240170.tif 000019
【図５】
2007240170.tif 000020
【図６】
2007240170.tif 000021
【図７】
2007240170.tif 000022
【図８】
2007240170.tif 000023