JPH0944186A - 雑音抑制装置 - Google Patents

Abstract

(57)【要約】 【目的】 雑音が混在した音響信号から雑音を抑制する
雑音抑制装置を提供する。 【構成】 信号の時間分割手段１００は入力信号を短時
間セグメントに分割しピッチ抽出手段２００は入力信号
のピッチを抽出する。周波数帯域分割手段４００は短時
間信号を周波数分割しエネルギーまたは実効値の算出手
段６００は周波数分割手段の出力のエネルギーまたは実
効値を算出する。特徴量の算出手段８００は特徴量を計
算し重み値の算出手段９００は信号の重み値を算出す
る。乗算手段９１０は入力信号の周波数と重み値を乗算
し総和手段９３０は乗算手段９１０の出力を時間と周波
数にわたり総和する。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、音声などの音響信号の
周波数成分を人間の聴覚特性をもちいて重み付けするこ
とにより、雑音が混在した音響信号から目的音を抽出装
置に関する。

【０００２】

【従来の技術】複数の信号が混合した音響信号や雑音が
混在する音響信号からひとつの音響信号だけを取り出す
技術は数種にわたる。例えば、複数のマイクロフォンを
配列しマイクロフォンに到達する信号の時間差や位相
差、レベル差を利用して信号を分離するものもある。特
開平1-118900はその一例である。図１はこの従来の雑音
抑圧装置の構成図を示すものであり、１は音声信号入力
用の接話マイク、２は雑音成分信号入力用のセンサマイ
ク、３及び４はローパスフィルタ、５はバンドパスフィ
ルタバンク、１０ｉから１０Ｎは雑音除去回路、１１ｉ
から１１Ｎは位相差検出補正回路、１２ｉから１２Ｎは
レベル差検出補正回路である。バンドパスフィルタバン
クにより分析された時間信号とそれを整流平滑したスペ
クトル信号を用いて雑音の除去を行うものである。しか
し、この手法ではモノーラルの信号に複数の音響信号が
混在する時には信号の分離ができない。また、雑音の性
質や信号の性質が既知の場合櫛形フィルタや自己相関を
用いた信号分離も試みられている。しかし、混在信号が
既知でない時この手法では信号の分離ができない。

【０００３】ところで、最近音の知覚的な側面に注目し
て音響信号の分離が試みられている。このような手法は
モノーラル信号に複数の信号が混在しているときでも有
効で信号の性質も未知でよく音声が楽器音といった調波
構造をもつ信号の分離に有効な手段である。知覚的な音
の統合のための音響信号の物理的な特徴はＢｒｅｇｍａ
ｎらによって心理学的に調べられている。人間が音を群
化する際周波数成分の高調波関係の有無、周波数成分の
振幅変化の共通性のずれ、音の時間的な連続性などの特
徴が重要なてがかりになっていることがわかっている
（A.S.Bregman Auditory Scene Analysis,Mit Press(19
90)）。

【０００４】Ｂｒｏｗｎらは、この特徴を利用して、信
号分離を試みている（G.J.Brown M.P.Cooke (1992) A
Computational model of auditory scene analysis, P
roceedings of International Conference on Spoken L
anguage Processing(ICSLP)）。システムは、まず人間
の聴覚的特性を模擬したモデルによって周波数解析に相
当する処理を行う。次に周波数成分を抽出し、それらの
立ち上がり時刻のずれ、および高調波という特徴に基づ
いて、ヒューリスティクスを用いて主要な領域をグルー
プ化し、これを目的音声とみなして音響信号を再合成す
る仕組みになっている。しかし、この手法では抽出する
目的信号と除去する信号が周波数上で重なるとき信号が
欠落するため不明瞭な音声信号しか抽出できなかった。
また、その信号の欠落を補償するためＤＰを用いた信号
の補間（M.P.Cooke(1993) Computational auditory sc
ene analysis: Exploiting principles of perceived c
ontinuity Speech Communication 13 391-399）やＨＭ
Ｍをもちいた補間（M.P.Cooke(1994) Learning to rec
ognise speech in noisy environments, Proceedings o
f ATR Workshop on A Biological Framework for Speec
h Perception and Producton）が試みられているが、一
度落ちた情報を再現するためには膨大な計算が必要にも
かかわらず効果は小さい。

【０００５】

【発明が解決しようとする課題】しかしながら、上記の
従来の構成では、 （１）抽出すべき信号の性質が既知でなければならな
い。 （２）複数のマイクロフォンが必要である。 （３）Ｓ／Ｎが改善されても、音声の認識にとって重要
な音響信号も欠落する。 （４）信号の特徴がロバストに抽出されなければならな
いので、計算量が膨大になる。という問題点を有してい
た。

【０００６】本発明は上記従来の課題を解決するもので
パラメータの抽出が不確かであっても抽出すべき信号が
含まれている見込みの大きさに従って入力信号の重みづ
けをすることで完全に抽出すべき信号を欠落されること
がなく音声の明瞭さを損なわず雑音の低減を実現するこ
とを目的としている。

【０００７】

【課題を解決するための手段】この目的を達成するため
に本発明の雑音抑制装置は、音響信号を時間分割する信
号の時間分割手段と、前記信号の時間分割手段による短
時間セグメントの信号のピッチを抽出するピッチ抽出手
段と、前記ピッチ抽出手段によって求められたピッチを
格納する第１の記憶手段と、前記信号の時間分割手段の
出力信号を複数の周波数帯域に分割する周波数帯域分割
手段と、前記周波数帯域分割手段の出力を格納する第２
の記憶手段と、前記第１の記憶手段の出力のエネルギー
を算出するエネルギー算出手段と、前記エネルギー算出
手段の出力を格納する第３の記憶手段と、前記第３の記
憶手段と前記第１の記憶手段の出力から特徴量の算出を
行う特徴量の算出手段と、ある時間区間のある周波数帯
域の信号の重み値の算出を行う重み値の算出手段と、前
記第３の記憶手段の出力と前記重み値の算出手段の積を
算出する乗算手段と、前記乗算手段の出力を格納する第
４の記憶手段と、前記第４の記憶手段の出力の総和を算
出する総和手段を有している。

【０００８】

【作用】本発明は前記構成により、人間の聴覚が聴感上
の雑音抑制のために利用している音の物理的な特徴を用
いて入力信号の周波数の重みづけを行うことにより雑音
が混在した音響信号から雑音を抑制することにより目的
の音響信号を抽出することができる。

【０００９】

【実施例】まず本発明の実施例に共通の構成図を図２に
示す。

【００１０】図２において、１００は入力信号を短時間
セグメントに分割する時間分割手段、２００は入力信号
のピッチを抽出するしピッチ抽出手段、４００は短時間
信号を周波数分割する周波数帯域分割手段、６００はエ
ネルギーまたは実効値の算出手段、８００は特徴量を計
算する特徴量算出手段、９００は重み値算出手段、９１
０は入力信号の周波数と重み値を乗算する乗算手段、９
３０は乗算手段９１０の出力を時間と周波数にわたり総
和する総和手段である。

【００１１】以下それぞれの具体的な実施例について説
明する。 （実施例１）図３は本発明の第１の実施例における雑音
抑制装置の構成を示す構成図である。

【００１２】図３において、１０００は信号の時間切り
出しを行う信号の時間分割手段、１００１はピッチ抽出
手段、１００２は第１の記憶手段、１００３は周波数帯
域分割フィルタ、１００４は第２の記憶手段、１００５
はエネルギーの算出手段、１００６は第３の記憶手段、
１００７は特徴量の算出手段、１００８は重み値の算出
手段、１００９は乗算手段、１０１０は第４の記憶手
段、１０１１は総和手段である。

【００１３】以上のように構成された雑音抑制装置の動
作について説明する。信号の時間分割手段１０００は音
響信号を時間Ｔ毎にδｔの時間セグメントに分割する。
入力信号が音声信号であるので、δｔは１３ｍｓ、Ｔは
２６ｍｓとする。前記信号の時間分割手段１０００の形
状をハニング窓とすると前記信号の時間分割手段Ｗ（ｔ
＋ｋＴ）は（数１）で表すことができる。

【００１４】

【数１】

【００１５】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とするとＳｉ
ｇ（ｔ＋ｋＴ）は（数２）で表される。

【００１６】

【数２】

【００１７】ピッチ抽出手段１００１はＳｉｇ（ｔ＋ｋ
Ｔ）のピッチを抽出する。ピッチは、ケプストラム法、
自己相関法などを用いて抽出される。抽出されたピッチ
をｆ _p（ｋＴ）とする。第１の記憶手段１００２は前記
ピッチ抽出手段１００１で求められたｆ_p（ｋＴ）を格
納する。また、周波数帯域分割フィルタ１００３はＳｉ
ｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯域にＮ分割する。
前記周波数帯域分割フィルタ１００３の出力をＳｉｇ
（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝１，２，・・・・
Ｎ）。ｆ_mは帯域の中心周波数である。Ｎ個の狭帯域信
号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２の記憶手段１００４
に格納される。エネルギーの算出手段１００５では各Ｓ
ｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギーを算出する。前記
エネルギーの算出手段１００５のＮ個の出力をＥｎｅ
（ｋＴ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，ｆ_m）は（数
３）で算出される。

【００１８】

【数３】

【００１９】Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の記憶
手段１００６に格納される。特徴量の算出手段１００７
はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。本実
施例では特徴量として、ピッチ、オンセット、オンセッ
トの同期性、オフセット、オフセットの同期性、連続
性、調波性を算出する。このうちピッチを除く特徴量に
ついて前記特徴量の算出手段１００７で算出する。

【００２０】図４は本実施例における特徴量の算出手段
１００７の詳細図である。１１０１はオンセットの算出
手段、１１０２はオフセットの算出手段、１１０３は連
続性の算出手段、１１０４は調波性の算出手段、１１０
５はオンセットの同期性の算出手段、１１０６はオフセ
ットの同期性の算出手段である。

【００２１】オンセットの算出手段１１０１では前記第
３の記憶手段１００６の出力の時間的な変化に基づきオ
ンセットを定める。つまり、現時刻のエネルギーが前記
時刻のエネルギーに対しある閾値以上に増加したときオ
ンセットであると判断する。

【００２２】オフセットの算出手段１１０２では前記第
３の記憶手段の出力の時間的な変化に基づきオフセット
を定める。つまり、現時刻のエネルギーが前記時刻のエ
ネルギーに対しある閾値以上に減少したときオフセット
であると判断する。

【００２３】連続性の算出手段１１０３では前記第３の
記憶手段１００６の出力の時間的な変化に基づき前時刻
の信号と現時刻の信号が連続であるかどうかを決定す
る。つまり、現時刻のエネルギーが前記時刻のエネルギ
ーに対しある閾値以下で増減したとき連続であると判断
する。

【００２４】調波性の算出手段１１０４ではＳｉｇ（ｔ
＋ｋＴ，ｆ_m）がピッチｆ_p（ｋＴ）に対して調波関係で
あるとき調波的であるとする。つまりピッチの整数倍の
周波数に相当するｆ_mがピッチに対し調波的であると判
断される。オンセットの同期性の算出手段１１０５では
Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセットがピッチｆ
_p（ｋＴ）のオンセット時刻を基準に同時であるとき同
期的であるとする。オフセットの同期性の算出手段１１
０６ではＳｉｇ（ｔ＋ｋＴ，ｆ_m）でのオフセットがピ
ッチｆ_p（ｋＴ）のオフセット時刻を基準に同時である
とき同期的であるとする。

【００２５】図３の重み値の算出手段１００８では、各
特徴量（ピッチ、オンセット、オフセット、オンセット
の同期性、オフセットの同期性、調波性、連続性）を組
み合わせてＳｉｇ（ｔ＋ｋＴ，ｆ_m）での重み値ｗ（ｋ
Ｔ,ｆ_m）を算出する。各特徴量のうち継時的特徴に関す
るものと周波数的特徴に関するものに分類しそれぞれの
統合を行う。前記オンセット、前記オンセットの同期
性、前記オフセット、前記オフセットの同期性、前記連
続性を統合して前記継時的特徴とする。また、前記ピッ
チ、前記調波性を統合して前記周波数的特徴とする。

【００２６】以下にＳｉｇ（ｔ＋ｋＴ，ｆ_m）における
各特徴量の統合方法の一例を示す。前記周波数的特徴及
び前記継時的特徴は論理式（数４）（数５）によって表
現できる。

【００２７】

【数４】

【００２８】

【数５】

【００２９】このとき、Ｓ（ｔ＋ｋＴ，ｆ_m）の重み値
ｗ（ｋＴ，ｆ_m）は（数６）となる。

【００３０】

【数６】

【００３１】論理値を（表１）のように割り当てる。

【００３２】

【表１】

【００３３】あるＳ（ｔ＋ｋＴ，ｆ_m）に対し、（表
２）のような論理値が割り当てられたとき、重み値ｗ
（ｋＴ，ｆ_m）は（数７）となる。

【００３４】

【表２】

【００３５】

【数７】

【００３６】図３の乗算手段１００９はＳｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）と重み値ｗ（ｋＴ，ｆ_m）を乗算する。前記乗
算手段１００９の出力は第４の記憶手段１０１０に格納
される。総和手段１０１１は全てのｋ，ｍについて、ｗ
（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を総和する。
前記総和手段１０１１の出力は雑音が抑制された音響信
号となっている。

【００３７】なお、各特徴量のうち調波性かまたは連続
性が含まれていれば全ての特徴量を用いる必要はない。
また、前記特徴量（ピッチ、オンセット、オフセット、
オンセットの同期性、オフセットの同期性、調波性、連
続性）以外の特徴量を新たに組み合わせてもよい。ま
た、前記信号の時間分割手段１０００の時間幅δｔや分
割の時間周期Ｔは前記δｔ＝１３ｍｓ、Ｔ＝２６ｍｓに
限ったものではなく、形状もハニング窓に限ったもので
はない。また、前記ピッチ抽出手段１００１は、ケプス
トラム法、自己相関法に限ったものではない。また、前
記周波数帯域分割フィルタ１００３の周波数帯域は２０
Ｈｚに限ったものではない。また、前記調波性の算出手
段１１０４において、調波的と判断するものはピッチの
整数倍に限ったものでなく、偶数倍や奇数倍でもよく、
また、それらの周波数を中心に幅をもたせてもよい。ま
た、前記オンセットの同期性の算出手段１１０５におい
て、同期的と判断する時間のずれは０でなくてもよい。
つまり、ピッチ帯域でオンセットが検出された時刻の前
後でも同期的であると判断してもよい。また、前記オフ
セットの同期性の算出手段１１０６において、同期的と
判断する時間のずれは０でなくてもよい。つまり、ピッ
チ帯域でオフセットが検出された時刻の前後でも同期的
であると判断してもよい。また、前記重み値の算出手段
１００８において、各特徴量を論理値（０，１）以外で
表現してもよく、また、特徴量の統合による重み値の算
出方法も論理積や論理和に限ったものではない。また、
オンセットやオフセットや連続性の判断に使うエネルギ
ーの差の算出は隣あう時間セグメントの比較に限らず時
間が離れていたり複数の時間セグメントの比較を行って
もよい。

【００３８】なお、後述する実施例９のように、前記周
波数帯域分割フィルタ１００３はフーリエ変換をもちい
てもよく、また、エネルギーの算出手段１００６は実効
値算出手段をもちいてもよい。

【００３９】（実施例２）以下本発明の信号分離装置に
おける一実施例について図面を参照しながら説明する。

【００４０】本実施例の主たる構成図は、第１の実施例
と同様、図３である。以上のように構成された雑音抑制
装置の動作について説明する。

【００４１】信号の時間分割手段１０００は音響信号を
時間Ｔ毎にδｔの時間セグメントに分割する。入力信号
が音声信号であるので、δｔは１３ｍｓ、Ｔは２６ｍｓ
とする。前記信号の時間分割手段の形状をハニング窓と
すると前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）は（数
１）で表すことができる。

【００４２】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とするとＳｉ
ｇ（ｔ＋ｋＴ）は（数２）で表される。

【００４３】ピッチ抽出手段１００１は複数の抽出方法
を用いてＳｉｇ（ｔ＋ｋＴ）のピッチを抽出する。図５
は本実施例における前記ピッチ抽出手段１００１の詳細
図である。図５において、１２００は第１のピッチ抽出
手段、１２０１は第２のピッチ抽出手段、１２０２は第
ｎのピッチ抽出手段（ｎ＝２，３，・・・・・）、１２
０３はピッチの推定値の算出手段である。第１のピッチ
抽出手段１２００によってピッチｆ_p1（ｋＴ）が、第２
のピッチ抽出手段１２０１によってピッチｆ_p2（ｋＴ）
が、第ｎのピッチ抽出手段１２０２によってピッチｆ_pn
（ｋＴ）が計算されたとする。このときピッチの推定値
の算出手段１２０３によってピッチの推定値ｆ_p（ｋ
Ｔ）は（数８）で求められる。

【００４４】

【数８】

【００４５】図３において、第１の記憶手段１００２
は、前記ピッチ抽出手段１００１で求められたピッチｆ
_p（ｋＴ）を格納する。

【００４６】また、周波数帯域分割フィルタ１００３は
Ｓｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯域にＮ分割す
る。前記周波数帯域分割フィルタ１００３の出力をＳｉ
ｇ（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝１，２，・・・
・Ｎ）。ｆ_mは帯域の中心周波数である。Ｎ個の狭帯域
信号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２の記憶手段１００
４に格納される。エネルギーの算出手段１００５では各
Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギーを算出する。前
記エネルギーの算出手段１００５のＮ個の出力をＥｎｅ
（ｋＴ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，ｆ_m）は（数
３）で算出される。

【００４７】Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の記憶
手段１００６に格納される。特徴量の算出手段１００７
はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。本実
施例では特徴量として、ピッチ、オンセット、オンセッ
トの同期性、オフセット、オフセットの同期性、連続
性、調波性を算出する。このうちピッチを除く特徴量に
ついて前記特徴量の算出手段１００７で算出する。

【００４８】図４は、前記特徴量の算出手段１００７の
詳細図である。１１０１はオンセットの算出手段、１１
０２はオフセットの算出手段、１１０３は連続性の算出
手段、１１０４は調波性の算出手段、１１０５はオンセ
ットの同期性の算出手段、１１０６はオフセットの同期
性の算出手段である。

【００４９】オンセットの算出手段１１０１では前記第
３の記憶手段１００６の出力の時間的な変化に基づきオ
ンセットを定める。つまり、現時刻のエネルギーが前記
時刻のエネルギーに対しある閾値以上に増加したときオ
ンセットであると判断する。

【００５０】オフセットの算出手段１１０２では前記第
３の記憶手段の出力の時間的な変化に基づきオフセット
を定める。つまり、現時刻のエネルギーが前記時刻のエ
ネルギーに対しある閾値以上に減少したときオフセット
であると判断する。

【００５１】連続性の算出手段１１０３では前記第３の
記憶手段１００６の出力の時間的な変化に基づき前時刻
の信号と現時刻の信号が連続であるかどうかを決定す
る。つまり、現時刻のエネルギーが前記時刻のエネルギ
ーに対しある閾値以下で増減したとき連続であると判断
する。

【００５２】調波性の算出手段１１０４ではＳｉｇ（ｔ
＋ｋＴ，ｆ_m）がピッチｆ_p（ｋＴ）に対して調波関係で
あるとき調波的であるとする。つまり、ピッチの整数倍
の周波数に相当するｆ_mがピッチに対し調波的であると
判断される。

【００５３】オンセットの同期性の算出手段１１０５で
はＳｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセットがピッチｆ_p
（ｋＴ）のオンセット時刻を基準に同時であるとき同期
的であるとする。

【００５４】オフセットの同期性の算出手段１１０６で
はＳｉｇ（ｔ＋ｋＴ，ｆ_m）でのオフセットがピッチｆ_p
（ｋＴ）のオフセット時刻を基準に同時であるとき同期
的であるとする。

【００５５】図３の重み値の算出手段１００８では、各
特徴量（ピッチ、オンセット、オフセット、オンセット
の同期性、オフセットの同期性、調波性、連続性）を組
み合わせてＳｉｇ（ｔ＋ｋＴ，ｆ_m）での重み値ｗ（ｔ
＋ｋＴ,ｆ_m）を算出する。各特徴量のうち継時的特徴に
関するものと周波数的特徴に関するものに分類しそれぞ
れの統合を行う。前記オンセット、前記オンセットの同
期性、前記オフセット、前記オフセットの同期性、前記
連続性を統合して前記継時的特徴とする。また、前記ピ
ッチ、前記調波性を統合して前記周波数的特徴とする。

【００５６】以下にＳｉｇ（ｔ＋ｋＴ，ｆ_m）における
前記各特徴量の統合方法の一例を示す。

【００５７】前記周波数的特徴及び前記継時的特徴は論
理式（数４）（数５）によって表現できる。

【００５８】このとき、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）の重
み値ｗ（ｋＴ，ｆ_m）は（数６）となる。

【００５９】論理値を（表１）のように割り当てる。あ
るＳｉｇ（ｔ＋ｋＴ，ｆ_m）に対し、（表２）のような
論理値が割り当てられたとき、重み値ｗ（ｋＴ，ｆ_m）
は（数７）となる。

【００６０】図３の乗算手段１００９はＳｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）と重み値ｗ（ｋＴ，ｆ_m）を乗算する。前記乗
算手段１００９の出力は第４の記憶手段１０１０に格納
される。総和手段１０１１は全てのｋ，ｍについて、ｗ
（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を総和する。
前記総和手段１０１１の出力は雑音が抑制された音響信
号となっている。

【００６１】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また、前記特徴量（ピッチ、オンセット、
オフセット、オンセットの同期性、オフセットの同期
性、調波性、連続性）以外の特徴量を新たに組み合わせ
てもよい。また、前記信号の時間分割手段１０００の時
間幅δｔや分割の時間周期Ｔは前記δｔ＝１３ｍｓ、Ｔ
＝２６ｍｓに限ったものではなく形状もハニング窓に限
ったものではない。また、前記ピッチ抽出手段１００１
は、ケプストラム法、自己相関法に限ったものではな
い。また、前記周波数帯域分割フィルタ１００３の周波
数帯域は２０Ｈｚに限ったものではない。また、前記調
波性の算出手段１１０４において調波的と判断するもの
はピッチの整数倍に限ったものでなく、ピッチの偶数倍
や奇数倍でもよく、また、それらの周波数を中心に幅を
もたせてもよい。また、前記オンセットの同期性の算出
手段１１０５において同期的と判断する時間のずれは０
でなくてもよい。つまり、ピッチ帯域でオンセットが検
出された時刻の前後でも同期的であると判断してもよ
い。また、前記オフセットの同期性の算出手段１１０６
において同期的と判断する時間のずれは０でなくてもよ
い。つまり、ピッチ帯域でオフセットが検出された時刻
の前後でも同期的であると判断してもよい。また、前記
重み値の算出手段１００８において各特徴量を論理値
（０，１）以外で表現してもよく、また、前記特徴量の
統合による重み値の算出方法も論理積や論理和に限った
ものではない。

【００６２】また、オンセットやオフセットや連続性の
判断に使うエネルギーの差の算出は隣あう時間セグメン
トの比較に限らず時間が離れていたり複数の時間セグメ
ントの比較を行ってもよい。

【００６３】なお、後述する実施例９のように、前記周
波数帯域分割フィルタ１００３はフーリエ変換をもちい
てもよく、また前記エネルギーの算出手段１００５は実
効値算出手段をもちいてもよい。

【００６４】（実施例３）以下本発明の信号分離装置に
おける一実施例について図面を参照しながら説明する。

【００６５】本実施例においても主の構成図は図３であ
る。以下、本実施例の動作について説明する。

【００６６】信号の時間分割手段１０００は音響信号を
時間Ｔ毎にδｔの時間セグメントに分割する。入力信号
が音声信号であるのでδｔは１３ｍｓ、Ｔは２６ｍｓと
する。前記信号の時間分割手段１０００の形状はハニン
グ窓とすると信号の時間分割手段Ｗ（ｔ＋ｋＴ）は（数
１）で表すことができる。

【００６７】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とすると、Ｓ
ｉｇ（ｔ＋ｋＴ）は（数２）で表される。

【００６８】ピッチ抽出手段１００１はＳｉｇ（ｔ＋ｋ
Ｔ）のピッチを抽出する。ピッチはケプストラム法、自
己相関法などを用いて抽出される。第１の記憶手段１０
０２は前記ピッチ抽出手段１００１で求められたピッチ
ｆ_p（ｋＴ）を格納する。また、周波数帯域分割フィル
タ１００３はＳｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯
域にＮ分割する。前記周波数帯域分割フィルタ１００３
の出力をＳｉｇ（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝
１，２，・・・・Ｎ）。ｆ_mは帯域の中心周波数であ
る。Ｎ個の狭帯域信号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２
の記憶手段１００４に格納される。エネルギーの算出手
段１００５では各Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギ
ーを算出する。前記エネルギーの算出手段１００５のＮ
個の出力をＥｎｅ（ｋＴ，ｆ_m）と記述する。Ｅｎｅ
（ｋＴ，ｆ_m）は（数３）で算出される。

【００６９】Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の記憶
手段１００６に格納される。特徴量の算出手段１００７
はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。本実
施例では特徴量として、ピッチ、オンセット、オンセッ
トの同期性、オフセット、オフセットの同期性、連続
性、調波性を算出する。

【００７０】図６は本実施例における特徴量の算出手段
１００７の詳細図である。図６において、１３０１はピ
ッチの評価測度の算出手段、１３０２はオンセットの評
価測度の算出手段、１３０３はオフセットの評価測度の
算出手段、１３０４は連続性の評価測度の算出手段、１
３０５は調波性の評価測度の算出手段、１３０６はオン
セットの同期性の評価測度の算出手段、１３０７はオフ
セットの同期性の評価測度の算出手段である。また、評
価測度としてＳｉｇ（ｔ＋ｋＴ，ｆ_m）で各特徴量が存
在する確からしさを用いる。

【００７１】以下で各特徴量の評価測度の算出手段につ
いて詳細に説明する。ピッチの評価測度の算出手段１３
０１において、ｆ_p（ｋＴ）に対しピッチが存在しうる
範囲を確率的に定義しＳｉｇ（ｋＴ，ｆ_m）にピッチが
存在する確からしさを算出する。図７はピッチが存在し
うる周波数の確率分布の一例である。ｆ_p（ｋＴ）に対
してピッチが存在する確からしさがガウス分布をすると
するとＳｉｇ（ｋＴ，ｆ_m）にピッチが存在する確から
しさＲｅｌ_{_pitch}（ｋＴ，ｆ_m）は分布の分散をσとす
ると（数９）によって与えられる。

【００７２】

【数９】

【００７３】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【００７４】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。オンセットの評価測度の算出手段１３０２ではこ
の変化が一つの音の時間変化ではなく図９（ｂ）のよう
に異なるもう一つの信号が加わったために起こったもの
であるという確からしさを人間の知覚に基づき定めこれ
をオンセットの確からしさとする。オンセットの確から
しさは前記第３の記憶手段１００６の出力によって算出
される。以下に算出方法の一つの例を挙げる。心理学的
に連続にこえるか新たに音が加わったかという判断基準
が丁度変化する音の大きさが調べられている（1994年日
本音響学会秋期研究発表会稿論集2-7-18）。前記エネル
ギーの算出手段１００５の出力を用いて現観測時刻のエ
ネルギーとそれより前の時刻までのエネルギーの平均か
ら前記オンセットの確からしさを定める。図１０はオン
セットの確からしさの算出方法のための説明図である。
前記前時刻までのエネルギーの平均は現観測周波数とそ
の周波数に近接する２つの周波数帯域における約２００
ｍｓ区間（１６個のセグメント）のセグメントの平均と
する。Ｓｉｇ（ｋＴ，ｆ_m）の前記オンセットの確から
しさＲｅｌ_{_on}（Ｅ_Dif）は（数１０）で算出される。

【００７５】

【数１０】

【００７６】図１１は（数１０）を図式化したものであ
る。前記オンセットの確からしさはピッチ帯域でオンセ
ットが検出されるまではピッチ帯域でのみ計算される。
図１０及び図１１は図８の信号に対する前記オンセット
の確からしさを説明する説明図である。図１２において
時刻ｋ＝２でピッチ帯域の前記オンセットの確からしさ
が１．０となり閾値０．８を越えたので時刻ｋ＝２から
全周波数帯域で前記オンセットの確からしさの計算を始
める。次に、（ｋ，ｍ）＝（２，１）と（２，３）で前
記オンセットの確からしさ１．０が算出できたとする。
この時図１３のように（ｋ，ｍ）＝（２，１）と（２，
３）に前記オンセットの確からしさ１．０を割り当て
る。

【００７７】オフセットの評価測度の算出手段１３０３
においてオフセットの確からしさが算出される。二つの
信号から一方の信号が除かれたとき音の大きさには変化
がある。そこで前記エネルギーの算出手段１００５の出
力を用いてオフセットを定める。前記オフセットの確か
らしさＲｅｌ_{_off}（Ｅ_Dif）は前記オンセットの確から
しさの算出方法と同様に（数１１）で与えられる。

【００７８】

【数１１】

【００７９】前記オフセットの確からしさはオンセット
が検出されてからピッチ帯域でオフセットが検出される
までピッチ周波数帯域のみ計算される。図１４及び図１
５は図８の信号に対する前記オフセットの確からしさを
説明する説明図である。図１４において時刻ｋ＝６でピ
ッチ帯域の前記オフセットの確からしさが１．０となり
ある閾値０．９を越えたので時刻ｋ＝６において前記オ
フセットの確からしさを全周波数帯域について算出す
る。（ｋ，ｍ）＝（６，１）と（６，３）で前記オフセ
ットの確からしさ１．０を算出できたとすると図１５の
ように前記オフセットの確からしさを割り当てる。

【００８０】連続性の評価測度の算出手段１３０４では
音の大きさが時間的に変化したときそれが一つの信号の
時間変化である確からしさを指標とする。連続性の確か
らしさは前記第３の記憶手段１００６の出力によって算
出される。図１６は前記連続性の確からしさの算出方法
のための説明図である。図１６において、現観測時刻を
含む現観測時刻以前約２００ｍｓ区間（１６個のセグメ
ント）の現観測セグメントのエネルギーの平均と現観測
周波数とその周波数に隣接する２つの周波数帯域におけ
る現観測時刻を含まない現観測時刻より前の約２００ｍ
ｓ区間（１６個のセグメント）のエネルギーの平均から
前記連続性の確からしさを決定する。前記連続性の確か
らしさＲｅｌ_{_con}（Ｅ_Dif）は（数１２）で定義する。

【００８１】

【数１２】

【００８２】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさはオンセットが検出された
後オフセットが検出されるまで計算される。図１８は図
８の信号に対する前記連続性の確からしさを説明する説
明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６でオ
フセットが検出されたとき、時刻ｋ＝３から５までの間
前記連続性の確からしさが計算され各（ｋ，ｍ）での連
続性の確からしさが図１８のように算出される。

【００８３】調波性の評価測度の算出手段１３０５では
Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）がｆ_p（ｋＴ）に対して知覚上
調波関係である確からしさをもって調波性の確からしさ
とする。ピッチの整数倍の周波数区間を中心に主観的に
調波性を感じる周波数分布を人間の心理測定関数から定
める。心理測定関数は電子情報通信学会論文誌'94/5Vo
l.J77-ANo.5 731-740に掲載のデータを参考に決める。
周波数帯域の中心周波数ｆ_mでの前記ピッチｆ_p（ｋＴ）
のｐ倍（ｐは整数）の周波数からのずれｕが（数１３）
で求められる。

【００８４】

【数１３】

【００８５】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【００８６】

【数１４】

【００８７】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【００８８】オンセットの同期性の評価測度の算出手段
１３０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5 731-740に掲載のデ
ータを参考に定義する。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でオ
ンセットが検出されたときそれと最も近いピッチのオン
セット時刻の差をｖ（ｍｓ）とすると、前記オンセット
の同期性の確からしさＲｅｌ_sync（ｋＴ，ｆ_m）は（数
１５）で記述できる。

【００８９】

【数１５】

【００９０】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【００９１】オフセットの同期性の評価測度の算出手段
１３０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【００９２】図３の重み値の算出手段１００８では評価
関数に前記各特徴（ピッチ、オンセット、オフセット、
オンセットの同期性、オフセットの同期性、調波性、連
続性）についての確からしさを入力し、Ｓｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算出する。

【００９３】各特徴量のうち継時的確からしさに関する
ものと周波数的確からしさに関するものに分類しそれぞ
れの確からしさの統合を行う。前記オンセットの確から
しさ、前記オンセットの同期性の確からしさ、前記オフ
セットの確からしさ、前記オフセットの同期性の確から
しさ、前記連続性の確からしさを統合して前記継時的確
からしさとする。また、前記ピッチの確からしさ、前記
調波性の確からしさを統合して前記周波数的確からしさ
とする。

【００９４】図３において乗算手段１００９はＳｉｇ
（ｔ＋ｋＴ，ｆ_m）と重み値ｗ（ｋＴ，ｆ_m）を乗算す
る。前記乗算手段１００９の出力は第４の記憶手段１０
１０に格納される。総和手段１０１１は全てのｋ，ｍに
ついて、ｗ（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を
総和する。前記総和手段１０１１の出力は雑音が抑制さ
れた音響信号となっている。

【００９５】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また、前記特徴量以外の特徴量を新たに組
み合わせてもよい。また、前記信号の時間分割手段１０
００の時間幅δｔや分割の時間周期Ｔは前記δｔ＝１３
ｍｓ、Ｔ＝２６ｍｓに限ったものではなく形状もハニン
グ窓に限ったものではない。また、前記ピッチ抽出手段
１００１はケプストラム法、自己相関法に限ったもので
はない。また、前記周波数帯域分割フィルタ１００３の
周波数帯域は２０Ｈｚに限ったものではない。また、ピ
ッチが存在する確からしさの算出においてガウス分布以
外の分布で定義してもよい。また、オンセット時刻の判
断に用いる閾値は０．８以外でもよい。また、オフセッ
ト時刻の判断に用いる閾値は０．９以外でもよい。ま
た、前記オンセットの確からしさの算出方法は（数１
０）に限らない。また、前記オフセットの確からしさの
算出方法は（数１１）に限らない。また、前記連続性の
確からしさの算出方法は（数１２）に限らない。また、
前記調波性の確からしさの算出方法は（数１３）および
（数１４）に限らない。また、オンセットの同期性の確
からしさの算出方法やオフセットの同期性の確からしさ
の算出方法は（数１５）に限らない。

【００９６】なお、後述する実施例９のように、前記周
波数帯域分割フィルタ１００３はフーリエ変換をもちい
てもよく、またエネルギーの算出手段１００６は実効値
算出手段をもちいてもよい。

【００９７】（実施例４）以下本発明の信号分離装置に
おける一実施例について、図面を参照しながら説明す
る。

【００９８】本実施例も主なる構成図は図３である。本
実施例の動作に付いて説明する。

【００９９】信号の時間分割手段１０００は音響信号を
時間Ｔ毎にδｔの時間セグメントに分割する。入力信号
が音声信号であるのでδｔは１３ｍｓ、Ｔは２６ｍｓと
する。前記信号の時間分割手段１０００の形状をハニン
グ窓とすると前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）は
（数１）で表すことができる。

【０１００】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とすると、Ｓ
ｉｇ（ｔ＋ｋＴ）は（数２）で表される。

【０１０１】ピッチ抽出手段１００１は複数の抽出方法
を用いてＳｉｇ（ｔ＋ｋＴ）のピッチを抽出する。図２
４は前記ピッチ抽出手段１００１の詳細図である。図２
４において、３８１０は第１のピッチ抽出手段、３８２
０は第２のピッチ抽出手段、３８３０はピッチの推定値
の算出手段である。第１のピッチ抽出手段３８１０によ
ってピッチｆ_p1（ｋＴ）が、第２のピッチ抽出手段３８
２０によってピッチｆ _p2（ｋＴ）が計算されたとする。
このとき、ピッチの推定値の算出手段３８３０によって
ピッチの推定値ｆ_p（ｋＴ）は（数１６）で求められ
る。

【０１０２】

【数１６】

【０１０３】図３において、第１の記憶手段１００２は
前記ピッチ抽出手段１００１で求められたｆ_p1（ｋＴ）
とｆ_p2（ｋＴ）とｆ_p（ｋＴ）を格納する。

【０１０４】また、周波数帯域分割フィルタ１００３は
Ｓｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯域にＮ分割す
る。前記周波数帯域分割フィルタ１００３の出力をＳｉ
ｇ（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝１，２，・・・
・Ｎ）。ｆ_mは帯域の中心周波数である。Ｎ個の狭帯域
信号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２の記憶手段１００
４に格納される。エネルギーの算出手段１００５では各
Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギーを算出する。前
記エネルギーの算出手段１００５のＮ個の出力をＥｎｅ
（ｋＴ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，ｆ_m）は（数
３）で算出される。

【０１０５】前記Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の
記憶手段１００６に格納される。特徴量の算出手段１０
０７はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。
本実施例では特徴量として、ピッチ、オンセット、オン
セットの同期性、オフセット、オフセットの同期性、連
続性、調波性を算出する。

【０１０６】図６は本実施例における特徴量の算出手段
１００７の詳細図である。図６において、１３０１はピ
ッチの評価測度の算出手段、１３０２はオンセットの評
価測度の算出手段、１３０３はオフセットの評価測度の
算出手段、１３０４は連続性の評価測度の算出手段、１
３０５は調波性の評価測度の算出手段、１３０６はオン
セットの同期性の評価測度の算出手段、１３０７はオフ
セットの同期性の評価測度の算出手段である。また、評
価測度としてＳｉｇ（ｔ＋ｋＴ，ｆ_m）で各特徴量が存
在する確からしさを用いる。

【０１０７】以下で前記特徴量の評価測度の算出手段に
ついて詳細に説明する。ピッチの評価測度の算出手段１
３０１は、前記ピッチ抽出手段１００１で計算されたｆ
_p1（ｋＴ）とｆ_p2（ｋＴ）の差からピッチｆ_p（ｋＴ）
の確からしさを決定する。ｆ_p1（ｋＴ）とｆ_p2（ｋＴ）
の差が０の時ピッチは正確に抽出されたとしピッチｆ_p
（ｋＴ）の確からしさは１とする。ｆ_p1（ｋＴ）とｆ_p2
（ｋＴ）の差が大きくなると確からしさは減少する。ピ
ッチの確からしさが図３９００のようにｆ_p1（ｋＴ）−
ｆ_p2（ｋＴ）の値によってガウス分布するとするとピッ
チの確からしさＲｅｌ（ｆ_p1（ｋＴ）−ｆ_p2（ｋＴ））
が（数１７）で決定される。

【０１０８】

【数１７】

【０１０９】同時に、図２６のように得られたピッチの
推定値ｆ_p（ｋＴ）に対しあるピッチが存在しうる範囲
を確率的に定義しＳｉｇ（ｋＴ，ｆ_m）にピッチが存在
する確からしさを算出する。ピッチの推定値ｆ_p（ｋ
Ｔ）に対してピッチが存在する確からしさがガウス分布
をするとすると、Ｓｉｇ（ｋＴ，ｆ_m）にピッチが存在
する確からしさＲｅｌ_{_pitch}（ｋＴ，ｆ_m）は分散をσ
とすると（数９）によって与えられる。

【０１１０】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【０１１１】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。オンセットの評価測度の算出手段１３０２ではこ
の変化が一つの音の時間変化ではなく図９（ｂ）のよう
に異なるもう一つの信号が加わったために起こったもの
であるという確からしさを人間の知覚に基づき定めこれ
をオンセットの確からしさとする。前記オンセットの確
からしさは前記第３の記憶手段１００６の出力によって
算出される。以下に算出方法の一つの例を挙げる。心理
学的に連続にきこえるか新たに音が加わったかという判
断基準が丁度変化する音の大きさが調べられている（19
94年日本音響学会秋期研究発表会稿論集2-7-18）。前記
エネルギーの算出手段１００５の出力を用いて現観測時
刻のエネルギーとそれより前の時刻までのエネルギーの
平均から前記オンセットの確からしさを定める。図１０
はオンセットの確からしさの算出方法のための説明図で
ある。前記前時刻までのエネルギーの平均は現観測周波
数とその周波数に近接する２つの周波数帯域における約
２００ｍｓ区間（１６個のセグメント）のセグメントの
平均とする。Ｓｉｇ（ｋＴ，ｆ_m）の前記オンセットの
確からしさＲｅｌ_{_on}（Ｅ_Dif）は（数１０）で算出され
る。

【０１１２】図１１は（数１０）を図式化したものであ
る。オンセットの確からしさはピッチ帯域でオンセット
が検出されるまではピッチ帯域でのみ計算される。図１
２及び図１３は図８の信号に対する前記オンセットの確
からしさを説明する説明図である。図１２において時刻
ｋ＝２でピッチ帯域の前記オンセットの確からしさが
１．０となり閾値０．８を越えたので時刻ｋ＝２から全
周波数帯域で前記オンセットの確からしさの計算を始め
る。（ｋ，ｍ）＝（２，１）と（２，３）で前記オンセ
ットの確からしさ１．０が算出できたとする。この時図
２７００のように（ｋ，ｍ）＝（２，１）と（２，３）
に前記オンセットの確からしさ１．０を割り当てる。

【０１１３】オフセットの評価測度の算出手段１３０３
においてオフセットの確からしさが算出される。二つの
信号から一方の信号が除かれたとき音の大きさには変化
がある。そこで前記エネルギーの算出手段１００５の出
力を用いてオフセットを定める。前記オフセットの確か
らしさＲｅｌ_{_off}（Ｅ_Dif）は前記オンセットの確から
しさの算出方法と同様に（数１１）で与えられる。

【０１１４】前記オフセットの確からしさはオンセット
が検出されてからピッチ帯域でオフセットが検出される
までピッチ周波数帯域のみ計算される。図１４及び図１
５は図８の信号に対する前記オフセットの確からしさを
説明する説明図である。図１４において時刻ｋ＝６でピ
ッチ帯域の前記オフセットの確からしさが１．０となり
ある閾値０．９を越えたので時刻ｋ＝６において前記オ
フセットの確からしさを全周波数帯域について算出す
る。（ｋ，ｍ）＝（６，１）と（６，３）で前記オフセ
ットの確からしさ１．０を算出できたとすると図１５の
ように前記オフセットの確からしさを割り当てる。

【０１１５】連続性の評価測度の算出手段１３０４で
は、音の大きさが時間的に変化したときそれが一つの信
号の時間変化である確からしさを指標とする。前記連続
性の確からしさは前記第３の記憶手段１００６の出力に
よって算出される。図１６はオンセットの確からしさの
算出方法のための説明図である。図１６において、現観
測時刻を含む現観測時刻以前約２００ｍｓ区間（１６個
のセグメント）の現観測セグメントのエネルギーの平均
と現観測周波数とその周波数に隣接する２つの周波数帯
域における現観測時刻を含まない現観測時刻より前の約
２００ｍｓ区間（１６個のセグメント）のエネルギーの
平均から前記連続性の確からしさを決定する。前記連続
性の確からしさＲｅｌ_{_con}（Ｅ_Dif）は（数１２）で定
義する。

【０１１６】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさは、オンセットが検出され
た後オフセットが検出されるまで計算される。図１８は
図８の信号に対する前記連続性の確からしさを説明する
説明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６で
オフセットが検出されたとき、時刻ｋ＝３から５までの
間前記連続性の確からしさが計算され、各（ｋ，ｍ）で
の連続性の確からしさが図１８のように算出される。

【０１１７】調波性の評価測度の算出手段１３０５で
は、前記Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）が前記ｆ_p（ｋＴ）に
対して知覚上調波関係である確からしさをもって調波性
の確からしさとする。ピッチの整数倍の周波数区間を中
心に主観的に調波性を感じる周波数分布を人間の心理測
定関数から定める。心理測定関数は電子情報通信学会論
文誌'94/5Vol.J77-ANo.5 731-740に掲載のデータを参
考に決める。周波数帯域の中心周波数ｆ_mでの前記ピッ
チｆ_p（ｋＴ）のｐ倍（ｐは整数）の周波数からのずれ
ｕが（数１３）で求められる。

【０１１８】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【０１１９】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【０１２０】オンセットの同期性の評価測度の算出手段
１３０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5 731-740に掲載のデ
ータを参考に定義する。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でオ
ンセットが検出されたときそれと最も近いピッチのオン
セット時刻の差をｖ（ｍｓ）とすると、前記オンセット
の同期性の確からしさＲｅｌ_sync（ｋＴ，ｆ_m）は、
（数１５）で記述できる。

【０１２１】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【０１２２】オフセットの同期性の評価測度の算出手段
１３０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【０１２３】図３の重み値の算出手段１００８では、評
価関数に各特徴（ピッチ、オンセット、オフセット、オ
ンセットの同期性、オフセットの同期性、調波性、連続
性）についての確からしさを入力しＳｉｇ（ｔ＋ｋＴ，
ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算出する。

【０１２４】各特徴量のうち継時的確からしさに関する
ものと周波数的確からしさに関するものに分類しそれぞ
れの確からしさの統合を行う。前記オンセットの確から
しさ、前記オンセットの同期性の確からしさ、前記オフ
セットの確からしさ、前記オフセットの同期性の確から
しさ、前記連続性の確からしさを統合して前記継時的確
からしさとする。また、前記ピッチの確からしさ、前記
ピッチの変動の範囲の確からしさ、前記調波性の確から
しさを統合して前記周波数的確からしさとする。

【０１２５】図３において、乗算手段１００９はＳｉｇ
（ｔ＋ｋＴ，ｆ_m）と重み値ｗ（ｋＴ，ｆ_m）を乗算す
る。前記乗算手段１００９の出力は第４の記憶手段１０
１０に格納される。総和手段１０１１は全てのｋ，ｍに
ついて、ｗ（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を
総和する。前記総和手段１０１１の出力は雑音が抑制さ
れた音響信号となっている。

【０１２６】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また前記特徴量以外の特徴量を新たに組み
合わせてもよい。また、前記信号の時間分割手段１００
０の時間幅δｔや分割の時間周期Ｔは前記δｔ＝１３ｍ
ｓ、Ｔ＝２６ｍｓに限ったものではなく形状もハニング
窓に限ったものではない。また、前記ピッチ抽出手段１
００１はケプストラム法、自己相関法に限ったものでは
ない。また、前記周波数帯域分割フィルタ１００３の周
波数帯域は２０Ｈｚに限ったものではない。また、ピッ
チが存在する確からしさの算出においてガウス分布以外
の分布で定義してもよい。また、オンセット時刻の判断
に用いる閾値は０．８以外でもよい。また、オフセット
時刻の判断に用いる閾値は０．９以外でもよい。また、
前記オンセットの確からしさの算出方法は（数１０）に
限らない。また、前記オフセットの確からしさの算出方
法は（数１１）に限らない。また、前記連続性の確から
しさの算出方法は（数１２）に限らない。また、前記調
波性の確からしさの算出方法は（数１３）および（数１
４）に限らない。また、オンセットの同期性の確からし
さの算出方法やオフセットの同期性の確からしさの算出
方法は（数１５）に限らない。

【０１２７】なお、実施例９のように、前記周波数帯域
分割フィルタ１００３はフーリエ変換を用いてもよく、
またエネルギーの算出手段１００６は実効値算出手段を
用いてもよい。

【０１２８】（実施例５）以下本発明の信号分離装置に
おける一実施例について、図面を参照しながら説明す
る。

【０１２９】本発明の第５の実施例における主な構成図
は図３である。以上のように構成された雑音抑制装置の
動作について説明する。

【０１３０】信号の時間分割手段１０００は音響信号を
時間Ｔ毎にδｔの時間セグメントに分割する。入力信号
が音声信号であるのでδｔは１３ｍｓ、Ｔは２６ｍｓと
する。前記信号の時間分割手段１０００の形状はハニン
グ窓とすると前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）は
（数１）で表すことができる。

【０１３１】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とするとＳｉ
ｇ（ｔ＋ｋＴ）は（数２）で表される。

【０１３２】ピッチ抽出手段１００１は複数の抽出方法
を用いてＳｉｇ（ｔ＋ｋＴ）のピッチを抽出する。図５
は本実施例における前記ピッチ抽出手段１００１の詳細
図である。図５において、１２００は第１のピッチ抽出
手段、１２０１は第２のピッチ抽出手段、１２０２は第
ｎのピッチ抽出手段、１２０３はピッチの推定値ｆ
_p（ｔ）の算出手段である。第１のピッチ抽出手段１２
００によってピッチｆ_p1（ｋＴ）が、第２のピッチ抽出
手段１２０１によってピッチｆ_p2（ｋＴ）が、第ｎのピ
ッチ抽出手段１２０２によってピッチｆ_pn（ｋＴ）が計
算されたとする。このとき、ピッチの推定値ｆ_p（ｔ）
の算出手段１２０３によってピッチの推定値ｆ_p（ｋ
Ｔ）は（数８）で求められる。

【０１３３】また、ｆ_p1（ｋＴ）からｆ_pn（ｋＴ）の分
散Ｖ_pitch（ｋＴ）を（数１８）で求める。

【０１３４】

【数１８】

【０１３５】図３において、第１の記憶手段１００２は
前記ピッチ抽出手段１００１で求められたｆ_p（ｋＴ）
と分散Ｖ_pitch（ｋＴ）を格納する。

【０１３６】また、周波数帯域分割フィルタ１００３は
Ｓｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯域にＮ分割す
る。前記周波数帯域分割フィルタ１００３の出力をＳｉ
ｇ（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝１，２，・・・
・Ｎ）。ｆ_mは帯域の中心周波数である。Ｎ個の狭帯域
信号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２の記憶手段１００
４に格納される。エネルギー算出手段１００５では各Ｓ
ｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギーを算出する。前記
エネルギー算出手段１００５のＮ個の出力をＥｎｅ（ｋ
Ｔ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，ｆ_m）は（数３）
で算出される。

【０１３７】前記Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の
記憶手段１００６に格納される。特徴量の算出手段１０
０７はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。
本実施例では、特徴量として、ピッチ、オンセット、オ
ンセットの同期性、オフセット、オフセットの同期性、
連続性、調波性を算出する。

【０１３８】図６は図３の特徴量の算出手段１００７の
詳細図である。図６において、１３０１はピッチの評価
測度の算出手段、１３０２はオンセットの評価測度の算
出手段、１３０３はオフセットの評価測度の算出手段、
１３０４は連続性の評価測度の算出手段、１３０５は調
波性の評価測度の算出手段、１３０６はオンセットの同
期性の評価測度の算出手段、１３０７はオフセットの同
期性の評価測度の算出手段である。また、評価測度とし
てＳｉｇ（ｔ＋ｋＴ，ｆ_m）で各特徴量が存在する確か
らしさを用いる。

【０１３９】以下で前記特徴量の評価測度の算出手段に
ついて詳細に説明する。ピッチの評価測度の算出手段
は、前記ピッチ抽出手段１００１で計算された分散Ｖ
_pitch（ｋＴ）から推定されたピッチｆ_p（ｋＴ）の確か
らしさを決定する。Ｖ_pitch（ｋＴ）が０の時ピッチは
正確に抽出されたとしピッチｆ_p（ｋＴ）の確からしさ
は１とする。Ｖ_pitch（ｋＴ）が大きくなると確からし
さは減少する。ピッチの確からしさが図４１００のよう
にＶ_pitch（ｋＴ）の値によってガウス分布するとする
とピッチの確からしさＲｅｌ（Ｖ_pitch（ｋＴ））が
（数１９）で決定される。

【０１４０】

【数１９】

【０１４１】同時に、図２８のように得られたピッチの
推定値ｆ_p（ｋＴ）に対しあるピッチが存在しうる範囲
を確率的に定義しＳｉｇ（ｋＴ，ｆ_m）がピッチである
確からしさを算出する。ピッチの推定値ｆ_p（ｋＴ）に
対してピッチが存在する確からしさがガウス分布をする
とするとＳｉｇ（ｋＴ，ｆ_m）にピッチが存在する確か
らしさＲｅｌ_{_pitch}（ｋＴ，ｆ_m）は（数２０）で与え
られる。

【０１４２】

【数２０】

【０１４３】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【０１４４】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。オンセットの評価測度の算出手段１３０２ではこ
の変化が一つの音の時間変化ではなく図９（ｂ）のよう
に異なるもう一つの信号が加わったために起こったもの
であるという確からしさを人間の知覚に基づき定めこれ
をオンセットの確からしさとする。前記オンセットの確
からしさは前記第３の記憶手段１００６の出力によって
算出される。以下に算出方法の一つの例を挙げる。心理
学的に連続にきこえるか新たに音が加わったかという判
断基準が丁度変化する音の大きさが調べられている（19
94年日本音響学会秋期研究発表会稿論集2-7-18）。前記
エネルギー算出手段１００５の出力を用いて現観測時刻
のエネルギーとそれより前の時刻までのエネルギーの平
均から前記オンセットの確からしさを定める。図１０は
オンセットの確からしさの算出方法のための説明図であ
る。前記前時刻までのエネルギーの平均は現観測周波数
とその周波数に近接する２つの周波数帯域における約２
００ｍｓ区間（１６個のセグメント）のセグメントの平
均とする。Ｓｉｇ（ｋＴ，ｆ_m）の前記オンセットの確
からしさＲｅｌ_{_on}（Ｅ_Dif）は（数１０）で算出され
る。

【０１４５】図１１は（数１０）を図式化したものであ
る。前記オンセットの確からしさはピッチ帯域でオンセ
ットが検出されるまでピッチ帯域でのみ計算される。図
１２及び図１３は図８の信号に対する前記オンセットの
確からしさを説明する説明図である。図１２において時
刻ｋ＝２でピッチ帯域の前記オンセットの確からしさが
１．０となり閾値０．８を越えたので時刻ｋ＝２から全
周波数帯域で前記オンセットの確からしさの計算を始め
る。（ｋ，ｍ）＝（２，１）と（２，３）で前記オンセ
ットの確からしさ１．０が算出できたとする。この時、
図１３のように（ｋ，ｍ）＝（２，１）と（２，３）に
前記オンセットの確からしさ１．０を割り当てる。

【０１４６】前記オフセットの評価測度の算出手段１３
０３においてオフセットの確からしさが算出される。二
つの信号から一方の信号が除かれたとき音の大きさには
変化がある。そこで前記エネルギーの算出手段１００５
の出力を用いてオフセットを定める。前記オフセットの
確からしさＲｅｌ_{_off}（Ｅ_Dif）は前記オンセットの確
からしさの算出方法と同様に（数１１）で与えられる。

【０１４７】前記オフセットの確からしさはオンセット
が検出されてからピッチ帯域でオフセットが検出される
までピッチ周波数帯域のみ計算される。図１４及び図１
５は図８の信号に対する前記オフセットの確からしさを
説明する説明図である。図１４において時刻ｋ＝６でピ
ッチ帯域の前記オフセットの確からしさが１．０となり
ある閾値０．９を越えたので時刻ｋ＝６において前記オ
フセットの確からしさを全周波数帯域について算出す
る。（ｋ，ｍ）＝（６，１）と（６，３）で前記オフセ
ットの確からしさ１．０を算出できたとすると図１５の
ように前記オフセットの確からしさを割り当てる。

【０１４８】連続性の評価測度の算出手段１３０４で
は、音の大きさが時間的に変化したときそれが一つの信
号の時間変化である確からしさを指標とする。前記連続
性の確からしさは前記第３の記憶手段１００６の出力に
よって算出される。図１６は前記連続性の確からしさの
算出方法のための説明図である。図１６において、現観
測時刻を含む現観測時刻以前約２００ｍｓ区間（１６個
のセグメント）の現観測セグメントのエネルギーの平均
と現観測周波数とその周波数に隣接する２つの周波数帯
域における現観測時刻を含まない現観測時刻より前の約
２００ｍｓ区間（１６個のセグメント）のエネルギーの
平均から前記連続性の確からしさを決定する。前記連続
性の確からしさＲｅｌ_{_con}（Ｅ_Dif）は（数１２）で定
義する。

【０１４９】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさは、オンセットが検出され
た後オフセットが検出されるまで計算される。図１８は
図８の信号に対する前記連続性の確からしさを説明する
説明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６で
オフセットが検出されたとき時刻ｋ＝３から５までの間
前記連続性の確からしさが計算され各（ｋ，ｍ）での連
続性の確からしさが図１８のように算出される。

【０１５０】調波性の評価測度の算出手段１３０５で
は、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）がｆ_p（ｋＴ）に対して知
覚上調波関係である確からしさをもって調波性の確から
しさとする。ピッチの整数倍の周波数区間を中心に主観
的に調波性を感じる周波数分布を人間の心理測定関数か
ら定める。心理測定関数は電子情報通信学会論文誌'94/
5Vol.J77-ANo.5 731-740に掲載のデータを参考に決め
る。周波数帯域の中心周波数ｆ_mでの前記ピッチｆ_p（ｋ
Ｔ）のｐ倍（ｐは整数）の周波数からのずれｕが（数１
３）で求められる。

【０１５１】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【０１５２】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【０１５３】オンセットの同期性の評価測度の算出手段
１３０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5 731-740に掲載のデ
ータを参考に前記同期性の確からしさの分布Ｒｅｌ_syn
c（ｋＴ，ｆ_m）を（数１５）で定義する。Ｓｉｇ（ｔ＋
ｋＴ，ｆ_m）でオンセットが検出されたときそれと最も
近いピッチのオンセット時刻の差をｖ（ｍｓ）とする
と、前記オンセットの同期性の確からしさＲｅｌ_sync
（ｋＴ，ｆ_m）は、（数１５）で記述できる。

【０１５４】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【０１５５】オフセットの同期性の評価測度の算出手段
１３０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【０１５６】図３の重み値の算出手段１００８では、評
価関数に各特徴（ピッチ、オンセット、オフセット、オ
ンセットの同期性、オフセットの同期性、調波性、連続
性）についての確からしさを入力し、Ｓｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算出する。

【０１５７】前記各特徴量のうち継時的確からしさに関
するものと周波数的確からしさに関するものに分類しそ
れぞれの確からしさの統合を行う。前記オンセットの確
からしさ、前記オンセットの同期性の確からしさ、前記
オフセットの確からしさ、前記オフセットの同期性の確
からしさ、前記パワーの連続性の確からしさを統合して
前記継時的確からしさとする。また、前記ピッチの確か
らしさ、前記ピッチの変動の範囲の確からしさ、前記調
波性の確からしさを統合して前記周波数的確からしさと
する。

【０１５８】図３において、前記乗算手段１００９はＳ
ｉｇ（ｔ＋ｋＴ，ｆ_m）と重み値ｗ（ｋＴ，ｆ_m）を乗算
する。前記乗算手段１００９の出力は第４の記憶手段１
０１０に格納される。総和手段１０１１は全てのｋ，ｍ
について、ｗ（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）
を総和する。前記総和手段１０１１の出力は雑音が抑制
された音響信号となっている。

【０１５９】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また前記特徴量以外の特徴量を新たに組み
合わせてもよい。また、前記信号の時間分割手段１００
０の時間幅δｔや分割の時間周期Ｔは前記δｔ＝１３ｍ
ｓ、Ｔ＝２６ｍｓに限ったものではなく、形状もハニン
グ窓に限ったものではない。また、前記ピッチ抽出手段
１００１はケプストラム法、自己相関法に限ったもので
はない。また、前記周波数帯域分割フィルタ１００３の
周波数帯域は２０Ｈｚに限ったものではない。また、ピ
ッチが存在する確からしさの算出においてガウス分布以
外の分布で定義してもよい。また、オンセット時刻の判
断に用いる閾値は０．８以外でもよい。また、オフセッ
ト時刻の判断に用いる閾値は０．９以外でもよい。ま
た、前記オンセットの確からしさの算出方法は（数１
０）に限らない。また、前記オフセットの確からしさの
算出方法は（数１１）に限らない。また、前記連続性の
確からしさの算出方法は（数１２）に限らない。また、
前記調波性の確からしさの算出方法は（数１３）および
（数１４）に限らない。また、オンセットの同期性の確
からしさの算出方法やオフセットの同期性の確からしさ
の算出方法は（数１５）に限らない。

【０１６０】なお、実施例９のように、前記周波数帯域
分割フィルタ１００３はフーリエ変換を用いてもよく、
またエネルギーの算出手段１００６は実効値算出手段を
用いてもよい。

【０１６１】（実施例６）以下本発明の信号分離装置に
おける一実施例について、図面を参照しながら説明す
る。

【０１６２】第６の実施例における主たる構成図は図３
である。以上のように構成された雑音抑制装置の動作に
ついて説明する。信号の時間分割手段１０００は音響信
号を時間Ｔ毎にδｔの時間セグメントに分割する。入力
信号が音声信号であるのでδｔは１３ｍｓ、Ｔは２６ｍ
ｓとする。前記信号の時間分割手段１０００の形状はハ
ニング窓とすると前記信号の時間分割手段Ｗ（ｔ＋ｋ
Ｔ）は（数１）で表すことができる。

【０１６３】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とするとＳｉ
ｇ（ｔ＋ｋＴ）は（数２）で表される。

【０１６４】ピッチ抽出手段１００１はＳｉｇ（ｔ＋ｋ
Ｔ）のピッチを抽出する。ピッチはケプストラム法、自
己相関法などを用いて抽出される。第１の記憶手段１０
０２は前記ピッチ抽出手段１００１で求められたピッチ
ｆ_p（ｋＴ）を格納する。また、前記周波数帯域分割フ
ィルタ１００３はＳｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波
数帯域にＮ分割する。前記周波数帯域分割フィルタ１０
０３の出力をＳｉｇ（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ
＝１，２，・・・・Ｎ）。ｆ_mは帯域の中心周波数であ
る。Ｎ個の狭帯域信号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２
の記憶手段１００４に格納される。エネルギーの算出手
段１００５では各Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギ
ーを算出する。前記エネルギーの算出手段１００５のＮ
個の出力をＥｎｅ（ｋＴ，ｆ_m）と記述する。Ｅｎｅ
（ｋＴ，ｆ_m）は（数３）で算出される。

【０１６５】前記Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の
記憶手段１００６に格納される。特徴量の算出手段１０
０７はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。
本実施例では、特徴量として、ピッチ、オンセット、オ
ンセットの同期性、オフセット、オフセットの同期性、
連続性、調波性を算出する。

【０１６６】図６は本実施例における特徴量の算出手段
１００７の詳細図である。１３０１はピッチの評価測度
の算出手段、１３０２はオンセットの評価測度の算出手
段、１３０３はオフセットの評価測度の算出手段、１３
０４は連続性の評価測度の算出手段、１３０５は調波性
の評価測度の算出手段、１３０６はオンセットの同期性
の評価測度の算出手段、１３０７はオフセットの同期性
の評価測度の算出手段である。また、評価測度としてＳ
ｉｇ（ｔ＋ｋＴ，ｆ_m）で各特徴量が存在する確からし
さを用いる。

【０１６７】以下で前記特徴量の評価測度の算出手段に
ついて詳細に説明する。ピッチの評価測度の算出手段１
３０１において、ｆ_p（ｋＴ）に対しピッチが存在しう
る範囲を確率的に定義しＳｉｇ（ｋＴ，ｆ_m）がピッチ
である確からしさを算出する。図７はピッチが存在しう
る周波数の確率分布の一例である。ｆ_p（ｋＴ）に対し
てピッチが存在する確からしさがガウス分布をするとす
るとＳｉｇ（ｋＴ，ｆ_m）にピッチが存在する確からし
さＲｅｌ_{_pitch}（ｋＴ，ｆ_m）は（数９）によって与え
られる。

【０１６８】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【０１６９】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。オンセットの評価測度の算出手段１３０２ではこ
の変化が一つの音の時間変化ではなく図９（ｂ）のよう
に異なるもう一つの信号が加わったために起こったもの
であるという確からしさを人間の知覚に基づき定めこれ
をオンセットの確からしさとする。前記オンセットの確
からしさは前記第３の記憶手段１００６の出力によって
算出される。以下に算出方法の一つの例を挙げる。心理
学的に連続にきこえるか新たに音が加わったかという判
断基準が丁度変化する音の大きさが調べられている（19
94年日本音響学会秋期研究発表会稿論集2-7-18）。前記
エネルギーの算出手段１００５の出力を用いて現観測時
刻のエネルギーとそれより前の時刻までのエネルギーの
平均から前記オンセットの確からしさを定める。図１０
はオンセットの確からしさの算出方法のための説明図で
ある。前記前時刻までのエネルギーの平均は、現観測周
波数とその周波数に近接する２つの周波数帯域における
約２００ｍｓ区間（１６個のセグメント）のセグメント
の平均とする。Ｓｉｇ（ｋＴ，ｆ_m）の前記オンセット
の確からしさＲｅｌ_{_on}（Ｅ_Dif）は（数１０）で算出さ
れる。

【０１７０】図１１は（数１０）を図式化したものであ
る。前記オンセットの確からしさはピッチ帯域でオンセ
ットが検出されるまでピッチ帯域でのみ計算される。図
１２及び図１３は図８の信号に対する前記オンセットの
確からしさを説明する説明図である。図１２において時
刻ｋ＝２でピッチ帯域の前記オンセットの確からしさが
１．０となり閾値０．８を越えたので時刻ｋ＝２から全
周波数帯域で前記オンセットの確からしさの計算を始め
る。（ｋ，ｍ）＝（２，１）と（２，３）で前記オンセ
ットの確からしさ１．０が算出できたとする。この時、
図１３のように（ｋ，ｍ）＝（２，１）と（２，３）に
前記オンセットの確からしさ１．０を割り当てる。

【０１７１】オフセットの評価測度の算出手段１３０３
においてオフセットの確からしさが算出される。二つの
信号から一方の信号が除かれたとき音の大きさには変化
がある。そこで前記エネルギーの算出手段１００５の出
力を用いてオフセットを定める。前記オフセットの確か
らしさＲｅｌ_{_off}（Ｅ_Dif）は前記オンセットの確から
しさの算出方法と同様に（数１１）で与えられる。

【０１７２】オフセットの確からしさはオンセットが検
出されてからピッチ帯域でオフセットが検出されるまで
ピッチ周波数帯域のみ計算される。図１４及び図１５は
図８の信号に対する前記オフセットの確からしさを説明
する説明図である。図１４において時刻ｋ＝６でピッチ
帯域の前記オフセットの確からしさが１．０となりある
閾値０．９を越えたので時刻ｋ＝６において前記オフセ
ットの確からしさを全周波数帯域について算出する。
（ｋ，ｍ）＝（６，１）と（６，３）で前記オフセット
の確からしさ１．０を算出できたとすると図１６のよう
に前記オフセットの確からしさを割り当てる。

【０１７３】連続性の評価測度の算出手段１３０４で
は、音の大きさが時間的に変化したときそれが一つの信
号の時間変化である確からしさを指標とする。前記連続
性の確からしさは前記第３の記憶手段１００６の出力に
よって算出される。図１６は連続性の確からしさの算出
方法のための説明図である。図１６において、現観測時
刻を含む現観測時刻以前約２００ｍｓ区間（１６個のセ
グメント）の現観測セグメントのエネルギーの平均と現
観測周波数とその周波数に隣接する２つの周波数帯域に
おける現観測時刻を含まない現観測時刻より前の約２０
０ｍｓ区間（１６個のセグメント）のエネルギーの平均
から前記連続性の確からしさを決定する。前記連続性の
確からしさＲｅｌ_{_con}（Ｅ_Dif）は（数１２）で定義す
る。

【０１７４】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさは、オンセットが検出され
た後オフセットが検出されるまで計算される。図１８は
図８の信号に対する前記連続性の確からしさを説明する
説明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６で
オフセットが検出されたとき、時刻ｋ＝３から５までの
間前記連続性の確からしさが計算され、各（ｋ，ｍ）で
の連続性の確からしさが図１８のように算出される。

【０１７５】調波性の評価測度の算出手段１３０５で
は、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）がｆ_p（ｋＴ）に対して知
覚上調波関係である確からしさをもって調波性の確から
しさとする。ピッチの整数倍の周波数区間を中心に主観
的に調波性を感じる周波数分布を人間の心理測定関数か
ら定める。心理測定関数は電子情報通信学会論文誌'94/
5Vol.J77-ANo.5 731-740に掲載のデータを参考に決め
る。周波数帯域の中心周波数ｆ_mでのピッチｆ_p（ｋＴ）
のｐ倍（ｐは整数）の周波数からのずれｕが（数１３）
で求められる。

【０１７６】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【０１７７】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【０１７８】オンセットの同期性の評価測度の算出手段
１３０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5 731-740掲載のデー
タを参考に前記同期性の確からしさの分布Ｒｅｌ_sync
（ｋＴ，ｆ_m）を（数１５）で定義する。Ｓｉｇ（ｔ＋
ｋＴ，ｆ_m）でオンセットが検出されたときそれと最も
近いピッチのオンセット時刻の差をｖ（ｍｓ）とする
と、前記オンセットの同期性の確からしさＲｅｌ_sync
（ｋＴ，ｆ_m）は、（数１５）で記述できる。

【０１７９】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【０１８０】オフセットの同期性の評価測度の算出手段
１３０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【０１８１】図３の重み値の算出手段１００８では、評
価関数に各特徴（ピッチ、オンセット、オフセット、オ
ンセットの同期性、オフセットの同期性、調波性、連続
性）についての確からしさを入力し、Ｓｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算出する。

【０１８２】前記各特徴量のうち継時的確からしさに関
するものと周波数的確からしさに関するものに分類しそ
れぞれの確からしさの統合を行う。前記オンセットの確
からしさ、前記オンセットの同期性の確からしさ、前記
オフセットの確からしさ、前記オフセットの同期性の確
からしさ、前記パワーの連続性の確からしさを統合して
前記継時的確からしさとする。また、前記ピッチの確か
らしさ、前記ピッチの変動の範囲の確からしさ、前記調
波性の確からしさを統合して前記周波数的確からしさと
する。

【０１８３】図２９は評価値の計算に関する流れ図であ
る。オンセットが検出された時刻では前記継時的確から
しさは前記オンセットの確からしさと前記同期性の確か
らしさの積で与えられ（数２１）で算出される。

【０１８４】

【数２１】

【０１８５】オンセットが検出された後オフセットが検
出されるまでの時間では、前記継時的確からしさＲｅｌ
_{_t}（ｋＴ，ｆ_m）は前時刻の前記継時的確からしさと前
時刻との前記連続性の確からしさとの積で与えられ（数
２２）で算出される。

【０１８６】

【数２２】

【０１８７】図３０は前記継時的確からしさの分布の一
例である。また前記周波数的確からしさＲｅｌ_{_f}（ｋ
Ｔ，ｆ_m）は（数２３）で表現できる。

【０１８８】

【数２３】

【０１８９】前記継時的確からしさと前記周波数的確か
らしさの重視度をｇ_t（ｋＴ，ｆ_m），ｇ_f（ｋＴ，
ｆ_m）とする。

【０１９０】本実施例では評価関数にλ−ファジィ積分
を用いる。このとき、総合評価値ｗ（ｋＴ，ｆ_m）は
（数２４）で計算できる。

【０１９１】

【数２４】

【０１９２】図３１は総合評価値の算出結果の一例であ
る。ただし、ｇ_t（ｋＴ，ｆ_m）＝ｇ_f（ｋＴ，ｆ_m）、
λ＝０である。

【０１９３】図３の乗算手段１００９はＳｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）と前記評価値ｗ（ｋＴ，ｆ_m）を乗算する。前
記乗算手段１９０の出力は第４の記憶手段１０１０に格
納される。総和手段１０１１は全てのｋ，ｍについてｗ
（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を総和する。
前記総和手段１０１１の出力は雑音が抑制された音響信
号となっている。

【０１９４】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また、前記特徴量以外の特徴量を新たに組
み合わせてもよい。また、前記信号の時間分割手段１０
００の時間幅δｔや分割の時間周期Ｔは前記δｔ＝１３
ｍｓ、Ｔ＝２６ｍｓに限ったものではなく、形状もハニ
ング窓に限ったものではない。また、前記ピッチ抽出手
段１００１はケプストラム法、自己相関法に限ったもの
ではない。また、前記周波数帯域分割フィルタ１００３
の周波数帯域は２０Ｈｚに限ったものではない。また、
前記オンセットの確からしさの閾値は０．８に限らな
い。また、前記オフセットの確からしさの閾値は０．９
に限らない。また、前記オンセットの確からしさの算出
方法は（数７）に限らない。また、前記オフセットの確
からしさの算出方法は（数１１）に限らない。また、前
記連続性の確からしさの算出方法は（数１２）に限らな
い。また、前記調波性の確からしさの算出方法は（数１
３）および（数１４）に限らない。また、オンセットの
同期性の確からしさの算出方法やオフセットの同期性の
確からしさの算出方法は（数１５）に限らない。また、
総合評価値の算出方法はλ−ファジィ積分に限らず、菅
野積分、ショケ積分、逆菅野積分など他の手段でもよ
い。

【０１９５】なお、実施例９のように、前記周波数帯域
分割フィルタ１００３はフーリエ変換をもちいてもよ
く、またエネルギーの算出手段１００６は実効値算出手
段をもちいてもよい。

【０１９６】（実施例７）以下本発明の信号分離装置に
おける一実施例について、図面を参照しながら説明す
る。

【０１９７】第７の実施例における主たる構成図は図３
である。以上のように構成された雑音抑制装置の動作に
ついて説明する。信号の時間分割手段１０００は音響信
号を時間Ｔ毎にδｔの時間セグメントに分割する。入力
信号が音声信号であるので、δｔは１３ｍｓ、Ｔは２６
ｍｓとする。前記信号の時間分割手段１０００の形状は
ハニング窓とすると、前記信号の時間分割手段Ｗ（ｔ＋
ｋＴ）は（数１）で表すことができる。

【０１９８】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とすると、Ｓ
ｉｇ（ｔ＋ｋＴ）は（数２）で表される。

【０１９９】ピッチ抽出手段１００１は複数の抽出方法
を用いてＳｉｇ（ｔ＋ｋＴ）のピッチを抽出する。図２
４は前記ピッチ抽出手段１００１の詳細図である。図２
４において、３８１０は第１のピッチ抽出手段、３８２
０は第２のピッチ抽出手段、３８３０はピッチの推定値
ｆ_p（ｔ）の算出手段である。第１のピッチ抽出手段３
８１０によってピッチｆ_p1（ｋＴ）が、第２のピッチ抽
出手段３８２０によってピッチｆ_p2（ｋＴ）が計算され
たとする。このとき、ピッチの推定値ｆ_p（ｔ）の算出
手段３８３０によって、ピッチの推定値ｆ_p（ｋＴ）は
（数８）で求められる。

【０２００】図３において、第１の記憶手段１００２
は、前記ピッチ抽出手段１００１で求められたｆ_p1（ｋ
Ｔ）とｆ_p2（ｋＴ）とｆ_p（ｋＴ）を格納する。

【０２０１】また、周波数帯域分割フィルタ１００３は
Ｓｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯域にＮ分割す
る。周波数帯域分割フィルタ１００３の出力をＳｉｇ
（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝１，２，・・・・
Ｎ）。ｆ_mは帯域の中心周波数である。Ｎ個の狭帯域信
号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は、第２の記憶手段１００
４に格納される。エネルギー算出手段１００５では各Ｓ
ｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギーを算出する。前記
エネルギー算出手段１００５のＮ個の出力をＥｎｅ（ｋ
Ｔ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，ｆ_m）は（数３）
で算出される。

【０２０２】前記Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の
記憶手段１００６に格納される。特徴量の算出手段１０
０７はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。
本実施例では、特徴量として、ピッチ、オンセット、オ
ンセットの同期性、オフセット、オフセットの同期性、
連続性、調波性を算出する。

【０２０３】図６は図本実施例の特徴量の算出手段１０
０７の詳細図である。１３０１はピッチの評価測度の算
出手段、１３０２はオンセットの評価測度の算出手段、
１３０３はオフセットの評価測度の算出手段、１３０４
は連続性の評価測度の算出手段、１３０５は調波性の評
価測度の算出手段、１３０６はオンセットの同期性の評
価測度の算出手段、１３０７はオフセットの同期性の評
価測度の算出手段である。また、評価測度としてＳｉｇ
（ｔ＋ｋＴ，ｆ_m）で各特徴量が存在する確からしさを
用いる。以下で各特徴量の評価測度の算出手段について
詳細に説明する。

【０２０４】ピッチの評価測度の算出手段は前記ピッチ
抽出手段１００１で計算されたｆ_p1（ｋＴ）とｆ_p2（ｋ
Ｔ）の差からピッチｆ_p（ｋＴ）の確からしさを決定す
る。ｆ_p1（ｋＴ）とｆ_p2（ｋＴ）の差が０の時ピッチは
正確に抽出されたとしピッチｆ_p（ｋＴ）の確からしさ
は１とする。ｆ_p1（ｋＴ）とｆ_p2（ｋＴ）の差が大きく
なると確からしさは減少する。ピッチの確からしさが図
２５のようにｆ_p1（ｋＴ）−ｆ_p2（ｋＴ）の値によって
ガウス分布するとするとピッチの確からしさＲｅｌ（ｆ
_p1（ｋＴ）−ｆ_p2（ｋＴ））が（数１７）で決定され
る。

【０２０５】同時に、図２６のように、得られたピッチ
の推定値ｆ_p（ｋＴ）に対しあるピッチが存在しうる範
囲を確率的に定義しＳｉｇ（ｋＴ，ｆ_m）がピッチであ
る確からしさを算出する。ピッチの推定値ｆ_p（ｋＴ）
に対してピッチが存在する確からしさがガウス分布をす
るとするとＳｉｇ（ｋＴ，ｆ_m）にピッチが存在する確
からしさＲｅｌ_{_pitch}（ｋＴ，ｆ_m）は（数９）で与え
られる。

【０２０６】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【０２０７】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。前記オンセットの評価測度の算出手段１３０２で
はこの変化が一つの音の時間変化ではなく図９（ｂ）の
ように異なるもう一つの信号が加わったために起こった
ものであるという確からしさを人間の知覚に基づき定め
これをオンセットの確からしさとする。前記オンセット
の確からしさは前記第３の記憶手段１００６の出力によ
って算出される。以下に算出方法の一つの例を挙げる。
心理学的に連続にきこえるか新たに音が加わったかとい
う判断基準が丁度変化する音の大きさが調べられている
（1994年日本音響学会秋期研究発表会稿論集2-7-18）。
前記エネルギー算出手段１００５の出力を用いて現観測
時刻のエネルギーとそれより前の時刻までのエネルギー
の平均から前記オンセットの確からしさを定める。図１
０はオンセットの確からしさの算出方法のための説明図
である。前記前時刻までのエネルギーの平均は現観測周
波数とその周波数に近接する２つの周波数帯域における
約２００ｍｓ区間（１６個のセグメント）のセグメント
の平均とする。Ｓｉｇ（ｋＴ，ｆ_m）の前記オンセット
の確からしさＲｅｌ_{_on}（Ｅ_Dif）は（数１０）で算出さ
れる。

【０２０８】図１１は（数１０）を図式化したものであ
る。オンセットの確からしさはピッチ帯域でオンセット
が検出されるまでピッチ帯域でのみ計算される。図１２
及び図１３は図８の信号に対する前記オンセットの確か
らしさを説明する説明図である。図１２において時刻ｋ
＝２でピッチ帯域の前記オンセットの確からしさが１．
０となり閾値０．８を越えたので時刻ｋ＝２から全周波
数帯域で前記オンセットの確からしさの計算を始める。
（ｋ，ｍ）＝（２，１）と（２，３）で前記オンセット
の確からしさ１．０が算出できたとする。この時、図１
３のように（ｋ，ｍ）＝（２，１）と（２，３）に前記
オンセットの確からしさ１．０を割り当てる。

【０２０９】オフセットの評価測度の算出手段１３０３
においてオフセットの確からしさが算出される。二つの
信号から一方の信号が除かれたとき音の大きさには変化
がある。そこで前記エネルギーの算出手段１００５の出
力を用いてオフセットを定める。前記オフセットの確か
らしさＲｅｌ_{_off}（Ｅ_Dif）は前記オンセットの確から
しさの算出方法と同様に（数１１）で与えられる。

【０２１０】前記オフセットの確からしさはオンセット
が検出されてからピッチ帯域でオフセットが検出される
までピッチ周波数帯域のみ計算される。図１４及び図１
５は図８の信号に対する前記オフセットの確からしさを
説明する説明図である。図１４において時刻ｋ＝６でピ
ッチ帯域の前記オフセットの確からしさが１．０となり
ある閾値０．９を越えたので時刻ｋ＝６において前記オ
フセットの確からしさを全周波数帯域について算出す
る。（ｋ，ｍ）＝（６，１）と（６，３）で前記オフセ
ットの確からしさ１．０を算出できたとすると図１５の
ように前記オフセットの確からしさを割り当てる。

【０２１１】連続性の評価測度の算出手段１３０４で
は、音の大きさが時間的に変化したときそれが一つの信
号の時間変化である確からしさを指標とする。前記連続
性の確からしさは前記第３の記憶手段１００６の出力に
よって算出される。図１６は連続性の確からしさの算出
方法のための説明図である。図１６において、現観測時
刻を含む現観測時刻以前約２００ｍｓ区間（１６個のセ
グメント）の現観測セグメントのエネルギーの平均と現
観測周波数とその周波数に隣接する２つの周波数帯域に
おける現観測時刻を含まない現観測時刻より前の約２０
０ｍｓ区間（１６個のセグメント）のエネルギーの平均
から前記連続性の確からしさを決定する。前記連続性の
確からしさＲｅｌ_{_con}（Ｅ_Dif）は（数１２）で定義す
る。

【０２１２】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさは、オンセットが検出され
た後オフセットが検出されるまで計算される。図１８は
図８の信号に対する前記連続性の確からしさを説明する
説明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６で
オフセットが検出されたとき、時刻ｋ＝３から５までの
間前記連続性の確からしさが計算され、各（ｋ，ｍ）で
の連続性の確からしさが図１８のように算出される。

【０２１３】調波性の評価測度の算出手段１３０５で
は、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）がｆ_p（ｋＴ）に対して知
覚上調波関係である確からしさをもって調波性の確から
しさとする。ピッチの整数倍の周波数区間を中心に主観
的に調波性を感じる周波数分布を人間の心理測定関数か
ら定める。心理測定関数は電子情報通信学会論文誌'94/
5Vol.J77-ANo.5 731-740に掲載のデータを参考に決め
る。周波数帯域の中心周波数ｆ_mでの前記ピッチｆ_p（ｋ
Ｔ）のｐ倍（ｐは整数）の周波数からのずれｕが（数１
３）で求められる。

【０２１４】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【０２１５】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【０２１６】オンセットの同期性の評価測度の算出手段
１３０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5731-740に掲載のデー
タを参考に前記同期性の確からしさの分布Ｒｅｌ_sync
（ｋＴ，ｆ_m）を（数１５）で定義する。Ｓｉｇ（ｔ＋
ｋＴ，ｆ_m）でオンセットが検出されたときそれと最も
近いピッチのオンセット時刻の差をｖ（ｍｓ）とすると
前記オンセットの同期性の確からしさＲｅｌ_sync（ｋ
Ｔ，ｆ_m）は（数１５）で記述できる。

【０２１７】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【０２１８】オフセットの同期性の評価測度の算出手段
１３０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【０２１９】図３の重み値の算出手段１００８では、評
価関数に各特徴（ピッチ、オンセット、オフセット、オ
ンセットの同期性、オフセットの同期性、調波性、連続
性）についての確からしさを入力し、Ｓｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算出する。
前記各特徴量のうち継時的確からしさに関するものと周
波数的確からしさに関するものに分類しそれぞれの確か
らしさの統合を行う。前記オンセットの確からしさ、前
記オンセットの同期性の確からしさ、前記オフセットの
確からしさ、前記オフセットの同期性の確からしさ、前
記パワーの連続性の確からしさを統合して前記継時的確
からしさとする。また、前記ピッチの確からしさ、前記
ピッチの変動の範囲の確からしさ、前記調波性の確から
しさを統合して前記周波数的確からしさとする。

【０２２０】図２９は評価値の計算に関する流れ図であ
る。オンセットが検出された時刻では前記継時的確から
しさは前記オンセットの確からしさと前記同期性の確か
らしさの積で与えられ（数２１）で算出される。

【０２２１】オンセットが検出された後オフセットが検
出されるまでの時間では、前記継時的確からしさＲｅｌ
_{_t}（ｋＴ，ｆ_m）は前時刻の前記継時的確からしさと前
時刻との前記連続性の確からしさとの積で与えられ（数
２２）で算出される。

【０２２２】図３０は前記継時的確からしさの分布の一
例である。また前記周波数的確からしさＲｅｌ_{_f}（ｋ
Ｔ，ｆ_m）は（数２３）で表現できる。

【０２２３】前記継時的確からしさと前記周波数的確か
らしさの重視度をｇ_t（ｋＴ，ｆ_m），ｇ_f（ｋＴ，
ｆ_m）とする。

【０２２４】本実施例では、評価関数にλ−ファジィ積
分を用いる。このとき、総合評価値ｗ（ｋＴ，ｆ_m）は
（数２４）で計算できる。

【０２２５】図３１は総合評価値の算出結果の一例であ
る。ただし、ｇ_t（ｋＴ，ｆ_m）＝ｇ_f（ｋＴ，ｆ_m）、
λ＝０である。

【０２２６】図３の乗算手段１００９はＳｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）と前記評価値ｗ（ｋＴ，ｆ_m）を乗算する。前
記乗算手段１９０の出力は第４の記憶手段１０１０に格
納される。総和手段１０１１は全てのｋ，ｍについて、
ｗ（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を総和す
る。前記総和手段１０１１の出力は雑音が抑制された音
響信号となっている。

【０２２７】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また、前記特徴量以外の特徴量を新たに組
み合わせてもよい。また、前記信号の時間分割手段１０
００の時間幅δｔや分割の時間周期Ｔは前記δｔ＝１３
ｍｓ、Ｔ＝２６ｍｓに限ったものではなく、形状もハニ
ング窓に限ったものではない。また、前記ピッチ抽出手
段１００１は、ケプストラム法、自己相関法に限ったも
のではない。また、前記周波数帯域分割フィルタ１００
３の周波数帯域は２０Ｈｚに限ったものではない。ま
た、前記オンセットの確からしさの閾値は０．８に限ら
ない。また、前記オフセットの確からしさの閾値は０．
９に限らない。また、前記オンセットの確からしさの算
出方法は（数１０）に限らない。また、前記オフセット
の確からしさの算出方法（数１１）に限らない。また前
記連続性の確からしさの算出方法は（数１２）に限らな
い。また、前記調波性の確からしさの算出方法は（数１
３）（数１４）に限らない。またオンセットの同期性の
確からしさの算出方法やオフセットの同期性の確からし
さの算出方法は（数１５）に限らない。

【０２２８】なお、実施例９のように、前記周波数帯域
フィルタ１００３はフーリエ変換をもちいてもよく、ま
たエネルギーの算出手段１００６は実効値算出手段をも
ちいてもよい。

【０２２９】（実施例８）以下本発明の信号分離装置に
おける一実施例について、図面を参照しながら説明す
る。

【０２３０】第８の実施例における主たる構成図は図３
である。以上のように構成された雑音抑制装置の動作に
ついて説明する。信号の時間分割手段１０００は音響信
号を時間Ｔ毎にδｔの時間セグメントに分割する。入力
信号が音声信号であるので、δｔは１３ｍｓ、Ｔは２６
ｍｓとする。前記信号の時間分割手段１０００の形状は
ハニング窓とすると、前記信号の時間分割手段Ｗ（ｔ＋
ｋＴ）は（数１）で表すことができる。

【０２３１】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とすると、Ｓ
ｉｇ（ｔ＋ｋＴ）は（数２）で表される。

【０２３２】ピッチ抽出手段１００１は複数の抽出方法
を用いてＳｉｇ（ｔ＋ｋＴ）のピッチを抽出する。図５
は前記ピッチ抽出手段１００１の詳細図である。図５に
おいて、１２００は第１のピッチ抽出手段、１２０１は
第２のピッチ抽出手段、１２０２は第ｎのピッチ抽出手
段、１２０３はピッチの推定値ｆ_p（ｔ）の算出手段で
ある。第１のピッチ抽出手段１２００によってピッチｆ
_p1（ｋＴ）が、第２のピッチ抽出手段１２０１によって
ピッチｆ_p2（ｋＴ）が、第ｎのピッチ抽出手段１２０２
によってピッチｆ_pn（ｋＴ）が計算されたとする。この
とき、ピッチの推定値ｆ_p（ｔ）の算出手段１２０３に
よって、ピッチの推定値ｆ_p（ｋＴ）は（数８）で求め
られる。

【０２３３】また、ｆ_p1（ｋＴ）からｆ_pn（ｋＴ）の分
散Ｖ_pitch（ｋＴ）を（数１８）で求める。

【０２３４】図３において、第１の記憶手段１００２は
前記ピッチ抽出手段１００１で求められたｆ_p（ｋＴ）
と分散Ｖ_pitch（ｋＴ）を格納する。

【０２３５】また、周波数帯域分割フィルタ１００３は
Ｓｉｇ（ｔ＋ｋＴ）を２０Ｈｚの周波数帯域にＮ分割す
る。前記周波数帯域分割フィルタ１００３の出力をＳｉ
ｇ（ｔ＋ｋＴ，ｆ_m）と記述する（ｍ＝１，２，・・・
・Ｎ）。ｆ_mは帯域の中心周波数である。Ｎ個の狭帯域
信号Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）は第２の記憶手段１００
４に格納される。エネルギー算出手段１００５では各Ｓ
ｉｇ（ｔ＋ｋＴ，ｆ_m）のエネルギーを算出する。前記
エネルギー算出手段１００５のＮ個の出力をＥｎｅ（ｋ
Ｔ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，ｆ_m）は（数３）
で算出される。

【０２３６】前記Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の
記憶手段１００６に格納される。特徴量の算出手段１０
０７はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。
本実施例では特徴量としてピッチ、オンセット、オンセ
ットの同期性、オフセット、オフセットの同期性、連続
性、調波性を算出する。

【０２３７】図６は図３の特徴量の算出手段１００７の
詳細図である。１３０１はピッチの評価測度の算出手
段、１３０２はオンセットの評価測度の算出手段、１３
０３はオフセットの評価測度の算出手段、１３０４は連
続性の評価測度の算出手段、１３０５は調波性の評価測
度の算出手段、１３０６はオンセットの同期性の評価測
度の算出手段、１３０７はオフセットの同期性の評価測
度の算出手段である。また、評価測度としてＳｉｇ（ｔ
＋ｋＴ，ｆ_m）で各特徴量が存在する確からしさを用い
る。

【０２３８】以下で各特徴量の評価測度の算出手段につ
いて詳細に説明する。前記ピッチの評価測度の算出手段
は、前記ピッチ抽出手段１００１で計算された分散Ｖ
_pitch（ｋＴ）から推定されたピッチｆ_p（ｋＴ）の確か
らしさを決定する。Ｖ_pitch（ｋＴ）が０の時ピッチは
正確に抽出されたとしピッチｆ_p（ｋＴ）の確からしさ
は１とする。Ｖ_pitch（ｋＴ）が大きくなると確からし
さは減少する。ピッチの確からしさが図２７のようにＶ
_pitch（ｋＴ）の値によってガウス分布するとすると、
ピッチの確からしさＲｅｌ（Ｖ_pitch（ｋＴ））が（数
１９）で決定される。

【０２３９】同時に、図２８のように得られたピッチの
推定値ｆ_p（ｋＴ）に対しあるピッチが存在しうる範囲
を確率的に定義しＳｉｇ（ｋＴ，ｆ_m）がピッチである
確からしさを算出する。ピッチの推定値ｆ_p（ｋＴ）に
対してピッチが存在する確からしさがガウス分布をする
とするとＳｉｇ（ｋＴ，ｆ_m）にピッチが存在する確か
らしさＲｅｌ_{_pitch}（ｋＴ，ｆ_m）は（数２０）で与え
られる。

【０２４０】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【０２４１】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。オンセットの評価測度の算出手段１３０２ではこ
の変化が一つの音の時間変化ではなく図９（ｂ）のよう
に異なるもう一つの信号が加わったために起こったもの
であるという確からしさを人間の知覚に基づき定めこれ
をオンセットの確からしさとする。前記オンセットの確
からしさは前記第３の記憶手段１００６の出力によって
算出される。以下に算出方法の一つの例を挙げる。心理
学的に連続にきこえるか新たに音が加わったかという判
断基準が丁度変化する音の大きさが調べられている（19
94年日本音響学会秋期研究発表会稿論集2-7-18）。前記
エネルギー算出手段１００５の出力を用いて現観測時刻
のエネルギーとそれより前の時刻までのエネルギーの平
均から前記オンセットの確からしさを定める。図１０は
オンセットの確からしさの算出方法のための説明図であ
る。前記前時刻までのエネルギーの平均は現観測周波数
とその周波数に近接する２つの周波数帯域における約２
００ｍｓ区間（１６個のセグメント）のセグメントの平
均とする。Ｓｉｇ（ｋＴ，ｆ_m）の前記オンセットの確
からしさＲｅｌ_{_on}（Ｅ_Dif）は（数１０）で算出され
る。

【０２４２】図１１は（数１０）を図式化したものであ
る。前記オンセットの確からしさはピッチ帯域でオンセ
ットが検出されるまでピッチ帯域でのみ計算される。図
１２及び図１３は図８の信号に対する前記オンセットの
確からしさを説明する説明図である。図１２において時
刻ｋ＝２でピッチ帯域の前記オンセットの確からしさが
１．０となり閾値０．８を越えたので時刻ｋ＝２から全
周波数帯域で前記オンセットの確からしさの計算を始め
る。（ｋ，ｍ）＝（２，１）と（２，３）で前記オンセ
ットの確からしさ１．０が算出できたとする。この時、
図１３のように（ｋ，ｍ）＝（２，１）と（２，３）に
前記オンセットの確からしさ１．０を割り当てる。

【０２４３】オフセットの評価測度の算出手段１３０３
においてオフセットの確からしさが算出される。二つの
信号から一方の信号が除かれたとき音の大きさには変化
がある。そこで前記エネルギーの算出手段１００５の出
力を用いてオフセットを定める。前記オフセットの確か
らしさＲｅｌ_{_off}（Ｅ_Dif）は前記オンセットの確から
しさの算出方法と同様に（数１１）で与えられる。

【０２４４】前記オフセットの確からしさはオンセット
が検出されてからピッチ帯域でオフセットが検出される
までピッチ周波数帯域のみ計算される。図２１４及び図
１５は図８の信号に対する前記オフセットの確からしさ
を説明する説明図である。

