JP2001086585A - スピーカの二次歪み除去装置 - Google Patents
スピーカの二次歪み除去装置Info
- Publication number
- JP2001086585A JP2001086585A JP25547999A JP25547999A JP2001086585A JP 2001086585 A JP2001086585 A JP 2001086585A JP 25547999 A JP25547999 A JP 25547999A JP 25547999 A JP25547999 A JP 25547999A JP 2001086585 A JP2001086585 A JP 2001086585A
- Authority
- JP
- Japan
- Prior art keywords
- convolution
- processing
- primary
- speaker
- eigenvectors
- 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.)
- Withdrawn
Links
Landscapes
- Circuit For Audible Band Transducer (AREA)
Abstract
(57)【要約】
【課題】 畳み込み処理の演算量を低減する。
【解決手段】 固有値分解処理部５０により歪み除去対
象のスピーカ３１の二次ボルテラ時間応答５１を固有ベ
クトルｕiと固有値λiに分解して二次畳み込み処理にと
って寄与が大きいもののみを選択して記憶部５４に格納
する。Ａ／Ｄ変換器１１により変換されたデジタルオー
ディオ信号ｘと記憶部５４に格納されていたＭ組の固有
ベクトルｕiがそれぞれ一次畳み込み器５５−１〜５５
−Ｍにより乗算され、次いで自乗器５６−１〜５６−Ｍ
によりそれぞれ一次畳み込み器５５−１〜５５−Ｍの乗
算結果が自乗され、次いで係数乗算器５７−１〜５７−
Ｍによりそれぞれ自乗器５６−１〜５６−Ｍの各乗算結
果と記憶部５４からの固有値λiが乗算される。
象のスピーカ３１の二次ボルテラ時間応答５１を固有ベ
クトルｕiと固有値λiに分解して二次畳み込み処理にと
って寄与が大きいもののみを選択して記憶部５４に格納
する。Ａ／Ｄ変換器１１により変換されたデジタルオー
ディオ信号ｘと記憶部５４に格納されていたＭ組の固有
ベクトルｕiがそれぞれ一次畳み込み器５５−１〜５５
−Ｍにより乗算され、次いで自乗器５６−１〜５６−Ｍ
によりそれぞれ一次畳み込み器５５−１〜５５−Ｍの乗
算結果が自乗され、次いで係数乗算器５７−１〜５７−
Ｍによりそれぞれ自乗器５６−１〜５６−Ｍの各乗算結
果と記憶部５４からの固有値λiが乗算される。
Description
【０００１】
【発明の属する技術分野】本発明は、オーディオ信号を
再生するスピーカの二次歪み除去装置に関する。
再生するスピーカの二次歪み除去装置に関する。
【０００２】
【従来の技術】スピーカの高調波歪みや混変調歪みを除
去する方法として、図４に示すようにボルテラフィルタ
と呼ばれる二次畳み込み演算を利用したシステムが提案
されている。なお、この種の従来例としては、例えば特
開平８−３１７４９０号公報、特開平９−１３９９９５
号公報、特開平１０−５１８８１号公報、特開平１１−
５５７８２号公報などに示されているものがある。
去する方法として、図４に示すようにボルテラフィルタ
と呼ばれる二次畳み込み演算を利用したシステムが提案
されている。なお、この種の従来例としては、例えば特
開平８−３１７４９０号公報、特開平９−１３９９９５
号公報、特開平１０−５１８８１号公報、特開平１１−
５５７８２号公報などに示されているものがある。
【０００３】図４に示すようにこのシステムでは、ま
ず、歪み除去対象のスピーカ３１を一次伝達関数Ｈ1で
表現される一次システム３７と、二次伝達関数Ｈ2で表
現される二次システム４１の和（加算器３８）と考え
る。ここで、一次システム３７では次式（１）
ず、歪み除去対象のスピーカ３１を一次伝達関数Ｈ1で
表現される一次システム３７と、二次伝達関数Ｈ2で表
現される二次システム４１の和（加算器３８）と考え
る。ここで、一次システム３７では次式（１）
【０００４】
【数１】
【０００５】の畳み込み演算が行われ、二次システム４
１では次式（２）
１では次式（２）
【０００６】
【数２】
【０００７】の二次畳み込み演算が行われる。ここで、
式（１）におけるｈ1(ｋ1)と式（２）におけるｈ
2（ｋ1，ｋ2）は、それぞれ一次伝達関数Ｈ1、二次伝達
関数Ｈ2の時間応答である。これを図５、図６に示す。
ただしｋ1、ｋ2は離散化された時間軸上における時刻を
表す整数値であり、ｈ1（ｋ1）はＮ、ｈ2（ｋ1，ｋ2）
はＮ×Ｎの応答長を持つものとする。
式（１）におけるｈ1(ｋ1)と式（２）におけるｈ
2（ｋ1，ｋ2）は、それぞれ一次伝達関数Ｈ1、二次伝達
関数Ｈ2の時間応答である。これを図５、図６に示す。
ただしｋ1、ｋ2は離散化された時間軸上における時刻を
表す整数値であり、ｈ1（ｋ1）はＮ、ｈ2（ｋ1，ｋ2）
はＮ×Ｎの応答長を持つものとする。
【０００８】そして、スピーカ３１において発生する二
次成分を除去するためにスピーカ３１の前段に二次歪み
除去装置３０が設けられる。二次歪み除去装置３０では
アナログオーディオ信号をＡ／Ｄ変換器１１によりデジ
タル信号に変換して、Ａ／Ｄ変換器１１の出力ｘを遅延
部１２により二次畳み込み処理部１７及び一次畳み込み
処理部１８の処理時間だけ遅延して加算器１３に出力す
るとともに、Ａ／Ｄ変換器１１の出力ｘを二次畳み込み
処理部１７により二次畳み込み処理して二次畳み込み処
理部１７の出力＾ｙ2を一次畳み込み処理部１８により
一次畳み込み処理して加算器１３に出力する（なお
「＾」の記号は、図中に示されているように、本来その
後の文字の上部に位置するが、便宜上分けて表記してい
る）。そして、加算器１３により遅延部１２の出力ｘ’
から一次畳み込み処理部１８の出力＾ｙ2’を減算し、
次いで加算器１３の出力ｘ’−＾ｙ2’をＤ／Ａ変換器
１４によりアナログ信号に変換してスピーカ３１に出力
する。
次成分を除去するためにスピーカ３１の前段に二次歪み
除去装置３０が設けられる。二次歪み除去装置３０では
アナログオーディオ信号をＡ／Ｄ変換器１１によりデジ
タル信号に変換して、Ａ／Ｄ変換器１１の出力ｘを遅延
部１２により二次畳み込み処理部１７及び一次畳み込み
処理部１８の処理時間だけ遅延して加算器１３に出力す
るとともに、Ａ／Ｄ変換器１１の出力ｘを二次畳み込み
処理部１７により二次畳み込み処理して二次畳み込み処
理部１７の出力＾ｙ2を一次畳み込み処理部１８により
一次畳み込み処理して加算器１３に出力する（なお
「＾」の記号は、図中に示されているように、本来その
後の文字の上部に位置するが、便宜上分けて表記してい
る）。そして、加算器１３により遅延部１２の出力ｘ’
から一次畳み込み処理部１８の出力＾ｙ2’を減算し、
次いで加算器１３の出力ｘ’−＾ｙ2’をＤ／Ａ変換器
１４によりアナログ信号に変換してスピーカ３１に出力
する。
【０００９】二次畳み込み処理部１７の伝達係数は、遅
延部１２→加算器１３→Ｄ／Ａ変換器１４→二次システ
ム４１の二次経路「１」を通過する信号と、二次畳み込
み処理部１７→一次畳み込み処理部１８→加算器１３→
Ｄ／Ａ変換器１４→一次システム３７の二次経路「２」
を通過する二次成分が相殺されるように設定され、この
例ではスピーカ３１の二次伝達関数Ｈ2の推定値＾Ｈ2に
設定される。二次畳み込み処理部１７では次式（３）
延部１２→加算器１３→Ｄ／Ａ変換器１４→二次システ
ム４１の二次経路「１」を通過する信号と、二次畳み込
み処理部１７→一次畳み込み処理部１８→加算器１３→
Ｄ／Ａ変換器１４→一次システム３７の二次経路「２」
を通過する二次成分が相殺されるように設定され、この
例ではスピーカ３１の二次伝達関数Ｈ2の推定値＾Ｈ2に
設定される。二次畳み込み処理部１７では次式（３）
【００１０】
【数３】
【００１１】の二次畳み込み演算が行われる。ここで、
式（３）における＾ｈ2（ｋ1，ｋ2）は推定値＾Ｈ2に対
応する二次時間応答である。
式（３）における＾ｈ2（ｋ1，ｋ2）は推定値＾Ｈ2に対
応する二次時間応答である。
【００１２】ところで、二次畳み込み演算式（３）に示
される処理は膨大な演算を要するので、十分なタップ長
をＤＳＰ（デジタル・シグナル・プロセッサ）によりリ
アルタイムで畳み込み処理することは困難である。そこ
で、例えば特開平１０−５１８８１号公報、特開平１１
−５５７８２号公報に示されるように、式（３）をフー
リエ変換して次式（４）
される処理は膨大な演算を要するので、十分なタップ長
をＤＳＰ（デジタル・シグナル・プロセッサ）によりリ
アルタイムで畳み込み処理することは困難である。そこ
で、例えば特開平１０−５１８８１号公報、特開平１１
−５５７８２号公報に示されるように、式（３）をフー
リエ変換して次式（４）
【００１３】
【数４】
【００１４】で表される処理を実施する手法が提案され
ている。式（４）における二次畳み込み演算では、式
（３）のように時間領域で行っていた積和演算が周波数
領域における乗算になるので、演算量を大幅に低減する
ことができる。また、畳み込む周波数を選択することが
できるので、不要な周波数帯域を畳み込まないように畳
み込む周波数を選択すれば、更に演算量を大幅に低減す
ることができる。
ている。式（４）における二次畳み込み演算では、式
（３）のように時間領域で行っていた積和演算が周波数
領域における乗算になるので、演算量を大幅に低減する
ことができる。また、畳み込む周波数を選択することが
できるので、不要な周波数帯域を畳み込まないように畳
み込む周波数を選択すれば、更に演算量を大幅に低減す
ることができる。
【００１５】
【発明が解決しようとする課題】しかしながら、式
（４）における周波数領域における処理の場合にも演算
量は依然多いという問題点がある。また、フーリエ変
換、逆フーリエ変換を行うのでオーバヘッドが大きく、
このため畳み込みの規模をある程度大きくしないと、演
算量低減の効果は少ない。
（４）における周波数領域における処理の場合にも演算
量は依然多いという問題点がある。また、フーリエ変
換、逆フーリエ変換を行うのでオーバヘッドが大きく、
このため畳み込みの規模をある程度大きくしないと、演
算量低減の効果は少ない。
【００１６】本発明は上記従来例の問題点に鑑み、畳み
込み処理の演算量を低減することができるスピーカの二
次歪み除去装置を提供することを目的とする。
込み処理の演算量を低減することができるスピーカの二
次歪み除去装置を提供することを目的とする。
【００１７】
【課題を解決するための手段】本発明は上記目的を達成
するために、歪み除去対象のスピーカの二次ボルテラ時
間応答を固有ベクトル及び固有値に分解し、この内、歪
み除去に寄与が大きいものを選択して一次畳み込み処理
を行うようにしたものである。
するために、歪み除去対象のスピーカの二次ボルテラ時
間応答を固有ベクトル及び固有値に分解し、この内、歪
み除去に寄与が大きいものを選択して一次畳み込み処理
を行うようにしたものである。
【００１８】すなわち本発明によれば、歪み除去対象の
スピーカの二次ボルテラ時間応答を固有ベクトル及び固
有値に分解し、前記固有ベクトル及び固有値の内、歪み
除去に寄与が大きいものを選択する固有値分解手段と、
アナログオーディオ信号をデジタル信号に変換するＡ／
Ｄ変換手段と、前記Ａ／Ｄ変換手段により変換されたデ
ジタルオーディオ信号と前記固有値分解手段により分解
及び選択された複数の固有ベクトルをそれぞれ一次畳み
込み処理する複数の第１の一次畳み込み手段と、前記複
数の第１の一次畳み込み手段の処理結果をそれぞれ自乗
する複数の自乗手段と、前記複数の自乗手段の処理結果
と前記固有値分解手段により分解及び選択された複数の
固有値をそれぞれ乗算する複数の乗算手段と、前記複数
の乗算手段の各処理結果を加算する第１の加算手段と、
前記第１の加算手段の処理結果を一次畳み込み処理する
第２の一次畳み込み手段と、前記Ａ／Ｄ変換手段により
変換されたデジタルオーディオ信号から前記第２の一次
畳み込み手段の処理結果を差し引く第２の加算手段と、
前記第２の加算手段の処理結果をアナログオーディオ信
号に変換するＤ／Ａ変換手段とを備え、前記Ｄ／Ａ変換
手段の出力を前記歪み除去対象のスピーカに印加するよ
うにしたスピーカの二次歪み除去装置が提供される。
スピーカの二次ボルテラ時間応答を固有ベクトル及び固
有値に分解し、前記固有ベクトル及び固有値の内、歪み
除去に寄与が大きいものを選択する固有値分解手段と、
アナログオーディオ信号をデジタル信号に変換するＡ／
Ｄ変換手段と、前記Ａ／Ｄ変換手段により変換されたデ
ジタルオーディオ信号と前記固有値分解手段により分解
及び選択された複数の固有ベクトルをそれぞれ一次畳み
込み処理する複数の第１の一次畳み込み手段と、前記複
数の第１の一次畳み込み手段の処理結果をそれぞれ自乗
する複数の自乗手段と、前記複数の自乗手段の処理結果
と前記固有値分解手段により分解及び選択された複数の
固有値をそれぞれ乗算する複数の乗算手段と、前記複数
の乗算手段の各処理結果を加算する第１の加算手段と、
前記第１の加算手段の処理結果を一次畳み込み処理する
第２の一次畳み込み手段と、前記Ａ／Ｄ変換手段により
変換されたデジタルオーディオ信号から前記第２の一次
畳み込み手段の処理結果を差し引く第２の加算手段と、
前記第２の加算手段の処理結果をアナログオーディオ信
号に変換するＤ／Ａ変換手段とを備え、前記Ｄ／Ａ変換
手段の出力を前記歪み除去対象のスピーカに印加するよ
うにしたスピーカの二次歪み除去装置が提供される。
【００１９】
【発明の実施の形態】以下、図面を参照して本発明の実
施の形態を説明する。図１は本発明に係るスピーカの二
次歪み除去装置の一実施形態を示すブロック図、図２は
図１の一次畳み込み器の構成を詳しく示すブロック図、
図３は畳み込み係数の固有値を示す説明図である。
施の形態を説明する。図１は本発明に係るスピーカの二
次歪み除去装置の一実施形態を示すブロック図、図２は
図１の一次畳み込み器の構成を詳しく示すブロック図、
図３は畳み込み係数の固有値を示す説明図である。
【００２０】二次畳み込み演算式（３）を変形して行列
で表すと、次式（５）
で表すと、次式（５）
【００２１】
【数５】
【００２２】のような二次形式で表現することができ
る。ここで、式（５）における畳み込み係数行列＾Ｈ2
は対称行列であるので、畳み込み係数の固有値λi、固
有ベクトルｕiは実数であり、また、固有ベクトルｕiは
お互いに直交する。また、畳み込み係数行列＾Ｈ2は、
固有ベクトルｕiを並べた正規行列Ｕと、固有値λiを対
角線上に並べた対角行列Λにより次式（６）
る。ここで、式（５）における畳み込み係数行列＾Ｈ2
は対称行列であるので、畳み込み係数の固有値λi、固
有ベクトルｕiは実数であり、また、固有ベクトルｕiは
お互いに直交する。また、畳み込み係数行列＾Ｈ2は、
固有ベクトルｕiを並べた正規行列Ｕと、固有値λiを対
角線上に並べた対角行列Λにより次式（６）
【００２３】
【数６】
【００２４】のように対角化される。そして、式（５）
は次式（７）
は次式（７）
【００２５】
【数７】
【００２６】のように展開されるので、二次畳み込みは
固有ベクトルｕiで構成されるモードの重ね合わせと解
釈することができる。そして、式（７）の各項λi（ｕi
Tｘ）2におけるｕi Tｘ（Ｔは転置を表す）は、固有ベク
トルｕiを係数とする一次畳み込みであるので、式
（７）を処理する回路ブロックは、図２に示すようにオ
ーディオ信号ｘを遅延する複数段の遅延器１０１
（Ｚ-1）と、各遅延器１０１の出力ｘに対して固有ベク
トルｕiを乗算する複数段の係数乗算器１０２と、この
乗算結果ｘｕiを加算する加算器１０３と、加算器１０
３の出力を自乗する自乗器５６と、自乗器５６の出力に
固有値λiを乗算する係数乗算器５７により実現するこ
とができる。遅延器１０１と、係数乗算器１０２と加算
器１０３は、図１に示す一次畳み込み演算器５５を構成
する。
固有ベクトルｕiで構成されるモードの重ね合わせと解
釈することができる。そして、式（７）の各項λi（ｕi
Tｘ）2におけるｕi Tｘ（Ｔは転置を表す）は、固有ベク
トルｕiを係数とする一次畳み込みであるので、式
（７）を処理する回路ブロックは、図２に示すようにオ
ーディオ信号ｘを遅延する複数段の遅延器１０１
（Ｚ-1）と、各遅延器１０１の出力ｘに対して固有ベク
トルｕiを乗算する複数段の係数乗算器１０２と、この
乗算結果ｘｕiを加算する加算器１０３と、加算器１０
３の出力を自乗する自乗器５６と、自乗器５６の出力に
固有値λiを乗算する係数乗算器５７により実現するこ
とができる。遅延器１０１と、係数乗算器１０２と加算
器１０３は、図１に示す一次畳み込み演算器５５を構成
する。
【００２７】さらに、Ｎ×Ｎの二次畳み込み係数を固有
値分解すると、Ｎ組の固有ベクトルｕiとＮ個の固有値
λiが得られるが、固有値λiは各モードを重ね合わせる
際の重みと考えることができるので、式（７）において λ1＞λ2・・・＞λM・・・＞λN として、λMが十分小さく第Ｍ項以降は演算を打ち切る
ことができる場合には式（７）は次式（８）
値分解すると、Ｎ組の固有ベクトルｕiとＮ個の固有値
λiが得られるが、固有値λiは各モードを重ね合わせる
際の重みと考えることができるので、式（７）において λ1＞λ2・・・＞λM・・・＞λN として、λMが十分小さく第Ｍ項以降は演算を打ち切る
ことができる場合には式（７）は次式（８）
【００２８】
【数８】
【００２９】で表すことができる。ここで、図３は実際
のスピーカの二次畳み込み係数を固有値分解したときの
固有値の一例を示し、固有値λiの大半は畳み込みに寄
与していないことがわかる。したがって、ノイズ除去に
貢献しない固有ベクトルｕiと固有値λiを除いて一次畳
み込みすれば、大幅に演算量を低減することができる。
のスピーカの二次畳み込み係数を固有値分解したときの
固有値の一例を示し、固有値λiの大半は畳み込みに寄
与していないことがわかる。したがって、ノイズ除去に
貢献しない固有ベクトルｕiと固有値λiを除いて一次畳
み込みすれば、大幅に演算量を低減することができる。
【００３０】次に図１を参照して本発明の二次歪み除去
装置について詳しく説明する。図１に示す装置は、図４
に示す二次畳み込み処理部１７の代わりに、固有値分解
処理部５０と二次畳み込み処理部１７ａが設けられ、他
の構成（Ａ／Ｄ変換器１１、遅延部１２、加算器１３、
Ｄ／Ａ変換器１４、一次畳み込み処理部１８）は図１と
同じである。まず、実際の歪み除去時（オーディオ再生
時）の前に、固有値分解処理部５０では、歪み除去対象
のスピーカ３１（図４参照）の二次ボルテラ時間応答５
１が固有値分解部５２により固有ベクトルｕiと固有値
λiに分離される。
装置について詳しく説明する。図１に示す装置は、図４
に示す二次畳み込み処理部１７の代わりに、固有値分解
処理部５０と二次畳み込み処理部１７ａが設けられ、他
の構成（Ａ／Ｄ変換器１１、遅延部１２、加算器１３、
Ｄ／Ａ変換器１４、一次畳み込み処理部１８）は図１と
同じである。まず、実際の歪み除去時（オーディオ再生
時）の前に、固有値分解処理部５０では、歪み除去対象
のスピーカ３１（図４参照）の二次ボルテラ時間応答５
１が固有値分解部５２により固有ベクトルｕiと固有値
λiに分離される。
【００３１】そして、この固有ベクトルｕiと固有値λi
の内、前述したように二次畳み込み処理にとって寄与が
大きいもののみが選択部５３により選択されて記憶部５
４に格納される。例えば二次ボルテラ時間応答５１の規
模がＮ×Ｎの場合、固有値分解部５２ではＮ組の固有ベ
クトルｕiとＮ個の固有値λiが算出されるが、選択部５
３では固有値λiの絶対値に基づいて二次畳み込み処理
にとって寄与が大きいＭ（＜Ｎ）組の固有ベクトルｕi
とＭ個の固有値λiが選択される。
の内、前述したように二次畳み込み処理にとって寄与が
大きいもののみが選択部５３により選択されて記憶部５
４に格納される。例えば二次ボルテラ時間応答５１の規
模がＮ×Ｎの場合、固有値分解部５２ではＮ組の固有ベ
クトルｕiとＮ個の固有値λiが算出されるが、選択部５
３では固有値λiの絶対値に基づいて二次畳み込み処理
にとって寄与が大きいＭ（＜Ｎ）組の固有ベクトルｕi
とＭ個の固有値λiが選択される。
【００３２】実際の歪み除去時（オーディオ再生時）に
は、Ａ／Ｄ変換器１１により変換されたデジタルオーデ
ィオ信号ｘと記憶部５４に格納されていたＭ組の固有ベ
クトルｕiがそれぞれ一次畳み込み器５５−１〜５５−
Ｍにより乗算され、次いで自乗器５６−１〜５６−Ｍに
より、それぞれ一次畳み込み器５５−１〜５５−Ｍの乗
算結果が自乗され、次いで係数乗算器５７−１〜５７−
Ｍにより、それぞれ自乗器５６−１〜５６−Ｍの各乗算
結果と記憶部５４からの固有値λiが乗算される。そし
て、係数乗算器５７−１〜５７−Ｍの各算出結果が加算
器５８により加算され、この加算結果＾ｙ2が一次畳み
込み処理部１８に送られる。
は、Ａ／Ｄ変換器１１により変換されたデジタルオーデ
ィオ信号ｘと記憶部５４に格納されていたＭ組の固有ベ
クトルｕiがそれぞれ一次畳み込み器５５−１〜５５−
Ｍにより乗算され、次いで自乗器５６−１〜５６−Ｍに
より、それぞれ一次畳み込み器５５−１〜５５−Ｍの乗
算結果が自乗され、次いで係数乗算器５７−１〜５７−
Ｍにより、それぞれ自乗器５６−１〜５６−Ｍの各乗算
結果と記憶部５４からの固有値λiが乗算される。そし
て、係数乗算器５７−１〜５７−Ｍの各算出結果が加算
器５８により加算され、この加算結果＾ｙ2が一次畳み
込み処理部１８に送られる。
【００３３】一次畳み込み処理部１８では、この加算結
果＾ｙ2に対して＾Ｈ1 -1が畳み込まれ、この処理結果＾
ｙ2’が加算器１３に送られる。加算器１３では遅延部
１２からのオーディオ信号ｘ’からこの処理結果＾
ｙ2’が差し引かれ、この処理結果ｘ’−＾ｙ2’がＤ／
Ａ変換器１４によりアナログ信号に変換されてスピーカ
３１に送られる。ここで、一次畳み込み処理部１８の処
理は、スピーカ３１により付加される一次成分を事前に
取り除くためであって、この処理結果＾ｙ2’がスピー
カ３１を通過する際に発生する二次歪みｙ2’を除去す
る。
果＾ｙ2に対して＾Ｈ1 -1が畳み込まれ、この処理結果＾
ｙ2’が加算器１３に送られる。加算器１３では遅延部
１２からのオーディオ信号ｘ’からこの処理結果＾
ｙ2’が差し引かれ、この処理結果ｘ’−＾ｙ2’がＤ／
Ａ変換器１４によりアナログ信号に変換されてスピーカ
３１に送られる。ここで、一次畳み込み処理部１８の処
理は、スピーカ３１により付加される一次成分を事前に
取り除くためであって、この処理結果＾ｙ2’がスピー
カ３１を通過する際に発生する二次歪みｙ2’を除去す
る。
【００３４】
【発明の効果】以上説明したように本発明によれば、歪
み除去対象のスピーカの二次ボルテラ時間応答を固有ベ
クトル及び固有値に分解し、この内、歪み除去に寄与が
大きいものを選択して一次畳み込み処理を行うようにし
たので、畳み込み処理の演算量を低減することができ
る。
み除去対象のスピーカの二次ボルテラ時間応答を固有ベ
クトル及び固有値に分解し、この内、歪み除去に寄与が
大きいものを選択して一次畳み込み処理を行うようにし
たので、畳み込み処理の演算量を低減することができ
る。
【図１】本発明に係るスピーカの二次歪み除去装置の一
実施形態を示すブロック図である。
実施形態を示すブロック図である。
【図２】図１の一次畳み込み器の構成を詳しく示すブロ
ック図である。
ック図である。
【図３】畳み込み係数の固有値を示す説明図である。
【図４】従来のスピーカの二次歪み除去装置を示すブロ
ック図である。
ック図である。
【図５】一次伝達関数の時間応答ｈ1（ｋ1）の一例であ
る。
る。
【図６】二次伝達関数の時間応答ｈ2（ｋ1，ｋ2）の一
例である。
例である。
１１ Ａ／Ｄ変換器（Ａ／Ｄ変換手段） １２ 遅延部 １３ 加算器（第２の加算手段） １４ Ｄ／Ａ変換器（Ｄ／Ａ変換手段） １８ 一次畳み込み処理部（第２の一次畳み込み手段） ３１ スピーカ ５０ 固有値分解処理部（固有値分解手段） ５２ 固有値分解部 ５３ 選択部 ５４ 記憶部 ５５，５５−１〜５５−Ｍ 一次畳み込み器（第１の一
次畳み込み手段） ５６，５６−１〜５６−Ｍ 自乗器（自乗手段） ５７，５７−１〜５７−Ｍ 係数乗算器（乗算手段） ５８ 加算器（第１の加算手段）
次畳み込み手段） ５６，５６−１〜５６−Ｍ 自乗器（自乗手段） ５７，５７−１〜５７−Ｍ 係数乗算器（乗算手段） ５８ 加算器（第１の加算手段）
Claims (1)
- 【請求項１】 歪み除去対象のスピーカの二次ボルテラ
時間応答を固有ベクトル及び固有値に分解し、前記固有
ベクトル及び固有値の内、歪み除去に寄与が大きいもの
を選択する固有値分解手段と、 アナログオーディオ信号をデジタル信号に変換するＡ／
Ｄ変換手段と、 前記Ａ／Ｄ変換手段により変換されたデジタルオーディ
オ信号と前記固有値分解手段により分解及び選択された
複数の固有ベクトルをそれぞれ一次畳み込み処理する複
数の第１の一次畳み込み手段と、 前記複数の第１の一次畳み込み手段の処理結果をそれぞ
れ自乗する複数の自乗手段と、 前記複数の自乗手段の処理結果と前記固有値分解手段に
より分解及び選択された複数の固有値をそれぞれ乗算す
る複数の乗算手段と、 前記複数の乗算手段の各処理結果を加算する第１の加算
手段と、 前記第１の加算手段の処理結果を一次畳み込み処理する
第２の一次畳み込み手段と、 前記Ａ／Ｄ変換手段により変換されたデジタルオーディ
オ信号から前記第２の一次畳み込み手段の処理結果を差
し引く第２の加算手段と、 前記第２の加算手段の処理結果をアナログオーディオ信
号に変換するＤ／Ａ変換手段とを備え、 前記Ｄ／Ａ変換手段の出力を前記歪み除去対象のスピー
カに印加するようにしたスピーカの二次歪み除去装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25547999A JP2001086585A (ja) | 1999-09-09 | 1999-09-09 | スピーカの二次歪み除去装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP25547999A JP2001086585A (ja) | 1999-09-09 | 1999-09-09 | スピーカの二次歪み除去装置 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JP2001086585A true JP2001086585A (ja) | 2001-03-30 |
Family
ID=17279343
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP25547999A Withdrawn JP2001086585A (ja) | 1999-09-09 | 1999-09-09 | スピーカの二次歪み除去装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP2001086585A (ja) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2014138226A (ja) * | 2013-01-16 | 2014-07-28 | Onkyo Corp | 音声再生装置 |
| US8996365B2 (en) | 2009-03-19 | 2015-03-31 | Yugengaisya Cepstrum | Howling canceller |
-
1999
- 1999-09-09 JP JP25547999A patent/JP2001086585A/ja not_active Withdrawn
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8996365B2 (en) | 2009-03-19 | 2015-03-31 | Yugengaisya Cepstrum | Howling canceller |
| JP2014138226A (ja) * | 2013-01-16 | 2014-07-28 | Onkyo Corp | 音声再生装置 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Panicker et al. | Parallel-cascade realizations and approximations of truncated Volterra systems | |
| JPH11261376A (ja) | 少ない乗算器によるディジタルｉｉｒフィルタ | |
| WO2007029536A1 (ja) | 雑音抑圧の方法及び装置並びにコンピュータプログラム | |
| US8108210B2 (en) | Apparatus and method to eliminate noise from an audio signal in a portable recorder by manipulating frequency bands | |
| JP5065784B2 (ja) | 同相成分抽出方法及び装置 | |
| US6519342B1 (en) | Method and apparatus for filtering an audio signal | |
| JP2001086585A (ja) | スピーカの二次歪み除去装置 | |
| JP6070572B2 (ja) | ディジタルフィルタ回路およびディジタルフィルタ処理方法 | |
| KR100454886B1 (ko) | 독립 성분 분석을 이용한 여파기 적응 알고리즘의 필터뱅크 접근 방법 | |
| US8463833B2 (en) | Filtering method and apparatus of low complexity fir filter, and recording medium thereof | |
| EP2809086B1 (en) | Method and device for controlling directionality | |
| JP2002175099A (ja) | 雑音抑制方法および雑音抑制装置 | |
| US20210183357A1 (en) | Method and system for implementing a modal processor | |
| JP3927701B2 (ja) | 音源信号推定装置 | |
| JP3831220B2 (ja) | 雑音抑圧方法及びその装置、雑音抑圧プログラム並びにそのプログラム記録媒体 | |
| JP3986457B2 (ja) | 入力信号推定方法、及び装置、入力信号推定プログラムならびにその記録媒体 | |
| JP3932960B2 (ja) | 信号成分抽出方法及び装置 | |
| CN107222808B (zh) | 一种高保真扬声器重放***设计方法 | |
| JP2006148453A (ja) | 信号推定方法、信号推定装置、信号推定プログラム及びその記録媒体 | |
| JP3602658B2 (ja) | スピーカの歪み除去装置を有する装置 | |
| JP4514153B2 (ja) | 音響装置 | |
| JPH06175691A (ja) | 音声強調装置と音声強調方法 | |
| JP2625015B2 (ja) | ゴースト除去装置 | |
| KR100863184B1 (ko) | 간섭 및 반향신호 제거를 위한 다단계 암묵 디콘볼루션방법 | |
| JPH05199071A (ja) | フィルタ設計方法および音声フィルタ |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A300 | Withdrawal of application because of no request for examination |
Free format text: JAPANESE INTERMEDIATE CODE: A300 Effective date: 20061205 |