(57)【要約】本発明は、送信器と受信器とを備え、材料が充填された共振器で誘電体材料の湿度及び／又は密度を決定するための方法に関する。課題を解決するために本方法によれば、送信器が信号を発し、充填された共振器の共振曲線を段階的に走査し、その際に異なる周波数（ｆｉ）でそれぞれ対応する信号強度（Ｕｉ）を測定し、充填された共振器の共振周波数（ｆｍ）と帯域幅（ＢＷｍ）を測定された点（ｆｉ／Ｕｉ）から求め、空の共振器と充填された共振器の共振周波数（ｆｒ０，ｆｒｍ）と帯域幅（ＢＷ０，ＢＷｍ）及び前記共振器の既知の較正係数（ａｆｒ１，ａｆｒ２，ｂｆｒ１，ｂｆｒ２，ｃｆｒ１，ｃｆｒ２，ａｆｂｗ，ａｆｂｗ，ｂｆｂｗ１，ｂｆｂｗ２，ｃｆｂｗ１，ｃｆｂｗ２）を含む第２の方程式系（Ｇ２）を解くことによって材料の湿度（ψ）及び／又は密度（σ）を計算する。本発明の目的は、密度に依存することなく迅速且つ正確に湿度を決定できるようにするための方法を提供することである。【選択図】図１
2007509328.tif 000002
【特許請求の範囲】
【請求項１】
送信器と受信器とを備えた共振器に充填された誘電体材料の湿度（ψ）及び／又は密度（ρ）を決定するための方法において、
送信器が信号を発し、
充填された共振器の共振曲線を段階的に走査し、その際に異なる周波数（ｆｉ）でそれぞれ付属する受信信号の信号強度（Ｕｉ）を測定し、
充填された共振器について測定された点（ｆｉ／Ｕｉ）から共振周波数（ｆｒｍ）と帯域幅（ＢＷｍ）を求め、
空の共振器と充填された共振器の共振周波数（ｆｒ０，ｆｒｍ）と帯域幅（ＢＷ０，ＢＷｍ）及び前記共振器の既知の較正係数（ａｆｒ１，ａｆｒ２，ｂｆｒ１，ｂｆｒ２，ｃｆｒ１，ｃｆｒ２，ａｆｂｗ１，ａｆｂｅ２，ｂｆｂｗ１，ｂｆｂｗ２，ｃｆｂｗ１，ｃｆｂｗ２）を含む第２の方程式系（Ｇ２）を解くことによって材料の湿度（ψ）及び／又は密度（ρ）を計算することを特徴とする方法。
【請求項２】
充填された共振器の帯域幅（ＢＷｍ）を求めるための点（ｆｉ／Ｕｉ）から、共振周波数（ｆｒｍ）、共振器品質（Ｑｍ）及び最大共振値（Ｕｒｍ）の値を求めてそこから帯域幅（ＢＷｍ）を計算するか、または遮断周波数（ｆａｍ，ｆｂｍ）を求めてそこから共振周波数（ｆｒｍ）及び帯域幅（ＢＷｍ）を計算することを特徴とする、請求項１に記載の方法。
【請求項３】
下方閾値（Ｕａ）を算出し、信号強度（Ｕｉ）が当該閾値より高い範囲で第２の走査過程をより小さい段階幅で行うことを特徴とする、請求項１又は２に記載の方法。
【請求項４】
共振曲線の走査を等間隔の段階（Δｆ１，Δｆ２）で行うことを特徴とする、請求項１、２又は３に記載の方法。
【請求項５】
送信器を一定の強度で運転することを特徴とする、前記請求項のいずれか１項に記載の方法。
【請求項６】
共振器の遮断周波数（ｆａ，ｆｂ）を求めるために、
受信器の信号強度が最大値（Ｕｍａｘ）となる測定点（ｆｉ／Ｕｉ）を決定し、これから出発して閾値（Ｕｇ）を計算し、
正及び負の勾配部分について、信号値（Ｕｉ，Ｕｉ＋１）がそれぞれ閾値（Ｕｇ）の下方及び上方にある、隣接する２個の点（ｆｉ／Ｕｉ，ｆｉ＋１／Ｕｉ＋１）を決定し、そこからそれぞれ隣接する点（ｆｉ／Ｕｉ，ｆｉ＋１／Ｕｉ＋１）の間を補間することによって第１の遮断周波数（ｆａ）と第２の遮断周波数（ｆｂ）を計算することを特徴とする、請求項２〜５のいずれか１項に記載の方法。
【請求項７】
閾値（Ｕｇ）が最大信号値（Ｕｍａｘ）に対して３ｄＢ減衰した値に相当することを特徴とする、請求項６に記載の方法。
【請求項８】
共振器の共振周波数（ｆｒ）、共振器品質（Ｑ）及び最大共振値（Ｕｒ）の値を決定するために、
３個の点（ｆｉ／Ｕｉ）を任意的及び／又はランダムに選択し、
これら３個の点（ｆｉ／Ｕｉ）に対して有効な分析的共振曲線の３つの方程式からなる第１の方程式系（Ｇ１）を解いてこれらの値（ｆｒ，Ｑ，Ｕｒ）を求めることを特徴とする、請求項２〜５のいずれか１項に記載の方法。
【請求項９】
共振器の共振周波数（ｆｒ）、共振器品質（Ｑ）及び最大共振値（Ｕｒ）の値を決定するために、
３の整数倍で６個以上の点（ｆｉ／Ｕｉ）を任意的及び／又はランダムに選択し、それらの点を大きさの等しい３つのグループ（Ｍ１，Ｍ２，Ｍ３）に分割し、
それぞれ異なるグループ（Ｍ１，Ｍ２，Ｍ３）から１個の点（ｆｉ／Ｕｉ）を取り出した３個の点（ｆｉ／Ｕｉ）の各々の組合せについて、これら３個の点（ｆｉ／Ｕｉ）に対して有効な分析的共振曲線の３つの方程式からなる第１の方程式系（Ｇ１）を解いて値（ｆｒｋ，Ｑｋ，Ｕｒｋ）を求め、
各々の値（ｆｒ，Ｑ，Ｕｒ）について、それぞれの組合せで算出された値（ｆｒｋ，Ｑｋ，Ｕｒｋ）の平均値を形成することを特徴とする、請求項２〜５のいずれか１項に記載の方法。
【請求項１０】
点（ｆｉ／Ｕｉ）を任意的及び／又はランダムに選択するための条件として、選択されるべき点（ｆｉ／Ｕｉ）の信号値（Ｕｉ）が３ｄＢ減衰された最大信号値（Ｕｍａｘ）より大きいことを特徴とする、請求項８又は９に記載の方法。
【請求項１１】
第２の方程式系（Ｇ２）が、湿度（ψ）及び密度（ρ）と、所定の湿度及び密度範囲における共振周波数（ｆｒ）及び共振器品質（Ｑ）又は共振周波数（ｆｒ）及び帯域幅（ＢＷ）との相関を高い近似性で表すことを特徴とする、前記請求項のいずれか１項に記載の方法。
【請求項１２】
第２の方程式系（Ｇ２）が非直線的であることを特徴とする、前記請求項のいずれか１項に記載の方法。
【請求項１３】
送信器による走査をマイクロ波の範囲にまで行うことを特徴とする、前記請求項のいずれか１項に記載の方法。
【請求項１４】
受信信号の測定のために受信器の電圧値又は電流値を利用することを特徴とする、前記請求項のいずれか１項に記載の方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、送信器と受信器とを備えた共振器に充填された誘電体材料の湿度及び／又は密度を決定するための方法に関する。
【０００２】
複素比誘電率εｒ＝ε’ｒ−ｉε”ｒによって表される材料の誘電特性は、多孔性材料の場合、湿度と密度によって影響される。湿度と密度は、空気で満たされた空の共振器のスカラーパラメータ（ｆｒ０及びＱ０）とは異なり、材料を充填した共振器のスカラーパラメータである共振周波数ｆｒと共振器のＱ値は
【数１】
2007509328.tif 000003
及び
【数２】
2007509328.tif 000004
のように湿度と密度によって影響を受ける。ここで、１≦ε’ｒ≦ε’ｒｍａｘ及び０≦ε”ｒ≦ε”ｒｍａｘであり、ε’ｒｍａｘ及びε”ｒｍａｘは、それぞれの材料について指定された湿度範囲及び密度範囲での最大値を表している。
【背景技術】
【０００３】
従来の技術において、材料を充填した共振器の共振特性を利用して粒状材料の湿度や密度を測定するための種々の方法が知られている。
【０００４】
例えばＵＳ５６６６０６１により知られている、マイクロ波を用いて粒状材料の湿度を測定するための方法では、ある周波数範囲を上昇及び下降する間に閾値に電子的に反応し、それによって形成されたパルスの時間特性から誘電定数の成分を推定する。
【０００５】
従来の方法はいずれも比較的遅く、材料の密度に依存するか、または湿度を決定する際に比較的大きい誤りを示す。
【０００６】
【特許文献１】ＵＳ５６６６０６１
【発明の開示】
【発明が解決しようとする課題】
【０００７】
本発明の目的は、密度に依存することなく湿度を迅速且つ正確に決定することが可能な、冒頭で述べた形式の方法を提供することである。
【課題を解決するための手段】
【０００８】
上記課題は、本発明により請求項１に記載された特徴を有する方法によって解決される。
【０００９】
本発明の有利な実施形態が、サブクレームに記載されている。
【００１０】
本発明についての以下の説明において、共振幅を表す任意の量を共振器の帯域と考えることができるが、それぞれの定義、特に閾値の定義及び方程式でなされている相応の規定を優先すべきである。
【００１１】
共振曲線をディジタル記録することにより測定値を迅速に獲得することが可能であり、それに基づいて決定された湿度を、共振器内に動的に案内される材料と時間的・空間的に密接に対応させることが可能である。更に本発明では、送信器が信号を発し、充填された共振器の共振曲線を段階的に走査し、その際に異なる周波数でそれぞれ対応する受信信号の信号強度を測定し、測定された点から充填された共振器の共振周波数と帯域幅を求め、空の共振器と充填された共振器の共振周波数と帯域幅及び前記共振器の既知の較正係数を含む第２の方程式系を解くことによって材料の湿度及び／又は密度を計算するようにした。この方法により、密度に依存することなく材料の湿度を決定できる。しかも材料の密度を容易に求めることができる。
【００１２】
好適な実施形態においては、充填された共振器の帯域幅を求めるための前記点から、共振周波数、共振器品質及び最大共振値の値を求めてそこから帯域幅を計算するか、または遮断周波数を求めてそこから共振周波数及び帯域幅を計算する。
【００１３】
別の実施形態では、下方閾値を算出し、信号強度が当該閾値より高い範囲で、第２の走査過程をより小さい段階幅で行う。この場合、２パス法により精度は著しく高められるにもかかわらず、時間的コストは少なく抑えられる。特に第２の走査過程で共振ピークの範囲においてのみより小さい段階幅を適用する場合はそうである。
【００１４】
共振曲線の走査を等間隔段階で行うと有利である。走査の最も単純で速い形式は、等間隔段階からなる。第２の走査段階で段階幅を可変にすれば、その時間を短縮できる。
【００１５】
送信器を一定の強度で運転することが好ましい。送信器の信号強度が一定であると、受信器での測定値を調整又はスケーリングを行うことなく使用できる。
【００１６】
可能な実施形態において、共振器の遮断周波数を求めるために、受信器の信号強度が最大となる点を決定し、ここから出発して閾値を計算し、正及び負の勾配部分について、信号強度がそれぞれ閾値の下方及び上方にある、隣接する２個の点を決定し、それぞれ隣接する点の間を補間することによって第１の遮断周波数と第２の遮断周波数を計算する。閾値を挟む点対の間を補間することによって、湿度及び／又は密度を決定するための出発パラメータとして遮断周波数を求めるための迅速且つ単純な方法である。
【００１７】
有利な実施形態においては、閾値が最大信号値に対して３ｄＢ減衰した値に対応する。信号強度の最大測定値から出発して３ｄＢの減衰に相当する閾値を選択すると、解かれるべき方程式が非常に単純な形式になる。
【００１８】
別の実施形態においては、共振器の共振周波数、共振器品質及び最大共振値の値を求めるために、３個の点を任意的及び／又はランダムに選択して、これら３点に対して有効な分析的共振曲線の３つの方程式からなる第１の方程式系を解いてこれらの値を求める。３つの方程式からなり、従って完全に規定された第１の方程式系を解くことにより、補間法よりも迅速で誤りが少なく、湿度及び／又は密度を決定するための出発パラメータとして直接共振周波数、共振器品質及び共振振幅を求めることができる。この場合、３個の点から成るグループの各点について３つの方程式のいずれか１つが有効である。このようにすることによって特にノイズに影響された共振曲線をより正確に分析できる。
【００１９】
別の実施形態において、共振器の共振周波数、共振器品質及び最大共振値の値を決定するために、３の整数倍で６個以上の点を任意的及び／又はランダムに選択し、それらの点を大きさの等しい３つのグループに分け、それぞれ異なるグループから１個の点を取り出した３個の点の各々の組合せについて、これら３個の点に対して有効な分析的共振曲線の３つの方程式からなる第１の方程式系を解いて上記の値を求め、各々の値について、それぞれの組合せで算出された値の平均値を求める。存在する点集合からより多くの３点グループを形成し、すべてのグループについて出発パラメータの平均を求めれば、３個の点だけ用いた場合よりも正確な値が得られる。
【００２０】
有利な実施形態において、点を任意的及び／又はランダムに選択するための条件として、選択されるべき点の信号値が３ｄＢ減衰された最大信号値より大きいことである。３個の点又は３の倍数の個数の点を用いた方法でも、あらゆる測定値の最大値から３ｄＢ減衰した値より信号強度が大きい点のみを共振曲線から選択することが好ましい。なぜならば、このようにすると有意な点のみが採用されるからである。
【００２１】
解かれるべき第２の方程式系は、湿度及び密度と、共振周波数及び共振器品質の変化若しくは共振周波数及び帯域幅の変化との相関を高い近似性で表すように選択することが好ましい。
【００２２】
第２の方程式系が非直線的であることが好都合である。
【００２３】
好適な実施形態において、送信器による走査をマイクロ波の範囲にまで行う。
【００２４】
受信信号の測定のために受信器の電圧値又は電流値を利用すると有利である。共振曲線を求めるために、受信器に発生する電圧を用いることが好ましい。なぜならば電圧はバックカップリングなしで容易に測定できるからである。しかし受信回路での電流の強さを測定してもよい。
【００２５】
以下に簡単にするために測定値はＵｉで表すが、それは電圧のみを使用することを意味するものではない。
【００２６】
離散的な測定点から再構成される共振曲線は、実際の共振曲線にできるだけ近似したものにすぎない。共振曲線は、供給された周波数ｆに対する受信器における信号振幅Ｕの推移によって表される。ここでＵｒは共振最大値、Ｑは共振の品質である。
【数３】
2007509328.tif 000005
遮断周波数ｆａ及びｆｂは、信号強度が所定の値を上回るか、または下回ることになる周波数である。このために最大値から３ｄＢ減衰した値を選択することが好ましい。
【数４】
2007509328.tif 000006
遮断周波数と曲線パラメータの共振周波数ｆｒ及び共振器品質Ｑとの間には、
【数５】
2007509328.tif 000007
の関係が成り立つ。遮断周波数の間隔は帯域幅
【数６】
2007509328.tif 000008
として定義される。
【発明を実施するための最良の形態】
【００２７】
以下に、本発明を実施例に基づいて更に詳細に説明する。
【００２８】
図１には、本発明による方法によって記録された２つの共振曲線が重ねて示されている。右側の共振曲線は空の共振器、左側の共振曲線は材料を充填した共振器によるものである。
【００２９】
スタート周波数ｆｓｔａｒｔ１及びストップ周波数ｆｓｔｏｐ１との間で共振器を走査（掃引）することにより、共振曲線が離散的に得られる。スタート周波数は最大限変位した周波数
【数７】
2007509328.tif 000009
最大限変化した品質値
【数８】
2007509328.tif 000010
及び正規化された電圧比
【数９】
2007509328.tif 000011
から、次のように計算される。
【数１０】
2007509328.tif 000012
ストップ周波数はｆｒ＝ｆｒ０により、次式のようになる。
【数１１】
2007509328.tif 000013
【００３０】
掃引速度は走査点の数ｎ１等に依存しており、図２に示すように２パス法によって高めることができる。そのために最初に比較的少ない数の走査点で走査、即ち、より大きい周波数幅
【数１２】
2007509328.tif 000014
で走査し、次に第２段階において予めスタート周波数及びストップ周波数
【数１３】
2007509328.tif 000015
又は
【数１４】
2007509328.tif 000016
を求めた後で、これらの周波数の間で小さい周波数幅
【数１５】
2007509328.tif 000017
で走査する。
【００３１】
走査によって測定された共振曲線から共振器パラメータ、ｆｒ，Ｑ及びＵｒを、例えばｆａｆｂ法、３点法及び３ｋ点法で求めることができる。これ以外の方法も使用可能である。
【００３２】
図３は、ｆａｆｂ法を示している。これは測定された共振曲線から第１の遮断周波数ｆａと第２の遮断周波数ｆｂを直接求めることに基づいている。このために最初に３ｄＢ閾値ライン
【数１６】
2007509328.tif 000018
を計算するために最大電圧Ｕ＝Ｕｍａｘを有する点（走査点）（ｆ＝ｆｍａｘ／Ｕ＝Ｕｍａｘ）を決定する。次いで３ｄＢ閾値ラインのすぐ近傍にあるそれぞれ２個の点を求める。走査点ａ１とａ２及びｂ１とｂ２の間で直線的に補間することにより、第１及び第２の３ｄＢ遮断周波数ｆａ、ｆｂが次式のように得られる。
【数１７】
2007509328.tif 000019
【００３３】
【数１８】
2007509328.tif 000020
上式によりｆｒとＱを計算できる。次に共振電圧Ｕｒから次式が得られる。
【数１９】
2007509328.tif 000021
【００３４】
共振曲線を離散化したために、遮断周波数の計算は程度の差はあれ誤りを含んでいる。それは第１に３ｄＢ閾値ライン
【数２０】
2007509328.tif 000022
をＵｍａｘから導き出すことにより、第２に走査点ａ１とａ２及びｂ１とｂ２の間で直線的に補間することによる。走査率を高めれば誤りは減少するが、同時に走査時間も長くなり、掃引速度も低下する。この方法において非直線的な補間を行えば、精度を改善することが可能である。
【００３５】
図４に示す３点法は３つの走査点Ｕ１、Ｕ２、Ｕ３を任意的に、またはランダム生成器によって、好ましくは線Ｕｓ＝ｓ・Ｕｍａｘの上方で選択し、ここでｓ＜１／√２である。
【００３６】
第１の方程式系Ｇ１
【数２１】
2007509328.tif 000023
を解くことにより、共振パラメータｆｒ、Ｑ及びＵｒが次のように得られる。
【数２２】
2007509328.tif 000024
【００３７】
ｆａｆｂ法と比べた３点法の利点は、共振曲線の離散化による誤りが発生せず、共振曲線の走査に少数の点しか必要とされず、従って高い掃引速度が達成されることである。
【００３８】
共振曲線がノイズに影響されている場合は、３点法によっても誤りを含んだ共振パラメーターが求められる。ノイズの影響を低減するために、走査点が３個以上選択される。即ち、３・ｋ個（ｋ＝２、３、４…）の走査点が選択され、３つのグループＭ１、Ｍ２、Ｍ３（各グループは、Ｕｓ＝ｓ・Ｕｍａｘの上方にあるｋ個の点を含む）に分割される。図５は、この３ｋ点法を示している。３つのグループの各々から取り出したそれぞれ１個の点を互いに組み合わせることによって得られるｋ３個の組み合わせについて、上述した３点法の要領で第１の方程式系Ｇ１から共振パラメータを計算し、次いでｋ３個の部分結果の平均を出す。
【数２３】
2007509328.tif 000025
【００３９】
求められた共振パラメータから湿度及び密度を決定するために、利用する装置を既知の湿度及び密度の材料を用いて較正することが必要である。この較正は通常運転の枠外で行うことが有利である。
【００４０】
較正の際は、物質値である湿度ψ及び密度ρに共振器の電気値、即ち共振周波数ｆｒ、共振器品質Ｑ、及び帯域幅ＢＷ（ＢＷ＝ｆｒ／Ｑ）を対応させる。この物質値と電気値との関係は非直線的であり、定義された湿度及び密度範囲について、例えば次の第２の方程式系Ｇ２により十分な精度で表すことができる。
【数２４】
2007509328.tif 000026
【００４１】
Δｆｒ及びΔＢＷは、空の共振器と材料が充填された共振器との共振周波数及び帯域幅の差である。
【数２５】
2007509328.tif 000027
【００４２】
較正の目的は、それぞれの材料について十分に大きい数の較正値（Δｆｒ及びΔＢＷとこれらに対応している物質値ψ及びρ）から１２個の較正係数、即ちａｆｒ１、ａｆｒ２、ｂｆｒ１、ｂｆｒ２、ｃｆｒ１、ｃｆｒ２、ａｆＢＷ１、ａｆＢＷ２、ｂｆＢＷ１、ｂｆＢＷ２、ｃｆＢＷ１、ｃｆＢＷ２を決定することである。
【００４３】
このために測定された共振周波数の差Δｆｒ及び帯域幅の差ΔＢＷを、適切な参照法で求めた湿度及び密度の値に対応させる。湿度及び密度に依存して共振周波数及び帯域幅を表すこれらの較正値は、較正係数を求めるための基礎である。等しい湿度による較正値の回帰は、
【数２６】
2007509328.tif 000028
及び
【数２７】
2007509328.tif 000029
の形式で行われ、回帰曲線は原点を通らなければならない。なぜならば空の共振器（ρ＝０の空気）にとってΔｆｒ及びΔＢＷもゼロだからである。回帰により湿度に依存して密度の等しい較正点が得られる。これらの較正点により新たに
【数２８】
2007509328.tif 000030
及び
【数２９】
2007509328.tif 000031
の形式の回帰を行う。こうして求められた回帰係数ａｆｒ、ｂｆｒ、ｃｆｒ及びａＢＷ、ｂＢＷ、ｃＢＷの値を密度に対してプロットし、そこから二次回帰によって較正係数を決定する。
【００４４】
図６は、較正結果を曲線群の形式で示しているが、これは分析の目的にも用いられる。
【００４５】
測定された共振周波数の差Δｆｒ及び帯域幅の差ΔＢＷから、上記の第２の方程式系Ｇ２を解くことによって、当該材料の湿度ψ及び密度ρを計算する。
【００４６】
この場合、２個の実数根と２個の虚数根が生じる。対象とする湿度及び密度範囲に真の解が１個だけ存在するか否かは、Δｆｒ−ΔＢＷグラフにおける較正曲線の推移から求めることができる。Δｆｒ−ΔＢＷグラフでは、等しい密度及び湿度の曲線について、帯域幅の差が共振周波数の差との関係で示されている。図６の点Ａ、Ｂ、Ｃ及びＤで特徴づけられるこれらの曲線の推移が対象とする湿度及び密度範囲において連続していて一義的である場合、この範囲には真の解が１個しか存在しない。
【００４７】
第２の方程式系Ｇ２自体を解くために、反復法が用いられる。このために第２の方程式系Ｇ２はψについて解かれる。
【数３０】
2007509328.tif 000032
【００４８】
湿度・密度グラフにおける２つの方程式の交点から、ψ及びρについて求められた値が得られる。
【図面の簡単な説明】
【００４９】
【図１】走査された共振曲線を示す。
【図２】２パス法で走査された共振曲線を示す。
【図３】共振パラメータを決定するための第１の方法を示す。
【図４】共振パラメータを決定するための第２の方法を示す。
【図５】共振パラメータを決定するための第３の方法を示す。
【図６】較正曲線を平面状に示したグラフである。
【符号の説明】
【００５０】
ｆｒ                      一般的な場合における共振周波数
ｆｒ０                           空の共振器の共振周波数
ｆｒｍ                           充填された共振器の共振周波数
ｆｓｔａｒｔ１、ｆｓｔｏｐ１     第１の走査過程のスタート周波数及びストップ周波数
Ｕｍａｘ                         最大測定信号値
ａ、ｓ                           閾値ファクター
ｆｓｔａｒｔ２、ｆｓｔｏｐ２            第２の走査過程のスタート周波数とストップ周波数
ｆｍａｘ                         最大信号値が生じる周波数
（ｆａ１／Ｕａ１）、（ｆａ２／Ｕａ２）  第１の遮断周波数の隣接点
（ｆｂ１／Ｕｂ１）、（ｆｂ２／Ｕｂ２）  第２の遮断周波数の隣接点
（ｆ１／Ｕ１）…（ｆ３／Ｕ３）   選択された３個の点
Ｍ１、Ｍ２、Ｍ３                 点グループ
（ｆ１１／Ｕ１１）…（ｆ１４／Ｕ１４）  点グループＭ１の要素
（ｆ２１／Ｕ２１）…（ｆ２４／Ｕ２４）  点グループＭ２の要素
（ｆ３１／Ｕ３１）…（ｆ３４／Ｕ３４）  点グループＭ３の要素
Ａ、Ｂ、Ｃ、Ｄ                   対象となる湿度及び密度の範囲の限界
【図１】
2007509328.tif 000033
【図２】
2007509328.tif 000034
【図３】
2007509328.tif 000035
【図４】
2007509328.tif 000036
【図５】
2007509328.tif 000037
【図６】
2007509328.tif 000038
【国際調査報告】
2007509328.tif 000039
2007509328.tif 000040
2007509328.tif 000041
2007509328.tif 000042
2007509328.tif 000043
2007509328.tif 000044
2007509328.tif 000045
2007509328.tif 000046
2007509328.tif 000047
2007509328.tif 000048