JP3784361B2 - Auditory filter shape estimation method and apparatus - Google Patents

Auditory filter shape estimation method and apparatus Download PDF

Info

Publication number
JP3784361B2
JP3784361B2 JP2002286594A JP2002286594A JP3784361B2 JP 3784361 B2 JP3784361 B2 JP 3784361B2 JP 2002286594 A JP2002286594 A JP 2002286594A JP 2002286594 A JP2002286594 A JP 2002286594A JP 3784361 B2 JP3784361 B2 JP 3784361B2
Authority
JP
Japan
Prior art keywords
subject
shape
notch width
sound
auditory filter
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.)
Expired - Fee Related
Application number
JP2002286594A
Other languages
Japanese (ja)
Other versions
JP2004121338A (en
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.)
Rion Co Ltd
Original Assignee
Rion Co Ltd
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 Rion Co Ltd filed Critical Rion Co Ltd
Priority to JP2002286594A priority Critical patent/JP3784361B2/en
Priority to US10/345,041 priority patent/US7048692B2/en
Priority to DE10301677A priority patent/DE10301677B4/en
Publication of JP2004121338A publication Critical patent/JP2004121338A/en
Application granted granted Critical
Publication of JP3784361B2 publication Critical patent/JP3784361B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)

Description

【0001】
【発明の属する技術分野】
本発明は、人の聴覚フィルタの形状を推定する方法とその装置に関する。
【0002】
【従来の技術】
現在、最も頻繁に行われている難聴者の聴覚特性検査は、聴力検査(オージオグラムの測定)と語音明瞭度検査である。聴力検査では、難聴者の最小可聴閾値の周波数特性を知ることができ、語音明瞭度検査では、難聴者の言葉の聞き取り能力を知ることができる。
しかし、難聴者の聴覚特性は個人毎に様々であるので、これら2つの方法のみでは、その複雑な聴覚特性の一端しか把握できないと考えられている。
【0003】
一般に、難聴者は、聴力の低下及び言葉の聞き取り能力の低下に加えて、周波数分解能が低下していると言われている。ここで、周波数分解能とは、周波数が異なる2つの音を聞き分ける能力である。健聴者は、例えば、1kHzと1.2kHzといった周波数が近接した2つの音を聞き分けることができるが、周波数分解能が低下した難聴者は、これら2つの音を聞き分けることができない。
【0004】
この周波数分解能の低下の度合いが大きくなると、言葉の聞き取り能力が低下したり、雑音下での音声弁別能力が低下したりする。この周波数分解能の低下の度合いを知ることは、難聴の診断、難聴者の聴覚特性の把握、補聴器フィッティング等においては非常に有益である。
【0005】
また、近年、人間の聴覚の周波数分析のメカニズムを表現するためのモデルとして、聴覚フィルタが提案されている。これは、人間の内耳の周波数分析のメカニズムを複数の帯域フィルタバンクで表現する考え方である。このフィルタバンク内の個々のフィルタ(聴覚フィルタ)の形状は、通常、ノッチノイズマスキングを用いて測定される。これまでに、この聴覚フィルタの理論を利用した難聴者のための周波数分解能簡易測定法が、特開平6−327654号公報に開示されている。人間の聴覚フィルタの形状は、roex(p,r)フィルタによってモデル化できることが知られている。
【0006】
【発明が解決しようとする課題】
しかし、ノッチノイズマスキングによる聴覚フィルタの測定は、個々の被験者の聴覚フィルタ形状を高精度で測定することができると言われているが、測定に要する時間は極めて長く、実際に耳鼻科臨床の現場や補聴器フィッティングの現場で、個々の難聴者について随時測定することは事実上不可能であった。特開平6−327654号公報では、短時間で周波数分解能が劣化しているかどうかを測定する方法が提案されているが、この方法は劣化の有無を判定することは可能であるが、周波数分解能の劣化度合い、延いては聴覚フィルタの形状そのものを測定することはできなかった。
【0007】
更に、聴覚フィルタの形状は、入力信号のレベルに応じて変化するといわれている。健聴者の聴覚フィルタ測定には、過去の様々な知見から、40dB/Hz程度のレベルのマスカー(ノッチノイズ)が最適であると言われており、40dB/Hz以上のレベルで測定した聴覚フィルタはその形状がブロードになるといわれている。しかし、全ての難聴者においてもこれらのレベルに応じて、健聴者と同様の聴覚フィルタの形状変化が起こっているとは限らない。個々の難聴者の、レベルに応じた聴覚フィルタの形状変化特性を測定する方法はこれまでのところ皆無である。
【0008】
本発明は、従来の技術が有するこのような問題点に鑑みてなされたものであり、その目的とするところは、難聴の診断、難聴者の聴覚特性の把握、補聴器フィッティング等を効率的かつ正確に行うために有効な聴覚フィルタの形状を短時間で正確に推定する方法とその装置を提案しようとするものである。
【0009】
【課題を解決するための手段】
上記課題を解決すべく請求項1に係る発明は、聴覚フィルタの形状をモデル化したroex(p,r)フィルタの係数pを求めて聴覚フィルタの形状を推定する方法であって、周波数fにおける被験者の最小可聴閾値Tに任意の値xを加算した信号音Sから任意の値aを減じた信号音S’を生成し、次いで被験者の外耳と中耳の周波数特性を考慮した前記周波数fをノッチに含む、ノッチ幅g、レベルNX’のマスカーを生成し、次いで前記信号音S’に前記マスカーを重畳した検査音を被験者に提示し、ノッチ幅gを変化させながら、その被験者の最小ノッチ幅gX-aを測定し、次いで最小ノッチ幅gX-aから被験者に適したノッチ幅gの上限値gmaxと前記最小ノッチ幅gX-aと係数rに相当する前記値xと前記値aとから係数pを算出し、次いで算出した係数pと係数rに相当する前記値xとから聴覚フィルタの形状を推定するものである。
【0010】
請求項2に係る発明は、請求項1に記載の聴覚フィルタの形状推定方法において、前記値xをパラメータとして聴覚フィルタの形状を推定するものである。
【0011】
請求項3に係る発明は、請求項1又は2に記載の聴覚フィルタの形状推定方法において、周波数fにおける被験者の最小可聴閾値Tに任意の値xを加算した信号音Sに基づいて限界マスキングレベルNX’を決定するものである。
【0012】
請求項4に係る発明は、請求項1、2又は3に記載の聴覚フィルタの形状推定方法において、最小ノッチ幅gX-aから被験者に適したノッチ幅gの上限値gmaxを推定するものである。
【0013】
請求項5に係る発明は、請求項1、2、3又は4に記載の聴覚フィルタの形状推定方法において、限界マスキングレベルNX’及びまたは最小ノッチ幅gX-aを測定する際に、被験者に対するマスカーの提示を開始した後に、所定の時間間隔を置いてから信号音の提示を開始するものである。
【0014】
請求項6に係る発明は、聴覚フィルタの形状をモデル化したroex(p,r)フィルタの係数pを求めて聴覚フィルタの形状を推定する装置であって、所定周波数の信号音を生成する信号音生成部と、ノッチを持たない被験者の外耳と中耳の周波数特性を考慮したノイズを生成するノイズ生成部と、このノイズ生成部で生成されたノイズを所定レベルに増幅・減衰するノイズレベル変更部と、ノイズに前記信号音の周波数を含むノッチを与えるノッチ幅変更部と、このノッチ幅変更部が出力するノッチノイズを前記信号音レベル変更部が出力する信号音に重畳するノッチノイズ重畳部と、このノッチノイズ重畳部が出力する検査音を被験者に提示する検査音提示部と、被験者が検査音を知覚できた時のノッチ幅及びその上限値に基づいてroex(p,r)フィルタの係数pを算出すると共に、得られた係数p,rからフィルタ形状を推定する聴覚フィルタ算出部と、推定したフィルタ形状を表示する聴覚フィルタ表示部を備えるものである。
【0015】
請求項7に係る発明は、請求項6に記載の聴覚フィルタの形状推定装置において、前記信号音生成部で生成された信号音を所定レベルに増幅・減衰する信号音レベル変更部を備えるものである。
【0016】
請求項8に係る発明は、請求項6又は7に記載の聴覚フィルタの形状推定装置において、被験者が検査音を知覚できた時のノッチ幅に基づいて被験者に適したノッチ幅の上限値を算出する上限値算出部を備えるものである。
【0017】
請求項9に係る発明は、請求項6、7又8に記載の聴覚フィルタの形状推定装置において、限界マスキングレベルNX’及びまたは最小ノッチ幅gX-aを測定する際に、被験者に対するマスカーの提示を開始した後に、所定の時間間隔を置いてから信号音の提示を開始するものである。
【0018】
【発明の実施の形態】
以下に本発明の実施の形態を添付図面に基づいて説明する。ここで、図1は本発明に係る聴覚フィルタの形状推定装置の構成図、図2は検査音の提示と被験者の応答の関係を示すタイミングチャート、図3は本発明に係る聴覚フィルタの形状推定方法の説明図、図4は同じく手順を示すフローチャート、図5は本発明で求めた聴覚フィルタの形状を示す図、図6はマスカーの構成図である。なお、図3におけるマスカーは1Hz当たりの音圧レベルが一様となっているが、実際には外耳と中耳の周波数特性が考慮されたレベルとなる。
【0019】
本発明に係る聴覚フィルタの形状推定装置は、図1に示すように、信号音生成部1、信号音レベル変更部2、ノイズ生成部3、ノイズレベル変更部4、ノッチ幅変更部5、ノッチノイズ重畳部6、検査音提示部7、上限値算出部8、聴覚フィルタ算出部9、聴覚フィルタ表示部10、ノッチ幅表示部11、応答部12などを備えている。
【0020】
信号音生成部1は、信号音(純音)として所定の周波数fの正弦波信号を出力する。周波数fの値は、任意に設定することができる。信号音生成部1は、CPUで構成して、所定のプログラムによって信号音を生成してもよいし、メモリで構成し、予め信号音信号を格納するようにしてもよい。
【0021】
信号音レベル変更部2は、信号音生成部1で生成された信号音を所定レベルに増幅・減衰する。信号音レベル変更部2は、ある被験者の最小可聴閾値T[dBSPL]の信号音、最小可聴閾値T[dBSPL]に任意の値x[dB]を加算したレベルT+x[dBSPL]の信号音S、レベルT+x[dBSPL]から任意の値a[dB]を減算したレベルT+x−a[dBSPL]の信号音S’などを出力する。なお、x>aとする。
【0022】
ノイズ生成部3は、内耳に到達した時に1Hz当たりの音圧レベルが一様になるように個々の被験者の外耳と中耳の周波数特性を考慮したノッチを持たないノイズを生成する。この場合には、例えば被験者の外耳と中耳の周波数特性を補正するようなマスカーとなるUEN(uniform exciting noise)を生成すればよい。ノイズ生成部3は、CPUで構成して、所定のプログラムによってノイズを生成してもよいし、メモリで構成し、予めノイズ信号を格納しておいてもよい。
【0023】
ここで、健聴者では、外耳、中耳における周波数特性の個人差は小さいため、補正値を個々に変える必要はない。
なお、健聴者における補正には、最小可聴音圧(minimum audible pressure:MAP)、最小可聴野(minimum audible field:MAF)、等感曲線(equal loudness contour:ELC)を用いて補正値を決定してマスカーを生成することも可能である。
【0024】
他に妨害する音のない状態で検知することができる音の最小レベルを指す音の絶対閾の測定方法には二つの方法がある。一つ目の方法は、外耳道の入口付近又は外耳道の内部で周波数と最小可聴音圧の関係を測定するもので、測定された閾値は最小可聴音圧(MAP)と呼ばれる。また、もう一つの方法は、大きな無響室(吸音性の高い壁で囲まれた部屋)の中でスピーカから提示される音のレベルを聴取者の頭の中心位置で測定するもので、測定された閾値は最小可聴野(MAF)と呼ばれる。
【0025】
また、音のレベルを固定した1000Hzの純音と検査音を交互に提示し、検査音の強さを変化させて大きさを純音に合わせる方法があり、これをさまざまな周波数の純音についておこなった結果得られたのが、音の大きさの等感曲線(ELC)と呼ばれる。
【0026】
一方、難聴者では外耳、中耳に障害を持つ場合があり、健聴者の外耳、中耳の周波数特性と異なることが予測できるので、難聴者の場合には個々の外耳、中耳の周波数特性を考慮した補正を行ってマスカーを生成する必要がある。
【0027】
外耳、中耳に障害がある難聴者の場合は、骨導と気導のオージオグラムの値が異なるため、この違いを利用して補正を考慮し、マスカーを生成することができる。また、外耳の周波数特性は、外耳道の入口付近の周波数特性と、鼓膜付近にプローブマイクを設けて周波数特性を測定した結果から補正を決定することも可能である。
【0028】
ノイズレベル変更部4は、ノイズ生成部3で生成されたノイズを所定レベルに増幅・減衰する。ノイズレベル変更部4は、レベルT+x[dBSPL]の信号音SをマスクすることができるレベルNX’のノイズなどを出力する。
【0029】
ノッチ幅変更部5は、ノイズレベル変更部4より出力されるノイズに中心周波数fCが信号音の周波数fと同一(fC=f)であるノッチを与える。このノッチのノッチ幅gは、被験者の応答によって随時変更される。ノッチ幅変更部5は、所望のノッチを実現するようなフィルタとして構成してもよいし、複数のノッチノイズ(マスカー)を予めメモリに格納しておいて、各種のノッチを有するノイズを随時選択的に使用してもよい。
【0030】
ノッチノイズ重畳部6は、ノッチ幅変更部5が出力する中心周波数fC、ノッチ幅g、レベルNX’のノッチノイズを、信号音レベル変更部2が出力する信号音S’に重畳して検査音とする。
検査音提示部7は、ノッチノイズ重畳部6が出力する検査音を被験者に提示する。被験者は検査音を聴取して、信号音S’を知覚できたか否かを応答する。
【0031】
上限値算出部8は、被験者が検査音を知覚できた時のノッチ幅gに基づいて被験者に適したノッチ幅の上限値gmaxを算出する。
聴覚フィルタ算出部9は、被験者が信号音S’を知覚できた時のノッチ幅g及び上限値gmaxに基づいてroex(p,r)フィルタのフィルタ係数pを算出すると共に、得られたフィルタ係数p,rから聴覚フィルタの形状を推定する。
【0032】
roex(p,r)フィルタは、次に示す式(1)で定義される。
【0033】
【数1】

