JPH0792407B2 - 電子体温計 - Google Patents
電子体温計Info
- Publication number
- JPH0792407B2 JPH0792407B2 JP61124416A JP12441686A JPH0792407B2 JP H0792407 B2 JPH0792407 B2 JP H0792407B2 JP 61124416 A JP61124416 A JP 61124416A JP 12441686 A JP12441686 A JP 12441686A JP H0792407 B2 JPH0792407 B2 JP H0792407B2
- Authority
- JP
- Japan
- Prior art keywords
- value
- temperature
- time
- memory
- guess
- 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
Landscapes
- Measuring Temperature Or Quantity Of Heat (AREA)
Description
【発明の詳細な説明】 （イ）産業上の利用分野 この発明は、推量式の電子体温計であって、推量値の算
出に応答曲線の応答値を反映させ、個人差に対応し得る
ようにした電子体温計に関する。
出に応答曲線の応答値を反映させ、個人差に対応し得る
ようにした電子体温計に関する。
（ロ）従来の技術 推量式電子体温計は、通常、一定時間温度を測定し、得
られた複数の温度データの変化等に基づいて収束温度、
つまり体温推量値を演算し、この推量値を表示器に順次
更新表示して、測定時間の短縮を図ることが行われてい
る。
られた複数の温度データの変化等に基づいて収束温度、
つまり体温推量値を演算し、この推量値を表示器に順次
更新表示して、測定時間の短縮を図ることが行われてい
る。
この推量値は、予めCPU（セントラルプロセッシングユ
ニット）に設定された推量式により算出されている。
ニット）に設定された推量式により算出されている。
この推量式には、例えば、以下の式が用いられる。
Ｓ（ｔ）＝Ｔ＋Ｋ・dT/dt この式において、Ｔは測定開始からｔ時間経過した時点
の検知温度であり、この検知温度Ｔと測定時間ｔと、検
知温度の上昇率dT/dt及び乗数Ｋが算定要素となってい
る。
の検知温度であり、この検知温度Ｔと測定時間ｔと、検
知温度の上昇率dT/dt及び乗数Ｋが算定要素となってい
る。
ところで、検出温度の時間に対する応答曲線は、被測定
者により異なることが知られている。例えば、収束温度
がほぼ同一であっても、応答曲線は全く異なる場合があ
る。ところが、従来の推量演算の方式は、乗数（Ｋ）を
固定した一定値として算出するものである。つまり、被
測定者により種々異なる応答曲線の応答値を画一的に捉
えたもので、測定時間とその時点の検出温度の他、その
時点の上昇率のみで推量することが行われている。
者により異なることが知られている。例えば、収束温度
がほぼ同一であっても、応答曲線は全く異なる場合があ
る。ところが、従来の推量演算の方式は、乗数（Ｋ）を
固定した一定値として算出するものである。つまり、被
測定者により種々異なる応答曲線の応答値を画一的に捉
えたもので、測定時間とその時点の検出温度の他、その
時点の上昇率のみで推量することが行われている。
このため、応答曲線が緩やかな場合、温度上昇率が小さ
いので、推量値が収束温度、即ち体温より低くなり、推
量精度が適正でないという問題があった。
いので、推量値が収束温度、即ち体温より低くなり、推
量精度が適正でないという問題があった。
そこで、応答曲線の観念を推量値演算のための算定要素
に取り入れ、応答曲線の変化特性を推量値に反映させる
ことを企図した電子体温計が提案されている（特開昭57
−28225号公報）。
に取り入れ、応答曲線の変化特性を推量値に反映させる
ことを企図した電子体温計が提案されている（特開昭57
−28225号公報）。
この体温計は、推量値（Tm）算出のために次の数式を用
いる。
いる。
Tm＝Tq＋Ｋ・dT/dt この式において、Tqは実測温度、Ｋは測温乗数であり、
dT/dtは温度上昇率である。
dT/dtは温度上昇率である。
この推量方式の特徴は、従来、一定の乗数として固定的
・画一的に捉えていた測温乗数Ｋを、予め値の異なる３
種類を設定し、被測定者の応答曲線に応じて３種類の
K1、K2、K3のうちから、何れかの乗数Ｋを選択すること
により、加算値（補正値）を決定するものである。この
方式によれば、乗数Ｋを画一的な固定定数とする方式に
比し、幾分、被測定者によって異なる応答曲線の相違を
推量値に反映させることが可能となる。
・画一的に捉えていた測温乗数Ｋを、予め値の異なる３
種類を設定し、被測定者の応答曲線に応じて３種類の
K1、K2、K3のうちから、何れかの乗数Ｋを選択すること
により、加算値（補正値）を決定するものである。この
方式によれば、乗数Ｋを画一的な固定定数とする方式に
比し、幾分、被測定者によって異なる応答曲線の相違を
推量値に反映させることが可能となる。
（ハ）発明が解決しようとする問題点 上記提案されている電子体温計の推量演算方式は、予め
３つの異なる測温乗数を設定しておき、応答曲線に応じ
て、３つの乗数のうち、応答値に対応する乗数を選択す
るものである。
３つの異なる測温乗数を設定しておき、応答曲線に応じ
て、３つの乗数のうち、応答値に対応する乗数を選択す
るものである。
ところが、この推量方式では、K1、K2、K3の乗数の選択
をdT/dtの絶対値でとる方式であるため、周囲温度によ
る影響を受け易い不利がある。また、固定定数K1、K2、
K3では、連続推量を行う場合に時間の経過と共に応答曲
線が変化し、適正値でなくなる問題がある。更に、３種
類の固定定数では、個人差のバラツキに充分対応し得な
い等の種々の不利があった。
をdT/dtの絶対値でとる方式であるため、周囲温度によ
る影響を受け易い不利がある。また、固定定数K1、K2、
K3では、連続推量を行う場合に時間の経過と共に応答曲
線が変化し、適正値でなくなる問題がある。更に、３種
類の固定定数では、個人差のバラツキに充分対応し得な
い等の種々の不利があった。
この発明は、従来のものが持つ、以上のような問題点を
解消させ、応答曲線の変化特性を推量値に反映させ、個
人差に対応した測定精度の高い推量値が得られる電子体
温計を提供することを目的とする。
解消させ、応答曲線の変化特性を推量値に反映させ、個
人差に対応した測定精度の高い推量値が得られる電子体
温計を提供することを目的とする。
（ニ）問題点を解決するための手段及び作用 この目的を達成させるために、この発明の電子体温計
は、次のような構成としている。
は、次のような構成としている。
電子体温計は、温度を測定する温度測定手段と、この温
度測定手段により検出された実測値を複数のサンプリン
グ分記憶するメモリと、時間についての２次式で近似を
行う補正関数と温度の時間微分との積に体温の実測値を
加えることにより体温推量値を求める推量値演算手段
と、前記メモリ内の最新の実測値または前記推量値演算
手段が求めた推量値を表示する表示手段と、前記メモリ
に記憶されている最新の実測値を含む複数回分の実測値
に基づいて前記補正関数の係数を決定する要素決定手段
と、前記メモリに記憶されている最新の実測値を含む複
数回分の実測値に基づいて前記時間微分値を求める特性
抽出手段とを有し、前記推量値演算手段は、前記要素決
定手段が求めた係数と前記特性抽出手段が求めた温度の
時間微分値と前記メモリの実測値とから前記推量値を求
めるようにしている。
度測定手段により検出された実測値を複数のサンプリン
グ分記憶するメモリと、時間についての２次式で近似を
行う補正関数と温度の時間微分との積に体温の実測値を
加えることにより体温推量値を求める推量値演算手段
と、前記メモリ内の最新の実測値または前記推量値演算
手段が求めた推量値を表示する表示手段と、前記メモリ
に記憶されている最新の実測値を含む複数回分の実測値
に基づいて前記補正関数の係数を決定する要素決定手段
と、前記メモリに記憶されている最新の実測値を含む複
数回分の実測値に基づいて前記時間微分値を求める特性
抽出手段とを有し、前記推量値演算手段は、前記要素決
定手段が求めた係数と前記特性抽出手段が求めた温度の
時間微分値と前記メモリの実測値とから前記推量値を求
めるようにしている。
このような構成を有する電子体温計では、温度測定手段
により実測値を検出し、この実測値を順次メモリに記憶
する一方、応答曲線を分析してパラメータＰ、つまり検
出温度の時間に対する変化特性を、変化率の比や変化率
の変化の比等として抽出する。そして、この抽出値（パ
ラメータＰ）に基づいて、従来は固定の乗数（ａ、ｂ、
ｃ）とされていた算定要素を個々に決定し、更にこの算
定要素により補正関数値ｈ（ｔ）を算出して、初めて推
量値Ｓ（ｔ）を算出する。
により実測値を検出し、この実測値を順次メモリに記憶
する一方、応答曲線を分析してパラメータＰ、つまり検
出温度の時間に対する変化特性を、変化率の比や変化率
の変化の比等として抽出する。そして、この抽出値（パ
ラメータＰ）に基づいて、従来は固定の乗数（ａ、ｂ、
ｃ）とされていた算定要素を個々に決定し、更にこの算
定要素により補正関数値ｈ（ｔ）を算出して、初めて推
量値Ｓ（ｔ）を算出する。
かくて、推量値は応答曲線を基盤として演算されるか
ら、被測定者によってそれぞれ異なる応答曲線が、算出
される推量値に直接反映される結果となり、極めて精確
な推量値が算出できる。
ら、被測定者によってそれぞれ異なる応答曲線が、算出
される推量値に直接反映される結果となり、極めて精確
な推量値が算出できる。
（ホ）実施例 第２図は、この発明に係る電子体温計の具体的な回路構
成例を示すブロック図である。
成例を示すブロック図である。
計器の先細先端部にサーミスタ等の温度センサ１が臨出
させてあり、この温度センサ１は接触する舌下或いは脇
下の温度を検出し、温度に応じた電気信号を取り出す。
そして、随時検出された温度情報（アナログ量）は、A/
D変換器２によりCPU3が処理し易いデジタル値に変換さ
れ、CPU3に取り込まれる。
させてあり、この温度センサ１は接触する舌下或いは脇
下の温度を検出し、温度に応じた電気信号を取り出す。
そして、随時検出された温度情報（アナログ量）は、A/
D変換器２によりCPU3が処理し易いデジタル値に変換さ
れ、CPU3に取り込まれる。
CPU3は、電源スイッチ４がON動作された時、温度センサ
１を作動させ、随時送られる温度情報をメモリ５に記憶
させると共に、表示器６に表示する。またCPU3には、こ
の温度情報に基づき推量値を演算するための温度上昇曲
線（応答曲線）を解析してパラメータＰを算出するパラ
メータＰ算出機能と、このパラメータＰに基づいて推量
値の演算に使用される算定要素（乗数）を決定する算定
要素決定機能と、この算定要素に基づき、応答曲線の補
正関数値ｈ（ｔ）を求める補正関数算出機能と、一定の
数式に前記算定要素を代入し、推量値を算出する推量値
演算機能とを備えている。
１を作動させ、随時送られる温度情報をメモリ５に記憶
させると共に、表示器６に表示する。またCPU3には、こ
の温度情報に基づき推量値を演算するための温度上昇曲
線（応答曲線）を解析してパラメータＰを算出するパラ
メータＰ算出機能と、このパラメータＰに基づいて推量
値の演算に使用される算定要素（乗数）を決定する算定
要素決定機能と、この算定要素に基づき、応答曲線の補
正関数値ｈ（ｔ）を求める補正関数算出機能と、一定の
数式に前記算定要素を代入し、推量値を算出する推量値
演算機能とを備えている。
前記パラメータＰの算出及びパラメータＰにより決定さ
れる算定要素の算出、及び算定要素に基づき演算される
補正関数の算出、更に推量値算出式は、以下の数式理論
によって説明される。
れる算定要素の算出、及び算定要素に基づき演算される
補正関数の算出、更に推量値算出式は、以下の数式理論
によって説明される。
通常、体温計で測定した場合における測定時間と温度上
昇曲線とは、次式によって表される。
昇曲線とは、次式によって表される。
％Temp この式において、τaは体温計感温部の熱応答時の定数
であり、t00、t11は体表と体温計間の熱伝導、熱容量等
に関するパラメータであり、τa、tcoは人体深部と体表
間の熱伝導に関するパラメータである。
であり、t00、t11は体表と体温計間の熱伝導、熱容量等
に関するパラメータであり、τa、tcoは人体深部と体表
間の熱伝導に関するパラメータである。
そこで、例えば第３図に示すような応答例に対して、前
記（１）式の各パラメータを求めると、 τa＝1.83,t00＝−5.77,t11＝−3.14, τc＝464.7,tco＝−1147.9, となり、第４図で示すように、０≦ｔ≦10秒ではτaが
支配的であり、10≦ｔ≦240秒程度ではt00、t11が支配
的であり、240≦ｔではτc、tcoが支配的となっている
ことが分かる。
記（１）式の各パラメータを求めると、 τa＝1.83,t00＝−5.77,t11＝−3.14, τc＝464.7,tco＝−1147.9, となり、第４図で示すように、０≦ｔ≦10秒ではτaが
支配的であり、10≦ｔ≦240秒程度ではt00、t11が支配
的であり、240≦ｔではτc、tcoが支配的となっている
ことが分かる。
ここで、従来技術では、推量値Ｓ（ｔ）を求める場合、
Ｓ（ｔ）＝Tq＋K1・ｄTq/dtにより算出されていた。こ
の式において、Tqは実測温度であり、K1は測温定数であ
る。
Ｓ（ｔ）＝Tq＋K1・ｄTq/dtにより算出されていた。こ
の式において、Tqは実測温度であり、K1は測温定数であ
る。
そこで、（１）式に基づいて推量値Ｓ（ｔ）をＳ（ｔ）
＝Ｔ＋ｈ（ｔ）・dT/dtの形に変形させると、 だだし つまり、推量値Ｓ（ｔ）は、 Ｓ（ｔ）＝Ｔ＋ｈ（ｔ）・dT/dt……（３） と表される。
＝Ｔ＋ｈ（ｔ）・dT/dtの形に変形させると、 だだし つまり、推量値Ｓ（ｔ）は、 Ｓ（ｔ）＝Ｔ＋ｈ（ｔ）・dT/dt……（３） と表される。
前記（２）式にて、ｔ＞10にて ただし、 となる。
この（４）式を多項式近似すると、K00＝100より、２次
式にて近似すると、 となる。この式において、a0乃至a2は定数である。t00
は、体表と体温計の熱伝導・熱容量等に関するパラメー
タであり、個人差により特徴的に決定されるものであ
る。
式にて近似すると、 となる。この式において、a0乃至a2は定数である。t00
は、体表と体温計の熱伝導・熱容量等に関するパラメー
タであり、個人差により特徴的に決定されるものであ
る。
第５図は、第３図の応答例に対して（２）式の補正関数
を算出した例である。ここで、ｔ≦60秒として２次式に
よる近似を行うと、60≦ｔ≦200では、 h1（ｔ）＝−0.0054(t-200)2＋255……（６）また、200
≦ｔ≦600では、 h2（ｔ）＝−0.0014(t-200)2＋255……（７）となる。
但しh1（ｔ）、h2（ｔ）＜０の時は、h1（ｔ）＝０、h2
（ｔ）＝０とする。
を算出した例である。ここで、ｔ≦60秒として２次式に
よる近似を行うと、60≦ｔ≦200では、 h1（ｔ）＝−0.0054(t-200)2＋255……（６）また、200
≦ｔ≦600では、 h2（ｔ）＝−0.0014(t-200)2＋255……（７）となる。
但しh1（ｔ）、h2（ｔ）＜０の時は、h1（ｔ）＝０、h2
（ｔ）＝０とする。
尚、近似は、２次以上の多項式でも、テーブル関数を用
いても良い。
いても良い。
更に、前記（６）式、（７）式の補正関数近似値を用い
て、（３）式により推量値Ｓ（ｔ）を算出すると、下表
のようになる。
て、（３）式により推量値Ｓ（ｔ）を算出すると、下表
のようになる。
第６図は、この応答曲線の推量値算出例を示している。
上表に示すように、h1（ｔ）、h2（ｔ）を使用して推量
値を算出すると、600秒の値に対して±0.01℃の誤差で
算出が可能であり、実用上は２次式を用いての近似でも
問題のないことが分かる。また、簡単化のためh1（ｔ）
のみを使用しても、240秒程度までは精度上問題ではな
い。上表から明らかなように、h1（ｔ）のみを使用の場
合、240秒以後420秒までの間、僅かずつ温度が低下傾向
を示し、420秒後に再び上昇して600秒の時点ではh
1（ｔ）、h2（ｔ）を使用する場合と同じ37.27℃まで上
昇する（第６図参照）。従って、240秒までの間であれ
ば、h1（ｔ）のみを使用する簡便方式で充分な精度が得
られる。本来、短時間にて飽和値を知る推量式電子体温
計では、長時間測定を行う意味がなく、長時間測定を行
う場合は、実測測定の方がより信頼性があるため、実用
上は好ましいと言える。
上表に示すように、h1（ｔ）、h2（ｔ）を使用して推量
値を算出すると、600秒の値に対して±0.01℃の誤差で
算出が可能であり、実用上は２次式を用いての近似でも
問題のないことが分かる。また、簡単化のためh1（ｔ）
のみを使用しても、240秒程度までは精度上問題ではな
い。上表から明らかなように、h1（ｔ）のみを使用の場
合、240秒以後420秒までの間、僅かずつ温度が低下傾向
を示し、420秒後に再び上昇して600秒の時点ではh
1（ｔ）、h2（ｔ）を使用する場合と同じ37.27℃まで上
昇する（第６図参照）。従って、240秒までの間であれ
ば、h1（ｔ）のみを使用する簡便方式で充分な精度が得
られる。本来、短時間にて飽和値を知る推量式電子体温
計では、長時間測定を行う意味がなく、長時間測定を行
う場合は、実測測定の方がより信頼性があるため、実用
上は好ましいと言える。
次に応答曲線を解析し、応答値（パラメータＰ）を算出
する理論を説明する。（５）式の補正関数の近似は、t
00のパラメータにより決定される。
する理論を説明する。（５）式の補正関数の近似は、t
00のパラメータにより決定される。
従って、次に（１）式からt00のパラメータを抽出する
手段を考えると、（１）式よりｔ＞10にて、 ただし、 となる。更に、この（８）式をｔで微分すると、 となる。ここで、ｔ＝ta、tb、tc（ta＜tb＜tc）の３点
の近似値を考えると、τc＞400より、 となり、前記（９）式は、 と近似できる。
手段を考えると、（１）式よりｔ＞10にて、 ただし、 となる。更に、この（８）式をｔで微分すると、 となる。ここで、ｔ＝ta、tb、tc（ta＜tb＜tc）の３点
の近似値を考えると、τc＞400より、 となり、前記（９）式は、 と近似できる。
かくて、（Ta−Tb）÷（Tb−Tc）より、 が得られる。ここで、ta、tb、tcをΔｔの等間隔でとる
と、 ta……（15） tb＝ta＋Δｔ……（16） tc＝ta＋２Δｔ……（17） となり、これらを（14）式に代入すると、 が得られる。この（18）式は近似的に、 となり、これにより、 となる。ここにおいて、t00が３点の温度変化率の差の
比率と関連してくることがわかる。
と、 ta……（15） tb＝ta＋Δｔ……（16） tc＝ta＋２Δｔ……（17） となり、これらを（14）式に代入すると、 が得られる。この（18）式は近似的に、 となり、これにより、 となる。ここにおいて、t00が３点の温度変化率の差の
比率と関連してくることがわかる。
ところで、この式は実際のデータに於いては、取り得る
値が限られているため、次式のように近似できる。
値が限られているため、次式のように近似できる。
ここにおいて、Tab／Tbcは、等間隔３点と各点の少なく
とも１点の温度データがあれば算出できることとなる。
例えば、第７図に示す３点ａ、ｂ、ｃと各点の前後点ａ
−１、ａ＋１、ｂ−１、ｂ＋１、ｃ−１、ｃ＋１より、 となり、ここで、 a-1−a-1＝b-1−b-1＝c-1−c-1＝Δｔ′と取ると、 となる。更に、a+1＝b-1、b+1＝c-1になるようにΔｔ′
を取ると、４点のデータT0、T1、T2、T3となり、 更に、簡略化のために、T2＝T1＝T12と取ると、３点デ
ータにて、 ここにおいて、多数の実際データを解析した結果、（2
7）式のパラメータと補正関数式には、下式のような相
関が得られた。
とも１点の温度データがあれば算出できることとなる。
例えば、第７図に示す３点ａ、ｂ、ｃと各点の前後点ａ
−１、ａ＋１、ｂ−１、ｂ＋１、ｃ−１、ｃ＋１より、 となり、ここで、 a-1−a-1＝b-1−b-1＝c-1−c-1＝Δｔ′と取ると、 となる。更に、a+1＝b-1、b+1＝c-1になるようにΔｔ′
を取ると、４点のデータT0、T1、T2、T3となり、 更に、簡略化のために、T2＝T1＝T12と取ると、３点デ
ータにて、 ここにおいて、多数の実際データを解析した結果、（2
7）式のパラメータと補正関数式には、下式のような相
関が得られた。
ｔ＝10、20、30秒の３点の温度データをT10、T20、T30
とすると、パラメータＰは次式となる。
とすると、パラメータＰは次式となる。
更に、補正関数ｈ（ｔ）は下式で表される。
ｈ（ｔ）＝ａ(t+b)2＋ｃ……（29） また、算定要素たる乗数ａ、ｂ、ｃは次式で求められ
る。
る。
ａ＝−0.00282＋0.0024・Ｐ……（30） ｂ＝−320＋110・Ｐ……（31） ｃ＝192−120・Ｐ……（32） そして、推量値Ｓ（ｔ）は次式となる。
以上のような数式理論に基づき、上記のような推量値演
算式が求められた。
算式が求められた。
かくして、推量値を演算するために行う温度カーブ（応
答曲線）の解析は、感温部の応答が飽和してから行うた
め、周囲温度の影響を受け難く、また、（28）式のパラ
メータは、個人にて特有のものであり、推量値算出にこ
のパラメータを導入することで、個人差による対応が可
能となる。
答曲線）の解析は、感温部の応答が飽和してから行うた
め、周囲温度の影響を受け難く、また、（28）式のパラ
メータは、個人にて特有のものであり、推量値算出にこ
のパラメータを導入することで、個人差による対応が可
能となる。
尚、このCPU3には、ブザー７及び電源８がそれぞれ連繋
させてあり、ブザー７は推量値Ｓ（ｔ）がほぼ一定値に
なると作動し、その旨を報知するように設定されてい
る。
させてあり、ブザー７は推量値Ｓ（ｔ）がほぼ一定値に
なると作動し、その旨を報知するように設定されてい
る。
もっとも、第２図に示すハード構成は、従来より周知の
一般的な電子体温計と変わらない。この実施例電子体温
計は、CPU3の保有する機能構成に特徴がある。
一般的な電子体温計と変わらない。この実施例電子体温
計は、CPU3の保有する機能構成に特徴がある。
それゆえ、次に第１図に示すフローにより、この実施例
電子体温計のソフト構成及び動作について説明する。
電子体温計のソフト構成及び動作について説明する。
電源スイッチ４がONされると、スタートフラグ、ステイ
トフラグ及びタイマがそれぞれイニシャライズされ、計
器の初期化が行われる〔ステップ（以下『ST』という）
１〕。
トフラグ及びタイマがそれぞれイニシャライズされ、計
器の初期化が行われる〔ステップ（以下『ST』という）
１〕。
スタートフラグとは、電源スイッチ４がONされた後、セ
ンサ温度が上昇し、検温開始状態を検知したか否かを判
断するためのフラグである。
ンサ温度が上昇し、検温開始状態を検知したか否かを判
断するためのフラグである。
また、ステイトフラグとは、実測が開始された後、パラ
メータＰの算出処理を実行させるか、推量値の演算処理
を実行させるかを選択するためのフラグである。
メータＰの算出処理を実行させるか、推量値の演算処理
を実行させるかを選択するためのフラグである。
更に、タイマとは、電源スイッチ４がONされた測定開始
時点から終了までを計時するために使用される。
時点から終了までを計時するために使用される。
電源スイッチ４がONされた後、次のST2では、サンプリ
ングタイムか否かを判定している。
ングタイムか否かを判定している。
実施例では、１秒毎に実測温度を計測する。従って、こ
のサンプリングタイム毎にセンサ１が現在温度情報をCP
U3に送り、順次上昇する実測値がメモリ５に更新記憶さ
れる一方、このサンプリングタイム毎に以下の動作が行
われ、ST2へ戻ることになる。
のサンプリングタイム毎にセンサ１が現在温度情報をCP
U3に送り、順次上昇する実測値がメモリ５に更新記憶さ
れる一方、このサンプリングタイム毎に以下の動作が行
われ、ST2へ戻ることになる。
今、サンプリングタイムが到来したとすると、このST2
の判定が“YES"となり、測定開始から時間を計時する時
間ｔがサンプリングタイム毎に１歩進される（ST3）。
そして、現在温度Ttが検出される（ST4）。
の判定が“YES"となり、測定開始から時間を計時する時
間ｔがサンプリングタイム毎に１歩進される（ST3）。
そして、現在温度Ttが検出される（ST4）。
ST5では、現在のタイムカウントが奇数であるか否かを
判定している。奇数秒でなければ、このST5の判定が“N
O"となってST7へ進み、逆に奇数秒であれば判定が“YE
S"となり、次のST6でこの実測値がメモリ５に更新記録
される（ST6）。つまり、実施例では温度データの記録
は、奇数秒のみ実行するように設定されている。
判定している。奇数秒でなければ、このST5の判定が“N
O"となってST7へ進み、逆に奇数秒であれば判定が“YE
S"となり、次のST6でこの実測値がメモリ５に更新記録
される（ST6）。つまり、実施例では温度データの記録
は、奇数秒のみ実行するように設定されている。
また、実施例では、過去４度の実測値（T0、T1、T2、
T3）が常時、レジスタに記憶保持されるようになってい
る（ST6）。従って、ST7において、今回の実測値Ttが、
メモリ５に記憶されるT3よりも高いか否か判定され、高
い場合にT3をT2とし、今回の実測値TtをT3として記憶す
る（ST8）。この処理が奇数秒ごとに実行される。
T3）が常時、レジスタに記憶保持されるようになってい
る（ST6）。従って、ST7において、今回の実測値Ttが、
メモリ５に記憶されるT3よりも高いか否か判定され、高
い場合にT3をT2とし、今回の実測値TtをT3として記憶す
る（ST8）。この処理が奇数秒ごとに実行される。
ST9では、スタートフラグが何の状態にあるかを判定し
ている。今、スタートフラグは０、つまり電源スイッチ
４がONされただけであり、未だ検出温度は室温状態にあ
り、検温開始が検知されていない状態とすると、次のST
10で今回の実測値T3と前回の実測値T2との差値が、所定
値（例えば一定値の１℃或いは１℃に相当する温度上昇
率）より高いか否かが判断される。勿論、このケースで
は検温開始状態にないから、このST10の判定が“NO"と
なり、ST2へ戻る。
ている。今、スタートフラグは０、つまり電源スイッチ
４がONされただけであり、未だ検出温度は室温状態にあ
り、検温開始が検知されていない状態とすると、次のST
10で今回の実測値T3と前回の実測値T2との差値が、所定
値（例えば一定値の１℃或いは１℃に相当する温度上昇
率）より高いか否かが判断される。勿論、このケースで
は検温開始状態にないから、このST10の判定が“NO"と
なり、ST2へ戻る。
今、検温が開始され実測値Ttが上昇し、この上昇値のタ
イマカウントが奇数秒であるとすると、ST5及びST7の判
定がいずれも“YES"となり、この実測値がピーク値T3と
してメモリ５に記憶される（ST7）。ここでは、スター
トフラグは未だ０の状態にある。従って、ST9を介してS
T10へ進み、このST10の判定が“YES"となって、ST11で
スタートフラグ０が１にセットされる。つまり、ここに
おいて検温開始と判定され、検温の開始時点から始まる
タイマが改めて０にセットされ（ST11）、以後タイマの
計時と共に同様な温度測定が繰返し実行される。
イマカウントが奇数秒であるとすると、ST5及びST7の判
定がいずれも“YES"となり、この実測値がピーク値T3と
してメモリ５に記憶される（ST7）。ここでは、スター
トフラグは未だ０の状態にある。従って、ST9を介してS
T10へ進み、このST10の判定が“YES"となって、ST11で
スタートフラグ０が１にセットされる。つまり、ここに
おいて検温開始と判定され、検温の開始時点から始まる
タイマが改めて０にセットされ（ST11）、以後タイマの
計時と共に同様な温度測定が繰返し実行される。
そして、以後のサンプリングタイムでは、スタートフラ
グが１となっているから、ST9を介してST13へ進む。ST1
3では、ステイトフラグが何の状態にあるかを判定して
いる。今、ステイトフラグは０である。つまり、実測開
始後40秒が経過していないため、推量値を演算するため
の温度データが検出されていない状態である。従って、
ST14へ移行する。
グが１となっているから、ST9を介してST13へ進む。ST1
3では、ステイトフラグが何の状態にあるかを判定して
いる。今、ステイトフラグは０である。つまり、実測開
始後40秒が経過していないため、推量値を演算するため
の温度データが検出されていない状態である。従って、
ST14へ移行する。
ST14乃至ST16は、第８図で示すように応答曲線を分析
し、パラメータＰを算出するための処理動作を実行す
る。ST14では時間がta、つまり、例えば20秒経過したか
否かを判定し、経過しておれば、このST14の判定が“YE
S"とり、この時点の検知温度T3（第８図では36.31℃）
をTaとする（ST17）。その後、このTa（36.31℃）を表
示器６に表示し（ST12）、ST2に戻る。そして、以後同
様に時間がtbつまり30秒経過（ta時点から10秒経過）し
た時、ST15の判定が“YES"となり、この時点の検知温度
T3（36.49℃）をTbとする（ST18）。更に、時間がtcつ
まり40秒経過したとすると、ST16の判定が“YES"とな
り、この時点の検知温度T3（36.60℃）がTcとされる（S
T19）。
し、パラメータＰを算出するための処理動作を実行す
る。ST14では時間がta、つまり、例えば20秒経過したか
否かを判定し、経過しておれば、このST14の判定が“YE
S"とり、この時点の検知温度T3（第８図では36.31℃）
をTaとする（ST17）。その後、このTa（36.31℃）を表
示器６に表示し（ST12）、ST2に戻る。そして、以後同
様に時間がtbつまり30秒経過（ta時点から10秒経過）し
た時、ST15の判定が“YES"となり、この時点の検知温度
T3（36.49℃）をTbとする（ST18）。更に、時間がtcつ
まり40秒経過したとすると、ST16の判定が“YES"とな
り、この時点の検知温度T3（36.60℃）がTcとされる（S
T19）。
ここにおいて、パラメータＰ算出のための温度情報が得
られたこととなる。従って、ST20において、前述のパラ
メータＰの算出式にこのデータ（Ta、Tb、Tc）を代入す
ることで、パラメータＰが0.61と算出される（ST20）。
つまり、温度曲線を多点比率により解析することとな
り、個人差の応答曲線に対応し得、且つ多点比率によっ
て周囲の温度による悪影響も防止できる。
られたこととなる。従って、ST20において、前述のパラ
メータＰの算出式にこのデータ（Ta、Tb、Tc）を代入す
ることで、パラメータＰが0.61と算出される（ST20）。
つまり、温度曲線を多点比率により解析することとな
り、個人差の応答曲線に対応し得、且つ多点比率によっ
て周囲の温度による悪影響も防止できる。
その後、ステイトフラグが１にセットされ（ST21）、ピ
ーク値を表示器６に表示して（ST12）、ST2へ戻る。
ーク値を表示器６に表示して（ST12）、ST2へ戻る。
次のサンプリングタイムでは、実測開始から40秒が経過
しており、推量可能となっている。つまり、ステイトフ
ラグが１になっている。
しており、推量可能となっている。つまり、ステイトフ
ラグが１になっている。
従って、ST13の判定が“YES"となり、ST22へ進む。ST22
では、現在のサンプリングタイムが奇数秒であるか否か
を判定している。
では、現在のサンプリングタイムが奇数秒であるか否か
を判定している。
この実施例ではST5と同様、奇数秒の実測値（ピーク
値）のみを対象として推量値の算出を行うように設定し
ている。
値）のみを対象として推量値の算出を行うように設定し
ている。
今、サンプリングタイムが奇数秒であるとすると、この
ST22の判定が“YES"となり、ST23へ進み、ここにおいて
先ず、温度上昇率dT/dtが、T3−T0/6により算出され
る。
ST22の判定が“YES"となり、ST23へ進み、ここにおいて
先ず、温度上昇率dT/dtが、T3−T0/6により算出され
る。
その後、推量値演算のための算定要素である乗数ａ、
ｂ、ｃを、パラメータＰより決定する（ST24）。つま
り、前述の（30）乃至（32）式により決定される。そし
て、乗数ａ、ｂ、ｃが求まると、次に補正関数値ｈ
（ｔ）を求める（ST25）。補正関数ｈ（ｔ）は前述の
（29）式にａ、ｂ、ｃを代入することで算出される。つ
まり、実施例では、補正係数に多次の時間関数を使用す
るため、連続推量によった場合であっても、適正値が得
られる。
ｂ、ｃを、パラメータＰより決定する（ST24）。つま
り、前述の（30）乃至（32）式により決定される。そし
て、乗数ａ、ｂ、ｃが求まると、次に補正関数値ｈ
（ｔ）を求める（ST25）。補正関数ｈ（ｔ）は前述の
（29）式にａ、ｂ、ｃを代入することで算出される。つ
まり、実施例では、補正係数に多次の時間関数を使用す
るため、連続推量によった場合であっても、適正値が得
られる。
ここにおいて、推量値演算のための全ての算定要素が決
定された。かくして、推量値Ｓ（ｔ）が前述の（33）式
に基づいて算出され（ST26）、この推量値が表示器６に
表示される（ST27）。
定された。かくして、推量値Ｓ（ｔ）が前述の（33）式
に基づいて算出され（ST26）、この推量値が表示器６に
表示される（ST27）。
尚、上記の推量演算は奇数秒のサンプリングタイム毎に
実行される。そして、このフローチャートには示してい
ないが、ここにおいて表示される推量値が精度高く、信
用できる旨をブザー７で報知する。その後、電源スイッ
チ４のOFFで測定が終了する。
実行される。そして、このフローチャートには示してい
ないが、ここにおいて表示される推量値が精度高く、信
用できる旨をブザー７で報知する。その後、電源スイッ
チ４のOFFで測定が終了する。
（ヘ）発明の効果 この発明では、以上のように、応答曲線を分析してパラ
メータＰを算出し、このパラメータＰから推量値を演算
するための算定要素を決定し、この算定要素から補正関
数値を求め、推量値を算出することとした。
メータＰを算出し、このパラメータＰから推量値を演算
するための算定要素を決定し、この算定要素から補正関
数値を求め、推量値を算出することとした。
この発明によれば、推量値の算定要素を検知温度の時間
に対する変化特性に対応して定める。従って、変化特性
が異なる場合には算定要素が異なることとなるので、変
化特性が推量値に反映されることとなり、個人差に対応
した一層精確な推量演算が可能となる。
に対する変化特性に対応して定める。従って、変化特性
が異なる場合には算定要素が異なることとなるので、変
化特性が推量値に反映されることとなり、個人差に対応
した一層精確な推量演算が可能となる。
また、この発明では、パラメータＰにより乗数ａ、ｂ、
ｃを求め、これから時間の補正関数ｈ（ｔ）を求めるこ
ととしたから、連続推量によっても適正値を得ることが
可能となる。更に、この発明では応答曲線の分析に多点
比率を用いることとしたから、周囲温度の影響を受ける
ことが少なく、より精度の高い測定が実行できる等、発
明目的を達成した優れた効果を有する。
ｃを求め、これから時間の補正関数ｈ（ｔ）を求めるこ
ととしたから、連続推量によっても適正値を得ることが
可能となる。更に、この発明では応答曲線の分析に多点
比率を用いることとしたから、周囲温度の影響を受ける
ことが少なく、より精度の高い測定が実行できる等、発
明目的を達成した優れた効果を有する。
第１図は、実施例電子体温計の処理動作を示すフローチ
ャート、第２図は、この電子体温計の回路構成例を示す
ブロック図、第３図は、時間に対する検知温度の応答曲
線図、第４図は、第３図に於ける100％比率による応答
曲線図、第５図は、第３図の応答曲線に対する補正関数
値を示す説明図、第６図は、第３図の応答曲線に対する
推量値算出例を示す説明図、第７図は、パラメータＰを
求める算出例を示す説明図、第８図は、パラメータＰの
具体的な算出例を示す説明図である。 1:温度センサ、3:CPU、5:メモリ、6:表示器。
ャート、第２図は、この電子体温計の回路構成例を示す
ブロック図、第３図は、時間に対する検知温度の応答曲
線図、第４図は、第３図に於ける100％比率による応答
曲線図、第５図は、第３図の応答曲線に対する補正関数
値を示す説明図、第６図は、第３図の応答曲線に対する
推量値算出例を示す説明図、第７図は、パラメータＰを
求める算出例を示す説明図、第８図は、パラメータＰの
具体的な算出例を示す説明図である。 1:温度センサ、3:CPU、5:メモリ、6:表示器。
フロントページの続き (56)参考文献 特開 昭60−157031（ＪＰ，Ａ) 特開 昭60−147621（ＪＰ，Ａ) 特開 昭56−46440（ＪＰ，Ａ)
Claims (1)
- 【請求項１】温度を測定する温度測定手段と、この温度
測定手段により検出された実測値を複数のサンプリング
分記憶するメモリと、時間についての２次式で近似を行
う補正関数と温度の時間微分との積に体温の実測値を加
えることにより体温推量値を求める推量値演算手段と、
前記メモリ内の最新の実測値または前記推量値演算手段
が求めた推量値を表示する表示手段と、前記メモリに記
憶されている最新の実測値を含む複数回分の実測値に基
づいて前記補正関数の係数を決定する要素決定手段と、
前記メモリに記憶されている最新の実測値を含む複数回
分の実測値に基づいて前記時間微分値を求める特性抽出
手段とを有し、前記推量値演算手段は、前記要素決定手
段が求めた係数と前記特性抽出手段が求めた温度の時間
微分値と前記メモリの実測値とから前記推量値を求める
こと特徴とする電子体温計。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP61124416A JPH0792407B2 (ja) | 1986-05-28 | 1986-05-28 | 電子体温計 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP61124416A JPH0792407B2 (ja) | 1986-05-28 | 1986-05-28 | 電子体温計 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPS62280629A JPS62280629A (ja) | 1987-12-05 |
| JPH0792407B2 true JPH0792407B2 (ja) | 1995-10-09 |
Family
ID=14884935
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP61124416A Expired - Lifetime JPH0792407B2 (ja) | 1986-05-28 | 1986-05-28 | 電子体温計 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0792407B2 (ja) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006308540A (ja) * | 2005-03-29 | 2006-11-09 | Citizen Watch Co Ltd | 電子体温計 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4873274A (en) * | 1987-09-24 | 1989-10-10 | Morton Thiokol, Inc. | In-mold coating powders with two initiators or 1,1-di(t-butylperoxy)-3,3,5-trimethylcyclohexane as a single initiator |
Family Cites Families (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH0238895B2 (ja) * | 1979-09-21 | 1990-09-03 | Terumo Corp | Denshitaionkei |
| JPS6093930A (ja) * | 1983-10-28 | 1985-05-25 | Seiko Instr & Electronics Ltd | 予測機能付電子体温計 |
| JPS60114731A (ja) * | 1983-11-25 | 1985-06-21 | Omron Tateisi Electronics Co | 電子体温計 |
| JPS60147621A (ja) * | 1984-01-12 | 1985-08-03 | Seiko Instr & Electronics Ltd | 電子体温計 |
| JPS60157031A (ja) * | 1984-01-27 | 1985-08-17 | Terumo Corp | 電子体温計 |
-
1986
- 1986-05-28 JP JP61124416A patent/JPH0792407B2/ja not_active Expired - Lifetime
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2006308540A (ja) * | 2005-03-29 | 2006-11-09 | Citizen Watch Co Ltd | 電子体温計 |
| US7637657B2 (en) | 2005-03-29 | 2009-12-29 | Citizen Holdings Co., Ltd. | Electronic thermometer |
Also Published As
| Publication number | Publication date |
|---|---|
| JPS62280629A (ja) | 1987-12-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| EP1183508B1 (en) | Predictive temperature measurement system | |
| CN101124464B (zh) | 温度预测系统和方法 | |
| US6439768B1 (en) | Electronic thermometer and temperature prediction method therefor | |
| EP0332700A1 (en) | Electronic thermometer | |
| JPS62165132A (ja) | 電子体温計 | |
| JPS61193037A (ja) | 電子式体温計 | |
| JP2005098982A (ja) | 電子体温計 | |
| JPH0792407B2 (ja) | 電子体温計 | |
| CA1177277A (en) | Electronic thermometer | |
| JPH073368B2 (ja) | 電子体温計 | |
| JPS62192625A (ja) | 電子体温計 | |
| JPS62175627A (ja) | 電子体温計 | |
| JPH047456B2 (ja) | ||
| JPS60114731A (ja) | 電子体温計 | |
| JPS62185133A (ja) | 電子体温計 | |
| JPH0792406B2 (ja) | 電子体温計 | |
| JPH0412409B2 (ja) | ||
| JPS62185132A (ja) | 電子体温計 | |
| JPS62267627A (ja) | 電子体温計 | |
| JPS60187830A (ja) | 電子体温計 | |
| JPS60114730A (ja) | 電子体温計 | |
| JPS62195530A (ja) | 電子体温計 | |
| JPS62161029A (ja) | 電子体温計 | |
| JPS63144223A (ja) | 電子体温計 | |
| JPS62190425A (ja) | 電子体温計 |