JP5945520B2 - 圧縮信号復元装置、圧縮信号復元方法、プログラム - Google Patents

Description

この発明は時系列信号およびその他の信号を保存或いは伝達するための符号化および復号化に関し、特に圧縮信号を復元する圧縮信号復元装置、圧縮信号復元方法、プログラムに関する。

多数のセンサを無線ネットワークで接続したセンサネットワークでは、各センサで取得したデータの送受信が頻繁に行われる。センサは処理能力が低く、外部電源からの電力供給が得られないようなものが用いられることも多いので、センサの省電力化のために、センサが送信する情報量を削減すること、およびセンサ自身で行う演算量を低減する手法が重要とされている。

その一つが、圧縮センシングという手法である。これは、用途が同じであり地理的に近接したセンサが計測した結果は類似性が高いことを利用して、各センサの情報を、他のセンサと共通する情報とそれ以外の情報とで表すことにより、情報量を圧縮するものである。非特許文献１では、単一センサの圧縮センシングの方法が提案されている。ここでは、送信装置は、測定の対象となる標本数Ｎの信号ｘを、Ｍ×Ｎの観測行列Ψを用いて線形変換することにより、ｘの次元に比して少ない次元Ｍの情報ｙ＝Ψｘに圧縮し、送信する。受信装置は、元の信号ｘの基底Φでの表現ｓ＝Φ^−１ｘの非０要素が非常に少ない（以下、「疎である」と表現）という仮定のもとに復元を行うことで、元の信号ｘと厳密に等しい信号を復号しようとするものである。

また、多数のセンサからなる無線ネットワーク上で圧縮センシングを行う方法として、非特許文献２では、各第ｉ番目のセンサの観測信号ｘ_ｉが

で表される信号の圧縮センシングが提案されている。ここで、

である。ａ_ｉは第ｉ番目のセンサの観測信号の振幅、θ_ｉは第ｉ番目のセンサの観測信号の位相を表す係数であり、ｚ_ｃは全てのチャネル間の共通信号である。また、

であり、ｊは虚数単位（ｊ^２＝−１）であり、ＤＦＴは離散フーリエ変換、ＩＤＦＴはＤＦＴの逆変換である。つまり、このモデルは、各センサが観測した信号の共通部分ｚ_ｃが群遅延および全体の振幅の減衰があるが形状が同じである信号ｚ_ｃ，ｉであるようなモデルである。

非特許文献２では、（１）ｘ_ｉを推定し、（２）推定したｘ_ｉと観測信号を用いて振幅ａ_ｉと位相θ_ｉを推定し、（３）推定した振幅ａ_ｉと位相θ_ｉを用いて共通信号ｚ_ｃ，ｉを推定し、（４）推定した共通信号と観測信号を用いて残りの信号ｚ_ｄ，ｉを推定する、という順で各信号を求めていく。（１）〜（４）の最適化では、それぞれ個別の最適化の式に従って最適化を行う。

E. Hale, W. Yin, and Y. Zhang, "Fixed-point continuation for `1-minimization: Methodology and convergence," SIAM Journal on Optimization, vol. 19, no. 3, pp. 1107-1130, 2008. 白木善史、鎌本優、守谷健弘、「減衰および遅延のある共通要素を持った多チャンネル信号の復元」、電子情報通信学会大会講演論文集、電子情報通信学会、平成24年8月28日、第2012巻、p.62

非特許文献１は、単一センサのための圧縮センシングの方法であるため、多数のセンサからなるセンサネットワークにそのまま適用することができない。

非特許文献２によれば、多数のセンサの圧縮センシングを行うことができるが、必ずしも元の信号ｘ_ｉ’と厳密に一致する信号ｘ_ｉを復号できないという問題がある。これは非特許文献２の方法では各信号の推定を、既に推定した他の信号の推定結果を用いて最適化するため、先に推定した信号の推定精度が低いと、残りの信号の推定精度にも影響することによる。加えて、スパース項ｚ_ｄ，ｉの推定は共通信号を求めた後に行っていたため式（１）で得られた推定結果が元の観測信号に最小二乗の意味で最適である保証はなかった。そこで、本発明では、元の信号を従来よりも精度良く復元することが出来る圧縮信号復元装置を提供することを目的とする。

本発明の圧縮信号復元装置は、Ｎｃ個のセンサからの圧縮信号を復元する圧縮信号復元装置であって、復号部を含む。

復号部は、インデクスｉを１以上Ｎｃ以下の整数とし、全てのセンサに共通の信号を信号ｚ_ｃとし、第ｉ番目のセンサが観測した信号の共通部分を共通信号ｚ_ｃ，ｉと表すものとし、第ｉ番目のセンサが観測した信号のうち共通信号ｚ_ｃ，ｉを除いた残りの信号を信号ｚ_ｄ，ｉと表すものとし、

で表される第ｉ番目のセンサの観測信号ｘ_ｉを所定の観測行列Ψ_ｉにより圧縮した信号ｙ_ｉ＝Ψ_ｉｘ_ｉを入力とし、ｅ_ｃ，ｉを第ｉ番目のセンサの補助共通信号とし、ａ_ｃ，ｉを補助共通信号ｅ_ｃ，ｉに対する第ｉ番目のセンサの共通信号要素の振幅を表す係数とし、θ_ｃ，ｉをｚ_{ｃ，ｉ＋１}に対する第ｉ番目のセンサの共通信号要素の位相を表す係数とし、Φを予め定められた基底行列とし、共通信号ｚ_ｃの基底Φでの表現ｓ_ｃ＝Φ^−１ｚ_ｃが疎ではなく、信号ｚ_ｄ，ｉの基底Φでの表現ｓ_ｄ，ｉ＝Φ^−１ｚ_ｄ，ｉが全てのセンサのインデクスｉに対して疎であるものとし、‖・‖_２をベクトルの各要素の二乗和であるＬ２ノルムとし、‖・‖_１をベクトルの各要素の絶対値和であるＬ１ノルムとし、Ωをインデクスｉの集合、すなわちΩ＝｛１，２，・・・，Ｎｃ｝とし、 ｗ_ｉを、予め定めた個別信号ｓ_ｄ，ｉに対する重みとし、∀を全称記号とし、ｊを虚数単位、ＤＦＴを離散フーリエ変換、ＩＤＦＴをＤＦＴの逆変換とし、記号Λを

と定義し、

の制約条件のもとで、

が最小となる観測信号ｘ_ｉを求め、元の信号の復号信号として出力する。

本発明によれば、元の信号を従来よりも精度良く復元することが出来る。

多数のセンサによる圧縮センシングの概要を示す図。 本発明の実施例１の圧縮信号復元装置の構成を示すブロック図。 本発明の実施例１の圧縮信号復元装置の動作を示すフローチャート。 本発明の実施例１の圧縮信号復元装置の動作を示すフローチャート。 本発明の実施例１の圧縮信号復元装置の信号推定部の構成を示すブロック図。 本発明の実施例１の圧縮信号復元装置の信号推定部の動作を示すフローチャート。 本発明の実施例２の圧縮信号復元装置の構成を示すブロック図。 本発明の実施例２の圧縮信号復元装置の動作を示すフローチャート。 本発明の実施例２の圧縮信号復元装置の信号推定部の構成を示すブロック図。 本発明の実施例２の圧縮信号復元装置の信号推定部の動作を示すフローチャート。 本発明の実施例３の圧縮信号復元装置の構成を示すブロック図。 本発明の実施例３の圧縮信号復元装置の動作を示すフローチャート。 本発明の実施例３の圧縮信号復元装置の信号推定部の構成を示すブロック図。 本発明の実施例３の圧縮信号復元装置の信号推定部の動作を示すフローチャート。 本発明の実施例４の圧縮信号復元装置の構成を示すブロック図。 本発明の実施例４の圧縮信号復元装置の動作を示すフローチャート。 本発明の実施例４の圧縮信号復元装置の信号推定部の構成を示すブロック図。 本発明の実施例４の圧縮信号復元装置の信号推定部の動作を示すフローチャート。 本発明の実施例５の圧縮信号復元装置の構成を示すブロック図。 本発明の実施例５の圧縮信号復元装置の動作を示すフローチャート。 本発明の実施例５の圧縮信号復元装置の信号推定部の構成を示すブロック図。 本発明の実施例５の圧縮信号復元装置の信号推定部の動作を示すフローチャート。

以下、本発明の実施の形態について、詳細に説明する。なお、同じ機能を有する構成部には同じ番号を付し、重複説明を省略する。

＜圧縮センシングの概要＞
以下、図１を参照して多数のセンサによる圧縮センシングについて説明する。図１は多数のセンサによる圧縮センシングの概要を示す図である。図１に示すように、Ｎｃ個（Ｎｃは１以上の整数）のセンサ（第１センサ９−１、第２センサ９−２、…、第Ｎｃセンサ９−Ｎｃ）と１個の中央受信機８がある。第１センサ９−１、第２センサ９−２、…、第Ｎｃセンサ９−Ｎｃはそれぞれ異なった信号を計測し、中央受信機８へ送信する。ここで各センサがある単位時間に取得する情報を複素ベクトルで表すこととし、このベクトルの長さはＮ（Ｎは正の整数）であるとする。すなわちある単位時間にＮｃ個のセンサが取得した情報は一般に各々ｘ’_ｉ∈Ｃ＾Ｎ、ｉ＝１、２、…、Ｎｃ（インデクスｉは１以上Ｎｃ以下の整数）と書けるとする。

なお、本発明の観測信号のモデルは、非特許文献２でも用いられている

で表せるモデルである。ここで、第ｉ番目のセンサが観測した信号の共通部分である共通信号ｚ_ｃ，ｉは、全てのセンサに共通の信号ｚ_ｃに群遅延および全体の振幅の減衰が加わったような信号

である。ｚ_ｄ，ｉは、第ｉ番目のセンサが観測した信号のうち共通信号ｚ_ｃ，ｉを除いた残りの信号である。ここで、Φを予め定められた基底行列とし、共通信号ｚ_ｃの基底Φでの表現ｓ_ｃ＝Φ^−１ｚ_ｃが疎ではなく、共通信号以外の残りの信号ｚ_ｄ，ｉの基底Φでの表現ｓ_ｄ，ｉ＝Φ^−１ｚ_ｄ，ｉが全てのセンサのインデクスｉに対して疎であることを前提とする。以下、ｓ_ｄ，ｉを「個別信号」と呼ぶ。

図２に示すように、各センサは測定部９１と圧縮部９２とで構成される観測装置９とみなすことができる。測定部９１は、単位時間ごとに観測した長さＮの信号ｘ＾_ｉを圧縮部９２へ送信する。圧縮部９２は、長さＮの信号ｘ＾_ｉを長さＭの圧縮された信号ｙ_ｉ∈Ｃ^Ｍに変換し出力する。ただし、ＭはＮより小さな正の整数である。また、この変換は線形であるものとする。すなわち、圧縮部９２は、予め定められた観測行列Ψ_ｉ∈Ｃ^Ｍ×Ｎを用いて

により圧縮された信号ｙ_ｉを求めて送信する。

中央受信機８は入力された信号ｙ_ｉから元の信号ｘ＾_ｉを復元する圧縮信号復元装置１を有する。図２に示すように、圧縮信号復元装置１は復号部１１を含む。本発明は、復号部１１で各々の信号ｘ＾_ｉを推定する技術に関する。

＜本発明のポイント＞
式（１）で表される観測信号のモデルにおいて、各センサ同士の関係に着目すれば

のように書くことができる。ここで、ｅ_ｃ，ｉは第ｉ番目のセンサの補助共通信号であり、ａ_ｃ，ｉは補助共通信号ｅ_ｃ，ｉに対する第ｉ番目のセンサの共通信号要素の振幅を表す係数であり、θ_ｃ，ｉはｚ_{ｃ，ｉ＋１}に対する第ｉ番目のセンサの共通信号要素の位相を表す係数である。この関係を利用して、本発明の復号部１１は、入力された信号ｙ_ｉを用いて、

を、

のような制約条件のもとで最小化する（最小となる観測信号ｘ_ｉを求める）ことにより、元の信号ｘ＾_ｉの推定値ｘ_ｉを復号する（Ｓ１１）。この最小化は、ADMM(Alternating Direction Method of Multiplier,参考非特許文献１)を応用することにより行うことができる。すなわち、

で表される目的関数Ｈを、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉのそれぞれについて最小化することにより、元の信号ｘ＾_ｉの推定値ｘ_ｉを復号する。ここで、‖・‖_２はベクトルの各要素の二乗和（Ｌ２ノルム）であり、‖・‖_１はベクトルの各要素の絶対値和（Ｌ１ノルム）であり、Ωはインデクスｉの集合、すなわちΩ＝｛１，２，・・・，Ｎｃ｝である。Ω⁻はＮｃを除くインデクスｉの集合、すなわちΩ⁻＝｛１，２，・・・，Ｎｃ−１｝である。Φは予め定められた基底行列とする。また、λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉはそれぞれ、式（１０），（１１），（１２），（１３）の制約条件に対する変数（未定乗数）である。ｗ_ｉは、個別信号ｓ_ｄ，ｉに対する重みであり、予め定めた定数である。α，β，γ，σは、それぞれ最小化問題を解く際の収束速度を調整するための重みであり、予め設定された定数（パラメータ）である。α，β，γ，σのかかる項は、オーバーフィッティングを避けるためのペナルティ項の役割を果たす。また、Ｔは行列またはベクトルの転置を表す。
（参考非特許文献１）D. Gabay and B. Mercier, “A dual algorithm for the solution of nonlinear variational problems via finite element approximation,” Computers & Mathematics with Applications, vol. 2, no. 1, pp. 17 - 40, 1976.

図２、図３、図４を参照して実施例１の圧縮信号復元装置１について説明する。図２は本実施例の圧縮信号復元装置１の構成を示すブロック図である。図３、図４は本実施例の圧縮信号復元装置１の動作を示すフローチャートである。前述したように、各センサは測定部９１と圧縮部９２とで構成される観測装置９とみなすことができる。圧縮信号復元装置１は、復号部１１を含み、復号部１１は、信号推定部１１１と、パラメータ記憶部１１２と、乗数更新部１１３と、終了判定部１１４とを含む。

＜パラメータ記憶部１１２＞
パラメータ記憶部１１２には、予め定められた観測行列Ψ_ｉ、基底行列Φ、センサの総数Ｎｃを記憶しておく。また、予め定められた目的関数のパラメータα，β，γ，σと、個別信号に対する重みｗ_ｉとを記憶しておく。Ψ_ｉ，Φ，Ｎｃ，α，β，γ，σ，ｗ_ｉは固定値である。また、パラメータ記憶部１１２には、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉの初期値と、目的関数における制約条件に対する未定乗数λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉの初期値をあらかじめ記憶しておく。これらの初期値には、ランダムな値を設定しておけば良い。

＜信号推定部１１１＞
信号推定部１１１は、上記式（１４）の目的関数Ｈを、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉのそれぞれについて最小化することにより、各信号および係数の推定値を求め、パラメータ記憶部１１２に記憶された信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉの値を更新する（Ｓ１１１）。信号推定部１１１の具体的な処理は後述する。

＜乗数更新部１１３＞
乗数更新部１１３は、信号推定部１１１で推定された信号および係数を用いて、式（１４）の目的関数における制約条件に対する未定乗数λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉを更新する（Ｓ１１３）。具体的には、下記式により各乗数の値を更新する。

＜終了判定部１１４＞
終了判定部１１４は、所定の終了条件を充たしたか否かを判定し、所定の条件を充たしていない場合は、信号推定部１１１、乗数更新部１１３を制御して、各信号および係数の推定と乗数の更新処理の繰り返し制御を実行する（Ｓ１１４）。所定の条件を充たした場合は、求めた信号ｘ_ｉを出力し、処理を終了する。

所定の終了条件は、例えば「信号推定部１１１の繰り返し回数ｋが予め定めた回数Ｉｔ_ｍａｘに到達したか否か」を用いる。終了判定部１１４は、ｋがＩｔ_ｍａｘを超えた場合に「所定の条件を充たした」と判定する。あるいは、「信号推定部１１１で更新される信号ｘ_ｉの更新量が所定の閾値εより小さくなったか否か」を用いることができる。この場合、終了判定部１１４は、ｘ_ｉを信号推定部１１１で更新された後の全体信号とし、更新前の全体信号を^〜ｘ_ｉとしたとき、

を充たす場合に、「所定の条件を充たした」と判定する。

＜信号推定部１１１の具体処理＞
以下、図５、図６を参照して、信号推定部１１１の詳細について説明する。図５は本実施例の圧縮信号復元装置１の信号推定部１１１の構成を示すブロック図である。図６は本実施例の圧縮信号復元装置１の信号推定部１１１の動作を示すフローチャートである。信号推定部１１１は、全体信号推定部１１１１と、要素信号推定部１１１２と、共通信号推定部１１１３と、個別信号推定部１１１４と、振幅位相推定部１１１５とを含む構成である。なお、これらの構成部はどんな順序で各処理を実行しても良いが、例えば、全体信号推定部１１１１、要素信号推定部１１１２、共通信号推定部１１１３、個別信号推定部１１１４、振幅位相推定部１１１５の順に各処理を実行すれば好適である。

＜全体信号推定部１１１１＞
全体信号推定部１１１１は入力された信号ｙ_ｉと、パラメータ記憶部１１２に記憶されたｚ_ｃ，ｉ，ｚ_ｄ，ｉ，λ_ｉ，α，Ψ_ｉ，Ｎｃを入力として受け取り、ｙ_ｉに対応する復号信号の推定値ｘ_ｉを求める（Ｓ１１１１）。具体的には、全体信号推定部１１１１は、

をｘ_ｉについて最小化することによりｘ_ｉを求め、パラメータ記憶部１１２に記憶する（Ｓ１１１１）。この最小化はたとえば、

のように実行される。全体信号推定部１１１１は、ｉ＝１，２，…，Ｎｃそれぞれについて最小化を行う（Ｓ１１１１）。また、Ｉは大きさＮの単位行列である。なお、式（２０）は、式（１４）のうち、推定したい変数ｘ_ｉに関連する項のみを取り出したものである。

＜要素信号推定部１１１２＞
要素信号推定部１１１２は、パラメータ記憶部１１２に記憶されたｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，λ_ｉ，α，μ_ｃ，ｉ，ａ_ｃ，ｉ，ｅ_ｃ，ｉ，βを用いて、センサ毎の要素信号の推定値ｚ_ｃ，ｉとｚ_ｄ，ｉを求め、パラメータ記憶部１１２に記憶する（Ｓ１１１２）。

具体的には、要素信号推定部１１１２は、ｉ＝１については

をｚ_ｃ，ｉについて最小化することによりｚ_ｃ，ｉを求める（Ｓ１１１２）。この最小化は、

のように行われる。また、要素信号推定部１１１２は、ｉ＝２，３，…，Ｎｃ−１については

をｚ_ｃ，ｉについて最小化することによりｚ_ｃ，ｉを求める（Ｓ１１１２）。この最小化は、

のように行われる。ここで、＊はエルミート転置行列を表す。
また、要素信号推定部１１１２は、ｉ＝Ｎｃについては

をｚ_ｃ、ｉについて最小化することによりｚ_ｃ，ｉを求める（Ｓ１１１２）。この最小化は、

のように行われる。なお、式（２２）（２４）（２６）は、式（１４）のうち、推定したい変数ｚ_ｃ、ｉに関連する項のみを取り出したものである。

また、要素信号推定部１１１２はｉ＝１，２，…，Ｎｃについて、

をｚ_ｄ，ｉについて最小化し、ｚ_ｄ，ｉを出力する（Ｓ１１１２）。この最小化はたとえば、

のように行われる。なお、式（２８）は、式（１４）のうち、推定したい変数ｚ_ｄ，ｉに関連する項のみを取り出したものである。

要素信号推定部１１１２は、パラメータ記憶部１１２に記憶されたｚ_ｃ，ｉ及びｚ_ｄ，ｉの値を、上述の最小化により求めたｚ_ｃ，ｉ及びｚ_ｄ，ｉの値にそれぞれ更新する（Ｓ１１１２）。

＜共通信号推定部１１１３＞
共通信号推定部１１１３は、パラメータ記憶部１１２に記憶されたｚ_ｃ，ｉ，θ_ｃ，ｉ，μ_ｃ，ｉ，ａ_ｃ，ｉ，ｅ_ｃ，ｉ，β，ν_ｃ，ｉ，Λ，ｚ_{ｃ，ｉ＋１}，γを用いて、補助共通信号ｅ_ｃ，ｉを推定する（Ｓ１１１３）。具体的には、共通信号推定部１１１３は、

をｅ_ｃ，ｉについて最小化することにより、ｅ_ｃ，ｉを求める（Ｓ１１１３）。この最小化は、

を求めることにより行われる。共通信号推定部１１１３は、これを、ｉ＝１，２，…，Ｎｃ−１それぞれについて求める（Ｓ１１１３）。なお、式（３０）は、式（１４）のうち、推定したいｅ_ｃ，ｉに関連する項のみを取り出したものである。

共通信号推定部１１１３は、パラメータ記憶部１１２に記憶されたｅ_ｃ，ｉの値を、上述の最小化により求めたｅ_ｃ，ｉの値に更新する（Ｓ１１１３）。

＜個別信号推定部１１１４＞
個別信号推定部１１１４は、パラメータ記憶部１１２に記憶されたξ_ｄ，ｉ，ｚ_ｄ，ｉ，Φ，ｓ_ｄ，ｉ，σ，ｗ_ｉを用いて、個別信号ｓ_ｄ，ｉを推定する（Ｓ１１１４）。具体的には、個別信号推定部１１１４は、

をｓ_ｄ，ｉについて最小化することにより個別信号ｓ_ｄ，ｉを求める（Ｓ１１１４）。この最小化は、

を求めることにより行われる。ここで、ｓｇｎは要素ごとの符号を出力する関数であり、○は要素ごとの掛け算を表す記号である。ｍａｘ｛｜・｜，０｝は要素ごとに計算されるものとする。この計算は、非特許文献１のSoft thresholding法を応用することで計算することができる。

個別信号推定部１１１４は、式（３３）をｉ＝１，２，…，Ｎｃそれぞれについて求める（Ｓ１１１４）。式（３２）は、式（１４）のうち、推定したい信号ｓ_ｄ，ｉに関連する項のみを取り出したものである。

個別信号推定部１１１４は、パラメータ記憶部１１２に記憶されたｓ_ｄ，ｉの値を、上述の最小化により求めたｓ_ｄ，ｉの値に更新する（Ｓ１１１４）。

＜振幅位相推定部１１１５＞
振幅位相推定部１１１５は、パラメータ記憶部１１２に記憶されたｚ_ｃ，ｉ，ｚ_{ｃ，ｉ＋１}，μ_ｃ，ｉ，ａ_ｃ，ｉ，ｅ_ｃ，ｉ，β，ν_ｃ，ｉ，Λ，γを用いて、 振幅係数ａ_ｃ，ｉと位相係数θ_ｃ，ｉを推定する（Ｓ１１１５）。具体的には、振幅位相推定部１１１５は、

をａ_ｃ，ｉについて最小化することによりａ_ｃ，ｉを求める（Ｓ１１１５）。この最小化は、

