(57)【要約】【課題】リアルタイム且つ高精度に二相正弦波信号の振幅、位相、オフセット誤差を検出し、内挿誤差の補正を行なう。【解決手段】所定のサンプリング周期でサンプリングして得られた位相差９０°の二相正弦波信号φA(n)、φB(n)から算出される位置に含まれる内挿誤差を補正する際、前記二相正弦波信号φA(n)、φB(n)に対して漸化的最小二乗法を適用して求めた振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)に基づいて、所定の更新周期毎に振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを更新する補正を行ない、更新された前記各補正誤差を用いて得られた二相正弦波信号φ’A(n)、φ’B(n)に、前記サンプリング周期毎に再度漸化的最小二乗法を適用して振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を推定し、前記二相正弦波信号φ’A(n)、φB’(n)の補正を行なう。【選択図】図２
2005172669.tif 000002
【特許請求の範囲】
【請求項１】
所定のサンプリング周期でサンプリングして得られた位相差９０°の二相正弦波信号φA(n)、φB(n)から算出される位置に含まれる内挿誤差を補正する内挿誤差補正方法において、
前記二相正弦波信号φA(n)、φB(n)に対して漸化的最小二乗法を適用して求めた振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)に基づいて、所定の更新周期Ｎ毎に振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを更新する補正を行ない、
更新された前記各補正誤差を用いて得られた二相正弦波信号φ’A(n)、φ’B(n)に、前記サンプリング周期毎に再度漸化的最小二乗法を適用して、振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を推定し、前記二相正弦波信号φ’A(n)、φ’B(n)の補正を行なうことを特徴とする内挿誤差補正方法。
【請求項２】
前記サンプリング周期に対して、前記二相正弦波信号φA(n)、φB(n)の周期が非常に長い場合、漸化的最小二乗法による前期振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)の推定と、更新周期Ｎ毎に更新していた前記振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bの補正を停止し、停止以前に求めた前記振幅、位相、オフセット誤差を用いて二相正弦波信号の補正を行なうことを特徴とする請求項１に記載の内挿誤差補正方法。
【請求項３】
所定のサンプリング周期でサンプリングして得られた位相差９０°の二相正弦波信号φA(n)、φB(n)から算出される位置に含まれる内挿誤差を補正する内挿誤差補正装置において、
前記二相正弦波信号φA(n)、φB(n)に対して漸化的最小二乗法を適用して求めた振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)に基づいて、所定の更新周期Ｎ毎に振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを更新する補正を行なう手段と、
更新された前記各補正誤差を用いて得られた二相正弦波信号φ’A(n)、φ’B(n)に、前記サンプリング周期毎に再度漸化的最小二乗法を適用して、振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を推定し、前記二相正弦波信号φ’A(n)、φ’B(n)の補正を行なう手段と、を備えたことを特徴とする内挿誤差補正装置。
【請求項４】
前記サンプリング周期に対して、前記二相正弦波信号φA(n)、φB(n)の周期が非常に長い場合、漸化的最小二乗法による前期振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)の推定と、更新周期Ｎ毎に更新していた前記振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bの補正を停止し、停止以前に求めた前記振幅、位相、オフセット誤差を用いて二相正弦波信号の補正を行なうことを特徴とする請求項３に記載の内挿誤差補正装置。
【請求項５】
請求項１又は２に記載の内挿誤差補正方法を実施するためのコンピュータプログラム。
【請求項６】
請求項３又は４に記載の内挿誤差補正装置を実現するためのコンピュータプログラム。
【請求項７】
請求項５又は６に記載のコンピュータプログラムが記録された、コンピュータ読取り可能な記録媒体。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、内挿誤差補正方法及び装置に係り、特に位置検出に使用される位相差９０°の二相正弦波信号を内挿して、高分解能の位置を検出する際に適用して好適な内挿誤差補正方法及び装置に関する。
【背景技術】
【０００２】
位相差９０°の二相正弦波信号を内挿して、高分解能の位置を検出する位置検出装置において、位置を検出する際に使用される二相正弦波信号は、理想的な正弦波であることが仮定されている。しかし、実際には、検出部を構成する素子のノイズや温度ドリフト等により、二相正弦波信号には振幅、位相、オフセット誤差が含まれている。その結果、検出部からサンプリング入力される二相正弦波信号から求めた位置には、理想的な正弦波の場合に対して内挿誤差が含まれてしまう。
【０００３】
この内挿誤差の補正には、（１）別の手段を用いて前もって内挿誤差を測定しておくテーブル補正方式（例えば、特許文献１参照）、（２）二相正弦波信号の最大値、最小値等から、振幅、位相、オフセット誤差を演算し内挿誤差を補正する最大振幅値方式（例えば、特許文献２参照）、（３）二相正弦波信号の１周期分のデータを用いて、振幅、位相、オフセット誤差を算出し内挿誤差を補正する最小二乗円方式（例えば、特許文献３参照）の３つが、主に行なわれている。
【０００４】
【特許文献１】特開平５−２３１８７９号公報
【特許文献２】特開平６−１６７３５４号公報
【特許文献３】特開２００３−２５４７８４号公報
【発明の開示】
【発明が解決しようとする課題】
【０００５】
しかしながら、前記（１）テーブル補正方式は、内挿誤差の時間的変化に対応しておらず、又、内挿誤差を測定するために、より高精度な基準が必要であるという欠点がある。
【０００６】
前記（２）最大振幅値方式は、計算時間が短くリアルタイム補正が可能であるが、補正を行なう際に特定のデータしか用いていないため、ノイズの影響を受け易いという欠点がある。
【０００７】
前記（３）最小二乗円方式は、ノイズの影響を受け難く、内挿誤差もかなり小さくできるが、計算時間が長くリアルタイム補正が難しいという欠点がある。
【０００８】
本発明は、前記従来の問題点を解決するべくなされたもので、リアルタイム且つ高精度に二相正弦波信号の振幅、位相、オフセット誤差を算出し、内挿誤差を補正することができる内挿誤差補正方法及び装置を提供することを課題とする。
【課題を解決するための手段】
【０００９】
本発明は、所定のサンプリング周期でサンプリングして得られた位相差９０°の二相正弦波信号φA(n)、φB(n)から算出される位置に含まれる内挿誤差を補正する内挿誤差補正方法において、前記二相正弦波信号φA(n)、φB(n)に対して漸化的最小二乗法を適用して求めた振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)に基づいて、所定の更新周期Ｎ毎に振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを更新する補正を行ない、更新された前記各補正誤差を用いて得られた二相正弦波信号φ’A(n)、φ’B(n)に、前記サンプリング周期毎に再度漸化的最小二乗法を適用して、振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を推定し、前記二相正弦波信号φ’A(n)、φ’B(n)の補正を行なうことにより、前記課題を解決したものである。
【００１０】
本発明は、又、所定のサンプリング周期でサンプリングして得られた位相差９０°の二相正弦波信号φA(n)、φB(n)から算出される位置に含まれる内挿誤差を補正する内挿誤差補正装置において、前記二相正弦波信号φA(n)、φB(n)に対して漸化的最小二乗法を適用して求めた振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)に基づいて、所定の更新周期Ｎ毎に振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを更新する補正を行なう手段と、更新された前記各補正誤差を用いて得られた二相正弦波信号φ’A(n)、φ’B(n)に、前記サンプリング周期毎に再度漸化的最小二乗法を適用して、振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を推定し、前記二相正弦波信号φ’A(n)、φ’B(n)の補正を行なう手段と、を備えたことにより、同様に前記課題を解決したものである。
【００１１】
又、前記サンプリング周期に対して、前記二相正弦波信号φA(n)、φB(n)の周期が非常に長い場合、漸化的最小二乗法による前期振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)の推定と、更新周期Ｎ毎に更新していた前記振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bの補正を停止し、停止以前に求めた前記振幅、位相、オフセット誤差を用いて二相正弦波信号の補正を行なうようにしたものである。
【００１２】
本発明は又、前記の内挿誤差補正方法を実施するためのコンピュータプログラムを提供するものである。
【００１３】
又、前記の内挿誤差補正装置を実現するためのコンピュータプログラムを提供するものである。
【００１４】
又、前記のコンピュータプログラムが記録された、コンピュータ読取り可能な記録媒体を提供するものである。
【発明の効果】
【００１５】
本発明によれば、リアルタイム且つ高精度に二相正弦波信号の振幅、位相、オフセット誤差を算出し、内挿誤差を補正することができる。従って、高分解能の位置をリアルタイム且つ高精度に検出することができる。
【発明を実施するための最良の形態】
【００１６】
以下、図面を参照して、本発明の実施の形態について詳細に説明する。
【００１７】
本発明に係る一実施形態の内挿誤差補正方法が適用される位置検出装置の概要を図１のブロック図に、そのアルゴリズムを図２のフローチャートに示す。
【００１８】
本実施形態は、以下に、数式に基づき具体的に説明するように、所定のサンプリング周期毎にサンプリングされた二相正弦波信号に対して、漸化的最小二乗法（Ｒecursive  Ｌeast  Ｓquares  method、以下ＲＬＳ法と略記する）を適用して、振幅、位相、オフセットを補正することにより、内挿誤差を補正する際、該ＲＬＳ法を適用する前に、サンプリングされた前記二相正弦波信号に対して、ＲＬＳ法によって既に得られている補正値に基づいた補正を行なうことにより、より高精度な補正を実現しようとするものである。なお、このＲＬＳ法については、例えば、足立修一著：“ＭＡＴＬＡＢによる制御のためのシステム同定”、東京電機大学出版局（１９９６）に説明されている。
【００１９】
いま、図１の位置検出装置において、検出部１０から所定のサンプリング周期でサンプリングされ、Ａ／Ｄ変換されてコンピュータからなるデータ処理部１２に入力される、位相差９０°の二相正弦波信号φA(n)、φB(n)を、式（１）で表わす。
【００２０】
【数１】
2005172669.tif 000003
【００２１】
ここで、ΔA(n)、ΔB(n)は振幅誤差、ΔαA(n)、ΔαB(n)は位相誤差、ΔＶA(n)、ΔＶB(n)はオフセット誤差、ｎはサンプリング番号である。
【００２２】
前記データ処理部１２において、プログラムにより構築された各演算手段（図示せず）により、図２に示すアルゴリズムに従って式（１）に対してＲＬＳ法が適用されるようになっている。具体的なＲＬＳ法の適用方法を以下に説明する。
【００２３】
最初に、式（１）に対してＲＬＳ法が適用できるように該式の線形化を行なう。
【００２４】
まず、二相正弦波信号φA(n)、φB(n)から、式（２）を用いて角度θ’(n)を求める（ステップ１０１）。但し、最初の補正処理以外は、二相正弦波信号φA(n)、φB(n)に対しては、後に詳述するステップ１００の補正を予め行ない、補正された二相正弦波信号φ’A(n)、φ’B(n)を得る。そして、二相正弦波信号φ’A(n)、φ’B(n)を二相正弦波信号φA(n)、φB(n)と考えて以下の演算を行う。
【００２５】
【数２】
2005172669.tif 000004
【００２６】
そして、ｘ(n)を式（３）のように定義する（ステップ１０２）。すると、式（１）は、下記の如く式（４）、式（５）のように線形化することができる。ここでは、式（１）のθ(n)と式（２）のθ’(n)は等しいものと考えて式（１）の線形化を行なう。なお、以下の式で、Ｔは転置行列を表わす。
【００２７】
ｘ(n)＝［ｘ1(n)  ｘ2(n)  ｘ3(n)］T＝［sinθ’(n)  cosθ’(n)  １］T
…（３）
φA(n)＝（１＋ΔA(n)）cos（θ’(n)＋ΔαA(n)）＋ΔＶA(n)
＝（１＋ΔA(n)）sin（θ’(n)＋ΔαA(n)＋π／２）＋ΔＶA(n)
＝（１＋ΔA(n)）sinθ’(n)cos（ΔαA(n)＋π／２）
＋（１＋ΔA(n)）cosθ’(n)sin（ΔαA(n)＋π／２）＋ΔＶA(n)
＝ｃA1(n)ｘ1(n)＋ｃA2(n)ｘ2(n)＋ｃA3(n)ｘ3(n)
＝ｃTA(n)ｘ(n)                                      …（４）
φB(n)＝（１＋ΔB(n)）sin（θ’(n)＋ΔαB(n)）＋ΔＶB(n)
＝（１＋ΔB(n)）sinθ’(n)cosΔαB(n)
＋（１＋ΔB(n)）cosθ’(n)sinΔαB(n)＋ΔＶB(n)
＝ｃB1(n)ｘ1(n)＋ｃB2(n)ｘ2(n)＋ｃB3(n)ｘ3(n)
＝ｃTB(n)ｘ(n)                                      …（５）
【００２８】
ここで、
【数３】
2005172669.tif 000005
である。
【００２９】
この式（４）、式（５）に対してＲＬＳ法を適用して、ｃA(n)、ｃB(n)を算出する。但し、予め下記の式（１２）の関係を求めておく（ステップ１０３）。
【００３０】
具体的には、ｘ(n)と二相正弦波信号φA(n)、φB(n)を用いて、式（１０）、式（１１）によりｃA (n)、ｃB (n)を算出する（ステップ１０４、１０５）。
【００３１】
ｃA(n)＝ｃA(n-1)＋Ｑ(n)ｘ(n)（φA(n)−ｘT(n)ｃA(n-1)）        …（１０）
ｃB(n)＝ｃB(n-1)＋Ｑ(n)ｘ(n)（φB(n)−ｘT(n)ｃB(n-1)）        …（１１）
ｃA(0)＝［０  １  ０］T
ｃB(0)＝［１  ０  ０］T
【００３２】
ここで、Ｑ(n)は式（１２）から算出される行列である。
【００３３】
【数４】
2005172669.tif 000006
【００３４】
なお、Ｑ(0)の値としては、最初は単位行列等を用いる。Ｑ(n)が十分に収束した場合は、次回以降からはこのＱ(n)の値をＱ(0)として用いる。又、γは忘却要素であり、１以下の正数である。
【００３５】
式（１０）、式（１１）により求めたｃA(n)、ｃB(n)を、式（６）、式（８）、式（７）、式（９）に適用して、振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を算出し（ステップ１０６、１０７）、各誤差を用いて二相正弦波信号φA(n)、φB(n)の補正を行ない、補正された二相正弦波信号φ”A(n)、φ”B(n)を得る（ステップ１０８）。そして、φ”A(n)、φ”B(n)から位置を算出することで内挿誤差を補正する。
【００３６】
通常のＲＬＳ法では、前述した如くθ(n)＝θ’(n)と考え、式（１）から式（４）、式（５）の線形化を行なっている。しかし、実際には二相正弦波信号φA(n)、φB(n)に含まれている振幅、位相、オフセット誤差のため、θ(n)＝θ’(n)とはならない。その結果、式（４）、式（５）の線形化には誤差が含まれてしまい、式（４）、式（５）にＲＬＳ法を適用して算出した振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)の推定結果にも誤差が含まれてしまうことになる。
【００３７】
そこで、本実施形態では、式（２）を適用してθ’(n)を算出する前に、ＲＬＳ法によって既に得られている補正値から求めた振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを用いて、サンプリング入力される二相正弦波信号φA(n)、φB(n)の補正を行なう。
【００３８】
その結果、二相正弦波信号φA(n)、φB(n)に含まれている振幅、位相、オフセット誤差は小さくなり、ＲＬＳ法を適用して求めた振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)の推定結果に含まれている誤差も小さくできる。
【００３９】
即ち、ここでは、前記ステップ１０１において式（２）を適用する前に、ステップ１００で行なう補正に用いる振幅誤差Δ’A、Δ’B、位相誤差Δα’A、Δα’B、オフセット誤差ΔＶ’A、ΔＶ’Bを、それ以前のＲＬＳ法による推定結果から予め求めておき、更新周期Ｎで式（１３）に示す更新を行なっている（ステップ１０９）。
【００４０】
【数５】
2005172669.tif 000007
【００４１】
例えば、更新周期Ｎ＝２００とした場合、振幅誤差Δ’Aは式（１４）のようになる。
【００４２】
【数６】
2005172669.tif 000008
【００４３】
なお、更新周期Ｎ毎に更新を行なっている理由は、ＲＬＳ法による推定結果が十分に収束したときの振幅誤差ΔA(n)、ΔB(n)、位相誤差ΔαA(n)、ΔαB(n)、オフセット誤差ΔＶA(n)、ΔＶB(n)を用いるためである。但し、各補正値の初期値は０とする。又、忘却要素γはγ＝０．９９、更新周期ＮはＮ＝１００〜２００程度が好ましい。
【００４４】
又、本実施形態では、振幅、位相、オフセット誤差の時間的変化にも対応できるように忘却要素を用いている。そのため、サンプリング周期に比べて二相正弦波信号の周期が非常に長い場合、振幅、位相、オフセット誤差の推定結果が発散してしまう場合がある。どのくらいの周期で発散してしまうかに関してはシミュレーション結果を図３に示す。この図から、例えば更新周期Ｎ＝１００のときは、サンプリング周期と比較して二相正弦波信号の周期が約３５０倍以下の場合は振幅、位相、オフセット誤差の推定結果は収束するが、約３５０倍以上の場合は振幅、位相、オフセット誤差の推定結果は発散してしまう。
【００４５】
そこで、そうした場合は式（２）から式（１３）による振幅、位相、オフセット誤差の計算を停止し、停止以前に求めた振幅、位相、オフセット誤差を用いて二相正弦波信号の補正を行なうことにより、前記問題に対処することが可能となる。
【００４６】
上記計算を停止するか否かの判断は、任意のサンプリング時に式（２）から求まる角度θ’(n)と、１サンプリング前の角度θ’(n-1)の差の絶対値を求め、その値が基準値より小さいか大きいかで行う。
【００４７】
例えば更新周期Ｎ＝１００のときは、前記図３から、サンプリング周期と比較して二相正弦波信号の周期が約３５０倍以上のときに振幅、位相、オフセット誤差の推定結果が発散してしまうことがわかる。そこで、この例では、角度θ’(n)と１サンプリング前の角度θ’(n-1)の差の絶対値が（２π／３５０）以下の時は計算を停止する。
【００４８】
即ち、上記計算を停止するか否かを判断する基準値を、（二相正弦波信号の周期）／（サンプリング周期）の値が、更新周期に応じて変化する上記推定誤差の発散と収束の境界線上に一致する境界倍率値に基づいて設定し、具体的には（２π／境界倍率値）で設定する。
【００４９】
次に、本発明に係る実施形態の効果をシミュレーションにより確認した結果を示す。適用したシミュレーション条件としては、二相正弦波信号の振幅誤差２％、位相誤差２°、オフセット誤差２％、ノイズを振幅の±２％、１周期のデータ数は１００個とした。又、内挿誤差の時間的変化として、オフセット誤差が１周期当たり０．２％ずつ増加していく時間的変化を与えた。
【００５０】
このシミュレーション条件のときのリサージュ波形を図４に示す。この図４から本実施形態による補正を行なうことにより、リサージュ波形は理想値である真円に近づいていることが分かる。又、図５に、波長λで正規化した位置Ｐについて、このときの内挿誤差ΔＰを示す。この図５から補正を行なうことにより内挿誤差が小さくなっていることが分かる。
【００５１】
又、図５の内挿誤差の空間周波数解析結果を図６に示す。この図６から本実施形態による補正を行なうことにより、内挿誤差であるλ、λ／２の周期の誤差が減少していることが分かる。
【００５２】
以上詳述した如く、本実施形態においては、まず、二相正弦波信号に対して、ある振幅、位相、オフセット誤差を用いて１回目の補正を行なう。そして、補正された二相正弦波信号に漸化的最小二乗法（ＲＬＳ法）を適用して振幅、位相、オフセット誤差を推定し、再度補正を行なう。更にＲＬＳ法による推定結果を１回目の補正に利用する。このようにして、本実施形態によれば、リアルタイム且つ高精度な内挿誤差の補正を実現することができる。
【００５３】
なお、前記実施形態では特に明示しなかったが、サンプリングして得られる位相差９０°の二相正弦波信号は、対応する二相信号をサンプリングする場合に限らず、多相正弦波信号を所定のサンプリング周期でサンプリングして得られた信号から求めるようにしてもよい。
【００５４】
又、サンプリング周期、更新周期Ｎは任意に設定可能であることは言うまでもない。
【産業上の利用可能性】
【００５５】
以上詳述したとおり、本発明によれば、位置検出装置により位相差９０°の二相正弦波信号を使って位置を検出する際に、リアルタイム且つ高精度に二相正弦波信号の振幅、位相、オフセット誤差を検出し、内挿誤差の補正を行なうことができることから、高分解能の位置検出をリアルタイム且つ高精度に行なうことができる。
【図面の簡単な説明】
【００５６】
【図１】本発明に係る一実施形態に適用される位置検出装置の概要を示すブロック図
【図２】上記実施形態のアルゴリズムを示すフローチャート
【図３】補正誤差のシミュレーション結果の一例を示す線図
【図４】本発明の効果を示す線図
【図５】本発明の効果を示す他の線図
【図６】本発明の効果を示す更に他の線図
【符号の説明】
【００５７】
１０…検出部
１２…データ処理部
【図１】
2005172669.tif 000009
【図２】
2005172669.tif 000010
【図３】
2005172669.tif 000011
【図４】
2005172669.tif 000012
【図５】
2005172669.tif 000013
【図６】
2005172669.tif 000014