【０２４５】図１４において時刻ｋ＝６でピッチ帯域の
前記オフセットの確からしさが１．０となりある閾値
０．９を越えたので時刻ｋ＝６において前記オフセット
の確からしさを全周波数帯域について算出する。（ｋ，
ｍ）＝（６，１）と（６，３）で前記オフセットの確か
らしさ１．０を算出できたとすると図１５のように前記
オフセットの確からしさを割り当てる。

【０２４６】連続性の評価測度の算出手段１３０４で
は、音の大きさが時間的に変化したときそれが一つの信
号の時間変化である確からしさを指標とする。前記連続
性の確からしさは前記第３の記憶手段１００６の出力に
よって算出される。図１６は連続性の確からしさの算出
方法のための説明図である。図１６において、現観測時
刻を含む現観測時刻以前約２００ｍｓ区間（１６個のセ
グメント）の現観測セグメントのエネルギーの平均と現
観測周波数とその周波数に隣接する２つの周波数帯域に
おける現観測時刻を含まない現観測時刻より前の約２０
０ｍｓ区間（１６個のセグメント）のエネルギーの平均
から前記連続性の確からしさを決定する。前記連続性の
確からしさＲｅｌ_{_con}（Ｅ_Dif）は（数１２）で定義す
る。

【０２４７】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさは、オンセットが検出され
た後オフセットが検出されるまで計算される。図１８は
図８の信号に対する前記連続性の確からしさを説明する
説明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６で
オフセットが検出されたとき、時刻ｋ＝３から５までの
間前記連続性の確からしさが計算され、各（ｋ，ｍ）で
の連続性の確からしさが図１８のように算出される。

【０２４８】調波性の評価測度の算出手段１３０５で
は、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）が前記ｆ_p（ｋＴ）に対し
て知覚上調波関係である確からしさをもって調波性の確
からしさとする。ピッチの整数倍の周波数区間を中心に
主観的に調波性を感じる周波数分布を人間の心理測定関
数から定める。心理測定関数は電子情報通信学会論文
誌'94/5Vol.J77-ANo.5 731-740に掲載のデータを参考
に決める。周波数帯域の中心周波数ｆ_mでの前記ピッチ
ｆ_p（ｋＴ）のｐ倍（ｐは整数）の周波数からのずれｕ
が（数１３）で求められる。

【０２４９】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【０２５０】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【０２５１】オンセットの同期性の評価測度の算出手段
１３０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5 731-740に掲載のデ
ータを参考に前記同期性の確からしさの分布Ｒｅｌ_syn
c（ｋＴ，ｆ_m）を（数１５）で定義する。Ｓｉｇ（ｔ＋
ｋＴ，ｆ_m）でオンセットが検出されたときそれと最も
近いピッチのオンセット時刻の差をｖ（ｍｓ）とする
と、前記オンセットの同期性の確からしさＲｅｌ_sync
（ｋＴ，ｆ_m）は、（数１５）で記述できる。

【０２５２】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【０２５３】オフセットの同期性の評価測度の算出手段
１３０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【０２５４】図３の重み値の算出手段１００８では、評
価関数に各特徴（ピッチ、オンセット、オフセット、オ
ンセットの同期性、オフセットの同期性、調波性、連続
性）についての確からしさを入力し、Ｓｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算出する。

【０２５５】各特徴量のうち継時的確からしさに関する
ものと周波数的確からしさに関するものに分類しそれぞ
れの確からしさの統合を行う。前記オンセットの確から
しさ、前記オンセットの同期性の確からしさ、前記オフ
セットの確からしさ、前記オフセットの同期性の確から
しさ、前記パワーの連続性の確からしさを統合して前記
継時的確からしさとする。また、前記ピッチの確からし
さ、前記ピッチの変動の範囲の確からしさ、前記調波性
の確からしさを統合して前記周波数的確からしさとす
る。

【０２５６】図２９は評価値の計算に関する流れ図であ
る。オンセットが検出された時刻では前記継時的確から
しさは前記オンセットの確からしさと前記同期性の確か
らしさの積で与えられ（数２１）で算出される。

【０２５７】オンセットが検出された後オフセットが検
出されるまでの時間では前記継時的確からしさＲｅｌ_{_t}
（ｋＴ，ｆ_m）は前時刻の前記継時的確からしさと前時
刻との前記連続性の確からしさとの積で与えられ（数２
２）で算出される。

【０２５８】図４４００は前記継時的確からしさの分布
の一例である。また前記周波数的確からしさＲｅｌ
_{_f}（ｋＴ，ｆ_m）は（数２３）で表現できる。

【０２５９】前記継時的確からしさと前記周波数的確か
らしさの重視度をｇ_t（ｋＴ，ｆ_m），ｇ_f（ｋＴ，
ｆ_m）とする。

【０２６０】本実施例では、評価関数にλ−ファジィ積
分を用いる。このとき、総合評価値ｗ（ｋＴ，ｆ_m）は
（数２４）で計算できる。

【０２６１】図３１は総合評価値の算出結果の一例であ
る。ただし、ｇ_t（ｋＴ，ｆ_m）＝ｇ_f（ｋＴ，ｆ_m）、
λ＝０である。

【０２６２】図３の乗算手段１００９はＳｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）と前記評価値ｗ（ｋＴ，ｆ_m）を乗算する。前
記乗算手段１００９の出力は第４の記憶手段１０１０に
格納される。総和手段１０１１は全てのｋ，ｍについ
て、ｗ（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）を総和
する。前記総和手段１０１１の出力は雑音が抑制された
音響信号となっている。

【０２６３】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また、前記特徴量以外の特徴量を新たに組
み合わせてもよい。また、前記ピッチ抽出手段１００１
はケプストラム法、自己相関法に限ったものではない。
また、前記周波数帯域分割フィルタ１００３の周波数帯
域は２０Ｈｚに限ったものではない。また、前記オンセ
ットの確からしさの閾値は０．８に限らない。また、前
記オフセットの確からしさの閾値は０．９に限らない。
また、前記オンセットの確からしさの算出方法は（数１
０）に限らない。また、前記オフセットの確からしさの
算出方法は（数１１）に限らない。また、前記調波性の
確からしさの算出方法は（数１３）（数１４）に限らな
い。また、オンセットの同期性の確からしさの算出方法
やオフセットの同期性の確からしさの算出方法は（数１
５）に限らない。また、ピッチが存在する確からしさの
算出においてガウス分布以外の分布で定義してもよい。

【０２６４】なお、実施例９のように、前記周波数帯域
フィルタ１００３はフーリエ変換をもちいてもよく、ま
たエネルギーの算出手段１００６は実効値算出手段を用
いてもよい。

【０２６５】（実施例９）以下本発明の信号分離装置に
おける一実施例について、図面を参照しながら説明す
る。

【０２６６】図３２は本発明の第９の実施例における雑
音抑制装置の構成を示す構成図である。図３２におい
て、９０００は信号の時間切り出しを行う信号の時間分
割手段、９００１はピッチ抽出手段、９００２は第１の
記憶手段、９００３はフーリエ変換手段、９００４は第
２の記憶手段、９００５は実効値算出手段、９００６は
第３の記憶手段、９００７は特徴量の算出手段、９００
８は重み値の算出手段、９００９は乗算手段、９０１０
はフーリエ逆変換手段、９０１１は第４の記憶手段、９
０１２は総和手段である。

【０２６７】以上のように構成された雑音抑制装置の動
作について説明する。信号の時間分割手段９０００は音
響信号を時間Ｔ毎にδｔの時間セグメントに分割する。
入力信号が音声信号であるので、δｔは１３ｍｓ、Ｔは
２６ｍｓとする。前記信号の時間分割手段９０００の形
状はハニング窓とすると、前記信号の時間分割手段Ｗ
（ｔ＋ｋＴ）は（数１）で表すことができる。

【０２６８】前記信号の時間分割手段Ｗ（ｔ＋ｋＴ）の
出力信号をＳｉｇ（ｔ＋ｋＴ）（ｔ＝０〜δｔ）とす
る。入力音声をＳｐｅｅｃｈ（ｔ＋ｋＴ）とすると、Ｓ
ｉｇ（ｔ＋ｋＴ）は（数２）で表される。

【０２６９】ピッチ抽出手段９００１はＳｉｇ（ｔ＋ｋ
Ｔ）のピッチを抽出する。ピッチは、ケプストラム法、
自己相関法などを用いて抽出される。第１の記憶手段９
００２は前記ピッチ抽出手段９００１で求められたピッ
チｆ_p（ｋＴ）を格納する。また、フーリエ変換手段９
００３はＳｉｇ（ｔ＋ｋＴ）をフーリエ変換する。この
時の周波数分解能は２０Ｈｚとする。前記フーリエ変換
手段９００３の出力をＳｉｇ（ｔ＋ｋＴ，ｆ_m）と記述
する（ｍ＝１，２，・・・・Ｎ）。ｆ_mは周波数であ
る。Ｎ個の周波数サンプル点Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）
は、第２の記憶手段９００４に格納される。実効値算出
手段９００５では、各Ｓｉｇ（ｔ＋ｋＴ，ｆ _m）の実効
値を算出する。前記実効値算出手段９００５のＮ個の出
力をＥｎｅ（ｋＴ，ｆ_m）と記述する。Ｅｎｅ（ｋＴ，
ｆ_m）は（数２５）で算出される。

【０２７０】

【数２５】

【０２７１】前記Ｎ個のＥｎｅ（ｋＴ，ｆ_m）は第３の
記憶手段９００６に格納される。特徴量の算出手段９０
０７はＳｉｇ（ｔ＋ｋＴ，ｆ_m）の特徴量を算出する。
本実施例では特徴量としてピッチ、オンセット、オンセ
ットの同期性、オフセット、オフセットの同期性、連続
性、調波性を算出する。

【０２７２】

【数２６】

【０２７３】図５は本実施例の特徴量の算出手段９００
７の詳細図である。１２０１はピッチの評価測度の算出
手段、１２０２はオンセットの評価測度の算出手段、１
２０３はオフセットの評価測度の算出手段、１２０４は
連続性の評価測度の算出手段、１２０５は調波性の評価
測度の算出手段、１２０６はオンセットの同期性の評価
測度の算出手段、１２０７はオフセットの同期性の評価
測度の算出手段、また、評価測度としてＳｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）で各特徴量が存在する確からしさを用いる。
以下で各特徴量の評価測度の算出手段について詳細に説
明する。

【０２７４】ピッチの評価測度の算出手段１２０１にお
いて、ｆ_p（ｋＴ）に対しピッチが存在しうる範囲を確
率的に定義しＳｉｇ（ｋＴ，ｆ_m）がピッチである確か
らしさを算出する。図７はピッチが存在しうる周波数の
確率分布の一例である。前記ｆ_p（ｋＴ）に対してピッ
チが存在する確からしさがガウス分布をするとすると前
記Ｓｉｇ（ｋＴ，ｆ_m）にピッチが存在する確からしさR
el_pitch（ｋＴ，ｆｍ）は（数９）によって与えられ
る。

【０２７５】図８は２つの成分から成る調波複合音が時
刻ｋ＝２からｋ＝６まで持続した時のスペクトル図であ
る。この信号に関して以下の各確からしさの算出手段に
ついて説明する。

【０２７６】図９は音の大きさの変化に対する聴覚知覚
に関する模式図である。図９（ａ）のように一つの信号
にもう一つの信号が加わったとき音の大きさには変化が
ある。前記オンセットの評価測度の算出手段１２０２で
はこの変化が一つの音の時間変化ではなく図９（ｂ）の
ように異なるもう一つの信号が加わったために起こった
ものであるという確からしさを人間の知覚に基づき定め
これをオンセットの確からしさとする。前記オンセット
の確からしさは第３の記憶手段９００６の出力によって
算出される。以下に算出方法の一つの例を挙げる。心理
学的に連続にきこえるか新たに音が加わったかという判
断基準が丁度変化する音の大きさが調べられている（19
94年日本音響学会秋期研究発表会稿論集2-7-18）。実効
値算出手段９００５の出力を用いて現観測時刻の実効値
とそれより前の時刻までの実効値の平均から前記オンセ
ットの確からしさを定める。図１０はオンセットの確か
らしさの算出方法のための説明図である。前記前時刻ま
での実効値の平均は、現観測周波数とその周波数に近接
する２つの周波数帯域における約２００ｍｓ区間（１６
個のセグメント）のセグメントの平均とする。Ｓｉｇ
（ｋＴ，ｆ_m）の前記オンセットの確からしさＲｅｌ_{_on}
（Ｅ_Dif）は（数１０）で算出される。

【０２７７】図１１は（数１０）を図式化したものであ
る。オンセットの確からしさはピッチ帯域でオンセット
が検出されるまでピッチ帯域でのみ計算される。図１２
及び図１３は図８の信号に対する前記オンセットの確か
らしさを説明する説明図である。図１２において時刻ｋ
＝２でピッチ帯域の前記オンセットの確からしさが１．
０となり閾値０．８を越えたので時刻ｋ＝２から全周波
数帯域で前記オンセットの確からしさの計算を始める。
（ｋ，ｍ）＝（２，１）と（２，３）で前記オンセット
の確からしさ１．０が算出できたとする。この時図２７
００のように（ｋ，ｍ）＝（２，１）と（２，３）に前
記オンセットの確からしさ１．０を割り当てる。

【０２７８】オフセットの評価測度の算出手段１２０３
においてオフセットの確からしさが算出される。二つの
信号から一方の信号が除かれたとき音の大きさには変化
がある。そこで前記実効値の算出手段９００５の出力を
用いてオフセットを定める。前記オフセットの確からし
さRel__off（E_Dif）は前記オンセットの確からしさの算
出方法と同様に（数１１）で与えられる。

【０２７９】前記オフセットの確からしさはオンセット
が検出されてからピッチ帯域でオフセットが検出される
までピッチ周波数帯域のみ計算される。図１４及び図１
５は図８の信号に対する前記オフセットの確からしさを
説明する説明図である。図１４において時刻ｋ＝６でピ
ッチ帯域の前記オフセットの確からしさが１．０となり
ある閾値０．９を越えたので時刻ｋ＝６において前記オ
フセットの確からしさを全周波数帯域について算出す
る。（ｋ，ｍ）＝（６，１）と（６，３）で前記オフセ
ットの確からしさ１．０を算出できたとすると図１５の
ように前記オフセットの確からしさを割り当てる。

【０２８０】連続性の評価測度の算出手段１２０４で
は、音の大きさが時間的に変化したときそれが一つの信
号の時間変化である確からしさを指標とする。前記連続
性の確からしさは前記第３の記憶手段９００６の出力に
よって算出される。図１６は連続性の確からしさの算出
方法のための説明図である。図１６において、現観測時
刻を含む現観測時刻以前約２００ｍｓ区間（１６個のセ
グメント）の現観測セグメントの実効値の平均と現観測
周波数とその周波数に隣接する２つの周波数帯域におけ
る現観測時刻を含まない現観測時刻より前の約２００ｍ
ｓ区間（１６個のセグメント）の実効値の平均から前記
連続性の確からしさを決定する。前記連続性の確からし
さRel_{_con}(E_Dif)は（数１２）で定義する。

【０２８１】図１７は（数１２）を図式化したものであ
る。前記連続性の確からしさは、オンセットが検出され
た後オフセットが検出されるまで計算される。図１８は
図８の信号に対する前記連続性の確からしさを説明する
説明図である。時刻ｋ＝２でオンセットが時刻ｋ＝６で
オフセットが検出されたとき、時刻ｋ＝３から５までの
間前記連続性の確からしさが計算され、各（ｋ，ｍ）で
の連続性の確からしさが図１８のように算出される。

【０２８２】調波性の評価測度の算出手段１２０５で
は、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）がｆ_p（ｋＴ）に対して知
覚上調波関係である確からしさをもって調波性の確から
しさとする。ピッチの整数倍の周波数区間を中心に主観
的に調波性を感じる周波数分布を人間の心理測定関数か
ら定める。心理測定関数は電子情報通信学会論文誌'94/
5Vol.J77-ANo.5731-740に掲載のデータを参考に決め
る。周波数帯域の中心周波数ｆ_mでの前記ピッチｆ_p（ｋ
Ｔ）のｐ倍（ｐは整数）の周波数からのずれｕが（数１
３）で求められる。