を求めることにより行われる。振幅位相推定部１１１５は、これを、ｉ＝１，２，…，Ｎｃ−１それぞれについて行う（Ｓ１１１５）。なお、式（３４）は、式（１４）のうち、推定したい振幅係数ａ_ｃ，ｉに関連する項のみを取り出したものである。
また、振幅位相推定部１１１５は、

をθ_ｃ，ｉについて最小化することにより、位相係数θ_ｃ，ｉを求める（Ｓ１１１５）。この最小化は数値計算法などによって行われる。振幅位相推定部１１１５は、ｉ＝１，２，…，Ｎｃ−１それぞれについて最小化を行う（Ｓ１１１５）。なお、式（３６）は、式（１４）のうち、推定したい位相係数θ_ｃ，ｉに関連する項のみを取り出したものである。

振幅位相推定部１１１５は、パラメータ記憶部１１２に記憶されたａ_ｃ，ｉ及びθ_ｃ，ｉの値を、上述の最小化により求めたａ_ｃ，ｉ及びθ_ｃ，ｉの値にそれぞれ更新する（Ｓ１１１５）。

このように、本実施例の圧縮信号復元装置１によれば、従来法と異なり、パラメータの推定を逐次的にかつペナルティを設けることでオーバーフィッティングを避けるように復元を実行することができるため、元の信号を従来よりも精度良く復元することが出来る。

次に、図７から図１０を参照して実施例２の圧縮信号復元装置２について説明する。図７は本実施例の圧縮信号復元装置２の構成を示すブロック図である。図８は本実施例の圧縮信号復元装置２の動作を示すフローチャートである。図９は本実施例の圧縮信号復元装置２の信号推定部２１１の構成を示すブロック図である。図１０は本実施例の圧縮信号復元装置２の信号推定部２１１の動作を示すフローチャートである。本実施例の圧縮信号復元装置２の信号推定部２１１は、図９に示すように実施例１の圧縮信号復元装置１の信号推定部１１１が備える振幅位相推定部１１１５を省略した構成である。従って、図１０に示すように本実施例の信号推定部２１１の動作（Ｓ２１１）は、ステップＳ１１１５が省略されること以外は、実施例１と同じであるため、説明を割愛する。本実施例においては、ａ_ｃ，ｉおよびθ_ｃ，ｉに予め定めた値を設定しておくことで、ステップＳ１１１５を省略可能である。

このように、本実施例の圧縮信号復元装置２によれば、あらかじめ振幅や位相を考慮する必要がない場合、あるいは振幅や位相が分かっている場合に、振幅位相推定処理（ステップＳ１１１５）を省略して実施例１と同様の効果を得ることができる。

次に、図１１から図１４を参照して実施例３の圧縮信号復元装置３について説明する。図１１は本実施例の圧縮信号復元装置３の構成を示すブロック図である。図１２は本実施例の圧縮信号復元装置３の動作を示すフローチャートである。図１３は本実施例の圧縮信号復元装置３の信号推定部３１１の構成を示すブロック図である。図１４は本実施例の圧縮信号復元装置３の信号推定部３１１の動作を示すフローチャートである。本実施例の圧縮信号復元装置３の信号推定部３１１は、実施例１の信号推定部１１１から共通信号推定部１１１３、振幅位相推定部１１１５を省略した構成である。さらに、本実施例の信号推定部３１１は、実施例１の要素信号推定部１１１２の代わりに要素信号推定部３１１２を備える。また、実施例１のパラメータ記憶部１１２は、本実施例においてパラメータ記憶部３１２に変更されている。また、実施例１の乗数更新部１１３は、本実施例において乗数更新部３１３に変更されている。

前述した実施例２の場合のように、振幅と位相が既知の場合は、ｅ_ｃ，ｉを導入する必要がない。そこで、本実施例においては、信号推定部３１１において、（式１４）の代わりに、

で表される目的関数Ｊを、信号ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｓ_ｄ，ｉのそれぞれについて最小化することにより、信号の値を求める例について説明する。

上記の処理を実行するため、本実施例では下記のように実施例１の各構成部が変更される。

＜パラメータ記憶部３１２＞
実施例１のパラメータ記憶部１１２との違いは、式（１４’）では用いないパラメータγ，ｅ_ｃ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉ，ν_ｃ，ｉが記憶されない点のみである。

＜信号推定部３１１＞
信号推定部３１１は、上記式（１４’）で表される目的関数Ｊを、信号ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｓ_ｄ，ｉのそれぞれについて最小化することにより、信号の値を求める（Ｓ３１１）。そして、パラメータ記憶部３１２に記憶された信号ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｓ_ｄ，ｉの値を、最小化することにより求めた各信号の値に更新する（Ｓ３１１）。

＜全体信号推定部１１１１＞
全体信号推定部１１１１の構成、入出力、機能は実施例１と同様である。

＜要素信号推定部３１１２＞
要素信号推定部３１１２は、パラメータ記憶部３１２に記憶されたｘ_ｉ，ｚ_ｃ，ｉ，ｚ_{ｃ，ｉ＋１}，ｚ_ｄ，ｉ，λ_ｉ，α，μ_ｃ，ｉ，βを用いて、センサ毎の要素信号の推定値ｚ_ｃ，ｉとｚ_ｄ，ｉを求め、パラメータ記憶部３１２に記憶する（Ｓ３１１２）。

具体的には、要素信号推定部３１１２は、ｉ＝１については

をｚ_ｃ，ｉについて最小化し、ｚ_ｃ，ｉを出力する（Ｓ３１１２）。この最小化は、

のように行われる。
また、要素信号推定部３１１２は、ｉ＝２，３，…，Ｎｃ−１については

をｚ_ｃ、ｉについて最小化し、ｚ_ｃ、ｉを出力する（Ｓ３１１２）。この最小化は、

のように行われる。
また、要素信号推定部３１１２は、ｉ＝Ｎｃについては

をｚ_ｃ、ｉについて最小化し、ｚ_ｃ，ｉを出力する（Ｓ３１１２）。この最小化は、

のように行われる。
また、要素信号推定部３１１２はｉ＝１，２，…，Ｎｃについて、上記式（２８）をｚ_ｄ，ｉについて最小化することによりｚ_ｄ，ｉを求め、パラメータ記憶部３１２に記憶されたｚ_ｄ，ｉの値を最小化して得られたｚ_ｄ，ｉの値に更新する（Ｓ３１１２）。この最小化はたとえば、式（２９）のように行われる。

＜個別信号推定部１１１４＞
個別信号推定部１１１４の構成、入出力、機能は実施例１と同様である。

＜乗数更新部３１３＞
乗数更新部３１３は、実施例１の式（１５）〜（１８）の代わりに、下記式により各乗数の値を更新する（Ｓ３１３）。

つまり、本実施例の乗数更新部３１３は、式（１７）は用いず、式（１６）が式（１６’）に置き換わる点を除いては、実施例１の乗数更新部１１３と同じである。

＜終了判定部１１４＞
終了判定部１１４の構成、入出力、機能は実施例１と同様である。

このように、本実施例の圧縮信号復元装置３によれば、あらかじめ振幅や位相を考慮する必要がない場合、あるいは振幅や位相が分かっている場合に、共通信号推定処理（ステップＳ１１１３）、振幅位相推定処理（ステップＳ１１１５）を省略し、式（１４）の代わりに式（１４’）を用いて実施例１と同様の効果を得ることができる。

次に、図１５から図１８を参照して実施例４の圧縮信号復元装置４について説明する。図１５は本実施例の圧縮信号復元装置４の構成を示すブロック図である。図１６は本実施例の圧縮信号復元装置４の動作を示すフローチャートである。図１７は本実施例の圧縮信号復元装置４の信号推定部４１１の構成を示すブロック図である。図１８は本実施例の圧縮信号復元装置４の信号推定部４１１の動作を示すフローチャートである。実施例１の信号推定部１１１、パラメータ記憶部１１２、乗数更新部１１３は、本実施例においてそれぞれ信号推定部４１１、パラメータ記憶部４１２、乗数更新部４１３に変更されている。実施例１の全体信号推定部１１１１、要素信号推定部１１１２、共通信号推定部１１１３、振幅位相推定部１１１５は、本実施例においてそれぞれ全体信号推定部４１１１、要素信号推定部４１１２、共通信号推定部４１１３、振幅位相推定部４１１５に変更されている。

実施例１〜３では共通信号は１つ（ｚ_ｃ）であることを仮定していたが、共通信号が複数である場合にも実施例１を拡張することができる。共通信号が複数の場合とは、例えば各センサが複数の観測対象（ソース）を観測して信号を取得する場合等が考えられる。
すなわち、ｊ＝１，２，…，Ｎｓ（Ｎｓは共通信号の総数）として、第ｉ番目のセンサの信号ｘ_ｉが、

のように記述できる場合である。Ωｓは共通信号のインデクスの集合Ωｓ＝｛１，２，・・・，Ｎｓ｝である。ここでｚ_{ｃ，ｉ，ｋ}は第ｉ番目のセンサにおけるｋ番目の共通信号である。それぞれの共通信号は

のような関係で書くことができると仮定する。ここでａ_ｉ，ｋを共通信号ｚ_ｃ，ｋに対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の振幅を表す係数とし、θ_ｃ，ｋをｚ_ｃ，ｋに対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の位相を表す係数とする。ａ_ｉ，ｋ、θ_ｃ，ｋはそれぞれ実数のスカラー量である。

上記の処理を実行するため、本実施例では下記のように実施例１の各構成部が変更される。

＜パラメータ記憶部４１２＞
パラメータ記憶部４１２には、予め定められた圧縮行列Ａ_ｉ、観測行列Ψ_ｉ、基底行列Φを記憶しておく。また、復号部４１の処理の過程で、信号の推定値を得るために必要なパラメータα，β，γ，σ，ｗと、信号および係数ｘ_ｉ，ｚ_{ｃ，ｉ，ｋ}，ｚ_ｄ，ｉ，ｅ_{ｃ，ｉ，ｋ}，ｓ_ｄ，ｉ，ａ_{ｃ，ｉ，ｋ}，θ_{ｃ，ｉ，ｋ}，λ_ｉ，μ_{ｃ，ｉ，ｋ}，ν_{ｃ，ｉ，ｋ}，ξ_ｄ，ｉの初期値をあらかじめ記憶しておく。なお、初期値は予めランダムな値を設定しておけば良い。なお、ａ_{ｃ，ｉ，ｋ}を補助共通信号ｅ_{ｃ，ｉ，ｋ}に対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の振幅を表す係数とし、θ_{ｃ，ｉ，ｋ}をｚ_{ｃ，ｉ＋１，ｋ}に対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の位相を表す係数とする。

＜信号推定部４１１＞
信号推定部４１１は、上記式（１４）の代わりに下記式（４０）を、信号および係数ｘ_ｉ，ｚ_{ｃ，ｉ，ｋ}，ｚ_ｄ，ｉ，ｅ_{ｃ，ｉ，ｋ}，ｓ_ｄ，ｉ，ａ_{ｃ，ｉ，ｋ}，θ_{ｃ，ｉ，ｋ}のそれぞれについて最小化することにより、信号および係数の値を求める（Ｓ４１１）。そして、パラメータ記憶部４１２に記憶された信号および係数ｘ_ｉ，ｚ_{ｃ，ｉ，ｋ}，ｚ_ｄ，ｉ，ｅ_{ｃ，ｉ，ｋ}，ｓ_ｄ，ｉ，ａ_{ｃ，ｉ，ｋ}，θ_{ｃ，ｉ，ｋ}の値を、最小化することにより求めた各信号および係数の値に更新する。

＜全体信号推定部４１１１＞
全体信号推定部４１１１は入力された信号ｙ_ｉと、パラメータ記憶部４１２に記憶されたｚ_{ｃ，ｉ，ｋ}，ｚ_ｄ，ｉ，λ_ｉ，α，Ψ_ｉ，Ｎｃ，Ｎｓを入力として受け取り、ｙ_ｉに対応する復号信号の推定値ｘ_ｉを求める（Ｓ４１１１）。具体的には、全体信号推定部４１１１は、

をｘ_ｉについて最小化することによりｘ_ｉを求め、パラメータ記憶部４１２に記憶する（Ｓ４１１１）。この最小化はたとえば、

のように実行される。全体信号推定部４１１１は、式（４２）を、ｉ＝１，２，…，Ｎｃそれぞれについて計算することにより、ｘ_ｉを求める（Ｓ４１１１）。

＜要素信号推定部４１１２＞
要素信号推定部４１１２は、パラメータ記憶部４１２に記憶されたｘ_ｉ，ｚ_{ｃ，ｉ，ｋ}，ｚ_ｄ，ｉ，λ_ｉ，α，μ_{ｃ，ｉ，ｋ}，ａ_{ｃ，ｉ，ｋ}，ｅ_{ｃ，ｉ，ｋ}，βを用いて、センサ毎の要素信号の推定値ｚ_{ｃ，ｉ，ｋ}とｚ_ｄ，ｉを求め、パラメータ記憶部４１２に記憶する（Ｓ４１１２）。
具体的には、要素信号推定部４１１２は、ｉ＝１については

をｚ_{ｃ，ｉ，ｋ}について最小化することによりｚ_{ｃ，ｉ，ｋ}を求める（Ｓ４１１２）。この最小化は、

のように行われる。ここで、Ωｓ＼ｋは、Ωｓからｋを除いた集合を表す。
また、要素信号推定部４１１２は、ｉ＝２，３，…，Ｎｃ−１については

をｚ_{ｃ、ｉ，ｋ}について最小化することによりｚ_{ｃ、ｉ，ｋ}を求める（Ｓ４１１２）。この最小化は、

のように行われる。
また、要素信号推定部４１１２は、ｉ＝Ｎｃについては

をｚ_{ｃ、ｉ，ｋ}について最小化することによりｚ_{ｃ，ｉ，ｋ}を求める（Ｓ４１１２）。この最小化は、

のように行われる。

また、要素信号推定部４１１２はｉ＝１，２，…，Ｎｃについて、

をｚ_ｄ，ｉについて最小化し、ｚ_ｄ，ｉを出力する（Ｓ４１１２）。この最小化はたとえば、

のように行われる。

要素信号推定部４１１２は、パラメータ記憶部４１２に記憶されたｚ_{ｃ，ｉ，ｋ}及びｚ_ｄ，ｉの値を、上述の最小化により求めたｚ_{ｃ，ｉ，ｋ}及びｚ_ｄ，ｉの値にそれぞれ更新する（Ｓ４１１２）。

＜共通信号推定部４１１３＞
共通信号推定部４１１３は、パラメータ記憶部４１２に記憶されたｚ_{ｃ，ｉ，ｋ}，θ_{ｃ，ｉ，ｋ}，μ_{ｃ，ｉ，ｋ}，ａ_{ｃ，ｉ，ｋ}，ｅ_{ｃ，ｉ，ｋ}，β，ν_{ｃ，ｉ，ｋ}，Λ，ｚ_{ｃ，ｉ＋１，ｋ}，γを用いて、補助共通信号ｅ_{ｃ，ｉ，ｋ}を推定する（Ｓ４１１３）。具体的には、共通信号推定部４１１３は、

をｅ_{ｃ，ｉ，ｋ}について最小化することにより、ｅ_{ｃ，ｉ，ｋ}を求める（Ｓ４１１３）。この最小化は、

