(57)【特許請求の範囲】
【請求項１】
信号を入力し、前記信号を再生するともに、非線形歪みを発生するスピーカと、
オーディオ信号を入力し、前記オーディオ信号をそのまま前記スピーカに入力した場合に発生する非線形歪みを除去するようにした前記オーディオ信号を生成し、前記スピーカに出力するスピーカの歪み除去装置と、
を有し、
前記スピーカの歪み除去装置は、
前記オーディオ信号を入力しディジタルの信号に変換するＡ／Ｄ変換手段と、
前記Ａ／Ｄ変換手段の出力信号を、部分的に重複させながら長さＮで分割して取り込むフレーム分割手段と、
前記フレーム分割手段で分割された時間領域の信号を高速フーリエ変換して周波数領域の信号に変換するＦＦＴ演算手段と、
周波数領域において、Ｎ個の第１の係数を記憶している第１の記憶手段と、
前記第１の係数と前記ＦＦＴ演算手段の出力信号とを用いて、（数１）の右辺第１項の乗算を行う第１の乗算器と、
【数１】
0003602658.tif 000032
（  但しＷ（ｍ）は周波数領域における歪み除去装置の出力信号のｍ成分、
Ｇ１（ｍ）は第１の係数、
Ｘ（ｍ）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ成分、
ｍ、ｍ１及びｍ２は離散化された周波数軸上のポイントの数を表す整数値、
Ｇ２（ｍ１、ｍ２）は第２の係数、
Ｘ（ｍ１）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ１成分、
Ｘ（ｍ２）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ２成分を表す、  ）
周波数領域において、２次元のＮ×Ｎ個の第２の係数を記憶する第２の記憶手段と、
前記第２の係数と前記ＦＦＴ演算手段の出力信号とを用いて、（数１）の右辺第２項の乗算及び加算を行う乗加算器と、
前記第１の乗算器の出力信号と前記乗加算器の出力信号とを加算する加算器と、
前記加算器の出力信号を逆高速フーリエ変換して時間領域の信号に変換するＩＦＦＴ演算手段と、
前記ＩＦＦＴ演算手段の出力信号の一部を順次連結して出力する、フレーム合成手段と、
前記フレーム合成手段の出力信号をアナログのオーディオ信号に変換し、前記アナログのオーディオ信号を前記スピーカに出力するＤ／Ａ変換手段と、を具備する、
ことを特徴とする、スピーカの歪み除去装置を有する装置。
【請求項２】
信号を入力し、前記信号を再生するともに、非線形歪みを発生するスピーカと、
オーディオ信号を入力し、前記オーディオ信号をそのまま前記スピーカに入力した場合に発生する非線形歪みを除去するようにした前記オーディオ信号を生成し、前記スピーカに出力するスピーカの歪み除去装置と、
を有し、
前記スピーカの歪み除去装置は、
前記オーディオ信号を入力しディジタルの信号に変換するＡ／Ｄ変換手段と、
前記Ａ／Ｄ変換手段の出力信号を、部分的に重複させながら長さＮで分割して取り込むフレーム分割手段と、
前記フレーム分割手段で分割された時間領域の信号を高速フーリエ変換して周波数領域の信号に変換するＦＦＴ演算手段と、
前記スピーカの１次のインパルス応答の群遅延量と実質的に同等の遅延量を有するＮタップの遅延器のインパルス応答をフーリエ変換して求められる周波数領域の係数を記憶した、第１の記憶手段と、
前記第１の記憶手段の係数を読み出して、前記高速フーリエ変換手段の出力信号との乗算を行う第１の乗算器と、
周波数領域において、Ｎ個の第１の係数を記憶する第２の記憶手段と、
前記第１の係数と、前記第１の乗算器の出力信号とを用いて、（数２）の右辺第１項の乗算を行う第２の乗算器と、
【数２】
0003602658.tif 000033
（  但しＷ（ｍ）は周波数領域における歪み除去装置の出力信号のｍ成分、
Ｇ１（ｍ）は第１の係数、
Ｄ１（ｍ）は第１の記憶手段に記憶された係数のｍ成分、
Ｘ（ｍ）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ成分、
ｍ、ｍ１及びｍ２は離散化された周波数軸上のポイントの数を表す整数値、
Ｇ２（ｍ１、ｍ２）は第２の係数、
Ｄ１（ｍ１）は第１の記憶手段から読出された係数のｍ１成分、
Ｘ（ｍ１）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ１成分、
Ｄ１（ｍ２）は第１の記憶手段から読出された係数のｍ２成分、
Ｘ（ｍ２）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ２成分を表す、  ）
周波数領域において２次元のＮ×Ｎ個の、第２の係数を記憶する第３の記憶手段と、
前記第２の係数と、前記第１の乗算器の出力信号とを用いて、（数２）の右辺第２項の乗算及び加算を行う乗加算器と、
前記第２の乗算器の出力信号と前記乗加算器の出力信号とを加算する加算器と、
前記加算器の出力信号を逆高速フーリエ変換して時間領域の信号に変換するＩＦＦＴ演算手段と、
前記ＩＦＦＴ演算手段の出力信号の一部を順次連結して出力するフレーム合成手段と、
前記フレーム合成手段の出力信号をアナログのオーディオ信号に変換し、前記アナログのオーディオ信号を前記スピーカに出力するＤ／Ａ変換手段と、を具備する、
ことを特徴とする、スピーカの歪み除去装置を有する装置。
【請求項３】
前記乗加算器は、周波数領域において表わされる２次元のＮ×Ｎ個の前記第２の係数Ｇ２（ｍ１，ｍ２）のうち、（数３）かつ（数４）で定義される領域および、（数５）で定義される領域のみの係数と、
【数３】
0003602658.tif 000034
【数４】
0003602658.tif 000035
【数５】
0003602658.tif 000036
入力される信号Ｘｉ（ｍ）とを用いて、（数６）
【数６】
0003602658.tif 000037
に示す演算を行うことを特徴とする、請求項１又は２に記載のスピーカの歪み除去装置を有する装置。
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】
本発明は、オーディオ信号を再生するスピーカにおいて生じる高調波歪み及び混変調歪みを除去し、入力信号に対して高忠実再生を実現するための、歪み除去装置を有する装置に関するものである。
【０００２】
【従来の技術】
従来のスピーカでは、ＣＤプレーヤーなどからアンプを通じてスピーカに入力されるオーディオ信号を音声として再生する際に、オーディオ信号を忠実に再生する音声とともに、一般に非線形歪みと呼ばれる高調波歪みおよび混変調歪みの音声を発生していた。
【０００３】
スピーカから再生された音楽信号に、これらの非線形歪みが付加されていると非常に不快に感じられるので、スピーカ設計に際しては、これらの歪みを低減することが大きな課題である。特に、２次の高調波歪みおよび混変調歪み（以下簡単のため、これらを合わせて２次歪みと呼ぶ）は不快に感じられる。そこで、本発明は２次の高調波歪み及び混変調歪みを低減することを目的としている。
【０００４】
従来、これらの非線形歪みを除去するには、大きく分類すると２つの手法があった。一つは、長年に渡り取られていた手法であり、スピーカの形状、寸法、材質を工夫し、最も歪みが少なくなるように、設計するという方法である。他の一つは、近年になって提案されている方法であり、ＣＤプレーヤなどの音源から出力されたオーディオ電気信号を、歪み除去装置によって補正してからスピーカに入力し、最終的にスピーカから放射される音を歪みの少ない音にするという方法である。発明者は後者に属する方法について研究し、図７に示す１つのシステム（未公開）の発明を本発明に先だって完成した。
【０００５】
図７は、発明者がさきに完成した上記の先の発明のシステムにおける歪み除去装置７０５とスピーカ７０１をモデル化して示したブロック図である。
スピーカ７０１は、歪みの無いスピーカの伝達関数Ｈ１の系７０２と、スピーカの２次歪みを表す２次の伝達関数Ｈ２の系７０３と、それら２つの系の出力信号を加算器７０４で加算するシステムとして表すことができる。
【０００６】
歪み除去装置７０５は、ＣＤプレーヤなどからのオーディオ電気信号である入力信号ｘ（ｔ）をディジタルに変換するＡ／Ｄ変換器７０９と、伝達関数Ｇ１で１次元の畳込み演算を行う第１のフィルタ７０６と、２次元の伝達関数Ｇ２で２次元の畳込み演算を行う第２のフィルタ７０７を有している。さらに第１のフィルタ７０６と第２のフィルタ７０７の出力信号を加算する加算器７０８と、加算器の出力信号をアナログ信号に変換するＤ／Ａ変換器７０９とを有している。
【０００７】
このような一般的なスピーカの歪み除去装置７０５を有するシステムでは、入力信号ｘ（ｔ）と、スピーカ７０１からの出力信号ｙ（ｔ）の関係は、周波数領域における表示をすると、次の（数７）で表される。
【０００８】
【数７】
0003602658.tif 000002
【０００９】
（数７）の右辺第１項は、入力信号ｘ（ｔ）が第１のフィルタ７０６とスピーカの伝達関数Ｈ１の系７０２とを通過する成分を示す。（数７）の第２項の中括弧内の第２項は、入力信号ｘ（ｔ）が第１のフィルタ７０６とスピーカの歪みを表す伝達関数Ｈ２の系７０３とを通過する成分を示す。
【００１０】
ただし、（数７）の第２項は、周波数領域に変換されたディジタル信号のアドレスを表すｍについて、ｍ＝ｍ１＋ｍ２または｜ｍ１−ｍ２｜を満たす全ての組み合わせについての和である。なお、入力信号ｘ（ｔ）が第２のフィルタ７０７とスピーカの歪みを表す伝達関数Ｈ２の系７０３とを通過する成分は、他の項に比べて微少であり、無視している。
【００１１】
スピーカ７０１の歪み成分を除去するためには（数７）の第２項の中括弧内の２つの項が相殺されてその値が０になればよい。
【００１２】
次に、第１のフィルタ７０６の伝達関数Ｇ１、および、第２のフィルタ７０７の伝達関数Ｇ２の決定方法について述べる。
【００１３】
まず、第１のフィルタ７０６の伝達関数Ｇ１は、（数７）の第１項で表される成分がスピーカから出力される所望の出力信号に等しくなるように、決定すればよい。例えば、出力信号Ｙ（ｍ）がオーディオ電気信号Ｘ（ｍ）に等しくなるようにするためには、（数８）より、（数９）が導かれる。
【００１４】
【数８】
0003602658.tif 000003
【００１５】
【数９】
0003602658.tif 000004
【００１６】
また、出力信号Ｙ（ｍ）について、オーディオ電気信号がスピーカの伝達関数Ｈ１の影響を受けた信号になるようにしたい場合は、（数１０）より、第１のフィルタ７０６の伝達関数Ｇ１を（数１１）のようにすればよい。
【００１７】
【数１０】
0003602658.tif 000005
【００１８】
【数１１】
0003602658.tif 000006
【００１９】
次に、第２のフィルタ７０７の伝達関数Ｇ２の決定方法について述べる。（数７）の第２項の中括弧内の２つの項が相殺されればよいので、括弧内が０に等しいとして、Ｇ２について解くと、（数１２）となる。
【００２０】
【数１２】
0003602658.tif 000007
【００２１】
伝達関数Ｇ２を決定するには、（数１２）に、スピーカ７０１の１次の系７０２の伝達関数Ｈ１と、スピーカ７０１の２次の系７０３の２次歪みの伝達関数Ｈ２と、上記の方法で決定した第１のフィルタ７０６の伝達関数Ｇ１とを代入すればよい。
【００２２】
【発明が解決しようとする課題】
しかしながら、従来の歪み除去装置７０５においては、入力信号ｘ（ｔ）がＣＤプレーヤなどから入力される時間的に非常に長い任意のオーディオ信号である場合には、タップ数Ｍの第１のフィルタ７０６、および、タップ数Ｍ×Ｍの第２のフィルタ７０７においては、（数１３）に示すような畳込み演算をリアルタイム処理する必要がある。このリアルタイム処理において、特に第２項で表されている第２のフィルタ７０７による２次元の畳込み演算は演算量が膨大であるために、これを実現するには演算装置の規模が極めて大きくになるという問題があった。
【００２３】
【数１３】
0003602658.tif 000008
【００２４】
（数１３）は時間領域における入力信号ｘと、出力信号ｗの関係を示している。また、ｎは離散化された信号の順番を示し、時間に相当する。（数７）より、アナログの入力信号ｘ（ｔ）をディジタルに変換した信号ｘ（ｎ）の１サンプル信号について、第１のフィルタ７０６において必要な乗算回数はＭ回、加算回数はＭ−１回である。また、同様に（数１３）より、第２のフィルタ７０７で必要な乗算回数はＭ×Ｍ×２回、加算回数はＭ×（Ｍ−１）回である。
【００２５】
例えば、Ｍが１２８（Ｍ＝１２８）とすると、第２のフィルタ７０６において必要な乗算回数は３２７６８回、加算回数は１６３６８回である。すなわち、オーディオ信号ｘ（ｎ）の１サンプル当たりに非常に多くの演算量が必要とされ、この装置を実現するには、装置が非常に大がかりになってしまう。
【００２６】
本発明はこのような従来の問題点に鑑みてなされたものであって、一般的なスピーカに無限長の一般的なオーディオ信号が入力された場合でも、高調波歪みや混変調歪みを除去することが可能な、スピーカの歪み除去装置を有する装置を実現することを目的とする。
【００２７】
【課題を解決するための手段】
本発明のスピーカの歪み除去装置を有する装置の歪み除去装置では、オーディオ信号を、フレーム分割手段によって長さＮの時間領域の信号に分割し、周波数領域で、（数１４）に示す演算を行った後、時間領域の信号に再変換し、フレーム合成手段によってその一部を順次連結して出力する。これによって、著しく演算量を減らすことが可能であり、従来に比べて歪み除去装置を極めて小型化することが可能になる。以下に詳しい手段を説明する。
【００２８】
まず、アナログのオーディオ信号をＡ／Ｄ変換手段に入力し、ディジタル信号に変換する。次に、フレーム分割手段によって、Ａ／Ｄ変換手段から出力されるディジタル信号を、部分的に重複させながら長さＮで分断して取り込む。分断されたディジタル信号を、ＦＦＴ演算手段で高速フーリエ変換して時間領域の信号から周波数領域の信号へ変換する。
【００２９】
第１の記憶手段はあらかじめＮ個の第１の係数を記憶している。ここで、第１の乗算器は、ＦＦＴ演算手段の出力側および第１の記憶手段に接続されている。よって、第１の乗算器は、Ｎ個の第１の係数と、先程ＦＦＴ演算手段によって周波数領域の信号へ変換された長さＮの信号を、掛け合わせる。この演算を（数１４）の右辺第１項に示す。
【００３０】
一方、第２の記憶手段は、あらかじめＮ×Ｎ個の２次元の第２の係数を記憶している。ここで、乗加算器は、ＦＦＴ演算手段の出力側及び第２の記憶手段に接続されている。よって、乗加算器は、第２の係数と、前記ＦＦＴ演算手段によって周波数領域に変換された信号とを用いて、（数１４）の右辺第２項に示すような乗算及び加算を行う。
【００３１】
第１の乗算器の出力信号および乗加算器の出力信号は、それぞれ加算器に入力され、加算される。この加算が、（数１４）の右辺の第１項と第２項の加算である。
【００３２】
加算された信号は、（数１４）の左辺であり、これをＩＦＦＴ演算手段により逆フーリエ変換して、周波数領域の信号から時間領域の信号に変換する。
【００３３】
時間領域の信号に変換されたＩＦＦＴの出力信号は、フレーム合成手段によって、その一部が順次連結されて出力される。最後にＤ／Ａ変換手段は、フレーム合成手段によって連結されたディジタルのオーディオ信号をアナログに変換して出力する。
【００３４】
【発明の実施の形態】
本発明のスピーカの歪み除去装置を有する装置は、信号を入力し、前記信号を再生するともに、非線形歪みを発生するスピーカと、オーディオ信号を入力し、前記オーディオ信号をそのまま前記スピーカに入力した場合に発生する非線形歪みを除去するようにした前記オーディオ信号を生成し、前記スピーカに出力するスピーカの歪み除去装置と、を有し、前記スピーカの歪み除去装置は、前記オーディオ信号を入力しディジタルの信号に変換するＡ／Ｄ変換手段と、前記Ａ／Ｄ変換手段の出力信号を、部分的に重複させながら長さＮで分割して取り込むフレーム分割手段と、前記フレーム分割手段で分割された時間領域の信号を高速フーリエ変換して周波数領域の信号に変換するＦＦＴ演算手段と、周波数領域において、Ｎ個の第１の係数を記憶している第１の記憶手段と、前記第１の係数と前記ＦＦＴ演算手段の出力信号とを用いて、（数１４）の右辺第１項の乗算を行う第１の乗算器と、
【００３５】
【数１４】
0003602658.tif 000009
【００３６】
（  但しＷ（ｍ）は周波数領域における歪み除去装置の出力信号のｍ成分、Ｇ１（ｍ）は第１の係数、Ｘ（ｍ）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ成分、ｍ、ｍ１及びｍ２は離散化された周波数軸上のポイントの数を表す整数値、Ｇ２（ｍ１、ｍ２）は第２の係数、Ｘ（ｍ１）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ１成分、Ｘ（ｍ２）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ２成分を表す、  ）
周波数領域において、２次元のＮ×Ｎ個の第２の係数を記憶する第２の記憶手段と、前記第２の係数と前記ＦＦＴ演算手段の出力信号とを用いて、（数１４）の右辺第２項の乗算及び加算を行う乗加算器と、前記第１の乗算器の出力信号と前記乗加算器の出力信号とを加算する加算器と、前記加算器の出力信号を逆高速フーリエ変換して時間領域の信号に変換するＩＦＦＴ演算手段と、前記ＩＦＦＴ演算手段の出力信号の一部を順次連結して出力するフレーム合成手段と、前記フレーム合成手段の出力信号をアナログのオーディオ信号に変換し、前記アナログのオーディオ信号を前記スピーカに出力するＤ／Ａ変換手段と、を具備する。
【００３７】
また、スピーカの歪み除去装置を有する装置は、信号を入力し、前記信号を再生するともに、非線形歪みを発生するスピーカと、オーディオ信号を入力し、前記オーディオ信号をそのまま前記スピーカに入力した場合に発生する非線形歪みを除去するようにした前記オーディオ信号を生成し、前記スピーカに出力するスピーカの歪み除去装置と、を有し、前記スピーカの歪み除去装置は、前記オーディオ信号を入力しディジタルの信号に変換するＡ／Ｄ変換手段と、前記Ａ／Ｄ変換手段の出力信号を、部分的に重複させながら長さＮで分割して取り込むフレーム分割手段と、前記フレーム分割手段で分割された時間領域の信号を高速フーリエ変換して周波数領域の信号に変換するＦＦＴ演算手段と、前記スピーカの１次のインパルス応答の群遅延量と実質的に同等の遅延量を有するＮタップの遅延器のインパルス応答をフーリエ変換して求められる周波数領域の係数を記憶した、第１の記憶手段と、前記第１の記憶手段の係数を読み出して、前記高速フーリエ変換手段の出力信号との乗算を行う第１の乗算器と、周波数領域において、Ｎ個の第１の係数を記憶する第２の記憶手段と、前記第１の係数と、前記第１の乗算器の出力信号とを用いて、（数１５）の右辺第１項の乗算を行う第２の乗算器と、
【００３８】
【数１５】
0003602658.tif 000010
【００３９】
（  但しＷ（ｍ）は周波数領域における歪み除去装置の出力信号のｍ成分、Ｇ１（ｍ）は第１の係数、Ｄ１（ｍ）は第１の記憶手段に記憶された係数のｍ成分、Ｘ（ｍ）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ成分、ｍ、ｍ１及びｍ２は離散化された周波数軸上のポイントの数を表す整数値、Ｇ２（ｍ１、ｍ２）は第２の係数、Ｄ１（ｍ１）は第１の記憶手段から読出された係数のｍ１成分、Ｘ（ｍ１）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ１成分、Ｄ１（ｍ２）は第１の記憶手段から読出された係数のｍ２成分、Ｘ（ｍ２）は入力信号を、離散化した後に高速フーリエ変換によって周波数領域へ変換した信号のｍ２成分を表す、  ）
周波数領域において２次元のＮ×Ｎ個の、第２の係数を記憶する第３の記憶手段と、前記第２の係数と、前記第１の乗算器の出力信号とを用いて、（数１５）の右辺第２項の乗算及び加算を行う乗加算器と、前記第２の乗算器の出力信号と前記乗加算器の出力信号とを加算する加算器と、前記加算器の出力信号を逆高速フーリエ変換して時間領域の信号に変換するＩＦＦＴ演算手段と、前記ＩＦＦＴ演算手段の出力信号の一部を順次連結して出力するフレーム合成手段と、前記フレーム合成手段の出力信号をアナログのオーディオ信号に変換し、前記アナログのオーディオ信号を前記スピーカに出力するＤ／Ａ変換手段と、を具備する。
【００４０】
前記乗加算器は、周波数領域において表わされる２次元のＮ×Ｎ個の前記第２の係数Ｇ２（ｍ１，ｍ２）のうち、（数１６）かつ（数１７）で定義される領域および、（数１８）で定義される領域のみの係数と、
【００４１】
【数１６】
0003602658.tif 000011
【００４２】
【数１７】
0003602658.tif 000012
【００４３】
【数１８】
0003602658.tif 000013
【００４４】
入力される信号Ｘｉ（ｍ）とを用いて、（数１９）
【００４５】
【数１９】
0003602658.tif 000014
【００４６】
に示す演算を行う。
【００４７】
【実施例】
以下、図１ないし図６を参照して本発明の実施例を説明する。
《実施例１》本発明の実施例１について、図１を参照して説明する。本実施例は従来に比べて、大幅に演算量を削減する事の出来る歪み除去装置を有する装置を提供するものである。この歪み除去装置１０は、オーディオ信号が入力されるＡ／Ｄ変換手段１１と、数値を表すＮとＮ１の関係がＮ＞Ｎ１であって前記Ａ／Ｄ変換手段のディジタル出力信号を（Ｎ−Ｎ１＋１）個毎に長さＮで分割し取り込むフレーム分割手段１２とを有する。さらに前記フレーム分割手段で分割された時間領域の信号を高速フーリエ変換して周波数領域の信号に変換するＦＦＴ演算手段１３と、周波数領域において、Ｎ個の、第１の係数を記憶する第１の記憶手段１５と、前記第１の係数と前記ＦＦＴ演算手段１３の出力信号とを用いて、（数２０）の第１項の乗算を行う乗算器１４とを有する。
【００４８】
【数２０】
0003602658.tif 000015
【００４９】
さらに、周波数領域において２次元のＮ×Ｎ個の第２の係数を記憶する第２の記憶手段１７と、前記第２の係数と前記ＦＦＴ演算手段の出力信号とを用いて（数２０）の第２項の乗算と加算を行う乗加算器１６と、前記第１の乗算器１４の出力信号と前記乗加算器１６の出力信号とを加算する加算器１８とを有している。さらに前記加算器の出力信号を逆高速フーリエ変換するＩＦＦＴ演算手段１９と、前記ＩＦＦＴ演算手段の出力信号のＮ１個めからＮ個めのデータを順次連結して出力するフレーム合成手段２０と、前記フレーム合成手段の出力信号をアナログのオーディオ信号に変換するＤ／Ａ変換手段２１と、を具備する。
【００５０】
このような構成において、ＣＤプレーヤなどからアナログのオーディオ信号ｘ（ｔ）が歪み除去装置１０に入力されると、Ａ／Ｄ変換手段１１はオーディオ信号ｘ（ｔ）をディジタルの信号ｘ（ｎ）に変換する。フレーム分割手段１２は、前記Ａ／Ｄ変換手段の出力信号ｘ（ｎ）が（Ｎ−Ｎ１＋１）個入力される毎に、出力信号ｘ（ｎ）をそれ以前に入力された長さＮの信号ｘｉ（ｎ）に分割し、出力する。ただし、ここで、Ｎ＞Ｎ１であり、Ｎ１は、図７の従来の歪み除去装置７０５の第１のフィルタ７０６のタップ数Ｍに対応する。分割された長さＮの信号は、高速フーリエ変換手段１３に入力され、周波数領域の信号Ｘｉ（ｍ）に変換される。「ｍ」は離散化された周波数軸上のポイントの数を表す整数値である。
【００５１】
信号Ｘｉ（ｍ）は、第１の乗算器１４、第１の記憶手段１５、乗加算器１６、第２の記憶手段１７及び加算器１８によって、（数２１）に示す演算が施され、信号Ｗｉ（ｍ）となる。
【００５２】
【数２１】
0003602658.tif 000016
【００５３】
（数２１）の演算について詳しく述べる。第１の乗算器１４は、第１の記憶手段１５に格納されているＮ個の第１の係数Ｇ１（ｍ）と、高速フーリエ変換手段１３の出力信号Ｘｉ（ｍ）をそれぞれのｍについて掛け合わせる。これが、（数２１）の第１項の演算である。一方、乗加算器１６は、第２の記憶手段１７に格納されているＮ×Ｎ個の第２の係数Ｇ２（ｍ１、ｍ２）と、信号Ｘｉ（ｍ１）、Ｘｉ（ｍ２）とを用いて、（数２１）の第２項の演算を行う。「ｍ１」、「ｍ２」はそれぞれ前記の整数値ｍの特定の値を表す。加算器１８は、乗算器１４の出力信号と、乗加算器１６の出力信号とをそれぞれのｍについて足し合わせて、信号Ｗｉ（ｍ）を出力する。
【００５４】
加算器１８の出力信号Ｗｉ（ｍ）は逆高速フーリエ変換を行うＩＦＦＴ演算手段１９に入力され、時間領域の信号ｗｉ（ｎ）に変換される。フレーム合成手段２０は、信号ｗｉ（ｎ）のうち、信号ｗｉ（Ｎ１＋１）から信号ｗｉ（Ｎ）を切り出し、ｗ１（Ｎ１＋１）、・・・ｗ１（Ｎ）、ｗ２（Ｎ１＋１）、・・・ｗ２（Ｎ）、ｗ３（Ｎ１＋１）、・・・、の順に連続して出力する。ここで、「ｉ」はフレーム分割手段１２で分割されたフレームの番号である。Ｄ／Ａ変換手段２１は、フレーム合成手段２０の出力信号をアナログ信号ｗ（ｔ）に変換し、出力する。ｗ（ｔ）は、スピーカ２２に入力される。
【００５５】
入力されるオーディオ信号とともに高調波歪みおよび混変調歪みを発生するスピーカ２２は、スピーカの線形の伝達関数Ｈ１の系２３と、歪みを発生する非線形の２次の伝達関数Ｈ２の系２４と、加算器２５とによって表わされる。スピーカ２２の入力信号と出力音圧の関係は、周波数領域での表示を行うと、（数２２）のように表わされる。
【００５６】
【数２２】
0003602658.tif 000017
【００５７】
ここで、ＣＤプレーヤなどからのオーディオ信号ｘ（ｔ）と、スピーカからの出力音圧ｙ（ｔ）の関係を、周波数領域で表わすと、（数２０）の信号Ｗ（ｍ）を（数２２）に代入して消去する事により、（数７）となる。
【００５８】
（数７）の第１項は、オーディオ信号ｘ（ｔ）が、乗算器１４とスピーカ２２の伝達関数Ｈ１の系２３とを通過する成分を示す。（数７）の右辺第２項の中括弧内の第１項は、オーディオ信号ｘ（ｔ）が乗加算器１６とスピーカ２２の伝達関数Ｈ１の系２３とを通過する成分を示す。（数７）の右辺第２項の中括弧内の第２項は、オーディオ信号ｘ（ｔ）が乗算器１４とスピーカ２２の歪みを表わす伝達関数Ｈ２の系２４とを通過する成分を示す。
【００５９】
ただし、（数７）の右辺第２項のｍは、ｍ＝ｍ１＋ｍ２または｜ｍ１−ｍ２｜を満たす全ての組み合わせについての和である。なお、オーディオ信号ｘ（ｔ）が乗加算器１６とスピーカ２２の歪みを表わす伝達関数Ｈ２の系２４とを通過する成分は、他の項に比べて微少であり、無視している。
【００６０】
スピーカ２２の歪み成分を除去するためには（数７）の第２項の括弧内の２つの項が相殺されてその値が０になれば良い。
【００６１】
次に、第１の記憶手段１５の係数Ｇ１（ｍ）および、第２の記憶手段１７の係数Ｇ２（ｍ１、ｍ２）の決定方法について述べる。
【００６２】
まず、第１の記憶手段１５の係数Ｇ１（ｍ）は、（数７）の第１項で表される成分が、スピーカから出力される所望の出力音圧に等しくなるように決定すれば良い。例えば、出力信号Ｙ（ｍ）がオーディオ信号Ｘ（ｍ）に等しくなるようにするためには、（数８）より、（数９）が導かれる。
【００６３】
あるいはまた、出力信号Ｙ（ｍ）について、オーディオ信号がスピーカ２２の線形の系２３の伝達特性Ｈ１の影響を受けた信号になるようにしたい場合は、（数１０）より、第１の記憶手段１５の係数Ｇ１を（数１１）のようにすればよい。
【００６４】
次に、第２の記憶手段１７の係数Ｇ２（ｍ１、ｍ２）の決定方法について述べる。（数７）の右辺第２項の中括弧内の２つの項が相殺されれば良いので、中括弧内が０に等しいとして、Ｇ２（ｍ１、ｍ２）について解くと、（数１２）となる。
【００６５】
Ｇ２（ｍ１、ｍ２）を決定するには、（数１２）に、スピーカ２２の１次の系２３の伝達関数Ｈ１と、高調波歪みおよび混変調歪みを表わす２次の系２４の伝達関数Ｈ２と、上記の方法で決定した第１の記憶手段１５の係数Ｇ１（ｍ）とを代入すれば良い。
【００６６】
以上のように、第１の記憶手段１５の係数Ｇ１（ｍ）および、第２の記憶手段１７の係数Ｇ２（ｍ１、ｍ２）を決定すれば、歪み除去装置１０によってスピーカ２２の高調波歪み及び混変調歪みを除去する事ができる。
【００６７】
本発明の実施例１における歪み除去装置１０の、第１の乗算器１４および乗加算器１６における、乗算回数および加算回数について以下に述べる。
【００６８】
第１の乗算器１４において必要な乗算回数は、ディジタルに変換されたオーディオ信号ｘ（ｎ）の１個につき（数２３）に示す回数となる。
【００６９】
【数２３】
0003602658.tif 000018
【００７０】
また、乗加算器１６において必要な乗算回数は、ディジタルに変換されたオーディオ信号ｘ（ｎ）の１個につき（数２４）に示す回数となり、加算回数は（数２５）に示す回数となる。
【００７１】
【数２４】
0003602658.tif 000019
【００７２】
【数２５】
0003602658.tif 000020
【００７３】
第１の乗算器１４と乗加算器１６と加算器１８とを合わせた系において必要な乗算回数および加算回数をまとめると、乗算回数は（数２６）に示す回数となり、加算回数は（数２７）に示す回数となる。
【００７４】
【数２６】
0003602658.tif 000021
【００７５】
【数２７】
0003602658.tif 000022
【００７６】
一方、従来の時間領域における畳込み演算による方法では、畳込み演算の部分で必要な乗算回数は（数２８）に示す回数となり、加算回数は（数２９）に示す回数となる。
【００７７】
【数２８】
0003602658.tif 000023
【００７８】
【数２９】
0003602658.tif 000024
【００７９】
例えば、上記の（数２６）、（数２７）、（数２８）、（数２９）に、Ｎ＝２５６、Ｎ１＝１２８を代入すると、本発明による方法では、乗算回数が約１０１８回となり、加算回数が５０８回となる。一方、従来の時間領域における畳込み演算による方法では、乗算回数が３２７６８回となり、加算回数が１６５１１回となる。以上のように本発明によれば乗算回数、加算回数とも従来の方法の約３０分の１になり本発明の有効性が分かる。
【００８０】
以上より、本発明の歪み除去装置１０においては演算量が大幅に削減され、装置の規模を縮小する事が出来る。
【００８１】
《実施例２》
本発明の実施例２について、図２ないし図５を参照して説明する。
【００８２】
実施例２の構成は、以下の点で実施例１で説明したものと異なっている。図１の第２の記憶手段は、Ｎ×Ｎ個の第２の係数Ｇ２（ｍ１、ｍ２）のうち、図２に斜線で示す領域２０１および領域２０２の係数のみを記憶していれば良い。図２に示した領域２０１を数式で表わすと（数３０）かつ（数３１）で定義される領域となる。また、領域２０２を数式で表わすと（数３２）で定義される領域となる。
【００８３】
【数３０】
0003602658.tif 000025
【００８４】
【数３１】
0003602658.tif 000026
【００８５】
【数３２】
0003602658.tif 000027
【００８６】
前記の実施例１の歪み除去装置１０では、歪み除去装置１０内部の乗加算器１６における（数２１）の右辺第２項の演算の乗算回数及び加算回数が装置１０全体の演算量の多くを占めていたが、実施例２では、乗加算器１６における演算量を減らし、装置１０の規模の縮小が実現できる。
【００８７】
第２の係数Ｇ２（ｍ１、ｍ２）の特長について述べる。第２の係数Ｇ２（ｍ１、ｍ２）においては、図３の斜線部の領域３０１の係数とそれ以外の領域３０２の係数がｍ１＝ｍ２を対称軸とする線対象の関係にある。また、図４の斜線部の領域４０１の係数と斜線部ではない領域４０２の係数が（Ｎ／２，Ｎ／２）を点対称の中心とする共役関係にある。また、図４の斜線部の領域４０３の係数と斜線部ではない領域４０４の係数が（Ｎ／２，Ｎ／２）を点対称の中心とする共役関係にある。また、図４の斜線部の領域４０５の係数と斜線部ではない領域４０６の係数が（Ｎ／２，Ｎ／２）を点対称の中心とする共役関係にある。また、サンプリング定理により考えると、（数２１）の第２項の積和演算において図５の斜線部の領域５０１、５０２とそれ以外の領域５０３、５０４では積和演算の結果が共役になるので、実際の計算では斜線部の領域のみを考慮すれば良い。そこで、図３、図４、図５の積集合を取ると、図２の斜線部の領域２０１および２０２となり、（数２１）の右辺第２項の演算を行う際には、乗加算器１６において、図２の斜線部の領域２０１及び２０２について（数３３）に示すように積和演算を行えば良い。
【００８８】
【数３３】
0003602658.tif 000028
【００８９】
これにより、ディジタルに変換されたオーディオ信号の１サンプルあたりに必要な、乗加算器１６における乗算回数は（数３４）に示す回数となり、また、加算回数は（数３５）に示す回数となる。
【００９０】
【数３４】
0003602658.tif 000029
【００９１】
【数３５】
0003602658.tif 000030
【００９２】
例えば、Ｎ＝２５６、Ｎ１＝１２８を（数３４）、（数３５）に代入すると、乗算回数は約１６１回、加算回数は約７９回となり、乗加算器１６における演算量が大幅に削減され、よって、歪み除去装置１０の規模を縮小する事ができる。
【００９３】
《実施例３》
本発明の実施例３について図６を参照して説明する。
【００９４】
図６の実施例３の構成において、実施例１と異なる部分は、高速フーリエ変換手段３３の後段に第２の乗算器４６を設け、係数Ｄ１（ｍ）を格納した第１の記憶手段４７を前記第２の乗算器４６に設けたことである。
【００９５】
スピーカ４２のインパルス応答ｈ１（ｔ）は、一般に、入力信号ｘ（ｔ）に対して群遅延を含んだ特性となっている。スピーカ４２の１次の系４３の伝達関数Ｈ１はインパルス応答ｈ１をフーリエ変換して求められる。しかし、インパルス応答ｈ１に群遅延がある場合に（数９）によって第１の係数Ｇ１（ｍ）を決定すると、時間的にランダムに変化する一般的なオーディオ信号を歪み除去装置３０に入力した場合に、歪み除去装置３０とスピーカ４２を含めた全体の系の入出力特性の因果律が犯され、歪みが所望通りに除去されないという場合がある。
【００９６】
本発明ではこの問題を解決するために、図６に示すように、前記第２の乗算器４６と第１の記憶手段４７を設ける。
【００９７】
以上の構成において、まず、第１の記憶手段４７には、インパルス応答ｈ１の群遅延量τとほぼ同等の遅延作用のあるＮタップの遅延器のインパルス応答特性をフーリエ変換した係数Ｄ１（ｍ）を格納する。
【００９８】
第２の乗算器４６は、信号Ｘ（ｍ）に対して、第１の記憶手段４７より読み出した係数Ｄ１（ｍ）を、ｍについてそれぞれ掛け合わせて出力する。これにより、第１の乗算器３４を通過する成分には、スピーカ４２の１次の系４３の伝達特性の群遅延量に見合った遅延が施されたこととなる。また、乗加算器３６を通る成分については、第２の乗算器４６において、第１の記憶手段４７より読み出された係数Ｄ１（ｍ１）が、信号Ｘ（ｍ１）に対して掛け合わされ、信号Ｘ（ｍ２）に対しては第１の記憶手段４７より読み出された係数Ｄ１（ｍ２）がそれぞれ掛け合わされて出力される。その結果スピーカ４２の高調波歪みや混変調歪みの群遅延量に見合った遅延が施された出力となる。これにより、第１の乗算器３４、第２の乗算器４６、乗加算器３６、加算器３８による信号Ｘ（ｍ）と信号Ｗ（ｍ）の関係は（数３６）となる。
【００９９】
【数３６】
0003602658.tif 000031
【０１００】
以上の作用により、歪み除去装置３０に、周期性の無い一般的なオーディオ信号が入力された場合でも、スピーカの歪みを除去する事が可能となる。
【０１０１】
なお、本実施例と同様の効果を得るには、以下の構成によっても可能である。まず、第１の記憶手段４７および第２の乗算器４６を設けず、図１に示した実施例１あるいは実施例２と同様の構成とする。第２の記憶手段３５には、前記第１の記憶手段４７に格納されている係数Ｄ１（ｍ）と、実施例３で第２の記憶手段３５に格納されていた係数Ｇ１（ｍ）とを、ｍについてそれぞれ掛け合わせた結果のデータＧ１（ｍ）Ｄ１（ｍ）をあらかじめ格納する。第３の記憶手段３７には、前記第１の記憶手段４７に格納されている係数Ｄ１（ｍ１）と、実施例３で第３の記憶手段３７に格納されていた係数Ｇ２（ｍ１、ｍ２）とを、ｍ１、ｍ２についてそれぞれ掛け合わせた結果のデータＧ２（ｍ１、ｍ２）Ｄ１（ｍ１）Ｄ１（ｍ２）をあらかじめ格納する。
【０１０２】
以上の構成により、第１の乗算器３４は、第２の記憶手段３５より読み出したデータＧ１（ｍ）Ｄ１（ｍ）を、信号Ｘ（ｍ）に対して、それぞれのｍについて掛け合わせて出力する。これにより、第１の乗算器３４を通過する成分には、スピーカの１次の系４３のインパルス応答ｈ１の群遅延量に見合った遅延が施されたこととなる。乗加算器３６は、第３の記憶手段３７より読み出したデータＧ２（ｍ１、ｍ２）Ｄ１（ｍ１）Ｄ１（ｍ２）を、信号Ｘ（ｍ）に対して、それぞれのｍ１、ｍ２について掛け合わせて出力する。ただしここで、ｍ、ｍ１、ｍ２の関係は、第１実施例あるいは第２実施例と同様に、ｍ＝ｍ１＋ｍ２あるいはｍ＝｜ｍ１−ｍ２｜を満たすものとする。
【０１０３】
これにより、第１の乗算器３４、乗加算器３６、加算器３８による信号Ｘ（ｍ）と信号Ｗ（ｍ）の関係は、（数３６）のようになる。
【０１０４】
以上の作用により、歪み除去装置４２に、周期性の無い一般的なオーディオ信号が入力された場合でも、スピーカの歪みを除去する事が可能となる。
【０１０５】
【発明の効果】
本発明によれば、Ａ／Ｄ変換されたオーディオ信号を所定の長さに分割して高速フーリエ変換を行い、周波数領域の信号に変換する。そして歪み除去装置における１次元および２次元の畳込み演算を、周波数領域において乗算器と乗加算器で行うことにより、少ない乗算回数と加算回数で演算処理が可能となり規模の小さな演算装置によりスピーカの歪み除去装置を有する装置が実現できる。
【０１０６】
また、上記の効果に加えて、前記乗算器および乗加算器の前段にさらにそれぞれ乗算器を設け、スピーカのインパルス応答の群遅延量とほぼ同等の遅延作用を有する遅延器のインパルス応答特性をフーリエ変換した係数Ｄ１（ｍ）の乗算を施す事により、スピーカの伝達特性に群遅延量が含まれる場合でも所望の歪み除去効果を得る事ができる。
【０１０７】
さらに、周波数領域において、２次元の所定領域の係数と高速フーリエ変換出力とを用いて乗算及び加算を行うことにより、さらに演算回数を削減することができる。
【図面の簡単な説明】
【図１】本発明の実施例１における、スピーカの歪み除去装置とスピーカの構成を示すブロック図
【図２】本発明の実施例２において、第２の記憶手段に格納すべき係数の領域を示す図
【図３】本発明の実施例２において、第２の記憶手段に格納すべき係数の対称性を示す図
【図４】本発明の実施例２において、第２の記憶手段に格納すべき係数の共役性を示す図
【図５】本発明の実施例２において、第２の記憶手段に格納すべき係数のうち、サンプリング定理により考えられる、有効な範囲を示す図
【図６】本発明の実施例３における、スピーカの歪み除去装置とスピーカの構成を示すブロック図
【図７】従来のスピーカの歪み除去装置とスピーカの構成を示すブロック図
【符号の説明】
１０、３０  歪み除去装置
１１、３１  Ａ／Ｄ変換手段
１２、３２  フレーム分割手段
１３、３３  高速フーリエ変換手段
１４、３４  第１の乗算器
１５、４７  第１の記憶手段
１６、３６  乗加算器
１７、３５  第２の記憶手段
１８、３８  加算器
１９、３９  逆高速フーリエ変換手段
２０、４０  フレーム合成手段
２１、４１  Ｄ／Ａ変換手段
２２、４２  スピーカ
２３、４３  線形の伝達関数Ｈ１の系
２４、４４  非線形歪みの２次の伝達関数Ｈ２の系
２５、４５  加算器
４６  第２の乗算器
３７  第３の記憶手段
２０１  実施例２において、第２の記憶手段に格納すべき係数の領域１
２０２  実施例２において、第２の記憶手段に格納すべき係数の領域２
【図１】
0003602658.tif 000038
【図２】
0003602658.tif 000039
【図３】
0003602658.tif 000040
【図４】
0003602658.tif 000041
【図５】
0003602658.tif 000042
【図６】
0003602658.tif 000043
【図７】
0003602658.tif 000044