【０２８３】前記調波性の確からしさの分布関数Ｒｅｌ
_harm（ｕ）はｕを用いて（数１４）のように定義す
る。

【０２８４】図１９は（数１４）を図式化したものであ
る。図２０は図８の信号に対する前記調波性の確からし
さを説明する説明図である。図２０ではピッチの確から
しさを１としている。

【０２８５】オンセットの同期性の評価測度の算出手段
１２０６では、Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）でのオンセッ
トがｆ_p（ｋＴ）のオンセット時刻を基準に主観的に同
時であると判断される時間分布に基づき同期性の確から
しさとする。知覚上のオンセットの同期性は電子情報通
信学会論文誌'94/5Vol.J77-ANo.5731-740に掲載のデー
タを参考に前記同期性の確からしさの分布Ｒｅｌ_sync
（ｋＴ，ｆ_m）を（数１５）で定義する。Ｓｉｇ（ｔ＋
ｋＴ，ｆ_m）でオンセットが検出されたときそれと最も
近いピッチのオンセット時刻の差をｖ（ｍｓ）とする
と、前記オンセットの同期性の確からしさＲｅｌ_sync
（ｋＴ，ｆ_m）は、（数１５）で記述できる。

【０２８６】図２１は（数１５）を図式化したものであ
る。図２２は図８の信号に対する前記オンセットの同期
性の確からしさを説明する説明図である。

【０２８７】オフセットの同期性の評価測度の算出手段
１２０７は前記オンセットの同期性の評価測度と同様な
方法により算出できる。Ｓｉｇ（ｔ＋ｋＴ，ｆ_m）での
オフセットがｆ_p（ｋＴ）のオフセット時刻を基準に主
観的に同時であると判断される時間分布に基づき前記オ
フセットの同期性の確からしさとする。図２３は図８の
信号に対する前記オフセットの同期性の確からしさを説
明する説明図である。

【０２８８】図３２の重み値の算出手段９００８では、
評価関数に前記各特徴（ピッチ、オンセット、オフセッ
ト、オンセットの同期性、オフセットの同期性、調波
性、連続性）についての確からしさを入力し、Ｓｉｇ
（ｔ＋ｋＴ，ｆ_m）での重み値ｗ（ｔ＋ｋＴ,ｆ_m）を算
出する。

【０２８９】各特徴量のうち継時的確からしさに関する
ものと周波数的確からしさに関するものに分類しそれぞ
れの確からしさの統合を行う。前記オンセットの確から
しさ、前記オンセットの同期性の確からしさ、前記オフ
セットの確からしさ、前記オフセットの同期性の確から
しさ、前記パワーの連続性の確からしさを統合して、前
記継時的確からしさとする。また、前記ピッチの確から
しさ、前記ピッチの変動の範囲の確からしさ、前記調波
性の確からしさを統合して、前記周波数的確からしさと
する。

【０２９０】図２９は評価値の計算に関する流れ図であ
る。オンセットが検出された時刻では前記継時的確から
しさは前記オンセットの確からしさと前記同期性の確か
らしさの積で与えられ（数２１）で算出される。

【０２９１】オンセットが検出された後オフセットが検
出されるまでの時間では、前記継時的確からしさＲｅｌ
_{_t}（ｋＴ，ｆ_m）は前時刻の前記継時的確からしさと前
時刻との前記連続性の確からしさとの積で与えられ、
（数２２）で算出される。

【０２９２】図３０は前記継時的確からしさの分布の一
例である。また前記周波数的確からしさＲｅｌ_{_f}（ｋ
Ｔ，ｆ_m）は、（数２３）で表現できる。

【０２９３】前記継時的確からしさと、前記周波数的確
からしさの重視度をｇ_t（ｋＴ，ｆ _m），ｇ_f（ｋＴ，
ｆ_m）とする。

【０２９４】本実施例では、評価関数にλ−ファジィ積
分を用いる。このとき、総合評価値ｗ（ｋＴ，ｆ_m）は
（数２４）で計算できる。

【０２９５】図３１は総合評価値の算出結果の一例であ
る。ただし、ｇ_t（ｋＴ，ｆ_m）＝ｇ_f（ｋＴ，ｆ_m）、
λ＝０である。

【０２９６】図３２の乗算手段９００９はＳｉｇ（ｔ＋
ｋＴ，ｆ_m）と前記評価値ｗ（ｋＴ，ｆ_m）を乗算する。
前記乗算手段１９０の出力は、第４の記憶手段９０１０
に格納される。フーリエ逆変換手段９０１１は全ての
ｋ，ｍについて、ｗ（ｋＴ，ｆ_m）×Ｓｉｇ（ｔ＋ｋ
Ｔ，ｆ_m）をフーリエ逆変換する。前記フーリエ逆変換
手段９０１１の出力は雑音が抑制された音響信号となっ
ている。

【０２９７】なお、前記各特徴量のうち前記調波性かま
たは前記連続性が含まれていれば全ての特徴量を用いる
必要はない。また、前記特徴量以外の特徴量を新たに組
み合わせてもよい。また。前記信号の時間分割手段９０
００の時間幅δｔや分割の時間周期Ｔは前記δｔ＝１３
ｍｓ、Ｔ＝２６ｍｓに限ったものではなく、形状もハニ
ング窓に限ったものではない。また、前記ピッチ抽出手
段９００１はケプストラム法、自己相関法に限ったもの
ではない。また、前記フーリエ変換手段９００３と前記
フーリエ逆変換手段９０１０の周波数分解能は２０Ｈｚ
に限ったものではない。また、前記オンセットの確から
しさの閾値は０．８に限らない。また、前記オンセット
の確からしさの閾値は０．９に限らない。また、前記オ
ンセットの確からしさの算出方法は（数１０）に限らな
い。また、前記オフセットの確からしさの算出方法は
（数１１）に限らない。また、前記連続性の確からしさ
の算出方法は（数１２）に限らない。また、前記調波性
の確からしさの算出方法は（数１３）および（数１４）
に限らない。また、オンセットの同期性の確からしさの
算出方法やオフセットの同期性の確からしさの算出方法
は（数１５）に限らない。また、総合評価値の算出方法
はλ−ファジィ積分に限らず、菅野積分、ショケ積分、
逆菅野積分など他の手段でもよい。

【０２９８】図３３は前記構成によって得られた出力信
号のＳ／Ｎ比と残差を示す図である。図３３の（Ａ）は
入力信号のＳ／Ｎ比と残差、（Ｂ）はλ＝−０．９のと
きのＳ／Ｎ比と残差、（Ｃ）はλ＝０．９のときのＳ／
Ｎ比と残差である。図１００００のように（Ｂ）（Ｃ）
とも（Ａ）に比べてＳ／Ｎ比が増加し残差が減少してお
り、前記構成の効果が確認できる。

【０２９９】

【発明の効果】以下のように本発明は、信号の時間分割
手段、ピッチ抽出手段、第１の記憶手段、周波数帯域分
割手段、第２の記憶手段、エネルギーまたは実効値の算
出手段、第３の記憶手段、特徴量の算出手段、重み値の
算出手段、乗算手段、第４の記憶手段、総和手段を設け
ることにより、人間の聴覚が聴感上の雑音抑制のためい
利用している音の物理的な特徴を用いて入力信号の周波
数の重みづけを行うことにより雑音が混在した音響信号
から雑音を抑制することにより目的の音響信号を抽出こ
とができる優れた雑音抑制装置を実現するものである。

【図面の簡単な説明】

【図１】従来例における雑音抑制装置の構成を示すブロ
ック図

【図２】本発明の雑音抑制装置の構成図

【図３】本発明の雑音抑制装置の実施例１及び実施例２
および実施例３および実施例４および実施例５および実
施例６および実施例７および実施例８における構成図

【図４】本発明の雑音抑制装置の実施例１および実施例
２における特徴量の算出手段の一例を示す図

【図５】本発明の雑音抑制装置の実施例２および実施例
５および実施例８および実施例９におけるピッチ抽出手
段の一例を示す図

【図６】本発明の雑音抑制装置の実施例３および実施例
４および実施例５および実施例６および実施例７および
実施例８における特徴量の算出手段の一例を示す図

【図７】本発明の雑音抑制装置の実施例３および実施例
６および実施例９におけるピッチが存在する周波数の確
率分布の一例を示す図

【図８】本発明の雑音抑制装置の実施例３および実施例
４およぴ実施例５および実施例６および実施例７および
実施例８および実施例９における動作の説明のための図

【図９】本発明の雑音抑制装置の実施例３および実施例
４および実施例５および実施例６および実施例７および
実施例８および実施例９においてオンセットである判断
の基準を説明するための音の大きさの変化に対する聴覚
知覚に関する模式図

【図１０】本発明の雑音抑制装置の実施例３および実施
例５および実施例６および実施例７および実施例８およ
び実施例９において、オンセットの確からしさの算出方
法の説明図

【図１１】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における（数１０）を図式化
した図

【図１２】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図１３】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図１４】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図１５】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図１６】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９において、連続性の確からし
さの算出方法の説明図

【図１７】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における（数１２）を図式化
した図

【図１８】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図１９】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における（数１４）を図式化
した図

【図２０】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図２１】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における（数１５）を図式化
した図

【図２２】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図２３】本発明の雑音抑制装置の実施例３および実施
例４および実施例５および実施例６および実施例７およ
び実施例８および実施例９における動作の説明のための
図

【図２４】本発明の雑音抑制装置の実施例４および実施
例７におけるピッチ抽出手段の一例を示す図

【図２５】本発明の雑音抑制装置の実施例４および実施
例７におけるピッチの確からしさの分布の一例を示す図

【図２６】本発明の雑音抑制装置の実施例４および実施
例７において、ピッチの推定値に対しピッチが存在する
確率分布の一例を示す図

【図２７】本発明の雑音抑制装置の実施例５および実施
例８において、ピッチの推定値に対する確からしさの確
率分布の一例を示す図

【図２８】本発明の雑音抑制装置の実施例５および実施
例８において、ピッチの推定値に対しピッチが存在する
確率分布の一例を示す図

【図２９】本発明の雑音抑制装置の実施例６および実施
例７および実施例８および実施例９において、評価値の
計算の一例を示す流れ図

【図３０】本発明の雑音抑制装置の実施例６および実施
例７および実施例８および実施例９における動作の説明
のための図

【図３１】本発明の雑音抑制装置の実施例６および実施
例７および実施例８および実施例９における動作の説明
のための図

【図３２】本発明の実施例９における雑音抑制装置の構
成を示す構成図

【図３３】本発明の実施例９における雑音抑制装置によ
る効果を示す図

【符号の説明】

１ 音声入力用の接話マイク ２ 雑音成分信号入力用のセンサマイク ３ ローパスフィルタ ４ ローパスフィルタ ５ バンドパスフィルタバンク １０ｉから１０Ｎ 雑音除去回路 １１ｉから１１Ｎ 位相差検出補正回路 １２ｉから１２Ｎ レベル差検出補正回路 １００ 信号の時間分割手段 ２００ ピッチ抽出手段 ３００ 第１の記憶手段 ４００ 周波数帯域分割手段 ５００ 第２の記憶手段 ６００ エネルギーまたは実効値の算出手段 ７００ 第３の記憶手段 ８００ 特徴量の算出手段 ９００ 重み値の算出手段 ９１０ 乗算手段 ９２０ 第４の記憶手段 ９３０ 総和手段 ３８１０ 第１のピッチ抽出手段 ３８２０ 第２のピッチ抽出手段 ３８３０ ピッチの推定値の算出手段 １０００ 信号の時間分割手段 １００１ ピッチ抽出手段 １００２ 第１の記憶手段 １００３ 周波数帯域分割フィルタ １００４ 第２の記憶手段 １００５ エネルギーの算出手段 １００６ 第３の記憶手段 １００７ 特徴量の算出手段 １００８ 重み値の算出手段 １００９ 乗算手段 １０１０ 第４の記憶手段 １０１１ 総和手段 １１０１ オンセットの算出手段 １１０２ オフセットの算出手段 １１０３ 連続性の算出手段 １１０４ 調波性の算出手段 １１０５ オンセットの同期性の算出手段 １１０６ オフセットの同期性の算出手段 １２００ 第１のピッチ抽出手段 １２０１ 第２のピッチ抽出手段 １２０２ 第ｎのピッチ抽出手段 １２０３ ピッチの推定値の算出手段 １３０１ ピッチの評価測度の算出手段 １３０２ オンセットの評価測度の算出手段 １３０３ オフセットの評価測度の算出手段 １３０４ 連続性の評価測度の算出手段 １３０５ 調波性の評価測度の算出手段 １３０６ オンセットの同期性の評価測度の算出手段 １３０７ オフセットの同期性の評価測度の算出手段 ９０００ 信号の時間分割手段 ９００１ ピッチ抽出手段 ９００２ 第１の記憶手段 ９００３ フーリエ変換手段 ９００４ 第２の記憶手段 ９００５ 実効値算出手段 ９００６ 第３の記憶手段 ９００７ 特徴量の算出手段 ９００８ 重み値の算出手段 ９００９ 乗算手段 ９０１０ フーリエ逆変換手段 ９０１１ 第４の記憶手段 ９０１２ 総和手段

Claims (10)

【特許請求の範囲】
1. 【請求項１】音響信号を時間分割する信号の時間分割手
   段と、前記信号の時間分割手段による短時間セグメント
   の信号のピッチを抽出するピッチ抽出手段と、前記ピッ
   チ抽出手段によって求められたピッチを格納する第２の
   記憶手段と、前記信号の時間分割手段の出力信号を複数
   の周波数帯域に分割する周波数帯域分割手段と、前記周
   波数帯域分割手段の出力を格納する第３の記憶手段と、
   前記第３の記憶手段の出力のエネルギーまたは実効値を
   算出するエネルギーまたは実効値算出手段と、前記エネ
   ルギーまたは実効値または実効値算出手段の出力を格納
   する第４の記憶手段と、前記第４の記憶手段と前記第２
   の記憶手段の出力から特徴量の算出を行う特徴量の算出
   手段と、ある時間区間のある周波数帯域の信号の重み値
   の算出を行う重み値の算出手段と、前記第３の記憶手段
   の出力と前記重み値の算出手段の積を算出する乗算手段
   と、前記乗算手段の出力を格納する第４の記憶手段と、
   前記第４の記憶手段の出力の総和を算出する総和手段と
   を備えたことを特徴とする雑音抑制装置。
2. 【請求項２】重み値の算出手段は特徴量の有無によって
   重み値を決定することを特徴とする請求項１記載の雑音
   抑制装置。
3. 【請求項３】特徴量の算出手段は特徴量の確からしさを
   特徴量の評価尺度とすることを特徴とする請求項１記載
   の雑音抑制装置。
4. 【請求項４】重み値の算出手段は特徴量の評価測度を乗
   算により統合することにより重み値を算出することを特
   徴とする請求項２記載の雑音抑制装置。
5. 【請求項５】重み値の算出手段は特徴量の評価測度をフ
   ァジィ積分により統合することにより重み値を算出する
   ことを特徴とする請求項２記載の雑音抑制装置。
6. 【請求項６】ピッチの抽出手段は複数のピッチ抽出方法
   を用い複数の出力の平均値をピッチとして用いることを
   特徴とする請求項２記載の雑音抑制装置。
7. 【請求項７】特徴量の算出手段は複数のピッチ抽出手段
   の出力からピッチの確からしさを算出することを特徴と
   する請求項３記載の雑音抑制装置。
8. 【請求項８】ピッチ抽出手段は異なる複数のピッチ推定
   法により推定されたピッチによりピッチを推定すること
   を特徴とする雑音抑制装置。
9. 【請求項９】特徴量の算出手段は異なる二つのピッチ推
   定法により推定されたピッチの差から推定されたピッチ
   の確からしさを求めることを特徴とする雑音抑制装置。
10. 【請求項１０】特徴量の算出手段は異なる複数のピッチ
    推定法により推定されたピッチの分散から推定されたピ
    ッチの確からしさを求めることを特徴とする雑音抑制装
    置。