(57)【要約】【課題】信号経路の周波数特性のばらつきが大きい場合でも十分な精度が得られるようにする。【解決手段】第１の推定手段２２は各Ａ／Ｄ変換器１２のサンプル値と次のサンプル値の間の値を、他のＡ／Ｄ変換器のサンプル値を用いて推定し、その推定値とサンプル値とを時系列に第１のイコライザ２３へ与えて各信号経路の差を補償し、第１の信号切換器２４から出力する。また、第１の信号切換器２４から出力された信号に対して第２の推定手段３０による推定処理と、第２のイコライザ３１による補正処理とを行い、その第２のイコライザ３１の出力値を順次選択するようにしたので、より精度の高いＡ／Ｄ変換処理が可能となる。【選択図】図２
2010141721.tif 000002
【特許請求の範囲】
【請求項１】
信号を入力するための入力端子（１０ａ）と、
複数のＡ／Ｄ変換器（１２）と、
前記入力端子からの入力信号を前記Ａ／Ｄ変換器にそれぞれ入力する信号分配器（１１）と、
前記各Ａ／Ｄ変換器に対し、サンプリングのための所定周期のクロックを、前記Ａ／Ｄ変換器数で前記所定周期を割って得られる時間にほぼ等しい時間差で所定順に且つ循環的に与えるサンプリング制御部（２１）と、
前記複数のＡ／Ｄ変換器の１つを基準とし、前記入力端子から前記各Ａ／Ｄ変換器の出力端子までのそれぞれの周波数特性と前記入力端子から前記基準のＡ／Ｄ変換器の出力端子までの周波数特性との差の特性を満たすインパルス応答を有するフィルタの係数を予め記憶しているＡＤ特性テーブル（２５）と、
前記Ａ／Ｄ変換器毎にそれぞれ設けられ、前記複数のＡ／Ｄ変換器によって変換出力されるサンプル値および前記ＡＤ特性テーブルに記憶されている係数に基づいて、前記クロックを受けたＡ／Ｄ変換器がサンプル値を更新するタイミングに他のＡ／Ｄ変換器が変換処理をおこなったと仮定して得られるサンプル値をそれぞれ推定する複数の第１の推定手段（２２）と、
前記基準のＡ／Ｄ変換器と各Ａ／Ｄ変換器の周波数特性の差をそれぞれ相殺する周波数特性を満たすインパルス応答を有するフィルタの係数を予め記憶しているイコライザ係数テーブル（２７）と、
前記各第１の推定手段にそれぞれ対応して設けられ、各第１の推定手段の出力値に対して、前記イコライザ係数テーブルに記憶されている係数に基づくフィルタリングをそれぞれ行って、誤差補正された値をそれぞれ出力する複数の第１のイコライザ（２３）と、
前記複数の第１のイコライザの出力値を受け、前記各Ａ／Ｄ変換器が前記クロックによってサンプリングする順に前記複数の第１のイコライザの出力値を選択して出力する第１の信号切換器（２４）と、
前記各第１のイコライザに対応して設けられ、前記ＡＤ特性テーブルに記憶されている係数に基づいて、前記第１の信号切換器から出力された信号に対して前記対応する各イコライザと逆の周波数補正処理を行うとともに、該補正処理で得られたデータのうち、前記対応するＡ／Ｄ変換器のサンプリングタイミングのデータを元のサンプル値に入れ換えて出力する複数の第２の推定手段（３０）と、
前記各第２の推定手段の出力値に対し、前記入力端子から前記各Ａ／Ｄ変換器までの信号経路の周波数特性の差異を補正するための処理を前記イコライザ係数テーブルに記憶された係数に基づいて行う複数の第２のイコライザ（３１）と、
前記各第２のイコライザの出力値を受け、前記各Ａ／Ｄ変換器が前記クロックによってサンプリングする順に前記複数の第２のイコライザの出力値を選択して出力する第２の信号選択手段（３２）とを備えたＡ／Ｄ変換装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、アナログの入力信号を複数のＡ／Ｄ変換器に共通に入力し各Ａ／Ｄ変換器のサンプリングタイミングをずらすことにより、入力信号に対して高速なサンプリングを等価的に行うインタリーブ方式のＡ／Ｄ変換装置において、その精度を高めるための技術に関する。
【背景技術】
【０００２】
図１１はインタリーブ方式のＡ／Ｄ変換装置１０の基本構成を示し、図１２はその動作を示している。
【０００３】
このＡ／Ｄ変換装置１０は、入力端子１０ａに入力される図１２の（ａ）のようなアナログの入力信号ｘ（ｔ）を、信号分配器１１によって複数Ｎ本の信号経路に分岐して、Ｎ個のＡ／Ｄ変換器１２０〜１２Ｎ−１にそれぞれ入力する。
【０００４】
サンプリング制御部１３は、図１２の（ｂ１）〜（ｂＮ）に示すように、それぞれが周期Ｔを持ち、位相がΔＴ（＝Ｔ／Ｎ）ずつシフトされたサンプリング用のクロックＣ０〜ＣＮ−１を生成してそれぞれＡ／Ｄ変換器１２０〜１２Ｎ−１に与えるとともに、図１２の（ｄ）のように、各Ａ／Ｄ変換器１２０〜１２Ｎ−１のうちサンプリングを行うＡ／Ｄ変換器を指定する指定信号ＡＤＮＵＭを信号切換器１４に与える。
【０００５】
各Ａ／Ｄ変換器１２０〜１２Ｎ−１は、クロックＣ０〜ＣＮ−１をそれぞれ受けたときの入力値ｘ（Ｐ）、ｘ（Ｐ＋１）、ｘ（Ｐ＋２）、…をサンプリングしてデジタル値に変換し、図１２の（ｃ１）〜（ｃＮ）のように、各サンプル値Ｘ０，Ｐ、Ｘ１，Ｐ＋１、Ｘ２，Ｐ＋２、…をそれぞれ信号切換器１４に出力する。
【０００６】
信号切換器１４は、各Ａ／Ｄ変換器１２０〜１２Ｎ−１のうち、指定信号ＡＤＮＵＭで指定されたＡ／Ｄ変換器から出力されるサンプル値Ｘ０，Ｐ、Ｘ１，Ｐ＋１、Ｘ２，Ｐ＋２、…を順次選択して、図１２の（ｅ）のように、サンプル値がそのサンプリング順に並んだデジタル信号列Ｙ（ｎ）を出力端子１０ｂに出力する。
【０００７】
このようにして得られるデジタル信号列Ｙ（ｎ）は、入力信号ｘ（ｔ）をクロック周期Ｔの１／Ｎのサンプリング周期ΔＴでサンプリングして得られるものと等価となり、低速なＡ／Ｄ変換器で高速なサンプリングが行える。
【０００８】
ところが、上記Ａ／Ｄ変換装置１０のように、入力信号ｘ（ｔ）を複数のＡ／Ｄ変換器１２０〜１２Ｎ−１に分配入力する場合、信号分配器１２自身の分配特性や分配経路の周波数特性の違い、各Ａ／Ｄ変換器１２０〜１２Ｎ−１の周波数特性の違いおよびサンプリングクロックの理想タイミングからのずれによって、得られたサンプル値を信号処理した結果に誤差を発生させる。
【０００９】
この問題を解決するための技術として、本願出願人は、入力端子から各Ａ／Ｄ変換器までの信号経路の周波数特性の差をなくして等しくするための等化処理（イコライズ処理）を行う技術を提案している（特許文献１）。
【００１０】
【特許文献１】特許第３７５２２３７号公報
【００１１】
この特許文献１の技術は、信号周波数に対して各Ａ／Ｄ変換器のサンプリング周期が長く、そのサンプル値だけで等価処理を行うことはできないので、各Ａ／Ｄ変換器の見かけ上のサンプリング周期をあげて等価処理を行うものであり、各Ａ／Ｄ変換器のサンプル値と次のサンプル値の間の値を、他のＡ／Ｄ変換器のサンプル値を用いて推定し、その推定値とサンプル値とを時系列にイコライザへ与えることで、各信号経路の差を補償するものである。
【発明の開示】
【発明が解決しようとする課題】
【００１２】
しかしながら、上記構成のＡ／Ｄ変換装置においても、各信号経路の周波数特性のばらつきが大きい場合等には十分な精度が得られない場合があり、より高精度な処理が望まれていた。
【００１３】
本発明は、この課題を解決するためになされたものであり、さらに高精度化したＡ／Ｄ変換装置を提供することを目的としている。
【課題を解決するための手段】
【００１４】
前記目的を達成するために、本発明のＡ／Ｄ変換装置は、
信号を入力するための入力端子（１０ａ）と、
複数のＡ／Ｄ変換器（１２）と、
前記入力端子からの入力信号を前記Ａ／Ｄ変換器にそれぞれ入力する信号分配器（１１）と、
前記各Ａ／Ｄ変換器に対し、サンプリングのための所定周期のクロックを、前記Ａ／Ｄ変換器数で前記所定周期を割って得られる時間にほぼ等しい時間差で所定順に且つ循環的に与えるサンプリング制御部（２１）と、
前記複数のＡ／Ｄ変換器の１つを基準とし、前記入力端子から前記各Ａ／Ｄ変換器の出力端子までのそれぞれの周波数特性と前記入力端子から前記基準のＡ／Ｄ変換器の出力端子までの周波数特性との差の特性を満たすインパルス応答を有するフィルタの係数を予め記憶しているＡＤ特性テーブル（２５）と、
前記Ａ／Ｄ変換器毎にそれぞれ設けられ、前記複数のＡ／Ｄ変換器によって変換出力されるサンプル値および前記ＡＤ特性テーブルに記憶されている係数に基づいて、前記クロックを受けたＡ／Ｄ変換器がサンプル値を更新するタイミングに他のＡ／Ｄ変換器が変換処理をおこなったと仮定して得られるサンプル値をそれぞれ推定する複数の第１の推定手段（２２）と、
前記基準のＡ／Ｄ変換器と各Ａ／Ｄ変換器の周波数特性の差をそれぞれ相殺する周波数特性を満たすインパルス応答を有するフィルタの係数を予め記憶しているイコライザ係数テーブル（２７）と、
前記各第１の推定手段にそれぞれ対応して設けられ、各第１の推定手段の出力値に対して、前記イコライザ係数テーブルに記憶されている係数に基づくフィルタリングをそれぞれ行って、誤差補正された値をそれぞれ出力する複数の第１のイコライザ（２３）と、
前記複数の第１のイコライザの出力値を受け、前記各Ａ／Ｄ変換器が前記クロックによってサンプリングする順に前記複数の第１のイコライザの出力値を選択して出力する第１の信号切換器（２４）と、
前記各第１のイコライザに対応して設けられ、前記ＡＤ特性テーブルに記憶されている係数に基づいて、前記第１の信号切換器から出力された信号に対して前記対応する各イコライザと逆の周波数補正処理を行うとともに、該補正処理で得られたデータのうち、前記対応するＡ／Ｄ変換器のサンプリングタイミングのデータを元のサンプル値に入れ換えて出力する複数の第２の推定手段（３０）と、
前記各第２の推定手段の出力値に対し、前記入力端子から前記各Ａ／Ｄ変換器までの信号経路の周波数特性の差異を補正するための処理を前記イコライザ係数テーブルに記憶された係数に基づいて行う複数の第２のイコライザ（３１）と、
前記各第２のイコライザの出力値を受け、前記各Ａ／Ｄ変換器が前記クロックによってサンプリングする順に前記複数の第２のイコライザの出力値を選択して出力する第２の信号選択手段（３２）とを備えている。
【発明の効果】
【００１５】
このように、本発明では、第１の信号切換器から出力された信号に対して第２の推定手段による推定処理と、第２のイコライザによる補正処理とを行い、その第２のイコライザの出力値を順次選択するようにしたので、より精度の高いＡ／Ｄ変換処理が可能となる。
【発明を実施するための最良の形態】
【００１６】
先ず、本発明の前提となる技術について説明する。
始めに、前記したＮ個のＡ／Ｄ変換器１２のうちの任意のものを基準Ａ／Ｄ変換器と定め、各Ａ／Ｄ変換器毎に、入力端子からＡ／Ｄ変換器までの入力特性や変換特性およびサンプリング系の位相誤差特性をまとめて周波数特性を算出し、その各周波数特性と基準のＡ／Ｄ変換器についての周波数特性との差を求めて、これをミスマッチ特性と定義する。
【００１７】
また、本発明で扱う入力信号ｘ（ｔ）は、Ｎ個のＡ／Ｄ変換器を用いて実現する高速サンプリング周波数をＦｓ（＝１／ΔＴ）としたとき、０〜Ｆｓ／２で帯域制限されているとする。
【００１８】
次に、各ミスマッチ特性を有するミスマッチ回路をそれぞれのＡ／Ｄ変換器の前段に挿入し、その周波数特性をＨｉ（ω）（ｉ＝０，１，…，Ｎ−１）と定義し、さらに、各ミスマッチ特性Ｈｉ（ω）をキャンセルする仮想等価器のイコライズ特性Ｇｉ（ω）を定義する。
【００１９】
ここで、入出力信号が０〜Ｆｓの周波数範囲に帯域制限されている条件下で、連続システムをサンプリング周期ΔＴ（＝１／Ｆｓ）で表される離散システムに置き換えた場合に、ミスマッチ特性Ｈｉ（ω）およびイコライズ特性Ｇｉ（ω）と等価な入出力特性を示すミスマッチ特性Ｈｉ＊（ω）およびイコライズ特性Ｇｉ＊（ω）を考え、これらの特性に対応するインパルス応答ｈｉ，ｕおよびｇｉ，ｋを次式によって算出する。なお、インパルス列の長さｕおよびｋは、必要精度で加減する。
【００２０】
Ｇｉ＊（ω）＝１／Ｈｉ＊（ω）      ……（１）
ｈｉ，ｕ＝Ｆ−１｛Ｈｉ＊（ω）｝      ……（２）
ｇｉ，ｋ＝Ｆ−１｛Ｇｉ＊（ω）｝      ……（３）
ただし、ｉ＝０，１，…，Ｎ
記号Ｆ−１は、離散フーリエ逆変換演算を示す
【００２１】
ここで、Ａ／Ｄ変換器１２０を基準として、図１の等価回路について考察する。
【００２２】
各Ａ／Ｄ変換器１２１〜１２Ｎ−１は、基準のＡ／Ｄ変換器１２０に対するミスマッチ成分がミスマッチ回路特性に換算されているので、図１の等価回路に示すように、入力信号ｘ（ｔ）を基準のＡ／Ｄ変換器１２０の変換特性１１０で離散システムに変換した信号ｘ（ｎ）を、各Ａ／Ｄ変換器についてのミスマッチ回路１１２０〜１１２Ｎ−１に通過させた後に、誤差が無い理想Ａ／Ｄ変換器１３００〜１３０Ｎ−１でＡ／Ｄ変換した場合と等価である。
【００２３】
さらに、各理想Ａ／Ｄ変換器１３００〜１３０Ｎ−１から順次出力されるデジタル値は、それぞれ仮想等価器１３１０〜１３１Ｎ−１に入力され、個々のＡ／Ｄ変換器毎に定義されたイコライザ（インパルス応答ｇｉ，ｋで定義される）で等価処理を実施した後、各仮想等価器１３１０〜１３１Ｎ−１からサンプル値Ｙ（ｎ）として出力されることになる。
【００２４】
なお、以下では説明を簡単化するために、基準のＡ／Ｄ変換特性１１０は、入力信号をそのまま出力に伝送しているものとするが、必要に応じて、この特性を補正してもよい。
【００２５】
上記等価回路において、各ミスマッチ回路１１２０〜１１２Ｎ−１の周波数特性を表すインパルス列の長さｕを等しくＵで表せば、理想Ａ／Ｄ変換器１３００〜１３０Ｎ−１の入力ｘｉ，ｎは、次式で表される。
【００２６】
ｘｉ，ｎ＝ｕΣｘ（ｎ−ｕ）・ｈｉ，ｕ      ……（４）
ただし、ｉ＝０，１，…，Ｎ−１
記号ｕΣは、ｕ＝−（Ｕ−１）〜（Ｕ−１）までの総和を示す
【００２７】
ここで、各Ａ／Ｄ変換器１２０〜１２Ｎ−１のサンプリングタイミングと理想Ａ／Ｄ変換器１３００〜１３０Ｎ−１のサンプリングタイミングを等しくすれば、理想Ａ／Ｄ変換器１３００〜１３０Ｎ−１は、入力された値ｘｉ，ｎを周期ＴでＡ／Ｄ変換処理した後、各Ａ／Ｄ変換器のサンプリングタイミングに合わせてサンプル値を仮想等価器１３１０〜１３１Ｎ−１に出力するから、理想Ａ／Ｄ変換器１３００がＰ番目のサンプル値を出力するとすれば、ｎ番目に出力されるサンプル値は次式で表されるＪ（ｎ）番目の理想Ａ／Ｄ変換器から出力されることになる。
【００２８】
ｘＪ（ｎ），ｎ＝ｕΣｘ（ｎ−ｕ）・ｈＪ（ｎ），ｕ……    （５）
記号ｕΣは、ｕ＝−（Ｕ−１）〜（Ｕ−１）までの総和を示す
【００２９】
ここで、Ｊ（ｎ）は、Ｎを法とする正の値であり、
Ｊ（ｎ）＝ｎ−Ｐ  ｍｏｄ（Ｎ）
と表す。
【００３０】
即ち、個々の理想Ａ／Ｄ変換器は、入力された値ｘｉ，ｎに対して、Ｎ個おき（周期Ｔ秒毎）にデータを仮想等価器に出力することになる。
【００３１】
今、仮に理想Ａ／Ｄ変換器がΔＴ毎にサンプル値を出力することにすれば、ミスマッチ回路から出力される値ｘｉ，ｎが、仮想等価器にそのまま入力されることになり、仮想等価器内部の対応するイコライザは、定義によりミスマッチ回路の特性を補正するように働くから、ミスマッチ回路およびイコライザの計算上の遅延が０となるように係数を定めれば、入力した値ｘ（ｎ）と同じ値のサンプル値Ｙ（ｎ）がＮ個の仮想等価器１３１０〜１３１Ｎ−１から出力されることになる。
【００３２】
理想Ａ／Ｄ変換器がΔＴ毎にサンプル値を出力したと仮定したときに、仮想等価器１３１０〜１３１Ｎ−１内部のイコライザによる処理は、対応するＡ／Ｄ変換器毎に定められるイコライザのインパルス応答ｇｉ，ｋを用いて次式で定められる。
【００３３】
Ｙ（ｎ）＝ｋΣｘＪ（ｎ），ｎ−ｋ・ｇＪ（ｎ），ｋ    ……（６）
ただし、Ｋはイコライザのインパルス列の長さを示し、記号ｋΣは、ｋ＝−（Ｋ−１）〜Ｋ−１までの総和を示す
【００３４】
ここで上式（６）が成立するためには、ｘＪ（ｎ），ｎ−ｋについて、ｋ＝−（Ｋ−１）〜Ｋ−１に対して全ての値が必要であるが、実際の各Ａ／Ｄ変換器は、前記したように、Ｎ個おきの値しか出力できない。
【００３５】
そこで、他のＡ／Ｄ変換器のサンプル値を用いて、イコライズに必要なサンプル値を推定し、その後に式（６）の等価演算処理を行う。
【００３６】
さらに、各仮想等価器１３１０〜１３１Ｎ−１が算出したｎ番目の出力候補のうち、最も誤差が少なくなるＪ（ｎ）番目（演算による遅延を０とした場合）の仮想等価器からの出力をサンプル値Ｙ（ｎ）として出力する。
【００３７】
ここで、Ｊ（ｎ）番目のＡ／Ｄ変換結果を推定するために、Ｊ（ｎ）番目以外のＡ／Ｄ変換出力
ｘＪ（ｎ−ｒ），ｎ−ｒ−ｋ
ただし、ｒ≠ｑ×Ｎ（ｑ：０，±１，±２，…）
の場合について考察する。
【００３８】
この場合、ｎ−ｒ番目の値をもっているのは、（ｎ−ｒ−Ｐ）  ｍｏｄ（Ｎ）番目のＡ／Ｄ変換器であり、定義によりｎ−ｒ番目の入力値ｘ（ｎ−ｒ）は、イコライズされた出力値Ｙ（ｎ−ｒ）と等しい値であるから、次式が成り立つ。
【００３９】
ｘ（ｎ−ｒ）＝Ｙ（ｎ−ｒ）
＝ｋΣｘＪ（ｎ−ｒ），ｎ−ｒ−ｋ・ｇＪ（ｎ−ｒ），ｋ    ……（７）
ただし、記号ｋΣは、ｋ＝−（Ｋ−１）〜Ｋ−１までの総和を示す
【００４０】
また、式（４）において、理想Ａ／Ｄ変換器がサンプリングタイミングをずらし、Ｊ（ｎ）番目のＡ／Ｄ変換器がｎ−ｒ番目のサンプリングを行ったと仮定して得られる推定サンプル値ｘＪ（ｎ），ｎ−ｒは、以下のように得られる。
【００４１】
ｘＪ（ｎ），ｎ−ｒ＝ｕΣｘ（ｎ−ｒ−ｕ）・ｈＪ（ｎ），ｕ  ……（８）
ただし、記号ｕΣは、ｕ＝−（Ｕ−１）〜Ｕ−１までの総和を示す
【００４２】
上記式（８）に式（７）を代入すれば、推定サンプル値ｘＪ（ｎ），ｎ−ｒが得られ、その得られた推定サンプル値に対して前記式（６）の処理を行うことで、Ｎ個のＡ／Ｄ変換器による出力値ｙ（ｎ）を得ることができる。
【００４３】
以下、図面に基づいて本発明の実施の形態を説明する。
図２は、上記前提技術に基づく実施形態のインタリーブ方式のＡ／Ｄ変換装置２０の構成を示している。
【００４４】
図２において、入力端子１０ａ、信号分配器１１、Ｎ個のＡ／Ｄ変換器１２０〜１２Ｎ−１および出力端子１０ｂは、前記したＡ／Ｄ変換装置１０と同一であるので同一符号を付して説明する。
【００４５】
このＡ／Ｄ変換装置２０では、前記したＡ／Ｄ変換器１０と同様に、入力端子１０ａに入力されるアナログの入力信号ｘ（ｔ）が、信号分配器１１によって複数Ｎ本の信号経路に分岐されて、周波数特性がほぼ等しいＮ個の信号ｘ０（ｔ）〜ｘＮ−１（ｔ）がＡ／Ｄ変換器１２０〜１２Ｎ−１にそれぞれ入力する。
【００４６】
また、サンプリング制御部２１から、周期ＴでΔＴ（＝Ｔ／Ｎ）時間ずつ位相がシフトしたサンプリング用のクロックＣ０〜ＣＮ−１を発生してそれぞれＡ／Ｄ変換器１２０〜１２Ｎ−１に与えて、入力信号に対するサンプリングを各Ａ／Ｄ変換器１２０〜１２Ｎ−１で行わせる。
【００４７】
このサンプリング制御部２１は、周期ΔＴのサンプリングタイミング信号（以下、単にタイミング信号と記す）Ｃｔを基に前記したクロックＣ０〜ＣＮ−１を生成するとともに、タイミング信号Ｃｔのタイミングに合わせて、Ａ／Ｄ変換結果（サンプル値）を更新するＡ／Ｄ変換器を指定する指定信号ＡＤＮＵＭを、後述する第１の推定手段２２０〜２２Ｎ−１、第１の信号切換器２４、第２の推定手段３００〜３０Ｎ−１および第２の信号切換器３２に出力する。
【００４８】
各Ａ／Ｄ変換器１２０〜１２Ｎ−１の出力は、それぞれＮ個の第１の推定手段２２０〜２２Ｎ−１に入力される。
【００４９】
各第１の推定手段２２０〜２２Ｎ−１は、それぞれがＮ個のＡ／Ｄ変換器１２０〜１２Ｎ−１の出力と、サンプリング制御部２１からの指定信号ＡＤＮＵＭを受けている。
【００５０】
各第１の推定手段２２０〜２２Ｎ−１は、タイミング信号Ｃｔで示されるタイミング毎に、入力されたＮ個のサンプル値、指定信号ＡＤＮＵＭおよび後述するＡＤ特性テーブル２５の係数とに基づいて、予め決定した推定値算出処理により定まる数Ｅ（３点のサンプリング点を用いて推定値を得る場合にＥ＝１以上となり、１点のサンプリング点を用いて推定値を得る場合にはＥ＝０以上となる）個前のサンプリングタイミングで、Ａ／Ｄ変換器がサンプリング動作したと仮定したときのサンプル値を推定する。
【００５１】
例えば、３点のサンプリング点を用いて推定を行う場合には、更新されたサンプル値をもつＡ／Ｄ変換器の番号をａ（ＡＤＮＵＭ＝ａ）とし、Ｎを法とする正の数ｂ、ｃを次式によって求める。
【００５２】
ｂ＝ａ−１  ｍｏｄ（Ｎ）      ……（１１ａ）
ｃ＝ａ−２  ｍｏｄ（Ｎ）      ……（１１ｂ）
【００５３】
そして、ｉ＝ｂのとき、推定サンプル値Ｗｉ，ｎを、
Ｗｉ，ｎ＝ｘｂ，ｎ        ……（１２ａ）
とする。
【００５４】
また、ｉ≠ｂのとき、推定サンプル値Ｗｉ，ｎを、次の演算で求める。
【００５５】
Ｗｉ，ｎ＝ｘｂ，ｎ・ｈｉ，０／ｈｂ，０
＋ｘａ，ｎ・（ｈｉ，０／ｈａ，０）
・｛（ｈｉ，−１／ｈｉ，０）−（ｈｂ，−１／ｈｂ，０）｝
＋ｘｃ，ｎ・（ｈｉ，０／ｈｃ，０）
・｛（ｈｉ，１／ｈｉ，０）−（ｈｂ，１／ｈｂ，０）｝
……（１２ｂ）
【００５６】
上記式で、ｈｉ，−１、ｈｉ，０、ｈｉ，１は、後述するＡＤ特性テーブル２５に予め記憶されている係数である。また、上記式（１２ｂ）の第１項は主に振幅誤差に関わる項、第２項および第３項は主に位相誤差に関わる項である。
【００５７】
各第１の推定手段２２から出力された推定サンプル値Ｗは、それぞれ第１のイコライザ２３０〜２３Ｎ−１に入力される。
【００５８】
各第１のイコライザ２３０〜２３Ｎ−１は、入力された推定サンプル値Ｗに対して、後述するイコライザ係数テーブル２７に記憶されている係数（フィルタ係数）を用いて等価演算処理を行って、その結果、即ち、基準のＡ／Ｄ変換器に対して誤差補正されたサンプル値ｙをタイミング信号Ｃｔで示されるタイミングでそれぞれ第１の信号切換器２４に出力する。
【００５９】
第１の信号切換器２４は、各イコライザ２３０〜２３Ｎ−１から出力されるサンプル値を受け、指定信号ＡＤＮＵＭで指定された値（ここではＡＤＮＵＭ＝ａ）、推定値算出処理によって定まる数Ｅおよびイコライザ係数テーブル２７を定義する際に定められるオフセット値ａ０を用いてイコライザを指定する値ｅを、
ｅ＝ａ−Ｅ−ａ０  ｍｏｄ（Ｎ）
の計算により求め、指定信号ＡＤＮＵＭで指定された値ａに対してｅ番目のイコライザ２３ｅの出力結果ｙｅ，ｎを選択して、最終のＡＤ変換結果Ｙ（ｎ）として出力する。
【００６０】
なお、得られるＡ／Ｄ変換結果は、推定値算出処理により理論計算よりＥ＋ａ０分のサンプリングタイミングだけ遅延して得られる。
【００６１】
一方、ＡＤ特性テーブル２５には、サンプリング周期ΔＴ（＝Ｔ／Ｎ）で表される離散システムで考慮した場合の入力端子１０ａから各Ａ／Ｄ変換器の出力端までの周波数特性に対する基準のＡ／Ｄ変換器との周波数特性の差Ｈｉ＊（ω）に対応したインパルス応答によって決まるイコライザ係数が必要なポイント数予め記憶されている。このイコライザ係数のポイント数は、上記第１の推定手段２２で用いる３ポイントだけでなく、後述する第２の推定手段３０で用いるそれ以上の分も含まれている。
【００６２】
このインパルス応答を求めるために、周波数特性の差の特性Ｈｉ＊（ω）を、基準Ａ／Ｄ変換器についての周波数特性ＨＯ＊（ω）および各Ａ／Ｄ変換器１２０〜１２Ｎ−１についての周波数特性ＨＯｉ＊（ω）から次式によって算出する。なお、差の特性は計算上では以下のように比となる。
【００６３】
Ｈｉ＊（ω）＝ＨＯｉ＊（ω）／ＨＯ＊（ω）……（１３）
【００６４】
次に、サンプリング定理を満たす範囲において、周波数特性Ｈｉ＊（ω）と等価なインパルス応答をもつＦＩＲフィルタを設計する。ただし、前記等価なインパルス応答をもつフィルタの設計に際しては、設計されるＮ個のフィルタ全てに共通する絶対遅延量τ０（秒）を任意に設定した後に、個々のフィルタ設計を行う。
【００６５】
得られるフィルタの係数を時系列順に、…、ｈｉ，−１、ｈｉ，０、ｈｉ，１、…（ただし、ｉ＝０，１，２，…，Ｎ−１）と表した場合、絶対遅延量τ０は、係数ｈｉ，０の絶対値が最大となり、かつ設計するＮ個のＦＩＲフィルタの係数を考慮した場合に、係数の２乗の総和Σ（ｈｉ，−１）２とΣ（ｈｉ，１）２がほぼ等しい値となるように絶対遅延量τ０を設定する。
【００６６】
次に、得られた係数の中から、｜ｈｉ，Ｕ１｜＜ε（ここでεは、予め定められた許容誤差）を満足する最小値Ｕ１を決定し、同様にして｜ｈｉ，Ｕ２｜＜εを満足する最大値Ｕ２を決定し、係数列ｈｉ，Ｕ１、…、ｈｉ，−１、ｈｉ，０、ｈｉ，１、…、ｈｉ，Ｕ２を用いて、図３に示すＡＤ特性テーブル２５を作成する。
【００６７】
このＡＤ特性テーブル２５は、例えばテーブル位置（ｉ，−１）にはｈｉ，−１を、テーブル位置（ｉ，０）にはｈｉ，０を、テーブル位置（ｉ，１）にはｈｉ，１を対応させる。
【００６８】
一方、イコライザ係数テーブル２７は、前記した式（１３）で算出した周波数特性の差Ｈｉ＊（ω）を基に、次式により周波数特性Ｇｉ＊（ω）を算出する。
【００６９】
Ｇｉ＊（ω）＝１／Ｈｉ＊（ω）      ……（１４）
ただし、Ｈｉ＊（ω）≠０
【００７０】
そして、サンプリング定理を満たす範囲では、周波数特性Ｇｉ＊（ω）と等価なインパルス応答をもつイコライザ（フィルタ）をｉ番目のＡ／Ｄ変換器に対応するイコライザと定義し、そのイコライザに要求されるフィルタ係数を求めてイコライザ係数テーブル２７に予め用意しておく。ただし、この等価なインパルス応答をもつフィルタの設計に際しては、設計されるＮ個のフィルタ全てに共通する絶対遅τ１（秒）を設定した後に、個々のフィルタ設計を行う。
【００７１】
得られるフィルタの係数を時系列順に、…、ｇｉ，−１、ｇｉ，０、ｇｉ，１、…と表した場合、全フィルタに共通する絶対遅延量τ１の設定値は任意であるが、イコライザ係数テーブル２７の設計においては、係数ｇｉ，０の絶対値が最大となり、かつ設計するＮ個のフィルタ係数の２乗の総和Σ（ｇｉ，−１）２とΣ（ｇｉ，１）２がほぼ等しくなるように絶対遅延量τ１を設定する。
【００７２】
次に、得られた係数の中から、｜ｇｉ，Ｍ１｜＜ε（ここでεは、予め定められた許容誤差）を満足する最小値Ｍ１を決定し、同様にして｜ｇｉ，Ｍ２｜＜εを満足する最大値Ｍ２を決定し、係数列ｇｉ，Ｍ１、…、ｇｉ，−１、ｇｉ，０、ｇｉ，１、…、ｇｉ，Ｍ２を用いて、図４のように、イコライザ係数テーブル２７を作成する。この場合、例えばテーブル位置（ｉ，Ｍ１）にはｇｉ，Ｍ１を、テーブル位置（ｉ，Ｍ１＋１）にはｈｉ，Ｍ１＋１を対応させ、以後同様に、テーブル位置（ｉ，Ｍ２）まで順に対応させる。
【００７３】
このとき、設計される第１の推定手段２２、第１のイコライザ２３の時間応答に合わせて、前記したオフセット値ａ０＝１（構成する回路の絶対遅延量により異なる）を決定する。
【００７４】
次に、このＡ／Ｄ変換装置２０の上記第１の信号切替器２４までの動作を図５、図６に基づいて説明する。
【００７５】
図５の（ａ）のように入力端子１０ａに入力された入力信号ｘ（ｔ）は、信号分配器１１によってＮ本の信号経路に分岐され，各Ａ／Ｄ変換器１２０〜１２Ｎ−１に入力される。
【００７６】
各Ａ／Ｄ変換器１２０〜１２Ｎ−１は、図５の（ｂ１）〜（ｂＮ）のように、サンプリング制御部２１から出力されるクロックＣ０〜Ｃ１２をそれぞれ受けて、それぞれの入力信号ｘ０（ｔ）〜ｘＮ−１（ｔ）に対するＡ／Ｄ変換処理をほぼΔＴ時間ずつ遅れたタイミングに順次行い、その変換処理によって得られたサンプル値Ｘ０，Ｐ、Ｘ１，Ｐ＋１、…、ＸＮ−１，Ｐ＋Ｎ−１を、図５の（ｃ１）〜（ｃＮ）のようにそれぞれ出力する。
【００７７】
ここで、サンプリングタイミング順に番号を付け、Ｐ番目のサンプリングでは、Ａ／Ｄ変換器１２０がＡ／Ｄ変換処理を行ってそのサンプル値を更新したと定義し、その更新されたサンプル値をＸ０，Ｐ、と表すとする。
【００７８】
このとき、サンプリング制御部２１は図５の（ｄ）、（ｅ）に示すように、Ａ／Ｄ変換結果の更新タイミングに合わせて、サンプル値を更新したＡ／Ｄ変換器１２０を指定する指定信号ＡＤＮＵＭ（例えばＡＤＮＵＭ＝０とする）と、入力信号に対するサンプリングタイミングを示すタイミング信号Ｃｔを出力する。
【００７９】
他のＡ／Ｄ変換器１２１〜１２Ｎ−１は変換結果を更新しないので、Ｐ番目のサンプリングが行われる前から保持している値を出力している。
【００８０】
即ち、
Ｘ１，Ｐ＝Ｘ１，Ｐ−１、Ｘ２，Ｐ＝Ｘ２，Ｐ−１、…、
ＸＮ−１，Ｐ＝ＸＮ−１，Ｐ−１
となる。
【００８１】
次のＰ＋１番目のサンプリングタイミングには、ＡＤＮＵＭ＝１となり、Ａ／Ｄ変換器１２１のサンプル値が更新され、他のＡ／Ｄ変換器１２０、１２２〜１２Ｎ−１は、Ｐ番目のサンプリングタイミングのときと同じ値を出力する。
【００８２】
以後同様に各Ａ／Ｄ変換器１２０〜１２Ｎ−１による変換処理が順番に行われ、Ｎ−１番目のＡ／Ｄ変換器１２Ｎ−１のサンプル値が更新された後に、再び０番目のＡ／Ｄ変換器１２０によるサンプル値の更新がなされ、上記動作が循環的に繰り返される。
【００８３】
各第１の推定手段２２０〜２２Ｎ−１は、前記したように、サンプル値が更新されていないＡ／Ｄ変換器がそのタイミングでサンプリング動作したと仮定したときのサンプル値を、更新されたサンプル値を用いて推定する。
【００８４】
例えば、Ｎが３以上の場合で、一つの第１の推定手段２２０についてみると、図６に示すように、Ａ／Ｄ変換器１２１によりＰ＋１番目のサンプル値が更新されたタイミングでは、各Ａ／Ｄ変換器について一つ前のサンプリングタイミングでＰ番目のサンプル値の推定が可能となる。第１の推定手段２２０のＰ番目の推定サンプル値Ｗ０，Ｐとしては、Ａ／Ｄ変換器１２０がサンプル値Ｘ０，Ｐを既にもっているから、この値をそのまま出力する。即ち、前記式（１２ａ）のｉ＝ｂ＝０の場合に相当する。
【００８５】
また、その次のＰ＋２番目のサンプリングタイミングにおける推定サンプル値Ｗ０，Ｐ＋１は、そのサンプリングタイミングに更新されたＡ／Ｄ変換器１２２のサンプル値Ｘ２，Ｐ＋２と、一つ前のサンプリングタイミングのサンプル値Ｘ１，Ｐ＋１と、さらにその一つ前のサンプリングタイミングのサンプル値Ｘ０，Ｐと、ＡＤ特性テーブル２５の係数とを用いて、前記式（１２ｂ）のｉ≠ｂの場合で示す演算式にしたがって算出する。
【００８６】
さらに、その次のＰ＋３番目のサンプリングタイミングにおける推定サンプル値Ｗ０，Ｐ＋２は、そのサンプリングタイミングに更新されたＡ／Ｄ変換器１２３のサンプル値Ｘ３，Ｐ＋３と、一つ前のサンプリングタイミングのサンプル値Ｘ２，Ｐ＋２と、さらにその一つ前のサンプリングタイミングのサンプル値Ｘ１，Ｐ＋１と、ＡＤ特性テーブル２５の係数とを用いて、前記式（１２ｂ）のｉ≠ｂの場合で示す演算式にしたがって算出する。
【００８７】
以下同様の推定処理がなされて、その推定サンプル値が時系列に並んだサンプル列Ｗ０，Ｐ、Ｗ０，Ｐ＋１、…が第１のイコライザ２３０に出力される。
【００８８】
他の第１の推定手段２２１〜２２Ｎ−１についても同様の推定処理がなされ、その推定サンプル値Ｗｍ，Ｐ、Ｗｍ，Ｐ＋１、…（ｍ＝１，２，…，Ｎ−１）がそれぞれ第１のイコライザ２３１〜２３Ｎ−１に出力される。
【００８９】
第１のイコライザ２３０〜２３Ｎ−１は、それぞれ入力される推定サンプル値Ｗに対して、イコライザ係数テーブル２７の係数による等価処理（フィルタリング）を行い、基準のＡ／Ｄ変換器について周波数特性に対して誤差補正されたサンプル値ｙｉ，Ｐ、ｙｉ，Ｐ＋１、…（ｉ＝０，１，…，Ｎ−１）を第１の信号切換器２４にそれぞれ出力する。
【００９０】
第１の信号切換器２４は、Ａ／Ｄ変換器を指定する指定信号ＡＤＮＵＭに対して前記したオフセット値ｅ分だけずれたタイミングにその指定信号ＡＤＮＵＭで指定されるＡ／Ｄ変換器に対応する第１のイコライザ２３の出力値を順次選択して、その選択値が時系列に並んだデジタル信号列Ｙ（ｎ）を出力する。
【００９１】
このようにして得られたＡ／Ｄ変換結果Ｙ（ｎ）は、各Ａ／Ｄ変換器１２０〜１２Ｎ−１の変換処理で実際に得られたサンプル値と各第１の推定手段２２で推定算出されたサンプル値からなるサンプル列を、それぞれ第１のイコライザ２３０〜２３Ｎ−１によって誤差補正しているため、信号分配器１１や配線等を含むＡ／Ｄ変換器間の周波数特性差による誤差の影響を格段に低減することができる。
【００９２】
また、各第１のイコライザ２３０〜２３Ｎ−１が出力するサンプル列のうち、同一サンプリングタイミングで得られる最も誤差の少ないサンプル値が第１の信号切換器２４によって選択されるようにしているので、時間波形解析や周波数スペクトラムによる解析誤差を大幅に改善することができる。
【００９３】
ただし、前記したように、各信号経路の周波数特性のばらつきが大きいような場合、十分な精度が得られない場合があり、より高精度な処理が望まれていた。
【００９４】
図７の（ａ）は、８個のＡ／Ｄ変換器をもつ従来のインタリーブ方式のＡ／Ｄ変換装置でサンプリングして得られたデータから生成したスペクトラムであり、横軸は装置全体としてのサンプリング周波数ｆｓの１／２で正規化した周波数であり、約０．１９の周波数の正弦波信号を与えたときの結果である。
【００９５】
これに対し、上記した第１の推定手段２２０〜２２Ｎ−１、第１のイコライザ２３０〜２３Ｎ−１および第１の信号切換器２４を用いた構成では、図７の（ｂ）のように、スプリアスを大幅に軽減することができているが、ＳＦＤＲ（スプリアスフリーダイナミックレンジ）は−７０ｄＢ程度にとどまっている。
【００９６】
そこでこの実施形態では、図２に示しているように、複数の第２の推定手段３００〜３０Ｎ−１、複数の第２のイコライザ３１０〜３１Ｎ−１および第２の信号切換器３２を設けて、より精度の高いＡ／Ｄ変換装置を実現している。
【００９７】
ここで、第２の推定手段３００〜３０Ｎ−１は、図８に示しているように、それぞれ各第１のイコライザ２３０〜２３Ｎ−１に対応して設けられ、ＡＤ特性テーブル２５に記憶されている係数ｈｉ，Ｕ１、…、ｈｉ，−１、ｈｉ，０、ｈｉ，１、…、ｈｉ，Ｕ２に基づいて、第１の信号切換器２４から出力された信号列Ｙ（ｎ）に対して対応する各第１のイコライザ２３０〜２３Ｎ−１と逆の周波数補正処理を行う補正処理部３０ａと、その補正処理で得られたデータ列Ｑのうち、対応するＡ／Ｄ変換器のサンプリングタイミングのデータＱ（ｋ）を指定信号ＡＤＮＵＭに基づいて元のサンプル値Ｘ（ｋ）に入れ換えるデータ置換部３０ｂとを有している。このような処理を行うことで、第１の推定手段２２０〜２２Ｎ−１による推定結果よりも精度の高い（誤差の少ない）、推定値を得ることができている。
【００９８】
そして、この誤差のより少ない推定値Ｗ′に対し、前記した第１のイコライザ２３０〜２３Ｎ−１と同様に構成された第２のイコライザ３１０〜３１Ｎ−１によって、入力端子１０ａから各Ａ／Ｄ変換器１２０〜１２Ｎ−１までの信号経路の周波数特性の差異を補正するための処理をイコライザ係数テーブル２７に記憶された係数に基づいて行い、それらのイコライザによる補正処理で得られた出力値ｚ０〜ｚＮ−１を第２の信号切換器３２に与え、各Ａ／Ｄ変換器１２０〜１２Ｎ−１がクロックによってサンプリングする順に複数の第２のイコライザ３１０〜３１Ｎ−１の出力値ｚ０〜ｚＮ−１を選択し、出力する。
【００９９】
このように構成した場合に得られたデジタル信号列Ｚ（ｎ）のスプリアス特性は、図７の（ｃ）となり、第２の推定手段３００〜３０Ｎ−１、第２のイコライザ３１０〜３１Ｎ−１および第２の信号切換器３２を用いたことにより、ＳＦＤＲを−９０ｄＢ程度まで改善することができた。
【０１００】
図９、図１０は、一定振幅の正弦波信号に対する一つのＡ／Ｄ変換器のサンプリング波形と各部の誤差を求めたシミュレーション結果である。
【０１０１】
図９の（ａ）はサンプリング値、図９の（ｂ）は第１の推定手段２２の出力誤差、図９の（ｃ）は第１のイコライザ２３の出力誤差、そして図９の（ｄ）は第１の信号切換器２４の出力誤差であり、この時点での誤差の最大値はほぼ１０−３Ｖ（１ｍＶ）である。
【０１０２】
また、図１０の（ａ）は第２の推定手段３０の出力誤差、図１０の（ｂ）は第２のイコライザ３１の出力誤差、そして図１０の（ｃ）は、第２の信号切換器３２の出力誤差、即ち本実施形態のＡ／Ｄ変換装置２０の最終出力誤差であり、その最大値はほぼ１０−４Ｖ（０．１ｍＶ）であり、前記したスペクトラムと同様にほぼ２０ｄＢ改善されていることがわかる。
【図面の簡単な説明】
【０１０３】
【図１】本発明の前提技術を説明するための図
【図２】本発明の実施形態の構成を示す図
【図３】実施形態の要部のテーブル図
【図４】実施形態の要部のテーブル図
【図５】実施形態の動作説明図
【図６】実施形態の動作説明図
【図７】実施形態の特性を示す図
【図８】実施形態の要部の構成図
【図９】実施形態の各部の出力誤差を示す図
【図１０】実施形態の各部の出力誤差を示す図
【図１１】インタリーブ方式の従来装置の基本構成を示す図
【図１２】従来装置の動作説明図
【符号の説明】
【０１０４】
１０ａ……入力端子、１０ｂ……出力端子、１１……信号分配器、１２……Ａ／Ｄ変換器、２０……Ａ／Ｄ変換装置、２１……サンプリング制御部、２２……第１の推定手段、２３……第１のイコライザ、２４……第１の信号切換器、２５……ＡＤ特性テーブル、２７……イコライザ係数テーブル、３０……第２の推定手段、３１……第２のイコライザ、３２……第２の信号切換器
【図１】
2010141721.tif 000003
【図２】
2010141721.tif 000004
【図３】
2010141721.tif 000005
【図４】
2010141721.tif 000006
【図５】
2010141721.tif 000007
【図６】
2010141721.tif 000008
【図７】
2010141721.tif 000009
【図８】
2010141721.tif 000010
【図９】
2010141721.tif 000011
【図１０】
2010141721.tif 000012
【図１１】
2010141721.tif 000013
【図１２】
2010141721.tif 000014