(57)【要約】【課題】ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質が含まれるか否かを精度良く測定する。【解決手段】本解析装置は、ある物質によって共振周波数に変化が生ずる測定系について、測定系本来の共振周波数を含む所定帯域内の特定の周波数における伝送特性、反射特性又はインピーダンス特性を表すデータを、ある物質の異なる複数の状態について取得する手段と、取得された複数回分の伝送特性、反射特性又はインピーダンス特性を表すデータから、所定の特徴量を算出する算出手段と、他の特定の物質の濃度若しくは含有量と所定の特徴量との関係又は他の物質が含まれているか否かと所定の特徴量との関係を表す関係データを格納したデータ格納部に格納された関係データに基づき、算出された所定の特徴量に対応する他の特定の物質の濃度若しくは含有量、又は他の物質が含まれているか否かを決定する手段とを有する。【選択図】図１
2020148713.tif 000002
【特許請求の範囲】
【請求項１】
ある物質に含まれる他の特定の物質の濃度若しくは含有量、又は前記ある物質に他の物質が含まれるか否かを解析する解析装置であって、
前記ある物質によって共振周波数に変化が生ずる測定系について、前記測定系本来の共振周波数を含む所定帯域内の特定の周波数における伝送特性、反射特性又はインピーダンス特性を表すデータを、前記ある物質の異なる複数の状態について取得する手段と、
取得された複数回分の前記伝送特性、反射特性又はインピーダンス特性を表すデータから、所定の特徴量を算出する算出手段と、
前記他の特定の物質の濃度若しくは含有量と所定の特徴量との関係又は前記他の物質が含まれているか否かと所定の特徴量との関係を表す関係データを格納したデータ格納部に格納された前記関係データに基づき、算出された前記所定の特徴量に対応する前記他の特定の物質の濃度若しくは含有量、又は前記他の物質が含まれているか否かを決定する手段と、
を有する解析装置。
【請求項２】
前記算出手段が、
複数回分の前記伝送特性、反射特性又はインピーダンス特性を表すデータで特定される振幅変化量及び位相変化量のデータから、振幅変化量及び位相変化量で張られる平面におけるフィッティングカーブを算出し、当該フィッティングカーブの係数を前記所定の特徴量として算出する
請求項１記載の解析装置。
【請求項３】
前記算出手段が、
複数回分の前記伝送特性、反射特性又はインピーダンス特性を表すデータで特定される実数部及び虚数部のデータから、実数部及び虚数部で張られる平面におけるフィッティングカーブを算出し、当該フィッティングカーブの係数を前記所定の特徴量として算出する
請求項１記載の解析装置。
【請求項４】
前記フィッティングカーブが、２次関数のフィッティングカーブである
請求項２記載の解析装置。
【請求項５】
前記フィッティングカーブが、２次関数のフィッティングカーブ又は円のフィッティングカーブである
請求項３記載の解析装置。
【請求項６】
ある物質に含まれる他の特定の物質の濃度若しくは含有量、又は前記ある物質に他の物質が含まれるか否かを解析する解析方法であって、
コンピュータが、
前記ある物質によって共振周波数に変化が生ずる測定系について、前記測定系本来の共振周波数を含む所定帯域内の特定の周波数における伝送特性、反射特性又はインピーダンス特性を表すデータを、前記ある物質の異なる複数の状態について取得するステップと、
取得された複数回分の前記伝送特性、反射特性又はインピーダンス特性を表すデータから、所定の特徴量を算出するステップと、
前記他の特定の物質の濃度若しくは含有量と所定の特徴量との関係又は前記他の物質が含まれているか否かと所定の特徴量との関係を表す関係データを格納したデータ格納部に格納された前記関係データに基づき、算出された前記所定の特徴量に対応する前記他の特定の物質の濃度若しくは含有量、又は前記他の物質が含まれているか否かを決定するステップと、
を実行する解析方法。
【請求項７】
ある物質に含まれる他の特定の物質の濃度若しくは含有量、又は前記ある物質に他の物質が含まれるか否かを解析する解析プログラムであって、
コンピュータに、
前記ある物質によって共振周波数に変化が生ずる測定系について、前記測定系本来の共振周波数を含む所定帯域内の特定の周波数における伝送特性、反射特性又はインピーダンス特性を表すデータを、前記ある物質の異なる複数の状態について取得するステップと、
取得された複数回分の前記伝送特性、反射特性又はインピーダンス特性を表すデータから、所定の特徴量を算出するステップと、
前記他の特定の物質の濃度若しくは含有量と所定の特徴量との関係又は前記他の物質が含まれているか否かと所定の特徴量との関係を表す関係データを格納したデータ格納部に格納された前記関係データに基づき、算出された前記所定の特徴量に対応する前記他の特定の物質の濃度若しくは含有量、又は前記他の物質が含まれているか否かを決定するステップと、
を実行させるための解析プログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質が含まれるか否かを解析する技術に関する。
【背景技術】
【０００２】
例えば、米等の粉状、粒状又は不定形のバラ状の対象物の水分量や材料の特性を、体積充填率の変動に影響されず、オンラインでリアルタイムに非接触且つ非破壊で測定する技術として、例えば以下に示す特許文献１には、周波数１０5乃至１０12［Ｈｚ］にある所定周波数の電磁波を照射する電磁波送信手段と、透過または反射した電磁波を受信する電磁波受信手段と、送受信波の位相と振幅を測定する測定器と、測定器での測定から振幅変化と位相差を検出し解析する解析装置とを備える水分量測定装置が開示されている。この解析装置は、予め既知の水分量の物質について振幅変化と位相差を検出し、振幅変化と位相差の関係を直線近似した直線の傾きを求め、さらに、傾きと水分量との関係を表す検量線を求めて記憶手段に記憶しておき、処理ラインを流れる物質に所定周波数の電磁波を照射して検出した振幅変化と位相差との組の複数の関係を直線近似した直線の傾きを求め、当該傾きを検量線に当てはめることにより水分量を検出するものである。
【０００３】
この技術は、直線近似できるような結果が得られ、検量線も線形性を維持できる周波数での測定が前提となっており、感度が悪いという問題点がある。
【先行技術文献】
【特許文献】
【０００４】
【特許文献１】特許第６２５３０９６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００５】
よって、本発明の目的は、一側面としては、ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質が含まれるか否かを精度良く測定できるようにするための技術を提供することである。
【課題を解決するための手段】
【０００６】
本発明に係る解析装置は、ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質が含まれるか否かを解析する解析装置であって、（Ａ）ある物質によって共振周波数に変化が生ずる測定系について、測定系本来の共振周波数を含む所定帯域内の特定の周波数における伝送特性、反射特性又はインピーダンス特性を表すデータを、ある物質の異なる複数の状態について取得する手段と、（Ｂ）取得された複数回分の伝送特性、反射特性又はインピーダンス特性を表すデータから、所定の特徴量を算出する算出手段と、（Ｃ）他の特定の物質の濃度若しくは含有量と所定の特徴量との関係又は他の物質が含まれているか否かと所定の特徴量との関係を表す関係データを格納したデータ格納部に格納された関係データに基づき、算出された所定の特徴量に対応する他の特定の物質の濃度若しくは含有量、又は他の物質が含まれているか否かを決定する手段とを有する。
【発明の効果】
【０００７】
一側面によれば、ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質が含まれるか否かを、精度良く測定できるようになる。
【図面の簡単な説明】
【０００８】
【図１】図１は、システム全体の概要を示す図である。
【図２】図２は、被測定物による共振周波数のシフトを説明するための図である。
【図３】図３（ａ）は、振幅変化量を説明するための図であり、図３（ｂ）は、位相変化量を説明するための図である。
【図４】図４は、解析装置の機能ブロック図である。
【図５】図５は、解析装置の処理フローを示す図である。
【図６】図６は、振幅変化量と位相変化量とで張られる平面における２次関数のフィッティングを説明するための図である。
【図７】図７は、複素平面における２次関数のフィッティングを説明するための図である。
【図８】図８は、検量線を説明するための図である。
【図９】図９は、複素平面における円フィッティングを説明するための図である。
【図１０】図１０は、コンピュータの機能ブロック図である。
【図１１】図１１（ａ）は、狭帯域アンテナの一例としてマイクロストリップのライン型アンテナを用いた場合におけるシステム全体の概要を示す図であり、（ｂ）は、狭帯域アンテナの一例としてパッチ型アンテナを用いた場合におけるシステム全体の概要を示す図である。
【図１２】図１２（ａ）は、狭帯域アンテナの一例としてマイクロストリップのライン型アンテナを用いた場合におけるシステム全体の概要を示す図であり、（ｂ）は、狭帯域アンテナの一例としてパッチ型アンテナを用いた場合におけるシステム全体の概要を示す図である。
【図１３】図１３（ａ）及び（ｂ）は、ベクトルネットワークアナライザを用いない場合の測定系の例を示す図である。
【発明を実施するための形態】
【０００９】
上で述べた特許文献１では、例えば米が処理ラインを空気と共に流れてくる際に、所定の周波数の電磁波を照射して、その透過電磁波または反射電磁波についての振幅変化及び位相差を複数時点で検出し、振幅変化及び位相差で張られる平面において直線近似した際の傾きにて、米の水分量を特定していた。
【００１０】
この技術においては、ＲＦ波、マイクロ波、ミリ波といった周波数帯の中で任意の周波数を選択できることになっているが、実際は、直線近似できるような測定結果が得られる周波数、そして水分量と直線の傾きとが線形関係を維持できるような周波数を用いて測定を行っている。
【００１１】
このような前提の下の測定では、感度が悪くて測定精度が落ちてしまうという問題がある。
【００１２】
そこで、本実施の形態では、共振を用いて測定系における感度を高めることで、測定精度を上げるものとする。但し、このような測定の結果は、非線形性を有しているので、それに対応した解析処理を行うことになる。共振を用いた測定系としては、干渉計を用いる場合や、狭帯域アンテナを用いる場合などがある。
【００１３】
なお、本実施の形態では、特許文献１のような粉状、粒状又は不定形のバラ状の物質の水分量だけではなく、ある物質（主たる物質）に含まれる他の物質の濃度又は含有量を測定する場合にも対応する。さらに、ある物質（主たる物質）に他の物質（糖分、塩分、アルコール、脂肪、二酸化炭素、尿素など。何らかの異物の場合もある。）が含まれているか否かを判定することにも対応する。
【００１４】
干渉計を用いる例を図１に示す。
【００１５】
本実施の形態に係るシステムは、干渉計１００と、測定器であるベクトルネットワークアナライザ２００と、解析装置３００とを有する。
【００１６】
干渉計１００は、ベクトルネットワークアナライザ２００の出力に接続された同軸ケーブル１０１と、Ｔ分岐回路１０３と、マイクロストリップライン１０４と、ベクトルネットワークアナライザ２００の入力に接続される同軸ケーブル１０２とを有する。
【００１７】
ベクトルネットワークアナライザ２００から出力された信号Ｓ１は、Ｔ分岐回路１０３において、マイクロストリップライン１０４を通らない信号Ｓ２と、マイクロストリップライン１０４への信号Ｓ３とに分岐される。マイクロストリップライン１０４へ分岐した信号Ｓ３は、マイクロストリップライン１０４を伝搬し、他端の短絡器で反射して反射信号Ｓ４がＴ分岐回路１０３へ戻ってくる。そして、Ｔ分岐回路１０３では、信号Ｓ２と信号Ｓ３とが合成され、合成された信号Ｓ５が、ベクトルネットワークアナライザ２００に入力される。
【００１８】
ベクトルネットワークアナライザ２００は、一般的なものであって、被測定物の伝送特性、反射特性、インピーダンス特性などを測定するものである。なお、これらの特性を表すデータとしては、振幅及び位相の組み合わせと、実部及び虚部の組み合わせとがある。両者は実質的に同じものを表しており、一方から他方に変換可能である。
【００１９】
Ｔ分岐回路１０３の分岐点からマイクロストリップライン１０４の短絡器までの電気長Ｌeは、物理的な長さＬと、基板及びケーブルやパターンなどの実効的な誘電率εrとから、Ｌe＝Ｌ／｛εr｝0.5で決まる。実効的な誘電率εrは、マイクロストリップライン１０４上に載せられた被測定物によっても変化する。
【００２０】
この電気長Ｌeに基づき、信号の波長λ／２＝２＊Ｌeの時に共振状態となり、信号の振幅が極小になる。このようすを図２に示す。
【００２１】
図２において、横軸は周波数を表し、縦軸は信号の振幅を表す。図２に示すように、電気長Ｌeに基づき、周期的に振幅極小のノッチが現れる。カーブＡは、マイクロストリップライン１０４上に被測定物を載せていない状態での特性を表しており、カーブＢは、マイクロストリップライン１０４上に被測定物を載せた場合の特性を表している。このように、被測定物によって実効誘電率が変化したため、共振周波数がずれて、ノッチがシフトしていることが分かる。
【００２２】
本実施の形態では、このような現象を用いて測定を行う。なお、図３（ａ）に示すように、振幅については、干渉計１００の共振周波数を含む所定帯域Ｗにおける特定の観測周波数ｆvにおけるカーブＢの値とカーブＡの値との差である振幅変化量をベースに解析を行う。また、図３（ｂ）に示すように、位相についても、特定の観測周波数ｆvにおける、マイクロストリップライン１０４上に被測定物を載せない場合のカーブＣの値とマイクロストリップライン１０４上に被測定物を載せた場合のカーブＤの値の差である位相変化量をベースに解析を行う。
【００２３】
なお、所定帯域Ｗは、共振を用いて感度を高められる範囲であって、測定系によっても異なるし、被測定物によっても異なる。観測周波数ｆvは、例えば共振周波数以外で、振幅変化量及び位相変化量が十分大きくなるような周波数として、所定帯域Ｗ内で指定される。
【００２４】
このような振幅変化及び位相変化の組み合わせを、実部及び虚部の組み合わせとして測定する場合もある。
【００２５】
図１に示したような測定系では、被測定物は、マイクロストリップライン１０４上に載せられて測定される。例えば、液体であれば袋などの容器に入れて載せることになる。この際、載せ方を変えたり、被測定物の体積や量を変えたりすることで、被測定物の状態を変化させて複数回測定する。
【００２６】
そして、解析装置３００は、ベクトルネットワークアナライザ２００から出力された複数回分の測定結果に基づき、ある物質に含まれる他の物質の濃度又は含有量などを決定する。
【００２７】
解析装置３００の機能構成例を、図４に示す。
【００２８】
解析装置３００は、データ取得部３０１と、測定データ格納部３０３と、特徴量算出部３０５と、特徴量格納部３０７と、決定部３０９と、検量線データ格納部３１１と、出力部３１３とを有する。
【００２９】
データ取得部３０１は、ベクトルネットワークアナライザ２００から測定データを取得し、測定データ格納部３０３に格納する。特徴量算出部３０５は、測定データ格納部３０３に格納されている測定データから、特徴量を算出する処理を実行し、算出された特徴量を特徴量格納部３０７に格納する。本実施の形態では、複数回分の測定データから、２次関数フィッティング又は円フィッティングを行って、フィッティングカーブの係数を算出する。
【００３０】
検量線データ格納部３１１は、例えば、フィッティングカーブの係数毎に、検量線のデータを格納している。検量線も、例えば２次関数である。決定部３０９は、検量線データ格納部３１１に格納されているデータに基づき、特徴量格納部３０７に格納されている特徴量に対応する他の物質の含有量若しくは濃度、又はある物質に他の物質が含まれるか否か等を決定する。出力部３１３は、例えば表示装置などの出力装置に、決定部３０９の処理結果を出力する。
【００３１】
次に、解析装置３００の処理内容について、図５乃至図９を用いて説明する。
【００３２】
まず、データ取得部３０１は、ベクトルネットワークアナライザ２００などの測定器から、被測定物の異なる状態について、複数回分の測定データを取得し、測定データ格納部３０３に格納する（図５：ステップＳ１０１）。
【００３３】
特徴量算出部３０５は、測定データ格納部３０３に格納されたデータから、特徴量として、２次関数のフィッティングカーブの係数を算出し、特徴量格納部３０７に格納する（ステップＳ１０３）。
【００３４】
なお、測定器から、振幅変化量及び位相変化量、又はその実部及び虚部とを得るものとして説明しているが、測定器からの測定データからこれらの値をデータ取得部３０１又は特徴量算出部３０５等により算出するようにしてもよい。
【００３５】
振幅変化量及び位相変化量のデータを取得した場合について、図６を用いて説明する。図６において、縦軸は振幅変化量［ｄＢ］を表し、横軸は位相変化量［ｄｅｇ］を表している。ここでは、周波数３００ＭＨｚにおける、高純度水と、塩素を不純物として混ぜた不純物水溶液との測定結果を表している。このように、振幅変化量と位相変化量とで張られる平面において測定データをプロットすると、いずれについても、それぞれ２次関数のカーブにフィッティングされる。図６の例では、高純度水の場合、カーブａにフィッティングされ、不純物水溶液についてはカーブｂにフィッティングされる。フィッティングの方法は、最小自乗法などの周知の方法を用いれば良い。
【００３６】
なお、高純度水の場合、位相変化量−５０°付近においては、変化が急峻でデータが欠落してしまっている。
【００３７】
そうすると、図６に示されているように、ｙ＝ａｘ2＋ｂｘ＋ｃの形で係数ａ、ｂ及びｃが、それぞれのカーブについて算出される。
【００３８】
なお、位相については、一般的に測定器の出力は−１８０°から＋１８０の範囲になるが、測定データによっては−１８０°付近と＋１８０°付近とに分離してしまう場合もある。この場合には、例えば−１７０°を＋１９０°に変換する折り返し処理を行うことで、データを一方にまとめてフィッティングを行いやすくする。
【００３９】
一方、実部及び虚部のデータを取得した場合については、図７を用いて説明する。図７において、縦軸は虚部（Imaginary）を表し、横軸は実部（Real）を表す。ここでは、周波数３００ＭＨｚにおける、高純度水と不純物水溶液との測定結果を表している。このように、実部及び虚部で張られる平面（すなわち複素平面）において測定データをプロットすると、いずれについても、それぞれ２次関数のカーブにフィッティングされる。図７の例では、高純度水の場合、カーブｃにフィッティングされ、不純物水溶液についてはカーブｄにフィッティングされる。
【００４０】
そうすると、図７の場合でも、ｙ＝ａｘ2＋ｂｘ＋ｃの形で係数ａ、ｂ及びｃが、それぞれのカーブについて算出される。
【００４１】
なお、実部及び虚部のデータの場合には、折り返し処理が不要であるという利点がある。
【００４２】
次に、決定部３０９は、検量線データ格納部３１１に格納されている検量線のデータに基づき、特徴量格納部３０７に格納されている特徴量のデータ（すなわちフィッティングカーブの係数）に対応する濃度又は含有量などを決定する（ステップＳ１０７）。
【００４３】
検量線データ格納部３１１には、例えば図８に示すようなデータが格納されている。図８の例では、横軸が濃度［ｐｐｍ］を表しており、左の縦軸は係数ａの値を表し、右の縦軸は係数ｂ又はｃ／１００の値を表している。なお、図８に示したものは、図６に示した振幅変化量と位相変化量についての検量線のデータである
【００４４】
本実施の形態では、各係数について２次関数でフィッティングされた検量線が引かれている。すなわち、濃度が既知のサンプルについて上記の測定系で測定し、上で述べたように振幅変化量と位相変化量についてのフィッティングカーブの係数を算出し、係数と濃度との関係を２次関数のカーブにフィッティングさせれば、検量線が得られる。
【００４５】
３つの係数があるので３本の検量線が得られるが、いずれか１つ信頼性の高い検量線を特定してその検量線から濃度を決定するようにしてもよい。
【００４６】
一方、３つの係数それぞれについて検量線から濃度を特定した上で、平均などによって濃度を決定しても良い。加重平均するようにしてもよい。
【００４７】
測定データが実部及び虚部であっても、同様に検量線を用意すれば、濃度を特定できる。また、濃度ではなく含有量であっても同様である。
【００４８】
なお、ある物質に他の物質が含まれているか否かについて判定する場合には、例えば各係数について、他の物質が含まれている場合の値域又は他の物質が含まれていない場合の値域を特定して、そのデータを検量線データ格納部３１１に格納しておく。そして、例えば、全ての係数又は特定の係数について、他の物質が含まれている場合の値域に入っているか否か又は他の物質が含まれていない場合の値域に入っているか否かを確認すれば良い。
【００４９】
場合によっては、物質Ａに、物質Ｂと物質Ｃが入っている可能性がある場合、物質Ａのみのケース、物質Ａに物質Ｂが入っているケース、物質Ａに物質Ｃが入っているケース、物質Ａに物質Ｂ及び物質Ｃが入っているケースとについて、各係数の値域を特定しておけば、上記のようなケースを区別できる場合もある。
【００５０】
出力部３０９は、決定部３０９による処理結果を、出力装置である表示装置などに出力する（ステップＳ１０７）。なお、ネットワークに接続されている他のコンピュータなどに出力する場合もある。
【００５１】
なお、図７では、２次関数でのフィッティングについて述べたが、被測定物の状態（配置状態、量、空気との粗密度合いなど）の変化が大きい場合等には、非線形性が強くなるため、円フィッティングが好ましい場合もある。
【００５２】
図９に円フィッティングの例を示す。図９において、縦軸は虚部を表し、横軸は実部を表しており、カーブｅは高純度水の円弧を表しており、カーブｆは不純物水溶液の円弧を表している。円フィッティングの場合には、ｘ2＋ｙ2＋Ａｘ＋Ｂｙ＋Ｃ＝０という関数における係数Ａ、Ｂ及びＣを算出する。図９の場合、カーブｅについてＡ＝−０．７２６、Ｂ＝０．８３０、Ｃ＝０．１２２であり、カーブｆについてはＡ＝−１．２２６、Ｂ＝１．４４８、Ｃ＝０．３０２である。
【００５３】
当然、円フィッティングの係数Ａ、Ｂ及びＣを算出する場合には、検量線データ格納部３１１には、係数Ａ、Ｂ及びＣのそれぞれについて検量線データを用意する。また、ある物質に他の物質が含まれているか否かについて判定する場合には、例えば各係数について、他の物質が含まれている場合の値域又は他の物質が含まれていない場合の値域を特定して、そのデータを検量線データ格納部３１１に格納しておく。
【００５４】
このようにすれば、被測定物の状態が変化しても、共振により感度を向上させて精度良く測定できるようになる。
【００５５】
以上本発明の実施の形態を説明したが、本発明はこれに限定されるものではない。
【００５６】
例えば、図４に示した機能ブロック構成は一例であって、プログラムモジュール構成とは一致しない場合もある。また、解析装置３００は、１台の装置ではなく複数台の装置にて実現される場合もある。解析装置３は、測定器と一体化されることもある。
【００５７】
また、干渉計１００を用いる測定系についても一例であって、被測定物によって共振周波数が変化することを利用して共振周波数付近の周波数にて測定できる他の測定系を採用しても良い。
【００５８】
図１に示したような干渉計１００に代わって、図１１及び図１２に示すような狭帯域アンテナを用いるようにしてもよい。図１１及び図１２の例で狭帯域アンテナ４００は、マイクロストリップのライン型アンテナ（図１１（ａ）及び図１２（ａ））及びパッチ型アンテナ（図１１（ｂ）及び図１２（ｂ））である。このような狭帯域アンテナ４００は、共振の原理に基づくものであり、特定の周波数（共振周波数に対応）を含む狭い帯域で電磁波を放出し、そのような帯域の電磁波を受信する。図１１（ａ）及び（ｂ）では、狭帯域アンテナ４００をベクトルネットワークアナライザ２００の１つのポートに接続する反射型の例を示す。すなわち狭帯域アンテナ４００の共振周波数付近の周波数の電磁波を被測定物に放射して、狭帯域アンテナ４００で被測定物からの反射波を受信して、ベクトルネットワークアナライザ２００が反射特性を表すデータを測定する。一方、図１２（ａ）及び（ｂ）の例では、狭帯域アンテナ４００をベクトルネットワークアナライザ２００の２つのポートに接続する干渉型の例を示す。これは、狭帯域アンテナ４００から被測定物に放射される電磁波と、被測定物から反射された電磁波との干渉の状態を、上で述べた干渉計と同様にベクトルネットワークアナライザ２００にて測定するものである。
【００５９】
また、被測定物の体積や量を変えたりすることで、被測定物の状態を変化させて複数回測定するが、袋などの容器に入れた被測定物を、例えばベルトコンベアなどの搬送器で移動させている途中で容器毎に測定を行うようにしてもよい。また、パイプ等に流れる被測定物についても、干渉計１００又は狭帯域アンテナ４００の配置を工夫することで、複数のタイミングで測定できるようになる。このようにすれば、オンラインでの解析が可能となる。
【００６０】
また、ベクトルネットワークアナライザ２００は、信号を出力し、伝送路又は物体を伝送した後の信号若しくはそれらからの反射信号の特性を、振幅及び位相、又は実数部及び虚数部として検出する装置である。上で述べた例では、ベクトルネットワークアナライザ２００を用いる例のみを示したが、図１３に示すように、ベクトルネットワークアナライザ２００に代わって、信号源、受信機及びそれに付随する装置の組み合わせを用いるようにしてもよい。図１３（ａ）では、受信機１にて信号源からの信号を受信すると共に、受信機２にてセンサまたは回路（共振回路など）からの反射信号を受信することで、それらの出力から振幅変化量及び位相変化量などを例えば解析装置３００にて算出するものである。同様に、図１３（ｂ）では、受信機１にて信号源からの信号を受信すると共に、受信機２にてセンサまたは回路からの出力である伝送信号を受信することで、それらの出力から振幅変化量及び位相変化量などを例えば解析装置３００にて算出するものである。
【００６１】
なお、上で述べた解析装置３００は、例えばコンピュータ装置であって、図１０に示すように、メモリ２５０１とＣＰＵ（Central Processing Unit）２５０３とハードディスク・ドライブ（ＨＤＤ：Hard Disk Drive）２５０５と表示装置２５０９に接続される表示制御部２５０７とリムーバブル・ディスク２５１１用のドライブ装置２５１３と入力装置２５１５とネットワークに接続するための通信制御部２５１７とがバス２５１９で接続されている。なお、ＨＤＤはソリッドステート・ドライブ（ＳＳＤ：Solid State Drive）などの記憶装置でもよい。オペレーティング・システム（ＯＳ：Operating System）及び本発明の実施の形態における処理を実施するためのアプリケーション・プログラムは、ＨＤＤ２５０５に格納されており、ＣＰＵ２５０３により実行される際にはＨＤＤ２５０５からメモリ２５０１に読み出される。ＣＰＵ２５０３は、アプリケーション・プログラムの処理内容に応じて表示制御部２５０７、通信制御部２５１７、ドライブ装置２５１３を制御して、所定の動作を行わせる。また、処理途中のデータについては、主としてメモリ２５０１に格納されるが、ＨＤＤ２５０５に格納されるようにしてもよい。本発明の実施例では、上で述べた処理を実施するためのアプリケーション・プログラムはコンピュータ読み取り可能なリムーバブル・ディスク２５１１に格納されて頒布され、ドライブ装置２５１３からＨＤＤ２５０５にインストールされる。インターネットなどのネットワーク及び通信制御部２５１７を経由して、ＨＤＤ２５０５にインストールされる場合もある。このようなコンピュータ装置は、上で述べたＣＰＵ２５０３、メモリ２５０１などのハードウエアとＯＳ及びアプリケーション・プログラムなどのプログラムとが有機的に協働することにより、上で述べたような各種機能を実現する。
【００６２】
なお、上で述べたような処理を実行することで用いられるデータは、処理途中のものであるか、処理結果であるかを問わず、メモリ２５０１又はＨＤＤ２５０５等の記憶装置に格納される。
【００６３】
以上述べた実施の形態をまとめると以下のようになる。
【００６４】
本実施の形態に係る解析装置は、ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質（単数の場合もあれば複数の場合もある）が含まれるか否かを解析する解析装置であって、（Ａ）ある物質によって共振周波数に変化が生ずる測定系について、測定系本来の共振周波数を含む所定帯域内の特定の周波数における伝送特性、反射特性又はインピーダンス特性を表すデータを、ある物質の異なる複数の状態について取得する手段と、（Ｂ）取得された複数回分の伝送特性、反射特性又はインピーダンス特性を表すデータから、所定の特徴量を算出する算出手段と、（Ｃ）他の特定の物質の濃度若しくは含有量と所定の特徴量との関係又は他の物質が含まれているか否かと所定の特徴量との関係を表す関係データを格納したデータ格納部に格納された関係データに基づき、算出された所定の特徴量に対応する他の特定の物質の濃度若しくは含有量、又は他の物質が含まれているか否かを決定する手段とを有する。
【００６５】
このように共振によって感度を向上させることで、精度良く、ある物質に含まれる他の特定の物質の濃度若しくは含有量、又はある物質に他の物質が含まれるか否かを特定できるようになる。
【００６６】
なお、上で述べた算出手段が、複数回分の伝送特性、反射特性又はインピーダンス特性を表すデータで特定される振幅変化量及び位相変化量のデータから、振幅変化量及び位相変化量で張られる平面におけるフィッティングカーブを算出し、当該フィッティングカーブの係数を所定の特徴量として算出するようにしてもよい。これによって、ある物質の特徴量を適切に抽出することができるようになる。なお、この場合には、２次関数フィッティングが好ましい。
【００６７】
また、上で述べた算出手段が、複数回分の伝送特性、反射特性又はインピーダンス特性を表すデータで特定される実数部及び虚数部のデータから、実数部及び虚数部で張られる平面におけるフィッティングカーブを算出し、当該フィッティングカーブの係数を所定の特徴量として算出するようにしてもよい。位相及び振幅ではなく、実部及び虚部に着目することで、位相に係る値域の変換を行わずに済む。この場合も、２次関数フィッティングが行われる場合もあるが、円フィッティングを行う方が好ましい場合もある。
【００６８】
以上述べた解析装置の処理をコンピュータに実行させるためのプログラムを作成することができて、そのプログラムは、様々な記憶媒体に記憶される。
【符号の説明】
【００６９】
１００  干渉計
２００  ベクトルネットワークアナライザ
３００  解析装置
３０１  データ取得部
３０３  測定データ格納部
３０５  特徴量算出部
３０７  特徴量格納部  
３０９  決定部
３１１  検量線データ格納部
３１３  出力部
【図１】
2020148713.tif 000003
【図２】
2020148713.tif 000004
【図３】
2020148713.tif 000005
【図４】
2020148713.tif 000006
【図５】
2020148713.tif 000007
【図６】
2020148713.tif 000008
【図７】
2020148713.tif 000009
【図８】
2020148713.tif 000010
【図９】
2020148713.tif 000011
【図１０】
2020148713.tif 000012
【図１１】
2020148713.tif 000013
【図１２】
2020148713.tif 000014
【図１３】
2020148713.tif 000015