JPS62180400A - 極零分析装置 - Google Patents

極零分析装置

Info

Publication number
JPS62180400A
JPS62180400A JP61023141A JP2314186A JPS62180400A JP S62180400 A JPS62180400 A JP S62180400A JP 61023141 A JP61023141 A JP 61023141A JP 2314186 A JP2314186 A JP 2314186A JP S62180400 A JPS62180400 A JP S62180400A
Authority
JP
Japan
Prior art keywords
value
zero
autocorrelation
pole
parameter
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP61023141A
Other languages
English (en)
Other versions
JPH0644198B2 (ja
Inventor
幸夫 三留
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
NEC Corp
Original Assignee
NEC Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by NEC Corp filed Critical NEC Corp
Priority to JP61023141A priority Critical patent/JPH0644198B2/ja
Priority to CA000528806A priority patent/CA1257002A/en
Priority to US07/010,597 priority patent/US4845754A/en
Publication of JPS62180400A publication Critical patent/JPS62180400A/ja
Publication of JPH0644198B2 publication Critical patent/JPH0644198B2/ja
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Landscapes

  • Filters That Use Time-Delay Elements (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】 〔産業上の利用分野〕 本発明は極零分析装置、特に音声等の信号のスペクトル
を近似する極零パラメータの値を抽出する極零分析装置
に関するものである。
r従来の技術〕 音声&成や音声認識の分野では、音声のスペクトルを近
似するパラメータの値を抽出することが重要である。ま
た、一般の信号についても、そのスペクトルを近似する
パラメータの値を抽出することが必要になる場合がある
音声等の信号(以下、単に音声等と表す)のスペクトル
を近似するパラメータとしては、極あるいは極零がしば
しば用いられる。これは、極や零はその物理的意味が明
確である等の特徴があり、応用上有利であるためである
従来、音声等のスペクトルを近似する極パラメータの値
を抽出するものとして、音声等から抽出された線形予測
係数を係数とする高次代数方程式を例えばニュー1〜ン
・ラプソン法等の逐次近似法を用いて解く方法が知られ
ている。この第1の従来例は、例えばジエイ・ディー・
マーケルとエイ・エイチ・グレイ(J、D、Marke
l、^、11.Gray)の著書[音声の線形予測(鈴
木訳) J  (Linear Prediction
 of 5peprb)7第7章に示されている。
極パラメータの値を抽出する第2の例として。
自己相関領域において逆フィルタリングするものが、伏
木田による論文「自己相関領域で逆フィルタリングを用
いたホルマンI・の多段推定方式」(日本音響学会音声
研究会資料581−41>に示されている。この第2の
従来例では、第1(段の2次の極回路の係数をαに、I
とα8.2とし、その入力の自己相関値をr−に−1,
1としたとき、出力の自己相関値r k、 Iを次式に
よって求める。
rk   i   ”   (If   α 2.  
鳳f  α 2に、  2)   rh−t、+−(α
に、1−αに、1  °αh、 2)rk−1,141
−(I  h、z(rb−t、l−2↑rk−1,1+
2>    −−−(1)第1段の入力には音声等の自
己相関を用いて。
その出力を次段の入力とし、式(1)によって出力の自
己相関値を逐次求め、最終段(これを第に段とする)の
出力の自己相関値rk、0を求める。これが自己相関領
域における逆フィルタリングと呼ばれる方法である。な
お、rk、0は予測残差のパワーあるいはエラーパワー
と呼ばれ、これを算出するには入力音声の自己相関値r
 o、 tは0次から2に次のものが必要である。
あらかじめ、極パラメータ値のテーブルを用意しておき
、それから各段の係数α3.1 ・α1.2を求めて上
記のrk、0を算出する。種々の極パラメータ値の候補
に対してr14.。の最小値を与えたものを最適な近似
を与える極パラメータの値として出力する。
一方、極零パラメータの値を抽出するものの一例(第3
の例)として、シー・ティー・マリスとアール・エイ・
ロバーツ(C,T、Mullis、R,A、Rober
Ls)によるアイイーイーイートラン′ザクジョン(1
1EE Transaction )エイニスニスと−
24(^S S +’−24>、N代3のページ226
から238に掲載の論文、[ザ ユース オブ セカン
ドオーダ イン・フォメーショ〉 イン ザ アブ口り
シ、メイショ〉 オブ ディスク−1−トタイム リニ
アー ジステムズ(The l1se of 5eco
nd−Or+ler 1nrorn+ati。
n  in  tl+e  Approximatio
n  of  0iscre!、e−Tin+e  l
、1near Systems)」に示されたものがあ
る。
これは、音声のスペクトル包絡H(ej6J)を、伝達
関数が次式で表される極零システムによって近似するも
のである。
ただしaO=1 ここにnは分母多項式の次数即ち極回路の次数で、1n
は分子多項式の次数即ち零回路の次数である。
この方法は次式の値を最小にするalおよびqkが最適
な近似を与えるとうい原理に基づいている。
ここに、rh−11は音声等の入力信号の自己相関値、
11−はそれに関連したインパルス応答である。
このとき、最適なatは次の連立方程式を解いて得られ
る。
ここに、 である。
また、最適なQbは、式(4)の解のa、を用いて、次
式によって与えられる。
mよりも高い次数の線形予測法、即ち第1の従来例の方
法によって容易に求めることができる。そのときには、
高次代数方程式の根を求める必要はなく、単にhlさえ
求めればよい。
本従来例では、伝達関数の分母多項式と分子多項式の係
数atとhが得られる。これから極零パラメータの値を
得るためには第1の従来例と同様に高次(0次もしくは
m次)代数方程式の根を求める必要がある6 〔発明が解決しようとする問題点〕 前記従来例のうち、第1の例即ち線形予測係数を係数と
する高次代数方程式を解くものは、高次代数方程式を解
く際に多くの演算を必要とすることや、極の周波数やバ
ンド幅を安定に求めることが困難であるといった問題点
があった。
また第2の例では、あらかじめ音声等のスペクトルを表
す極パラメータとして妥当と思われる候補の中から最適
な値を抽出するので安定に求められるが、多くの演算を
必要とするといった間1題点があった。
更にこれら2つの従来例では極パラメータの値しか得ら
れないので、音声等のスペクトルの形状によっては比較
的低い次数では最適な近似を与えるパラメータ値とはな
らず、逆にそのような場合には高い次数のパラメータを
求めなければならず、いっそう多くの演算が必要になる
といった問題点があった。
第3の従来例、即ちマリスとロバーツによるものは、極
パラメータと零パラメータの値が得ら五るので、比較的
低い次数で最適な近似を与えるパラメータ値となるが、
第1の例と同様に高次代数方程式を解かねばならず第1
の例と同様の問題点があった。この点は、音声の場合に
は零の数は少ないので問題とならないが、極の次数は少
なくとも10次程度必要であり問題となる。
本発明の目的は、音声等の信号のスペクトルを近似する
極零パラメータの値を比較的少ない演算量で抽出できる
極零分析装置を提供することにある。
L問題点を解決するための手段〕 本発明の極零分析装置は、入力信号から時間窓内の自己
相関値とインパルス応答を算出し一時記憶rる手段と、
極パラメータの値を記憶し入力コードに対する極パラメ
ータの値を出力する極パラメータ値テーブルと、その極
パラメータ値から自己相関領域における逆フィルタリ〉
′グを行なうための係数を算出する手段と、その係数を
用いて前記一時記憶された自己相関値を自己相関領域に
おいて逆フィルタリングする手段と、極パラメータ値か
ら時間領域における逆フィルタリングを行なうための係
数を算出する手段と、その係数を用いて前記一時記憶さ
れたインパルス応答を時間領域において逆フィルタリン
グする手段と、前記自己相関領域において逆フィルタリ
ングする手段の出力と前記時間領域において逆フィルタ
リングする手段の出力とから4次以下の極パラメータ値
とエラーパワーを算出する手段と、そのエラーパワーの
最小値を検出する手段と、前記時間領域において逆フィ
ルタリングする手段の出力と4次以下の極パラメータ値
とから零パラメータ値を算出する手段を有し、そのエラ
ーパワーの最小値を与えた時の前記極パラメータ値テー
ブルの出力と前記4次以下の極パラメータ値と前記零パ
ラメータ値を算出する手段の出力とを抽出された横倒パ
ラメータの値として出力するように制御する手段から構
成される。
[作用〕 本発明では、スペクI・ル含近似する極零システムを4
次以下(これをIく次とする)の部分と残り(L次とす
る)の部分に分けに次の部分に第3の従来例の方法を利
用するものである。4次以下の場合には、有限回の演算
で多項式の根を求められることが知られており、第3の
従来例のような問題はない。L次の部分に相当するのは
自己相関領域および時間領域における逆フィルタリング
を行なう部分である。零を求めるために時間領域におけ
る逆フィルタリングを行なう点が第2の従来例と異なる
極パラメータ値テーブルから候補を選び自己相関領域お
よび時間領域での逆フィルタリングを行なう手段の係数
を算出しそれにより逆フィルタリングを行なう。その出
力に対し第3の従来例の方法で■く次の極零パラメータ
の最適値とエラーパワーを算出する。このときの最適値
はに次の部分のみに対する最適値であり、全体の最適値
とは限らないので、極パラメータのさまざまな候補に対
してエラーパワーを算出してその最小値を検出し、その
ときの極パラメータの候補の値とKの部分から得られる
極パラメータを最適な極パラメータの値とする。零パラ
メータの値は1時間領域で逆フィルタリングしたしのを
さらにに次の部分から得られた極パラ、メータを用いて
さらに逆フィルタリングすれば、式(6)と等価なもの
が求まる7本発明において必要な時間領域での逆フィル
タリ〉′グの結果のサンプル数はmであり、このために
必要な逆フィルタリングの演算量はわずかでよい、。
〔実施例] 次に、本発明について図面を用いて説明する。
第1図は本発明の第1の実施例を示すブロック図である
あらかじめ定められた分析フレームの時間ごとに制御回
路101から信号線123,124,125.126,
127,128,129,130および131を介して
それぞれ自己相関インパルス応答算出回路102.バッ
ファ103.バッファ104.自己相関逆フィルタ10
7.波形逆フィルタ108.M寄算出回路109.バ・
ソファ110および比較回路111に吋して初期化を指
示する情報が送られて、初期化が行なわれる。
自己相関インパルス応答算出回路102では制御回路1
01から送られる指示に従い、信号線121から入力さ
れる信号の自己相関係数値とインパルス応答の始めの部
分が算出されそれぞれ信号線135および136を介し
てバッファ103およびバ・ソファ104へ送られる。
このとき必要なインパルス応答のサンプル数は、分析し
ようとする零の次駿に等しい。
バッファ103では、制御回路101からの初期化の指
示に従って内部がクリアされた後、自己相関インパルス
応答算出回路102から送られる自己相関係数値が蓄え
られ、それが信号線137を介して自己相関逆フィルタ
107に送られる。
ここに蓄えられた自己相関係数値は、以後、制御回路1
01から送出のt旨示があるごとにくり返し自己相関逆
フィルタ107に送られる。
一方バッファ104でも同様に、制御回路101からの
初期1ヒの指示に従って内部がクリアされた後、自己相
関インパルス応答算出回路102から送られるインパル
ス応答値が蓄えられ、信号線138を介して波形逆フィ
ルタ108にくり返し送られる。
制御回路101からは極値テーブル105に対して信号
線126を介して極パラメータ値の候補のコードが送ら
れる。それに対する極パラメータ値が極値テーブル10
5から信号線139を介して係数算出回路106に送ら
れ、そこで自己相関領域での逆フィルタリングのための
係数、即ち前記式(1)のかっこ内の係数が算出され信
号線140を介して自己相関逆フィルタ107に送られ
、同時に時間領域での逆フィルタリングのための係数、
即ち式(1)のα工、1とαに、jが算出され信号線1
41を介して波形逆フィルタlO8に送られる。この過
程は制御回路101から極パラメータの候補が送られる
たびにくり返される。
自己相関逆フィルタ107では、制御回路101からの
指示に従って、係数算出回路106から送られる係数値
を用いてバッファ103か15送られる自己相関係数値
が自己相関領域において逆フィルタリングされる。即ち
式(1)を逐次用いることによって出力の4次までの自
己相関係数値が算出される。この値は信号線142を介
して極零算出回路109に送られる。この処理も制御回
路101からの指示があるごとにくり返される。
同様に、波形逆フィルタ108では、制御回路101か
らの指示に従って、係数算出回路106から送られる係
数値を用いてバッファ104から送られるインパルス応
答の波形が時間領域において逆フィルタリングされ、信
号線143を介して極零算出回路109に送られる。こ
の処理も制御回路101からの指示があるごとにくり返
し行なわれる。
極零算出回路109では、自己相関逆フィルタ107か
ら送られた4次までの自己相関係数値と波形逆フィルタ
108から送られた逆フィルタリングされたインパルス
応答値に対する最適な4次の極零パラメータ値とそのと
きのエラーパワーが直達の第3の従来例の方法によって
算出され、エラーパワーは信号数145を介して比較回
路111に送られる。このあと、制御回路101からの
指示があれば、極零パラメータ値は信号線1 ll 4
を介してバッファ11(〕に送られる。
比較回路111では、制御回路101から初期(ヒの指
示が送られたときには、エラーパワーの最小値を検出す
る際の初期値が信号線133を介してメモリ112に送
られる。次に、比較回路111では制御回路101から
の指示に従って、極零算出回路109から送られるエラ
ーパワーの値とメモリ112から信号線134を介して
読み出された値とが比較され1.比較結果が信号線13
2を介して制御回路101に送られるとともに、極零算
出回路109から送られたエラーパワーの値の方が小さ
いときは、極零算出回路109から送られたエラーパワ
ーの値がメモリ112に送られてその内容が書き換えら
れる。この処理も制御回路101からの指示があるごと
にくり返される。
制御回路101はS比較回路111から送られた比較結
果が極零算出回路109から送ら・れなエラーパワーの
値の方が小さいことを示している場合は、極値テーブル
105に対してそのときの極パラメータ値を信号線14
6を介してバッファ110に送らせ、極零算出回路10
9に対してそのときの4次の極零パラメータ値を信号線
144を介してバッファ110に送らせる。バッファ1
10では制御回路101からの指示があるごとに、極値
テーブル105から送られる極パラメータ値と極零算出
回路109から送られる極零パラメータの値に書き換え
られ記憶される。
このようにして、極パラメータの種々の候補について逆
フィルタリングが行なわれた時点では、バッファ110
には最適な極零パラメータの値が記憶されている。各分
析フレームの処理の最後に制御回路101からバッファ
110に出力の指示が送られ、最適な極零パラメータの
値が信号線122から出力される。
以上の処理が分析フレームごとにくり返される。
なお、本実施例では極零算出回路109において抽出さ
れる極零パラメータの次数は4次としたが、4次以下の
ものでも良い。
第2図は本発明の第2の実施例を示すブロック図である
本実施例では、音声合成のための極零分析などにおいて
、音声合成モデルに合ったものを抽出できるように、極
パラメータ値に制限を加えられるようにしたものである
。そのために、第1の実施例に対し比較回路213とメ
モリ214が追加されている。
本実施例における自己相関インパルス応答算出回路20
2.バッファ203.バッファ204゜極値テーブル2
05.係数算出回路206.自己相関逆フィルタ207
.波形逆フィルタ208゜極零算出回路209.バッフ
ァ210.比較回路211およびメモリ212は、それ
ぞれ第1の実施例における自己相関インパルス応答算出
回路102、バッファ103.バッファ104.極値テ
ーブル105.係数算出回路106.自己相関逆フィル
タ107.波形逆フィルタ108.極零算出回路109
.バッファ110.比較回路111およびメモリ112
に対応し同様の動作をする。
制御回路201の行なう制御で第1の実施例と異なる点
について以下に説明する。
まず初期化の際には、極パラメータの値に対する制限の
値を信号線247を止して比較回路213に送り、比較
回路213はさらにその値を信号線250を介してメモ
リ214に送る。
極値テーブル205に対してモデルの制限に合った極パ
ラメータの候補のコードを送るごとに、極零算出回路2
09に対しても、そのつと4次の極零パラメータの値を
比較回路213に送らせる。
比較回路213では、その値がメモリ214から信号線
249を介して読み出された制限の値に適合するか否か
を判定し、その結果を信号線248を庁して制御回路2
01に送る。
制御回路201では比較回路213からの判定結果が、
制限に適合していることを示しているときだけ比較回路
211からのエラーパワーの比較結果を参照し、最適な
極パラメータ値の抽出を行なう、。
本実施例では、制限に合った極パラメータの値のみを抽
出できるので、音声のボルマン1−分析なとではより安
定したホルマント抽出が行なえるという効果がある。
〔発明の効果1 以上説明したように本発明は、第2の従来例より少ない
演算量で、より精度良くスペクトルを近似できる術零パ
ラメータ値を、高次の代数方程式を解く必要がないため
第3の従来例より安定に求められるという効果がある。
【図面の簡単な説明】
第1図および第2図はそれぞれ本発明の第1および第2
の実施例を示すブロック図である。 101.201・・・制御回路、102,202・・・
自己相関インパルス応答算出回路、103,104.1
10,203,204.210・・・バッファ、105
.205・・・極値テーブル、106,206・・・係
数算出回路、107.207・・・自己相関逆フィルタ
、108,208・・・波形逆フィルタ、109.20
9・・・極零算出回路、111,211.213・・・
比較回路、112,212.214・・・メモリ。 1−1゜

Claims (1)

    【特許請求の範囲】
  1. 音声等の信号から抽出された自己相関値とインパルス応
    答に基づいて極零パラメータの値を抽出する型の極零分
    析装置において、入力信号から時間窓内の自己相関値と
    インパルス応答を算出し一時記憶する手段と、極パラメ
    ータの値を記憶し入力コードに対する極パラメータの値
    を出力する極パラメータ値テーブルと、その極パラメー
    タ値から自己相関領域における逆フィルタリングを行な
    うための係数を算出する手段と、その係数を用いて前記
    一時記憶された自己相関値を自己相関領域において逆フ
    ィルタリングする手段と、極パラメータ値から時間領域
    における逆フィルタリングを行なうための係数を算出す
    る手段と、その係数を用いて前記一時記憶されたインパ
    ルス応答を時間領域において逆フィルタリングする手段
    と、前記自己相関領域において逆フィルタリングする手
    段の出力と前記時間領域において逆フィルタリングする
    手段の出力とから4次以下の極パラメータ値とエラーパ
    ワーを算出する手段と、そのエラーパワーの最小値を検
    出する手段と、前記時間領域において逆フィルタリング
    する手段の出力と4次以下の極パラメータ値とから零パ
    ラメータ値を算出する手段を有し、そのエラーパワーの
    最小値を与えた時の前記極パラメータ値テーブルの出力
    と前記4次以下の極パラメータ値と前記零パラメータ値
    を算出する手段の出力とを抽出された極零パラメータの
    値として出力するように制御する手段を備えることを特
    徴とする極零分析装置。
JP61023141A 1986-02-04 1986-02-04 極零分析装置 Expired - Lifetime JPH0644198B2 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP61023141A JPH0644198B2 (ja) 1986-02-04 1986-02-04 極零分析装置
CA000528806A CA1257002A (en) 1986-02-04 1987-02-03 Pole-zero analyzer
US07/010,597 US4845754A (en) 1986-02-04 1987-02-03 Pole-zero analyzer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP61023141A JPH0644198B2 (ja) 1986-02-04 1986-02-04 極零分析装置

Publications (2)

Publication Number Publication Date
JPS62180400A true JPS62180400A (ja) 1987-08-07
JPH0644198B2 JPH0644198B2 (ja) 1994-06-08

Family

ID=12102281

Family Applications (1)

Application Number Title Priority Date Filing Date
JP61023141A Expired - Lifetime JPH0644198B2 (ja) 1986-02-04 1986-02-04 極零分析装置

Country Status (1)

Country Link
JP (1) JPH0644198B2 (ja)

Also Published As

Publication number Publication date
JPH0644198B2 (ja) 1994-06-08

Similar Documents

Publication Publication Date Title
US6047254A (en) System and method for determining a first formant analysis filter and prefiltering a speech signal for improved pitch estimation
JPS63113600A (ja) 音声信号の符号化及び復号化のための方法及び装置
JP2778567B2 (ja) 信号符号化装置及び方法
CN109473091A (zh) 一种语音样本生成方法及装置
Bozkurt et al. Mixed-phase speech modeling and formant estimation, using differential phase spectrums
JPS62180400A (ja) 極零分析装置
US4845754A (en) Pole-zero analyzer
US6049814A (en) Spectrum feature parameter extracting system based on frequency weight estimation function
US20060229868A1 (en) Method for estimating resonance frequencies
Tremain Linear predictive coding systems
JPH0644197B2 (ja) 極零分析装置
JPH0644196B2 (ja) 極パラメ−タ値抽出装置
JPH0235990B2 (ja)
JP2539351B2 (ja) 音声合成方法
US4899386A (en) Device for deciding pole-zero parameters approximating spectrum of an input signal
JPH08328593A (ja) スペクトル分析方法
JPH0670750B2 (ja) 位相特性の推定を伴う極零分析装置
JPH0670751B2 (ja) 位相特性の推定を伴う極零分析装置
JPS63226699A (ja) 極零分析装置
JPS6367200B2 (ja)
Salor et al. Kalman filter approach for pitch determination of speech signals
JPH0117599B2 (ja)
JPH0863192A (ja) スペクトルパラメータ抽出装置
JPH0670754B2 (ja) 極零分析装置
JPH0237402A (ja) パラメータ推定方式