Figure 0003784361
【0034】
ここで、pはフィルタのバンド幅(傾斜角度)を表す係数、rはフィルタのダイナミックレンジを表す係数、gは聴覚フィルタの中心周波数からの距離を正規化した値でノッチ幅(=Δf/fC)である。聴覚フィルタを式(1)でモデル化するためには、ノッチノイズマスキングによって得られるノッチ雑音マスキングデータPs(g)が用いられる。これは、次に示す式(2)で定義される。
【0035】
【数2】
Figure 0003784361
【0036】
ここで、gmaxはノッチ幅gの上限値、Kは各個人の感度、fCはノッチノイズの中心周波数、NXは内耳に到達した時のマスカーMのレベルである。
【0037】
聴覚フィルタ表示部10は、聴覚フィルタ算出部9によって得られたフィルタ係数p,rから聴覚フィルタの形状を表示する。
ノッチ幅表示部11は、ノッチ幅変更部5が出力するマスカーのノッチ幅gを表示する。
【0038】
応答部12は、被験者の操作によって検査音提示部7により提示される検査音から信号音S’を知覚できた場合の応答信号と知覚できない場合の応答信号を出力する。そして、信号音S’を知覚できない場合の応答信号が、ノッチ幅変更部5へ出力される。
【0039】
なお、信号音S’を知覚できた場合の応答信号のみを使用し、所定時間知覚できた場合の応答信号が出力されない時には、知覚できなかったものと判断し、次のステップに進むようにしてもよい。
【0040】
以上のように構成した聴覚フィルタの形状推定装置の動作について説明する。信号音S’に、内耳に到達した時に1Hz当たりの音圧レベルが一様になるように個々の被験者の外耳と中耳の周波数特性を考慮したマスカーMをノッチノイズ重畳部6において重畳し、検査音を作成する。作成された検査音は、検査音提示部7から被験者に提示される。
被験者は、検査音を聴取し、信号音S’を知覚できたか否かを応答部12により応答する。
【0041】
ここで、図2に示すように、被験者に対して最初の検査音を提示する際には、検査音提示の開始時に所定時間tだけ信号音レベル変更部2が出力する信号音S’のレベルを0にする(信号音S’の提示をマスカーMよりも所定時間tだけ遅らせる)ことができる。
【0042】
ノッチ幅変更部5では、その応答に応じて新たなマスカーMを作成し、再びノッチノイズ重畳部6に出力する。ここで、マスカーMのノッチ幅gの変更は、被験者の応答に応じて自動的に行うようにしてもよいし(CPU等で構成し、専用プログラムを用意する)、測定者がその都度手動で指示してもよい。
【0043】
ノッチ幅gを増加させた後に、被験者に対して新たな検査音を提示する際にも、新たな検査音提示の開始時に所定時間tだけ信号音レベル変更部2が出力する信号音S’のレベルを0にする(信号音S’の提示をマスカーMよりも所定時間tだけ遅らせる)ことができる。
【0044】
測定中のノッチ幅gは、ノッチ幅表示部11に表示される。そして、ノッチ幅変更部5は、ノッチ幅gを徐々に大きくしていき、被験者の応答により信号音S’が知覚できる最小ノッチ幅(限界ノッチ幅)gX-aを測定する。最小ノッチ幅gX-aは、図3に示すように、ノッチ雑音マスキングデータ特性において、頂点からa[dB]減衰した点のバンド幅と等価である。
【0045】
上限値算出部8では、ノッチ幅変更部5が測定した最小ノッチ幅gX-aを用いて上限値gmaxを算出する。
聴覚フィルタ算出部9では、ノッチ幅変更部5が測定した最小ノッチ幅gX-aと上限値算出部8が算出した上限値gmaxを用いてroex(p,r)フィルタのフィルタ係数pを、次のようにして算出する。ここでは、式(2)においてPs(0)とPs(gX-a)の差分aを、次に示す式(3)とする。
【0046】
【数3】
Figure 0003784361
【0047】
更に、式(2)、(3)より差分aは、次に示す式(4)となる。
【0048】
【数4】
Figure 0003784361
【0049】
ここで、gX-aの値は測定値であり、rはダイナミックレンジを定める係数xを用いて10log10r=-xとし、aは測定条件によって決まる任意の定数であるので、gmaxの値が決まれば、フィルタ係数pの値を式(4)から算出することができる。
【0050】
また、gmaxの値の決定法としては、以下のような方法が考えられる。信号検出閾値と最小可聴閾値が一致するノッチ幅gの値(以下gtと記す)は、健聴者で1Hz当たりの音圧レベルが40dBSPLになるように作成されたマスカーを使用して測定を行った場合、概ねgt=0.6〜0.7であり、gmax>0.8の区間の積分はノッチ雑音マスキングデータPs(g)の値に殆ど影響を与えない。
【0051】
従って、聴覚フィルタの過去の測定データを参照すると、gmaxの値とgtの値の関係は、次に示す式(5)となる。
【0052】
【数5】
Figure 0003784361
【0053】
ここで、Cの値は、概ね1.1〜1.3となっている。一方、本発明の簡易測定法では、測定値がノッチ幅gX-aのみであるため、信号検出閾値と最小可聴閾値が一致するノッチ幅gtの値は未知である。加えて、信号検出閾値と最小可聴閾値が一致するノッチ幅gtの値は、測定結果(ノッチ幅gと信号検出閾値[dBSPL]との関係)からダイナミックレンジ(値xに相当)によっても変化することが分かっている。よって、gmaxの値はノッチ幅gX-aからxの値ごとに推定する必要があり、gmaxの値とgX-aの値の関係は、次に示す式(6)となる。
【0054】
【数6】
Figure 0003784361
【0055】
ここで、Axはxの値に応じて変化する係数である。このAxの値は、式(5)でC=1.2と仮定し、2名の健聴被験者の測定結果(gtの値)を用いてgmaxの値を算出し、このgmaxの値とgX-aの値から実験的に決定すると、次に示す表1のようになる。
【0056】
【表1】
Figure 0003784361
【0057】
なお、表1では、x=10,15,20,30[dB]におけるgtの値を測定し、Axの値を算出しているが、例えばx=40,50[dB]におけるAxの値についても予め用意しておくことができる。
【0058】
従って、聴覚フィルタ算出部9では、上限値gmaxと最小ノッチ幅gX-aと係数rに相当する値xと値aを式(4)に代入してroex(p,r)フィルタのフィルタ係数pを算出する。更に、フィルタ係数p,rから聴覚フィルタの形状を推定する。
そして、聴覚フィルタ表示部10では、フィルタ係数pとフィルタ係数rに相当する値xからフィルタ形状を表示する。
【0059】
このようにして、フィルタ係数rに相当する値xを、例えば10,20,30,40,50[dB]として、それぞれについてフィルタ係数pの推定値を求め、図5に示すように、値xをパラメータとした聴覚フィルタの形状(横軸がノッチ幅g、縦軸が減衰量[dB])を表示する。
【0060】
次に、本発明に係る聴覚フィルタの形状推定方法について、図4に示すフローチャートより説明する。
先ず、ステップSP1で、ノイズを付加しない状態で周波数fを含む気導、骨導のオージオグラムを測定する。ここで、周波数fにおける純音最小可聴閾値をT[dBSPL/Hz]とする。ステップSP2で、測定する周波数(信号音周波数)fを設定する。
【0061】
次いで、ステップSP3で、ノッチ雑音マスキングデータのダイナミックレンジを定める係数xの値を設定する。ここでは、x=10としている。更に、ステップSP4において、周波数f[Hz]、レベルT+x[dBSPL]の純音(信号音)Sを作成する。
【0062】
次いで、ステップSP5で、信号音Sにレベルが十分小さく内耳に到達した時、1Hz当たりの音圧レベルが一様になるよう補正されたノイズやステップSP1で測定した気骨導のオージオグラムの差から考慮したノイズを重畳し、このノイズのレベルを徐々に上昇させながら、被験者が信号音Sを知覚できなくなる最小のレベル(限界マスキングレベル)NX’を測定する。
【0063】
ステップSP6では、レベルNX’のノイズのノッチ幅gを所定の値に設定する。ここでは、g=0.05としている。また、ステップSP7では、信号音SのレベルT+x[dBSPL]から減じる値a[dB]を所定の値(x>a)に設定する。ここでは、a=5としている。
【0064】
次いで、ステップSP8で、周波数f[Hz]、レベルT+x−a[dBSPL]の純音(信号音)S’を作成する。
更に、ステップSP9で、ノッチの中心周波数fC(=f)、ノッチ幅g、レベルNX’のマスカーMを作成する。
【0065】
次いで、ステップSP10で、信号音S’にマスカーMを重畳し、検査音として被験者に提示する。更に、ステップSP11で、被験者に検査音から信号音S’を知覚できたか否か判断させ、信号音S’を知覚できたという応答がない場合には、ステップSP12へ進みノッチ幅gの値を徐々に増加させ(ここでは、増加幅を0.05としている)、ステップSP9〜ステップSP11を信号音S’を知覚できたという応答があるまで繰り返す。
【0066】
ここで、図2に示すように、被験者に対して最初に検査音を提示する際には、検査音提示の開始時に所定時間tだけ信号音レベル変更部2が出力する信号音S’のレベルを0にする(信号音S’の提示をマスカーMよりも所定時間tだけ遅らせる)ことができる。同様に、ステップSP12でノッチ幅gを増加させた後に、被験者に対して新たな検査音を提示する際にも、新たな検査音提示の開始時に所定時間tだけ信号音レベル変更部2が出力する信号音S’のレベルを0にすることができる。
【0067】
そして、信号音S’の提示開始のタイミングと応答のタイミングが著しく異なる場合には、再度同条件で測定を行い、安定した応答が得られるまで測定を繰り返してもよい。また、所定時間tは一定値である必要はなく、測定回毎にランダムに変更してもよい。
【0068】
このように、最初の検査音提示の際に信号音S’をマスカーMよりも所定時間tだけ遅らせて提示したり、マスカーMのノッチ幅gの値を変化させた時に信号音S’のみを一旦停止し、所定時間tが経過した後に信号音S’を被験者に提示したりして、その時の被験者の応答のタイミングを観察すれば、被験者の応答の正当性を判断することも可能であると考えられるからである。
【0069】
図2に示すように、信号音S’が被験者に提示された後に、信号音S’を知覚したとの応答が有り、信号音S’の停止と同時に知覚したとの応答が無くなれば、被験者の応答は正常であると判断できる。一方、信号音S’が被験者に提示される前に、信号音S’を知覚したとの応答が有り、しかも信号音S’の停止後も知覚したとの応答が維持される場合などは、被験者の応答は正常でないと判断できるので、再測定を行うか、または信号音S’のレベル及びマスカーMのレベルの設定を変更すればよい。
【0070】
ステップSP11で、信号音S’を知覚できたという応答がされると、ステップSP13へ進む。
ステップSP13では、その時のノッチ幅gを、その被験者の最小ノッチ幅gX-aと定める。更に、設定した値x(=10[dB])におけるAxの値を、表1を用いて決定する(Ax=3.15)。
【0071】
次いで、ステップSP14で、式(6)にAxの値(3.15)及び測定した最小ノッチ幅gX-aの値を代入して被験者に適したノッチ幅gの上限値gmaxの値を決定する。
【0072】
ステップSP15では、式(4)に測定した最小ノッチ幅gX-aと、この最小ノッチ幅gX-aから推定したノッチ幅の上限値gmaxと、予め設定した係数rに相当する値x及び値aを代入して、roex(p,r)フィルタのフィルタ係数pX(x=10)の推定値を算出する。
【0073】
次いで、ステップSP16で、x=50であるか否かを判断し、x=50でなければ、ステップSP17に進んでx=x+10とし、x=50になるまでステップSP4〜ステップSP17を繰り返し、新たなx(x=20,30,40,50[dB])についてのフィルタ係数pXの推定値を式(4)より算出する。ここでは、xの増加幅を10とし、x=50までの測定を行っている。
【0074】
ステップSP16で、x=50になったと判断されると、ステップSP18において、図5に示すように、値xをパラメータ(x=10,20,30,40,50[dB])とした聴覚フィルタの形状(横軸がノッチ幅g、縦軸が減衰量[dB])を表示する。
【0075】
なお、本発明の実施の形態では、図6(a)に示すように、マスカーMを被験者の外耳と中耳の周波数特性を考慮したノイズにノッチを付加する方法で作成しているが、このマスカーMは、図6(b)に示すように、2つのバンドノイズ(高周波数側と低周波数側)で構成してもよい。
【0076】
また、上述した本発明の実施の形態においては、信号音レベル変更部2が出力する信号音信号の周波数fと、ノッチ幅変更部5が出力するマスカーMの中心周波数fCを一致させた(f=fC)場合について説明したが、必ずしも、一致させる必要はなく、ノッチ幅内に信号音信号が含まれていればよい。
【0077】
更に、上述した本発明の実施の形態においては、被験者が知覚できる最小ノッチ幅(限界ノッチ幅)gX-aを求めるのに際し、はじめにマスカーMのノッチ幅gを狭く設定(g=0.05)して検査音を提示し、被験者が検査音から信号音S’を知覚できない場合に、徐々にマスカーMのノッチ幅gを広くして信号音S’を知覚できるまで再度検査音を提示する手順を説明した。
【0078】
しかし、被験者が知覚できる最小ノッチ幅gX-aが求まればよいのであるから、はじめにマスカーMのノッチ幅gを広く設定して検査音を提示し、被験者が検査音から信号音S’を知覚できる状態にし、徐々にマスカーMのノッチ幅gを狭くして信号音S’を知覚できなくなるまで再度検査音を提示する手順としてもよい。
【0079】
また、上述した本発明の実施の形態においては、信号音のレベルや限界マスキングレベルNX’の設定に際し、いわゆる極限法を用いた場合について説明したが、信号音を加算した検査音と信号音を加算していない検査音のいずれかから選択させる二者強制選択法(Two alternative forced choice)等の心理物理測定法に基づいて設定してもよい。
【0080】
更に、上述した本発明の実施の形態においては、被験者が知覚できる最小ノッチ幅gX-aを求めるに際し、極限法を用いた場合について説明したが、二者強制選択法等の心理物理測定法に基づいて求めてもよい。
【0081】
また、上述した本発明の実施の形態においては、聴覚フィルタのモデルをroex(p,r)フィルタとしてフィルタ形状の推定を行ったが、このモデル関数は必ずしもroex(p,r)フィルタである必要はなく、聴覚フィルタのモデルとして適正であれば、他の関数(例えば、roex(p,r)フィルタにおいて、10log10r=0としたroex(p)など)であっても差し支えない。
【0082】
更に、上述した本発明の実施の形態においては、最小ノッチ幅gX-a測定時のステップSP10における検査音提示の開始時に、信号音S’の提示をマスカーMよりも所定時間tだけ遅らせて被験者の応答の正当性を判断しているが、ステップSP5で限界マスキングレベルNX’測定時にマスカーのレベルを変更する際にも、同様の時間遅れを信号音Sに与えて、限界マスキングレベルNX測定における正当性も併せて判断してもよい。
【0083】
【発明の効果】
以上説明したように請求項1に係る発明によれば、被験者の外耳と中耳の周波数特性を考慮した周波数fをノッチに含む、ノッチ幅g、レベルNX’のマスカーを生成することにより、難聴の診断、難聴者の聴覚特性の把握、補聴器フィッティング等を効率的かつ正確に行うために有効な聴覚フィルタの形状を正確に推定することができる。
【0084】
請求項2に係る発明によれば、入力信号レベルに対応した聴覚フィルタの形状を推定することができ、入力信号レベルと聴覚フィルタとの関係を知ることができる。
【0085】
請求項3に係る発明によれば、周波数fにおける被験者の最小可聴閾値Tに任意の値xを加算した信号音Sに基づいて限界マスキングレベルNX’を決定することによって、より正確に被験者の外耳と中耳における周波数特性を考慮した聴覚フィルタの形状を推定することができる。
【0086】
請求項4に係る発明によれば、最小ノッチ幅gX-aから被験者に適したノッチ幅gの上限値gmaxを推定することによって、より正確に被験者の外耳と中耳における周波数特性を考慮した聴覚フィルタの形状を推定することができる。
【0087】
請求項5に係る発明によれば、信号音提示のタイミングと被験者の応答のタイミングから被験者の応答の正当性を容易に判断することができる。
【0088】
請求項6に係る発明によれば、被験者の外耳と中耳の周波数特性を考慮したノイズを生成するノイズ生成部を備えることにより、難聴の診断、難聴者の聴覚特性の把握、補聴器フィッティング等を効率的かつ正確に行うために有効な聴覚フィルタの形状を正確に推定することができる。
【0089】
請求項7に係る発明によれば、信号音生成部で生成された信号音を所定レベルに増幅・減衰する信号音レベル変更部を備えることによって、より正確に聴覚フィルタの形状を推定することができる。
【0090】
請求項8に係る発明によれば、被験者が検査音を知覚できた時のノッチ幅に基づいて被験者に適したノッチ幅の上限値を算出する上限値算出部を備えることによって、より正確に聴覚フィルタの形状を推定することができる。
【0091】
請求項9に係る発明によれば、信号音提示のタイミングと被験者の応答のタイミングから被験者の応答の正当性を容易に判断することができる。
【図面の簡単な説明】
【図1】本発明に係る聴覚フィルタの形状推定装置の構成図
【図2】検査音の提示と被験者の応答の関係を示すタイミングチャート
【図3】本発明に係る聴覚フィルタの形状推定方法の説明図
【図4】本発明に係る聴覚フィルタの形状推定方法の手順を示すフローチャート
【図5】本発明で求めた聴覚フィルタの形状を示す図
【図6】マスカーの構成図
【符号の説明】
1…信号音生成部、2…信号音レベル変更部、3…ノイズ生成部、4…ノイズレベル変更部、5…ノッチ幅変更部、6…ノッチノイズ重畳部、7…検査音提示部、8…上限値算出部、9…聴覚フィルタ算出部、10…聴覚フィルタ表示部、11…ノッチ幅表示部、12…応答部。[0001]
BACKGROUND OF THE INVENTION
The present invention relates to a method and apparatus for estimating the shape of a human auditory filter.
[0002]
[Prior art]
At present, the hearing characteristics of hearing-impaired people that are most frequently used are hearing tests (audiogram measurement) and speech intelligibility tests. In the hearing test, the frequency characteristics of the minimum hearing threshold of the deaf person can be known, and in the speech intelligibility test, the hearing ability of the deaf person can be known.
However, since the hearing characteristics of the hearing-impaired person vary from person to person, it is considered that only these two methods can grasp only one part of the complicated hearing characteristics.
[0003]
In general, it is said that a hearing-impaired person has a reduced frequency resolution in addition to a decrease in hearing ability and a ability to hear words. Here, the frequency resolution is the ability to distinguish two sounds having different frequencies. For example, a normal hearing person can distinguish between two sounds whose frequencies are close to each other such as 1 kHz and 1.2 kHz, but a hearing-impaired person whose frequency resolution is low cannot distinguish between these two sounds.
[0004]
If the degree of decrease in the frequency resolution increases, the ability to hear words decreases, and the ability to discriminate speech under noise decreases. Knowing the degree of the reduction in frequency resolution is very useful in diagnosis of hearing loss, grasping of hearing characteristics of a hearing impaired person, hearing aid fitting, and the like.
[0005]
In recent years, auditory filters have been proposed as models for expressing the mechanism of frequency analysis of human hearing. This is a concept of expressing the frequency analysis mechanism of the human inner ear with a plurality of band filter banks. The shape of the individual filters (auditory filters) in this filter bank is usually measured using notch noise masking. So far, a simple method for measuring frequency resolution for a hearing impaired person using this auditory filter theory has been disclosed in Japanese Patent Laid-Open No. 6-327654. It is known that the shape of a human auditory filter can be modeled by a roex (p, r) filter.
[0006]
[Problems to be solved by the invention]
However, it is said that auditory filter measurement by notch noise masking can measure the shape of an individual subject's auditory filter with high accuracy. However, the time required for measurement is extremely long, and it is actually an otolaryngology clinical site. It was virtually impossible to measure individual hearing loss at any time in the field of hearing aid fittings. Japanese Patent Laid-Open No. 6-327654 proposes a method for measuring whether or not the frequency resolution has deteriorated in a short time. This method can determine whether or not the frequency resolution has deteriorated. It was not possible to measure the degree of deterioration, and thus the shape of the auditory filter itself.
[0007]
Furthermore, the shape of the auditory filter is said to change according to the level of the input signal. It is said that a masker (notch noise) level of about 40 dB / Hz is optimal for auditory filter measurement of a normal hearing person based on various past knowledge. The shape is said to be broad. However, not all hearing-impaired people have the same change in the shape of the auditory filter as that of normal hearing, depending on these levels. To date, there is no method for measuring the shape change characteristics of an auditory filter according to the level of an individual deaf person.
[0008]
The present invention has been made in view of such problems of the prior art, and its purpose is to efficiently and accurately perform diagnosis of hearing loss, grasp of hearing characteristics of a hearing impaired person, fitting of hearing aids, and the like. Therefore, an object of the present invention is to propose a method and apparatus for accurately estimating the shape of an auditory filter effective in a short time.
[0009]
[Means for Solving the Problems]
In order to solve the above-mentioned problem, the invention according to claim 1 is a method for estimating the shape of an auditory filter by obtaining a coefficient p of a roex (p, r) filter that models the shape of the auditory filter, and at a frequency f. A signal sound S ′ obtained by subtracting an arbitrary value a from a signal sound S obtained by adding an arbitrary value x to the minimum audible threshold T of the subject is generated, and then the frequency f considering the frequency characteristics of the outer ear and the middle ear of the subject is generated. A masker having a notch width g and a level N X ′ included in the notch is generated, and then a test sound in which the masker is superimposed on the signal sound S ′ is presented to the subject. The notch width g Xa is measured, and then, from the minimum notch width g Xa, an upper limit value g max of the notch width g suitable for the subject, the minimum notch width g Xa , the value x corresponding to the coefficient r, and the value a calculate p, then calculate And estimates the shape of the auditory filter from said value x corresponding to the coefficient p and the coefficient r.
[0010]
According to a second aspect of the present invention, in the auditory filter shape estimation method according to the first aspect, the shape of the auditory filter is estimated using the value x as a parameter.
[0011]
The invention according to claim 3 is the method for estimating the shape of an auditory filter according to claim 1 or 2, wherein the limit masking level is based on the signal sound S obtained by adding an arbitrary value x to the minimum audible threshold value T of the subject at the frequency f. N X 'is determined.
[0012]
According to a fourth aspect of the present invention, in the auditory filter shape estimation method according to the first, second, or third aspect, the upper limit value g max of the notch width g suitable for the subject is estimated from the minimum notch width g Xa. .
[0013]
According to a fifth aspect of the present invention, in the method for estimating the shape of an auditory filter according to the first, second, third, or fourth aspect, when measuring the limit masking level N X ′ and / or the minimum notch width g Xa , After the start of presentation, the signal sound presentation is started after a predetermined time interval.
[0014]
The invention according to claim 6 is an apparatus for estimating the shape of an auditory filter by obtaining a coefficient p of a roex (p, r) filter that models the shape of the auditory filter, and a signal that generates a signal sound of a predetermined frequency. Sound generator, noise generator that generates noise considering the frequency characteristics of the outer and middle ears of subjects without notches, and noise level change that amplifies and attenuates the noise generated by this noise generator to a predetermined level A notch width changing unit that gives the noise a notch including the frequency of the signal sound, and a notch noise superimposing unit that superimposes the notch noise output by the notch width changing unit on the signal sound output by the signal sound level changing unit And a test sound presenting unit that presents the test sound output by the notch noise superimposing unit to the subject, and a roex (p based on the notch width when the test subject can perceive the test sound and its upper limit value. To calculate the coefficient of r) filter p, those having coefficients obtained p, and auditory filter calculating section for estimating the filter shape from r, the auditory filter display unit for displaying the estimated filter shape.
[0015]
The invention according to claim 7 is the auditory filter shape estimation device according to claim 6, further comprising a signal sound level changing unit that amplifies and attenuates the signal sound generated by the signal sound generating unit to a predetermined level. is there.
[0016]
The invention according to claim 8 is the auditory filter shape estimation apparatus according to claim 6 or 7, wherein the upper limit value of the notch width suitable for the subject is calculated based on the notch width when the subject can perceive the test sound. The upper limit calculation part to perform is provided.
[0017]
The invention according to claim 9 is the auditory filter shape estimation device according to claim 6, 7 or 8, wherein the masker is presented to the subject when the limit masking level N X 'and / or the minimum notch width g Xa is measured. After starting, the presentation of the signal sound is started after a predetermined time interval.
[0018]
DETAILED DESCRIPTION OF THE INVENTION
Embodiments of the present invention will be described below with reference to the accompanying drawings. Here, FIG. 1 is a block diagram of the auditory filter shape estimation apparatus according to the present invention, FIG. 2 is a timing chart showing the relationship between the presentation of the test sound and the response of the subject, and FIG. 3 is the auditory filter shape estimation according to the present invention. FIG. 4 is a flowchart showing the procedure, FIG. 5 is a diagram showing the shape of the auditory filter obtained by the present invention, and FIG. 6 is a block diagram of the masker. Note that the masker in FIG. 3 has a uniform sound pressure level per 1 Hz, but in reality, it is a level that takes into account the frequency characteristics of the outer ear and the middle ear.
[0019]
As shown in FIG. 1, an auditory filter shape estimation apparatus according to the present invention includes a signal sound generator 1, a signal sound level changer 2, a noise generator 3, a noise level changer 4, a notch width changer 5, a notch A noise superimposing unit 6, an examination sound presenting unit 7, an upper limit calculating unit 8, an auditory filter calculating unit 9, an auditory filter display unit 10, a notch width display unit 11, a response unit 12 and the like are provided.
[0020]
The signal sound generation unit 1 outputs a sine wave signal having a predetermined frequency f as a signal sound (pure sound). The value of the frequency f can be set arbitrarily. The signal sound generation unit 1 may be configured by a CPU and generate a signal sound by a predetermined program, or may be configured by a memory and store a signal sound signal in advance.
[0021]
The signal sound level changing unit 2 amplifies and attenuates the signal sound generated by the signal sound generating unit 1 to a predetermined level. The signal sound level changing unit 2 is a signal sound having a minimum audible threshold value T [dBSPL] of a subject, a signal sound S having a level T + x [dBSPL] obtained by adding an arbitrary value x [dB] to the minimum audible threshold value T [dBSPL]. A signal sound S ′ of level T + x−a [dBSPL] obtained by subtracting an arbitrary value a [dB] from level T + x [dBSPL] is output. Note that x> a.
[0022]
The noise generation unit 3 generates noise having no notch considering the frequency characteristics of the outer ear and the middle ear of each subject so that the sound pressure level per 1 Hz becomes uniform when reaching the inner ear. In this case, for example, a uniform exciting noise (UEN) that is a masker that corrects the frequency characteristics of the outer ear and the middle ear of the subject may be generated. The noise generation unit 3 may be configured by a CPU and generate noise by a predetermined program, or may be configured by a memory and store a noise signal in advance.
[0023]
Here, in a normal hearing person, since the individual difference of the frequency characteristic in an outer ear and a middle ear is small, it is not necessary to change a correction value separately.
For normal hearing, correction values are determined using minimum audible pressure (MAP), minimum audible field (MAF), and equal loudness contour (ELC). It is also possible to generate a masker.
[0024]
There are two methods for measuring the absolute threshold of sound, which indicates the minimum level of sound that can be detected without any other disturbing sound. The first method is to measure the relationship between the frequency and the minimum audible sound pressure near the entrance of the ear canal or inside the ear canal, and the measured threshold is called the minimum audible sound pressure (MAP). Another method is to measure the level of sound presented from the speaker in a large anechoic room (a room surrounded by walls with high sound absorption) at the center of the listener's head. The threshold value is called the minimum audible field (MAF).
[0025]
In addition, there is a method in which 1000Hz pure sound with fixed sound level and inspection sound are presented alternately, and the intensity of the inspection sound is changed to match the volume with the pure sound. What is obtained is called an isosensitive curve (ELC) of loudness.
[0026]
On the other hand, the hearing impaired may have disabilities in the outer and middle ears, which can be predicted to be different from the frequency characteristics of the outer ear and middle ear of a normal hearing person. It is necessary to generate a masker by performing correction considering the above.
[0027]
In the case of a hearing-impaired person with a disability in the outer ear or the middle ear, the audiogram values of the bone conduction and the air conduction differ, so that the difference can be used to generate a masker by taking correction into account. Further, the correction of the frequency characteristics of the outer ear can be determined from the frequency characteristics near the entrance of the ear canal and the result of measuring the frequency characteristics by providing a probe microphone near the eardrum.
[0028]
The noise level changing unit 4 amplifies and attenuates the noise generated by the noise generating unit 3 to a predetermined level. The noise level changing unit 4 outputs noise of level N X ′ that can mask the signal sound S of level T + x [dBSPL].
[0029]
The notch width changing unit 5 gives the noise output from the noise level changing unit 4 a notch whose center frequency f C is the same as the frequency f of the signal sound (f C = f). The notch width g of this notch is changed at any time according to the response of the subject. The notch width changing unit 5 may be configured as a filter that realizes a desired notch, or a plurality of notch noises (maskers) are stored in a memory in advance, and noise having various notches is selected as needed. May be used.
[0030]
The notch noise superimposing unit 6 superimposes the notch noise of the center frequency f C , the notch width g, and the level N X ′ output from the notch width changing unit 5 on the signal sound S ′ output from the signal sound level changing unit 2. Use inspection sound.
The inspection sound presenting unit 7 presents the inspection sound output by the notch noise superimposing unit 6 to the subject. The subject listens to the test sound and responds whether or not the signal sound S ′ can be perceived.
[0031]
The upper limit calculator 8 calculates the upper limit g max of the notch width suitable for the subject based on the notch width g when the subject can perceive the test sound.
The auditory filter calculation unit 9 calculates the filter coefficient p of the roex (p, r) filter based on the notch width g and the upper limit g max when the subject can perceive the signal sound S ′, and the obtained filter The shape of the auditory filter is estimated from the coefficients p and r.
[0032]
The roex (p, r) filter is defined by the following equation (1).
[0033]
[Expression 1]
Figure 0003784361
[0034]
Here, p is a coefficient representing the bandwidth (tilt angle) of the filter, r is a coefficient representing the dynamic range of the filter, g is a value obtained by normalizing the distance from the center frequency of the auditory filter, and the notch width (= Δf / f C ). In order to model the auditory filter by the equation (1), notch noise masking data Ps (g) obtained by notch noise masking is used. This is defined by the following equation (2).
[0035]
[Expression 2]
Figure 0003784361
[0036]
Here, g max is the upper limit value of the notch width g, K is the sensitivity of each individual, f C is the center frequency of the notch noise, and N X is the level of the masker M when it reaches the inner ear.
[0037]
The auditory filter display unit 10 displays the shape of the auditory filter from the filter coefficients p and r obtained by the auditory filter calculation unit 9.
The notch width display unit 11 displays the notch width g of the masker output from the notch width changing unit 5.
[0038]
The response unit 12 outputs a response signal when the signal sound S ′ can be perceived from the test sound presented by the test sound presenting unit 7 by the operation of the subject and a response signal when the signal sound S ′ cannot be perceived. Then, a response signal when the signal sound S ′ cannot be perceived is output to the notch width changing unit 5.
[0039]
It should be noted that only the response signal when the signal sound S ′ can be perceived is used, and when the response signal when the signal sound S ′ can be perceived for a predetermined time is not output, it is determined that the signal sound S ′ cannot be perceived and the process proceeds to the next step. .
[0040]
The operation of the auditory filter shape estimation apparatus configured as described above will be described. The notch noise superimposing unit 6 superimposes the masker M in consideration of the frequency characteristics of the outer ear and the middle ear of each subject so that the sound pressure level per 1 Hz is uniform when the signal sound S ′ reaches the inner ear. Create inspection sound. The created inspection sound is presented to the subject from the inspection sound presentation unit 7.
The subject listens to the test sound and responds by the response unit 12 as to whether or not the signal sound S ′ has been perceived.
[0041]
Here, as shown in FIG. 2, when the first test sound is presented to the subject, the level of the signal sound S ′ output by the signal sound level changing unit 2 for a predetermined time t at the start of the test sound presentation. Can be set to 0 (presentation of the signal sound S ′ is delayed from the masker M by a predetermined time t).
[0042]
The notch width changing unit 5 creates a new masker M according to the response and outputs it again to the notch noise superimposing unit 6. Here, the change of the notch width g of the masker M may be automatically performed according to the response of the subject (configured by a CPU or the like, and a dedicated program is prepared), or manually by the measurer each time. You may instruct.
[0043]
Even when a new examination sound is presented to the subject after increasing the notch width g, the signal sound S ′ output by the signal sound level changing unit 2 for a predetermined time t at the start of the presentation of the new examination sound. The level can be set to 0 (presentation of the signal sound S ′ can be delayed by a predetermined time t from the masker M).
[0044]
The notch width g being measured is displayed on the notch width display section 11. Then, the notch width changing unit 5 gradually increases the notch width g, and measures the minimum notch width (limit notch width) g Xa at which the signal sound S ′ can be perceived by the response of the subject. As shown in FIG. 3, the minimum notch width g Xa is equivalent to the bandwidth of a point attenuated by a [dB] from the apex in the notch noise masking data characteristic.
[0045]
The upper limit calculation unit 8 calculates the upper limit g max using the minimum notch width g Xa measured by the notch width changing unit 5.
The auditory filter calculation unit 9 uses the minimum notch width g Xa measured by the notch width changing unit 5 and the upper limit value g max calculated by the upper limit value calculation unit 8 to calculate the filter coefficient p of the roex (p, r) filter as follows: Calculate as follows. Here, the difference a between Ps (0) and Ps (g Xa ) in the equation (2) is defined as the following equation (3).
[0046]
[Equation 3]
Figure 0003784361
[0047]
Furthermore, the difference a is expressed by the following expression (4) from the expressions (2) and (3).
[0048]
[Expression 4]
Figure 0003784361
[0049]
Here, the value of g Xa is a measure, r is a 10 log 10 r = -x using coefficients x defining the dynamic range, since a is an optional constant determined by the measurement conditions, the value of g max is If determined, the value of the filter coefficient p can be calculated from the equation (4).
[0050]
Further, as a method for determining the value of g max , the following method can be considered. Signal detection threshold and minimum audible threshold (hereinafter referred to as g t) value matching notch width g is subjected to measurement using a masker sound pressure level per 1Hz is created so as to 40dBSPL in hearing person In this case, g t = 0.6 to 0.7, and integration in the interval of g max > 0.8 has little effect on the value of the notch noise masking data Ps (g).
[0051]
Therefore, referring to past measurement data of the auditory filter, the relationship between the value of g max and the value of g t is expressed by the following equation (5).
[0052]
[Equation 5]
Figure 0003784361
[0053]
Here, the value of C is approximately 1.1 to 1.3. On the other hand, in the simple measurement method of the present invention, since the measured value is only the notch width g Xa , the value of the notch width g t at which the signal detection threshold matches the minimum audible threshold is unknown. In addition, the value of the notch width g t to the signal detection threshold and minimum audible threshold matches the measurement result changes from (relationship between notch width g and the signal detection threshold [dBSPL]) by the dynamic range (corresponding to a value x) I know you will. Therefore, the value of g max needs to be estimated for each value of x from the notch width g Xa, and the relationship between the value of g max and the value of g Xa is expressed by the following equation (6).
[0054]
[Formula 6]
Figure 0003784361
[0055]
Here, A x is a coefficient that changes according to the value of x. The value of this A x assumes that C = 1.2 in equation (5), calculates the value of g max using 2 people hearing subjects the measurement result (the value of g t), the g max When experimentally determined from the value and the value of g Xa, the following Table 1 is obtained.
[0056]
[Table 1]
Figure 0003784361
[0057]
In Table 1, x = 10, 15, 20, 30 to measure the value of g t in [dB], but calculates the value of A x, for example x = 40, 50 A in [dB] x The value of can also be prepared in advance.
[0058]
Therefore, the auditory filter calculation unit 9 substitutes the upper limit value g max , the minimum notch width g Xa , the value x and the value a corresponding to the coefficient r into the equation (4), and the filter coefficient p of the roex (p, r) filter. Is calculated. Further, the shape of the auditory filter is estimated from the filter coefficients p and r.
The auditory filter display unit 10 displays the filter shape from the value x corresponding to the filter coefficient p and the filter coefficient r.
[0059]
In this way, the value x corresponding to the filter coefficient r is set to 10, 20, 30, 40, 50 [dB], for example, and an estimated value of the filter coefficient p is obtained for each, and as shown in FIG. The shape of the auditory filter with the parameter as a parameter (notch width g on the horizontal axis and attenuation [dB] on the vertical axis) is displayed.
[0060]
Next, the auditory filter shape estimation method according to the present invention will be described with reference to the flowchart shown in FIG.
First, in step SP1, the air conduction and bone conduction audiograms including the frequency f are measured without adding noise. Here, the pure tone minimum audible threshold value at the frequency f is T [dBSPL / Hz]. In step SP2, a frequency (signal sound frequency) f to be measured is set.
[0061]
Next, in step SP3, the value of the coefficient x that determines the dynamic range of the notch noise masking data is set. Here, x = 10. Further, in step SP4, a pure sound (signal sound) S having a frequency f [Hz] and a level T + x [dBSPL] is created.
[0062]
Next, in step SP5, when the level of the signal sound S is sufficiently small and reaches the inner ear, the noise corrected to make the sound pressure level per 1 Hz uniform and the difference in the air bone conduction audiogram measured in step SP1. The minimum level (limit masking level) N X ′ at which the subject cannot perceive the signal sound S is measured while superimposing the considered noise and gradually increasing the noise level.
[0063]
In step SP6, the noise notch width g at level N X 'is set to a predetermined value. Here, g = 0.05. In step SP7, the value a [dB] subtracted from the level T + x [dBSPL] of the signal sound S is set to a predetermined value (x> a). Here, a = 5.
[0064]
Next, in step SP8, a pure tone (signal tone) S ′ having a frequency f [Hz] and a level T + x−a [dBSPL] is created.
Further, in step SP9, a masker M having a notch center frequency f C (= f), a notch width g, and a level N X ′ is created.
[0065]
Next, in step SP10, the masker M is superimposed on the signal sound S ′ and presented as a test sound to the subject. Further, in step SP11, the subject is made to judge whether or not the signal sound S ′ can be perceived from the inspection sound. If there is no response that the signal sound S ′ can be perceived, the process proceeds to step SP12 and the value of the notch width g is set. It is gradually increased (in this case, the increase range is set to 0.05), and steps SP9 to SP11 are repeated until there is a response that the signal sound S ′ can be perceived.
[0066]
Here, as shown in FIG. 2, when the test sound is first presented to the subject, the level of the signal sound S ′ output by the signal sound level changing unit 2 for a predetermined time t at the start of the test sound presentation. Can be set to 0 (presentation of the signal sound S ′ is delayed from the masker M by a predetermined time t). Similarly, when the new test sound is presented to the subject after the notch width g is increased in step SP12, the signal sound level changing unit 2 outputs only the predetermined time t at the start of the new test sound presentation. The level of the signal sound S ′ to be performed can be reduced to zero.
[0067]
If the timing of starting the presentation of the signal sound S ′ and the response timing are significantly different, the measurement may be performed again under the same conditions, and the measurement may be repeated until a stable response is obtained. Further, the predetermined time t does not need to be a constant value, and may be changed randomly every measurement.
[0068]
In this way, when the first inspection sound is presented, the signal sound S ′ is presented after being delayed by a predetermined time t from the masker M, or only the signal sound S ′ is displayed when the value of the notch width g of the masker M is changed. It is also possible to judge the validity of the response of the subject by temporarily stopping and presenting the signal sound S ′ to the subject after the predetermined time t has elapsed and observing the response timing of the subject at that time. Because it is considered.
[0069]
As shown in FIG. 2, after the signal sound S ′ is presented to the subject, there is a response that the signal sound S ′ is perceived, and if there is no response that the signal sound S ′ is perceived simultaneously with the stop of the signal sound S ′, Can be determined to be normal. On the other hand, when there is a response that the signal sound S ′ is perceived before the signal sound S ′ is presented to the subject, and the response that the signal sound S ′ is perceived after the stop of the signal sound S ′ is maintained, etc. Since it can be determined that the response of the subject is not normal, remeasurement may be performed, or the settings of the level of the signal sound S ′ and the level of the masker M may be changed.
[0070]
If a response that the signal sound S ′ can be perceived in step SP11, the process proceeds to step SP13.
In step SP13, the notch width g at that time is determined as the minimum notch width g Xa of the subject. Further, the value of A x at the set value x (= 10 [dB]) is determined using Table 1 (A x = 3.15).
[0071]
Next, at step SP14, the value of A x (3.15) and the value of the measured minimum notch width g Xa are substituted into Equation (6) to determine the upper limit value g max of the notch width g suitable for the subject. To do.
[0072]
In step SP15, the minimum notch width g Xa measured in the equation (4), the upper limit value g max of the notch width estimated from the minimum notch width g Xa , the value x and the value a corresponding to the preset coefficient r are obtained. By substituting, an estimated value of the filter coefficient p X (x = 10) of the roex (p, r) filter is calculated.
[0073]
Next, in step SP16, it is determined whether or not x = 50. If x = 50 is not satisfied, the process proceeds to step SP17 to set x = x + 10, and steps SP4 to SP17 are repeated until x = 50. An estimated value of the filter coefficient p X for x (x = 20, 30, 40, 50 [dB]) is calculated from the equation (4). Here, the increase width of x is set to 10, and measurement up to x = 50 is performed.
[0074]
If it is determined in step SP16 that x = 50, as shown in FIG. 5, an auditory filter using the value x as a parameter (x = 10, 20, 30, 40, 50 [dB]) in step SP18. (The horizontal axis is the notch width g and the vertical axis is the attenuation [dB]).
[0075]
In the embodiment of the present invention, as shown in FIG. 6A, the masker M is created by a method of adding notches to noise considering the frequency characteristics of the outer ear and middle ear of the subject. The masker M may be composed of two band noises (high frequency side and low frequency side) as shown in FIG.
[0076]
In the above-described embodiment of the present invention, the frequency f of the signal sound signal output from the signal sound level changing unit 2 and the center frequency f C of the masker M output from the notch width changing unit 5 are matched ( The case of f = f C ) has been described, but it is not always necessary to match, and it is only necessary that the signal sound signal is included within the notch width.
[0077]
Further, in the above-described embodiment of the present invention, when determining the minimum notch width (limit notch width) g Xa that can be perceived by the subject, first, the notch width g of the masker M is set narrow (g = 0.05). The test sound is presented, and when the subject cannot perceive the signal sound S ′ from the test sound, the procedure of presenting the test sound again until the signal sound S ′ can be perceived by gradually increasing the notch width g of the masker M. explained.
[0078]
However, since it is only necessary to obtain the minimum notch width g Xa that can be perceived by the subject, first, the test sound is presented by setting the notch width g of the masker M wide, and the subject can perceive the signal sound S ′ from the test sound. The procedure may be such that the test sound is presented again until the notch width g of the masker M is gradually narrowed until the signal sound S ′ cannot be perceived.
[0079]
In the above-described embodiment of the present invention, the case where the so-called limit method is used for setting the signal sound level and the limit masking level N X ′ has been described. You may set based on psychophysical measurement methods, such as the two alternative forced choice method (Two alternative forced choice) which makes it select from either of the test | inspection sound which has not added.
[0080]
Furthermore, in the above-described embodiment of the present invention, the case where the limit method is used to determine the minimum notch width g Xa that can be perceived by the subject has been described. However, based on a psychophysical measurement method such as a two-party forced selection method. You may ask.
[0081]
In the above-described embodiment of the present invention, the filter shape is estimated using the model of the auditory filter as the roex (p, r) filter. However, the model function is not necessarily a roex (p, r) filter. However, as long as it is appropriate as a model of an auditory filter, other functions (for example, roex (p) with 10 log 10 r = 0 in the roex (p, r) filter) may be used.
[0082]
Further, in the above-described embodiment of the present invention, the presentation of the signal sound S ′ is delayed from the masker M by a predetermined time t at the start of the inspection sound presentation in step SP10 when measuring the minimum notch width g Xa . Although the legitimacy of the response is judged, when the masker level is changed at the time of measuring the limit masking level N X ′ at step SP5, a similar time delay is given to the signal sound S to measure the limit masking level N X. You may also judge the legitimacy.
[0083]
【The invention's effect】
As described above, according to the first aspect of the invention, by generating a masker with notch width g and level N X ′ including frequency f in consideration of frequency characteristics of the subject's outer ear and middle ear, It is possible to accurately estimate the shape of an auditory filter effective for efficiently and accurately performing hearing loss diagnosis, hearing characteristics of a hearing impaired person, hearing aid fitting, and the like.
[0084]
According to the second aspect of the present invention, the shape of the auditory filter corresponding to the input signal level can be estimated, and the relationship between the input signal level and the auditory filter can be known.
[0085]
According to the invention of claim 3, by determining the limit masking level N X ′ based on the signal sound S obtained by adding an arbitrary value x to the minimum audible threshold T of the subject at the frequency f, It is possible to estimate the shape of the auditory filter in consideration of frequency characteristics in the outer ear and the middle ear.
[0086]
According to the fourth aspect of the invention, the upper limit value g max of the notch width g suitable for the subject is estimated from the minimum notch width g Xa, thereby more accurately considering the frequency characteristics in the outer ear and the middle ear of the subject. The shape of the filter can be estimated.
[0087]
According to the invention which concerns on Claim 5, the validity of a test subject's response can be judged easily from the timing of signal sound presentation and the test subject's response timing.
[0088]
According to the invention of claim 6, by providing a noise generation unit that generates noise in consideration of the frequency characteristics of the outer ear and middle ear of the subject, diagnosis of hearing loss, grasp of hearing characteristics of the deaf person, hearing aid fitting, etc. It is possible to accurately estimate the shape of an effective auditory filter for efficient and accurate execution.
[0089]
According to the seventh aspect of the invention, the shape of the auditory filter can be estimated more accurately by including the signal sound level changing unit that amplifies and attenuates the signal sound generated by the signal sound generating unit to a predetermined level. it can.
[0090]
According to the eighth aspect of the present invention, by providing the upper limit value calculation unit that calculates the upper limit value of the notch width suitable for the subject based on the notch width when the subject can perceive the test sound, the auditory sense can be more accurately heard. The shape of the filter can be estimated.
[0091]
According to the invention which concerns on Claim 9, the validity of a test subject's response can be judged easily from the timing of signal sound presentation and the test subject's response timing.
[Brief description of the drawings]
FIG. 1 is a configuration diagram of an auditory filter shape estimating apparatus according to the present invention. FIG. 2 is a timing chart showing the relationship between the presentation of a test sound and the response of a subject. FIG. 4 is a flowchart showing the procedure of the method for estimating the shape of the auditory filter according to the present invention. FIG. 5 is a diagram showing the shape of the auditory filter obtained by the present invention.
DESCRIPTION OF SYMBOLS 1 ... Signal sound production | generation part, 2 ... Signal sound level change part, 3 ... Noise production | generation part, 4 ... Noise level change part, 5 ... Notch width change part, 6 ... Notch noise superimposition part, 7 ... Test sound presentation part, 8 ... upper limit calculation unit, 9 ... auditory filter calculation unit, 10 ... auditory filter display unit, 11 ... notch width display unit, 12 ... response unit.

Claims (9)

聴覚フィルタの形状をモデル化したroex(p,r)フィルタの係数pを求めて聴覚フィルタの形状を推定する方法であって、周波数fにおける被験者の最小可聴閾値Tに任意の値xを加算した信号音Sから任意の値aを減じた信号音S’を生成し、次いで被験者の外耳と中耳の周波数特性を考慮した前記周波数fをノッチに含む、ノッチ幅g、レベルNX’のマスカーを生成し、次いで前記信号音S’に前記マスカーを重畳した検査音を被験者に提示し、ノッチ幅gを変化させながら、その被験者の最小ノッチ幅gX-aを測定し、次いで最小ノッチ幅gX-aから被験者に適したノッチ幅gの上限値gmaxと前記最小ノッチ幅gX-aと係数rに相当する前記値xと前記値aとから係数pを算出し、次いで算出した係数pと係数rに相当する前記値xとから聴覚フィルタの形状を推定することを特徴とする聴覚フィルタの形状推定方法。A method for estimating the shape of an auditory filter by obtaining a coefficient p of a roex (p, r) filter that models the shape of an auditory filter, wherein an arbitrary value x is added to the minimum audible threshold T of a subject at a frequency f. A signal sound S ′ obtained by subtracting an arbitrary value a from the signal sound S, and then a masker having a notch width g and a level N X ′ including the frequency f considering the frequency characteristics of the outer ear and the middle ear of the subject. Next, the test sound in which the masker is superimposed on the signal sound S ′ is presented to the subject, and the minimum notch width g Xa of the subject is measured while changing the notch width g, and then the minimum notch width g Xa The coefficient p is calculated from the upper limit value g max of the notch width g suitable for the subject, the minimum notch width g Xa , the value x corresponding to the coefficient r, and the value a, and then the calculated coefficient p and coefficient r Listen from the corresponding value x Shape estimation method of auditory filter and estimates the shape of the filter. 請求項1に記載の聴覚フィルタの形状推定方法において、前記値xをパラメータとして聴覚フィルタの形状を推定することを特徴とする聴覚フィルタの形状推定方法。2. The method for estimating an auditory filter shape according to claim 1, wherein the shape of the auditory filter is estimated using the value x as a parameter. 請求項1又は2に記載の聴覚フィルタの形状推定方法において、周波数fにおける被験者の最小可聴閾値Tに任意の値xを加算した信号音Sに基づいて限界マスキングレベルNX’を決定することを特徴とする聴覚フィルタの形状推定方法。3. The method of estimating a shape of an auditory filter according to claim 1, wherein the limit masking level N X ′ is determined based on a signal sound S obtained by adding an arbitrary value x to the minimum audible threshold value T of the subject at the frequency f. A method for estimating the shape of a characteristic auditory filter. 請求項1、2又は3に記載の聴覚フィルタの形状推定方法において、最小ノッチ幅gX-aから被験者に適したノッチ幅gの上限値gmaxを推定することを特徴とする聴覚フィルタの形状推定方法。4. The method for estimating the shape of an auditory filter according to claim 1, wherein the upper limit value g max of the notch width g suitable for the subject is estimated from the minimum notch width g Xa. . 請求項1、2、3又は4に記載の聴覚フィルタの形状推定方法において、限界マスキングレベルNX’及びまたは最小ノッチ幅gX-aを測定する際に、被験者に対するマスカーの提示を開始した後に、所定の時間間隔を置いてから信号音の提示を開始することを特徴とする聴覚フィルタの形状推定方法。5. The shape estimation method for an auditory filter according to claim 1, 2, 3 or 4, wherein when the limit masking level N X ′ and / or the minimum notch width g Xa is measured, the masker is presented to the subject after starting. A method for estimating the shape of an auditory filter, wherein the presentation of a signal sound is started after a time interval of. 聴覚フィルタの形状をモデル化したroex(p,r)フィルタの係数pを求めて聴覚フィルタの形状を推定する装置であって、所定周波数の信号音を生成する信号音生成部と、被験者の外耳と中耳の周波数特性を考慮したノッチ幅を持たないノイズを生成するノイズ生成部と、このノイズ生成部で生成されたノイズを所定レベルに増幅・減衰するノイズレベル変更部と、ノイズに前記信号音の周波数を含むノッチを与えるノッチ幅変更部と、このノッチ幅変更部が出力するノッチノイズを前記信号音レベル変更部が出力する信号音に重畳するノッチノイズ重畳部と、このノッチノイズ重畳部が出力する検査音を被験者に提示する検査音提示部と、被験者が検査音を知覚できた時のノッチ幅及びその上限値に基づいてroex(p,r)フィルタの係数pを算出すると共に、得られた係数p,rからフィルタ形状を推定する聴覚フィルタ算出部と、推定したフィルタ形状を表示する聴覚フィルタ表示部を備えることを特徴とする聴覚フィルタの形状推定装置。A device for estimating a shape of an auditory filter by obtaining a coefficient p of a roex (p, r) filter that models the shape of an auditory filter, a signal sound generating unit that generates a signal sound of a predetermined frequency, and an outer ear of a subject A noise generation unit that generates noise not having a notch width in consideration of the frequency characteristics of the middle ear, a noise level change unit that amplifies and attenuates the noise generated by the noise generation unit to a predetermined level, and the signal to noise A notch width changing unit for providing a notch including the frequency of the sound, a notch noise superimposing unit for superimposing the notch noise output by the notch width changing unit on the signal sound output by the signal sound level changing unit, and the notch noise superimposing unit A test sound presenting unit that presents the test sound output by the subject to the test subject, and a coefficient of the roex (p, r) filter based on the notch width and the upper limit value when the test subject can perceive the test sound Calculates the coefficients resulting p, and auditory filter calculating section for estimating the filter shape from r, the shape estimation apparatus of the auditory filter characterized in that it comprises an auditory filter display unit for displaying the estimated filter shape. 請求項6に記載の聴覚フィルタの形状推定装置において、前記信号音生成部で生成された信号音を所定レベルに増幅・減衰する信号音レベル変更部を備えることを特徴とする聴覚フィルタの形状推定装置。The shape estimation device for an auditory filter according to claim 6, further comprising a signal sound level changing unit that amplifies and attenuates the signal sound generated by the signal sound generation unit to a predetermined level. apparatus. 請求項6又は7に記載の聴覚フィルタの形状推定装置において、被験者が検査音を知覚できた時のノッチ幅に基づいて被験者に適したノッチ幅の上限値を算出する上限値算出部を備えることを特徴とする聴覚フィルタの形状推定装置。The shape estimation device for an auditory filter according to claim 6 or 7, further comprising an upper limit value calculation unit that calculates an upper limit value of a notch width suitable for the subject based on the notch width when the subject can perceive the test sound. An auditory filter shape estimation device characterized by the above. 請求項6、7又8に記載の聴覚フィルタの形状推定装置において、限界マスキングレベルNX’及びまたは最小ノッチ幅gX-aを測定する際に、被験者に対するマスカーの提示を開始した後に、所定の時間間隔を置いてから信号音の提示を開始することを特徴とする聴覚フィルタの形状推定装置。9. The shape estimation apparatus for an auditory filter according to claim 6, 7 or 8, wherein when a limit masking level N X ′ and / or a minimum notch width g Xa is measured, a predetermined time has elapsed after the masker is presented to the subject. An auditory filter shape estimation apparatus, wherein presentation of a signal sound is started after an interval.
JP2002286594A 2002-01-22 2002-09-30 Auditory filter shape estimation method and apparatus Expired - Fee Related JP3784361B2 (en)

