(57)【要約】【課題】複数のＡ／Ｄ変換器を用いたＡ／Ｄ変換において、正確にＡ／Ｄ変換を行う。【解決手段】アナログ信号処理装置1が、複数のＡ／Ｄ変換器14a,14bと、Ａ／Ｄ変換器14a,14bの各々の出力を、サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換するフーリエ変換部16a,16bと、周波数領域の信号X1、X2と、周波数領域の信号の周波数f1と、サンプリングタイミングの各々の時間差δtとに基づき、アナログ信号の測定帯域内の測定周波数成分S1と、周波数領域の信号において測定周波数成分と同一の周波数を有する重畳成分R12の基となる、アナログ信号の測定帯域外の重畳要因成分S2とを導出する成分導出部18,19と、測定周波数成分S1と、重畳要因成分S2とを有するスペクトルを導出するスペクトル導出部20とを備える。【選択図】図１
2010283684.tif 000002
【特許請求の範囲】
【請求項１】
同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器と、
前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換部と、
前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出部と、
前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出部と、
を備えたアナログ信号処理装置。
【請求項２】
請求項１に記載のアナログ信号処理装置であって、
前記重畳成分は、前記重畳要因成分のレプリカである、
アナログ信号処理装置。
【請求項３】
請求項２に記載のアナログ信号処理装置であって、
前記成分導出部は、
前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差と、前記重畳要因成分の周波数とに基づき、前記測定周波数成分と、前記重畳要因成分とを導出する、
アナログ信号処理装置。
【請求項４】
請求項１に記載のアナログ信号処理装置であって、
前記成分導出部は、
前記重畳要因成分を前記重畳成分に基づいて導出する、
アナログ信号処理装置。
【請求項５】
請求項４に記載のアナログ信号処理装置であって、
前記成分導出部は、
前記重畳要因成分を、前記重畳成分そのもの、または前記重畳成分の共役複素数として導出する、
アナログ信号処理装置。
【請求項６】
請求項１に記載のアナログ信号処理装置であって、
前記サンプリングタイミングの位相が、２π[rad]を前記Ａ／Ｄ変換器の個数で割った値ずつ異なっているものとして、前記測定帯域が定められている、
アナログ信号処理装置。
【請求項７】
請求項１ないし６のいずれか一項に記載のアナログ信号処理装置であって、
前記スペクトル導出部の出力を時間領域に変換する時間領域変換部、
を備えたアナログ信号処理装置。
【請求項８】
同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器を備えたアナログ信号処理装置によるアナログ信号処理方法であって、
前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換工程と、
前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出工程と、
前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出工程と、
を備えたアナログ信号処理方法。
【請求項９】
同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器を備えたアナログ信号処理装置によるアナログ信号処理をコンピュータに実行させるためのプログラムであって、
前記アナログ信号処理は、
前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換工程と、
前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出工程と、
前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出工程と、
を備えたプログラム。
【請求項１０】
同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器を備えたアナログ信号処理装置によるアナログ信号処理をコンピュータに実行させるためのプログラムを記録したコンピュータによって読み取り可能な記録媒体であって、
前記アナログ信号処理は、
前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換工程と、
前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出工程と、
前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出工程と、
を備えた記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、複数のＡ／Ｄ変換器を用いたＡ／Ｄ変換に関する。
【背景技術】
【０００２】
従来より、複数のＡ／Ｄ変換器を用いた、インターリーブ方式によるＡ／Ｄ変換が知られている（例えば、特許文献１、２、３および４を参照）。インターリーブ方式によるＡ／Ｄ変換においては、各々のＡ／Ｄ変換器のサンプリングタイミングを、他のＡ／Ｄ変換器のサンプリングタイミングと、所定の位相差だけ正確にずらす必要がある。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０００−３４６９１３号公報
【特許文献２】特開２００２−７１７２３号公報
【特許文献３】特開２００２−２１７７３２号公報
【特許文献４】特開２００２−２４６９１０号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
しかしながら、各々のＡ／Ｄ変換器のサンプリングタイミングどうしの位相差を、正確に所定値にすることは、難しい。
【０００５】
そこで、本発明は、複数のＡ／Ｄ変換器を用いたＡ／Ｄ変換において、各々のＡ／Ｄ変換器のサンプリングタイミングどうしの位相差が、正確には所定値になっていない場合でも、正確にＡ／Ｄ変換を行うことを課題とする。
【課題を解決するための手段】
【０００６】
本発明にかかるアナログ信号処理装置は、同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器と、前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換部と、前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出部と、前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出部とを備えるように構成される。
【０００７】
上記のように構成されたアナログ信号処理装置によれば、複数のＡ／Ｄ変換器が、同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する。周波数領域変換部が、前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する。成分導出部が、前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する。スペクトル導出部が、前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出する。
【０００８】
なお、本発明にかかるアナログ信号処理装置は、前記重畳成分は、前記重畳要因成分のレプリカであるようにしてもよい。
【０００９】
なお、本発明にかかるアナログ信号処理装置は、前記成分導出部は、前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差と、前記重畳要因成分の周波数とに基づき、前記測定周波数成分と、前記重畳要因成分とを導出するようにしてもよい。
【００１０】
なお、本発明にかかるアナログ信号処理装置は、前記成分導出部が、前記重畳要因成分を前記重畳成分に基づいて導出するようにしてもよい。
【００１１】
なお、本発明にかかるアナログ信号処理装置は、前記成分導出部が、前記重畳要因成分を、前記重畳成分そのもの、または前記重畳成分の共役複素数として導出するようにしてもよい。
【００１２】
なお、本発明にかかるアナログ信号処理装置は、前記サンプリングタイミングの位相が、２π[rad]を前記Ａ／Ｄ変換器の個数で割った値ずつ異なっているものとして、前記測定帯域が定められているようにしてもよい。
【００１３】
なお、本発明にかかるアナログ信号処理装置は、前記スペクトル導出部の出力を時間領域に変換する時間領域変換部を備えるようにしてもよい。
【００１４】
本発明は、同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器を備えたアナログ信号処理装置によるアナログ信号処理方法であって、前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換工程と、前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出工程と、前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出工程とを備えたアナログ信号処理方法である。
【００１５】
本発明は、同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器を備えたアナログ信号処理装置によるアナログ信号処理をコンピュータに実行させるためのプログラムであって、前記アナログ信号処理は、前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換工程と、前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出工程と、前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出工程とを備えたプログラムである。
【００１６】
本発明は、同じ周波数で位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号をサンプリングしてデジタル信号に変換する複数のＡ／Ｄ変換器を備えたアナログ信号処理装置によるアナログ信号処理をコンピュータに実行させるためのプログラムを記録したコンピュータによって読み取り可能な記録媒体であって、前記アナログ信号処理は、前記Ａ／Ｄ変換器の各々の出力を、前記サンプリングタイミングの周波数の２分の１以下の幅の測定帯域の内の周波数領域の信号に変換する周波数領域変換工程と、前記周波数領域の信号と、前記周波数領域の信号の周波数と、前記サンプリングタイミングの各々の位相差とに基づき、前記アナログ信号の前記測定帯域内の測定周波数成分と、前記周波数領域の信号において前記測定周波数成分と同一の周波数を有する重畳成分の基となる、前記アナログ信号の前記測定帯域外の重畳要因成分とを導出する成分導出工程と、前記測定周波数成分と、前記重畳要因成分とを有するスペクトルを導出するスペクトル導出工程とを備えた記録媒体である。
【図面の簡単な説明】
【００１７】
【図１】本発明の第一の実施形態にかかるアナログ信号処理装置１の構成を示すブロック図である。
【図２】第一の実施形態にかかるアナログ信号の周波数スペクトルである。
【図３】第一の実施形態にかかるＡ／Ｄ変換器１４ａが使用するサンプリングタイミング（図３（ａ））、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミング（図３（ｂ））である。
【図４】第一の実施形態にかかる出力X1およびX2を説明するために、測定帯域Band1内の信号の周波数f1の成分および帯域Band2内の信号の周波数f2の成分を図示した図である。
【図５】第一の実施形態にかかる周波数f2の成分S2が折り返して、周波数f1の成分S1に重なった状態を示す図である。
【図６】第一の実施形態にかかる測定周波数成分S1（図６（ａ））、重畳要因成分S2（図６（ｂ））、スペクトル導出部２０の導出するスペクトルS（図６（ｃ））を示す図である。
【図７】本発明の第二の実施形態にかかるアナログ信号処理装置１の構成を示すブロック図である。
【図８】第二の実施形態にかかるアナログ信号の周波数スペクトルである。
【図９】第二の実施形態にかかるＡ／Ｄ変換器１４ａが使用するサンプリングタイミング（図９（ａ））、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミング（図９（ｂ））、Ａ／Ｄ変換器１４ｃが使用するサンプリングタイミング（図９（ｃ））である。
【図１０】第二の実施形態にかかる出力X1、X2およびX3を説明するために、測定帯域Band1内の信号の周波数f1の成分、帯域Band2内の信号の周波数f2の成分および、帯域Band3内の信号の周波数f3の成分を図示した図である。
【図１１】第二の実施形態にかかる周波数f2の成分S2が折り返し、さらに周波数f3の成分S3が折り返して、周波数f1の成分S1に重なった状態を示す図である。
【図１２】第二の実施形態にかかる測定周波数成分S1（図１２（ａ））、重畳要因成分S2（図１２（ｂ））、重畳要因成分S3（図１２（ｃ））、スペクトル導出部２０の導出するスペクトルS（図１２（ｄ））を示す図である。
【発明を実施するための形態】
【００１８】
以下、本発明の実施形態を図面を参照しながら説明する。
【００１９】
第一の実施形態
図１は、本発明の第一の実施形態にかかるアナログ信号処理装置１の構成を示すブロック図である。第一の実施形態にかかるアナログ信号処理装置１は、バンドパスフィルタ１２、Ａ／Ｄ変換器１４ａ、１４ｂ、フーリエ変換部（周波数領域変換部）１６ａ、１６ｂ、測定周波数成分導出部１８、共役複素数変換部１９、スペクトル導出部２０、逆フーリエ変換部（時間領域変換部）２２を備える。
【００２０】
バンドパスフィルタ１２は、アナログ信号を受け、周波数がFs以上２Fs以下の成分を取り出して、Ａ／Ｄ変換器１４ａ、１４ｂに与える。なお、周波数Fsは、Ａ／Ｄ変換器１４ａ、１４ｂのサンプリング周波数である。
【００２１】
図２は、第一の実施形態にかかるアナログ信号の周波数スペクトルである。アナログ信号は、周波数がFs以上２Fs以下の成分も、それ以外の成分（バンドパスフィルタ１２を通過しない）も有する。ここで、周波数がFs以上２Fs以下の帯域をBand、周波数がFs以上３Fs／２以下の帯域をBand1、周波数が３Fs／２以上２Fs以下の帯域をBand2という。
【００２２】
複数（２個）のＡ／Ｄ変換器１４ａ、１４ｂは、同じ周波数Ｆｓで位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号（すなわち、バンドパスフィルタ１２を通過したアナログ信号）をサンプリングしてデジタル信号に変換する。
【００２３】
図３は、第一の実施形態にかかるＡ／Ｄ変換器１４ａが使用するサンプリングタイミング（図３（ａ））、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミング（図３（ｂ））である。なお、図３においては、縦の実線がサンプリングタイミングを示す。
【００２４】
Ａ／Ｄ変換器１４ａが使用するサンプリングタイミングは、周期がＴ（＝１／Fs）であり、時間が０、Ｔ、２Ｔ、…において、アナログ信号をサンプリングしてデジタル信号に変換する。
【００２５】
Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミングは、周期がＴであり、時間がＴ／２、３Ｔ／２、…において、アナログ信号をサンプリングしてデジタル信号に変換する。すなわち、２π[rad]をＡ／Ｄ変換器の個数（２個）で割った値（π）だけ、サンプリングタイミングの位相が、Ａ／Ｄ変換器１４ａの場合と異なっている。
【００２６】
ただし、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミングには誤差Δがあり、周期がＴであるものの、実際には、時間がＴ／２＋Δ、３Ｔ／２＋Δ、…において、アナログ信号をサンプリングしてデジタル信号に変換する。
【００２７】
誤差Δが無ければ、Ａ／Ｄ変換器１４ａ、１４ｂの出力結果を並べ替えることで、アナログ信号をサンプリング周波数２Fsでサンプリングした信号を得ることができる。
【００２８】
フーリエ変換部（周波数領域変換部）１６ａは、Ａ／Ｄ変換器１４ａの出力を、測定帯域Band1の内の周波数領域の信号に変換する。フーリエ変換部（周波数領域変換部）１６ｂは、Ａ／Ｄ変換器１４ｂの出力を、測定帯域Band1の内の周波数領域の信号に変換する。ただし、測定帯域Band1は、サンプリングタイミングの周波数Fsの２分の１以下の幅の帯域である。
【００２９】
なお、第一の実施形態では、測定帯域Band1の幅は、サンプリングタイミングの周波数Fsの２分の１であるFs／２となっている。
【００３０】
フーリエ変換部１６ａの出力X1と、フーリエ変換部１６ｂの出力X2とを以下に説明する。
【００３１】
図４は、第一の実施形態にかかる出力X1およびX2を説明するために、測定帯域Band1内の信号の周波数f1の成分および帯域Band2内の信号の周波数f2の成分を図示した図である。ただし、(f1+f2)/2 = 3Fs/2である。
【００３２】
ある時間t1における測定帯域Band1内の成分S1（図４においては、周波数f1の成分となる）は、式（１）のように表される。
【００３３】
【数１】
2010283684.tif 000003
ある時間t1における帯域Band2内の成分S2（図４においては、周波数f2の成分となる）は、式（２）のように表される。
【００３４】
【数２】
2010283684.tif 000004
ただし、A1およびA2は、成分S1およびS2の振幅である。
【００３５】
ここで、サンプリングタイミングの周波数Fsなので、折り返し（エイリアシング）が生じないでＡ／Ｄ変換器１４ａ、１４ｂにより測定できる帯域の幅はFs／２である。ここで、折り返しが生じないで、Ａ／Ｄ変換器１４ａが測定できる帯域が測定帯域Band1（Fs〜３Fs／２）である。よって、周波数f2の成分S2が折り返して、周波数f1の成分S1に重なって、フーリエ変換部１６ａの出力となる。
【００３６】
図５は、第一の実施形態にかかる周波数f2の成分S2が折り返して、周波数f1の成分S1に重なった状態を示す図である。
【００３７】
周波数f2の成分S2が折り返した成分をレプリカR12という。レプリカR12は、式（３）のように表される。なお、(f1+f2)/2 = 3Fs/2であるため、レプリカR12は周波数f1の成分S1に重なる。
【００３８】
【数３】
2010283684.tif 000005
ここで、測定帯域Band1内の周波数f1の成分S1を、測定周波数成分という。レプリカR12を重畳成分という。レプリカR12（重畳成分）は、周波数領域の信号において測定周波数成分S1と同一の周波数f1を有する。周波数f2の成分S2を、重畳要因成分という。周波数f2の成分S2（重畳要因成分）は、レプリカR12（重畳成分）の基となっている。すなわち、レプリカR12は、成分S2に基づいて生じている。レプリカR12（重畳成分）は、成分S2（重畳要因成分）のレプリカである。重畳要因成分S2の周波数f2は、測定帯域Band1の外である。
【００３９】
なお、式（２）および式（３）から、レプリカR12（重畳成分）の共役複素数は、成分S2（重畳要因成分）であることがわかる。
【００４０】
また、Ａ／Ｄ変換器１４ａ、１４ｂに用いられるサンプリングタイミングの時間差δtは、式（４）のように表される。なお、δtは、Ａ／Ｄ変換器１４ａ、１４ｂに同じ正弦波（掃引してもよい）を与え、Ａ／Ｄ変換器１４ａ、１４ｂの出力に基づき求めることができる。
【００４１】
【数４】
2010283684.tif 000006
すると、フーリエ変換部１６ａの出力X1は、式（５）のように表される。
【００４２】
【数５】
2010283684.tif 000007
また、フーリエ変換部１６ｂの出力X2は、式（５）におけるS1およびR12のt1をt1+δtに置き換えれば得ることができ、式（６）のように表される。
【００４３】
【数６】
2010283684.tif 000008
測定周波数成分導出部１８および共役複素数変換部１９は、共に、成分導出部を構成する。
【００４４】
測定周波数成分導出部１８は、フーリエ変換部１６ａ、１６ｂの出力する周波数領域の信号X1、X2と、周波数領域の信号の周波数f1と、重畳要因成分の周波数f2と、サンプリングタイミングの時間差δt（サンプリングタイミングの位相差に対応する）とに基づき、測定周波数成分S1と、レプリカR12（重畳成分）とを導出する。
【００４５】
共役複素数変換部１９は、測定周波数成分導出部１８からレプリカR12を受け、レプリカR12の共役複素数として成分S2（重畳要因成分）を導出する。
【００４６】
よって、測定周波数成分導出部１８および共役複素数変換部１９は、測定周波数成分S1と、重畳要因成分S2とを導出する。
【００４７】
なお、式（５）および式（６）から、測定周波数成分S1と、レプリカR12とが、それぞれ式（７）および式（８）のように導出される。
【００４８】
【数７】
2010283684.tif 000009
【００４９】
【数８】
2010283684.tif 000010
測定周波数成分導出部１８は、式（７）および式（８）に、X1、X2、f1、f2、δtを代入して、測定周波数成分S1と、レプリカR12とを導出するようにしてもよい。
【００５０】
また、測定周波数成分導出部１８は、式（７）から測定周波数成分S1を導出し、その後、式（９）に、X1およびS1を代入して、レプリカR12を導出するようにしてもよい。
【００５１】
【数９】
2010283684.tif 000011
また、測定周波数成分導出部１８は、式（８）からレプリカR12を導出し、その後、式（１０）に、X1およびR12を代入して、測定周波数成分S1を導出するようにしてもよい。
【００５２】
【数１０】
2010283684.tif 000012
スペクトル導出部２０は、測定周波数成分S1と、重畳要因成分S2とを有するスペクトルSを導出する。なお、スペクトル導出部２０は、測定周波数成分S1を測定周波数成分導出部１８から受け、重畳要因成分S2を共役複素数変換部１９から受ける。
【００５３】
図６は、第一の実施形態にかかる測定周波数成分S1（図６（ａ））、重畳要因成分S2（図６（ｂ））、スペクトル導出部２０の導出するスペクトルS（図６（ｃ））を示す図である。
【００５４】
測定周波数成分S1および重畳要因成分S2は、測定周波数成分導出部１８および共役複素数変換部１９から出力される。ただし、図６においては、周波数f1、f2成分以外の成分も図示している。
【００５５】
スペクトル導出部２０の導出するスペクトルSは、測定周波数成分S1と、重畳要因成分S2とを有する。すなわち、スペクトルSは、測定帯域Band1内では測定周波数成分S1であり、帯域Band2内では重畳要因成分S2であるようなスペクトルである。
【００５６】
逆フーリエ変換部（時間領域変換部）２２は、スペクトル導出部２０の出力を逆フーリエ変換して、時間領域に変換する。これにより、アナログ信号をサンプリング周波数２Fsでサンプリングしたものが得られることになる。
【００５７】
次に、第一の実施形態の動作を説明する。
【００５８】
まず、アナログ信号（図２参照）がバンドパスフィルタ１２に与えられ、周波数がFs以上２Fs以下の成分を取り出され、Ａ／Ｄ変換器１４ａ、１４ｂに与えられる。Ａ／Ｄ変換器１４ａ、１４ｂのサンプリングタイミング（図３参照）は、ほぼ位相がπずれている。しかし、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミングには誤差Δがあり、実際には、Ａ／Ｄ変換器１４ａが使用するサンプリングタイミングとは、時間がＴ／２＋Δだけずれている（時間差δt）。
【００５９】
フーリエ変換部１６ａ、１６ｂは、Ａ／Ｄ変換器１４ａ、１４ｂの出力を、測定帯域Band1の内の周波数領域の信号に変換する。すると、フーリエ変換部１６ａの出力は、図５を参照して、式（５）のようになる。すなわち、測定周波数成分S1と、重畳要因成分S2のレプリカR12（重畳成分）とが加算されたものとなる。なお、フーリエ変換部１６ｂの出力は、時間差δtを考慮して、式（６）のようになる。
【００６０】
測定周波数成分導出部１８は、式（５）および式（６）から、X1、X2、f1、f2、δtに基づき、測定周波数成分S1とレプリカR12とを導出する（式（７）〜式（１０）参照）。
【００６１】
共役複素数変換部１９は、測定周波数成分導出部１８からレプリカR12を受け、レプリカR12の共役複素数として成分S2（重畳要因成分）を導出する。
【００６２】
スペクトル導出部２０は、測定周波数成分S1と、重畳要因成分S2とを有するスペクトルSを導出する（図６参照）。逆フーリエ変換部（時間領域変換部）２２は、スペクトル導出部２０の出力を逆フーリエ変換して、時間領域に変換する。
【００６３】
第一の実施形態によれば、Ａ／Ｄ変換器１４ａ、１４ｂに用いられるサンプリングタイミングの時間差δtが正確にＴ／２になっていなくても（位相差が正確にはπになっていない）、正確なスペクトルSを導出することができる。
【００６４】
第二の実施形態
第二の実施形態にかかるアナログ信号処理装置１は、Ａ／Ｄ変換器を３個（Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃ）備えた点が、第一の実施形態にかかるアナログ信号処理装置１と異なる。
【００６５】
図７は、本発明の第二の実施形態にかかるアナログ信号処理装置１の構成を示すブロック図である。第二の実施形態にかかるアナログ信号処理装置１は、バンドパスフィルタ１２、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃ、フーリエ変換部（周波数領域変換部）１６ａ、１６ｂ、１６ｃ、測定周波数成分導出部１８、共役複素数変換部１９、スペクトル導出部２０、逆フーリエ変換部（時間領域変換部）２２を備える。
【００６６】
バンドパスフィルタ１２は、アナログ信号を受け、周波数が３Fs／２以上３Fs以下の成分を取り出して、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃに与える。なお、周波数Fsは、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃのサンプリング周波数である。
【００６７】
図８は、第二の実施形態にかかるアナログ信号の周波数スペクトルである。アナログ信号は、周波数が３Fs／２以上３Fs以下の成分も、それ以外の成分（バンドパスフィルタ１２を通過しない）も有する。ここで、周波数が３Fs／２以上３Fs以下の帯域をBand、周波数が３Fs／２以上２Fs以下の帯域をBand1、周波数が２Fs以上５Fs／２以下の帯域をBand2、周波数が５Fs／２以上３Fs以下の帯域をBand3という。
【００６８】
複数（３個）のＡ／Ｄ変換器１４ａ、１４ｂ、１４ｃは、同じ周波数Ｆｓで位相がそれぞれ異なるサンプリングタイミングで、同じアナログ信号（すなわち、バンドパスフィルタ１２を通過したアナログ信号）をサンプリングしてデジタル信号に変換する。
【００６９】
図９は、第二の実施形態にかかるＡ／Ｄ変換器１４ａが使用するサンプリングタイミング（図９（ａ））、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミング（図９（ｂ））、Ａ／Ｄ変換器１４ｃが使用するサンプリングタイミング（図９（ｃ））である。なお、図９においては、縦の実線がサンプリングタイミングを示す。
【００７０】
Ａ／Ｄ変換器１４ａが使用するサンプリングタイミングは、周期がＴ（＝１／Fs）であり、時間が０、Ｔ、２Ｔ、…において、アナログ信号をサンプリングしてデジタル信号に変換する。
【００７１】
Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミングは、周期がＴであり、時間がＴ／３、４Ｔ／３、…において、アナログ信号をサンプリングしてデジタル信号に変換する。すなわち、２π[rad]をＡ／Ｄ変換器の個数（３個）で割った値（２π／３）だけ、サンプリングタイミングの位相が、Ａ／Ｄ変換器１４ａの場合と異なっている。
【００７２】
ただし、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミングには誤差Δ１があり、周期がＴであるものの、実際には、時間がＴ／３＋Δ１、４Ｔ／３＋Δ１、…において、アナログ信号をサンプリングしてデジタル信号に変換する。
【００７３】
Ａ／Ｄ変換器１４ｃが使用するサンプリングタイミングは、周期がＴであり、時間が２Ｔ／３、５Ｔ／３、…において、アナログ信号をサンプリングしてデジタル信号に変換する。すなわち、２π[rad]をＡ／Ｄ変換器の個数（３個）で割った値（２π／３）だけ、サンプリングタイミングの位相が、Ａ／Ｄ変換器１４ｂの場合と異なっている。
【００７４】
ただし、Ａ／Ｄ変換器１４ｃが使用するサンプリングタイミングには誤差Δ２があり、周期がＴであるものの、実際には、時間が２Ｔ／３＋Δ２、５Ｔ／３＋Δ２、…において、アナログ信号をサンプリングしてデジタル信号に変換する。
【００７５】
誤差Δ１、Δ２が無ければ、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃの出力結果を並べ替えることで、アナログ信号をサンプリング周波数３Fsでサンプリングした信号を得ることができる。
【００７６】
フーリエ変換部（周波数領域変換部）１６ａは、Ａ／Ｄ変換器１４ａの出力を、測定帯域Band1の内の周波数領域の信号に変換する。フーリエ変換部（周波数領域変換部）１６ｂは、Ａ／Ｄ変換器１４ｂの出力を、測定帯域Band1の内の周波数領域の信号に変換する。フーリエ変換部（周波数領域変換部）１６ｃは、Ａ／Ｄ変換器１４ｃの出力を、測定帯域Band1の内の周波数領域の信号に変換する。ただし、測定帯域Band1は、サンプリングタイミングの周波数Fsの２分の１以下の幅の帯域である。
【００７７】
なお、第二の実施形態では、測定帯域Band1の幅は、サンプリングタイミングの周波数Fsの２分の１であるFs／２となっている。
【００７８】
フーリエ変換部１６ａの出力X1と、フーリエ変換部１６ｂの出力X2と、フーリエ変換部１６ｃの出力X3とを以下に説明する。
【００７９】
図１０は、第二の実施形態にかかる出力X1、X2およびX3を説明するために、測定帯域Band1内の信号の周波数f1の成分、帯域Band2内の信号の周波数f2の成分および、帯域Band3内の信号の周波数f3の成分を図示した図である。ただし、(f1+f2)/2 = 3Fs/2であり、(f2+f3)/2
= 5Fs/2である。
【００８０】
ある時間t1における測定帯域Band1内の成分S1（図１０においては、周波数f1の成分となる）は、式（１１）のように表される。
【００８１】
【数１１】
2010283684.tif 000013
ある時間t1における帯域Band2内の成分S2（図１０においては、周波数f2の成分となる）は、式（１２）のように表される。
【００８２】
【数１２】
2010283684.tif 000014
ある時間t1における帯域Band3内の成分S3（図１０においては、周波数f3の成分となる）は、式（１３）のように表される。
【００８３】
【数１３】
2010283684.tif 000015
ただし、A1、A2およびA3は、成分S1、S2およびS3の振幅である。
【００８４】
ここで、サンプリングタイミングの周波数Fsなので、折り返し（エイリアシング）が生じないでＡ／Ｄ変換器１４ａ、１４ｂ、１４ｃにより測定できる帯域の幅はFs／２である。ここで、折り返しが生じないで、Ａ／Ｄ変換器１４ａが測定できる帯域が測定帯域Band1（３Fs／２〜２Fs）である。よって、周波数f2の成分S2が折り返して、周波数f1の成分S1に重なり、さらに、周波数f3の成分S3が折り返して、周波数f1の成分S1に重なり、フーリエ変換部１６ａの出力となる。
【００８５】
図１１は、第二の実施形態にかかる周波数f2の成分S2が折り返し、さらに周波数f3の成分S3が折り返して、周波数f1の成分S1に重なった状態を示す図である。
【００８６】
周波数f2の成分S2が折り返した成分をレプリカR12という。周波数f3の成分S3が折り返した成分をレプリカR13という。
【００８７】
レプリカR12は、式（１４）のように表される。なお、(f1+f2)/2 = 3Fs/2であるため、レプリカR12は周波数f1の成分S1に重なる。
【００８８】
【数１４】
2010283684.tif 000016
レプリカR13は、式（１５）のように表される。なお、(f1+f2)/2 = 3Fs/2であり、(f2+f3)/2
= 5Fs/2であるため、レプリカR13は周波数f1の成分S1に重なる。
【００８９】
【数１５】
2010283684.tif 000017
ここで、測定帯域Band1内の周波数f1の成分S1を、測定周波数成分という。レプリカR12、R13を重畳成分という。レプリカR12、R13（重畳成分）は、周波数領域の信号において測定周波数成分S1と同一の周波数f1を有する。
【００９０】
周波数f2の成分S2を、重畳要因成分という。周波数f2の成分S2（重畳要因成分）は、レプリカR12（重畳成分）の基となっている。すなわち、レプリカR12は、成分S2に基づいて生じている。レプリカR12（重畳成分）は、成分S2（重畳要因成分）のレプリカである。重畳要因成分S2の周波数f2は、測定帯域Band1の外である。
【００９１】
周波数f3の成分S3を、重畳要因成分という。周波数f3の成分S3（重畳要因成分）は、レプリカR13（重畳成分）の基となっている。すなわち、レプリカR13は、成分S3に基づいて生じている。レプリカR13（重畳成分）は、成分S3（重畳要因成分）のレプリカである。重畳要因成分S3の周波数f3は、測定帯域Band1の外である。
【００９２】
なお、式（１２）および式（１４）から、レプリカR12（重畳成分）の共役複素数は、成分S2（重畳要因成分）であることがわかる。また、式（１３）および式（１５）から、レプリカR13（重畳成分）そのものが、成分S3（重畳要因成分）であることがわかる。
【００９３】
また、Ａ／Ｄ変換器１４ａに用いられるサンプリングタイミングと、Ａ／Ｄ変換器１４ｂに用いられるサンプリングタイミングとの時間差δt1は、式（１６）のように表される。Ａ／Ｄ変換器１４ａに用いられるサンプリングタイミングと、Ａ／Ｄ変換器１４ｃに用いられるサンプリングタイミングとの時間差δt2は、式（１７）のように表される。
【００９４】
【数１６】
2010283684.tif 000018
【００９５】
【数１７】
2010283684.tif 000019
すると、フーリエ変換部１６ａの出力X1は、式（１８）のように表される。
【００９６】
【数１８】
2010283684.tif 000020
また、フーリエ変換部１６ｂの出力X2は、式（１８）におけるS1、R12およびR13のt1をt1+δt1に置き換えれば得ることができ、式（１９）のように表される。
【００９７】
【数１９】
2010283684.tif 000021
さらに、フーリエ変換部１６ｃの出力X3は、式（１８）におけるS1、R12およびR13のt1をt1+δt2に置き換えれば得ることができ、式（２０）のように表される。
【００９８】
【数２０】
2010283684.tif 000022
測定周波数成分導出部１８および共役複素数変換部１９は、共に、成分導出部を構成する。
【００９９】
測定周波数成分導出部１８は、フーリエ変換部１６ａ、１６ｂ、１６ｃの出力する周波数領域の信号X1、X2、X3と、周波数領域の信号の周波数f1と、重畳要因成分の周波数f2、f3と、サンプリングタイミングの時間差δt1、δt2（サンプリングタイミングの位相差に対応する）とに基づき、測定周波数成分S1と、レプリカR12、R13（重畳成分）とを導出する。
【０１００】
共役複素数変換部１９は、測定周波数成分導出部１８からレプリカR12を受け、レプリカR12の共役複素数として成分S2（重畳要因成分）を導出する。
【０１０１】
なお、レプリカR13（重畳成分）そのものが成分S3（重畳要因成分）である。
【０１０２】
よって、測定周波数成分導出部１８および共役複素数変換部１９は、測定周波数成分S1と、重畳要因成分S2と、重畳要因成分S3とを導出する。
【０１０３】
式（１８）、（１９）、（２０）において、X1、X2、X3、f1、f2、f3、δt1、δt2がわかるので、未知の変数S1、R12、R13に対して、式が３個あることになる。よって、変数S1、R12、R13を導出することができる。変数S1、R12、R13を導出する方法は、周知の三元連立一次方程式の解法によればよい。
【０１０４】
スペクトル導出部２０は、測定周波数成分S1と、重畳要因成分S2とを有するスペクトルSを導出する。なお、スペクトル導出部２０は、測定周波数成分S1を測定周波数成分導出部１８から受け、重畳要因成分S2を共役複素数変換部１９から受ける。
【０１０５】
図１２は、第二の実施形態にかかる測定周波数成分S1（図１２（ａ））、重畳要因成分S2（図１２（ｂ））、重畳要因成分S3（図１２（ｃ））、スペクトル導出部２０の導出するスペクトルS（図１２（ｄ））を示す図である。
【０１０６】
測定周波数成分S1および重畳要因成分S2、S3は、測定周波数成分導出部１８および共役複素数変換部１９から出力される。ただし、図１２においては、周波数f1、f2、f3成分以外の成分も図示している。
【０１０７】
スペクトル導出部２０の導出するスペクトルSは、測定周波数成分S1と、重畳要因成分S2、S3とを有する。すなわち、スペクトルSは、測定帯域Band1内では測定周波数成分S1であり、帯域Band2内では重畳要因成分S2であり、帯域Band3内では重畳要因成分S3であるようなスペクトルである。
【０１０８】
逆フーリエ変換部（時間領域変換部）２２は、スペクトル導出部２０の出力を逆フーリエ変換して、時間領域に変換する。これにより、アナログ信号をサンプリング周波数３Fsでサンプリングしたものが得られることになる。
【０１０９】
次に、第二の実施形態の動作を説明する。
【０１１０】
まず、アナログ信号（図８参照）がバンドパスフィルタ１２に与えられ、周波数がFs以上２Fs以下の成分を取り出され、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃに与えられる。Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃのサンプリングタイミング（図９参照）は、ほぼ位相が２π／３ずれている。しかし、Ａ／Ｄ変換器１４ｂが使用するサンプリングタイミングには誤差Δ１があり、実際には、Ａ／Ｄ変換器１４ａが使用するサンプリングタイミングとは、時間がＴ／３＋Δ１だけずれている（時間差δt1）。また、Ａ／Ｄ変換器１４ｃが使用するサンプリングタイミングには誤差Δ２があり、実際には、Ａ／Ｄ変換器１４ａが使用するサンプリングタイミングとは、時間が２Ｔ／３＋Δ２だけずれている（時間差δt2）。
【０１１１】
フーリエ変換部１６ａ、１６ｂ、１６ｃは、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃの出力を、測定帯域Band1の内の周波数領域の信号に変換する。すると、フーリエ変換部１６ａの出力は、図１１を参照して、式（１８）のようになる。すなわち、測定周波数成分S1と、重畳要因成分S2のレプリカR12（重畳成分）と、重畳要因成分S3のレプリカR13（重畳成分）とが加算されたものとなる。なお、フーリエ変換部１６ｂの出力は、時間差δt1を考慮して、式（１９）のようになる。また、フーリエ変換部１６ｃの出力は、時間差δt2を考慮して、式（２０）のようになる。
【０１１２】
測定周波数成分導出部１８は、X1、X2、X3、f1、f2、f3、δt1、δt2に基づき、測定周波数成分S1とレプリカR12、R13(=S3)とを導出する。
【０１１３】
共役複素数変換部１９は、測定周波数成分導出部１８からレプリカR12を受け、レプリカR12の共役複素数として成分S2（重畳要因成分）を導出する。
【０１１４】
スペクトル導出部２０は、測定周波数成分S1と、重畳要因成分S2、S3とを有するスペクトルSを導出する（図１２参照）。逆フーリエ変換部（時間領域変換部）２２は、スペクトル導出部２０の出力を逆フーリエ変換して、時間領域に変換する。
【０１１５】
第二の実施形態によれば、Ａ／Ｄ変換器１４ａ、１４ｂ、１４ｃに用いられるサンプリングタイミングの時間差δt1、δt2が正確にＴ／３、２Ｔ／３になっていなくても（位相差が正確には２π／３になっていない）、正確なスペクトルSを導出することができる。
【０１１６】
なお、Ａ／Ｄ変換器が二個の場合（第一の実施形態）と、三個の場合（第二の実施形態）を説明したが、四個以上であっても本願発明を実施できる。
【０１１７】
また、上記の実施形態は、以下のようにして実現できる。ＣＰＵ、ハードディスク、メディア（フロッピー（登録商標）ディスク、ＣＤ−ＲＯＭなど）読み取り装置を備えたコンピュータに、上記の各部分、例えばフーリエ変換部（周波数領域変換部）１６ａ、１６ｂ、１６ｃ、測定周波数成分導出部１８、共役複素数変換部１９、スペクトル導出部２０、逆フーリエ変換部（時間領域変換部）２２を実現するプログラムを記録したメディアを読み取らせて、ハードディスクにインストールする。このような方法でも、上記の機能を実現できる。
【符号の説明】
【０１１８】
１  アナログ信号処理装置
１２  バンドパスフィルタ
１４ａ、１４ｂ、１４ｃ  Ａ／Ｄ変換器
１６ａ、１６ｂ、１６ｃ  フーリエ変換部（周波数領域変換部）
１８  測定周波数成分導出部
１９  共役複素数変換部
２０  スペクトル導出部
２２  逆フーリエ変換部（時間領域変換部）
Band1  測定帯域
S1  測定周波数成分
R12、R13  レプリカ（重畳成分）
S2、S3  重畳要因成分
【図１】
2010283684.tif 000023
【図２】
2010283684.tif 000024
【図３】
2010283684.tif 000025
【図４】
2010283684.tif 000026
【図５】
2010283684.tif 000027
【図６】
2010283684.tif 000028
【図７】
2010283684.tif 000029
【図８】
2010283684.tif 000030
【図９】
2010283684.tif 000031
【図１０】
2010283684.tif 000032
【図１１】
2010283684.tif 000033
【図１２】
2010283684.tif 000034