を求めることにより行われる。共通信号推定部４１１３は、これを、ｉ＝１，２，…，Ｎｃ−１それぞれについて求める（Ｓ４１１３）。共通信号推定部４１１３は、パラメータ記憶部４１２に記憶されたｅ_{ｃ，ｉ，ｋ}の値を、上述の最小化により求めたｅ_{ｃ，ｉ，ｋ}の値に更新する（Ｓ４１１３）。

＜個別信号推定部１１１４＞
個別信号推定部１１１４の構成、入出力、機能は実施例１と同様である。

＜振幅位相推定部４１１５＞
振幅位相推定部４１１５は、パラメータ記憶部４１２に記憶されたｚ_{ｃ，ｉ，ｋ}，θ_{ｃ，ｉ，ｋ}，μ_{ｃ，ｉ，ｋ}，ａ_{ｃ，ｉ，ｋ}，ｅ_{ｃ，ｉ，ｋ}，β，ν_{ｃ，ｉ，ｋ}，Λ，ｚ_{ｃ，ｉ＋１，ｋ}，γを用いて、振幅ａ_{ｃ，ｉ，ｋ}と位相係数θ_{ｃ，ｉ，ｋ}を推定する（Ｓ４１１５）。具体的には、振幅位相推定部４１１５は、

をａ_{ｃ，ｉ，ｋ}について最小化することによりａ_{ｃ，ｉ，ｋ}を求める（Ｓ４１１５）。この最小化は、

を求めることにより行われる。振幅位相推定部４１１５は、これを、ｉ＝１，２，…，Ｎｃ−１それぞれについて行う（Ｓ４１１５）。
また、振幅位相推定部４１１５は、

をθ_{ｃ，ｉ，ｋ}について最小化することにより、位相係数θ_{ｃ，ｉ，ｋ}を求める（Ｓ４１１５）。この最小化は数値計算法などによって行われる。振幅位相推定部４１１５は、ｉ＝１，２，…，Ｎｃ−１それぞれについて最小化を行う（Ｓ４１１５）。

振幅位相推定部４１１５は、パラメータ記憶部４１２に記憶されたａ_{ｃ，ｉ，ｋ}及びθ_{ｃ，ｉ，ｋ}の値を、上述の最小化により求めたａ_{ｃ，ｉ，ｋ}及びθ_{ｃ，ｉ，ｋ}の値にそれぞれ更新する（Ｓ４１１５）。

＜乗数更新部４１３＞
乗数更新部４１３は、信号推定部４１１で推定したλ_ｉ，α，ｘ_ｉ，ｚ_{ｃ，ｉ，ｋ}，ｚ_ｄ，ｉ，μ_{ｃ，ｉ，ｋ}，β，ν_{ｃ，ｉ，ｋ}，γ，ｅ_{ｃ，ｉ，ｋ}，ｅ_{ｃ，ｉ−１，ｋ}，ｓ_ｄ，ｉ，Λ，ａ_{ｃ，ｉ，ｋ}，θ_{ｃ，ｉ−１，ｋ}，ξ_ｄ，ｉ，σ，Φを用いて、式（４０）の各制約条件に対応する未定乗数λ_ｉ，μ_{ｃ，ｉ，ｋ}，ν_{ｃ，ｉ，ｋ}，ξ_ｄ，ｉを更新する（Ｓ４１３）。具体的には、乗数更新部４１３は、下記式により各乗数の値を更新する（Ｓ４１３）。

＜終了判定部１１４＞
終了判定部１１４の構成、入出力、機能は実施例１と同様である。

このように、本実施例の圧縮信号復元装置４によれば、共通信号が複数ある場合にも実施例１と同様の効果を奏する。

次に、図１９から図２２を参照して実施例５の圧縮信号復元装置５について説明する。図１９は本実施例の圧縮信号復元装置５の構成を示すブロック図である。図２０は本実施例の圧縮信号復元装置５の動作を示すフローチャートである。図２１は本実施例の圧縮信号復元装置５の信号推定部５１１の構成を示すブロック図である。図２２は本実施例の圧縮信号復元装置５の信号推定部５１１の動作を示すフローチャートである。実施例１の信号推定部１１１、乗数更新部１１３は、本実施例においてそれぞれ信号推定部５１１、乗数更新部５１３に変更されている。また、本実施例の圧縮信号復元装置５は実施例１の圧縮信号復元装置１が備えないインデクス入れ替え部５１５、対応テーブル記憶部５１６を備える。実施例１の要素信号推定部１１１２、共通信号推定部１１１３、振幅位相推定部１１１５は、本実施例においてそれぞれ要素信号推定部５１１２、共通信号推定部５１１３、振幅位相推定部５１１５に変更されている。

実施例４においては複数の共通信号がある場合を仮定していたが、場合によっては共通信号間の区別をつけることが困難である（あいまいさがある）ために最適化が上手くいかない場合がある。実施例５ではこれを回避するために、共通信号のパワーが大きい順になるように共通信号のインデクスを入れ替えることを行う。

上記の処理を実行するため、本実施例では下記のように実施例１の各構成部が変更される。

＜対応テーブル記憶部５１６＞
対応テーブル記憶部５１６は、変換前の共通信号のインデクスと、変換後（パワーが大きい順に並べ替えた後の）インデクスとの対応関係を表す対応テーブルを記憶する。対応テーブルの初期値は、変換前のインデクスと変換後のインデクスが同じ（変化しない）状態としておく。

＜インデクス入れ替え部５１５＞
インデクス入れ替え部５１５は、要素信号推定部５１１２からｚ_{ｃ，ｉ，ｋ}を受け取り、予め定めたいずれかの第ｉ番目のセンサに着目して、各ｋ＝１，２，…，Ｎｓについてｚ_{ｃ，ｉ，ｋ}のＬ２ノルムを計算する（Ｓ５１５）。そして、インデクス入れ替え部５１５は、変換後のインデクスｋ’＝１，２，・・・，Ｎｓが、ｚ_{ｃ，ｉ，ｋ’}のＬ２ノルムの大きい順になるように、並べ替え前のインデクスｋと変換後のインデクスｋ’との対応関係を求め、対応テーブル記憶部５１６の対応テーブルの内容を更新する（Ｓ５１５）。

要素信号推定部５１１２、共通信号推定部５１１３、振幅位相推定部５１１５、乗数更新部５１３は、基本的には実施例４と同じ動作をするが、共通信号のインデクスｋについては、対応テーブル記憶部５１６に記憶された対応テーブルを参照して、変換後のインデクスの順に各変数を入れ替えてから計算を行う。他の部の構成、入出力、機能については実施例１と同様である。

このように、本実施例の圧縮信号復元装置５によれば、共通信号間の区別をつけることが困難な場合であっても、共通信号のパワーが大きい順になるように共通信号のインデクスを入れ替えることで、最適化を実行できる。

実施例１〜４については個別信号のスパース性を仮定していた。本実施例ではこの制約をなくす。具体的には、信号推定部６１１として、

または

を目的関数として最小化する構成とすることができる。この場合は各実施例から個別信号推定部１１１４を取り除くとともに、上記式（６０）〜（６２）にそれぞれ対応した各変数について最小化を行えば良い。

実施例１〜４において重みｗ_ｉは固定された値であったが、ＦＰＣ（Ｆｉｘｅｄ−Ｐｏｉｎｔ Ｃｏｎｔｉｎｕａｔｉｏｎ，非特許文献１）のような方法でｗ_ｉを更新してもよい。乗数更新部７１３において、さらに重み係数ｗ_ｉの更新処理を行う。

＜乗数更新部７１３＞
乗数更新部７１３は、実施例１〜４のいずれかに記載の乗数更新部１１３〜４１３の処理に加え、下記の更新式に基づいて、パラメータ記憶部に記憶された重みｗ_ｉを更新する（Ｓ７１３）。

ここで、ηは予め定めた定数である。ただしｗ_ｉが予め定めた下限値ｗ_ｍｉｎよりも小さい場合はｗ_ｍｉｎを出力する構成としても良い。

なお、上述の重みの更新処理を行うか否かを