Priority Applications (3)

Application Number Priority Date Filing Date Title
JP2002286594A JP3784361B2 (en) 2002-09-30 2002-09-30 Auditory filter shape estimation method and apparatus
US10/345,041 US7048692B2 (en) 2002-01-22 2003-01-15 Method and apparatus for estimating auditory filter shape
DE10301677A DE10301677B4 (en) 2002-01-22 2003-01-17 A method of measuring frequency selectivity and method and apparatus for estimating forms of acoustic filter by a frequency selectivity measurement method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2002286594A JP3784361B2 (en) 2002-09-30 2002-09-30 Auditory filter shape estimation method and apparatus

Publications (2)

Publication Number Publication Date
JP2004121338A JP2004121338A (en) 2004-04-22
JP3784361B2 true JP3784361B2 (en) 2006-06-07

Family

ID=32279616

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002286594A Expired - Fee Related JP3784361B2 (en) 2002-01-22 2002-09-30 Auditory filter shape estimation method and apparatus

Country Status (1)

Country Link
JP (1) JP3784361B2 (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009276248A (en) * 2008-05-15 2009-11-26 Mitsubishi Electric Corp Laser radar device
US8571245B2 (en) 2009-06-16 2013-10-29 Panasonic Corporation Hearing assistance suitability determining device, hearing assistance adjustment system, and hearing assistance suitability determining method

Also Published As

Publication number Publication date
JP2004121338A (en) 2004-04-22

Similar Documents

Publication Publication Date Title
JP3448586B2 (en) Sound measurement method and system considering hearing impairment
KR20090060598A (en) Method and apparatus of the hair cell stimulation using acoustic signal
Kimlinger et al. High-frequency audibility: The effects of audiometric configuration, stimulus type, and device
EP3384686A1 (en) Self-fitting of a hearing device
CN107320109B (en) Frequency identification test method
TW201019908A (en) Method for detecting hearing impairment and device thereof
JP2532007B2 (en) Hearing aid fitting device
Ho et al. Reference equivalent threshold sound pressure levels for Apple EarPods
JP2004179965A (en) Automatic amplification characteristic adjusting device for hearing aid
JP3784361B2 (en) Auditory filter shape estimation method and apparatus
US7048692B2 (en) Method and apparatus for estimating auditory filter shape
JP3784360B2 (en) Auditory filter shape estimation method and apparatus
JP4909263B2 (en) Judgment method of subjective characteristics of binaural sound signal
WO2010139760A2 (en) System and method for conducting an alternative forced choice hearing test
JP3741651B2 (en) Auditory filter shape estimation method and apparatus
JP3621028B2 (en) Auditory filter shape estimation device
US10258260B2 (en) Method of testing hearing and a hearing test system
US9474473B2 (en) Method to measure real-ear-to-coupler difference
KR20090065749A (en) Hearing aid and method for audiometry thereof
US10842418B2 (en) Method and apparatus for tinnitus evaluation with test sound automatically adjusted for loudness
JP4164035B2 (en) Auditory filter sensitivity determination device
Navarro Effects of ear canal occlusion and masking on the perception of voice
KR100917714B1 (en) Observation device for hearing and control method thereof
JP3579866B2 (en) Frequency resolution measurement device
Beynon When is a decibel not a decibel?: The application of decibel scales and calibration in clinical audiology

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20040202

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20050228

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20051206

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20060314

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20060314

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110324

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120324

Year of fee payment: 6

LAPS Cancellation because of no payment of annual fees