(57)【要約】【課題】再生する目的信号と周囲雑音とで、制限される周波数帯域が異なっているために信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりする場合でも、明瞭度を向上させることが可能な信号帯域拡張装置を提供する。【解決手段】  第１の周波数範囲に帯域制限された入力信号に対して音量または周波数特性を変化させる信号処理装置であって、前記第１の周波数範囲とは異なる第２の周波数範囲に帯域制限された集音信号から周波数特性情報を抽出する情報抽出手段と、前記情報抽出手段によって抽出された周波数特性情報に対して前記第１の周波数範囲を含んだ周波数特性情報を推定する周波数特性情報拡張手段と、前記周波数特性情報拡張手段によって得られた周波数特性情報に応じて、前記入力信号の音量または周波数特性を変化させる信号補正手段とを具備して構成するようにした。【選択図】      図２
2010224321.tif 000002
【特許請求の範囲】
【請求項１】
第１の周波数範囲に帯域制限された入力信号に対して音量または周波数特性を変化させる信号処理装置であって、
前記第１の周波数範囲とは異なる第２の周波数範囲に帯域制限された集音信号から周波数特性情報を抽出する情報抽出手段と、
前記情報抽出手段によって抽出された周波数特性情報に対して前記第１の周波数範囲を含んだ周波数特性情報を推定する周波数特性情報拡張手段と、
前記周波数特性情報拡張手段によって得られた周波数特性情報に応じて、前記入力信号の音量または周波数特性を変化させる信号補正手段とを備えることを特徴とする信号処理装置。
【請求項２】
第１の周波数範囲に帯域制限された入力信号に対して音量または周波数特性を変化させる信号処理装置であって、
前記第１の周波数範囲とは異なる第２の周波数範囲に帯域制限されて集音される集音信号に含まれる周囲雑音を抽出する周囲雑音抽出手段と、
前記周囲雑音抽出手段によって抽出された周囲雑音から周波数特性情報を抽出する情報抽出手段と、
前記情報抽出手段によって抽出された周波数特性情報に対して前記第１の周波数範囲を含んだ周波数特性情報を推定する周波数特性情報拡張手段と、
前記周波数特性情報拡張手段によって得られた周波数特性情報に応じて、前記入力信号の音量または周波数特性を変化させる信号補正手段とを備えることを特徴とする信号処理装置。
【請求項３】
前記情報抽出手段は、前記第２の周波数範囲よりも狭い第３の周波数範囲に制限された周波数特性情報を抽出することを特徴とする請求項１または請求項２に記載の信号処理装置。
【請求項４】
事前に取得した信号の第２の周波数範囲の周波数特性情報と第１の周波数範囲の周波数特性情報とを対応づけて記憶する記憶手段を更に有し、
前記周波数特性情報拡張手段は、前記記憶手段に記憶された第２の周波数範囲の周波数特性情報と第１の周波数範囲の周波数特性情報との対応を用いて周波数特性情報の拡張を行うことを特徴とする請求項１乃至請求項３のいずれか１項に記載の信号処理装置。
【請求項５】
前記周波数特性情報拡張手段は、前記第１の周波数範囲に含まれる周波数帯域において、連続性または微分連続性を保つように補正を行って前記第１の周波数範囲を含んだ周波数特性情報を推定することを特徴とする請求項１乃至請求項４のいずれか１項に記載の信号処理装置。
【請求項６】
前記情報抽出手段によって抽出される周波数特性情報は、周波数ごとのマスキングレベルであることを特徴とする請求項１乃至請求項５のいずれか１項に記載の信号処理装置。
【請求項７】
前記情報抽出手段によって抽出される周波数ごとのマスキングレベルは、多項式で近似表現されることを特徴とする請求項６に記載の信号処理装置。
【請求項８】
第１の周波数範囲に帯域制限された入力信号に対して前記第１の周波数範囲とは異なる第４の周波数範囲に周波数帯域を拡張する拡張手段と、
前記第４の周波数範囲と異なる第２の周波数範囲に帯域制限されて集音される集音信号から周波数特性情報を抽出する情報抽出手段と、
前記情報抽出手段によって抽出された周波数特性情報に対して前記第４の周波数範囲を含んだ周波数特性情報を推定する周波数特性情報拡張手段と、
前記周波数特性情報拡張手段によって取得された周波数特性情報と前記第４の周波数範囲とに応じて前記拡張手段によって第４の周波数範囲に拡張された入力信号の音量または周波数特性を変化させる信号補正手段とを備えることを特徴とする信号処理装置。
【請求項９】
第１の周波数範囲に帯域制限された入力信号に対して前記第１の周波数範囲とは異なる第４の周波数範囲に周波数帯域を拡張する拡張手段と、
前記第２の周波数範囲に帯域制限されて集音される集音信号に含まれる周囲雑音を抽出する周囲雑音抽出手段と、
前記周囲雑音抽出手段によって抽出された周囲雑音から周波数特性情報を抽出する情報抽出手段と、
前記情報抽出手段によって抽出された周囲雑音に対して前記第４の周波数範囲を含んだ周波数特性情報を推定する周波数特性情報拡張手段と、
前記周波数特性情報拡張手段によって取得された周波数特性情報と前記第４の周波数範囲とに応じて前記拡張手段によって第４の周波数範囲に拡張された入力信号の音量または周波数特性を変化させる信号補正手段とを備えることを特徴とする信号処理装置。
【請求項１０】
第１の周波数範囲に帯域制限された入力信号に対して音量または周波数特性を変化させる信号処理装置であって、
前記第１の周波数範囲とは異なる第２の周波数範囲に帯域制限された集音信号に含まれる周囲雑音を抽出する周囲雑音抽出手段と、
前記周囲雑音抽出手段によって抽出された周囲雑音に対して前記第１の周波数範囲を含む周波数帯域に拡張する周囲雑音拡張手段と、
前記周囲雑音拡張手段によって帯域が拡張された周囲雑音の音量または周波数特性に応じて、前記入力信号の音量または周波数特性を変化させる信号補正手段とを備えることを特徴とする信号処理装置。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、音声や音楽・オーディオなどの信号に対して明瞭度を向上させる信号処理装置に関する。
【背景技術】
【０００２】
音声や音楽・オーディオなどの信号を再生するときに、音声や音楽・オーディオなどの所望の信号（以降、目的信号と称する）以外の周囲雑音などの影響を受けて目的信号の明瞭度が低下する場合がある。そこで、目的信号の明瞭度を向上させるために、集音した信号に含まれる周囲雑音に応じた信号処理を施す必要がある。従来、このような信号処理方法としては、周囲雑音の音量を用いる手法、周囲雑音の周波数特性を用いる手法（例えば、特許文献１）があった。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２００１−１８８５９９号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、目的信号と周囲雑音とで、制限される周波数帯域が異なっているために信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりする場合がある。このような場合、従来の信号処理装置では、周囲雑音の音量や周波数特性が高精度に求まらないために音質劣化を招き、明瞭度を向上させることができないという課題があった。
【０００５】
また、音声信号や音楽・オーディオ信号などの目的信号に対して、エイリアシングを用いたり非線形関数を用いたり線形予測分析を用いたりするような帯域を拡張する従来技術をそのまま用いて、集音した周囲雑音の帯域を拡張しても、周囲雑音の周波数特性を高精度に推定することはできないという課題があった。
【０００６】
この発明は上記の問題を解決すべくなされたもので、再生する目的信号と周囲雑音とで、制限される周波数帯域が異なっているために信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりする場合でも、明瞭度を向上させることが可能な信号帯域拡張装置を提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記の目的を達成するために、この発明は、第１の周波数範囲に帯域制限された入力信号に対して音量または周波数特性を変化させる信号処理装置であって、前記第１の周波数範囲とは異なる第２の周波数範囲に帯域制限された集音信号から周波数特性情報を抽出する情報抽出手段と、前記情報抽出手段によって抽出された周波数特性情報に対して前記第１の周波数範囲を含んだ周波数特性情報を推定する周波数特性情報拡張手段と、前記周波数特性情報拡張手段によって得られた周波数特性情報に応じて、前記入力信号の音量または周波数特性を変化させる信号補正手段とを具備して構成するようにした。
【発明の効果】
【０００８】
本発明によれば、再生する目的信号と周囲雑音で、制限される周波数帯域が異なっているために信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりする場合でも、明瞭度を向上させることが可能な信号処理装置を提供することができる。  
【図面の簡単な説明】
【０００９】
【図１】この発明に係わる信号処理装置の第１の実施例を適用した通信装置の構成を示す回路ブロック図。
【図２】この発明に係わる信号処理部の第１の実施例の構成を示す回路ブロック図。
【図３】図２に示した信号処理部の周囲雑音推定部の構成例を示す回路ブロック図。
【図４】図２に示した信号処理部の周囲雑音情報帯域拡張部の構成例を示す回路ブロック図。
【図５】図４に示した周囲雑音情報帯域拡張部の辞書格納部における辞書の生成方法の動作を説明するための処理フロー図。
【図６】図２に示した信号処理部の信号特性補正部の構成例を示す回路ブロック図。
【図７】この発明に係わる信号処理装置の第１の実施例を適用した通信装置およびディジタルオーディオプレイヤの構成を示す回路ブロック図。
【図８】この発明に係わる信号処理部の変形例１の構成を示す回路ブロック図。
【図９】図８に示した信号処理部の周囲雑音情報帯域拡張部の構成例を示す回路ブロック図。
【図１０】図９に示した周囲雑音情報帯域拡張部の辞書格納部における辞書の生成方法の動作を説明するための処理フロー図。
【図１１】広帯域マスキング閾値の例を示す図。
【図１２】図８に示した信号処理部の信号特性補正部の構成例を示す回路ブロック図。
【図１３】この発明に係わる信号処理部の変形例３の構成を示す回路ブロック図。
【図１４】図１３に示した信号処理装置の周囲雑音情報帯域拡張部の構成例を示す回路ブロック図。
【図１５】図１４に示した周囲雑音情報帯域拡張部の辞書格納部における辞書の生成方法の動作を説明するための処理フロー図。
【図１６】図１４に示した周囲雑音情報帯域拡張部の閾値補正部の動作を説明するための例を示す図。
【図１７】図１４に示した周囲雑音情報帯域拡張部の辞書格納部における辞書の他の生成方法の動作を説明するための処理フロー図。
【図１８】図１４に示した周囲雑音情報帯域拡張部の辞書格納部における辞書の他の生成方法の動作を説明するための処理フロー図。
【図１９】この発明に係わる信号処理部の第２の実施例を適用した通信装置およびディジタルオーディオプレイヤの構成を示す回路ブロック図。
【図２０】この発明に係わる信号処理装置の第２の実施例の構成を示す回路ブロック図。
【図２１】図２０に示した信号処理装置の周囲雑音推定部と周囲雑音抑圧処理部の構成例を示す回路ブロック図。
【発明を実施するための形態】
【００１０】
以下、図面を参照して、この発明の実施形態について説明する。  
（第１の実施例）
図１は、この発明の一実施形態である通信装置の構成を示すものである。この図に示す通信装置は、例えば携帯電話などの無線通信装置の受信系を示すものであって、無線通信部１と、デコーダ２と、信号処理部３と、ディジタル・アナログ（D/A）変換器４と、ス
ピーカ５と、マイク６と、アナログ・ディジタル（A/D）変換器７と、ダウンサンプリン
グ部８と、エコー抑圧処理部９と、エンコーダ１０とを備えている。本実施形態では、再生する目的信号は、受信した入力信号に含まれる遠端話者の音声信号であるとして説明する。
【００１１】
無線通信部１は、移動通信網に収容される無線基地局と無線通信し、そしてこの無線基地局および移動通信網を通じて通信相手局との間に通信リンクを確立して通信する。
【００１２】
デコーダ２は、無線通信部１が通信相手局から受信した受信データを、事前に決められた時間単位である1フレーム（=20[ms]）ごとに復号して、ディジタルの入力信号x[n] (n=0,1,…2N-1)を得て、フレーム単位で信号処理部３に出力する。ただし、この入力信号x[n]は、サンプリング周波数はfs’[Hz]でfs_wb_low[Hz]からfs_wb_high[Hz]までに帯域制限された広帯域の信号である。ここでは、後述する集音信号z[n]のサンプリング周波数fs[Hz]との関係を、fs’=2fsとする。また、サンプリング周波数fs’[Hz]のときの１フレームのデータ長は2Nサンプルする。つまり、N=20[ms]×fs[Hz]÷1000とする。
【００１３】
信号処理部３は、後述するエコー抑圧処理部８においてエコー低減された集音信号z[n] (n=0,1,…N-1)に応じて、１フレーム単位で入力信号x[n] (n=0,1,…2N-1)に対して信号
補正処理を施し、音量または周波数特性を変化させて、その出力信号をy[n] (n=0,1,…2N-1)としてD/A変換器４とダウンサンプリング部８に出力する。なお、信号処理部３の具体的な構成例については後に詳述する。
【００１４】
D/A変換器４は、上記信号補正された出力信号y[n]をアナログ信号y(t)に変換して、ス
ピーカ５に出力する。スピーカ５は、アナログ信号である出力信号y(t)を音響空間へ出力する。
【００１５】
マイク６は、音を集音してアナログ信号である集音信号z(t)を取得し、A/D変換器７に
出力する。このアナログ信号には、近端話者の音声信号と、それ以外の周囲環境に起因するノイズ成分、出力信号y(t)と音響空間に起因するエコー成分などが混在する。例えばノイズ成分としては、電車などが出す騒音、車などのカーノイズ、人ごみでのストリートノイズなどが挙げられる。本実施形態では、通信装置として近端話者の音声信号は通信相手局との間での通信にて所望される必要な信号であるため、近端話者の音声信号以外の成分を周囲雑音として扱う。
【００１６】
A/D変換器７は、アナログ信号である集音信号z(t)をディジタル信号に変換して、ディ
ジタルの集音信号z’[n] (n=0,1,…N-1)を得て、Nサンプル単位でエコー抑圧処理部８に
出力する。ただし、この集音信号z[n]は、サンプリング周波数はfs[Hz]でfs_nb_low[Hz]
からfs_nb_high[Hz]までに帯域制限された狭帯域の信号である。ただし、fs_wb_low ≦ fs_nb_low < fs_nb_high < fs/2 ≦ fs_wb_high < fs’/2 を満たすものとする。
【００１７】
ダウンサンプリング部８は、出力信号y[n]をサンプリング周波数fs’[Hz]からfs[Hz]にダウンサンプリングして、fs_nb_low[Hz]からfs_nb_high[Hz]までに帯域制限した信号をy’[n] (n=0,1,…N-1)としてエコー抑圧処理部９に出力する。
【００１８】
エコー抑圧処理部９は、ダウンサンプリングされた出力信号y’[n] (n=0,1,…N-1)を利用して、集音信号z’[n] (n=0,1,…N-1)に含まれるエコー成分を低減する処理を行い、そのエコー低減された信号をz[n] (n=0,1,…N-1)として、信号処理部３とエンコーダ１０に出力する。ここでエコー抑圧処理部９は、例えば、特登４０４７８６７号公報や、特開２００６−２０３３５８号公報や、特開２００７−６０６４４号公報などに記載される既存の技術で実施してよい。
【００１９】
エンコーダ１０は、エコー抑圧処理部８においてエコー低減された集音信号z[n] (n=0,1,…N-1)をNサンプルごとに符号化して無線通信部１に出力し、無線通信部１によって送
信データとして通信相手局へ送信される。
【００２０】
次に、信号処理部３の実施例について説明する。以下の説明では、例えば、fs=8000[Hz]、fs’=16000[Hz]、fs_nb_low=340[Hz]、fs_nb_high=3950[Hz]、fs_wb_low=50[Hz]、fs_wb_high=7950[Hz]とする。帯域制限の周波数帯域やサンプリング周波数については、これに限らない。また、ここではN=160とする。
【００２１】
図２は、信号処理部３の構成例を示すものである。信号処理部３は、周囲雑音推定部３１と、周囲雑音情報帯域拡張部３２と、信号特性補正部３３とを備える。これらは、１つのプロセッサと、図示しない記憶媒体に記録されたソフトウェアによって実現することも可能である。
【００２２】
周囲雑音推定部３１は、エコー抑圧処理部８においてエコー低減された信号から近端話者の音声信号以外の信号を周囲雑音と推定し、この周囲雑音を特徴付ける特徴量を抽出する。なお、集音信号z[n]が狭帯域の信号であるため、周囲雑音も狭帯域の信号である。そこで、周囲雑音を特徴付ける特徴量を、狭帯域信号情報と称する。狭帯域信号情報は、パワースペクトル、振幅スペクトルや位相スペクトル、PARCOR係数や反射係数、線スペクトル周波数、ケプストラム係数、メルケプストラム係数など、周囲雑音を特徴付ける特徴量であればどのようなものでも構わない。
【００２３】
周囲雑音情報帯域拡張部３２は、狭帯域信号情報を用いて、周囲雑音を入力信号x[n]の周波数帯域と同じ周波数帯域（広帯域）に拡張した場合にこの周囲雑音を特徴付ける特徴量を推定する。この特徴量を、広帯域信号情報と称する。
【００２４】
信号特性補正部３３は、周囲雑音情報帯域拡張部３２を用いて、目的信号の信号特性を補正する。
【００２５】
このように、周囲雑音が狭帯域の信号であっても、広帯域に拡張した場合の特徴量を推定することによって、信号特性補正部３３での補正処理によって明瞭度を向上させることができる。
【００２６】
以下の説明では、信号処理部３の具体的な構成について説明する。なお、以下の説明では、狭帯域信号情報は周囲雑音のパワースペクトル、広帯域信号情報は周囲雑音を広帯域の信号に拡張した場合のパワー値（広帯域パワー値）であるとして説明する。
【００２７】
図３に周囲雑音推定部３１の構成例を示す。周囲雑音推定部３１は、周波数領域変換部３１１と、パワー算出部３１２と、周囲雑音区間判定部３１３と、周波数スペクトル更新部３１４とを備える。
【００２８】
周囲雑音推定部３１は、エコー抑圧処理部８においてエコー低減された集音信号z[n] (n=0,1,…N-1)から近端話者の音声信号以外である周囲雑音を推定してこの信号のパワースペクトル|N[f,w]|2 を抽出して、周囲雑音情報帯域拡張部３２へ出力する。
【００２９】
周波数領域変換部３１１は、現在のフレームｆの集音信号z[n] (n=0,1,…N-1)が入力される。そして、このフレームの１フレーム前の集音信号から窓掛けによるオーバーラップサンプル数分のサンプルを抽出し、現在のフレームの入力信号と時間方向に結合し、適宜零詰めなどを行って、周波数領域変換に必要なサンプル分の信号を取り出す。次のフレームでの集音信号z[n]のシフト幅と集音信号z[n]のデータ長の比であるオーバーラップは５０%である場合が考えられるが、ここでは例として、１フレーム前とのオーバーラップの
サンプル数をL=48として、１フレーム前の集音信号Ｌサンプルと当該フレームの集音信号z[n]のN=160サンプル分とLサンプル分の零詰めから、２Ｍ=256サンプルを用意するとする。この２Mサンプルに対して正弦波窓による窓関数を乗じることで窓掛けを行う。そして
、窓掛けを行った２Ｍサンプルの信号に対して、周波数領域変換を行う。周波数領域への変換は、例えばＦＦＴの次数を2MとしＦＦＴによって行うことができる。なお、周波数領域変換を施す信号に零詰めすることによってデータ長を２のべき乗（２Ｍ）にし、周波数領域変換の次数を２のべき乗（２Ｍ）にするとしたが、周波数領域変換の次数はこれに限らない。
【００３０】
集音信号z[n]が実信号である場合は、周波数領域変換を施して得られた信号から冗長なＭ=128ビンを除くと、周波数スペクトルZ[f,w]  (w=0,1,…M-1)が得られる、これを出力
する。ただし、ωは、周波数ビンを表す。なお、実信号のとき冗長なのは本来Ｍ-1 (=127)ビンであり、最高域の周波数ビンw=Ｍ (=128)を考慮するべきである。しかしながら、ここで周波数領域変換する信号は、帯域制限された音声信号を含むディジタル信号を前提としており、帯域制限によって最高域の周波数ビンw=Mを考慮しなくても音質に影響を及ぼ
さない。そこで、これ以降説明の簡略化のために、最高域の周波数ビンw=Mを考慮しない
記述にする。勿論、最高域の周波数ビンw=Mを考慮しても構わない。その際、最高域の周
波数ビンw=Mは、w=M-1と同等に扱うか、単独で扱うようにする。
【００３１】
なお、窓掛けに用いる窓関数は、ハミング窓に限定せず、他の対称窓（ハニング窓、ブラックマン窓、正弦波窓など）あるいは音声符号化処理で用いられるような非対称窓などに適宜変更してよい。また、周波数領域変換は、DFT(Discrete Fourier Transform)や離
散コサイン変換（ＤＣＴ: Discrete Cosine Transform）などの周波数領域に変換する他
の直交変換を代用することも可能である。
【００３２】
パワー算出部３１２は、周波数領域変換部３１１から出力された周波数スペクトルZ[f,w] (w=0,1,…M-1)における実部と虚部の２乗和であるパワースペクトル|Z[f,w]|2  (w=0,1,…M-1)を算出して出力する。
【００３３】
周囲雑音区間判定部３１３は、集音信号z[n] (n=0,1,…N-1)と、パワー算出部３１２から出力されるパワースペクトル|Z[f,w]|2  (w=0,1,…M-1)と、周波数スペクトル更新部３１４から出力される１フレーム前の各周波数帯域の周囲雑音のパワースペクトル|N[f-1,w]|2 を用いて、集音信号z[n]に周囲雑音が支配的に含まれている区間（周囲雑音区間）であるか、周囲雑音には含まない近端話者の音声信号と周囲雑音が混在している区間（音声区間）のどちらであるかの判別をフレーム毎に行い、フレーム毎に判定結果を表すフレーム判定情報vad[f]を出力する。ここでは、周囲雑音区間であるときフレーム判定情報vad[f]=0とし、音声区間であるときvad[f]=1とする。なおこれ以降、当該成分のみしか存在しないか、あるいは当該成分が他の成分よりも非常に多く含まれる場合（所定の閾値以上含まれる場合）を「支配的に含まれる」と表現する。
【００３４】
具体的には、集音信号z[n] (n=0,1,…N-1)とパワースペクトル|Z[f,w]|2 と１フレーム前の周囲雑音のパワースペクトル|N[f-1,w]|2 を用いて、複数の特徴量を算出し、フレーム判定情報vad[f]を出力する。ここでは複数の特徴量として、１次自己相関係数Acorr[f,1]、自己相関係数最大値Acorr_max[f]、周波数別ＳＮ比総和snr_sum[f]、周波数別ＳＮ比分散snr_var[f]を例に挙げて説明する。
【００３５】
まず、式（１）に示すように、フレーム単位でのパワーで正規化されて絶対値をとったk次自己相関係数Acorr[f,k] (k=1,…N-1)を計算する。
【数１】
2010224321.tif 000003
このとき併せて、k=1である１次自己相関係数Acorr[f,1]も計算する。１次自己相関係数Acorr[f,1]は０から１の値をとり、０に近づくほどノイズ性が強い。つまり、１次自己相
関係数の値が小さいほど、集音信号に周囲雑音が多く含まれ、周囲雑音には含まない音声信号が少ないと判断される。そして、正規化されたk次自己相関係数Acorr[f,k] (k=1,…N-1)から式（２）に示すように、最大となる自己相関係数Acorr[f,k]を計算して、自己相
関係数最大値Acorr_max[f]とする。自己相関係数最大値Acorr_max[f]は０から１の値をとり、０に近づくほどノイズ性が強い。つまり、自己相関係数の値が小さいほど、集音信号に周囲雑音が多く含まれ、周囲雑音には含まない音声信号が少ないと判断される。
【数２】
2010224321.tif 000004
次に、パワースペクトル|Z[f,w]|2と周囲雑音のパワースペクトル|N[f,w]|2とを入力として、それらの比である各周波数帯域のＳＮ比を、ここではｄＢ表現したsnr[f,w] (w=0,1,…M-1)として式（３）で算出する。
【数３】
2010224321.tif 000005
そして、各周波数帯域のＳＮ比snr[f,w] (w=0,1,…M-1)の和を式（４）で算出し、周波数別ＳＮ比総和値snr_sum[f]とする。周波数別ＳＮ比総和値snr_sum[f]は０以上の値をとり、この値が小さいほど集音信号中にノイズ成分である周囲雑音が多く含まれ、周囲雑音には含まない音声信号が少ないと判断される。
【数４】
2010224321.tif 000006
また、各周波数帯域のＳＮ比snr[f,w] (w=0,1,…M-1)の分散を式（５）で算出し、周波数別ＳＮ比分散値snr_var[f]とする。周波数別ＳＮ比分散値snr_var[f]は０以上の値をとり、この値が小さいほどノイズ成分である周囲雑音が多く含まれ、周囲雑音には含まない音声信号が少ないと判断される。
【数５】
2010224321.tif 000007
最後に、複数の特徴量である、１次自己相関係数Acorr[f,１]、自己相関係数最大値Acorr_max[f]、周波数別ＳＮ比総和値snr_sum[f]、周波数別ＳＮ比分散値snr_var[f]を用い
て、これらにそれぞれ所定の重み付けによる重み付けを行い、複数の特徴量の重み付け和として周囲雑音度合type[f]を算出する。ここでは、周囲雑音度合type[f]が小さいほど周囲雑音が支配的であるとし、大きいほど周囲雑音には含まない音声信号が支配的であるとしているので、例えば、線形識別関数による判定を用いた学習アルゴリズムなどで重みw1、w2、w3、w4（ただしw1≧０、w2≧０、w3≧０、w4≧０）を設定して、式（６）で算出する。そして、周囲雑音度合type[f]が所定の閾値THRよりも大きければvad[f]=1とし、周囲雑音度合type[f]が所定の閾値THR以下であればvad[f]=0とする。
【数６】
2010224321.tif 000008
以上の説明では、複数の特徴量を求める際に、周波数ビンごとに処理するとして説明したが、周波数領域変換による隣接する複数の周波数ビンをまとめてグループを作り、そのグループ単位で処理を行っても構わない。また、フィルタバンクなどの帯域分割フィルタなどの周波数領域変換によって実現してもよい。
【００３６】
なお、前述した複数の特徴量を全て使わなくてもよいし、他の特徴量を追加して用いてもよい。また、無線通信部１あるいはデコーダ２から出力されるコーデック情報、例えば、無音挿入記述子（ＳＩＤ）や音声検出器（ＶＡＤ）による音声であるか音声でないかを表す音声検出情報や擬似背景雑音を生成したかどうかの情報などを用いてもよい。
【００３７】
周波数スペクトル更新部３１４は、周囲雑音区間判定部３１３から出力されるフレーム判定情報vad[f]と、パワー算出部３１２から出力されるパワースペクトル|Z[f,w]|2  (w=0,1,…M-1)を用いて、各周波数帯域の周囲雑音のパワースペクトルである|N[f,w]|2  (w=0,1,…M-1)を推定して出力する。例えば、フレーム判定情報vad[f]を0として周囲雑音が
支配的に含まれる区間（周囲雑音区間）であると判別されたフレームのパワースペクトル|Z[f,w]|2 をフレーム単位で忘却させて平均的なパワースペクトルを算出し、これを各周波数帯域の周囲雑音のパワースペクトル|N[f,w]|2  (w=0,1,…M-1)として出力する。なお、具体的には、各周波数帯域の周囲雑音のパワースペクトル|N[f,w]|2 の算出は、式（８）に示すように１フレーム前の各周波数帯域の周囲雑音のパワースペクトル|N[f-1,w]|2 を用いて再帰的に行う。ただし、式（７）の忘却係数αN[ω]は１以下の係数であって、
好ましくは０．７５〜０．９５程度である。
【数７】
2010224321.tif 000009
周囲雑音情報帯域拡張部３２は、各周波数帯域の周囲雑音のパワースペクトル|N[f,w]|2 を用いて、入力信号x[n]には存在して集音信号z[n]には存在しない周波数帯域成分を含めた信号のパワー値を生成する。
【００３８】
図４は周囲雑音情報帯域拡張部３２の構成例を示す図である。周囲雑音情報帯域拡張部３２は、パワー正規化部３２１と、辞書格納部３２２と、広帯域パワー算出部３２３とを備える。
【００３９】
周囲雑音情報帯域拡張部３２では、狭帯域信号情報から狭帯域特徴量データを算出し、狭帯域信号情報から算出される狭帯域特徴量データと広帯域特徴量データとの対応を事前にモデル化しておき、このモデルと取得した狭帯域特徴量データとの対応を用いて広帯域特徴量データを算出し、広帯域特徴量データから広帯域信号情報を生成する。前述のとおり、ここでは、狭帯域信号情報は周囲雑音のパワースペクトルである。またここでは、広帯域特徴量データと広帯域信号情報は同じであるとし、広帯域信号情報は広帯域パワー値N_wb_level[f]で示される音量である。狭帯域特徴量データと広帯域特徴量データとの対
応のモデル化には、GMM(Gaussian mixture model)を利用する手法を用いる。ここでは、
狭帯域パワー値Pow_N[f]と周囲雑音の正規化したパワースペクトル|Nn[f,w]|2 (w=0,1,…M-1)を次数方向に連結してDnb次の狭帯域特徴量データとして用い、広帯域パワー値N_wb_level[f]をDwb次の広帯域特徴量データとして用いる（Dnb=M+1、Dwb=1）。
【００４０】
まず狭帯域信号情報から狭帯域特徴量データを算出するために、パワー正規化部３２１には、周囲雑音推定部３１から出力される周囲雑音のパワースペクトル|N[f,w]|2 (w=0,1,…M-1)が入力され、この周囲雑音のパワースペクトルを用いて狭帯域特徴量データを算
出する。狭帯域特徴量データの１つは、式（８）に基づいて算出される、パワースペクトルの各周波数ビンの総和である狭帯域パワー値Pow_N[f]である。
【数８】
2010224321.tif 000010
また他の狭帯域特徴量データとしては、狭帯域パワー値Pow_N[f]を用いて式（９）に従って各周波数ビンのパワースペクトル|N[f,w]|2を正規化したパワースペクトル|Nn[f,w]|2
を算出する。
【数９】
2010224321.tif 000011
辞書格納部３２２は、事前に集音した周囲雑音に基づいてDnb次の狭帯域特徴量データ
とDwb次の広帯域特徴量データとの対応をモデル化して学習された混合数Ｑ（ここではQ=64）のGMMの辞書  λ1q={ｗq,μq,Σq}（q=1,…,Q）を格納する。なお、ｗqはq番目の混合正規分布の混合重みを示し、μqはq番目の混合正規分布の平均ベクトル、Σqはq次数目の混合正規分布の共分散行列（対角共分散行列または全共分散行列）を表している。なお、平均ベクトルμqと共分散行列Σqの成分の数である次数は、Dnb＋Dwbである。
【００４１】
辞書格納部３２２における事前の辞書λ1qの学習生成方法について、フローチャートを図５に示し、説明する。
【００４２】
GMMの生成に用いる信号は、入力信号x[n]と同様のサンプリング周波数fs’[Hz]でfs_wb_low[Hz]からfs_wb_high[Hz]までに帯域制限された広帯域な信号を別途事前に集音した信号群である。この信号群は、多数の様々な環境、様々な音量であることが望ましい。以下では、GMMの生成に用いる広帯域信号の信号群をまとめて広帯域信号データwb[n]と表記する。ｎは時刻（サンプル）を表す。
【００４３】
まず、広帯域信号データwb[n]を入力として、ダウンサンプリングフィルタによってサ
ンプリング周波数fs[Hz]にダウンサンプリングし、fs_nb_low[Hz]からfs_nb_high[Hz]ま
での狭帯域に帯域制限された狭帯域信号データnb[n]を得る（ステップＳ１０１）。この
ようにして、集音信号z[n]と同じように帯域制限された信号群を生成する。なお、図示しないが、上記ダウンサンプリングフィルタや帯域制限処理でアルゴリズム遅延が生じる場合には、狭帯域信号データnb[n]を広帯域信号データwb[n]と同期を合わせる処理を行う。
【００４４】
次に、上記狭帯域信号データnb[n]からフレームf単位で狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）を抽出する（ステップＳ１０２）。狭帯域特徴量データPnd[f,d]は、所定の次数の狭帯域信号情報を表す特徴量データである。ステップＳ１０２ではまず、狭帯域信号データnb[n]からフレーム毎に前述の周波数領域変換部３１１における処理と同様に周
波数領域変換処理を行い、M次の狭帯域信号データnb[n]のパワースペクトルを得る（ステップＳ１０２１）。次に、前述のパワー正規化部３２１における処理と同様の処理によって、狭帯域信号データnb[n]からフレーム毎にパワー算出を行い、1次のパワー値を得る（ステップＳ１０２２）。そして、これらのパワースペクトルとパワー値からM次の狭帯域
信号データnb[n]の正規化されたパワースペクトルを得る（ステップＳ１０２３）。そし
て、M次の正規化されたパワースペクトルと1次のパワー値をフレーム単位で次数方向（次元方向）に連結して、次数Dnb（=M+1）の狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）を
生成する（ステップＳ１０２４）。
【００４５】
一方、上記に並行して、広帯域信号データwb[n]からフレームf単位で広帯域特徴量データPwb[f,d]（d=1,…,Dwb）を抽出する（ステップＳ１０３）。広帯域特徴量データPwb[f,d]は、所定の次数の広帯域信号情報を表す特徴量データである。ステップＳ１０３ではまず、広帯域信号データwb[n]からフレーム毎に前述の周波数領域変換部３１１における処
理のＦＦＴ点数を倍の4M点にして、同様に周波数領域変換処理を行い、2M次の広帯域信号データwb[n]のパワースペクトルを得る（ステップＳ１０３１）。次に、前述のパワー正
規化部３２１における処理と同様の処理によって、広帯域信号データwb[n]からフレーム
毎にパワー算出を行って1次のパワー値を得る。このパワー値を次数Dwb（=1）の広帯域特徴量データPwb[f,d]とする（ステップＳ１０３２）。
【００４６】
次に、狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）と広帯域特徴量データPwb[f,d]（d=1,…,Dwb）の時間的に同期が取れた２つの特徴量データをフレーム単位で次数方向（次元方向）に連結して、次数Dnb+Dwbの連結特徴量データP[f,d]（d=1,…,Dnb+Dwb）を生成す
る（ステップＳ１０４）。
【００４７】
そして、連結特徴量データＰ[f,d]から混合数Q=1の初期GMMを生成し、各GMMの平均ベクトルをわずかにずらして別の混合分布を生成することで混合数Qを2倍に増やす処理と、連結特徴量データＰ[f,d]を用いてＥＭアルゴリズムにより収束するまでGMMの尤度最大化学習を行う処理とを交互に繰り返し行い、混合数Ｑ（ここではQ=64）のGMM  λ1q={ｗq,μq,Σq}（q=1,…,Q）を生成する（ステップＳ１０５）。EMアルゴリズムについては、D.A.Reynols and R.C.Rose,“Robust text-independent speaker identification using Gaussian mixture models”,IEEE Trans. Speech and Audio Processing, Vol.3, no.1, pp.72-83, Jan.1995. などの文献に詳細な記述がある。  
図４の説明に戻る。広帯域パワー算出部３２３には、パワー正規化部３２１から出力された狭帯域パワー値Pow_N[f]と周囲雑音の正規化したDnb次のパワースペクトル|Nn[f,w]|2 (w=0,1,…M-1)が連結されて、狭帯域特徴量データPn_nb[f]（d=1,…,Dnb）として入力
される。また、広帯域パワー算出部３２３は、辞書格納部３２２からGMMの辞書λ1q={ｗi,μq,Σq}（q=1,…,Q）を読み出して、最小平均２乗誤差（MMSE：Minimam Mean Square Error）推定に従って、式（１０）に示すように、複数の正規分布モデルによるソフトクラスタリングと連続的な線形回帰によって、周波数帯域が拡張された広帯域に対応する特徴量データへの変換を行い、狭帯域特徴量データPn_nb[f]から広帯域特徴量データである広帯域パワー値N_wb_level[f]を算出して出力する。式（１０）は、次元（d=1,…,Dnb+Dwb
）方向のベクトルとして記載している。また、平均ベクトルμq（d=1,…,Dnb+Dwb）は次
元方向で、μqN（d=1,…,Dnb）とμqW（d=Dnb,…,Dnb+Dwb）に分割し、（Dn+Dw）×（Dn+Dw）行列である共分散行列Σqも以下のように、Dn×Dn行列であるΣqNNとDn×Dw行列であるΣqNWとDw×Dn行列であるΣqWNとDw×Dw行列であるΣqWWとに分割する。
【数１０】
2010224321.tif 000012
周囲雑音情報帯域拡張部３２では、広帯域特徴量データと広帯域信号情報は同じであるとしたため、このようにして、狭帯域信号情報である周囲雑音のパワースペクトル|N[f,w]|2から、広帯域信号情報である広帯域パワー値N_wb_level[f]が得られる。
【００４８】
図６に信号特性補正部３３の構成例を示す。信号特性補正部３３は、周波数領域変換部３３１と、補正度合決定部３３２と、補正処理部３３３と、時間領域変換部３３４とを備える。信号特性補正部３３には、入力信号x[n] (n=0,1,…2N-1)と広帯域パワー値N_wb_level[f]が入力され、入力信号x[n]が集音信号に含まれる周囲雑音に埋もれてしまわないよう明瞭化する信号補正処理を行い、その補正後の出力信号y[n] (n=0,1,…2N-1)を出力す
る。
【００４９】
周波数領域変換部３３１には、周波数領域変換部３１１における集音信号z[n] (n=0,1,…N-1)の代わりに入力信号x[n] (n=0,1,…2N-1)が入力される。周波数領域変換部３３１
は、周波数領域変換部３１１と同様の処理によって、入力信号x[n]の周波数スペクトルX[f,w]を出力する。例えば、周波数領域変換部３３１は、１フレーム前とのオーバーラップのサンプル数をL=96とし、１フレーム前の入力信号Lサンプルと当該フレームの入力信号x[n]の2N=320サンプル分とLサンプル分の零詰めから、4M=512サンプルを用意する。そして、この4Mサンプルに対して正弦波窓による窓関数を乗じることで窓掛けを行った信号に対して、ＦＦＴの次数を4MとしてＦＦＴによる周波数領域変換を行う。
【００５０】
補正度合決定部３３２には、周囲雑音情報帯域拡張部３２から出力された広帯域パワー値N_wb_level[f]を入力される。そして、式（１１）で補正ゲインG[f,w] (w=0,1,…2M-1)を算出して出力する。
【数１１】
2010224321.tif 000013
式（１１）のＮ０は、通常の利用環境における周囲雑音のパワーを入力信号x[n]と同じサンプリング周波数・同じ帯域制限で事前に測定しておいて設定した周囲雑音の基準パワー値である。このようにすることで、通常の利用環境よりも周囲雑音のパワー値が大きい環境（すなわち、周囲雑音が多く含まれる環境）でも、補正ゲインG[f,w]をその分大きく設定することで、入力信号x[n]を明瞭化することができる。
【００５１】
補正処理部３３３には、入力信号x[n]の周波数スペクトルX[f,w] (w=0,1,…2M-1)と補
正度合決定部３３２から出力された補正ゲインG[f,w] (w=0,1,…2M-1)が入力される。そ
して、式（１２）によって入力信号x[n]の周波数スペクトルX[f,w]を補正し、その補正結果である出力信号y[n]の周波数スペクトルY[f,w] (w=0,1,…2M-1)を出力する。
【数１２】
2010224321.tif 000014
時間領域変換部３３４は、補正処理部３３３から出力された周波数スペクトルY[f,w] (w=0,1,…2M-1)に対して時間領域変換（周波数逆変換）を行って、周波数領域変換部３３
１における窓掛けを考慮しオーバーラップを戻す処理を適宜行い、補正された信号である出力信号y[n]を算出する。例えば、周波数スペクトルY[f,w] (w=0,1,…2M-1)に対して、
入力信号x[n]が実信号であることを考慮に入れて周波数スペクトルY[f,w]をw=0,1,…4M-1まで復元した上で、4M点のＩＦＦＴ（Inverse Fast Fourier Transform）を行い、窓掛けを考慮し１フレーム前の補正された信号である出力信号y[n]を用いてオーバーラップを戻し、出力信号y[n]を算出する。
【００５２】
以上のように、再生される入力信号と集音信号で、信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりしていても、集音信号の音量について入力信号の周波数帯域を加味して拡張し推定することで、集音信号の音量が高精度に求まり、入力信号の明瞭度を向上させることができる。
【００５３】
なお、上述の説明では、本発明を通信装置に適用した場合について説明したが、図７（ａ）に示すように、ディジタルオーディオプレイヤに本発明を適用することも可能である。このディジタルオーディオプレイヤは、フラッシュメモリやHDD(Hard Disk Drive)を用いた記憶部１１を備え、この記憶部１１から読み出した音楽・オーディオデータをデコーダ２が復号する。このとき復号して再生したい所望の信号である目的信号は音楽・オーディオ信号である。このディジタルオーディオプレイヤのマイク６で集音される集音信号z(t)は、近端話者の音声や周囲環境に起因するノイズ成分、出力信号y(t)と音響空間に起因するエコー成分などで構成され、音楽・オーディオ信号は含まれない。この場合は、
通信装置と異なって近端話者の音声は不要であるため、近端話者の音声を含んだこれら全ての成分を周囲雑音として取り扱う。
【００５４】
また、図７（ｂ）に示すように、本発明を通信装置に適用して、音声帯域拡張通話装置に適用することも可能である。この音声帯域拡張通話装置は、デコーダ２Ａを備え、デコーダ２Ａと信号処理部３の間に信号帯域拡張処理部１２を備えた構成である。そして、この場合の信号処理部３は、帯域が拡張された入力信号x’[n]に対して前述の処理をする。
【００５５】
なお、信号帯域拡張処理部１２で行われる処理は、fs_nb_low[Hz]からfs_nb_high[Hz]ま
でに帯域制限された狭帯域の入力信号をfs_wb_low[Hz]からfs_wb_high[Hz]までの広帯域
の信号に帯域を拡張する処理であって、例えば、特登３１８９６１４号公報や特登３２４３１７４号公報や特開平９−５５７７８号公報などに記載される既存の技術で実施してよい。
【００５６】
（信号処理部の変形例１）
次に、信号処理部にて用いる狭帯域信号情報は周囲雑音のパワースペクトル、広帯域信号情報は周囲雑音を広帯域の信号に拡張した場合のマスキング閾値（広帯域マスキング閾値）である場合を例にして説明する。
【００５７】
図８は、その構成を示すものである。信号処理部３０は、信号処理部３で用いていた周囲雑音情報帯域拡張部３２および信号特性補正部３３に代わって、周囲雑音情報帯域拡張部３４と、信号特性補正部３５とを備えて構成される。
【００５８】
図９に周囲雑音情報帯域拡張部３４の構成例を示す。周囲雑音情報帯域拡張部３４は、パワー正規化部３２１と、辞書格納部３４２と、広帯域パワースペクトル算出部３４３と、広帯域マスキング閾値算出部３４４と、パワー制御部３４５とを備える。
【００５９】
周囲雑音情報帯域拡張部３４は、周囲雑音情報帯域拡張部３２と同様に、周囲雑音のパワースペクトルを入力として、入力信号x[n]には存在して集音信号z[n]には存在しない周波数帯域成分を含めた情報（広帯域信号情報）を生成する。つまり、周囲雑音情報帯域拡張部３４では、狭帯域信号情報から狭帯域特徴量データを算出し、狭帯域信号情報から算出される狭帯域特徴量データと広帯域特徴量データとの対応を事前にモデル化しておき、このモデルと取得した狭帯域特徴量データとの対応を用いて広帯域特徴量データを算出し、広帯域特徴量データから広帯域信号情報を生成する。ただし、周囲雑音情報帯域拡張部３４では、狭帯域特徴量データと広帯域特徴量データとの対応のモデル化に、ベクトル量子化によるコードブックを利用する手法を用いる。ここでは、周囲雑音の正規化したパワースペクトル|Nn[f,w]|2 (w=0,1,…M-1)をDnb次の狭帯域特徴量データとして用い、周囲
雑音の正規化した広帯域のパワースペクトル|Nw[f,w]|2 (w=0,1,…2M-1)をDwb次の広帯域特徴量データとして用いる（Dnb=M、Dwb=2M）。具体的には、周囲雑音情報帯域拡張部３
４は、周囲雑音のパワースペクトル|N[f,w]|2 (w=0,1,…M-1)を入力として、周囲雑音の
パワースペクトル|N[f,w]|2について入力信号x[n]には存在して集音信号z[n]には存在し
ない周波数帯域成分のパワースペクトルを周波数帯域拡張によって生成して、その帯域拡張されたパワースペクトルに対してマスキング閾値を求め、その結果である広帯域マスキング閾値N_wb_th[f,w] (w=0,1,…2M-1)を出力する。
【００６０】
辞書格納部３４２は、Dnb次の狭帯域特徴量データとDwb次の広帯域特徴量データとの対応をモデル化して事前に学習されたサイズＱ（ここではQ=64）のコードブックの辞書λ2q={μxq,μyq}（q=1,…,Q）を格納している。なお、μxqはq番目のコードブックにおける
狭帯域特徴量データのセントロイドベクトル、μyqはq番目のコードブックにおける広帯
域特徴量データのセントロイドベクトルを表している。なお、コードブックのコードベクトルの次数は、狭帯域特徴量データのセントロイドベクトルμxqと広帯域特徴量データのセントロイドベクトルμyqの成分の和であるDnb＋Dwbである。
【００６１】
辞書格納部３４２における事前の辞書λ2qの学習生成方法について、フローチャートを図１０に示し、説明する。以下の説明では、前述の辞書λ1qの学習生成方法と同じ処理については同じ番号を付番し、説明を簡明にするために必要に応じて重複する説明を省略する。
【００６２】
コードブックの辞書生成に用いる信号は、入力信号x[n]と同様でサンプリング周波数fs’[Hz]でfs_wb_low[Hz]からfs_wb_high[Hz]までに帯域制限された広帯域な信号を別途事
前に集音した信号群である。この信号群は、多数の様々な環境、様々な音量であることが望ましい。以下では、コードブックの辞書生成に用いる広帯域信号の信号群をまとめて広帯域信号データwb[n]と表記する。また、ｎは時刻（サンプル）を表す。
【００６３】
まず、広帯域信号データwb[n]を入力として、サンプリング周波数fs[Hz]にダウンサン
プリングし狭帯域信号データnb[n]を得る（ステップＳ１０１）。そして、狭帯域信号デ
ータnb[n]から狭帯域信号情報を表す特徴量データである狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）を抽出する（ステップＳ２０２）。このステップＳ２０２では、狭帯域信号
データnb[n]のパワースペクトル（M次）を得て（ステップＳ１０２１）、狭帯域信号データnb[n]のパワー値を得て（ステップＳ１０２２）、これらのパワースペクトルとパワー
値から狭帯域信号データnb[n]の正規化されたパワースペクトルを得て（ステップＳ１０
２３）、これを次数Dnb（=M）の狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）とすること
によって狭帯域特徴量データの抽出を行う。
【００６４】
一方、広帯域信号データwb[n]から広帯域信号情報を表す特徴量データである広帯域特
徴量データPwb[f,d]（d=1,…,Dwb）を抽出する（ステップＳ２０３）。このステップＳ２０３では、広帯域信号データwb[n]のパワースペクトルを得て（ステップＳ１０３１）、
広帯域信号データwb[n]から広帯域信号データwb[n]のパワー値をフレーム単位で得て（ステップＳ２０３２）、これらのパワースペクトルとパワー値から広帯域信号データwb[n]
の正規化されたパワースペクトルをフレーム単位で得て（ステップＳ２０３３）、これを次数Dwb（=2M）の広帯域特徴量データPwb[f,d]（d=1,…,Dwb）とすることによって広帯域特徴量データの抽出を行う。
【００６５】
次に、狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）と広帯域特徴量データPwb[f,d]（d=1,…,Dwb）を連結して、次数Dnb+Dwbの連結特徴量データP[f,d]（d=1,…,Dnb+Dwb）を生
成する（ステップＳ１０４）。
【００６６】
上記連結特徴量データＰ[f,d]からサイズＱ（ここではQ=64）のコードブックの辞書λ2q={μxq,μyq}（q=1,…,Q）をｋ−ｍｅａｎｓアルゴリズムやＬＢＧアルゴリズムなどに
よるクラスタリング手法を用いて生成する（ステップＳ２０５）。ステップＳ２０５では、まず狭帯域セントロイドベクトルμx1を狭帯域特徴量データの全部の平均とし、広帯域セントロイドベクトルμy1を広帯域特徴量データの全部の平均としてサイズQ=1の初期コ
ードブックを生成する（ステップＳ２０５１）。そして、コードブックのサイズQが所定
数（ここでは64）に達したかどうかを判定する（ステップＳ２０５２）。コードブックのサイズQが所定数に達していなければ、コードブックλ2qの各コードベクトルにおける狭
帯域セントロイドベクトルμxqと広帯域セントロイドベクトルμyqをわずかにずらして別のコードベクトルを生成することでコードブックのサイズQを２倍に増やす処理を行う（
ステップＳ２０５３）。そして、次数Dnb+Dwbの連結特徴量データP[f,d]について、コー
ドブックλ2qの各コードベクトルにおける狭帯域セントロイドベクトルμxqとの所定の距離尺度（例えばユークリッド距離やマハラノビス距離）が最小となるコードベクトルを求めて、連結特徴量データP[f,d]をその該当するコードベクトルに割り当てる。その後、コードブックλ2qのコードベクトルごとに割り当てられた連結特徴量データP[f,d]を用いて、コードベクトルごとに新しい狭帯域セントロイドベクトルμxqと広帯域セントロイドベクトルμyqを求めて、コードブックλ2qを更新する（ステップＳ２０５４）。コードブックのサイズQが所定数に達していれば、そのコードブックλ2q={μxq,μyq}（q=1,…,Q）
を出力する。
【００６７】
広帯域パワースペクトル算出部３４３は、パワー正規化部３２１から出力された周囲雑音の正規化したパワースペクトル|Nn[f,w]|2 (w=0,1,…M-1)をDnb次の特徴量データとし
て入力し、辞書格納部３４２からコードブックの辞書λ2q={μxq,μyq}（q=1,…,Q）を読み出して、Dnb次の狭帯域特徴量データとDwb次の広帯域特徴量データとの対応から広帯域パワースペクトル|Nw[f,w]|2 (w=0,1,…2M-1)を求める。具体的には、Q個ある狭帯域セントロイドベクトルμxq（q=1,…,Q）から、周囲雑音の正規化したパワースペクトル|Nn[f,w]|2 (w=0,1,…M-1)と所定の距離尺度で一番距離が近いものを求めて、一番距離が近いコードベクトルにおける広帯域セントロイドベクトルμyqを広帯域パワースペクトル|Nw[f,w]|2 (w=0,1,…2M-1)とする。
【００６８】
広帯域マスキング閾値算出部３４４は、広帯域パワースペクトル算出部３４３から出力される広帯域パワースペクトル|Nw[f,w]|2 (w=0,1,…2M-1)を入力として、周波数成分ご
とに周囲雑音のマスキング閾値である広帯域マスキング閾値N_wb_th1[f,w] (w=0,1,…2M-1)を算出する。
【００６９】
一般にマスキング閾値は、spreading functionと呼ばれる関数を信号のパワースペクトルに畳み込むことで算出することができる。すなわち、周囲雑音の広帯域マスキング閾値N_wb_th1[f,w] (w=0,1,…2M-1)は、spreading functionを関数sprdngf()として、式（１
３）の式で算出される。周囲雑音の広帯域パワースペクトル|Nw[f,w]|2が広帯域マスキング閾値N_wb_th1[f,w]以下であるならば、周波数ビンω以外の周波数帯域の周囲雑音の広
帯域パワースペクトルによってマスクされる。図１１に、横軸を周波数[Hz]、縦軸をパワー[dB]として、屋外など様々な環境で採取した周囲雑音の広帯域マスキング閾値の例を示す。
【数１３】
2010224321.tif 000015
ここで、bark[w]は周波数ビンωをバーク尺度に変換したバーク値を表し、spreading functionでは、バーク尺度bark[w]に適宜変換する。バーク尺度は、聴覚の分解能を考慮して、低域ほど細かく、高域ほど粗く設定された尺度である。
【００７０】
ここでは、spreading functionを関数sprdngf()として、ISO/IEC13818-7で定義されて
いる方式を用いるとする。spreading functionは、例えばITU-R1387、3GPP TS 26.403と
いった文献で説明されている他の方式を用いても良い。なお、バーク尺度でなくても、メル尺度、ＥＲＢ尺度など人間の音の高さの知覚特性や聴覚フィルタから得られた尺度を用いたspreading functionを適宜用いても構わない。
【００７１】
パワー制御部３４５は、パワー正規化３２１から出力された狭帯域パワー値Pow_N[f]と広帯域マスキング閾値算出部３４５から出力された広帯域マスキング閾値N_wb_th1[f,w] (w=0,1,…2M-1)を入力として、広帯域マスキング閾値N_wb_th1[f,w]のfs_nb_low[Hz]からfs_nb_high[Hz]におけるパワーを狭帯域パワー値Pow_N[f]と同じになるように、広帯域マスキング閾値N_wb_th1[f,w]を増幅あるいは減衰させることで制御し、このパワー制御さ
れたN_wb_th1[f,w]を広帯域マスキング閾値N_wb_th[f,w]として出力する。
【００７２】
このようにして、周囲雑音情報帯域拡張部３４では、狭帯域信号情報である周囲雑音のパワースペクトル|N[f,w]|2から、広帯域信号情報である広帯域マスキング閾値N_wb_th[f,w]を求める。
【００７３】
図１２に信号特性補正部３５の構成例を示す。信号特性補正部３５は、周波数領域変換部３３１と、パワー算出部３５２と、マスキング閾値算出部３５３と、マスキング判定部３５４と、パワー平滑化部３５５と、補正度合決定部３５６と、補正処理部３３３と、時間領域変換部３３４とを備える。
【００７４】
信号特性補正部３５は、入力信号x[n] (n=0,1,…2N-1)と広帯域マスキング閾値N_wb_th[f,w]を入力とし、入力信号x[n]が集音信号に含まれる周囲雑音に埋もれてしまわないよ
う明瞭化する信号補正処理を行い、その補正後の出力信号y[n] (n=0,1,…2N-1)を出力す
る。
【００７５】
パワー算出部３５２は、周波数領域変換部３３１から出力された入力信号x[n]の周波数スペクトルX[f,w] (w=0,1,…2M-1)における実部と虚部の２乗和であるパワースペクトル|X[f,w]|2  (w=0,1,…2M-1)を算出して出力する。
【００７６】
マスキング閾値算出部３５３は、パワー算出部３５２から出力された入力信号x[n]のパワースペクトル|X[f,w]|2  (w=0,1,…2M-1)を入力として、spreading functionを関数sprdngf()として、式（１４）の式で入力信号x[n]の広帯域マスキング閾値X_th[f,w] (w=0,1,…2M-1)を算出し出力する。広帯域マスキング閾値X_th[f,w]は、入力信号x[n]のパワー
スペクトル|X[f,w]|2 が入力信号x[n]の広帯域マスキング閾値X_th[f,w]以下であるなら
ば、周波数ビンω以外の周波数帯域の入力信号x[n]のパワースペクトル|X[f,w]|2 によってマスクされることを表す。
【数１４】
2010224321.tif 000016
マスキング判定部３５４は、パワー算出部３５２から出力されたパワースペクトル|X[f,w]|2  (w=0,1,…2M-1)とマスキング閾値算出部３５３から出力された広帯域マスキング
閾値X_th[f,w]とを入力とし、周波数帯域ごとに入力信号x[n]自身によってマスクされる
か否かを表すマスキング判定情報X_flag[f,w] (w=0,1,…2M-1)を出力する。具体的には、パワースペクトル|X[f,w]|2 と広帯域マスキング閾値X_th[f,w]の大小比較を行い、パワ
ースペクトル|X[f,w]|2 が広帯域マスキング閾値X_th[f,w]以上ならば、その周波数成分
は入力信号x[n]中の他の周波数成分にマスクされないとしてX_flag[f,w]=0とする。また
、パワースペクトル|X[f,w]|2 が広帯域マスキング閾値X_th[f,w]未満ならば、その周波
数成分は入力信号x[n]中の他の周波数成分にマスクされるとしてX_flag[f,w]=1とする
パワー平滑化部３５５は、パワー算出部３５２から出力されたパワースペクトル|X[f,w]|2  (w=0,1,…2M-1)とマスキング判定部３５４から出力されたマスキング判定情報X_flag[f,w]とを入力として、パワースペクトル|X[f,w]|2 を式（１５）の式による三角窓による移動平均によって平滑化して、平滑化されたパワースペクトル|XS[f,w]|2 を出力する
。なお、Kは平滑化を計算する範囲であり、αX[j]は、ｊが０に近いほど係数が大きくな
るようなスムージング係数である。例えば、K=３で、αX[j]は［０．１、０．２、０．４、０．８、０．４、０．２、０．１］とする。
【数１５】
2010224321.tif 000017
補正度合決定部３５６は、パワー平滑化部３５５から出力された平滑化されたパワースペクトル|XS[f,w]|2 (w=0,1,…2M-1)とマスキング判定部３５４から出力されたマスキン
グ判定情報X_flag[f,w] (w=0,1,…2M-1)と周囲雑音情報帯域拡張部３２から出力されたN_wb_th[f,w] (w=0,1,…2M-1)とを入力として、補正ゲインG[f,w] (w=0,1,…2M-1)を算出して出力する。補正ゲインG[f,w]の具体的な算出は、まずマスキング判定情報X_flag[f,w]
により入力信号x[n]中の他の周波数成分にマスクされる（X_flag[f,w]=1）と判定された
周波数帯域であれば、G[f,w]=1とし補正による増幅も減衰も行わないようにする。そして、マスキング判定情報X_flag[f,w]により入力信号x[n]中の他の周波数成分にマスクされ
ない（X_flag[f,w]=0）と判定された周波数帯域について、パワースペクトル|X[f,w]|2 
と周囲雑音の広帯域マスキング閾値N_wb_th[f,w]との大小比較を行う。ここで、パワースペクトル|X[f,w]|2 が周囲雑音の広帯域マスキング閾値N_wb_th[f,w]以上ならば、その周波数成分は集音信号z[n]中の他の周波数成分にマスクされないのでG[f,w]=1とし、補正による増幅を行わないようにする。一方で、パワースペクトル|X[f,w]|2 が周囲雑音の広帯域マスキング閾値N_wb_th[f,w]未満ならば、集音信号z[n]中の周囲雑音が少なければ知覚できるにも関わらず、周囲雑音があるためにマスクされていると判断し、式（１６）の式のように補正ゲインG[f,w]を周囲雑音の広帯域マスキング閾値N_wb_th[f,w]と平滑化されたパワースペクトル|XS[f,w]|2との比に基づいて算出する。なお、関数Ｆは、平滑されたパワースペクトル|XS[f,w]|2のスペクトル傾斜を周囲雑音の広帯域マスキング閾値N_wb_th[f,w]の形状と平行に近くなるように増幅するような関数である。ここで、α、βは正の定数であり、γは正負いずれかの定数である。これらの定数は、入力信号x[n]の増幅度合いを調整するために用いられる。
【数１６】
2010224321.tif 000018
【数１７】
2010224321.tif 000019
補正度合決定部３５６において、このように求めた補正ゲインG[f,w]をさらに式（２２）の式による三角窓による移動平均によって平滑化して、平滑化された補正ゲインGS[f,w]を用いてもよい。なお、Kは平滑化を計算する範囲であり、αX[j]は、ｊが０に近いほど係数が大きくなるようなスムージング係数である。例えば、K=３で、αG[j]は［０．１、０．２、０．４、０．８、０．４、０．２、０．１］とする。
【数１８】
2010224321.tif 000020
以上のように、再生される入力信号と集音信号で、信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりしていても、集音信号の周波数特性であるパワースペクトルについて入力信号の周波数帯域を加味して帯域拡張して推定することで、集音信号の周波数特性が高精度に求まり、入力信号の明瞭度を向上させることができる。
【００７７】
なお、図７（ｂ）に示す音声帯域拡張通話装置に本変形例を適用する場合は、信号帯域拡張処理部１２において事前に設定した周波数f_limit（f_limitは500〜1200[Hz]程度で
、例えばf_limit=1000[Hz]とする）以下の低い周波数帯域が拡張されるとき、つまりfs_wb_low < fs_nb_lowかつfs_wb_low < f_limitであるときは、信号特性補正部３５でf_limit以下の周波数帯域について信号補正処理をしないようにする。低域（f_limit以下の周波数）においては、集音する環境やノイズ成分の種類によって、周囲雑音のバラツキが大きいため、このようにすることで、信号帯域拡張処理部１２において拡張した低い周波数帯域での周囲雑音のバラツキによって信号補正処理が不安定になることを防止できる。
【００７８】
（信号処理部の変形例２）
本変形例では、図８に示す信号処理部３０にて用いる狭帯域信号情報を周囲雑音のパワースペクトルとし、広帯域信号情報を周囲雑音の広帯域パワースペクトル（周囲雑音を広帯域の信号に拡張した場合のパワースペクトル）とした場合を例にして説明する。この場合、周囲雑音情報帯域拡張部３４では、狭帯域信号情報である周囲雑音のパワースペクトルを入力として、狭帯域特徴量データとして周囲雑音の正規化されたパワースペクトルを算出し、広帯域特徴量データである周囲雑音の正規化された広帯域パワースペクトルを事前にモデル化された狭帯域特徴量データと広帯域特徴量データとの対応を用いて算出し、この広帯域特徴量データから広帯域信号情報である周囲雑音の広帯域パワースペクトルを生成するようにする。なお、狭帯域特徴量データと広帯域特徴量データとの対応のモデル化には、図５に示すGMMを利用する手法を用いる。これによれば、再生される入力信号と
集音信号で、信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりしていても、集音信号の周波数特性であるパワースペクトルについて入力信号の周波数帯域を加味して帯域拡張して推定することで、集音信号の周波数特性が高精度に求まり、入力信号の明瞭度を向上させることができる。
【００７９】
（信号処理部の変形例３）
次に、信号処理部にて用いる狭帯域信号情報は周囲雑音のパワースペクトル、広帯域信号情報は周囲雑音を広帯域の信号に拡張した場合のマスキング閾値（広帯域マスキング閾値）である場合を例にして説明する。
【００８０】
図１３は、その構成を示すものである。信号処理部３００では、信号処理部３０で用いていた周囲雑音情報帯域拡張部３４に代わって、周囲雑音情報帯域拡張部３６を用いた構成になっている。
【００８１】
図１４に周囲雑音情報帯域拡張部３６の構成例を示す。周囲雑音情報帯域拡張部３６は、パワー正規化部３２１と、狭帯域マスキング閾値算出部３６２と、帯域制御部３６３と、辞書格納部３６４と、広帯域マスキング閾値算出部３６５と、閾値補正部３６６と、パワー制御部３４５とを備える。
【００８２】
周囲雑音情報帯域拡張部３６は、周囲雑音情報帯域拡張部３４と同様に、集音信号z[n]の周波数帯域成分における情報（狭帯域信号情報）を入力として、入力信号x[n]には存在して集音信号z[n]には存在しない周波数帯域成分を含めた情報（広帯域信号情報）を生成する。つまり、周囲雑音情報帯域拡張部３６では、狭帯域信号情報から狭帯域特徴量データを算出し、狭帯域特徴量データと広帯域特徴量データとの対応を事前にモデル化しておき、このモデルと取得した狭帯域特徴量データとの対応を用いて広帯域特徴量データを算出し、広帯域特徴量データから広帯域信号情報を生成する。この際、周囲雑音情報帯域拡張部３６では、狭帯域特徴量データと広帯域特徴量データとの対応のモデル化に、ベクトル量子化によるコードブックを利用する手法を用いる。ここでは、周囲雑音の帯域制御された狭帯域マスキング閾値N_th[f,w] (w=0,1,…MC-1)のDnb次の狭帯域特徴量データとし
て用い、周囲雑音の広帯域マスキング閾値N_wb_th1[f,w](w=0,1,…2M-1)のDwb次の広帯域特徴量データとして用いる（Dnb=MC、Dwb=2M）。具体的には、周囲雑音情報帯域拡張部３６は、周囲雑音のパワースペクトル|N[f,w]|2 (w=0,1,…M-1)を入力として、周囲雑音の
マスキング閾値を求め、このマスキング閾値を帯域制限し、帯域制限されたマスキング閾値について入力信号x[n]には存在して集音信号z[n]には存在しない周波数帯域成分を周波数帯域拡張して生成して、この帯域拡張されたマスキング閾値である広帯域マスキング閾値N_wb_th[f,w] (w=0,1,…2M-1)を出力する。
【００８３】
狭帯域マスキング閾値算出部３６２は、パワー正規化部３２１から出力される周囲雑音の正規化したパワースペクトル|Nn[f,w]|2(w=0,1,…M-1)を入力として、周波数成分ごと
に周囲雑音のマスキング閾値である狭帯域マスキング閾値N_th1[f,w] (w=0,1,…M-1)を算出する。前述した広帯域マスキング閾値算出部３４４と同様にして、データ長である2MをMで置き換え、周囲雑音の狭帯域マスキング閾値N_th1[f,w] (w=0,1,…M-1)は、spreading functionを関数sprdngf()として、式（１９）の式で算出される。狭帯域マスキング閾値N_th1[f,w]は、周囲雑音の正規化したパワースペクトル|Nn[f,w]|2が狭帯域マスキング閾値N_th1[f,w]以下であるならば、周波数ビンω以外の周波数帯域の周囲雑音の正規化したパワースペクトルによってマスクされることを示す。
【数１９】
2010224321.tif 000021
帯域制御部３６３は、狭帯域マスキング閾値算出部３６２から出力された周囲雑音の狭帯域マスキング閾値N_th1[f,w] (w=0,1,…M-1)を入力として、帯域制御する下限周波数limit_low[Hz]から帯域制御する上限周波数limit_high[Hz]までの周波数帯域の信号情報の
みを用いるように制御し、帯域制御された狭帯域マスキング閾値であるN_th[f,w]を出力
する。ただし、fs_nb_low ≦ limit_low < limit_high ≦ fs_nb_high < fs/2とする。例えばlimit_low=1000[Hz]、limit_high=3400[Hz]とするとき、これらの周波数帯域を式（
２４）で周波数ビンωに変換して考慮すると、狭帯域マスキング閾値N_th1[f,w] (w=0,1,…M-1)のうちw=32,33,…108のみを用いるようにする。ＭCをN_th[f,w]の配列の個数とし
て、帯域制御された狭帯域マスキング閾値N_th[f,w] (w=0,1,…MC-1)は、狭帯域マスキング閾値N_th1[f,w] (w=32,…108)そのものを代入する。この場合はMC=108-32+1=77である
。
【００８４】
図１１にも示す通り、低域においては、集音する環境やノイズ成分の種類によって、周囲雑音のマスキング閾値の分散・バラツキが大きいことが分かる。周囲雑音の主要な成分はノイズ成分であるため、狭帯域マスキング閾値N_th1[f,w]も低域においては、分散・バラツキが大きくなる。そこで、狭帯域特徴量データと広帯域特徴量データとの対応のモデル化をベクトル量子化によるコードブックを利用する手法を用いて、広帯域マスキング閾値を高精度に求めるために、分散・バラツキが大きい低域を用いないように帯域制御する。つまり、ここで、帯域制御する下限周波数limit_low[Hz]は、狭帯域マスキング閾値の
分散・バラツキが所定の値よりも小さいような周波数帯域の下限に設定することが望ましい。こうすることによって、広帯域マスキング閾値を高精度に求めることができ、入力信号の明瞭度を向上させることができる。
【００８５】
また、マスキング閾値は、その周波数帯域のパワースペクトルだけではなくて周囲の周波数帯域のパワースペクトルを加味して算出される。そのため、マスキング閾値を求める元々の信号の帯域制限されている周波数帯域付近では、正確にマスキング閾値が算出できない。つまり、帯域制御する上限周波数limit_high[Hz]は、帯域制限を加味してもマスキング閾値が正確に求まる周波数帯域の上限に設定することが望ましい。こうすることによって、広帯域マスキング閾値を高精度に求めることができ、入力信号の明瞭度を向上させることができる。
【００８６】
辞書格納部３６４は、Dnb次の狭帯域特徴量データとDwb次の広帯域特徴量データとの対応をモデル化して事前に学習されたサイズＱ（ここではQ=64）のコードブックの辞書λ3q={μxq,μyq}（q=1,…,Q）を格納している。なお、μxqはq番目のコードブックにおける
狭帯域特徴量データのセントロイドベクトル、μyqはq番目のコードブックにおける広帯
域特徴量データのセントロイドベクトルを表している。なお、コードブックのコードベクトルの次数は、狭帯域信号情報のセントロイドベクトルμxqと広帯域信号情報のセントロイドベクトルμyqの成分の和であるDnb＋Dwbである。
【００８７】
辞書格納部３６４における事前の辞書λ3qの学習生成方法の一手法について、フローチャートを図１５に示し、説明する。以下の説明では、上述した変形例１における辞書λ2qの学習生成方法と同じ処理については同じ番号を付番し、説明を簡明にするために必要に応じて重複する説明を省略する。
【００８８】
まず、広帯域信号データwb[n]を入力として、サンプリング周波数fs[Hz]にダウンサン
プリングし狭帯域信号データnb[n]を得る（ステップＳ１０１）。そして、狭帯域信号デ
ータnb[n]から狭帯域信号情報を表す特徴量データである狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）を抽出する（ステップＳ２０２）。このステップＳ２０２では、狭帯域信号
データnb[n]のパワースペクトル（M次）を得て（ステップＳ１０２１）、狭帯域信号データnb[n]のパワー値を得て（ステップＳ１０２２）、これらのパワースペクトルとパワー
値から狭帯域信号データnb[n]の正規化されたパワースペクトルを得て（ステップＳ１０
２３）、式（２３）と同様にして狭帯域信号データnb[n]のマスキング閾値を算出する（
ステップＳ３０２４）。そして、狭帯域信号データnb[n]のマスキング閾値に対して、帯
域制御部３６３での処理と同様に帯域制御する（ステップＳ３０２５）。これを次数Dnb
（=MC）の狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）とすることによって狭帯域特徴量
データの抽出を行う。
【００８９】
一方、広帯域信号データwb[n]から広帯域信号情報を表す特徴量データである広帯域特
徴量データPwb[f,d]（d=1,…,Dwb）を抽出する（ステップＳ３０３）。このステップＳ３０３では、広帯域信号データwb[n]のパワースペクトル（2M次）を得て（ステップＳ１０
３１）、広帯域信号データwb[n]から広帯域信号データwb[n]のパワー値を得て（ステップＳ２０３２）、これらのパワースペクトルとパワー値から広帯域信号データwb[n]の正規
化されたパワースペクトルをフレーム単位で得て（ステップＳ２０３３）、式（２３）の次数をMから2Mにして同様にして広帯域信号データwb[n]のマスキング閾値を算出する（ステップＳ３０３４）。これを次数Dwb（=2M）の広帯域特徴量データPwb[f,d]（d=1,…,Dwb）とすることによって広帯域特徴量データの抽出を行う。
【００９０】
次に、狭帯域特徴量データPnb[f,d]（d=1,…,Dnb）と広帯域特徴量データPwb[f,d]（d=1,…,Dwb）を連結して、次数Dnb+Dwbの連結特徴量データP[f,d]（d=1,…,Dnb+Dwb）を生
成する（ステップＳ１０４）。
【００９１】
そして、連結特徴量データＰ[f,d]からコードブックの各コードベクトルにおける狭帯
域セントロイドベクトルμxqと広帯域セントロイドベクトルμyqを求め、サイズＱ（ここではQ=64）のコードブックをｋ−ｍｅａｎｓアルゴリズムやＬＢＧアルゴリズムなどによるクラスタリング手法を用いて生成する（ステップＳ２０５）。コードブックの各コードベクトルにおける広帯域セントロイドベクトルμyqである広帯域信号データwb[n]のマス
キング閾値を近似多項式係数で表現して、近似多項式係数を広帯域セントロイドベクトルμ’yqとして辞書に格納して、辞書λ3q={μxq,μ’yq}（q=1,…,Q）を生成する（ステップＳ３０７）。近似多項式係数ｍp（p=0,…,P）とはここでは、縦軸をパワー値X[dB]、横軸を周波数Y[Hz]として、式（２０）のようにマスキング閾値を所定の次数（ここではPとし、例えばP=6とする）の多項式で近似した、その多項式の係数のことであり、これ以降
そのように呼ぶ。
【数２０】
2010224321.tif 000022
このように、マスキング閾値を近似多項式係数で表現して辞書として格納しておくことで、マスキング閾値を辞書として格納しておくよりも、辞書の格納に掛かるメモリ量を削減することができ、辞書の配列の数を小さくなるため辞書の利用時の処理量を削減することができる。
【００９２】
広帯域マスキング閾値算出部３６５は、帯域制御部３６３から出力された帯域制御された狭帯域マスキング閾値N_th[f,w] (w=0,1,…MC-1)をDnb次の特徴量データとして入力し
、辞書格納部３６４からコードブックの辞書λ3q={μxq,μ’yq}（q=1,…,Q）を読み出して、Dnb次の狭帯域特徴量データとDwb次の広帯域特徴量データとの対応から周囲雑音の広帯域マスキング閾値N_wb_th1[f,w](w=0,1,…2M-1)を求める。具体的には、Q個ある狭帯域セントロイドベクトルμxq（q=1,…,Q）から、帯域制御された狭帯域マスキング閾値N_th[f,w] (w=0,1,…MC-1)と所定の距離尺度で一番距離が近いものを求めて、一番距離が近いコードベクトルにおける広帯域セントロイドベクトルμ’yqをそのまま広帯域マスキング閾値の近似多項式係数として設定し、式（２０）と同様にして広帯域マスキング閾値N_wb_th1[f,w](w=0,1,…2M-1)を算出する。
【００９３】
閾値補正部３６６は、狭帯域マスキング閾値算出部３６２から出力された周囲雑音の狭帯域マスキング閾値N_th1[f,w] (w=0,1,…M-1)と広帯域マスキング閾値算出部３６５から出力された周囲雑音の広帯域マスキング閾値N_wb_th1[f,w] (w=0,1,…2M-1)を入力として、狭帯域と広帯域における境界帯域付近での不連続性あるいは微分不連続性を解消するように補正し、その補正された広帯域マスキング閾値N_wb_th2[f,w] (w=0,1,…2M-1)を出力する。図１６（ａ）に、境界帯域fs/2[Hz]前後の周波数において、狭帯域マスキング閾値N_th[f,w]と広帯域マスキング閾値N_wb_th1[f,w]とに不連続性が生じ、それを解消するように補正された広帯域マスキング閾値N_wb_th2[f,w]の例を示す。図１６（ｂ）に、境界
帯域fs/2[Hz]前後の周波数において、狭帯域マスキング閾値N_th[f,w]と広帯域マスキン
グ閾値N_wb_th1[f,w]とに不連続性と微分不連続性の両方が生じ、それを解消するように
補正された広帯域マスキング閾値N_wb_th2[f,w]の例を示す。両図共に、実線は狭帯域マ
スキング閾値N_th[f,w]を、破線は広帯域マスキング閾値N_wb_th2[f,w]を、太実線は補正された広帯域マスキング閾値N_wb_th2[f,w]における補正箇所を表す。ただし、adjust_low[Hz] < fs/2 < adjust_high[Hz]とする。ここで、adjust_lowは周波数ビンωL−１に対応する周波数以上で周波数ビンωLに対応する周波数未満であり、adjust_highは周波数ビンωHに対応する周波数以上で周波数ビンωH＋１に対応する周波数未満であるとする。例えばfs=8000[Hz]であるとき、adjust_low=3600[Hz]、adjust_high=4400[Hz]とする。具体的には、少なくとも境界帯域fs/2[Hz]前後の周波数において不連続あるいは微分不連続が検出された場合に、adjust_low[Hz]以上かつadjust_high[Hz]以下であるような境界帯域
付近について、周波数ビンωL、ωL＋１…、ωL＋SとωH、ωH−１…、ωH−Sにおける広帯域マスキング閾値N_wb_th1[f,w]を用いて、周波数ビンωL＋S＋１からωH−S−１まで
の広帯域マスキング閾値を（2S-1）次関数で模擬し、スプライン補間を行うことで、補正された広帯域マスキング閾値N_wb_th2[f,w]を求める。ここで、狭帯域マスキング閾値N_th1[f,M-1]と広帯域マスキング閾値N_wb_th1[f,M]との中点を通過するように模擬する関数を設定してスプライン補間を行ってもよい。
【００９４】
このように閾値補正部３６６において広帯域マスキング閾値を補正することで、広帯域マスキング閾値における不連続性あるいは微分不連続性が解消され、信号補正においても周波数方向の不連続性が無くなって違和感のない自然な信号補正にすることができ、高い明瞭感を得ることができる。
【００９５】
以上のように、再生される入力信号と集音信号で、信号成分が存在する周波数帯域が異なっていたり、サンプリング周波数が異なっていたりしていても、集音信号のマスキング閾値について入力信号の周波数帯域を加味して帯域拡張して推定することで、集音信号のマスキング閾値が高精度に求まり、入力信号の明瞭度を向上させることができる。
【００９６】
（信号処理部の変形例４）
信号処理部３００の辞書格納部３６４における事前の辞書λ3qの学習生成方法の他の手法について、フローチャートを図１７に示し、説明する。ここでは、狭帯域信号データnb[n]を生成しないで広帯域信号データwb[n]のみから辞書λ3qを学習生成する方法について説明する。以下の説明では、上述した変形例２における辞書λ3qの学習生成方法と同じ処理については同じ番号を付番し、説明を簡明にするために必要に応じて重複する説明を省略する。
【００９７】
まず、ステップＳ３０３で広帯域信号データwb[n]から広帯域信号情報を表す特徴量デ
ータ（ここではマスキング閾値）である広帯域特徴量データPwb[f,d]（d=1,…,Dwb）を抽出する。この広帯域特徴量データPwb[f,d]（d=1,…,Dwb）のみを用いて、ステップＳ２０５でサイズＱのコードブックを作成する。そして、コードブックの各コードベクトルにおける広帯域セントロイドベクトルμyqである広帯域信号データwb[n]の広帯域マスキング
閾値に対して、帯域制御する下限周波数limit_low[Hz]から帯域制御する上限周波数limit_high[Hz]までの周波数帯域の広帯域マスキング閾値のみを用いるように制御する（ステ
ップＳ３０２５）。これにより狭帯域に帯域制御された狭帯域マスキング閾値が求まり、これをコードブックの各コードベクトルにおける狭帯域セントロイドベクトルμxq（q=1,…,Q）とする（ステップＳ３０６）。その後、ステップＳ３０７で広帯域信号データwb[n]のマスキング閾値の近似多項式係数である広帯域セントロイドベクトルμ’yqと併せて
辞書に格納して、辞書λ3q={μxq,μ’yq}を生成する。
【００９８】
狭帯域特徴量データを併用してクラスタリングする図１５における手法では、狭帯域特徴量データに狭帯域と広帯域における境界帯域付近で誤差を含む。このように、広帯域特徴量データのみを用いてクラスタリングして、広帯域セントロイドベクトルを帯域制限して狭帯域セントロイドベクトルを求めることで、理想的なデータである広帯域特徴量データのみを用いてクラスタリングするため、図１５における手法よりも、高精度にクラスタリングを行うことができる。
【００９９】
（信号処理部の変形例５）
信号処理部３００の辞書格納部３６４における事前の辞書λ3qの学習生成方法の他の手法について、フローチャートを図１８に示し、説明する。以下の説明では、上述した変形例２における辞書λ3qの学習生成方法と同じ処理については同じ番号を付番し、説明を簡明にするために必要に応じて重複する説明を省略する。
【０１００】
ステップＳ２０５でサイズＱのコードブックを作成した後、コードブックの各コードベクトルにおける狭帯域セントロイドベクトルμxqである狭帯域信号データnb[n]のマスキ
ング閾値を式（２０）のように近似多項式で表現して、近似多項式係数を狭帯域セントロイドベクトルμ’xq（q=1,…,Q）とする（ステップＳ３０６Ａ）。その後、ステップＳ３０７で広帯域信号データwb[n]のマスキング閾値の近似多項式係数である広帯域セントロ
イドベクトルμ’yqと併せて辞書に格納して、辞書λ3q={μ’xq,μ’yq}を生成する。
【０１０１】
一方で、この手法においては、広帯域マスキング閾値算出部３６５では、帯域制御部３６３から出力された帯域制御された狭帯域マスキング閾値N_th[f,w] (w=0,1,…MC-1)をDnb次の特徴量データとして入力し、辞書格納部３６４からコードブックの辞書λ3q={μ’xq,μ’yq}（q=1,…,Q）を読み出して、Dnb次の狭帯域特徴量データとDwb次の広帯域特徴
量データとの対応から周囲雑音の広帯域マスキング閾値N_wb_th1[f,w](w=0,1,…2M-1)を
求めるようにする。具体的には、Q個ある狭帯域セントロイドベクトルμ’xq（q=1,…,Q
）の近似多項式から、帯域制御された狭帯域マスキング閾値N_th[f,w] (w=0,1,…MC-1)と所定の距離尺度で一番距離が近いものを近似多項式に代入していくことで求めて、一番距離が近いコードベクトルにおける広帯域セントロイドベクトルμ’yqをそのまま広帯域マスキング閾値の近似多項式係数として設定し、式（２０）と同様にして広帯域マスキング閾値N_wb_th1[f,w](w=0,1,…2M-1)を算出する。
【０１０２】
このように、狭帯域マスキング閾値も近似多項式係数で表現して辞書として格納しておくことで、マスキング閾値を辞書として格納しておくよりも、図１５における手法と比較しても、辞書の格納に掛かるメモリ量を削減することができ、辞書の配列の数を小さくなるため辞書の利用時の処理量を削減することができる。
【０１０３】
（第２の実施例）
図１９（ａ）は、本発明の第２の実施形態に係わる通信装置の構成を示すものである。
【０１０４】
この図に示す通信装置は、例えば携帯電話などの無線通信装置の受信系を示すものであって、無線通信部１と、デコーダ２と、信号処理部３Ａと、ディジタル・アナログ（D/A）
変換器４と、スピーカ５と、マイク６と、アナログ・ディジタル（A/D）変換器７と、ダ
ウンサンプリング部８と、エコー抑圧処理部９と、エンコーダ１０とを備えている。
【０１０５】
なお、第１の実施例と同様に、本発明は、図１９（ａ）のような通信装置だけでなく、図１９（ｂ）に示すディジタルオーディオプレイヤに適用することも可能である。また、図１９（ｃ）に示す音声帯域拡張通話装置に適用することも可能である。
【０１０６】
次に、信号処理部３Ａについて説明する。図２０は、その構成を示すものである。信号処理部３Ａは、第１の実施例にて説明した信号処理部３に周囲雑音抑圧処理部３７を追加して構成される。以下の説明では、上述した実施例と同じ構成については同じ番号を付番し、必要に応じて重複する説明を省略する。
【０１０７】
図２１に周囲雑音抑圧処理部３７の構成例を示す。周囲雑音抑圧処理部３７は、
抑圧ゲイン算出部３７１と、スペクトル抑圧部３７２と、パワー算出部３７３と、時間領域変換部３７４とを備える。
【０１０８】
周囲雑音抑圧処理部３７は、周囲雑音推定部３１から出力される周囲雑音のパワースペクトルと集音信号z[n]のパワースペクトルと集音信号z[n]の周波数スペクトルを用いて、集音信号z[n]に含まれる周囲雑音であるノイズ成分を抑圧して、周囲雑音であるノイズ成分が抑圧された信号s[n]をエンコーダ１０に出力する。エンコーダ１０では周囲雑音抑圧処理部３７から出力された信号s[n]を符号化して無線通信部１に出力する。
【０１０９】
抑圧ゲイン算出部３７１は、パワー算出部３１２から出力される集音信号z[n]のパワースペクトル|Z[f,w]|2  (w=0,1,…M-1)と、周波数スペクトル更新部３１４から出力される周囲雑音のパワースペクトル|N[f,w]|2  (w=0,1,…M-1)と、パワー算出部３７３から出力される１フレーム前の抑圧処理された信号のパワースペクトル|S[f-1,w]|2  (w=0,1,…M-1)とを用いて、各周波数帯域の抑圧ゲインG[f,w]  (w=0,1,…M-1)を出力する。例えば、
抑圧ゲインG[f,w]の算出は、以下のアルゴリズムまたはそれらの組み合わせによって行う。すなわち、一般のノイズキャンセラであるスペクトル・サブトラクション（Spectral Subtraction）法（S. F. Boll, “Suppression of acoustic noise in speech using spectral subtraction”, IEEE Trans. Acoustics, Speech, and Signal Processing, vol.ASSP-29, pp.113-120 (1979).）、ウィナー・フィルター（Wiener Filter）法（J. S. Lim, A. V. Oppenheim, “Enhancement and bandwidth compression of noisy speech”, Proc. IEEE Vol.67, No.12, pp.1586-1604, Dec.1979.）及び最尤推定（Maximum Likelihood）法（R. J. McAulay, M. L. Malpass, “Speech enhancement using a soft-decision noise suppression filter”, IEEE Trans. on Acoustics, Speech, and Signal Processing, vol.ASSP-28, no.2, pp.137-145, Apr.1980.）などである。ここでは一例としてウィナー・フィルター法を用いて、抑圧ゲインG[f,w]を算出するとする。
【０１１０】
スペクトル抑圧部３７２は、周波数領域変換部３１１から出力された集音信号z[n]の周波数スペクトルZ[f,w] と、抑圧ゲイン算出部３７１から出力された抑圧ゲインG[f,w]と
を入力として、集音信号z[n]の周波数スペクトルZ[f,w]を集音信号z[n]の振幅スペクトル|Z[f,w]|  (w=0,1,…M-1)と位相スペクトルθZ[f,w] (w=0,1,…M-1)に分け、集音信号z[n]の振幅スペクトル|Z[f,w]| に抑圧ゲインG[f,w]を乗じることで周囲雑音であるノイズ成分を抑圧し、その抑圧処理された信号の振幅スペクトル|S[f-1,w]|とし、位相スペクトルθZ[f,w]をそのまま抑圧処理された信号の位相スペクトルθS[f,w]として、抑圧処理された信号の周波数スペクトルS[f,w]  (w=0,1,…M-1)を算出する。
【０１１１】
パワー算出部３７３は、スペクトル抑圧部３７２から出力された抑圧処理された信号の周波数スペクトルS[f,w]  (w=0,1,…M-1)から抑圧処理された信号のパワースペクトル|S[f,w]|2  (w=0,1,…M-1)を算出し出力する。
【０１１２】
時間領域変換部３７４は、スペクトル抑圧部３７２から出力された抑圧処理された信号の周波数スペクトルS[f,w]  (w=0,1,…M-1)を入力として、周波数領域を時間領域に変換
する処理（例えば、ＩＦＦＴ）を施し、周波数領域変換部３１１における窓掛けによるオーバーラップ分を考慮して１フレーム前の抑圧処理された信号s[n]を適宜加算して、抑圧処理された時間領域の信号s[n] (n=0,1,…N-1)を算出する。
【０１１３】
以上のように、周囲雑音推定処理に周囲雑音抑圧処理を併用することで処理量の増加を抑えつつ、入力信号を明瞭化すると同時に、集音信号における周囲雑音成分を抑圧して高音質な集音信号を得ることができる。
【０１１４】
なお、この発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また上記実施形態に開示されている複数の構成要素を適宜組み合わせることによって種々の発明を形成できる。また例えば、実施形態に示される全構成要素からいくつかの構成要素を削除した構成も考えられる。さらに、異なる実施形態に記載した構成要素を適宜組み合わせてもよい。
【０１１５】
例えば、入力信号（あるいは目的信号）のサンプリング周波数は、集音信号（あるいは周囲雑音）のサンプリング周波数の２倍に限定されるものではなく、整数倍でも、非整数倍でもよい。また、入力信号（あるいは目的信号）のサンプリング周波数は、集音信号（あるいは周囲雑音）のサンプリング周波数と等しい上で、入力信号（あるいは目的信号）の周波数帯域制限の範囲と集音信号（あるいは周囲雑音）の周波数帯域制限の範囲が異なっている場合であっても構わない。さらに、入力信号（あるいは目的信号）の周波数帯域制限の範囲が、集音信号（あるいは周囲雑音）の周波数帯域制限の範囲を包含していなくても構わない。さらにまた、入力信号（あるいは目的信号）の周波数帯域制限の範囲は、集音信号（あるいは周囲雑音）の周波数帯域制限の範囲と隣接していなくても構わない。
【０１１６】
また、入力信号がモノラル信号ではなくステレオ信号であったとしても、例えばＬ（左）チャネルとＲ（右）チャネルにそれぞれ上記信号処理部３における信号処理を施したり、和信号（ＬチャネルとＲチャネルの信号の和）と差信号（ＬチャネルからＲチャネルの信号の差）にそれぞれ上記の信号処理を施したりすることで同様の効果が得られる。勿論、マルチチャネル信号であったとしても例えば同様にそれぞれのチャネル信号に対して上記の信号処理を施したりすることで同様の効果が得られる。
【０１１７】
その他、この発明の要旨を逸脱しない範囲で種々の変形を施しても同様に実施可能であることはいうまでもない。
【符号の説明】
【０１１８】
１…無線通信部、２、２Ａ…デコーダ、３、３０、３００、３Ａ…信号処理部、４…ディジタル・アナログ（D/A）変換器、５…スピーカ、６…マイク、７…アナログ・ディジ
タル（A/D）変換器、８…ダウンサンプリング部、９…エコー抑圧処理部、１０…エンコ
ーダ、１１…記憶部、１２…信号帯域拡張処理部、３１…周囲雑音推定部、３２、３４、３６…周囲雑音情報帯域拡張部、３３、３５…信号特性補正部、３７…周囲雑音抑圧処理部、３１１、３３１…周波数領域変換部、３１２、３５２、３７３…パワー算出部、３１３…周囲雑音区間判定部、３１４…周波数スペクトル更新部、３２１…パワー正規化部、３２２、３４２、３６４…辞書格納部、３２３…広帯域パワー算出部、３３２、３５６…補正度合決定部、３３３…補正処理部、３３４、３７４…時間領域変換部、３４３…広帯域パワースペクトル算出部、３４４、３６５…広帯域マスキング閾値算出部、３４５…パワー制御部、３５３…マスキング閾値算出部、３５４…マスキング判定部、３５５…パワー平滑化部、３６２…狭帯域マスキング閾値算出部、３６３…帯域制御部、３６６…閾値補正部、３７１…抑圧ゲイン算出部、３７２…スペクトル抑圧部。
【図１】
2010224321.tif 000023
【図２】
2010224321.tif 000024
【図３】
2010224321.tif 000025
【図４】
2010224321.tif 000026
【図５】
2010224321.tif 000027
【図６】
2010224321.tif 000028
【図７】
2010224321.tif 000029
【図８】
2010224321.tif 000030
【図９】
2010224321.tif 000031
【図１０】
2010224321.tif 000032
【図１１】
2010224321.tif 000033
【図１２】
2010224321.tif 000034
【図１３】
2010224321.tif 000035
【図１４】
2010224321.tif 000036
【図１５】
2010224321.tif 000037
【図１６】
2010224321.tif 000038
【図１７】
2010224321.tif 000039
【図１８】
2010224321.tif 000040
【図１９】
2010224321.tif 000041
【図２０】
2010224321.tif 000042
【図２１】
2010224321.tif 000043