"第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換手段と、所定周波数の前記位相成分を回転させる第1制御手段と、前記第1制御手段が施した回転による前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御手段と、前記第1制御手段により回転された位相成分と、前記第2制御手段により補正された振幅成分又はパワー成分と、を合成する合成手段と、を含む信号処理装置。前記第1制御手段は、前記位相成分の回転量を生成する回転量生成手段を含み、前記第2制御手段は、前記回転量生成手段から提供された前記回転量に基づいて前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正する、請求項1に記載の信号処理装置。前記第2制御手段は、前記位相成分を遅延させる位相成分遅延手段を含み、遅延前後の前記位相成分を比較することにより、前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正する、請求項1に記載の信号処理装置。前記第2制御手段は、前記合成手段からの出力と前記混在信号とを比較することにより前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正する、請求項1に記載の信号処理装置。前記第2制御手段は、前記合成手段からの出力の平均値と前記混在信号の平均値とを比較することにより前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正する、請求項4に記載の信号処理装置。前記振幅成分又はパワー成分を遅延させる振幅成分遅延手段と、をさらに含み、前記第2制御手段は、前記合成手段からの出力と前記混在信号とを比較することにより前記変化量を算出し、算出した前記変化量に応じて、前記振幅成分遅延手段により遅延した前記振幅成分又はパワー成分を補正する、請求項1に記載の信号処理装置。前記変換手段は、前記混在信号を、互いに一部オーバーラップした複数フレームに分割する分割手段を含み、前記第2制御手段は、さらに、フレームのオーバーラップ率に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する、請求項1乃至6のいずれか1項に記載の信号処理装置。前記第1制御手段は、前記位相成分の回転量を制限する回転量制限手段を含む、請求項1乃至7のいずれか1項に記載の信号処理装置。第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換ステップと、所定周波数の前記位相成分を回転させる第1制御ステップと、前記第1制御ステップにおける回転によって生じた前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御ステップと、前記第1制御ステップにより回転された位相成分と、前記第2制御ステップにより補正された振幅成分又はパワー成分と、を合成する合成ステップと、を含む信号処理方法。第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換ステップと、所定周波数の前記位相成分を回転させる第1制御ステップと、前記第1制御ステップにおける回転によって生じた前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御ステップと、前記第1制御ステップにより回転された位相成分と、前記第2制御ステップにより補正された振幅成分又はパワー成分と、を合成する合成ステップと、をコンピュータに実行させる信号処理プログラムを格納したプログラム記録媒体。しかしながら、特許文献1や非特許文献1に記載の技術のように、位相スペクトルをランダムにシフトさせると、隣接するフレーム間で位相スペクトルの不一致が生じ、フレーム合成時に信号レベルが低下する。その結果、雑音抑圧区間と雑音非抑圧区間で出力信号レベルに差が生じ、出力信号の品質が低下する。本発明は、上述の課題を解決する信号処理技術を提供することを目的とする。上記目的を達成するため、本発明に係る装置は、第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換手段と、所定周波数の前記位相成分を回転させる第1制御手段と、前記第1制御手段が施した回転による前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御手段と、前記第1制御手段により回転された位相成分と、前記第2制御手段により補正された振幅成分又はパワー成分と、を合成する合成手段と、を含む。上記目的を達成するため、本発明に係る方法は、第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換ステップと所定周波数の前記位相成分を回転させる第1制御ステップと、前記第1制御ステップにおける回転によって生じた前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御ステップと、前記第1制御ステップにより回転された位相成分と、前記第2制御ステップにより補正された振幅成分又はパワー成分と、を合成する合成ステップと、を含む。上記目的を達成するため、本発明に係るプログラム記録媒体に格納されるプログラムは、第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換ステップと所定周波数の前記位相成分を回転させる第1制御ステップと、前記第1制御ステップにおける回転によって生じた前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御ステップと、前記第1制御ステップにより回転された位相成分と、前記第2制御ステップにより補正された振幅成分又はパワー成分と、を合成する合成ステップと、をコンピュータに実行させる。本発明によれば、出力信号の品質を低下させずに入力信号の位相成分を制御する信号処理技術を提供することができる。以下に、図面を参照して、本発明の実施の形態について例示的に詳しく説明する。ただし、以下の実施の形態に記載されている構成要素はあくまで例示であり、本発明の技術範囲をそれらのみに限定する趣旨のものではない。第1実施形態図1は、本発明の第1実施形態に係る信号処理装置100の概略構成を示す図である。図1において、信号処理装置100は、変換部101と第1制御部102と第2制御部103と合成部104とを含む。変換部101は、第1信号と第2信号とが混在した混在信号110を、周波数ごとの位相成分120及び振幅成分又はパワー成分130に変換する。また、第1制御部102は、所定周波数の位相成分120を回転させる。さらに、第2制御部103は、第1制御部102が施した回転による振幅成分又はパワー成分の変化量に応じて、所定周波数の振幅成分又はパワー成分を補正する。一方、合成部104は、第1制御部102により回転された位相成分140と、第2制御部103により補正された振幅成分又はパワー成分150と、を合成する。以上の構成により、第1実施形態に係る信号処理装置100は、出力信号の品質を低下させずに入力信号の位相成分を制御することができる。第2実施形態《全体構成》本発明の第2実施形態としての雑音抑圧装置200について図2乃至図11を用いて説明する。図2は、雑音抑圧装置200の全体構成を示すブロック図である。本実施形態の雑音抑圧装置200は、たとえばデジタルカメラ、ノートパソコン、携帯電話などといった装置の一部としても機能するが、本発明はこれに限定されるものではない。雑音抑圧装置200は、入力信号からのノイズ除去を要求されるあらゆる情報処理装置に適用可能である。入力端子206には、劣化信号所望信号と雑音の混在する信号が、サンプル値系列として供給される。入力端子206に劣化信号が供給されると、変換部201は、供給された劣化信号にフーリエ変換などの変換を施して、複数の周波数成分に分割する。変換部201は、複数の周波数成分を各周波数で独立に処理する。ここでは、特定の周波数成分に注目して説明を続ける。変換部201は、複数の周波数成分のうち振幅スペクトル振幅成分230を雑音抑圧部205に供給する。変換部201は、複数の周波数成分のうち位相スペクトル位相成分220を位相制御部202に供給する。なお、ここでは、変換部201は、雑音抑圧部205に劣化信号振幅スペクトル230を供給しているが、本発明はこれに限定されるものではない。変換部201は、劣化信号振幅スペクトル230の二乗に相当するパワースペクトルを雑音抑圧部205に供給しても良い。雑音抑圧部205は、変換部201から供給される劣化信号振幅スペクトル230を用いて、雑音を推定し、推定雑音スペクトルを生成する。また、雑音抑圧部205は、変換部201から供給された劣化信号振幅スペクトル230と、生成した推定雑音スペクトルとを用いて雑音を抑圧する。雑音抑圧部205は、雑音抑圧結果としての強調信号振幅スペクトルを振幅制御部203に伝達する。位相制御部202は、変換部201から供給された劣化信号位相スペクトル220を回転シフトさせ、強調信号位相スペクトル240として逆変換部204へ供給する。また、位相制御部202は、位相の回転量シフト量を、振幅制御部203へ伝達する。振幅制御部203は、位相制御部202から位相の回転量シフト量を受け取って振幅補正量を算出する。振幅制御部203は、その振幅補正量を用いて、強調信号振幅スペクトルを各周波数で補正し、補正振幅スペクトル250を生成する。振幅制御部203は、補正振幅スペクトル250を逆変換部204へ供給する。逆変換部204は、位相制御部202から供給された強調信号位相スペクトル240と、振幅制御部203から供給された補正振幅スペクトル250とを合成して逆変換を行い、強調信号として、出力端子207に供給する。《変換部201の構成》図3は、変換部201の構成を示すブロック図である。図3に示すように、変換部201はフレーム分割部301、窓がけ処理部windowingunit302、及びフーリエ変換部303を含む。劣化信号サンプルは、フレーム分割部301に供給され、K/2サンプル毎のフレームに分割される。ここで、Kは偶数とする。フレームに分割された劣化信号サンプルは、窓がけ処理部302に供給され、窓関数windowfunctionであるwtとの乗算が行なわれる。第nフレームの入力信号yntt=0,1,...,K/2-1に対するwtで窓がけwindowingされた信号は、次式1で与えられる。また、窓がけ処理部302は、連続する2フレームの一部を重ね合わせオーバーラップして窓がけしてもよい。オーバーラップ長としてフレーム長の50を仮定すれば、t=0,1,...,K/2-1に対して、以下の式2で得られる左辺が、窓がけ処理部302の出力となる。窓がけ処理部302は、実数信号に対しては、左右対称窓関数を用いても良い。また、窓関数は、MMSESTSA法における抑圧係数を1に設定したとき、又はSS法においてゼロを減算したときの入力信号と出力信号が計算誤差を除いて一致するように設計される。これは、wtwtK/2=1となることを意味する。以後、連続する2フレームの50をオーバーラップして窓がけする場合を例として説明を続ける。窓がけ処理部302は、wtとして、たとえば、次式3に示すハニング窓を用いても良い。このほかにも、ハミング窓、三角窓など、様々な窓関数が知られている。窓がけされた出力はフーリエ変換部303に供給され、劣化信号スペクトルYnkに変換される。劣化信号スペクトルYnkは位相と振幅に分離され、劣化信号位相スペクトルargYnkは、位相制御部202に、劣化信号振幅スペクトル|Ynk|は、雑音抑圧部205に供給される。既に説明したように、振幅スペクトルの代わりにパワースペクトルが利用されても良い。《逆変換部204の構成》図4は、逆変換部204の構成を示すブロック図である。図4に示すように、逆変換部204は逆フーリエ変換部401、窓がけ処理部402及びフレーム合成部403を含む。逆フーリエ変換部401は、振幅制御部203から供給された補正振幅スペクトル250と位相制御部202から供給された強調信号位相スペクトルargXnkとを乗算して、強調信号以下の式4の左辺を求める。逆フーリエ変換部401は、得られた強調信号に逆フーリエ変換を施す。逆フーリエ変換された強調信号は、1フレームがKサンプルを含む時間領域サンプル値系列xntt=0,1,...,K-1として、窓がけ処理部402に供給され、窓関数wtとの乗算が行なわれる。第nフレームの入力信号xntt=0,1,...,K/2-1に対してwtで窓がけされた信号は、次式5の左辺で与えられる。また、窓がけ処理部302は、連続する2フレームの一部を重ね合わせオーバーラップして窓がけしてもよい。フレーム長の50をオーバーラップ長として仮定すれば、t=0,1,...,K/2-1に対して、以下の式の左辺が、窓がけ処理部402の出力となり、フレーム合成部403に伝達される。フレーム合成部403は、窓がけ処理部402からの隣接する2フレームの出力を、K/2サンプルずつ取り出して重ね合わせ、以下の式7によって、t=0,1,...,K-1における出力信号式7の左辺を得る。得られた出力信号は、フレーム合成部403から出力端子207に伝達される。なお、図3と図4において変換部201と逆変換部204における変換をフーリエ変換として説明したが、変換部201、逆変換部204は、フーリエ変換に代えて、アダマール変換、ハール変換、ウェーブレット変換など、他の変換を用いても良い。変換部201、逆変換部204がハール変換を用いた場合には、乗算が不要となり、LSI化したときの面積を小さくすることができる。変換部201、逆変換部204がウェーブレット変換を用いた場合には、周波数によって時間解像度を異なったものに変更できるために、雑音抑圧効果の向上が期待できる。また、変換部201において得られる周波数成分を複数統合してから、雑音抑圧部205で実際の抑圧を行っても良い。その際、聴覚特性の弁別能力が高い低周波領域から、能力が低い高周波領域に向かって、よりたくさんの周波数成分を統合して、高い音質を達成することができる。このように、複数の周波数成分を統合してから雑音抑圧を実行すると、雑音抑圧を適用する周波数成分の数が少なくなり、全体の演算量を削減することができる。《雑音抑圧部205の構成》雑音抑圧部205は、変換部201から供給される劣化信号振幅スペクトルを用いて雑音を推定し、推定雑音スペクトルを生成する。雑音抑圧部205は、変換部201からの劣化信号振幅スペクトルと生成した推定雑音スペクトルとを用いて抑圧係数を求める。雑音抑圧部205は、求めた抑圧係数を劣化信号振幅スペクトルに乗じ、強調信号振幅スペクトルとして、振幅制御部203へ供給する。雑音の推定には、M.Kato,A.Sugiyama,andM.Serizawa,′′NoisesuppressionwithhighspeechqualitybasedonweightednoiseestimationandMMSESTSA,′′IEICETrans.FundamentalsJapaneseEdition,vol.J87-A,no.7,pp.851-860,July2004.文献1やR.Martin,′′Spectralsubtractionbasedonminimumstatistics,′′EUSPICO-94,pp.1182-1185,Sept.1994文献2に示されるように、様々な推定方法が利用できる。たとえば、推定雑音スペクトルを、目的音が発生していないフレームの劣化信号振幅スペクトルの平均値とする方法文献1がある。この方法では目的音の発生を検出する必要がある。目的音の発生している区間は、強調信号のパワーで判断することができる。理想的な動作状態として、強調信号は雑音以外の目的音となっている。また、目的音や雑音のレベルは、隣接フレーム間で大きく変化しない。これらのことから、1フレーム過去の強調信号レベルが雑音区間判定の指標とされる。1フレーム過去の強調信号パワーが一定値以下の時には、現フレームは雑音区間と判定される。雑音スペクトルは、雑音区間と判定されたフレームの劣化信号振幅スペクトルを平均化することで推定される。また、推定雑音スペクトルを、劣化信号振幅スペクトルが供給され始めた推定初期の平均値とする方法文献1もある。この場合、推定が開始された直後には目的音が含まれないという条件を満たす必要がある。条件が満たされる場合、推定初期の劣化信号振幅スペクトルが推定雑音スペクトルとされる。さらに、推定雑音スペクトルを、統計的な劣化信号振幅スペクトルの最小値から求める方法文献2もある。この方法では、統計的に一定時間における劣化信号振幅スペクトルの最小値を保持し、その最小値から雑音スペクトルを推定する。劣化信号振幅スペクトルの最小値は、雑音スペクトルのスペクトル形状と似ているため、雑音スペクトル形状の推定値として用いることができる。しかし、最小値では、本来の雑音レベルより小さくなる。そのため、最小値を適切に増幅させたものを推定雑音スペクトルとして用いる。雑音抑圧部205は、様々な抑圧を行うことが可能である。抑圧方法には、代表的なものとして、SSSpectrumSubtraction:スペクトル減算法とMMSESTSAMinimumMean-SquareErrorShort-TimeSpectralAmplitudeEstimator:最小二乗平均誤差短時間振幅スペクトル推定法とがある。雑音抑圧部205がSS法を用いる場合は、雑音抑圧部205は変換部201から供給された劣化信号振幅スペクトルから推定雑音スペクトルを減算する。雑音抑圧部205がMMSESTSA法を用いる場合は、雑音抑圧部205は変換部201から供給された劣化信号振幅スペクトルと生成した推定雑音スペクトルとを用いて、抑圧係数を計算し、この抑圧係数を劣化信号振幅スペクトルに乗算する。この抑圧係数は、強調信号の平均二乗パワーを最小化するように決定される。雑音抑圧部205は、様々な雑音推定を行うことが可能である。推定方法には、たとえば目的音が発生していないフレームの劣化信号振幅スペクトルを利用し、その平均値を推定雑音スペクトルとする方法がある。また、劣化信号振幅スペクトルが供給され始めた初期段階の平均値を推定雑音スペクトルとする方法もある。さらには、統計的に劣化信号振幅スペクトルの最小値を観測し、その最小値を推定雑音スペクトルとする方法などもある。《位相制御部202及び振幅制御部203の構成》図5は位相制御部202及び振幅制御部203の構成を示すブロック図である。図5に示すように、位相制御部202は位相回転部501と回転量生成部502とを含む。振幅制御部203は、補正量算出部503と振幅補正部504とを含む。回転量生成部502は、劣化信号位相スペクトルの回転量を生成し、位相回転部501と補正量算出部503に供給する。位相回転部501は、変換部201から劣化信号位相スペクトル220が供給され、回転量生成部502から回転量が供給される。位相回転部501は、劣化信号位相スペクトル220を供給された回転量だけ回転シフトさせ、強調信号位相スペクトル240として逆変換部204へ供給する。補正量算出部503は、回転量生成部502から供給される回転量に基づいて、振幅の補正係数を決定し、振幅補正部504に供給する。回転量生成部502は、たとえば乱数によって回転量を生成する。位相回転部501は、回転量生成部502が乱数により生成した回転量を用いて劣化信号位相スペクトルを各周波数で回転させる。位相回転部501が行う回転により、劣化信号位相スペクトル220の形状が変化する。この形状の変化により、雑音の特徴を弱めることができる。乱数には、その発生確率が一様な一様乱数や、発生確率が正規分布を示す正規乱数などがある。ここでは一様乱数による回転量の生成法を説明する。一様乱数は線形合同法などで発生する。ここでは、線形合同法での一様乱数を例に挙げ、説明を続ける。線形合同法で発生させた一様乱数は、0~2^M-1の範囲に一様に分布する。ここで、Mは任意の整数であり、^はべき乗を表している。位相の回転量φは、0~2πの範囲に分布させる必要がある。そこで、発生させた一様乱数は、0~2πの範囲に分布されるように変換される。変換は、以下の式8で行う。ここで、Rは一様乱数である。Rmaxはその一様乱数が発生しうる最大の値である。前述した線形合同法で発生させる場合には、Rmax=2^M-1となる。計算の簡単化のために、Rの値は、そのまま回転量とされても良い。回転量であるので2πはちょうど一回転を表す。位相を2π回転させた場合は回転させていない場合と同一である。よって、2παという回転量は、回転量がαの場合と同じになる。ここでは、線形合同法により一様乱数を発生させた場合を説明したが、それ以外の方法で一様乱数を発生させた場合でも、式8により回転量φを求めればよい。位相回転部501は、回転量生成部502から回転量を受け取り、劣化信号位相スペクトルを回転させる。もし、劣化信号位相スペクトルが角度で表現されている場合には、位相回転部501は、その角度に回転量φの値を加算することで回転させる。劣化信号位相スペクトルが複素数の正規ベクトルで表現されている場合には、位相回転部501は、回転量φの正規ベクトルを求め、劣化信号位相スペクトルに乗算することで回転させる。回転量φの正規ベクトルは、式9で求められる。ここで、Φは回転ベクトルである。jはsqrt-1を示す。なお、sqrtは二乗根を表す。補正量算出部503による、補正係数算出方法について説明する。まず、位相の回転により出力レベルの低下が生じることを、図6及び図7を用いて説明する。図6及び図7は、劣化信号が図2に示すブロック図で処理された場合の信号を示している。図6と図7の違いは位相回転の有無である。図6は位相回転を行わない場合の信号を、図7は位相回転をフレーム3から行った場合の信号を示している。まず、位相を回転しない場合の信号について図6を用いて説明する。図6の最上部に描かれているのは、劣化信号である。劣化信号はフレーム分割部301においてフレーム分割される。点線で区切られた、上から2番目の信号が、フレーム分割後の信号である。図6では、連続した4フレーム分の信号を図示した。また、フレームのオーバーラップ率は50としている。窓掛け処理部302は、フレームに分割された信号に窓掛けを行う。点線で区切られた、上から3番目の信号が窓掛け処理後の信号である。図6では位相回転による影響を明確に示すため、矩形窓による重み付けをしている。次に、フーリエ変換部303によって、周波数領域の信号に変換されるが、図6では周波数領域での信号は省略した。位相回転の点線より下部は、逆変換部204の逆フーリエ変換部401により、時間領域へ変換された信号を図示している。点線で区切られた、上から4番目の信号が位相回転後の信号である。ただし、図6では位相回転を行っていないため、窓掛け処理後の信号から変化していない。逆変換部204の逆フーリエ変換部401から出力された強調信号は、窓掛け処理を再度実施される。図6では、矩形窓による重み付けを実施した場合を示している。窓掛け処理された信号は、フレーム合成部403において、合成される。この時、フレーム間の時間を揃える必要がある。フレームのオーバーラップ率が50であるので、ちょうど半分ずつフレームが重なる。位相の回転を実施していない場合、図6のように入力信号と出力信号は一致する。一方、位相を回転する場合の信号について図7を用いて説明する。図7に示したのは位相回転をフレーム3から実施した場合の信号である。最上部に描かれているのは、図6と同じ劣化信号である。フレーム分割後及び窓掛け処理後の信号も図6と同様である。図7では、フレーム3から一定の位相回転を実施した場合を描いている。位相回転処理の点線下部に示した右向き三角形の区間に注目する。位相回転処理により、フレーム3及び4の信号が時間方向にシフトしている。位相回転を施した信号は、再度窓掛け処理が行われ、フレーム合成される。このとき、フレーム2とフレーム3が重なるiiの区間で、フレーム2とフレーム3の信号に違いが生じる。これにより、フレーム合成後の出力信号レベルがiiの区間において小さくなる。つまり、位相の回転を実施する場合、図7のiiの区間で出力信号レベルが低下する。この位相回転による出力信号レベルの低下は、時間領域における加算を周波数領域の加算に置き換え、周波数領域のベクトル合成でも説明できる。図8に、フレーム分割及び窓掛け処理後の連続した2フレームの劣化信号を、x1[n]及びx2[m]として示す。なお、オーバーラップ率は50としている。ここで、nはx1の離散時間を表す。mはx2の離散時間を表す。オーバーラップ率が50の場合には、以下の式10が成立する。また、x1とx2の関係は、以下の式11のようになる。まず、時間領域信号から周波数領域信号への変換及び逆変換の式を示す。周波数領域信号X[k]は、時間領域信号x[n]のフーリエ変換によって、以下の式12のように表現される。ここで、kは離散周波数を表し、Lはフレーム長である。また、周波数領域信号X[k]を逆変換により、時間領域信号x[n]に戻すと、以下の式13のように表現される。これに基づいて、時間領域信号x1[n],x2[m]を、周波数領域信号X1[k],X2[k]に変換すると、以下の式14、式15のように表現される。周波数領域信号X1[k],X2[k]を、逆変換によってそれぞれ時間領域信号x1[n],x2[m]に戻すと、式13より、以下の式16、式17のように表現される。逆変換部は、逆フーリエ変換により周波数領域信号を時間領域信号に変換する。その後、フレーム合成部403は、前フレームと現フレームの強調音声をオーバーラップ加算する。たとえば、図示した例のオーバーラップ率50では、フレーム合成部403は、離散時間m=L/2~L-1の区間で隣接フレームの加算を行う。この加算区間m=L/2~L-1を考える。時間領域の加算に、式16及び式17を代入すると、以下の式18のように表現される。さらに、式18中の周波数領域信号X1[k],X2[k]に式14,式15を代入すると、以下の式19のように表現される。さらに、式19を展開すると、以下の式20のように表現される。ここで、式20の各項に含まれる総和演算を考える。任意の整数gを導入し、以下の式21が成立する。デルタ関数δ[g]のフーリエ逆変換式は、式22で示される。ここで、デルタ関数δ[g]は、以下の式23で示される。式22より、式21は、以下の式24に変形できる。式24の関係から、式20は、以下の式25で表わされる。よって、式20は、以下の式26となる。さて、ここで周波数領域信号X2[k]に対し、位相回転を行った場合を考える。このときの時間領域信号は、図9のようになる。X2[k]の位相スペクトルを、φ[k]回転させたとき、その逆変換は以下の式27となる。これを、式18に代入すると、以下の式28が成立する。これを展開すると、以下の式29が成立する。ここで、オーバーラップ率50と仮定し、オーバーラップ区間のn=L/2~L-1について考える。オーバーラップ区間では、式11より、以下の式30のように展開できる。ここで、それぞれの項にある括弧内の項式31に示すは、ベクトル合成であるから、特定の周波数kに注目すると、図10のように描ける。もし、位相回転が行われていないとき、つまりφ[k]=0の場合は、図11のようになる。式31の絶対値を求めると、以下の式32となる。よって、式31の絶対値が最大になる条件は、φ[k]=0の場合であり、その値は2である。つまり、位相回転が行われると、出力信号の大きさが小さくなることがわかる。この出力信号レベルの低下量を補正するように、補正量算出部503は、強調信号振幅スペクトルの振幅補正量を決定する。ここでは、位相回転量が一様乱数によって決定されると仮定し、補正量の算出方法を具体的に説明する。ここでは、問題簡単化のため、位相の回転による大きさの変動に着目し、それぞれの周波数成分が単位ベクトルに正規化されているものとする。まず、位相回転を行わない場合を考える。連続するフレーム間で位相が同じ場合の合成ベクトルは、図11に示されるベクトルSのようになり、そのベクトルの大きさ|S|は、以下の式33で表わされる。一方、一様乱数により位相回転を行う場合、連続するフレーム間位相差φは、-πからπまでの範囲に一様に分布する。この連続するフレーム間で位相が異なる場合の合成ベクトルは、図10に示されるベクトルS′のようになる。そのベクトルの大きさ|S′|は、以下の式34で表わされる。さて、ここで、期待値E|S′|^2を求めると、以下の式35のようになる。ここで、φが-πからπまで一様分布することから、以下の式36のようになる。このため、|S′|^2の期待値は、以下の式37で表わされる。さて、位相の回転をしない場合の、期待値E|S|^2は、式33より以下の式38で表わされる。式37と式38の比を取ると、以下の式39が成立する。つまり、位相を一様乱数で回転させる場合、出力信号のパワー平均値は入力と比較して、1/2だけ小さくなる。振幅補正部504は振幅値に対し補正を行う。そのため、補正量算出部503は補正係数をsqrt2とし、振幅補正部504に伝達する。ここでは、一様乱数による回転量生成の場合を例に挙げて説明したが、正規乱数などでもその分散と平均値が決まれば、補正係数は一意に定められる。正規乱数を使う場合の補正係数の導出を以下に説明する。正規乱数の場合φの発生確率が正規分布により決定される。そのため、正規乱数による位相回転を実施した場合のパワー期待値を求めるには、φの発生確率に基づいた重み付けが行われる必要がある。具体的には、φの発生確率に基づいた、重み関数fφを導入する。その重み関数fφにより、cosφは、重みづけされる。さらに、重み付け関数fφの積分値で正規化することにより、パワー期待値を求めることができる。正規乱数による位相回転を行った場合の出力パワー期待値ES′′^2は、一様乱数の出力パワー期待値である式35に、重み付け関数fφ及びその積分値を導入することで、以下の式40のように表現される。重み関数fφは、正規分布により表現できるので、以下の式41が成立する。ここで、σは分散をμは平均を表している。たとえば、平均値μ=0,分散σ=1の標準正規分布では、以下の式42となるので、となるので、これを式40に代入すると、式43のようになる。ここで、式43の右辺、第2項を数値計算すると、式44が成立するので、位相の回転をしない場合の、E|S^2|との比は、式45で表わされる。補正量算出部503は、位相を標準正規分布の正規乱数で回転させる場合、振幅補正部504に補正係数をsqrt1/0.805として、伝達する。位相制御部202は、位相回転をすべての周波数に対して行っても良いし、一部の周波数に対して行っても良い。振幅制御部203は、振幅補正を、位相回転が行われた周波数に対してのみ行う。よって、位相回転を行わない周波数の補正係数は、1.0とする。位相回転を行った周波数の補正係数のみ導出した値とする。《振幅補正部504の説明》振幅制御部203は、位相制御部202から伝達された位相の回転量を用いて振幅補正係数を算出する。振幅制御部203は、算出した振幅補正係数を雑音抑圧部205から供給された強調信号振幅スペクトルに乗じ、逆変換部204に供給する。強調信号振幅スペクトルに振幅補正係数を乗じることにより、雑音抑圧装置200は、劣化信号位相スペクトルを回転させ、強調信号位相スペクトルを得た場合の出力レベル低下を解消することができる。以上のように、本実施形態に係る雑音抑圧装置200は、位相スペクトルを回転操作することによる出力信号レベルへの影響を振幅制御部203により取り除くことができる。そのため、雑音抑圧装置200は、高品質な強調信号を得ることができる。第3実施形態本発明の第3実施形態について、図12を用いて説明する。本実施形態では、位相の回転量の上限を制限する点で第2実施形態と異なる。それ以外の構成及び動作については第2実施形態と同様であるためここでは詳しい説明を省略する。図12は、本実施形態に係る位相制御部1202及び振幅制御部203の構成を示すブロック図である。図12に示すように、本実施形態における位相制御部1202は、第2実施形態で説明した回転量生成部502と位相回転部501に加え、回転量制限部1221を含む。回転量生成部502は、回転量制限部1221からの制限を受けつつ、劣化信号位相スペクトルの回転量を生成し、位相回転部501に供給する。回転量制限部1221は、回転量生成部502によって生成される回転量を一定範囲内に制限する。つまり、回転量制限部1221は、φの分布を、0~2πから任意の範囲に制限する。たとえば、回転量制限部1221は、φの分布を、0~π/2のように制限する。これにより、強調信号位相スペクトルに劣化信号位相スペクトルの特徴がある程度残ることになる。完全に位相をランダムに回転させる場合に比較し、劣化信号の特徴がある程度保持されるため、目的音への影響が少なくなる。そのため、目的音の歪が軽減される。本実施形態に係る本発明は、第2実施形態の効果に加えて、位相の回転量を制限することにより、目的音の劣化を低減させることができる。第4実施形態本発明の第4実施形態について、図13を用いて説明する。本実施形態に係る本発明は、位相成分を遅延させ、フレーム間での位相成分の差分を求めた上でそこから補正量を算出する点で第2実施形態と異なる。つまり、位相制御部1302と振幅制御部1303の内部構成において第2実施形態と異なる。他の構成及び動作は第2実施形態と同様であるためここではその説明を省略する。図13は、本実施形態における位相制御部1302及び振幅制御部1303の構成を示すブロック図である。図13に示すように、第4実施形態における位相制御部1302は、第2実施形態に含まれる回転量生成部502と位相回転部501を含む。位相制御部1302は、回転後の位相を振幅制御部1303に供給する。また振幅制御部1303は、位相成分遅延部1331、補正量算出部1333及び振幅補正部504を含む。位相成分遅延部1331は、位相回転部501から供給された強調信号位相スペクトルを1フレーム分保持し、補正量算出部1333に供給する。補正量算出部1333は位相成分遅延部1331からの1フレーム前の強調信号位相スペクトルと位相回転部501からの現在の強調信号位相スペクトルから振幅補正量を算出し、振幅補正部504に伝達する。位相成分遅延部1331は、強調信号位相スペクトルを1フレーム保持し、1フレーム分遅延させ、補正量算出部1333に供給する。位相成分遅延部1331は、単純に1フレーム分強調信号位相スペクトルを遅延することで実現される。本実施形態に係る本発明は、第2実施形態の効果に加えて、位相回転量を数学的に出力レベルの期待値を導出できない場合でも出力レベルを補正することができる。たとえば、乱数テーブルによる乱数ではその発生確率に偏りがあるが、本実施形態に係る本発明は出力レベルを問題なく補正できる。補正量算出部1333は、前フレームと現フレームの強調信号位相スペクトルから各周波数で合成ベクトルの大きさを求め、その大きさから補正係数を決定する。前フレームの位相をα,現フレームの位相をβとすると、その合成ベクトルの大きさ|S′|は、以下の式46で表わされる。連続するフレームで位相が一致する場合の合成ベクトルの大きさ|S|は、式33において既に導出しているように|S|=2である。よって、振幅補正量は、以下の式47で表わされる。本実施形態における本発明は、この値を、振幅制御部203に供給し、強調信号振幅スペクトルを補正することで、出力レベルの低下を解消することができる。本実施形態において、位相回転部以外の構成及び動作については、第2実施形態と同様であるためここではその説明を省略する。第5実施形態本発明の第5実施形態について、図14を用いて説明する。図14は、本実施形態に係る位相制御部202及び振幅制御部1403の構成を示すブロック図である。本実施形態に係る本発明は、入出力比算出部1431を含む点で、第2実施形態と異なる。入出力比算出部1431は、入力端子206からの劣化信号及び逆変換部204からの強調信号を受け取り、その入出力レベル比を算出する。入出力比算出部1431は、入出力レベル比を補正量算出部1433に供給する。補正量算出部1433は、強調信号のレベルが劣化信号と同等になるように、補正量を算出する。振幅補正部504は、算出された補正量により強調信号振幅スペクトルを補正する。入出力比算出部1431は、劣化信号と強調信号の時間領域信号から、そのレベルの比を求める。第nフレームの劣化信号yntと第nフレームの強調信号xntのレベル比Rは、以下の式48で表わされる。ここで、tはサンプル時間を示す。Lはフーリエ変換のフレーム長を示す。補正量算出部1433は、この比の値Rと位相回転を行った周波数成分の数により、振幅補正量Gを求める。時間領域信号が変換部によりN個の周波数成分に分割され、そのうちM個の位相スペクトルに対し、位相回転を行った場合、振幅補正量Gは、以下の式49のように求められる。また、振幅制御部1403は、位相制御部202から伝達された位相回転の有無の情報より、位相回転が行われた周波数のみで振幅補正を行う。本実施形態において、入出力比算出部1431及び補正量算出部1433以外の構成及び動作については、第2実施形態と同様であるためここではその説明を省略する。本実施形態に係る本発明は、時間領域の信号から補正係数を求めているため、位相回転量がいかなる方法で決定されても、出力レベルを補正できる。第6実施形態本発明の第6実施形態について、図15を用いて説明する。図15は、本実施形態に係る位相制御部202及び振幅制御部1503の構成を示すブロック図である。図15に示すように、本実施形態における振幅制御部1503は、第5実施形態に含まれる入出力比算出部1431に加え、平均化処理部1531を含む。平均化処理部1531以外の構成及び動作については、第5実施形態と同様であるためここではその説明を省略する。平均化処理部1531は、入力端子206から劣化信号を受け取り、平均化処理を施した上で、その平均値を入出力比算出部1431に供給する。また、平均化処理部1531は逆変換部204から強調信号を受け取り、平均化処理を施した上で、その平均値を入出力比算出部1431に供給する。入出力比算出部1431は平均化処理部1531から劣化信号及び強調信号の平均値を受け取り、そのレベル比を算出する。平均化処理部1531は、任意の時間長で劣化信号及び強調信号のレベルを平均化する。具体的には、平均化処理部1531は、移動平均や漏れ積分などを用いて劣化信号及び強調信号のレベルを平均化する。本実施形態に係る本発明は、第5実施形態の効果に加えて、平均化されたレベルを使用するため、補正量の変動が抑制され、出力信号の品質向上が実現できる。第7実施形態本発明の第7実施形態について、図16及び図17を用いて説明する。図16は、本実施形態に係る雑音抑圧装置1600の構成を示す図である。本実施形態に係る雑音抑圧装置1600は、第2実施形態の図2の構成に加え、振幅成分遅延部1611、位相成分遅延部1612及び逆変換部1613を含む。また、振幅制御部1603の内部構成にも差異がある。本実施形態において、振幅成分遅延部1611、位相成分遅延部1612及び振幅制御部1603以外の動作については、第6実施形態と同様であるためここではその説明を省略する。入力端子206に供給された劣化信号は、変換部201と振幅制御部1603に供給される。変換部201は、劣化信号振幅スペクトル230を振幅成分遅延部1611及び逆変換部1613に供給する。また、変換部201は、劣化信号位相スペクトル220を位相制御部202に供給する。位相制御部202は、変換部201から供給された劣化信号位相スペクトル220を回転させ、強調信号位相スペクトルとして逆変換部1613及び位相成分遅延部1612に供給する。また、位相制御部202は、各周波数における位相回転の有無を、振幅制御部1603へ伝達する。逆変換部1613は、変換部201から供給された劣化信号振幅スペクトル230と位相制御部202から供給される劣化信号位相スペクトルを用いて、位相回転によるレベル低下の生じた信号を振幅制御部1603に伝達する。振幅成分遅延部1611は、変換部201からの劣化信号振幅スペクトル230を遅延させ、雑音抑圧部205に供給する。位相成分遅延部1612は、位相制御部202からの強調信号位相スペクトルを遅延させ、逆変換部204に供給する。雑音抑圧部205は、振幅成分遅延部1611から供給される劣化信号振幅スペクトルを用いて、雑音を推定し、推定雑音スペクトルを生成する。逆変換部204は、位相制御部202から位相成分遅延部1612を介して供給された、強調信号位相スペクトル240と、振幅制御部1603から供給された補正振幅スペクトル250とを合成して逆変換を行い、強調信号として、出力端子207に供給する。劣化信号位相スペクトル220は、位相制御部202で制御され、逆変換部1613において時間領域の信号に変換される。振幅制御部1603は、その信号と劣化信号210とを用いて、位相回転によるレベルの変動量を求める。この変動量は位相回転部501による回転処理のみの変動である。そのため、振幅制御部1603は、位相の回転によるレベル変動を正確に捉えることができる。振幅制御部1603は、このレベル比を用いて振幅補正を行うが、求まるレベル比は1フレーム前のものである。そこで、振幅成分遅延部1611と位相成分遅延部1612を導入し、1フレーム前の周波数成分に対して、振幅制御部1603において振幅補正を行う。図17は、本実施形態に係る位相制御部202及び振幅制御部1603の内部構成を説明するためのブロック図である。入出力比算出部1731は、入力端子206から供給される劣化信号と逆変換部1613から供給される、位相回転によるレベル低下分を含んだ信号から、レベル比を算出し、補正量算出部1233に供給する。補正量算出部1233は、位相制御部202から各周波数における位相回転の有無の情報を受け取り、振幅補正量を算出する。振幅補正部504は、その振幅補正量に基づいて、強調信号位相スペクトルを各周波数で補正し、逆変換部204へ供給する。なお、振幅成分遅延部1611は、雑音抑圧部205と振幅制御部1603の間に導入しても構わない。本実施形態に係る雑音抑圧装置1600は、第5実施形態の効果に加え、第5実施形態及び第6実施形態では不可避だった入出力比の遅延を回避でき、より正確な出力レベルの補正を実現できる。第8実施形態本発明の第8実施形態について、図18を用いて説明する。図18に示すように、本実施形態に係る雑音抑圧装置1800は、第2実施形態の構成に加え、フレームオーバーラップ制御部1808を含む。フレームオーバーラップ制御部1808は、変換部201及び逆変換部204において、フレーム分割、合成されるときのオーバーラップ率の制御を行う。フレームオーバーラップ制御部1808は、そのオーバーラップ率を振幅制御部203に供給する。既に説明したとおり、位相回転によるレベル低下はオーバーラップによって生じる。このレベル低下量はオーバーラップ率により変化し、オーバーラップ率が大きくなるほど、低下量も大きくなる。よって、オーバーラップ率が変化した場合には、振幅補正量を制御する必要がある。具体的には、オーバーラップ率50の場合の振幅補正量Gを基準に補正量を求める。オーバーラップ率が0のとき振幅補正は必要無い。また、オーバーラップ率50のとき振幅補正量はGである。そこで、フレーム長Lとオーバーラップ長Qの比を用いて、以下の式50のようになる。ここで、G′はオーバーラップ率による補正を行った振幅補正量である。たとえば、オーバーラップ率50の場合は、Q=L/2であるので、以下の式51が成立する。オーバーラップ率25の場合には、Q=L/4より以下の式52が成立する。振幅制御部203は、式50に基づき、位相制御部202から伝達された補正係数を補正し、強調信号振幅スペクトルを補正する。本実施形態において、フレームオーバーラップ制御部1808以外の構成及び動作については、第2実施形態と同様であるためここではその説明を省略する。本実施形態に係る雑音抑圧装置1800は、第2実施形態の効果に加えて、フレームのオーバーラップ率を自由に設定できる。第9実施形態本発明の第9実施形態について、図19を用いて説明する。図19では、雑音抑圧部205に雑音記憶部1909が接続されている。本実施形態に係る雑音抑圧装置1900は、既知の雑音スペクトルを用いて雑音を抑圧する。雑音記憶部1909は、抑圧すべき雑音のスペクトルを記憶する。その記憶された雑音スペクトルに基づき、雑音抑圧部205が雑音抑圧を行うことで、雑音抑圧装置1900は、抑圧対象の雑音のみを抑圧することができる。ここで、雑音記憶部1909は、雑音のスペクトルとして、その雑音の平均スペクトルや最大スペクトル、最小スペクトルなどを記憶していても良い。あるいは、雑音記憶部1909は、それらを組み合わせを記憶していても良い。本実施形態において、雑音記憶部1909以外の構成及び動作については、第2実施形態と同様であるためここではその説明を省略する。本実施形態に係る雑音抑圧装置1900は、第2実施形態の効果を、特定の雑音にのみ適用することができ、抑圧したい雑音が明確な場合には、より高品質な強調信号を得ることができる。他の実施形態以上説明してきた第1乃至第9実施形態では、それぞれ別々の特徴を持つ雑音抑圧装置について説明したが、それらの特徴を如何様に組み合わせた雑音抑圧装置も、本発明の範疇に含まれる。また、本発明は、複数の機器から構成されるシステムに適用しても良いし、単体の装置に適用しても良い。さらに、本発明は、実施形態の機能を実現するソフトウェアの信号処理プログラムが、システムあるいは装置に直接あるいは遠隔から供給される場合にも適用可能である。したがって、本発明の機能をコンピュータで実現するために、コンピュータにインストールされるプログラム、あるいはそのプログラムを格納した媒体、そのプログラムをダウンロードさせるWWWサーバも、本発明の範疇に含まれる。図20は、第1実施形態を信号処理プログラムにより構成する場合に、その信号処理プログラムを実行するコンピュータ2000の構成図である。コンピュータ2000は、入力部2001と、CPU2002と、出力部2003と、メモリ2004とを含む。CPU2002は、信号処理プログラムを読み込むことにより、コンピュータ2000の動作を制御する。すなわち、CPU2002は、メモリ2004に格納された信号処理プログラムを実行し、第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換するS2011。次に、CPU2002は、所定周波数の位相成分を回転させるS2012。位相成分の回転によって生じる、出力信号における振幅成分又はパワー成分の変化量に応じて、CPU2002は、所定周波数の振幅成分又はパワー成分を補正するS2013。CPU2002は、ステップS2012により回転された位相成分と、ステップS2013により補正された振幅成分又はパワー成分と、を合成するS2014。これにより、第1実施形態と同様の効果を得ることができる。[実施形態の他の表現]上記の実施形態の一部又は全部は、以下の付記のようにも記載されうるが、以下には限られない。付記1第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換手段と、所定周波数の前記位相成分を回転させる第1制御手段と、前記第1制御手段が施した回転による前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御手段と、前記第1制御手段により回転された位相成分と、前記第2制御手段により補正された振幅成分又はパワー成分と、を合成する合成手段と、を備えたことを特徴とする信号処理装置。付記2前記第1制御手段は、前記位相成分の回転量を生成する回転量生成手段を含み、前記第2制御手段は、前記回転量生成手段から提供された前記回転量に基づいて前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正することを特徴とする付記1に記載の信号処理装置。付記3前記第2制御手段は、前記位相成分を遅延させる位相成分遅延手段を含み、遅延前後の前記位相成分を比較することにより、前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正することを特徴とする付記1に記載の信号処理装置。付記4前記第2制御手段は、前記合成手段からの出力と前記混在信号とを比較することにより前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正することを特徴とする付記1に記載の信号処理装置。付記5前記第2制御手段は、前記合成手段からの出力の平均値と前記混在信号の平均値とを比較することにより前記変化量を算出し、算出した前記変化量に応じて前記振幅成分又はパワー成分を補正することを特徴とする付記4に記載の信号処理装置。付記6前記振幅成分又はパワー成分を遅延させる振幅成分遅延手段と、をさらに有し、前記第2制御手段は、前記合成手段からの出力と前記混在信号とを比較することにより前記変化量を算出し、算出した前記変化量に応じて、前記振幅成分遅延手段により遅延した前記振幅成分又はパワー成分を補正することを特徴とする付記1に記載の信号処理装置。付記7前記変換手段は、前記混在信号を、互いに一部オーバーラップした複数フレームに分割する分割手段を含み、前記第2制御手段は、さらに、フレームのオーバーラップ率に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正することを特徴とする付記1乃至6のいずれか1項に記載の信号処理装置。付記8前記第1制御手段は、前記位相成分の回転量を制限する回転量制限手段を含むことを特徴とする付記1乃至7のいずれか1項に記載の信号処理装置。付記9第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換ステップと所定周波数の前記位相成分を回転させる第1制御ステップと、前記第1制御ステップにおける回転によって生じた前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御ステップと、前記第1制御ステップにより回転された位相成分と、前記第2制御ステップにより補正された振幅成分又はパワー成分と、を合成する合成ステップと、を備えたことを特徴とする信号処理方法。付記10第1信号と第2信号とが混在した混在信号を、周波数ごとの位相成分及び振幅成分又はパワー成分に変換する変換ステップと所定周波数の前記位相成分を回転させる第1制御ステップと、前記第1制御ステップにおける回転によって生じた前記振幅成分又はパワー成分の変化量に応じて、前記所定周波数の前記振幅成分又はパワー成分を補正する第2制御ステップと、前記第1制御ステップにより回転された位相成分と、前記第2制御ステップにより補正された振幅成分又はパワー成分と、を合成する合成ステップと、をコンピュータに実行させることを特徴とする信号処理プログラムを格納したプログラム記録媒体。以上、実施形態を参照して本願発明を説明したが、本願発明は以上の実施形態に限定されるものではない。本願発明の構成や詳細には、本願発明のスコープ内で同業者が理解し得る様々な変更をすることができる。この出願は、2010年11月24日に出願された日本出願特願2010-261831を基礎とする優先権を主張し、その開示の全てをここに取り込む。"