により判定しても良い。つまり、式（６４）が真ならば式（６３）による重みの更新処理を行い、式（６４）が偽ならは重みの更新処理は行わない構成としても良い。ここで、ε_ｄ，ｉは予め定めた定数である。

＜終了判定部７１４＞
終了判定部７１４は、実施例１の終了判定部１１４における「所定の終了条件」として、式（１９）に加えて、下記式

を用いる。すなわち、「信号推定部の繰り返し回数ｋが予め定めた回数Ｉｔ_ｍａｘに到達し」かつ「更新後の重みｗ_ｉがｗ_ｍｉｎより小さい」という条件を充たした場合に、「所定の条件を充たした」と判定する。

また、上述の各種の処理は、記載に従って時系列に実行されるのみならず、処理を実行する装置の処理能力あるいは必要に応じて並列的にあるいは個別に実行されてもよい。その他、本発明の趣旨を逸脱しない範囲で適宜変更が可能であることはいうまでもない。

また、上述の構成をコンピュータによって実現する場合、各装置が有すべき機能の処理内容はプログラムによって記述される。そして、このプログラムをコンピュータで実行することにより、上記処理機能がコンピュータ上で実現される。

この処理内容を記述したプログラムは、コンピュータで読み取り可能な記録媒体に記録しておくことができる。コンピュータで読み取り可能な記録媒体としては、例えば、磁気記録装置、光ディスク、光磁気記録媒体、半導体メモリ等どのようなものでもよい。

また、このプログラムの流通は、例えば、そのプログラムを記録したＤＶＤ、ＣＤ−ＲＯＭ等の可搬型記録媒体を販売、譲渡、貸与等することによって行う。さらに、このプログラムをサーバコンピュータの記憶装置に格納しておき、ネットワークを介して、サーバコンピュータから他のコンピュータにそのプログラムを転送することにより、このプログラムを流通させる構成としてもよい。

このようなプログラムを実行するコンピュータは、例えば、まず、可搬型記録媒体に記録されたプログラムもしくはサーバコンピュータから転送されたプログラムを、一旦、自己の記憶装置に格納する。そして、処理の実行時、このコンピュータは、自己の記録媒体に格納されたプログラムを読み取り、読み取ったプログラムに従った処理を実行する。また、このプログラムの別の実行形態として、コンピュータが可搬型記録媒体から直接プログラムを読み取り、そのプログラムに従った処理を実行することとしてもよく、さらに、このコンピュータにサーバコンピュータからプログラムが転送されるたびに、逐次、受け取ったプログラムに従った処理を実行することとしてもよい。また、サーバコンピュータから、このコンピュータへのプログラムの転送は行わず、その実行指示と結果取得のみによって処理機能を実現する、いわゆるＡＳＰ（Application Service Provider）型のサービスによって、上述の処理を実行する構成としてもよい。

なお、本形態におけるプログラムには、電子計算機による処理の用に供する情報であってプログラムに準ずるもの（コンピュータに対する直接の指令ではないがコンピュータの処理を規定する性質を有するデータ等）を含むものとする。また、この形態では、コンピュータ上で所定のプログラムを実行させることにより、本装置を構成することとしたが、これらの処理内容の少なくとも一部をハードウェア的に実現することとしてもよい。

Claims (8)

1. Ｎｃ個のセンサからの圧縮信号を復元する圧縮信号復元装置であって、
   インデクスｉを１以上Ｎｃ以下の整数とし、全てのセンサに共通の信号を信号ｚ_ｃとし、第ｉ番目のセンサが観測した信号の共通部分を共通信号ｚ_ｃ，ｉと表すものとし、第ｉ番目のセンサが観測した信号のうち共通信号ｚ_ｃ，ｉを除いた残りの信号を信号ｚ_ｄ，ｉと表すものとし、
   

   

   
   で表される第ｉ番目のセンサの観測信号ｘ_ｉを所定の観測行列Ψ_ｉにより圧縮した信号ｙ_ｉ＝Ψ_ｉｘ_ｉを入力とし、
   ｅ_ｃ，ｉを第ｉ番目のセンサの補助共通信号とし、ａ_ｃ，ｉを補助共通信号ｅ_ｃ，ｉに対する第ｉ番目のセンサの共通信号要素の振幅を表す係数とし、θ_ｃ，ｉをｚ_{ｃ，ｉ＋１}に対する第ｉ番目のセンサの共通信号要素の位相を表す係数とし、Φを予め定められた基底行列とし、共通信号ｚ_ｃの基底Φでの表現ｓ_ｃ＝Φ^−１ｚ_ｃが疎ではなく、信号ｚ_ｄ，ｉの基底Φでの表現ｓ_ｄ，ｉ＝Φ^−１ｚ_ｄ，ｉが全てのセンサのインデクスｉに対して疎であるものとし、‖・‖_２をベクトルの各要素の二乗和であるＬ２ノルムとし、‖・‖_１をベクトルの各要素の絶対値和であるＬ１ノルムとし、Ωをインデクスｉの集合、すなわちΩ＝｛１，２，・・・，Ｎｃ｝とし、ｗ_ｉを、予め定めた個別信号ｓ_ｄ，ｉに対する重みとし、∀を全称記号とし、ｊを虚数単位、ＤＦＴを離散フーリエ変換、ＩＤＦＴをＤＦＴの逆変換とし、記号Λを
   

   

   
   と定義し、
   

   

   
   の制約条件のもとで、
   

   

   
   が最小となる観測信号ｘ_ｉを求め、元の信号の復号信号として出力する復号部
   を含む圧縮信号復元装置。
2. 請求項１に記載の圧縮信号復元装置であって、
   λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉをそれぞれ前記制約条件に対する未定乗数であるものとし、α，β，γ，σを、最小化問題を解く際の収束速度を調整するための重みとし、Ｔを行列またはベクトルの転置を表すものとし、
   前記復号部は、
   

   

   
   で表される目的関数Ｈを、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉのそれぞれについて最小化することにより、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉを推定する信号推定部と、
   前記推定された信号および係数に基づいて、前記未定乗数λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉを更新する乗数更新部と、
   所定の条件を充たすまで、前記信号推定部と、前記乗数更新部の処理を繰り返し制御する終了判定部と、
   を含む圧縮信号復元装置。
3. 請求項２に記載の圧縮信号復元装置であって、
   ＊をエルミート転置行列を表す記号とし、Ｉを大きさＮの単位行列とし
   前記信号推定部が、
   

   

   
   により信号ｘ_ｉを求める全体信号推定部と、
   ｉ＝１について、
   

   

   
   ｉ＝２，３，…，Ｎｃ−１について、
   

   

   
   ｉ＝Ｎｃについて、
   

   

   
   により、信号ｚ_ｃ，ｉを求め、
   ｉ＝１，２，…，Ｎｃについて、
   

   

   
   により信号ｚ_ｄ，ｉを求める要素信号推定部と、
   ｉ＝１，２，…，Ｎｃ−１について、
   

   

   
   により信号ｅ_ｃ，ｉを求める共通信号推定部と、
   ｓｇｎを要素ごとの符号を出力する関数とし、○を要素ごとの掛け算を表す記号とし、ｍａｘ｛｜・｜，０｝は要素ごとに計算されるものとし、
   ｉ＝１，２，…，Ｎｃについて、
   

   

   
   により信号ｓ_ｄ，ｉを求める個別信号推定部と、
   ｉ＝１，２，…，Ｎｃ−１について、
   

   

   
   により係数ａ_ｃ，ｉを求め、
   

   

   
   をθ_ｃ，ｉについて最小化することにより係数θ_ｃ，ｉを求める振幅位相推定部と、
   を含む圧縮信号復元装置。
