JP2012060295A - 伝達関数調整装置、伝達関数調整プログラムおよび伝達関数調整方法 - Google Patents
伝達関数調整装置、伝達関数調整プログラムおよび伝達関数調整方法 Download PDFInfo
- Publication number
- JP2012060295A JP2012060295A JP2010199773A JP2010199773A JP2012060295A JP 2012060295 A JP2012060295 A JP 2012060295A JP 2010199773 A JP2010199773 A JP 2010199773A JP 2010199773 A JP2010199773 A JP 2010199773A JP 2012060295 A JP2012060295 A JP 2012060295A
- Authority
- JP
- Japan
- Prior art keywords
- transfer function
- peak
- dip
- frequency
- frequency amplitude
- 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
Links
- 238000012546 transfer Methods 0.000 title claims abstract description 161
- 238000000034 method Methods 0.000 title claims description 21
- 238000012937 correction Methods 0.000 claims abstract description 98
- 238000000605 extraction Methods 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 153
- 238000012545 processing Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 2
- 239000000284 extract Substances 0.000 abstract description 2
- 230000004044 response Effects 0.000 description 19
- 238000005070 sampling Methods 0.000 description 4
- 238000013459 approach Methods 0.000 description 3
- 230000004807 localization Effects 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Landscapes
- Stereophonic System (AREA)
- Circuit For Audible Band Transducer (AREA)
Abstract
【解決手段】伝達関数の周波数振幅特性を算出する周波数振幅特性算出手段と、周波数振幅の平均値を算出する周波数振幅平均値算出手段と、閾値を超えるピーク(ディップ)を周波数振幅特性から抽出するピーク・ディップ抽出手段と、全ピーク(ディップ)のうち周波数振幅の平均値からの乖離が最大である1つのピーク(ディップ)を有する伝達関数を補正対象伝達関数と決定し、当該ピーク(ディップ)を補正対象として決定する補正対象決定手段と、補正対象に対応した周波数の正規化周波数との差が最小である偏角を有する伝達関数の極および零点の一方を制御対象に決定する制御対象決定手段と、制御対象を補正対象が閾値内に収まる最終位置まで移動させる制御対象移動手段と、を備える。
【選択図】図2
Description
すなわち、上記方法は時間領域あるいは周波数領域で平均化を行いピークおよびディップの除去をするものであり、伝達関数の極および零点の位置は考慮外であるためピークおよびディップを除去した結果伝達関数の極がZ平面の原点を中心とする単位円外に移動する可能性もあり、伝達関数の安定性は保証されない。
前記補正対象がピークであり、前記制御対象が極である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点から遠ざけ、
前記補正対象がピークであり、前記制御対象が零点である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点に近付け、
前記補正対象がディップであり、前記制御対象が極である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点に近付け、
前記補正対象がディップであり、前記制御対象が零点である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点から遠ざけるものである。
すなわち、汎用コンピュータ1は、プログラムを実行するCPU11と、CPU11が実行するプログラムを記憶するプログラムメモリ12と、CPU11の処理するデータを記憶するデータメモリ13と、周辺機器接続用のインターフェイス14と、CPU11、プログラムメモリ12、データメモリ13およびインターフェイス14を相互に接続するバス15と、から構成される。
すなわち、汎用コンピュータ1は、例えばCD−ROMである所定のメディアに記録された伝達関数調整プログラムを、ディスクドライブ16を使用して読み込み、プログラムメモリ12に記憶する。
以下、伝達関数調整プログラムのフローチャートを参照しつつ、伝達関数調整装置の動作を説明する。
この入力形式は特に限定されないが、例えば[数1]で表されるパルス伝達関数A(z)の次数(n,m)ならびに係数(a0,a1,・・・an−1,an)および(b0,b1,・・・bm−1,bm)を読み込む方法を適用できる。
なお、周波数応答を直接読み込むようにしても、時間領域におけるインパルス応答から直接パルス伝達関数A(z)を同定するようにしてもよい。
図3の場合は、4つのパルス伝達関数A1(z)、A2(z)、A3(z)およびA4(z)のそれぞれの周波数振幅の平均値A1av、A2av、A3avおよびA4avが算出される。
第1のパルス伝達関数A1(z)を例として、ピークおよびディップの抽出方法を説明する。
他のパルス伝達関数A2(z)、A3(z)およびA4(z)についても、同様にピークおよびディップの抽出を行う。
他のパルス伝達関数A2(z)、A3(z)およびA4(z)についても同様の判定を行い、いずれかのパルス伝達関数にピーク閾値Pthを超えるピークもしくはディップ閾値Dthを超えるディップがある場合はステップS26に進む。
CPU11は、ピーク閾値Pthを超えるピークならびにディップ閾値Dthを超えるディップの中から各パルス伝達関数の周波数振幅の平均値からの乖離が最大であるピークあるいはディップを補正対象に決定する(ステップS26)。
図4は、CPU11がステップS27で実行する制御対象決定処理のフローチャートである。
まずCPU11は、調整対象のパルス伝達関数が1つであるか否かを判定する(ステップS271)。
CPU11は、調整対象は1つでないと判定した場合には、全パルス伝達関数が極あるいは零点を共有しているか否かを判定する(ステップS273)。
[数2]で表される逆フィルタの例では、4つのパルス伝達関数A1(z)、A2(z)、A3(z)およびA4(z)は共通の分母Ad(z)を有しているので、Ad(z)=0の根を極として共有していることとなる。
CPU11は、全パルス伝達関数が極もしくは零点を共有していると判定した場合は、補正対象伝達関数であるパルス伝達関数以外の他のパルス伝達関数(補正対象外伝達関数)が、補正対象に対応する周波数(補正対象周波数fo)において補正対象がピークである場合はピークを、補正対象がディップである場合はディップを有するか否かを判定する(ステップS274)。
図3の場合には、第1のパルス伝達関数A1(z)のピークP12に対応する周波数f13において、他のパルス伝達関数A2(z)、A3(z)およびA4(z)も平均値A2av、A3avおよびA4avから予め定めた副ピーク閾値pth乖離したピークを有するか否かを判定する。
[数2]で表される逆フィルタの例では、探索範囲は、4つのパルス伝達関数A1(z)、A2(z)、A3(z)およびA4(z)に共通する極、すなわちAd(z)=0の根となる。
[数2]で表される逆フィルタの例では、探索範囲は、4つのパルス伝達関数A1(z)、A2(z)、A3(z)およびA4(z)に共通でない零点、例えばA1'(z)=0の根となる。
CPU11は、補正対象周波数foの正規化補正対象周波数Foを算出する。
図3の場合には、第1のパルス伝達関数A1(z)のピークP12に対応する周波数f13が補正対象周波数foに決定される。
図5は、制御対象の決定方法を説明するためのグラフであって、円41はZ平面の原点Oを中心とする単位円、点Qは正規化補正対象周波数Foを偏角とする単位円上の点である。
そして、点Rの偏角をωR、点Sの偏角をωSとすると、図5の場合には[数3]が成り立つ。
なお、制御対象の決定に際し、正規化補正対象周波数Foに対応する単位円上の点Qに対してZ平面の実軸(横軸)を挟んで位置する極または零点は偏角の差が詐称であるとしても制御対象としては決定しないものとする。
1−1:制御対象を、単位円を越えて移動させない。
制御対象は単位円の内部にあることが普通であるが、その場合は制御対象を単位円の外側に移動させない。これは伝達関数が不安定となることを防止するためである。
また、制御対象が単位円の外側にある場合であっても、制御対象を単位円の内側に移動させない。これは時間領域の特性が大きく変化することを防止するためである。
1−2:制御対象を、Z平面の実軸(横軸)を越えて移動させない。
すなわち、制御対象が第1または第2象限に存在する場合は、制御対象を第4または第3象限に移動させない。これは極および零点は実軸に対して線対称に配置されているので、実軸を越えて移動させる必要は通常ありえないからである。
なお、Z平面の虚軸(縦軸)を超える移動、すなわち、第1象限と第2象限間の移動、第3象限と第4象限間の移動は許容する。
2−1:補正対象がピークであって、制御対象が極Rであるときは、極Rを点Qから遠ざけるように移動する。
2−2:補正対象がピークであって、制御対象が零点Sであるときは、零点Sを点Qに近付けるように移動する。
2−3:補正対象がディップであって、制御対象が極Rであるときは、極Rを点Qに近付けるように移動する。
2−4:補正対象がディップであって、制御対象が零点Sであるときは、極Sを点Qから遠ざけるように移動する。
規則2においては制御対象を、点Qに近ける、または点Qから遠ざけるように移動させれば、その移動方向は特に規定されない。
なお、規則2に代えて規則3を適用してもよい。
3−1:補正対象がピークであって、制御対象が極Rであるときは、極Rを直線RQ上で点Qから遠ざけるように移動する。
3−2:補正対象がピークであって、制御対象が零点Sであるときは、零点Sを直線SQ上で点Qに近付けるように移動する。
3−3:補正対象がディップであって、制御対象が極Rであるときは、極Rを直線RQ上で点Qに近付けるように移動する。
3−4:補正対象がディップであって、制御対象が零点Sであるときは、零点Sを直線SQ上で点Qから遠ざけるように移動する。
ただし、補正対象がピーク閾値Pth以下となる極Rの移動量を解析的に求めることはできない。
そこで、実際には、極Rを所定の移動量ΔLだけ移動させて周波数振幅特性を計算し、ピークP12がピーク閾値Pth以下となったか否かを判定し、ピークP12がピーク閾値Pth以下となっていない場合には伝達関数調整プログラムを繰り返し実行するという再帰的な手法を適用する。
すなわち、本発明に係る伝達関数調整装置によれば、調整による周波数振幅特性および周波数位相特性への影響を局所的に限定し、かつ、安定性を維持しつつ、パルス伝達関数の周波数振幅特性のピークおよびディップを抑制できる。
11...CPU
12...プログラムメモリ
13...データメモリ
14...インターフェイス
15...バス
16...ディスクドライブ
17...キーボード
18...ポインティングデバイス
19...ディスプレイ
Claims (7)
- 少なくとも1つの伝達関数の周波数振幅特性を算出する周波数振幅特性算出手段と、
前記伝達関数ごとに、前記周波数振幅の平均値を算出する周波数振幅平均値算出手段と、
前記伝達関数ごとに、予め定められたピーク閾値を超える前記周波数振幅の平均値を基準としたピークおよび予め定められたディップ閾値を超える前記周波数振幅の平均値を基準としたディップの少なくとも一方を、前記周波数振幅特性から抽出するピーク・ディップ抽出手段と、
すべての前記伝達関数の中で、前記ピークおよび前記ディップのうち前記周波数振幅の平均値からの乖離が最大である1つのピークまたはディップを有する伝達関数を補正対象伝達関数として決定するとともに、当該ピークまたはディップを補正対象として決定する補正対象決定手段と、
前記補正対象に対応した周波数の正規化周波数との差が最小である偏角を有する前記少なくとも1つの伝達関数の極および零点の一方を制御対象に決定する制御対象決定手段と、
前記制御対象を前記補正対象がピーク閾値またはディップ閾値内に収まる最終位置まで予め定めた規則に従って移動させる制御対象移動手段と、を備える伝達関数調整装置。 - 前記制御対象決定手段が、前記補正対象伝達関数と前記補正対象伝達関数以外の伝達関数である補正対象外伝達関数とが極または零点を共有し、前記補正対象伝達関数と前記補正対象外伝達関数とが前記補正対象に対応する周波数においてピークまたはディップを有するときは、前記補正対象伝達関数および前記補正対象外伝達関数の極および零点の中から制御対象を決定するものである請求項1に記載の伝達関数調整装置。
- 前記制御対象移動手段が、前記制御対象を、前記単位円の円周を横切らず、かつ、前記Z平面の実軸を横切らない範囲で移動させるものである請求項1または請求項2に記載の伝達関数調整装置。
- 前記制御対象移動手段が、
前記補正対象がピークであり、前記制御対象が極である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点から遠ざけ、
前記補正対象がピークであり、前記制御対象が零点である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点に近付け、
前記補正対象がディップであり、前記制御対象が極である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点に近付け、
前記補正対象がディップであり、前記制御対象が零点である場合は、前記制御対象を前記補正対象の正規化周波数に対応する前記単位円上の点から遠ざけるものである請求項3に記載の伝達関数調整装置。 - 前記制御対象移動手段が、
前記制御対象を、前記制御対象と前記補正対象の正規化周波数に対応する前記単位円上の点とを結ぶ直線上で移動させるものである請求項4に記載の伝達関数調整装置。 - コンピュータで、
少なくとも1つの伝達関数の周波数振幅特性を算出する周波数振幅特性算出処理と、
前記伝達関数ごとに、前記周波数振幅の平均値を算出する周波数振幅平均値算出処理と、
前記伝達関数ごとに、予め定められたピーク閾値を超える前記周波数振幅の平均値を基準としたピークおよび予め定められたディップ閾値を超える前記周波数振幅の平均値を基準としたディップの少なくとも一方を、前記周波数振幅特性から抽出するピーク・ディップ抽出処理と、
すべての前記伝達関数の中から、前記ピークおよび前記ディップの中から前記周波数振幅の平均値からの乖離が最大である1つのピークまたはディップを有する伝達関数を補正対象伝達関数として決定するとともに、当該ピークまたはディップを補正対象として決定する補正対象決定処理と、
前記補正対象に対応した周波数の正規化周波数との差が最小である偏角を有する前記少なくとも1つの伝達関数の極および零点の一方を制御対象に決定する制御対象決定処理と、
前記制御対象を前記補正対象がピーク閾値またはディップ閾値内に収まる最終位置まで予め定めた規則に従って移動させる制御対象移動処理と、を行う伝達関数調整プログラム。 - 少なくとも1つの伝達関数の周波数振幅特性を算出する周波数振幅特性算出段階と、
前記伝達関数ごとに、前記周波数振幅の平均値を算出する周波数振幅平均値算出段階と、
前記伝達関数ごとに、予め定められたピーク閾値を超える前記周波数振幅の平均値を基準としたピークおよび予め定められたディップ閾値を超える前記周波数振幅の平均値を基準としたディップの少なくとも一方を、前記周波数振幅特性から抽出するピーク・ディップ抽出段階と、
すべての前記伝達関数の中から、前記ピークおよび前記ディップの中から前記周波数振幅の平均値からの乖離が最大である1つのピークまたはディップを有する伝達関数を補正対象伝達関数として決定するとともに、当該ピークまたはディップを補正対象として決定する補正対象決定段階と、
前記補正対象に対応した周波数の正規化周波数との差が最小である偏角を有する前記少なくとも1つの伝達関数の極および零点の一方を制御対象に決定する制御対象決定段階と、
前記制御対象を前記補正対象がピーク閾値またはディップ閾値内に収まる最終位置まで予め定めた規則に従って移動させる制御対象移動段階と、を含む伝達関数調整方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010199773A JP5514050B2 (ja) | 2010-09-07 | 2010-09-07 | 伝達関数調整装置、伝達関数調整プログラムおよび伝達関数調整方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010199773A JP5514050B2 (ja) | 2010-09-07 | 2010-09-07 | 伝達関数調整装置、伝達関数調整プログラムおよび伝達関数調整方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012060295A true JP2012060295A (ja) | 2012-03-22 |
JP5514050B2 JP5514050B2 (ja) | 2014-06-04 |
Family
ID=46056887
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010199773A Active JP5514050B2 (ja) | 2010-09-07 | 2010-09-07 | 伝達関数調整装置、伝達関数調整プログラムおよび伝達関数調整方法 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5514050B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016052082A (ja) * | 2014-09-02 | 2016-04-11 | 日本放送協会 | 制御器設計装置、制御器設計方法及びプログラム |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2755081B2 (ja) * | 1992-11-30 | 1998-05-20 | 日本ビクター株式会社 | 音像定位制御方法 |
JPH10509565A (ja) * | 1994-08-25 | 1998-09-14 | アダプティブ オーディオ リミテッド | 録音及び再生システム |
JPH11136790A (ja) * | 1997-10-29 | 1999-05-21 | Asahi Kasei Micro Syst Co Ltd | 帯域強調回路 |
JP2003248703A (ja) * | 2002-02-22 | 2003-09-05 | Univ Nihon | フィルタの設計方法及びプログラム |
JP2006071775A (ja) * | 2004-08-31 | 2006-03-16 | Aiphone Co Ltd | 音声再生方法およびその装置 |
JP2006279863A (ja) * | 2005-03-30 | 2006-10-12 | Clarion Co Ltd | 頭部伝達関数の補正方法 |
-
2010
- 2010-09-07 JP JP2010199773A patent/JP5514050B2/ja active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2755081B2 (ja) * | 1992-11-30 | 1998-05-20 | 日本ビクター株式会社 | 音像定位制御方法 |
JPH10509565A (ja) * | 1994-08-25 | 1998-09-14 | アダプティブ オーディオ リミテッド | 録音及び再生システム |
JPH11136790A (ja) * | 1997-10-29 | 1999-05-21 | Asahi Kasei Micro Syst Co Ltd | 帯域強調回路 |
JP2003248703A (ja) * | 2002-02-22 | 2003-09-05 | Univ Nihon | フィルタの設計方法及びプログラム |
JP2006071775A (ja) * | 2004-08-31 | 2006-03-16 | Aiphone Co Ltd | 音声再生方法およびその装置 |
JP2006279863A (ja) * | 2005-03-30 | 2006-10-12 | Clarion Co Ltd | 頭部伝達関数の補正方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016052082A (ja) * | 2014-09-02 | 2016-04-11 | 日本放送協会 | 制御器設計装置、制御器設計方法及びプログラム |
Also Published As
Publication number | Publication date |
---|---|
JP5514050B2 (ja) | 2014-06-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107925816B (zh) | 用于在波束成形的音频中重新创建方向提示的方法和装置 | |
CN110024419B (zh) | 增益相位均衡(gpeq)滤波器和调谐方法 | |
JP6337519B2 (ja) | 音声処理装置、雑音抑圧方法、およびプログラム | |
JP2011061422A (ja) | 情報処理装置、情報処理方法およびプログラム | |
EP3441965A1 (en) | Signal processing device, signal processing method, and program | |
US9559658B2 (en) | Filter coefficient group computation device and filter coefficient group computation method | |
JP2019047478A (ja) | 音響信号処理装置、音響信号処理方法および音響信号処理プログラム | |
US10008997B2 (en) | Filter generator, filter generation method, and filter generation program | |
KR101944758B1 (ko) | 스테레오 신호의 스테레오 이미지를 수정하기 위한 오디오 신호 처리 장치 및 방법 | |
JP5514050B2 (ja) | 伝達関数調整装置、伝達関数調整プログラムおよび伝達関数調整方法 | |
JP6171926B2 (ja) | 頭外音像定位装置、頭外音像定位方法、及び、プログラム | |
CN110858487A (zh) | 一种音频信号缩放处理方法及装置 | |
US8135137B2 (en) | Sound image localization apparatus | |
JP4413043B2 (ja) | 周期性ノイズ抑圧方法、周期性ノイズ抑圧装置、周期性ノイズ抑圧プログラム | |
EP3444946B1 (en) | Acoustic processing device and acoustic processing method | |
US10362396B2 (en) | Phase control signal generation device, phase control signal generation method, and phase control signal generation program | |
US8306242B2 (en) | Heyser spiral low frequency correction of FIR filters | |
JP2009065452A (ja) | 音像定位制御装置、音像定位制御方法、プログラム、および集積回路 | |
CN114067817A (zh) | 低音增强方法、装置、电子设备及存储介质 | |
WO2021251136A1 (ja) | 信号処理装置、信号処理方法、信号処理プログラム、信号処理用モデル製造方法及び音響出力機器 | |
JP2017212732A (ja) | チャンネル数変換装置およびプログラム | |
JP2017163448A (ja) | 音調整装置及び音調整方法 | |
JP6603725B2 (ja) | 音声信号生成装置、音声信号生成方法、及び、プログラム | |
JP2007101782A (ja) | 音響調整装置、音響調整プログラム、および音響調整方法 | |
KR20120022650A (ko) | 출력음향신호에 따라 필터의 계수를 적응적으로 조정하는 음향 신호 재생 방법 및 장치 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130108 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131024 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131029 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131220 |
|
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: 20140304 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140328 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5514050 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |