(57)【要約】【課題】入力音響信号のレベルが一定ではない、あるいは強い背景雑音が含まれる場合においても調波構造区間を適切に推定できる調波構造区間推定方法を提供する。【解決手段】連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定するステップと、前記音響信号中の瞬時周波数を推定するステップと、瞬時周波数に基づき不動点を推定するステップと、占有度を計算するステップと、基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算するステップと、調波構造占有度の値に対して閾値処理を行うことで調波構造区間を推定するステップと、を備える。【選択図】図１
2004109742.tif 000002
【特許請求の範囲】
【請求項１】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定するステップと、
前記音響信号中の瞬時周波数を推定するステップと、
瞬時周波数に基づき不動点を推定するステップと、
占有度を計算するステップと、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算するステップと、
調波構造占有度の値に対して閾値処理を行うことで調波構造区間を推定するステップと、を備えたことを特徴とする調波構造区間推定方法。
【請求項２】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定するステップと、
前記音響信号中の瞬時周波数を推定するステップと、
瞬時周波数に基づき不動点を推定するステップと、
占有度を計算するステップと、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算するステップと、
各時刻の調波構造占有度について、値ごとの頻度をプロットするヒストグラムを計算するステップと、
ヒストグラム中の二つの分布の山の境界の値を計算して閾値とするステップと、を備えたことを特徴とする調波構造区間の閾値決定方法。
【請求項３】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定する基本周波数推定部と、
前記音響信号中の瞬時周波数を推定する瞬時周波数推定部と、
瞬時周波数に基づき不動点を推定する不動点推定部と、
占有度を計算する占有度計算部と、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算する調波構造占有度推定部と、
調波構造占有度の値に対して閾値処理を行うことで調波構造区間を推定する閾値処理部と、を備えたことを特徴とする調波構造区間推定装置。
【請求項４】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定する基本周波数推定部と、
前記音響信号中の瞬時周波数を推定する瞬時周波数推定部と、
瞬時周波数に基づき不動点を推定する不動点推定部と、
占有度を計算する占有度計算部と、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算する調波構造占有度推定部と、
各時刻の調波構造占有度について、値ごとの頻度をプロットするヒストグラムを計算するヒストグラム計算部と、
ヒストグラム中の二つの分布の山の境界の値を計算して閾値とする分布境界抽出部と、を備えたことを特徴とする調波構造区間の閾値決定装置。
【請求項５】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定する処理と、
前記音響信号中の瞬時周波数を推定する処理と、
瞬時周波数に基づき不動点を推定する処理と、
占有度を計算する処理と、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算する処理と、
調波構造占有度の値に対して閾値処理を行うことで調波構造区間を推定する処理と、をコンピュータに実行させるために調波構造区間推定プログラム。
【請求項６】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定する処理と、
前記音響信号中の瞬時周波数を推定する処理と、
瞬時周波数に基づき不動点を推定する処理と、
占有度を計算する処理と、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算する処理と、
各時刻の調波構造占有度について、値ごとの頻度をプロットするヒストグラムを計算する処理と、
ヒストグラム中の二つの分布の山の境界の値を計算して閾値とする処理と、をコンピュータに実行させるための調波構造区間の閾値決定プログラム。
【請求項７】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定する処理と、
前記音響信号中の瞬時周波数を推定する処理と、
瞬時周波数に基づき不動点を推定する処理と、
占有度を計算する処理と、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算する処理と、
調波構造占有度の値に対して閾値処理を行うことで調波構造区間を推定する処理と、をコンピュータに実行させるための調波構造区間推定プログラムを記録した記録媒体。
【請求項８】
連続的に入力される音響信号に対し、各時刻ごとの基本周波数を推定する処理と、
前記音響信号中の瞬時周波数を推定する処理と、
瞬時周波数に基づき不動点を推定する処理と、
占有度を計算する処理と、
基本周波数の整数倍と不動点が一致する点における占有度の和をとる調波構造占有度を計算する処理と、
各時刻の調波構造占有度について、値ごとの頻度をプロットするヒストグラムを計算する処理と、
ヒストグラム中の二つの分布の山の境界の値を計算して閾値とする処理と、をコンピュータに実行させるための調波構造区間の閾値決定プログラムを記録した記録媒体。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
この発明は、複数の音や雑音下の音響信号から、目的音の調波構造（音声の場合は有声音）が含まれている区間（有声音区間）を推定する調波構造区間推定法に関する。
調波構造区間の推定は、音声合成、音声認識、音声符号化等の信号処理の前処理として利用されている。したがって、雑音下での高精度な調波構造区間の推定は、後処理として実施される信号処理装置の性能を向上させることに寄与する。このような信号処理装置には、以下のようなものが含まれる。
１．調波構造区間の情報をもとに複数音源の混合音から各構成音を分離する音源分離装置。
２．有声音区間の情報をもとに音声を符号化する音声符号化・復号化装置。
３．騒がしい環境で人が鼻歌で歌った音の有声音区間からメロディを推定して、楽曲を検索する音楽検索装置。
４．音楽演奏の音響信号を受け取って調波構造区間を推定し、楽譜もしくは楽譜相当の音楽情報を推定するための自動採譜装置。
５．人が発した声の有声音区間における基本周波数の高さで機械にコマンドをわたす機械制御インターフェース、および機械との対話装置。
【０００２】
【従来の技術】
図５に従来の調波構造区間推定装置の構成を示す。
図５を参照して従来の調波構造区間推定装置と調波構造区間推定方法を説明する。
調波構造区間推定装置は、入力された音響信号のケプストラムピークを抽出するケプストラムピーク抽出部と、入力された音響信号のパワーを抽出する信号パワー抽出部と、単位時間当たりのゼロクロス点の数を抽出するｚｅｒｏ　ｃｒｏｓｓｉｎｇ比抽出部と、抽出された各信号と予め設定された閾値とを比較する各閾値処理部と、閾値処理部の比較結果に基づき調波構造区間を判定する統合部とから構成される。
この従来例では、調波構造区間の判定は入力音響信号の以下の特徴を利用する。
（１）ケプストラム係数のピーク値が大きくなる。
（２）パワー値が大きくなる。
（３）ｚｅｒｏ　ｃｒｏｓｓｉｎｇ比が小さくなる。
上記３つの値毎に閾値を決定し、（１），（２）の値が閾値を超えて、（３）の値が閾値以下の場合のみ、当該入力区間を調波構造区間と判定する。（例えば、非特許文献１　参照。）
【０００３】
【非特許文献１】
Ａｈｍａｄｉ，Ｓ．，ａｎｄ　Ｓｐａｎｉａｓ，Ａ．Ｓ．，”Ｃｅｐｓｔｒｕｍ−ｂａｓｅｄ　ｐｉｔｃｈ　ｄｅｔｅｃｔｉｏｎ　ｕｓｉｎｇ　ａ　ｎｅｗ　ｓｔａｔｉｓｔｉｃａｌ　Ｖ／ＵＶ　ｃｌａｓｓｉｆｉｃａｔｉｏｎ　ａｌｇｏｒｉｔｈｍ”，Ｓｐｅｅｃｈ　ａｎｄ　Ａｕｄｉｏ　Ｐｒｏｃｅｓｓｉｎｇ，　ＩＥＥＥ　Ｔｒａｎｓａｃｔｉｏｎ　ｏｎ，Ｖｏｌｕｍｅ：７　Ｉｓｓｕｅ：３　，Ｍａｙ　１９９９　Ｐａｇｅ（ｓ）：３３３−３３８．
【０００４】
【発明が解決しようとする課題】
上述した従来の調波構造区間推定装置及び方法では、入力された音響信号の状態に応じて、閾値を変更する必要があった。例えば、入力音響信号のレベルが大きくなると、音響信号のケプストラム係数やパワー値が調波構造区間であるかないかに関わらず大きくなるため、閾値を変更しない場合は、調波構造区間ではない区間を調波構造区間であると誤判定してしまう場合が多くなってしまう。また、背景雑音が大きくなる場合にも、調波構造区間以外のケプストラム係数やパワー値が上昇するため、誤判定を減少させるためには、同様に閾値を上げる必要があった。ｚｅｒｏ　ｃｒｏｓｓｉｎｇ比自体は、入力レベルとは関係のない値であるが、調波構造区間を推定するためには信頼性の低い尺度であるため、他の二つの特徴と組み合わせて用いるのが通常であり、単独では用いることが出来なかった。
また、状況に応じて適切な閾値を決定することは容易ではなく、専門家による調整が必要であった。しかも、入力音や背景雑音のレベルは、一般に、一つの録音環境でも一定ではないため、推定精度を上げるためには、その都度調整をしなおす等の必要があった。
【０００５】
【課題を解決するための手段】
上記課題を解決するために、本発明では、入力音響信号の各周波数成分が雑音の影響を受けていない度合いを表す占有度を利用し、占有度を利用した調波構造区間推定方法及び装置、調波構造区間推定の閾値決定方法及び装置を構成する。
占有度による方法では、雑音下で占有的なパワーを持つ調波成分を特定し、その特徴量のみを用いて特徴抽出を行うことで、入力音の状況によらず安定した調波構造区間推定方法、調波構造区間推定の閾値決定方法を構成する。
【０００６】
占有度は、瞬時周波数に関する次の性質に基づき定義されている。瞬時周波数φ’とは、短時間フーリエ変換の各周波数ｂｉｎを等間隔にならんだ狭帯域通過フィルタ群とみなした場合の、その各出力波の位相φの時間微分のことである。ある時刻のある帯域に強いパワーを持った占有的な周波数成分があると、短時間フーリエ変換（以下、ＳＴＦＴという）におけるその周波数近傍のｂｉｎでは、瞬時周波数がほぼ一定値になることが知られている。このため、雑音が少ない入力音響信号中の調波構造を持った音の瞬時周波数を縦軸に、ＳＴＦＴの周波数ｂｉｎを横軸にとってプロットすると、図６（ａ）の細実線で示したような階段状になる。この階段の水平部分と各周波数ｂｉｎの中心周波数ωｃが一致する点（φ’＝ωｃ、以下、不動点と呼ぶ）が、各高調波成分の周波数とみなすことができる。一方、強い雑音がある入力信号中では、瞬時周波数は明確な階段状にならず、図６（ｂ）の細実線（６００［Ｈｚ］以上）で示したように、なだらかな右上がりの線になる。
【０００７】
上述の瞬時周波数の性質を用いて、調波構造が周波数ｂｉｎの出力をどの程度占めているかを評価するために、占有度（ｄｅｇｒｅｅ　ｏｆ　ｄｏｍｉｎａｎｃｅ）Ｄ０（ωｃ）は以下で定義される。
Ｄ０（ωｃ）＝１０ｌｏｇ１０（１／Ｂ（ωｃ）２）　　　　　　　　　（１）
【数１】
2004109742.tif 000003
Ｂ（ωｃ）２は、中心周波数ωｃを持つ周波数ｂｉｎの近傍（ωｃ−Δω／２＜ω＜ωｃ＋Δω／２）のｂｉｎにおいて、φ’（ω）とωｃの差分をパワースペクトルＳ（ω）２で重み付き和をとったものである。占有的な周波数成分に対応する不動点近傍では、φ’（ω）とωｃはほぼ同じ値をとるため、Ｂ（ωｃ）２は極小値をとると期待される。その逆数（の対数）をとって、同じ点で極大値をとるようにしたものがＤ０（ωｃ）である。
占有度は、パワースペクトルＳ（ω）２で正規化した値でもあるため、非調波成分に対する調波成分の相対的な強度のみを評価するものである。信号の絶対的なパワーに依存せず、ほぼ同じ範囲の値をとる（約−４０〜０ｄＢ）。
【０００８】
図７に、幼児音声の有声音／無声音区間の占有度を各周波数ｂｉｎごとにプロットした結果を示す。図（ａ）より、有声音区間では、占有度は各調波成分に相当する規則的な鋭いピークを持ち、瞬時周波数の不動点（＝瞬時周波数と各ＦＦＴ　ｂｉｎの中心周波数が一致する点）とも一致していることがわかる。また、対数スペクトル上では背景雑音に由来する調波成分以外のピークが表れているが、占有度においてはこれらのピークは消失している。一方、図（ｂ）より無声音区間ではピークの位置は不規則かつ不明確であり不動点の位置とも一致しない。占有度のこれらの性質は、入力信号中で占有的な調波成分の構造を抽出するのに極めて有効である。
【０００９】
占有度による方法では調波構造に対応する占有度のピーク値の和をとり、その大小で有声音／無声音区間（Ｖ／ＵＶ）を判定する。各フレーム（例えば、１［ｍｓｅｃ］シフト）ごとに以下の調波構造占有度を求め、メジアンフィルタ（例えば、６１サンプルポイント）で時間方向の平滑化処理を行った後に閾値処理する。
【数２】
2004109742.tif 000004
ここで、ｌは高調波の次数、ｆ０は基本周波数Ｆ０の推定値、Ｄ０，Ｆ（ｌ・２πｆ０）はｌ次高調波の近傍の不動点における占有度（不動点がない場合はＥ（Ｄ０（ωｃ）））を返す関数である。なお、Ｅ（Ｄ０（ωｃ））は占有度のバイアスを除去する項で占有度の周波数方向の平均値を返す関数である。
【００１０】
前述したように、有声音区間では各高調波に相当する周波数と不動点および占有度のピークが一致するため、無声音区間に比べてＤｔ０（ｆ０）は大きな値をとることになる。しかも、入力音声のパワーによらず、占有的な調波成分の占有度はほぼ近い値をとるため、有声音区間においてＤｔ０（ｆ０）はある一定の大きな値の範囲におさまる。これに対し、無声音区間の占有度は小さな値で一定の範囲におさまる。この有声音／無声音区間のそれぞれの占有度がとりうる範囲の境界に閾値を設定することでＶ／ＵＶ判定を行うことが出来る。
【００１１】
【発明の実施の形態】
（調波構造区間の推定方法、装置）
図１に本発明の調波構造区間推定装置の構成例を示す。
図１を用いて、調波構造区間推定装置の構成及び調波構造区間推定方法を説明する。
調波構造区間推定装置は、基本周波数推定部と、窓切り出し手段と瞬時周波数推定手段からなる瞬時周波数推定部と、信号パワー推定手段と占有度計算手段からなる占有度計算部と、不動点推定部と、調波構造占有度推定部と閾値処理部と、から構成される。
本実施例では、まず、基本周波数推定部は、入力された音響信号から有声音／無声音区間の区別なく基本周波数を推定する。基本周波数の推定には、特願２００２−６２５１３号に示されている方法などが利用できる。その他、既存の様々な基本周波数推定法が利用可能である。
【００１２】
次に、この基本周波数推定値ｆ０を用いて有声音区間の判定を行う。
まず、窓切り出し手段は、入力音響音声を短時間単位で分析するために、窓切り出し処理を行う。後段の処理で行う瞬時周波数推定を精度よく行うためには、入力音響信号の推定基本周波数ｆ０に応じて窓の切り出し幅を変えることが有効であることが知られている（詳しくは文献［１］阿竹　他，“瞬時周波数に基づく雑音環境下でのピッチ推定”，電子情報通信学会論文誌，Ｖｏｌ．Ｊ７９−Ｄ−ＩＩ，Ｎｏ．１１，ｐｐ．１７７１−１７８１，１９９６参照）。例えば、ある時刻の推定基本周波数をｆ０［Ｈｚ］としたとき、３．５／ｆ０［ｓｅｃ］程度の長さのハニング窓を利用すればよい。ハニング窓以外にもハミング窓やブラックマン窓など、様々な既存の窓関数を利用することが出来る。さらに、瞬時周波数推定精度は下がるが可変長時間窓ではなく固定長時間窓（例えば、４２［ｍｓｅｃ］ハニング窓など）を用いても、同様に以降の分析を進めることが出来る。
【００１３】
次に、占有度を求めるのに必要な量である、信号のパワースペクトルＳ（ω）２と瞬時周波数φ’（ωｃ）を求める。
信号パワー推定手段において、音響信号のパワースペクトルは、まず、切り出した各短時間入力音響信号に対して、短時間フーリエ変換を施して周波数領域の信号に変換すると、各周波数帯域（以後、中心周波数をωｃと表す）ごとの特徴を表す複素信号が得られる。この複素信号の平方を計算することで、各周波数ωｃごとの信号パワーＳ（ωｃ）２を求めることが出来る。周波数領域への変換は、ウォーブレット変換、余弦変換など他の手法を用いてもよい。
瞬時周波数推定手段において、瞬時周波数は、上記のようにして求められる周波数領域の複素信号に対して、各周波数ごとに位相を時間微分することで求めることが出来る。このためには、例えば、例えば、可変長窓切り出し処理で切り出すディジタル信号波形を切り出す位置を１サンプルずらして求められる二つの信号波形から得られる周波数領域の複素信号の位相φ（ｔ１），φ（ｔ２）に対して、各位相の差を１サンプル間の時間差△ｔ＝ｔ２−ｔ１で割ってやることで、近似的に微分を行うことが出来る。
【数３】
2004109742.tif 000005
もしくは、文献［２］（阿部　他“調波成分の瞬時周波数を用いた基本周波数推定方法”，電子情報通信学会論文誌，Ｖｏｌ．Ｊ８３−Ｄ−ＩＩ，Ｎｏ．１１，ｐｐ．２０７７−２０８６，２０００）に示されている方法などを用いて、ひとつの時間窓で切り出したディジタル信号波形から直接計算する方法も知られている。
【００１４】
占有度計算手段における占有度の推定は、上記のようにして周波数変換の各中心周波数ωｃごとに得られる信号パワーＳ（ωｃ）２、瞬時周波数φ’（ωｃ）に基づき、式（１），（２）を用いて行う。まず、各周波数帯の中心周波数ごとに、中心周波数ωｃと、その近傍の周波数の瞬時周波数φ’（ω）との２乗誤差（φ’（ω）−ωｃ）２を計算し、それに信号パワーＳ（ω）２を乗じた値を、その近傍の周波数全体で総和をとるとともに、同じ近傍周波数で信号パワーＳ（ω）２のみの総和をとった値で割ることで、各中心周波数ωｃに対する占有度を推定できる。総和をとる近傍の周波数の範囲△ωは、基本周波数ｆ０［Ｈｚ］を用いてより適切に決定することが出来、例えばｆ０［Ｈｚ］に０．９程度の値を乗じて得られる範囲などを用いればよい。
【００１５】
一方、不動点推定部において、不動点は各中心周波数ωｃとその点での瞬時周波数φ’（ωｃ）から求めることが出来る。二つのとなりあった中心周波数ωｃ１，ωｃ２およびその瞬時周波数φ’（ωｃ１），φ’（ωｃ２）の間に式（６）の不等式が成立する場合、ωｃ１，ωｃ２の間に不動点が存在し、不動点の周波数φ’は、式（５）で与えられる。
【数４】
2004109742.tif 000006
Ｗｈｅｒｅ　　φ’１＞ωｃ１　　ａｎｄ　　φ’２＜ωｃ２　　　　（６）
ただし、φ’はφ’（ωｃ１），φ’（ωｃ２）の間の値をとるため、以後のディジタル信号処理では、φ’（ωｃ１），φ’（ωｃ２）のうちφ’に近い方を不動点として扱うようにする。
【００１６】
次に、調波構造占有度推定部は、調波構造占有度の推定を式（３）に基づいて行う。さらに、各時刻の調波構造占有度に対し、時間方向の平滑化処理としてメジアンフィルターを用いることは有効である。メジアンフィルターは、各時刻を中心にして、前後一定数のサンプルを切り出すとともに、切り出したサンプルの値をその時刻の値とするフィルターである。すべての時刻に対して同じ処理を行うことで、平滑化を行える。前後一定数のサンプルの長さとしては、例えば３０［ｍｓｅｃ］（前後合わせて６０［ｍｓｅｃ］）程度の長さが有効である。
最後に、閾値処理部は、時間平滑化処理を行った調波構造占有度に対し、その値が予め決められた閾値より大きな値をとる時刻を調波構造区間とすることで、調波構造区間を推定することが出来る。
【００１７】
（調波構造区間の閾値推定法）
図２に調波構造区間の閾値決定装置の構成例を示す。
多くの音響信号から構成されるデータベースを用いて、調波構造区間を推定するための閾値の決定方法及び装置を説明する。
まず、データベースに含まれる各音響信号に対し、前節で示される方法に基づき、各時間窓ごとに調波構造占有度を求める。すなわち、図１に示された基本周波数推定部、窓切り出し手段と瞬時周波数推定手段からなる瞬時周波数推定部と、信号パワー推定手段と占有度計算手段からなる占有度計算部と、不動点推定部と、調波構造占有度推定部を用いて調波構造占有度を求める。
次にヒストグラム計算部では、占有度−１０〜１８０［ｄＢ］を細かいいくつかの区間に分け、得られた各調波構造占有度がどの区間に含まれるかを判定することで、各区間に含まれる値を調波構造占有度がとる総数を計算する。このようにして得られた各区間ごとの回数を、横軸に調波構造占有度の値をとってプロットすると、左右二つの分布の山をもつヒストグラムを生成する（図３　参照）。
分布境界抽出部では、この二つの山の境界の低くなっているところ（例えば、最小値を検出して）を閾値とすることで、調波構造区間推定の閾値を得ることができる。
なお、調波構造占有度の値は、入力音響信号の状態にあまり依存しないため、このようにして求められた閾値は、異なる環境で録音された音響信号に対しても利用することができる。
【００１８】
本発明の音響信号の調波構造区間推定装置及び調波構造区間の閾値決定装置はＣＰＵやメモリ等を有するコンピュータと、ユーザが利用する端末と、記録媒体とから構成される。記録媒体は、ＣＤ−ＲＯＭ、磁気ディスク装置、半導体メモリ等の機械読み取り可能な記録媒体であり、ここに記録された調波構造区間推定プログラム及び調波構造区間の閾値決定プログラムはコンピュータに読み取られ、コンピュータの動作を制御し、コンピュータ上に前述した各構成要素を実現する。
【００１９】
【発明の効果】
図３（ａ）に、幼児音声データからランダムに選んだ１７４９データから、また図３（ｂ）に雑音を含まない成人の音声（２８人×３０発話）と女性の音声から抽出した調波構造占有度のヒストグラムを示す。各図では、横軸上３４［ｄＢ］付近をはさんで左右に一つづつ分布の山が出来ている。左が無声音区間、右が有声音区間に相当する分布である。成人でも幼児でも同様な性質を持つ分布が得られ、しかも有声音／無声音の各分布の境界がほぼ同じ程度の値になる。これはＶ／ＵＶ判定の尺度として有効であるといえる。以上の考察により、音声の有声音／無声音の判定の閾値は、３４［ｄＢ］付近に設定すればよいことがわかる。
最後に発明の構成例を用いて、幼児音声の分析を行った結果の例を図４に示す。図の実線で示すように、ほぼ正しく基本周波数ｆ０および有声音／無声音区間（Ｖ／ＵＶ）を推定出来ていることがわかる。
【図面の簡単な説明】
【図１】本発明の調波構造区間推定装置の構成例を示す図。
【図２】本発明の調波構造区間の閾値決定装置の構成例を示す図。
【図３】（ａ）雑音下の幼児音声データの調波構造占有度のヒストグラムと（ｂ）非雑音下の成人男性と女性音声データの調波構造占有度のヒストグラム。
【図４】幼児音声のスペクトルグラムと基本周波数推定及び有声音／無声音区間の判定結果を示す図。
【図５】従来の調波構造区間推定装置の構成を示す図。
【図６】（ａ）非雑音下における有声音の瞬時周波数と対数パワースペクトルと占有度を示す図、及び（ｂ）白色雑音下における有声音の瞬時周波数と対数パワースペクトルと占有度を示す図。
【図７】（ａ）幼児音声有声音区間の占有度と対数パワースペクトルと不動点を示す図、及び（ｂ）幼児音声無声音区間の占有度と対数パワースペクトルと不動点を示す図。
【図１】
2004109742.tif 000007
【図２】
2004109742.tif 000008
【図３】
2004109742.tif 000009
【図４】
2004109742.tif 000010
【図５】
2004109742.tif 000011
【図６】
2004109742.tif 000012
【図７】
2004109742.tif 000013