(57)【特許請求の範囲】
【請求項１】
複数のチャネルｋ（ｋ＝１、２、…、Ｋ；Ｋは２以上の予め定められた正整数）の信号について、各チャネルの信号が全てのチャネルに共通の信号（以下、「共通信号」という）と当該共通信号以外の信号（以下、「残信号」という）との和で表されるとし、
ｋ＝１、２、…、Ｋについて、チャネルごとに元の信号が圧縮された圧縮信号ｙｋは、元の信号ｓｋの基底行列Ψによる表現ｘｋに対して観測行列Φｋを用いてｙｋ＝Ａｋｘｋ＝Φｋｓｋ、ｓｋ＝Ψｘｋ、Ａｋ＝ΦｋΨなる関係があるとして、
チャネルごとに元の信号が圧縮された圧縮信号ｙｋをチャネルごとに独立に元の信号（以下、「仮の復元信号」という）を復元する仮伸長部と、
各チャネルの上記仮の復元信号ｓ~ｋから、上記共通信号に対する相対的な振幅の
大きさと遅延量を表す振幅係数ａｋｃと群遅延ｂｋｃを求める差異推定部と、
上記圧縮信号ｙｋと、上記振幅係数ａｋｃと、上記群遅延ｂｋｃと（ｋ＝１、２、…、Ｋ）を用いて、Λを回転行列
【数１５】
0005785903.tif 000016
、ｊを虚数単位、Ｎを元の信号のベクトルの要素数として、上記共通信号の周波数領域表現
【数１６】
0005785903.tif 000017
を求める共通信号推定部と、
上記圧縮信号ｙｋと、上記振幅係数ａｋｃと、上記群遅延ｂｋｃと、上記共通信号の周波数領域表現ｘ'ｃを用いて、上記圧縮信号ｙｋと上記共通信号の周波数領域表現ｘ'ｃに基づく値との差分ｙｄｋを得て、さらに当該差分ｙｄｋについて圧縮センシングのアルゴリズムを適用して、時間領域表現であるｓ'ｄｋを得て、当該ｓ'ｄｋと上記共通信号の時間領域表現との和によって、元の信号の推定値ｓ'ｋを得る残信号推定部と
を含む圧縮信号復元装置。
【請求項２】
請求項１に記載の圧縮信号復元装置において、
さらに、上記共通信号の周波数領域表現ｘ'ｃを用いて、行列Ａｋから残信号を表現する基底行列または上記共通信号を表現する基底行列を推定する残信号基底推定部を含み、
上記共通信号推定部は、上記残信号基底推定部によって得られた基底行列に基づいて上記共通信号の周波数領域表現ｘ'ｃを更新する処理を行う
ことを特徴とする圧縮信号復元装置。
【請求項３】
請求項２に記載の圧縮信号復元装置において、
上記残信号基底推定部によって得られた基底行列に基づいて残信号が除去された圧縮信号を上記仮伸長部の入力とする
ことを特徴とする圧縮信号復元装置。
【請求項４】
複数のチャネルｋ（ｋ＝１、２、…、Ｋ；Ｋは２以上の予め定められた正整数）の信号について、各チャネルの信号が全てのチャネルに共通の信号（以下、「共通信号」という）と当該共通信号以外の信号（以下、「残信号」という）との和で表されるとし、
ｋ＝１、２、…、Ｋについて、チャネルごとに元の信号が圧縮された圧縮信号ｙｋは、元の信号ｓｋの基底行列Ψによる表現ｘｋに対して観測行列Φｋを用いてｙｋ＝Ａｋｘｋ＝Φｋｓｋ、ｓｋ＝Ψｘｋ、Ａｋ＝ΦｋΨなる関係があるとして、
仮伸長部が、チャネルごとに元の信号が圧縮された圧縮信号ｙｋをチャネルごとに独立に元の信号（以下、「仮の復元信号」という）を復元する仮伸長ステップと、
差異推定部が、各チャネルの上記仮の復元信号ｓ~ｋから、上記共通信号に対する相対的な振幅の大きさと遅延量を表す振幅係数ａｋｃと群遅延ｂｋｃを求める差異推定ステップと、
共通信号推定部が、上記圧縮信号ｙｋと、上記振幅係数ａｋｃと、上記群遅延ｂｋｃと（ｋ＝１、２、…、Ｋ）を用いて、Λを回転行列
【数１７】
0005785903.tif 000018
、ｊを虚数単位、Ｎを元の信号のベクトルの要素数として、上記共通信号の周波数領域表現
【数１８】
0005785903.tif 000019
を求める共通信号推定ステップと、
残信号推定部が、上記圧縮信号ｙｋと、上記振幅係数ａｋｃと、上記群遅延ｂｋｃと、上記共通信号の周波数領域表現ｘ'ｃを用いて、上記圧縮信号ｙｋと上記共通信号の周波数領域表現ｘ'ｃに基づく値との差分ｙｄｋを得て、さらに当該差分ｙｄｋについて圧縮センシングのアルゴリズムを適用して、時間領域表現であるｓ'ｄｋを得て、当該ｓ'ｄｋと上記共通信号の時間領域表現との和によって、元の信号の推定値ｓ'ｋを得る残信号推定ステップと
を含む圧縮信号復元方法。
【請求項５】
請求項１から請求項３のいずれかに記載の圧縮信号復元装置としてコンピュータを機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数のチャネルの信号が圧縮された圧縮信号の復元技術に関する。
【背景技術】
【０００２】
近年、多数のセンサーを用いて環境（温度、照度など）や人間の行動などを逐一センシングし、その情報を活用しようという動きが高まっている。多数のセンサーをネットワークによって結び、お互いに通信すること或いは特定の中央受信装置へ情報を送ることを可能にしたシステムはセンサーネットワークと呼ばれる。
【０００３】
このようなシステムでは、それぞれのセンサーは離れた位置に配置されており、その運用を簡便にするため無線通信を用いていることが多い。また、センサーの電源は電池を用い外部からの電源入力を得ていない場合も想定される。そのため、センサーが送信すべき情報をできうる限り削減すること、およびセンサー自身における演算量を低減することが重要であり多くの提案がなされている（たとえば非特許文献１）。
【０００４】
一方で、センサーネットワークは多くのセンサーを有し、かつ、用途によってはそれらが非常に近接している場合がある。このような構成では、各センサーが計測した結果は、或る観点からすれば非常に似通っていることが期待される（図７参照。図７は、３ chの生体信号（筋電信号）を重ねて表示したグラフであり、各チャネルの信号が非常に類似していることがわかる）。このため、例えば、分散情報源符号化や階層的に情報の相関を利用して送信すべき情報量をなるべく減らそうという試みがなされている。
【０００５】
このような試みのなかで近年注目を受けているのが圧縮センシングという手法である（例えば非特許文献２）。単一センサーでの圧縮センシングは測定の対象となる標本数Ｎの信号ｓに比して少ない数Ｍの情報ｙ＝Φｓ（ΦはＭ×Ｎの行列）のみをもって、元の信号ｓが或る基底行列Ψにおいての信号ｓの表現ｘ＝Ψｓの非0要素が非常に少ない（以下、「疎である」と表現する）という仮定の下に復元を行うことで元の信号ｓの厳密に等しいｓ’を得ようとするものである。また、この拡張タイプとして測定結果にノイズが含まれているケースについては同様なアルゴリズムで近似的な解ｓ”が得られることが分かっている。
【０００６】
圧縮センシング技術において、センサーネットワークのような多数のセンサーが存在する場合に応用する提案はすでにいくつか存在している。そのうちの一つが、各センサーによって得られた測定信号ｓｋ（ｋはセンサーの識別インデックス）を全てのセンサーに共通した信号ｓｃとそれ以外の信号ｓｄｋの和で表現できると仮定したものである。このモデルにもいくつかの派生型があるが、例えば共通信号の基底行列Ψでの表現ｘｃ＝Ψｓｃが疎ではなく、一方でｘｄｋ＝Ψｓｋが全てのｋに対して疎であるようなモデルがある。観測の対象となる信号がこのようなモデルで記述できる場合は、単独のセンサーでは信号が疎ではないために単一圧縮センシングの枠組みでは長さＭの情報ｙｋ＝Φｋｓｋから長さＮの情報ｓｋを正確に復元することはできない。しかしながら、多数のセンサーがあることで互いに情報を補い合い送信すべき情報を削減することができる場合がある。
【０００７】
各センサーはＣＰＵなどで演算を行うことが可能であり、またメモリなどにデータを保持することも可能である。また、通常の無線送受信機あるいは可視光通信の送受信機を用い双方向ないし一方向の通信を行うことができる。
【０００８】
以下に、通常の圧縮センシングについて説明する。通常の圧縮センシングでは、下記の条件に基づいて元の信号の推定値を得る。
［条件］ｙ＝Ａｘ＝Φｓを満たすような量||ｘ||１を最小化する      …（１）
【０００９】
ただし、||・||１はベクトル・の絶対値和（Ｌ１ノルム）を表し、ｙを観測ベクトル、Ａ、Φ、Ψはそれぞれ圧縮行列、観測行列、基底行列とする。またｓ＝Ψｘであり、推定された推定信号ｘから元の信号ｓを計算することができる。
【００１０】
また、上記条件に代えて、下記のような条件の問題に緩和することもできる。
||ｙ―Ａｘ||２＜εの下で量||ｘ||１を最小化する      …（２）
【００１１】
ただし、εはあらかじめ決められた微小なスカラー量であり、||・||２はベクトル・の二乗和である。
【００１２】
これ以外にも、条件によっては上記と同様の効果を持ついくつかのアルゴリズムが発見されており（例えば非特許参考文献３）、本明細書ではこれらを総称して「通常の圧縮センシングのアルゴリズム」と呼ぶことにする。
【先行技術文献】
【非特許文献】
【００１３】
【非特許文献１】D. Baron，M.B. Wakin，M.F. Duarte，S. Sarvotham，and R.G. Baraniuk. “Distributed compressed sensing,” Rice Univ. Technical report，2005.
【非特許文献２】Emmanuel J. Candes and Michael B. Wakin, "An Introduction ToCompressive Sampling" IEEE SIGNAL PROCESSING MAGAZINE, MARCH 2008.
【非特許文献３】J.A. Tropp, S.J. Wright, "Computational Methods for Sparse Solution of Linear Inverse Problems," Proceedings of the IEEE, vol.98, no.6, pp.948-958, June 2010.
【発明の概要】
【発明が解決しようとする課題】
【００１４】
各センサーが測定した信号ｓｋ（ｋはセンサーの識別インデックス）を全てのセンサーに共通した信号ｓｃとそれ以外の信号ｓｄｋの和で表現できると仮定したモデルにもいくつかの派生型があるが、例えば、共通信号の基底行列Ψでの表現ｘｃ＝Ψｓｃが疎ではなく、一方でｘｄｋ＝Ψｓｋが全てのインデックスｋに対して疎であるような信号モデルがある。このような信号モデルでも、多数のセンサーがある場合には、この信号モデルに基づいた復元を行うことによって、各々のセンサーについては長さＭ＜Ｎの情報ｙｋを用いて元の情報ｓｋを得ることが可能である場合がある。
【００１５】
しかしながら、上記の信号モデルで観測信号のモデル化を行った場合、共通信号ｓｃが全てのセンサーの間で厳密には共有されていないと、却って復元誤差を大きくしてしまう場合もありうる。
【００１６】
従来の圧縮センシングでは共通信号が厳密に一致していない場合には共通信号の基底行列Ψでの表現ｘｃ＝Ψｓｃも疎であるという仮定を導入して近似的に復元問題を解こうという試みが行われてきたが、得られた情報を十分活用しているとは言い難い。
【００１７】
本発明はこのような点に鑑みてなされたものであり、各チャネル間で共通信号が厳密に一致していな場合においても、複数のチャネルの信号が圧縮された圧縮信号の復元技術を提供することを目的とする。
【課題を解決するための手段】
【００１８】
本発明の圧縮信号復元技術は、共通信号ｓｃとともに振幅係数ａｋｃ及び群遅延βｋｃを考慮した信号モデルに基づいて、これら２つの因子をチャネルごとに推定することによって、より正確な信号復元を行う。
具体的には、複数のチャネルｋ（ｋ＝１、２、…、Ｋ；Ｋは２以上の予め定められた正整数）の信号について、各チャネルの信号が全てのチャネルに共通の信号（共通信号）と当該共通信号以外の信号（残信号）との和で表されるとし、ｋ＝１、２、…、Ｋについて、チャネルごとに元の信号が圧縮された圧縮信号ｙｋは、元の信号ｓｋの基底行列Ψによる表現ｘｋに対して観測行列Φｋを用いてｙｋ＝Ａｋｘｋ＝Φｋｓｋ、ｓｋ＝Ψｘｋ、Ａｋ＝ΦｋΨなる関係があるとして、チャネルごとに元の信号が圧縮された圧縮信号ｙｋをチャネルごとに独立に元の信号（仮の復元信号）を復元し［仮伸長処理］、各チャネルの仮の復元信号ｓ~ｋから、共通信号に対する相対的な振幅の大きさと遅延量を表す振幅係数ａｋｃと群遅延ｂｋｃを求め［差異推定処理］、圧縮信号ｙｋと、振幅係数ａｋｃと、群遅延ｂｋｃと（ｋ＝１、２、…、Ｋ）を用いて、Λを回転行列
【数１】
0005785903.tif 000002
、ｊを虚数単位、Ｎを元の信号のベクトルの要素数として、共通信号の周波数領域表現
【数２】
0005785903.tif 000003
を求め［共通信号推定処理］、圧縮信号ｙｋと、振幅係数ａｋｃと、群遅延ｂｋｃと、共通信号の周波数領域表現ｘ'ｃを用いて、圧縮信号ｙｋと共通信号との差分ｙｄｋを得て、さらに当該差分ｙｄｋを復元してｓ'ｄｋを得て、当該ｓ'ｄｋと共通信号の時間領域表現との和によって、元の信号の推定値ｓ'ｋを得る［残信号推定処理］。
【発明の効果】
【００１９】
本発明に拠れば、本来必要な情報より少ない数の線形観測によってえられた情報を用いて、各チャネルの振幅差や遅延量を推定することで、各チャネルに含まれる共通信号が全く同一でなくチャネルごとに振幅や群遅延量が異なっている場合においても原信号を従来よりも精度良く復元することができる。
【図面の簡単な説明】
【００２０】
【図１】センサーネットワークを示す図。
【図２】従来の圧縮センシングの枠組を示す図。
【図３】本発明の実施形態１の機能構成を示す図。
【図４】本発明の実施形態２の機能構成を示す図。
【図５】本発明の実施形態４の機能構成を示す図。
【図６】比較実験例：（ａ）元の信号、（ｂ）ＴＥＣＣ法、（ｃ）本発明に拠る手法。
【図７】３ｃｈの生体信号の例。
【発明を実施するための形態】
【００２１】
＜記号の説明＞
Ｒ：実数の全体（集合）
Ｃ：複素数の全体（集合）
Ｎ∈Ｒ：観測対象のベクトルの大きさ
Ｍ∈Ｒ：１chあたりの観測数
ｓｋ∈ＣＮ：ｋ chの観測対象の元のベクトル表現
ｘｋ∈ＣＮ：ｋ chの観測対象の基底行列Ψでの表現
ｙｋ∈ＣＭ：ｋ chの観測結果のベクトル表現
ｓｃ∈ＣＮ：全チャネルに含まれる共通信号のベクトル表現
ｘｃ∈ＣＮ：共通信号の基底行列Ψでの表現
ｓｄｋ∈ＣＮ：ｋ chの共通信号を取り除いた信号（個別信号）のベクトル表現
ｘｄｋ∈ＣＮ：個別信号の基底行列Ψでの表現
Ψ∈ＣＮ×Ｎ：ｓｋ＝Ψｘとなるような基底行列
Φｋ∈ＣＭ×Ｎ：ｙｋ＝Φｓｋとなるような観測行列
Ａｋ∈ＣＭ×Ｎ：Ａｋ＝ΦｋΨで表せられる行列
ａｋｍ∈Ｒ：ｍ chを基準としたｋ chの振幅の大きさ
ｂｋｍ∈Ｒ：ｍ chを基準としたｋ chの位相遅延の大きさ
【００２２】
＜従来の信号モデルによる圧縮センシング＞
実施形態の説明に先立ち、図１と図２を用いて従来の信号モデルによる圧縮センシングについて述べる。
【００２３】
図１のようにＫ個（Ｋは２以上の正整数）のセンサーと１個の中央受信機９５０で構成されるシステムを考える。センサー９１０−１、センサー９１０−２、…、センサー９１０−Ｋはそれぞれ異なった信号を計測し、中央受信機９２０へ送信する。ここで各センサーが或る単位時間に取得する情報を複素ベクトルで表すこととし、このベクトルの長さはＮ（Ｎは正の整数）であるとする。すなわち、或る単位時間にＫ個のセンサーが取得した情報は一般に各々ｓｋ∈ＣＮ、ｋ＝１、２、…、Ｋと書けるとする。
【００２４】
図２のように、各センサー９１０は観測装置９１１と送信装置９１３で構成されている。観測装置９１１はふたつの部分、測定部９１１１と圧縮部９１１２で構成されている。測定部９１１１は、信号（例えば音声信号や生体信号などであり、特に限定はない）を観測し、単位時間の長さＮの信号ｓｋを出力する。圧縮部９１１２は、信号ｓｋを長さＭの圧縮された信号ｙｋ∈ＣＭに変換する。ただし、ＭはＮより小さな正の整数である。また、この変換は線形変換であるとする。すなわち、観測行列Φｋ∈ＣＭ×Ｎによって式（５）のように書けるとする。
ｙｋ＝Φｋｓｋ      …（５）
【００２５】
送信装置９１３は、符号化部９１３１によって圧縮信号ｙｋを通信路符号化して中央受信機９５０に送信する。
【００２６】
中央受信機９５０は、各センサーから、それぞれ通信路符号化されたＫ個の符号列を受信する。中央受信機９５０は、受信装置９５１と復元装置９５３で構成されている。受信装置９５１の復号部９５１１は、各センサーから受信した符号列をそれぞれ復号して圧縮信号ｙｋ、ｋ＝１、２、…、Ｋを出力する。
【００２７】
復元装置９５３は各行列Φｋ、ｋ＝１、２、…、Ｋ及びΨを既知とする。
また、各圧縮行列Ａｋは、式（６）のような行列の積で表すことができ、Φｋをセンサーｋの観測行列、Ψを基底行列と呼ぶこととする。
Ａｋ＝ΦｋΨ、Φｋ∈ＣＭ×Ｎ、Ψ∈ＣＮ×Ｎ、ｋ＝１、２、…、Ｋ      …（６）
【００２８】
復元装置９５３は、これら圧縮行列Ａｋについても既知であるとする。
【００２９】
従来技術では、ここで下記のような特別な場合を考えていた。すなわち、或る単位時間に得られた信号ｓｋ、ｋ＝１、２、…、Ｋは、式（７）のように全てのセンサーに共通した信号ｓｃとそれ以外の部分ｓｄｋで表すことができ、かつ、基底行列Ψにおけるｓｄｋの表現ｘｄｋ（式（８）参照）は非ゼロ要素が少ない疎なベクトルであるという信号モデルである。
ｓｋ＝ｓｃ＋ｓｄｋ      …（７）
ｓｄｋ＝Ψｘｄｋ        …（８）
【００３０】
このとき復元装置９５３を構成する共通信号推定部９５３１によって共通信号の基底行列Φにおける表現ｘ'ｃは式（９）のように推定される（例えば非特許参考文献１）。［・］＊は行列［・］のエルミート共役行列を表す。
【数３】
0005785903.tif 000004
【００３１】
復元装置９５３を構成する残信号推定部９５３２は、圧縮行列Ａｋと、式（９）に従って得られたｘ'ｃと、観測信号ｙｋとから式（１０）に従って残信号ｙ'ｄｋを計算し、ｙ'ｄｋとＡｋに対して通常の圧縮センシングのアルゴリズムを適用することでｘ'ｄｋを再構成する。
ｙ'ｄｋ＝ｙｋ−Ａｋｘ'ｃ      …（１０）
【００３２】
残信号推定部９５３２は、すべてのチャネル（ｋ＝１、２、…、Ｋ）について残信号ｙ'ｄｋを推定する。さらに、ｘ'ｃとｘ'ｄｋの和ｘ'ｋから元の基底による数ベクトル空間への変換、すなわち
ｓ'ｋ＝Ψｘ'ｋ      …（１１）
を計算し、元の信号の推定値ｓ'ｋを各チャネルについて得る。
【００３３】
上述した従来の信号再構成の方法は、下記のような信号モデルを前提としている。すなわち、各チャネルの測定信号ｓｋは、共通信号ｓｃと基底行列Φにおいて疎な信号ｓｄｋの和
ｓｋ＝ｓｃ＋ｓｄｋ      …（１２）
と表現される信号モデルである。しかしながら、音響信号や生体信号などの場合においては観測地点の違いなどから各チャネル間において共通信号ｓｃに相当する部分に振幅の違いや位相の違いがありうる。この信号モデルではこれらの影響を無視してしまっているため、実際の信号をよくモデル化しているとは言い難く、結果として上述した再構成方法では単チャネルの圧縮センシングによる再構成方法に比しても必ずしも性能が向上しない場合がある。このような課題に対する解決方法はいまだ提示されたことがない。
【００３４】
＜本発明の原理＞
本発明は、下記に述べるように共通信号に振幅の違いや位相の違いがある場合でも、多チャネルの情報を利用することで効率的に信号を再構成することができるようにしたことを特徴としている。ここで、全実施形態に共通する信号モデルについて記述しておく。
【００３５】
本発明では、信号モデルを共通信号ｓｃとともに振幅係数ａｋｃ及び群遅延βｋｃも考慮した信号モデルに拡張した上でこれら２つの因子をチャネルごとに推定することによって、より正確な信号再構成を行うことを特徴としている。具体的には、各信号の離散フーリエ変換（DFT）が
ｘｋ＝ｘｃｋ＋ｘｄｋ      …（１３）
という和の形式で書けるとすると、式（１４）で近似することができるという信号モデルである。
【数４】
0005785903.tif 000005
【００３６】
ただし、Λは式（１５）で表される。πは円周率、ｊは虚数単位である。
【数５】
0005785903.tif 000006
【００３７】
上記の式を合わせると、元の信号の離散フーリエ変換は式（１６）で近似的に記述できることを前提としている。
【数６】
0005785903.tif 000007
【００３８】
［実施形態１］
図３を参照して、実施形態１について説明する。実施形態１において観測装置と送信装置と受信装置は既述の従来技術と同じ構成および処理であるから説明を省略する。本発明の実施形態は、従来の復元装置９５３に替えて復元装置１００を備える。
【００３９】
以下では、復元装置１００が、受信装置９５１から受けた信号ベクトルｙｋ、圧縮行列Ａｋ、観測行列Φｋ、基底行列Ψを知っているものとする（例えば、復元装置１００が事前に各センサーから各行列についての情報を受信している）。ただし、ｋ＝１、２、…、Ｋであり、下記では基底行列Ψを離散フーリエ変換に対応する行列とする。
【００４０】
実施形態１では、復元装置１００は、仮伸長部１０１、差異推定部１０３、共通信号推定部１０５、残信号推定部１０７を有する。
【００４１】
仮伸長部１０１は、信号ベクトルｙｋから仮の信号の再構成を行う。この処理には、通常の圧縮センシングのアルゴリズムが用いられる。これによって仮の復元信号ｓ~１、ｓ~２、…、ｓ~Ｋを得る。
【００４２】
差異推定部１０３は、ｓ~１、ｓ~２、…、ｓ~Ｋを用いて下記のように振幅及び位相の差異を推定する。具体的には、差異推定部１０３は、式（１７）で表される量を最小化することで各チャネル間の相対的な振幅の差ａｋｍと群遅延量ｂｋｍを推定する。
【数７】
0005785903.tif 000008
【００４３】
ここで[DFT]および、[IDFT]はそれぞれ離散フーリエ変換、逆離散フーリエ変換に対応する変換行列である。上記の量はａｋｍとｂｋｍを別々に変化させることで最適化することが可能である。
【００４４】
差異推定部１０３は、さらに、
ａｋ＝（ａｋｍ＋１／ａｍｋ）／２、    βｋ＝（βｋｍ−βｍｋ）／２    …（１８）
とする。この処理は推定誤差を軽減するためのものであるが、この処理を行わず、
аｋ＝аｋｍ、    βｋ＝βｋｍ
としてもよい。
また、аｋｍは、例えば
а１３＝ａ１２＋а２３
を近似的に満たすので、この観点からаｋｍに関するコスト関数
||ａｋｍ−（ａｋｊ＋ａｊｍ）||２      …（１９）
を最小化するようにａｋｍの組を調整し推定誤差を軽減してもよい。また上の関係式からさらに４項以上のаｋｍの組によるコスト関数を立てることができるので、このようなコスト関数を最小化するようにａｋｍの組を調整し推定誤差を軽減してもよい。
【００４５】
差異推定部１０３は、さらに、
Σｋａｋ＝１、Σｋｂｋ＝０      …（２０）
となるように正則化を行う。この正則化によって得られた係数をａｋｃ、ｂｋｃとする。ただし、上述の方法の正則化に限らず、或る定められたベクトルと各チャネル間の相対的な振幅の差と群遅延量の両方の情報を含めることができれば、どのような正則化を用いてもよい。例えば、特定のチャネルを共通ベクトルとみなして、そのベクトルからの振幅の差と群遅延量をａｋｃ、ｂｋｃとしてもよい。
【００４６】
次に、共通信号推定部１０５は、上記の係数ａｋｃ、ｂｋｃ及び観測された信号ｙ１、ｙ２、…、ｙＫを用いて式（２１）を計算する。ただし、［・］†は行列［・］の一般化逆行列を表す。
【数８】
0005785903.tif 000009
【００４７】
残信号推定部１０７は、上記共通信号推定部１０５で求めたｘ'ｃを用いて、式（２２）に従って観測信号ｙｋとの差分ｙｄｋを得る。
【数９】
0005785903.tif 000010
【００４８】
このｙｄｋとＡｋについて通常の圧縮センシングのアルゴリズムを適用することによってｓ'ｄ１、ｓ'ｄ２、…、ｓ'ｄＫを得る。さらにこれらとｘ'ｃ、ａｋｃ、ｂｋｃを用いることで、元の信号の推定値ｓ'ｋを得る（式（２３））。
【数１０】
0005785903.tif 000011
【００４９】
以上によって復元装置１００は元の信号よりも少ない情報ｙ１、ｙ２、…、ｙＫから元の信号の推定値ｓ'１、ｓ'２、…、ｓ'Ｋを得ることができる。
【００５０】
［実施形態２］
実施形態１では全ての処理を各チャネルについて一度だけ行うものであったが、一連の処理を反復的に行うことによって推定精度を向上させることができる。実施形態２を、図４を参照して説明する。なお、復元装置１００以外の構成は実施形態１と同様である。
【００５１】
仮伸長部１０１と差異推定部１０３の処理は実施形態１と同様である。
共通信号推定部１０５は実施形態１と異なり、式（２４）を計算することによって共通信号を求める。
【数１１】
0005785903.tif 000012
【００５２】
ただし、ここでＱｋは、初期値がサイズＭ×Ｍの単位行列であり、後述する残信号基底推定部１０９によって更新される。
【００５３】
残信号基底推定部１０９は、残信号推定部１０７と同様の手続きで個別信号ｘｄｋを推定し、さらに下記のような処理を行う（例えば非特許文献１）。この処理の説明に先立ち、前提事項について説明を加えておく。
【００５４】
行列Ｂｋ∈ＣＭ×Ｍを、当該行列の列ベクトル（或いは行ベクトル）によって張られる数ベクトル空間に信号ｙｋが含まれるような基底行列とする（換言すれば、信号ｙｋが行列Ｂｋの列ベクトル（或いは行ベクトル）の線形結合で表現される）。このとき基底行列Ｂｋは下記のように２つの部分に分けることができる。
Ｂｋ＝[ＡｋΩ(ｋ)  Ｑｋ]      …（２５）
【００５５】
ここでＡkΩ(k)は圧縮行列Ａｋの一部を列単位で切り取った得た行列で、Ω(k)は切り取った列を特定するインデックスを表している。なお、ＱｋはＡｋΩ(k)を構成するすべての列ベクトルに直交した列ベクトルからなる基底行列である。。
【００５６】
残信号基底推定部１０９は、予め決められた個数Ｄｋ（ｋ＝１、２、…、Ｋ）に従って圧縮行列Ａｋのうちｙｄｋを２乗誤差の観点から良く記述しているＤｋ個の列を特定するインデックスΩ(ｔｋ)を決定し、これと現在のΩ(k)との積集合を新たにΩ(k)として、行列Ｑｋを更新する。ただしΩ(k)の初期値は空集合である。また、Ｑｋの初期値はＭ×Ｍの単位行列である。下記に残信号基底推定部１０９の処理の具体例を示す。
１．ｙｄｋを用いて通常の圧縮センシングのアルゴリズムを実行してｘｄｋを取得する。すなわち、
||ｙｄｋ―Ａｋｘｄｋ||２＜εの下で量||ｘｄｋ||１を最小化する、という問題を解く。
２．ステップ１．で得られたｘｄｋの要素のうちそれら各要素の絶対値の大きさを比較する。このうち上位Ｄｋ個のｘｄｋの要素のインデックスをΩ(ｔｋ)として保持する。
３．Ω(ｔｋ)とΩ(k)との積集合を新たなΩ(k)として保持する。
４．ＡkΩ(k)に対する直交基底Ｑｋを計算する。
【００５７】
このようにして得られる行列Ｑｋは大きさがＭ×（Ｍ−|Ωｋ|）の行列となる。Ｄｋは反復のたびに変更してもよい。また、積集合の代わりに和集合を用いてΩ(k)を更新してもよい。この場合はΩ(k)の大きさが過大にならないように、上限Ｅｋを設けることが望ましい。
【００５８】
残信号基底推定部１０９は、行列Ｑｋを共通信号推定部１０５に渡す。以下、あらかじめ定めた反復回数に達するまで、残信号基底推定部１０９と共通信号推定部１０５との間で共通信号ｘ'ｃと基底行列Ｑｋの更新を行う。
【００５９】
反復終了後に残信号推定部１０７はａｋｃ、ｂｋｃ、ｘ'ｃ及び観測行列Φｋ、基底行列Ψ、観測結果ｙｋを用いて、実施形態１と同様に残差信号を推定し、最後に再構成した復元信号ｓ'１、ｓ'２、…、ｓ'Ｋを出力する。
【００６０】
［実施形態３］
上記実施形態２において共通信号推定部１０５の処理（式（２４））を下記のように置き換えることも可能である。すなわち、
ｙ~ｋ＝ｙｋ−ｙ~ｄｋ      …（２６ａ）
という減算によって残差信号の影響を除去した上で共通信号
【数１２】
0005785903.tif 000013
を推定する。ただし、ｙ~ｄｋは、実施形態２における残信号基底推定部１０９で得られたＡkΩ(k)の基底のみで残差ｙｄｋを近似したベクトルである。すなわち、ｙ~ｄｋは、
ｘ~ｄｋΩ（ｋ）＝ＡｋΩ（ｋ）†ｙｄｋ
として
ｙ~ｄｋ＝Φｋｘ~ｄｋ
で表される。
【００６１】
［実施形態４］
実施形態２または実施形態３では、仮伸長部１０１による処理は始めに一度だけ実行されるだけであるが、図５に示すように、個別信号を除去した信号を用いて仮伸長部１０１で再び通常の圧縮センシングのアルゴリズムで仮の復元信号ｓ~１、ｓ~２、…、ｓ~Ｋを得ることも可能である。このようにして、残信号の基底更新ごとに新たに仮の復元信号を用いて差異推定部１０３によってａｋcやｂｋcを推定することができる。
【００６２】
具体的には実施形態２で定義した残信号の基底を取り除く行列Ｑｋを用いて
Ｑｋｙｋ＝ＱｋＡｋｘｋ      …（２７）
を考えれば、ｙｋの代わりにＱｋｙｋを観測結果、ＱｋＡｋを観測の際に用いた圧縮行列であるとみなすことによって通常の圧縮センシングのアルゴリズムによって仮の復元信号ｓ~ｋの再構成を行うことができる。
【００６３】
また、式（２７）による上記の方法に代えて、単なる減算処理を行ってもよい。具体的には
ｙｃｋ＝ｙｋ−Ａｋｘ'ｄｋ      …（２８）
を観測結果として用い、Ａｋはそのまま用いて、仮伸長部１０１で通常の圧縮センシングのアルゴリズムを用いることで仮の復元信号ｓ~ｋを得ることができる。
再度得られたこれら仮の復元信号ｓ~１、ｓ~２、…、ｓ~Ｋは差異推定部１０３に送られる。
【００６４】
［実施形態５］
上記実施形態２〜４では反復処理の回数は適当な固定値であった。しかし、終了条件判定部（図示せず）を設けることで残差信号の基底更新が行われなくなったことをもって反復処理を終了することも可能である。すなわち、今回処理の残信号基底推定部１０９の出力（行列Ｑｋ）が前回処理の残信号基底推定部１０９の出力（行列Ｑｋ）と同じであれば反復を終了し残信号推定部１０７に結果を渡すようにすることが可能である。
【００６５】
［実施形態６］
上記実施形態１〜４について、送信装置や受信装置を用いるのではなく、観測装置で得た結果ｙｋを外部記憶装置に記録しておき、当該外部記憶装置から得た観測結果ｙｋを用いて復元装置で復元してもよい。また、各装置の間に適当な外部記憶装置を介在させ、各装置で得られた情報を外部記憶装置に保持させる処理を含めてもよい。
【００６６】
［実施形態７］
実施形態１〜６ではｓｄｋが疎になる基底が離散フーリエ変換であるとしてきた。しかし、これを一般的な基底に拡張しても、フーリエ領域（周波数領域）との間に変換行列を規定することができるので、下記のように僅かな変更を実施形態１〜６を適用することができる。
【００６７】
以下ではｚｄｋ＝Ψｓｄｋとし、ｚｄｋをｓｄｋの基底Ψでの表現とし、これが疎であるとする。ｘｋ＝[DFT]ｓｋは引き続き、ｓｋの離散フーリエ変換とする。実施形態１〜６の共通信号推定部１０５では、Ａｋ＝ΦｋΨのかわりにＯｋ＝Φｋ[IDFT]を使えばｘｃはフーリエ領域の表現になる。ｙｋはそのままでよい。残信号推定部１０７と残信号基底推定部１０９では残差を
【数１３】
0005785903.tif 000014
としたが、これを
【数１４】
0005785903.tif 000015
に置き換えればよい。これらの推定部で通常の圧縮センシングのアルゴリズムを適用する際には圧縮行列をＡｋ＝ΦｋΨとし、
||ｙｄｋ―Ａｋｚ||２＜εの下で量||ｚ||１を最小化する
という条件の問題を解いて得られるｚをｚｄｋとする。
【００６８】
図６は従来法と本発明による手法との性能差を表わしている。ＴＥＣＣアルゴリズムは非特許参考文献１による従来のアルゴリズムである。比較結果からわかるように実際の信号に対して本発明による手法は従来法よりもＳＮＲが良い。
【００６９】
＜復元装置のハードウェア構成例＞
上述の実施形態に関わる復元装置は、ＣＰＵ（Central Processing Unit）やＤＳＰ（Digital Sygnal Processor）〔キャッシュメモリなどを備えていてもよい。〕、メモリであるＲＡＭ（Random Access Memory）やＲＯＭ（Read Only Memory）と、ハードディスクである外部記憶装置、並びにこれらのＣＰＵやＤＳＰ、ＲＡＭやＲＯＭ、外部記憶装置間のデータのやり取りが可能なように接続するバスなどを備えている。また必要に応じて、音響再生装置に、ＣＤ−ＲＯＭなどの記憶媒体を読み書きできる装置（ドライブ）などを設けるとしてもよい。
【００７０】
復元装置の外部記憶装置には、上述の復元処理のためのプログラム並びにこのプログラムの処理において必要となるデータ（圧縮行列、観測行列、基底行列、回転行列）などが記憶されている〔外部記憶装置に限らず、例えばプログラムを読み出し専用記憶装置であるＲＯＭに記憶させておくなどでもよい。〕。また、これらのプログラムの処理によって得られるデータなどは、ＲＡＭや外部記憶装置などに適宜に記憶される場合がある。データやその格納領域のアドレスなどを記憶する記憶装置を単に「記憶部」と呼ぶことにする。
【００７１】
復元装置の記憶部には、圧縮信号ｙｋをチャネルごとに独立に元の信号（仮の復元信号）を復元するためのプログラム、仮の復元信号ｓ~ｋから、振幅係数ａｋｃと群遅延ｂｋｃを求めるためのプログラム、圧縮信号ｙｋと、振幅係数ａｋｃと、群遅延ｂｋｃとを用いて共通信号の周波数領域表現を求めるためのプログラム、圧縮信号ｙｋと、振幅係数ａｋｃと、群遅延ｂｋｃと、共通信号の周波数領域表現ｘ'ｃを用いて、圧縮信号ｙｋと共通信号との差分ｙｄｋを得て、さらに当該差分ｙｄｋを復元してｓ'ｄｋを得て、当該ｓ'ｄｋと共通信号の時間領域表現との和によって、元の信号の推定値ｓ'ｋを得るためのプログラムなどが記憶されている。
【００７２】
復元装置では、記憶部に記憶された各プログラムとこの各プログラムの処理に必要なデータが必要に応じてＲＡＭに読み込まれて、ＣＰＵで解釈実行・処理される。この結果、ＣＰＵが所定の機能（仮伸長部、差異推定部、共通信号推定部、残信号基底推定部、残信号推定部など）を実現することで上述の音響再生が実現される。
【００７３】
＜補記＞
本発明は上述の実施形態に限定されるものではなく、本発明の趣旨を逸脱しない範囲で適宜変更が可能である。また、上記実施形態において説明した処理は、記載の順に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されるとしてもよい。
【００７４】
また、上記実施形態において説明したハードウェアエンティティ（復元装置）における処理機能をコンピュータによって実現する場合、ハードウェアエンティティが有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記ハードウェアエンティティにおける処理機能がコンピュータ上で実現される。
【００７５】
この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。具体的には、例えば、磁気記録装置として、ハードディスク装置、フレキシブルディスク、磁気テープ等を、光ディスクとして、ＤＶＤ（Digital Versatile Disc）、ＤＶＤ−ＲＡＭ（Random Access Memory）、ＣＤ−ＲＯＭ（Compact Disc Read Only Memory）、ＣＤ−Ｒ（Recordable）／ＲＷ（ReWritable）等を、光磁気記録媒体として、ＭＯ（Magneto-Optical disc）等を、半導体メモリとしてＥＥＰ−ＲＯＭ（Electronically Erasable and Programmable-Read Only Memory）等を用いることができる。
【００７６】
また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。
【００７７】
このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。
【００７８】
また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、ハードウェアエンティティを構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。
【図１】
0005785903.tif 000020
【図２】
0005785903.tif 000021
【図３】
0005785903.tif 000022
【図４】
0005785903.tif 000023
【図５】
0005785903.tif 000024
【図６】
0005785903.tif 000025
【図７】
0005785903.tif 000026