4. 請求項１に記載の圧縮信号復元装置であって、
   λ_ｉ，μ_ｃ，ｉ，ξ_ｄ，ｉをそれぞれ前記制約条件に対する未定乗数であるものとし、α，β，σを、最小化問題を解く際の収束速度を調整するための重みとし、Ｔを行列またはベクトルの転置を表すものとし、
   前記復号部は、
   

   

   
   で表される目的関数Ｊを、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｓ_ｄ，ｉのそれぞれについて最小化することにより、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｓ_ｄ，ｉを推定する信号推定部と、
   前記推定された信号および係数に基づいて、前記未定乗数λ_ｉ，μ_ｃ，ｉ，ξ_ｄ，ｉを更新する乗数更新部と、
   所定の条件を充たすまで、前記信号推定部と、前記乗数更新部の処理を繰り返し制御する終了判定部と、
   を含む圧縮信号復元装置。
5. Ｎｃ個のセンサからの圧縮信号を復元する圧縮信号復元装置であって、
   インデクスｉを１以上Ｎｃ以下の整数とし、全てのセンサに共通の信号を信号ｚ_ｃとし、第ｉ番目のセンサが観測した信号の共通部分を共通信号ｚ_ｃ，ｉと表すものとし、第ｉ番目のセンサにおけるｋ番目の共通信号をｚ_{ｃ，ｉ，ｋ}と表すものとし、第ｉ番目のセンサが観測した信号のうち共通信号ｚ_ｃ，ｉを除いた残りの信号を信号ｚ_ｄ，ｉと表すものとし、Ｎｓを共通信号の総数とし、Ωｓを共通信号のインデクスの集合、すなわちΩｓ＝｛１，２，・・・，Ｎｓ｝とし、ｊ＝１，２，…，Ｎｓとし、ｊを虚数単位、ＤＦＴを離散フーリエ変換、ＩＤＦＴをＤＦＴの逆変換とし、記号Λを
   

   

   
   と定義し、
   

   

   
   で表される第ｉ番目のセンサの観測信号ｘ_ｉを所定の観測行列Ψ_ｉにより圧縮した信号ｙ_ｉ＝Ψ_ｉｘ_ｉを入力とし、ａ_ｉ，ｋを共通信号ｚ_ｃ，ｋに対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の振幅を表す実数のスカラー量である係数とし、θ_ｃ，ｋをｚ_ｃ，ｋに対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の位相を表す実数のスカラー量である係数とし、 λ_ｉ，μ_{ｃ，ｉ，ｋ}，ν_{ｃ，ｉ，ｋ}，ξ_ｄ，ｉをそれぞれ前記制約条件に対する未定乗数であるものとし、α，β，γ，σを、最小化問題を解く際の収束速度を調整するための重みとし、Ｔは行列またはベクトルの転置を表すものとし、
   ｅ_{ｃ，ｉ，ｋ}を第ｉ番目のセンサにおけるｋ番目の補助共通信号とし、ａ_{ｃ，ｉ，ｋ}を補助共通信号ｅ_{ｃ，ｉ，ｋ}に対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の振幅を表す係数とし、θ_{ｃ，ｉ，ｋ}をｚ_{ｃ，ｉ＋１，ｋ}に対する第ｉ番目のセンサにおけるｋ番目の共通信号要素の位相を表す係数とし、Φを予め定められた基底行列とし、共通信号ｚ_ｃの基底Φでの表現ｓ_ｃ＝Φ^−１ｚ_ｃが疎ではなく、信号ｚ_ｄ，ｉの基底Φでの表現ｓ_ｄ，ｉ＝Φ^−１ｚ_ｄ，ｉが全てのセンサのインデクスｉに対して疎であるものとし、‖・‖_２をベクトルの各要素の二乗和であるＬ２ノルムとし、‖・‖_１をベクトルの各要素の絶対値和であるＬ１ノルムとし、Ωをセンサのインデクスｉの集合、すなわちΩ＝｛１，２，・・・，Ｎｃ｝とし、Ω⁻をＮｃを除くインデクスｉの集合、すなわちΩ⁻＝｛１，２，・・・，Ｎｃ−１｝とし、ｗ_ｉを予め定めた個別信号ｓ_ｄ，ｉに対する重みとし、
   

   

   
   で表される目的関数Ｈ_２を、信号および係数ｘ_ｉ，ｚ_{ｃ，ｉ、ｋ}，ｚ_ｄ，ｉ，ｓ_ｄ，ｉのそれぞれについて最小化することにより、元の信号の復号信号として出力する復号部
   を含む圧縮信号復元装置。
6. Ｎｃ個のセンサからの圧縮信号を復元する圧縮信号復元方法であって、
   インデクスｉを１以上Ｎｃ以下の整数とし、全てのセンサに共通の信号を信号ｚ_ｃとし、第ｉ番目のセンサが観測した信号の共通部分を共通信号ｚ_ｃ，ｉと表すものとし、第ｉ番目のセンサが観測した信号のうち共通信号ｚ_ｃ，ｉを除いた残りの信号を信号ｚ_ｄ，ｉと表すものとし、
   

   

   
   で表される第ｉ番目のセンサの観測信号ｘ_ｉを所定の観測行列Ψ_ｉにより圧縮した信号ｙ_ｉ＝Ψ_ｉｘ_ｉを入力とし、
   ｅ_ｃ，ｉを第ｉ番目のセンサの補助共通信号とし、ａ_ｃ，ｉを補助共通信号ｅ_ｃ，ｉに対する第ｉ番目のセンサの共通信号要素の振幅を表す係数とし、θ_ｃ，ｉをｚ_{ｃ，ｉ＋１}に対する第ｉ番目のセンサの共通信号要素の位相を表す係数とし、Φを予め定められた基底行列とし、共通信号ｚ_ｃの基底Φでの表現ｓ_ｃ＝Φ^−１ｚ_ｃが疎ではなく、信号ｚ_ｄ，ｉの基底Φでの表現ｓ_ｄ，ｉ＝Φ^−１ｚ_ｄ，ｉが全てのセンサのインデクスｉに対して疎であるものとし、‖・‖_２をベクトルの各要素の二乗和であるＬ２ノルムとし、‖・‖_１をベクトルの各要素の絶対値和であるＬ１ノルムとし、Ωをインデクスｉの集合、すなわちΩ＝｛１，２，・・・，Ｎｃ｝とし、ｗ_ｉを、予め定めた個別信号ｓ_ｄ，ｉに対する重みとし、∀を全称記号とし、ｊを虚数単位、ＤＦＴを離散フーリエ変換、ＩＤＦＴをＤＦＴの逆変換とし、記号Λを
   

   

   
   と定義し、
   

   

   
   の制約条件のもとで、
   

   

   
   が最小となる観測信号ｘ_ｉを求め、元の信号の復号信号として出力する復号ステップ
   を含む圧縮信号復元方法。
7. 請求項６に記載の圧縮信号復元方法であって、
   λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉをそれぞれ前記制約条件に対する未定乗数であるものとし、α，β，γ，σを、最小化問題を解く際の収束速度を調整するための重みとし、Ｔを行列またはベクトルの転置を表すものとし、
   前記復号ステップは、
   

   

   
   で表される目的関数Ｈを、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉのそれぞれについて最小化することにより、信号および係数ｘ_ｉ，ｚ_ｃ，ｉ，ｚ_ｄ，ｉ，ｅ_ｃ，ｉ，ｓ_ｄ，ｉ，ａ_ｃ，ｉ，θ_ｃ，ｉを推定する信号推定サブステップと、
   前記推定された信号および係数に基づいて、前記未定乗数λ_ｉ，μ_ｃ，ｉ，ν_ｃ，ｉ，ξ_ｄ，ｉを更新する乗数更新サブステップと、
   所定の条件を充たすまで、前記信号推定サブステップと、前記乗数更新サブステップの処理を繰り返し制御する終了判定サブステップと、
   を含む圧縮信号復元方法。
8. 請求項６又は７に記載の圧縮信号復元方法を実行すべき指令をコンピュータに対してするプログラム。

Images