JP3461284B2 - 赤外線ガス分析計における検量線の作成方法 - Google Patents
赤外線ガス分析計における検量線の作成方法Info
- Publication number
- JP3461284B2 JP3461284B2 JP12422098A JP12422098A JP3461284B2 JP 3461284 B2 JP3461284 B2 JP 3461284B2 JP 12422098 A JP12422098 A JP 12422098A JP 12422098 A JP12422098 A JP 12422098A JP 3461284 B2 JP3461284 B2 JP 3461284B2
- Authority
- JP
- Japan
- Prior art keywords
- calibration curve
- gas analyzer
- gas
- infrared gas
- concentration
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
- 238000011088 calibration curve Methods 0.000 title claims description 36
- 238000000034 method Methods 0.000 claims description 21
- 239000007789 gas Substances 0.000 description 43
- 238000005259 measurement Methods 0.000 description 16
- 238000001514 detection method Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000004868 gas analysis Methods 0.000 description 3
- 230000009897 systematic effect Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 230000006835 compression Effects 0.000 description 1
- 238000007906 compression Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 239000003085 diluting agent Substances 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000002834 transmittance Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/27—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands using photo-electric detection ; circuits for computing concentration
- G01N21/274—Calibration, base line adjustment, drift correction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/25—Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
- G01N21/31—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
- G01N21/35—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
- G01N21/3504—Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing gases, e.g. multi-gas analysis
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Engineering & Computer Science (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Description
【発明の詳細な説明】
【０００１】
【発明の属する技術分野】この発明は、例えばＣＯやＣ
Ｏ2 などの濃度を測定する赤外線ガス分析計における検
量線の作成方法に関する。 【０００２】 【従来の技術】例えば自動車排ガスに含まれるＣＯやＣ
Ｏ2 の濃度測定には、非分散型赤外線ガス分析計（ＮＤ
ＩＲ）が用いられるが、エンジンが急速に低エミッショ
ン化する近年においては、より高度な赤外線ガス分析計
が望まれている。 【０００３】前記赤外線ガス分析計の分析精度を左右す
る要因の一つに濃度目盛校正用検量線の適正な維持管理
がある。この検量線の作成は、一般的には、図３に示す
ような構成によって行われる。この図において、１はガ
ス分析部で、例えばＣＯ2 を測定する赤外線ガス分析計
２、この赤外線ガス分析計２の出力を適宜波形成形する
プリアンプ３、このプリアンプ３を経た赤外線ガス分析
計出力をＡＤ変換するＡＤ変換器４などよりなる。 【０００４】５は前記ガス分析部１に対して所定濃度の
校正用標準ガスＳＧを供給するガス分割器で、予め濃度
が判っている成分ガス（スパンガス）ＣＧとこれを適宜
の濃度に希釈するための希釈ガス（例えば窒素ガスまた
はエアー）ＤＧを入力とし、所定濃度の校正用標準ガス
ＳＧを出力するものである。 【０００５】６は装置全体を統括制御するＭＣＵ（Ｍａ
ｉｎ Ｃｏｎｔｒｏｌ Ｕｎｉｔ）で、例えばパソコン
などのコンピュータである。このＭＣＵ６は、ガス分析
部１やガス分割器５をインタフェースコントローラ（Ｉ
ＦＣ）７，８を介して制御するとともに、ガス分析部１
からの出力に基づいて濃度計算などを行う演算処理機能
を備えている。９はＭＣＵ６と接続された表示装置で、
例えばタッチパネルディスプレイである。また、１０は
例えばハードディスクなどの記憶媒体で、測定データを
記憶したり、表示装置９に表示する内容を収めたファイ
ルなどが複数設けられている。 【０００６】上記構成において、校正用の成分ガスＣＧ
と希釈ガスＤＧとをガス分割器５に供給して、適宜の濃
度の成分ガスＳＧを得て、この適宜濃度に調整された成
分ガスＳＧをガス分析部１の赤外線ガス分析計２に供給
して、赤外線ガス分析計２の出力信号（測定値）を複数
点（測定点として４つ以上）得て、これをＭＣＵ６に入
力し、前記採取された４以上の測定値を基にして最小二
乗法を用いて、検量線を４次以下の多項式で近似するよ
うにしている。 【０００７】すなわち、任意適宜濃度をｘとするとき、
測定値との関係がｆ（ｘ）で表すことができるとする。
そのとき、ｆ（ｘ）は、 ｆ（ｘ）＝ａ0 ＋ａ1 ｘ＋ａ2 ｘ2 ＋ａ3 ｘ3 ＋ａ4 ｘ
4 また、ｉ番目に調整された適宜濃度をｘ i とし、そのと
きの測定値をｙ i とする。このとき、ｙ i とｆ（ｘ i ）
との差Δε i は、 Δε i ＝ｙ i −ｆ（ｘ i ） となり、Δε i 2 は、下記（１）式で表される。 Δεi 2 ＝｛ｙ i −ｆ（ｘ i ）｝ 2 ＝｛ｙi −（ａ0 ＋ａ1 ｘi ＋ａ2 ｘi 2 ＋ａ3 ｘi 3 ＋ａ4 ｘi 4 ）｝2 ……（１）この（１）式を展開し、最小二乗法を用いて 前記係数ａ
0 〜ａ4 を求める。 このとき、ε i 2 の和、すなわち、 が最小になるように、下記（２）式に示す偏微分方程式
を解く。 ……（２） 【０００８】ところで、前記赤外線ガス分析計２の測定
原理は、下記（３）式に示すランバート−ベールの法則
（Ｌａｍｂｅｒｔ−Ｂｅｅｒ’ｓ ｌａｗ）に則ってい
ることはよく知られている。 Ｉ＝Ｉ0 ｅｘｐ（−μｃＬ） ……（３） ここで、Ｉ0 ：入射光強度、Ｉ：透過光強度、μ：吸収
係数、ｃ：測定対象ガスの濃度、Ｌ：ガス層の厚さであ
る。 【０００９】前記（３）式は、透過光強度の変化の割合
であるところの透過率が低濃度領域では１次式近似が可
能で、高濃度になるに伴って高次式へ展開できることを
示している。 【００１０】 【発明が解決しようとする課題】そして、例えば自動車
排ガス測定においては、赤外線ガス分析計２で用いる検
量線を４次以下の多項式で近似する必要があるが、測定
範囲が広範囲の場合、つまり、前記１次式近似可能領域
と高次（２次、３次、４次）式近似可能領域とを一つに
含む場合、１次式近似可能領域の検量線精度が悪くなっ
てしまう。 【００１１】図４は、前記検量線精度を説明するもの
で、この図において、実線で示す曲線４１は検量線であ
り、仮想線で示す曲線４２は、実際の測定によってえら
れた濃度を表す曲線である。また、符号Ｉは１次式近似
可能領域であり、IIは高次（２次、３次、４次）式近似
可能領域である。この図から、検量線４１は測定曲線４
２に対して、高濃度領域である高次式近似可能領域IIに
おいてはほぼ一致しているが、低濃度領域である１次式
近似可能領域Ｉにおいては大きく外れて検量線精度が悪
くなっており、低濃度領域ＩでのＥrr％（誤差％）が規
格値（±２．０％ＰＴ）を満足しないことがある。ここ
で、％ＰＴは読み値に対する誤差のことである。 【００１２】上述のような赤外線ガス分析計における検
量線４１の測定曲線４２に対する誤差は、赤外線ガス分
析計の測定原理上生ずる系統誤差ともいうべきものであ
る。 【００１３】この発明は、上述の事柄に留意してなされ
たもので、その目的は、赤外線ガス分析計における測定
原理上生ずる系統誤差を可及的に小さくし、測定精度を
向上させることができる検量線を得るための赤外線ガス
分析計における検量線の作成方法（以下、単に検量線の
作成方法という）を提供することである。 【００１４】 【課題を解決するための手段】上記目的を達成するた
め、この発明では、赤外線ガス分析計に濃度が異なるガ
スを供給して、そのときのガス濃度に対応する出力を４
以上の複数点採取し、この採取された４以上の測定値を
基にして最小二乗法を用いて、検量線を４次の多項式で
近似するようにした検量線の作成方法において、前記多
項式近似を行う場合、検量線作成データの各点におい
て、濃度の二乗の逆数を重み付けとして用いるようにし
ている。 【００１５】従来の検量線の作成方法によって得られる
検量線は、低濃度領域において誤差が大きかったが、こ
の発明の検量線の作成方法によれば、この低濃度領域に
おいて大きな重み付けが行われるため、低濃度領域にお
ける誤差が問題にならない程度に圧縮される。また、従
来の検量線の作成方法によって得られる検量線は、高濃
度領域においてはそれほど大きな誤差は生じてなく、こ
の発明の検量線の作成方法においてもそれほど大きな重
み付けは行われないため、この領域における誤差はほと
んど生ずることがない。つまり、この発明の検量線の作
成方法によれば、低濃度領域および高濃度領域における
誤差を可及的に小さく抑えることができる。 【００１６】 【発明の実施の形態】発明の実施の形態を図面を参照し
ながら説明する。この発明では、赤外線ガス分析計の検
量線を作成するのに、従来において用いていた（１）式
に代えて、下記（４）式を用いるのである。 Δεi 2 ＝ｗi ｛ｙi −（ａ0 ＋ａ1 ｘi ＋ａ2 ｘi 2 ＋ａ3 ｘi 3 ＋ａ4 ｘi 4 ）｝2 ……（４） ここで、ｗi ＝１／ｘ i 2 （ｘ i ：ｉ番目に調整された
適宜濃度）である。 【００１７】つまり、前記（１）式の右辺に重み係数ｗ
i として濃度ｘ i の二乗の逆数を乗ずるのである。 【００１８】このように、カーブフィッティングを行う
際、濃度ｘ i の二乗の逆数を用いることに想到した過程
は、次の通りである。すなわち、従来の検量線において
は、式（１）を用いてカーブフィッティングを行うよう
にしていたため、図４に示すように、検量線が低濃度領
域Ｉにおいて実際の測定結果と大きな誤差を生じてい
た。この誤差を低減する手法としては、前記（１）式の
右辺に何らかの重み付けを行うことが考えられる。その
場合、重み付けとして、濃度の逆数、濃度の二乗の
逆数、濃度の三乗の逆数の３種類が考えられる。 【００１９】図１は上記各重み付けによる検量線の誤差
の違いを示すものであり、図２は図１の低濃度部分を拡
大して示すものである。いずれの図においても、横軸は
濃度（％）、縦軸は誤差の大きさ（％ＰＴ）を表してい
る。また、図中、符号〜を付した曲線は、前記各重
み付け〜と対応している。 【００２０】前記両図から次のことがわかる。すなわ
ち、低濃度領域Ｉにおける誤差は、よりも、より
もというように改善される。そして、高濃度領域IIに
おける誤差は、逆によりも、よりもというよう
に改善される。これは、前記の手法、重み付けとして
濃度の逆数を用いた場合には、ランバート−ベールの法
則に従う低濃度のリニア領域が４次式検量線からずれて
しまうためである。 【００２１】そして、前記の手法（重み付けとして濃
度の二乗の逆数を用いた場合）、の手法（重み付けと
して濃度の三乗の逆数を用いた場合）では、低濃度領域
Ｉでの重み付けがより大きくなるので、その誤差は小さ
くなる。しかし、の手法では、その重み付けが大きす
ぎて、高濃度領域IIでの誤差が大きくなる。 【００２２】これに対して、前記の手法によれば、低
濃度領域Ｉおよび高濃度領域IIのいずれにおいても誤差
は規格値（±２．０％ＰＴ）以下であり、低濃度から高
濃度にわたって広範囲の領域で系統誤差が小さくなる。
つまり、濃度と出力との直線性（リニアリティ）が向上
する。 【００２３】上述のように、検量線を多項式近似によっ
て行う場合、検量線作成データの各点において、濃度の
二乗の逆数を重み付けとして用いることにより、測定可
能レンジを広くすることができ、例えば最低レンジの５
０倍といった広範囲で、所望の測定を精度よく行うこと
ができる。 【００２４】 【発明の効果】この発明においては、赤外線ガス分析計
に濃度が異なるガスを供給して、広範囲のガス濃度に対
応する出力を４以上の複数点採取し、この採取された４
以上の測定値を基にして最小二乗法を用いて、検量線を
４次以下の多項式で近似するようにした赤外線ガス分析
計における検量線の作成方法において、前記多項式近似
を行う場合、検量線作成データの各点において、濃度の
二乗の逆数を重み付けとして用いるようにしているの
で、赤外線ガス分析計における測定原理上生ずる系統誤
差を可及的に小さくすることができる。したがって、こ
の発明によれば、低濃度領域から高濃度領域まで広範囲
にわたって精度よく測定できる赤外線ガス分析計を得る
ことができる。
Ｏ2 などの濃度を測定する赤外線ガス分析計における検
量線の作成方法に関する。 【０００２】 【従来の技術】例えば自動車排ガスに含まれるＣＯやＣ
Ｏ2 の濃度測定には、非分散型赤外線ガス分析計（ＮＤ
ＩＲ）が用いられるが、エンジンが急速に低エミッショ
ン化する近年においては、より高度な赤外線ガス分析計
が望まれている。 【０００３】前記赤外線ガス分析計の分析精度を左右す
る要因の一つに濃度目盛校正用検量線の適正な維持管理
がある。この検量線の作成は、一般的には、図３に示す
ような構成によって行われる。この図において、１はガ
ス分析部で、例えばＣＯ2 を測定する赤外線ガス分析計
２、この赤外線ガス分析計２の出力を適宜波形成形する
プリアンプ３、このプリアンプ３を経た赤外線ガス分析
計出力をＡＤ変換するＡＤ変換器４などよりなる。 【０００４】５は前記ガス分析部１に対して所定濃度の
校正用標準ガスＳＧを供給するガス分割器で、予め濃度
が判っている成分ガス（スパンガス）ＣＧとこれを適宜
の濃度に希釈するための希釈ガス（例えば窒素ガスまた
はエアー）ＤＧを入力とし、所定濃度の校正用標準ガス
ＳＧを出力するものである。 【０００５】６は装置全体を統括制御するＭＣＵ（Ｍａ
ｉｎ Ｃｏｎｔｒｏｌ Ｕｎｉｔ）で、例えばパソコン
などのコンピュータである。このＭＣＵ６は、ガス分析
部１やガス分割器５をインタフェースコントローラ（Ｉ
ＦＣ）７，８を介して制御するとともに、ガス分析部１
からの出力に基づいて濃度計算などを行う演算処理機能
を備えている。９はＭＣＵ６と接続された表示装置で、
例えばタッチパネルディスプレイである。また、１０は
例えばハードディスクなどの記憶媒体で、測定データを
記憶したり、表示装置９に表示する内容を収めたファイ
ルなどが複数設けられている。 【０００６】上記構成において、校正用の成分ガスＣＧ
と希釈ガスＤＧとをガス分割器５に供給して、適宜の濃
度の成分ガスＳＧを得て、この適宜濃度に調整された成
分ガスＳＧをガス分析部１の赤外線ガス分析計２に供給
して、赤外線ガス分析計２の出力信号（測定値）を複数
点（測定点として４つ以上）得て、これをＭＣＵ６に入
力し、前記採取された４以上の測定値を基にして最小二
乗法を用いて、検量線を４次以下の多項式で近似するよ
うにしている。 【０００７】すなわち、任意適宜濃度をｘとするとき、
測定値との関係がｆ（ｘ）で表すことができるとする。
そのとき、ｆ（ｘ）は、 ｆ（ｘ）＝ａ0 ＋ａ1 ｘ＋ａ2 ｘ2 ＋ａ3 ｘ3 ＋ａ4 ｘ
4 また、ｉ番目に調整された適宜濃度をｘ i とし、そのと
きの測定値をｙ i とする。このとき、ｙ i とｆ（ｘ i ）
との差Δε i は、 Δε i ＝ｙ i −ｆ（ｘ i ） となり、Δε i 2 は、下記（１）式で表される。 Δεi 2 ＝｛ｙ i −ｆ（ｘ i ）｝ 2 ＝｛ｙi −（ａ0 ＋ａ1 ｘi ＋ａ2 ｘi 2 ＋ａ3 ｘi 3 ＋ａ4 ｘi 4 ）｝2 ……（１）この（１）式を展開し、最小二乗法を用いて 前記係数ａ
0 〜ａ4 を求める。 このとき、ε i 2 の和、すなわち、 が最小になるように、下記（２）式に示す偏微分方程式
を解く。 ……（２） 【０００８】ところで、前記赤外線ガス分析計２の測定
原理は、下記（３）式に示すランバート−ベールの法則
（Ｌａｍｂｅｒｔ−Ｂｅｅｒ’ｓ ｌａｗ）に則ってい
ることはよく知られている。 Ｉ＝Ｉ0 ｅｘｐ（−μｃＬ） ……（３） ここで、Ｉ0 ：入射光強度、Ｉ：透過光強度、μ：吸収
係数、ｃ：測定対象ガスの濃度、Ｌ：ガス層の厚さであ
る。 【０００９】前記（３）式は、透過光強度の変化の割合
であるところの透過率が低濃度領域では１次式近似が可
能で、高濃度になるに伴って高次式へ展開できることを
示している。 【００１０】 【発明が解決しようとする課題】そして、例えば自動車
排ガス測定においては、赤外線ガス分析計２で用いる検
量線を４次以下の多項式で近似する必要があるが、測定
範囲が広範囲の場合、つまり、前記１次式近似可能領域
と高次（２次、３次、４次）式近似可能領域とを一つに
含む場合、１次式近似可能領域の検量線精度が悪くなっ
てしまう。 【００１１】図４は、前記検量線精度を説明するもの
で、この図において、実線で示す曲線４１は検量線であ
り、仮想線で示す曲線４２は、実際の測定によってえら
れた濃度を表す曲線である。また、符号Ｉは１次式近似
可能領域であり、IIは高次（２次、３次、４次）式近似
可能領域である。この図から、検量線４１は測定曲線４
２に対して、高濃度領域である高次式近似可能領域IIに
おいてはほぼ一致しているが、低濃度領域である１次式
近似可能領域Ｉにおいては大きく外れて検量線精度が悪
くなっており、低濃度領域ＩでのＥrr％（誤差％）が規
格値（±２．０％ＰＴ）を満足しないことがある。ここ
で、％ＰＴは読み値に対する誤差のことである。 【００１２】上述のような赤外線ガス分析計における検
量線４１の測定曲線４２に対する誤差は、赤外線ガス分
析計の測定原理上生ずる系統誤差ともいうべきものであ
る。 【００１３】この発明は、上述の事柄に留意してなされ
たもので、その目的は、赤外線ガス分析計における測定
原理上生ずる系統誤差を可及的に小さくし、測定精度を
向上させることができる検量線を得るための赤外線ガス
分析計における検量線の作成方法（以下、単に検量線の
作成方法という）を提供することである。 【００１４】 【課題を解決するための手段】上記目的を達成するた
め、この発明では、赤外線ガス分析計に濃度が異なるガ
スを供給して、そのときのガス濃度に対応する出力を４
以上の複数点採取し、この採取された４以上の測定値を
基にして最小二乗法を用いて、検量線を４次の多項式で
近似するようにした検量線の作成方法において、前記多
項式近似を行う場合、検量線作成データの各点におい
て、濃度の二乗の逆数を重み付けとして用いるようにし
ている。 【００１５】従来の検量線の作成方法によって得られる
検量線は、低濃度領域において誤差が大きかったが、こ
の発明の検量線の作成方法によれば、この低濃度領域に
おいて大きな重み付けが行われるため、低濃度領域にお
ける誤差が問題にならない程度に圧縮される。また、従
来の検量線の作成方法によって得られる検量線は、高濃
度領域においてはそれほど大きな誤差は生じてなく、こ
の発明の検量線の作成方法においてもそれほど大きな重
み付けは行われないため、この領域における誤差はほと
んど生ずることがない。つまり、この発明の検量線の作
成方法によれば、低濃度領域および高濃度領域における
誤差を可及的に小さく抑えることができる。 【００１６】 【発明の実施の形態】発明の実施の形態を図面を参照し
ながら説明する。この発明では、赤外線ガス分析計の検
量線を作成するのに、従来において用いていた（１）式
に代えて、下記（４）式を用いるのである。 Δεi 2 ＝ｗi ｛ｙi −（ａ0 ＋ａ1 ｘi ＋ａ2 ｘi 2 ＋ａ3 ｘi 3 ＋ａ4 ｘi 4 ）｝2 ……（４） ここで、ｗi ＝１／ｘ i 2 （ｘ i ：ｉ番目に調整された
適宜濃度）である。 【００１７】つまり、前記（１）式の右辺に重み係数ｗ
i として濃度ｘ i の二乗の逆数を乗ずるのである。 【００１８】このように、カーブフィッティングを行う
際、濃度ｘ i の二乗の逆数を用いることに想到した過程
は、次の通りである。すなわち、従来の検量線において
は、式（１）を用いてカーブフィッティングを行うよう
にしていたため、図４に示すように、検量線が低濃度領
域Ｉにおいて実際の測定結果と大きな誤差を生じてい
た。この誤差を低減する手法としては、前記（１）式の
右辺に何らかの重み付けを行うことが考えられる。その
場合、重み付けとして、濃度の逆数、濃度の二乗の
逆数、濃度の三乗の逆数の３種類が考えられる。 【００１９】図１は上記各重み付けによる検量線の誤差
の違いを示すものであり、図２は図１の低濃度部分を拡
大して示すものである。いずれの図においても、横軸は
濃度（％）、縦軸は誤差の大きさ（％ＰＴ）を表してい
る。また、図中、符号〜を付した曲線は、前記各重
み付け〜と対応している。 【００２０】前記両図から次のことがわかる。すなわ
ち、低濃度領域Ｉにおける誤差は、よりも、より
もというように改善される。そして、高濃度領域IIに
おける誤差は、逆によりも、よりもというよう
に改善される。これは、前記の手法、重み付けとして
濃度の逆数を用いた場合には、ランバート−ベールの法
則に従う低濃度のリニア領域が４次式検量線からずれて
しまうためである。 【００２１】そして、前記の手法（重み付けとして濃
度の二乗の逆数を用いた場合）、の手法（重み付けと
して濃度の三乗の逆数を用いた場合）では、低濃度領域
Ｉでの重み付けがより大きくなるので、その誤差は小さ
くなる。しかし、の手法では、その重み付けが大きす
ぎて、高濃度領域IIでの誤差が大きくなる。 【００２２】これに対して、前記の手法によれば、低
濃度領域Ｉおよび高濃度領域IIのいずれにおいても誤差
は規格値（±２．０％ＰＴ）以下であり、低濃度から高
濃度にわたって広範囲の領域で系統誤差が小さくなる。
つまり、濃度と出力との直線性（リニアリティ）が向上
する。 【００２３】上述のように、検量線を多項式近似によっ
て行う場合、検量線作成データの各点において、濃度の
二乗の逆数を重み付けとして用いることにより、測定可
能レンジを広くすることができ、例えば最低レンジの５
０倍といった広範囲で、所望の測定を精度よく行うこと
ができる。 【００２４】 【発明の効果】この発明においては、赤外線ガス分析計
に濃度が異なるガスを供給して、広範囲のガス濃度に対
応する出力を４以上の複数点採取し、この採取された４
以上の測定値を基にして最小二乗法を用いて、検量線を
４次以下の多項式で近似するようにした赤外線ガス分析
計における検量線の作成方法において、前記多項式近似
を行う場合、検量線作成データの各点において、濃度の
二乗の逆数を重み付けとして用いるようにしているの
で、赤外線ガス分析計における測定原理上生ずる系統誤
差を可及的に小さくすることができる。したがって、こ
の発明によれば、低濃度領域から高濃度領域まで広範囲
にわたって精度よく測定できる赤外線ガス分析計を得る
ことができる。
【図面の簡単な説明】
【図１】重み付けによる検量線の誤差の違いを示す図で
ある。 【図２】前記図１における低濃度領域を拡大して示す図
である。 【図３】検量線を得るための装置の構成を概略的に示す
図である。 【図４】従来の検量線の作成方法の問題点を説明するた
めの図である。 【符号の説明】 ２…赤外線ガス分析計、ＳＧ…ガス。
ある。 【図２】前記図１における低濃度領域を拡大して示す図
である。 【図３】検量線を得るための装置の構成を概略的に示す
図である。 【図４】従来の検量線の作成方法の問題点を説明するた
めの図である。 【符号の説明】 ２…赤外線ガス分析計、ＳＧ…ガス。
─────────────────────────────────────────────────────
フロントページの続き
(58)調査した分野(Int.Cl.7，ＤＢ名)
G01N 21/00 - 21/01
G01N 21/17 - 21/61
G01J 3/00 - 3/52
G01J 4/00 - 4/04
G01J 9/00 - 9/00
G01N 21/62 - 21/74
G01N 27/00 - 27/72
ＰＡＴＯＬＩＳ
Claims (1)
- (57)【特許請求の範囲】 【請求項１】 赤外線ガス分析計に濃度が異なるガスを
供給して、広範囲のガス濃度に対応する出力を４以上の
複数点採取し、この採取された４以上の測定値を基にし
て最小二乗法を用いて、検量線を４次以下の多項式で近
似するようにした赤外線ガス分析計における検量線の作
成方法において、前記多項式近似を行う場合、検量線作
成データの各点において、濃度の二乗の逆数を重み付け
として用いることを特徴とする赤外線ガス分析計におけ
る検量線の作成方法。
Priority Applications (3)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12422098A JP3461284B2 (ja) | 1998-04-18 | 1998-04-18 | 赤外線ガス分析計における検量線の作成方法 |
| US09/244,307 US6202032B1 (en) | 1998-04-18 | 1999-02-04 | Method for forming a calibration line in an infrared gas analyzer |
| DE1999115018 DE19915018B4 (de) | 1998-04-18 | 1999-04-01 | Verfahren zum Herstellen einer Kalibrierungskurve in einem Infrarot-Gasanalysator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP12422098A JP3461284B2 (ja) | 1998-04-18 | 1998-04-18 | 赤外線ガス分析計における検量線の作成方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH11304695A JPH11304695A (ja) | 1999-11-05 |
| JP3461284B2 true JP3461284B2 (ja) | 2003-10-27 |
Family
ID=14879977
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP12422098A Expired - Lifetime JP3461284B2 (ja) | 1998-04-18 | 1998-04-18 | 赤外線ガス分析計における検量線の作成方法 |
Country Status (3)
| Country | Link |
|---|---|
| US (1) | US6202032B1 (ja) |
| JP (1) | JP3461284B2 (ja) |
| DE (1) | DE19915018B4 (ja) |
Families Citing this family (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| GB2386415B (en) * | 2002-03-12 | 2004-04-21 | Sira Ltd | Method and apparatus for classifying emissions |
| JP4648739B2 (ja) * | 2005-03-29 | 2011-03-09 | 株式会社堀場製作所 | 分析装置及び校正方法 |
| US8148691B1 (en) * | 2009-08-21 | 2012-04-03 | Jacob Y Wong | Calibration methodology for NDIR dew point sensors |
| DE102011109000B4 (de) | 2011-07-29 | 2014-05-15 | Abb Technology Ag | Optische Gasanalysatoreinrichtung mit Mitteln zumVerbessern des Dynamikbereichs der Messung |
| CN103616345B (zh) * | 2013-12-16 | 2016-03-30 | 四川旌科仪器制造有限公司 | 红外碳硫分析仪的含碳量测定方法 |
| JP6571476B2 (ja) * | 2015-06-29 | 2019-09-04 | 旭化成エレクトロニクス株式会社 | ガス濃度測定装置 |
| US20170100057A1 (en) * | 2015-10-12 | 2017-04-13 | Oriental System Technology Inc. | Portable ndir breath acetone measurement apparatus with sub-ppm accuracy |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US3957372A (en) * | 1974-12-19 | 1976-05-18 | United Technologies Corporation | Vehicle exhaust gas analysis system |
| US4499378A (en) * | 1982-03-09 | 1985-02-12 | Horiba, Ltd. | Infrared radiation gas analyzer |
| JP3279756B2 (ja) * | 1993-08-24 | 2002-04-30 | 東ソー株式会社 | 定量計算装置 |
| US5804695A (en) * | 1993-11-02 | 1998-09-08 | Horiba Instruments Incorporated | Gas dividing method and apparatus |
| US5475223A (en) * | 1994-04-25 | 1995-12-12 | Ford Motor Company | System for monitoring exhaust gas composition |
| DE29504088U1 (de) * | 1995-03-10 | 1996-07-11 | Palocz-Andresen, Michael, Dr.-Ing.habil., 20459 Hamburg | On-Board-Diagnose-/OBD/-Vorrichtung im Mikromaßstab zur kontinuierlichen Messung des Schadstoffaustrages aus Kraftfahrzeugen |
| TW491961B (en) * | 1996-08-27 | 2002-06-21 | Nippon Oxygen Co Ltd | Method for analyzing impurities in gas and its analyzer |
-
1998
- 1998-04-18 JP JP12422098A patent/JP3461284B2/ja not_active Expired - Lifetime
-
1999
- 1999-02-04 US US09/244,307 patent/US6202032B1/en not_active Expired - Lifetime
- 1999-04-01 DE DE1999115018 patent/DE19915018B4/de not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| DE19915018B4 (de) | 2005-10-06 |
| DE19915018A1 (de) | 1999-10-28 |
| JPH11304695A (ja) | 1999-11-05 |
| US6202032B1 (en) | 2001-03-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Ho et al. | The UWO contribution to the NIST aerodynamic database for wind loads on low buildings: Part 1. Archiving format and basic aerodynamic data | |
| Driver et al. | The Millennium Galaxy Catalogue: the space density and surface-brightness distribution (s) of galaxies | |
| US20180321208A1 (en) | Determining the net emissions of air pollutants | |
| Liang et al. | Field comparison of electrochemical gas sensor data correction algorithms for ambient air measurements | |
| Bouwens et al. | Cloning Hubble Deep Fields. I. A model-independent measurement of galaxy evolution | |
| Horvath et al. | The relationship between atmospheric light scattering coefficient and visibility | |
| Barreira et al. | Response approach to the matter power spectrum covariance | |
| CN103175759A (zh) | 基于多种地基遥感技术获取城市气溶胶复折射指数的方法 | |
| US20070180936A1 (en) | Comprehensive particulate matter measurement system and method for using the same | |
| Kushner | On determining the statistical parameters for pollution concentration from a truncated data set | |
| CN111965140B (zh) | 基于特征峰的波长点重组方法 | |
| KR101733026B1 (ko) | 행정구역별 온실가스 분포 데이터 생성 장치 및 방법 | |
| JP3461284B2 (ja) | 赤外線ガス分析計における検量線の作成方法 | |
| Kovalev | Solutions in LIDAR Profiling of the Atmosphere | |
| JP3737257B2 (ja) | ２次元表現によるスペクトルデータの処理方法及び補正方法 | |
| CN105784672A (zh) | 一种基于双树复小波算法的毒品检测仪标准化方法 | |
| CN113111529A (zh) | 融合数值模式和卫星微波云反演数据的红外亮温模拟方法 | |
| Conti et al. | The star formation history of galaxies measured from individual pixels. I. The hubble deep field north | |
| CN108627433B (zh) | 一种环境气溶胶分粒径黑碳质量分布的测量方法及其*** | |
| CN108362697B (zh) | 一种增加分层数目的大气视宁度分层测量方法 | |
| US6507798B1 (en) | Time-frequency dependent damping via Hilbert damping spectrum | |
| Tsai et al. | Line profile reconstruction: validation and comparison of reconstruction methods | |
| Kerekes | An elementary error model for terrestrial laser scanning | |
| JP4233423B2 (ja) | 定量方法及びスペクトル測定装置 | |
| Koner et al. | A comparison of regularization techniques for atmospheric trace gases retrievals |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090815 Year of fee payment: 6 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120815 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120815 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120815 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120815 Year of fee payment: 9 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130815 Year of fee payment: 10 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |