(57)【要約】      （修正有）【課題】測定開始時における測定温度と体温との温度差が大きい場合における予測誤差を低減する。【解決手段】予測関数によって熱平衡状態に至った平衡温度を予測するものであり、測定開始時における測定温度と体温との温度差が大きい場合には、予測値を補正することによって高温側に生じる予測誤差を低減する。補正値算出部は、パラメータ算出部が算出するパラメータｂの値と、温度検出手段が検出する温度の変化から求めた時間長とに基づいて判定条件から補正値αを算出する。あるいは、温度検出手段が検出する測定開始時の測定温度に基づいて判定条件から補正値αを算出する。【選択図】図４
2010122164.tif 000002
【特許請求の範囲】
【請求項１】
被測定部位の温度を検出する温度検出手段と、
検出した温度に応じて平衡温度を算出する予測手段とを備え、
前記予測手段は、
２つのパラメータａ，ｂおよび補正値αを有する平衡予測温度を求める予測関数
Ｔu＝α・（Ｃ−ｔ）＋ａ・ｔ＾ｂ
（Ｔu：熱平衡予測温度
ｔ：測定開始からの測定時間
ａ：予測関数の曲線の倍率を定めるパラメータ
ｂ：予測関数の曲線の形状を定めるパラメータ
α：予測関数を直線状に補正する補正値
Ｃ：時間定数
＾：冪乗）
のパラメータａ，ｂおよび補正値αを求めて予測関数を決定する予測関数決定部と、
前記予測関数決定部により決定された予測関数を用いて熱平衡時間時における平衡温度を算出する熱平衡予測温度算出手段とを備え、
前記予測関数決定部は、
前記予測関数の曲線の倍率を定めるパラメータａおよび前記予測関数の曲線の形状を定めるパラメータｂを算出するパラメータ算出部と、
前記予測関数の測定開始から熱平衡時点までを直線状に補正する補正値αを、測定温度の増加特性を判定条件として算出する補正値算出部と、
を備えることを特徴とする電子体温計。
【請求項２】
前記増加特性を前記予測関数の曲線の形状を定めるパラメータｂの値と測定開始時から温度変化が安定するまでの時間長とによって特定し、
前記補正値算出部は、
前記パラメータｂの値と前記時間長とを判定条件とする補正値αを予め定めて備え、
前記パラメータ算出部が算出するパラメータｂの値と、前記温度検出手段が検出する温度の変化から求めた時間長とに基づいて前記判定条件から補正値αを算出することを特徴とする、請求項１に記載の電子体温計。
【請求項３】
前記測定温度の増加特性を測定温度の温度変化幅によって特定し、
前記補正値算出部は、
測定温度の温度変化幅を表す値として測定開始時の測定温度を用い、当該測定開始時の測定温度を判定条件とする補正値αを予め定めて備え、
前記温度検出手段が検出する測定開始時の測定温度に基づいて前記判定条件から補正値αを算出することを特徴とする、請求項１に記載の電子体温計。
【請求項４】
前記パラメータ算出部は、
前記補正値αを零とし、２測定点の検出温度および当該検出温度の測定時間を用いてパラメータｂを決定すると共に仮パラメータａ′を仮決定し、
さらに、前記パラメータｂ及び前記補正値α、並びに、２測定点の何れか一方の測定点の検出温度及び当該検出温度の測定時間を前記予測関数に代入してパラメータａを再決定することを特徴とする請求項１から３にいずれか一つに記載の電子体温計。
【請求項５】
前記パラメータ算出部は、
パラメータｂを、
ｂ＝ln（Ｔ１／Ｔ２）／ln（ｔ１／ｔ２）
（Ｔ１、Ｔ２は２測定点の検出温度、ｔ１、ｔ２は測定時間、lnは対数）
の演算式で算出し、
前記仮決定される仮パラメータａ′を、
lnａ′＝（lnＴ２・lnｔ１−lnＴ１・lnｔ２）／（lnｔ１−lnｔ２）
（Ｔ１、Ｔ２は２測定点の検出温度、ｔ１、ｔ２は測定時間、lnは対数）
の演算式で算出し、
前記再決定されるパラメータａを
予測関数Ｔu＝α・（Ｃ−ｔ）＋ａ・ｔ＾ｂ
の演算式で算出することを特徴とする請求項４に記載の電子体温計。
【請求項６】
前記予測手段は、
前記温度検出手段により検出された温度に基づいて温度勾配を算出する温度勾配算出部を備え、
前記温度勾配算出部は、検出温度の温度勾配が所定範囲内となった時点の検出温度および当該検出温度の測定時間を、前記２測定点の検出温度と当該検出温度の測定時間の内の最初の第１点目とすることを特徴とする請求項１から５のいずれか一つに記載の電子体温計。
【請求項７】
前記予測手段は計時部を備え、
前記計時部は、前記第１点目の温度検出からの時間を計時して、第１点目の測定時間から所定時間後の時間を第２点目の測定時間とすることを特徴とする請求項６に記載の電子体温計。
【請求項８】
前記熱平衡予測温度算出手段は、前記パラメータａ，パラメータｂ、および補正値αで定まる予測関数の測定時間ｔに熱平衡予測時間を代入して熱平衡予測温度を算出することを特徴とする請求項１から７のいずれか一つに記載の電子体温計。
【請求項９】
前記時定数Ｃは、熱平衡予測時間であることを特徴とする、請求項１から８のいずれか一つに記載記載の電子体温計。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、検出温度から平衡温度を予測する電子体温計に関する。
【背景技術】
【０００２】
一般に、電子体温計は、温度センサで体温を検出し、この温度センサの検出信号を信号処理し、その実測値を表示している。測定開始直後には温度センサの温度と体温との間に温度差があるため、温度センサが被測定部位と熱平衡状態となるまでに時間を要する。
【０００３】
そこで、電子体温計では、実測値の変化から熱平衡状態での平衡温度を予測し、この予測値を順次更新して表示することによって、測定時間を短縮している。
【０００４】
平衡温度の予測は、典型的には測定温度とその変化率を監視し、この測定温度と変化率を変数とする予測関数を演算することで行われる。この予測関数を用いて平衡温度を予測する電子体温計として種々のものが提案されている。
【０００５】
例えば、予測関数による平衡温度の予測では、温度予測として用いられた予測関数が適切でない場合には、予測精度が低下するという問題がある。この問題に対して、平衡温度を予測する予測関数を複数用意しておき、所定間隔で平衡温度を予測し、今回の平衡温度予測値と前回の平衡温度予測値との差が所定範囲外の場合には、新たな予測関数を選択し、所定範囲内の場合には平衡温度予測値を表示し、予測演算を続行する電子体温計が提案されている（特許文献１参照）。
【０００６】
また、個人的ばらつきによる測定精度へ影響を低減させる平衡温度の予測として、センサから出力される体温をサンプリングし、各検出出力の時間微分の対数値に基づいて直線式ＴL＝Ａ−τ´ｔ（ＴLは各検出出力の時間微分の対数値、ｔは時間、Ａ及びτ´はパラメータ部）のパラメータ部を回帰法で算出し、算出された値から熱平衡後の体温を予測演算するものが提案されている（特許文献２参照）。
【０００７】
また、センサの初期温度や熱時定数の影響を受けずに平衡温度を予測するものとして、センサと被測定体とが熱平衡に達する以前における３点の値Ｔ１，Ｔ２，Ｔ３を用い、Ｔu＝（Ｔ２＾２−Ｔ１Ｔ３）／（２Ｔ２−Ｔ１−Ｔ３）で平衡温度を予測するものが提案されている。なお、“＾”は冪乗を表している（特許文献３参照）。
【０００８】
また、本出願人は平衡温度の予測において、３つのパラメータ（ａ，ｂ，ｃ）を有する予測関数Ｔu＝ｃ・（Ｃ−ｔ）＋ａ・ｔ＾ｂを用いることによって、平衡予測温度を求めるサンプリング数を少なくし、外的要因や個人差等の影響を低減する電子体温計を提案している。
【０００９】
この電子体温計では、３つのパラメータ（ａ，ｂ，ｃ）を有する予測関数Ｔu＝ｃ・（Ｃ−ｔ）＋ａ・ｔ＾ｂを用いて平衡温度を予測する場合は、測定点から熱平衡点までの間を直線的に補正するパラメータｃを決定する際、予測関数の曲線の形状を定めるパラメータｂと一測定点の測定温度とパラメータｃとの関係を表したテーブルを用意しておき、パラメータｂと測定温度からパラメータｃを求めることができる。
【００１０】
図１０は、パラメータｂと測定点の測定温度に基づいてパラメータｃを定めるテーブルを説明するための図である。図１０において、縦軸をパラメータｂとして横軸を測定温度Ｔ１としている。パラメータｃは、パラメータｂと測定温度Ｔ１とを結ぶ線によって求めることができる（特許文献４参照）。
【００１１】
【特許文献１】特公平２−５９４１８号公報
【特許文献２】特公平４−７４５６号公報
【特許文献３】特開昭５５−７１９１９号公報
【特許文献４】特開２００６−３０８５４０号公報
【発明の開示】
【発明が解決しようとする課題】
【００１２】
前記した複数の予測関数の中から使用する予測関数を測定毎に選択する場合には、予測のための制御が複雑となる他、予め複数の予測関数を用意しておく必要があるといった構成上の問題がある他に、適切な予測関数が選択されない場合には、予測値が求められないという問題がある。
【００１３】
また、回帰法により算出を行う場合には、サンプリング数を増やすことで精度を向上させることができるが、予測値を決定するまでに時間がかかるという問題がある。これに対して、サンプリング数を減らして測定時間を短縮しようとした場合には、正確な予測値を得ることができなくなるという問題がある。
【００１４】
また、Ｔu＝（Ｔ２＾２−Ｔ１Ｔ３）／（２Ｔ２−Ｔ１−Ｔ３）という予測関数は、Ｔu−Ｔ＝（Ｔu−Ｔ０）ｅ＾（−ｔ／τ）という理論式を元としている（なお、Ｔuは熱平衡温度、Ｔ０は初期温度、Ｔは検出温度、ｔは時間、τは熱時定数、“＾”は冪乗である）。この理論式は理想的な体温上昇をモデルとするものであって、得られた予測関数は被測定対象に依存する自由度を持たない関数であるため、外的要因や個人差等により必ずしも正確な予測値を決定することができないという問題がある。
【００１５】
また、３つのパラメータ（ａ，ｂ，ｃ）を有する予測関数Ｔu＝ｃ・（Ｃ−ｔ）＋ａ・ｔ＾ｂを用いることによって平衡予測温度を求める場合には、測定時における外気温が低いなどの環境条件によって測定開始時における測定温度と体温との温度差が大きい場合には、パラメータを決定するための複数の測定点の間で急な温度上昇傾向となるため、予測関数の曲線の形状を決めるパラメータｂの値が大きくなって予測値が高くなる傾向があり、予測誤差が増大するという問題がある。
【００１６】
図１１は、測定開始時における測定温度と体温との温度差による測定誤差を説明するための図である。図１１において、実線は実測値を示し、破線は測定開始時における測定温度と体温との温度差が小さい場合の予測関数Ａを示し、一点鎖線は測定開始時における測定温度と体温との温度差が大きい場合の予測関数Ｂを示している。測定開始時における測定温度と体温との温度差が大きい場合の予測関数Ｂは、パラメータｂの値が大きくなる傾向にあるため予測値が高くなる。
【００１７】
そこで、本発明は従来の問題を解決し、測定開始時における測定温度と体温との温度差が大きい場合における予測誤差を低減することを目的とする。
【課題を解決するための手段】
【００１８】
本発明の電子体温計は、予測関数によって熱平衡状態に至った平衡温度を予測するものであり、測定開始時における測定温度と体温との温度差が大きい場合には、予測値を補正することによって高温側に生じる予測誤差を低減する。
【００１９】
本発明では、測定開始時における測定温度と体温との温度差が大きい場合に生じる予測誤差の要因は、温度差によって測定温度の増加特性に影響が生じるためにあるとして、測定温度の増加特性に基づいて補正値を設定することによって予測誤差を低減する。
【００２０】
本発明の電子体温計は、被測定部位の温度を検出する温度検出手段と、検出した温度に応じて平衡温度を算出する予測手段とを備える。
【００２１】
本発明の予測手段は、２つのパラメータａ，ｂおよび補正値αを有する平衡予測温度を求める予測関数Ｔu＝α・（Ｃ−ｔ）＋ａ・ｔ＾ｂ（Ｔu：熱平衡予測温度、ｔ：測定開始からの測定時間、ａ：予測関数の曲線の倍率を定めるパラメータ、ｂ：予測関数の曲線の形状を定めるパラメータ、α：予測関数を直線状に補正する補正値、Ｃ：時間定数、＾：冪乗）のパラメータａ，ｂおよび補正値補正値αを求めて予測関数を決定する予測関数決定部と、予測関数決定部により決定された予測関数を用いて熱平衡時間時における平衡温度を算出する熱平衡予測温度算出手段とを備える。
【００２２】
予測手段は、予測関数決定部によって測定開始から所定の時間が経過した測定点のデータを用いて予測関数を決定し、予測関数を用いて熱平衡予測温度算出手段によって熱平衡時間時における平衡温度を算出する。
【００２３】
本発明の予測関数決定部は、予測関数の曲線の倍率を定めるパラメータａおよび予測関数の曲線の形状を定めるパラメータｂを算出するパラメータ算出部と、予測関数の測定開始から熱平衡時点までを直線状に補正する補正値αを、測定温度の増加特性を判定条件として算出する補正値算出部とを備え、パラメータ算出部は算出したパラメータａ，ｂおよび補正値算出部で算出した補正値αに基づいて予測関数を決定する
【００２４】
測定温度の増加特性は、外気温と被測定者の体温の温度差によって影響を受ける。本発明は、この測定温度の増加特性を、予測関数の曲線の形状を定めるパラメータｂの値と測定開始時から温度変化が安定するまでの時間長によって特定する他、測定温度の増加特性を測定温度の温度変化幅によって特定する。
【００２５】
パラメータｂは予測関数の曲線の形状を定めるパラメータであり、このパラメータｂの値が大きいほど、測定温度が増加する割合が時間経過に対して大きいことを表している。また、測定開始時から温度変化が安定するまでの時間長は、予測関数の立ち上がりの急峻さを表す指標となり、この時間長は長いほど、温度差が大きいことを表している。
【００２６】
したがって、外気温と被測定者の体温の温度差が大きいほど、パラメータｂの値は大きく、かつ、安定するまでの時間長は長くなる。
【００２７】
補正値算出部は、測定温度の増加特性を予測関数の曲線の形状を定めるパラメータｂの値と測定開始時から温度変化が安定するまでの時間長とによって特定する場合には、パラメータｂの値と時間長とを判定条件とする補正値αを予め定めて備える。補正値算出部は、パラメータ算出部が算出するパラメータｂの値と、温度検出手段が検出する温度の変化から求めた時間長とに基づいて判定条件から補正値αを算出する。
【００２８】
一方、測定温度の温度変化幅は、測定開始時の測定温度と平衡状態に至ったあるいは平衡状態に至ったと見なせるときの測定温度との温度差である。平衡状態における測定温度はほぼ一定の範囲内にあるため、測定開始時における測定温度を測定温度の温度変化幅と見なすことができ、平衡状態における測定温度を測定することなく、測定開始時における測定温度を測定温度の温度変化幅を表す指標として用いることができる。
【００２９】
補正値算出部は、測定温度の増加特性を測定温度の温度変化幅によって特定する場合には、測定温度の温度変化幅を表す値として測定開始時の測定温度を用いる。補正値算出部は、この測定開始時の測定温度を判定条件とする補正値αを予め定めて備える。補正値算出部は、温度検出手段が検出する測定開始時の測定温度に基づいて判定条件から補正値αを算出する。
【００３０】
本発明のパラメータ算出部は、補正値αを零とし、２測定点の検出温度およびこの検出温度の測定時間を用いてパラメータｂを決定すると共に仮パラメータａ′を仮決定し、さらに、パラメータｂ及び補正値α、並びに、２測定点の何れか一方の測定点の検出温度及びこの検出温度の測定時間を予測関数に代入してパラメータａを再決定する。
【００３１】
パラメータ算出部において、パラメータｂは、ｂ＝ln（Ｔ１／Ｔ２）／ln（ｔ１／ｔ２）（Ｔ１、Ｔ２は２測定点の検出温度、ｔ１、ｔ２は測定時間、lnは対数）の演算式で算出することができ、仮決定される仮パラメータａ′は、lnａ′＝（lnＴ２・lnｔ１−lnＴ１・lnｔ２）／（lnｔ１−lnｔ２）（Ｔ１、Ｔ２は２測定点の検出温度、ｔ１、ｔ２は測定時間、lnは対数）の演算式で算出することができ、再決定されるパラメータａは予測関数Ｔu＝α・（Ｃ−ｔ）＋ａ・ｔ＾ｂの演算式で算出することができる。
【００３２】
本発明の予測手段は、温度検出手段により検出された温度に基づいて温度勾配を算出する温度勾配算出部を備える。この温度勾配算出部は、検出温度の温度勾配が所定範囲内となった時点の検出温度を２つもしくは３つの検出温度の内の最初の第１点目の検出温度を検出する測定時間とする。
【００３３】
また、本発明の予測手段は計時部を備える。この計時部は、第１点目の温度検出からの時間を計時して、第１点目の測定時間から所定時間後における第２点目乃至第３点目の測定時間とする。
【００３４】
本発明の熱平衡予測温度算出手段は、パラメータａ，パラメータｂ、及び補正値αとにより定まる予測関数の測定時間ｔに熱平衡予測時間を代入して熱平衡予測温度を算出する。予測関数Ｔuが備える演算時定数Ｃは熱平衡予測時間とすることができ、熱平衡予測時間時における平衡温度を求める際にｔ＝Ｃとすることで、熱平衡予測温度を算出することができる。時定数Ｃを熱平衡予測時間とし、例えば、Ｃ＝６００秒とすることによって測定開始から６００秒後の予測時間を算出することができる。
【００３５】
また、本発明の電子体温計は、予測関数を定めるパラメータを検出温度とその測定時間を用いて決定することによって、外的要因や個人差等に応じた予測関数となり、外的要因や個人差によるばらつきを低減することができる。
【００３６】
さらに、予測関数のａ・ｔ＾ｂ部で温度上昇曲線を近似し、一次関数であるα・（Ｃ−ｔ）部でその曲線を補正するという単純な式で熱平衡時点の温度を正確に予測でき、予測手段をシンプルに設計することができる。
【発明の効果】
【００３７】
本発明の電子体温計は、測定開始時における測定温度と体温との温度差が大きい場合における予測誤差を低減することができる。
【発明を実施するための最良の形態】
【００３８】
以下、本発明の電子体温計について図を用いて詳細に説明する。
【００３９】
図１は本発明の電子体温計の概略構成を説明するための図である。
【００４０】
電子体温計は、被測定部位の温度を検出する温度検出手段１と、温度検出手段１で検出した実測値を用いて熱平衡時の温度を算出する予測手段２と、予測手段２で予測した予測値を表示する表示手段３と、予測値表示を音で告知するブザー４と、前回の測定値あるいは予測値を記憶する前回値記憶手段５と、表示手段３やブザー４の他に電子体温計の駆動を制御する制御手段６を備える。
【００４１】
温度検出手段１は、被測定部位の温度を測定する温度センサ１ａと、温度センサ１ａから出力される検出信号を温度信号に変換して実測値として出力する温度測定部１ｂを備える。温度測定部１ｂはＡ／Ｄ変換器を備え、実測値をデジタル値で出力する。
【００４２】
表示手段３は、予測値や実測値を表示する表示部３ｂと、表示部３ｂに対して予測値の表示から実測値の表示に切り替え等の表示切り替えを制御する表示切替部３ａとを備える。
【００４３】
表示部３ｂは、予測手段２で予測した予測値を表示する他に、前回値記憶手段５に記憶する前回実測値あるいは前回予測値や、測定値が安定した後に予測値から実測値に切り替える場合にはその実測値を表示する。また、予測中であること、予測値の表示中であること、あるいは実測値の表示中であること等の表示状態についても、必要に応じて表示することができる。
【００４４】
制御手段６は、電子体温計全体の制御を行う他に、測定開始から所定時間（例えば、３分）が経過したか否かの判定を行い、表示切替部３ａに対して表示切替の制御を行う。このとき、予測表示から実測表示への切替制御は、測定値が安定したか否かの安定検出により行なっても良い。また、予測値表示時や予測値表示から実測値への表示切り替え時等に、ブザー４を駆動して音で告知を行う。なお、予測手段２については、図３を用いて後述する。
【００４５】
本発明の電子体温計の動作例を、図２のフローチャートを用いて説明する。なお、以下のフローチャートにおいて、S7を除く各工程は、予測温度を算出する従来の電子体温計の動作と同様とすることができる。また、各工程の制御は、ＣＰＵ及び制御プログラムが記憶されたメモリ等で構成される制御手段によって行われる。
【００４６】
電子体温計のスイッチがオンされると、表示手段３の表示セグメントを全点灯表示させる等の表示動作（S1）、及びブザーによる音動作を行って、スイッチがオン状態となったことを表示する（S2）。なお、この全点灯表示は一例であって、他の表示形態とすることもできる。表示手段３は、全点灯表示を行った後、前回値記憶手段５に記憶しておいた前回の測定値を読み出して表示する（S3）。
【００４７】
温度検出手段１は温度測定を開始し、検出した実測値を予測手段２及び制御手段６に送る（S4）。温度測定が開始されると、表示手段３は予測中であることを表示する（S5）。予測手段２は、温度検出手段１で検出された実測値を順次入力し、予め設定された時点での実測値を予測演算用データとして取り込んで（S6）、予測値の算出処理を行う（S7）。
【００４８】
なお、予測演算用データとして取り込むデータ数や、取り込むタイミングは種々設定することができる。図３は予測演算用データの一例を説明するための図である。図３（ａ）は予測関数のパラメータを決定するために用いる予測演算用データ（Ｔ１，Ｔ２）の算出を説明するための図である。図３（ａ）において、時間ｔ０で測定を開始した後、温度上昇率から第１番目のデータ点（時間ｔ１）を選択し、第１番目から所定時間経過した時点を第２のデータ点（時間ｔ２）を選択する等によって２点の予測演算用データ（［ｔ１，Ｔ１］，［ｔ２，Ｔ２］）を抽出し、これらの予測演算用データ（［ｔ１，Ｔ１］，［ｔ２，Ｔ２］）を所定の予測用の演算式に適用させることで予測値を算出する。
【００４９】
また、図３（ｂ）は予測関数の補正値を決定するために、測定温度の増加特性の時間変化の傾向を求めるために用いる予測演算用データ（［ｔＡ，ＴＡ］〜［ｔＤ，ＴＤ］）の算出を説明するための図である。図３（ｂ）において、例えば、［ｔＡ，ＴＡ］と［ｔＢ，ＴＢ］の間、［ｔＢ，ＴＢ］と［ｔＣ，ＴＣ］の間、［ｔＣ，ＴＣ］と［ｔＤ，ＴＤ］の間の測定温度の増加特性の時間変化の傾きをパラメータｂによって求め、求めたパラメータｂの増減から測定温度の増加特性の時間変化の傾向を求める。
【００５０】
求めた測定温度の増加特性の時間変化の傾向から、予測関数を直線状に補正する補正値を定める。
【００５１】
予測手段２は、S7の工程で算出した予測値を表示手段３に送って、予測値を表示すると共に（S8）、ブザー６を鳴らすことによって利用者に予測値が表示されたことを告知する（S9）（予測値表示モード）。
【００５２】
時間Ｔ０で測定を開始してから３分が経過したか否かの判定を行い（S10）、表示切替部３ａに対して表示切替の制御を行って、表示手段３に実測値を表示する（実測値表示モード）。ここで、表示手段３は実測値を表示し（S11）、ブザー６を鳴らすことによって利用者に安定状態となったことを告知する（S12）。
【００５３】
予測表示から実測表示への切替制御は、測定値が安定したか否かの安定検出により行なっても良い。この場合には、制御手段６は検出温度の変化が安定状態に達したことを検出することによって温度が熱平衡状態にあると判定する。
【００５４】
次に、本発明の電子体温計が備える予測手段について、図４を用いて予測手段の構成例を説明し、図５のフローチャート及び図６〜図９の温度特性曲線を用いて予測動作について説明する。
【００５５】
図４において、予測手段２は、予測演算用データ保持部２ａと、予測関数決定部２ｂと、温度勾配算出部２ｃと計時部２ｄと、熱平衡予測温度算出部２ｅを備える。予測関数決定部２ｂと熱平衡予測温度算出部２ｅは予測演算部を構成している。
【００５６】
温度勾配算出部２ｃは、温度検出手段１からの実測値を入力し、この上昇率から温度勾配を算出し、この温度勾配が所定範囲内となった時点を第１のデータ点を取得する時間ｔ１とし、この時の実測値Ｔ１と共に予測演算用データ保持部２ａに保持させる。また、計時部２ｄは、温度勾配算出部２ｃで定めた時間ｔ１から計時を開始し、予め設定した時間が経過した時点を第２のデータ点を取得する時間ｔ２とし、この時の実測値Ｔ２と共に予測演算用データ保持部２ａに保持させる。これによって、予測演算用データ保持部２ａには、２つの予測演算用データ［ｔ１，Ｔ１］，［ｔ２，Ｔ２］が記憶される。
【００５７】
また、予測関数決定部２ｂは、パラメータ算出部２ｂ１と補正値算出部２ｂ２とを備える。パラメータ算出部２ｂ１は、パラメータ決定部２ｂ１１とパラメータ再決定部２ｂ１２とを備える。
【００５８】
パラメータ決定部２ｂ１１は、２つの予測演算用データ（［ｔ１，Ｔ１］，［ｔ２，Ｔ２］）を用いて、予測関数が有するパラメータの内のパラメータａを仮決定すると共に、パラメータｂを決定する。なお、以下、仮決定したパラメータａをパラメータａ′で表す。
【００５９】
また、パラメータ再決定部２ｂ１２は、１つの予測演算用データ（［ｔ２，Ｔ２］）と、パラメータ決定部２ｂ１１で決定したパラメータｂと、補正値算出部２ｂ２で算出した補正値αと用いて、予測関数のパラメータａを再決定する。
【００６０】
補正値算出部２ｂ２は、測定温度の増加特性に基づく判定条件によって補正値αを算出する。判定条件は、予測関数の曲線の形状を定めるパラメータｂの値と測定開始時から温度変化が安定するまでの時間長ｔ１との組み合わせとする他、測定開始時の測定温度Ｔ１とすることができる。
【００６１】
補正値算出部２ｂ２は、予測関数の曲線の形状を定めるパラメータｂの値と測定開始時から温度変化が安定するまでの時間長ｔ１との組み合わせ、あるいは測定開始時の測定温度Ｔ１を判定条件とし、この判定条件に対応する補正値αを予め定めて記憶しておき、パラメータ決定部２ｂ１１から得られたパラメータｂや予測演算用データ保持部２ａから得られた測定温度Ｔ１に対応する補正値αを求める。
【００６２】
熱平衡予測温度算出部２ｅは、予測関数決定部２ｂで決定したパラメータａ，ｂおよび補正値αを用いて決定された予測関数を用いて熱平衡予測温度を算出し、予測温度Ｔｕを出力する。
【００６３】
予測関数決定部２ｂおよび熱平衡予測温度算出部２ｅからなる予測演算部による演算動作を、図５に示すフローチャートを用いて説明する。なお、図５のフローチャートにおいて、S101〜S108は、予測演算用データ保持部２ａ、温度勾配算出部２ｃ，計時部２ｄが行う処理である。
【００６４】
図５のフローチャートにおいて、計時部２ｄは、測定開始を検出すると（S101）、ｔ＝０として計時を開始する（S102）。また、温度勾配算出部２ｃは、実測値を入力してその上昇率を算出し、その上昇率が所定範囲であるか否かを判定する。上昇率は、例えば、所定時間毎に入力する実測値の差分値が、所定の温度範囲にあるか否かで判定することができる。この上昇率の判定は、前回入力した実測値をＴ（ｔ−１）とし今回入力した実測値をＴ(ｔ)としたとき、その差分Ｔ(ｔ)−Ｔ（ｔ−１）が、例えば０．００℃と０．０２℃の温度範囲内にあるか否かを判定することで行うことができる（S103）。
【００６５】
本例では温度上昇率が所定範囲内となった後、所定時間ｔｂ（例えば１秒後）後に（S104）、改めて入力された実測値Ｔ(ｔ)と前回入力された実測値Ｔ（ｔ−１）の差分Ｔ(ｔ)−Ｔ（ｔ−１）が再度所定温度範囲内（例えば０．００℃と０．０２℃）にあるか否かを判定することで、上昇率に大きな変動がないことを確認し、次のステップに移る（S105）。これは、正常な温度上昇によらず、何らかの原因で上昇率の変化が所定範囲に入ってしまうようなケースを除外することによって、温度上昇率が確実に安定の方向に移行していること確認するためである。仮に、差分Ｔ(ｔ)−Ｔ（ｔ−１）が再度所定温度範囲内と成らないときは、前回のS103での温度上昇の判定は誤りであったと判断してS103に戻り、再び所定範囲内に入るまで判定を繰り返す（S105）。
【００６６】
S105の工程において、温度差が所定範囲内と成った場合には、その時の時間ｔ１と実測値Ｔ１を予測演算用データ保持部２ａに保持して、第１のデータ点［ｔ１，Ｔ１］を決定する（S106）。
【００６７】
計時部２ｄは、時間の経過を監視し、時間ｔ１から所定時間ｔｃが経過すると、時間ｔ２を予測演算用データ保持部２ａに通知する（S107）。予測演算用データ保持部２ａは、時間ｔ２とそのときの実測値Ｔ２を保持して、第２のデータ点［ｔ２，Ｔ２］を決定する（S108）。第２のデータ点［ｔ２，Ｔ２］を定める所定時間ｔｃは、例えば１０秒とすることができる。図６は２つのデータ点［ｔ１，Ｔ１］，［ｔ２，Ｔ２］を示している。
【００６８】
パラメータ決定部２ｂ１１は、第１のデータ点［ｔ１，Ｔ１］と第２のデータ点［ｔ２，Ｔ２］を用いて、予測関数のパラメータａ及びパラメータｂを算出する。ここで、熱平衡予測温度Ｔuの予測関数は、
予測関数：Ｔu＝α・（Ｃ−ｔ）＋ａ・ｔ＾ｂ    …（１）
Ｔu：熱平衡予測温度
ｔ：測定開始からの測定時間
ａ：予測関数の曲線の倍率を定めるパラメータ
ｂ：予測関数の曲線の形状を定めるパラメータ
α：予測関数を直線状に補正する補正値、
Ｃ：時間定数
＾：冪乗
を備える。
【００６９】
予測関数Ｔuは、（ａ・ｔ＾ｂ）  で表される時間ｔと共に増加率が減少する項と、α・（Ｃ−ｔ）で表される時間ｔと“α・Ｃ”（時間ｔ＝０）の値から“０”（時間ｔ＝Ｃ）まで直線的に減少する補正項に相当する項とを備える演算式で表される。
【００７０】
この時のパラメータａとパラメータｂの算出は、補正値αを零とすることで得られ、パラメータｂを決定すると共に、パラメータａ´を仮決定する。仮決定したパラメータａ´は後に再度求めて、パラメータａを決定する。  
【００７１】
上記した予測関数に、第１のデータ点［ｔ１，Ｔ１］と第２のデータ点［ｔ２，Ｔ２］を代入すると、パラメータａ´は、
lnａ´＝（lnＴ２・lnｔ１−lnＴ１・lnｔ２）／（lnｔ１−lnｔ２）  …（２）
として仮決定される。ここで、仮決定のパラメータａをパラメータａ´で表している。“ln”は対数関数を表している。
【００７２】
また、パラメータｂは、
ｂ＝ln（Ｔ１／Ｔ２）／ln（ｔ１／ｔ２）  …（３）
により決定される。この工程によって、仮パラメータａ´とパラメータｂが定まる。図７（ａ）はこの段階における予測関数を示している（S109）。
【００７３】
補正値算出部２ｂ２は、パラメータ決定部２ｂ１１からパラメータｂの値（ｂ値）および第１のデータ点［ｔ１，Ｔ１］、あるいは測定開始時における測定温度Ｔ１を取得し、これらの値を判定条件と照合して判定を行い、判定条件に対応する補正値αを求める。
【００７４】
表１、２は、判定条件を説明するための表である。
【００７５】
【表１】
2010122164.tif 000003
【００７６】
表１は、パラメータｂの値（ｂ値）と第１のデータ点［ｔ１，Ｔ１］のｔ１とを組み合わせた判定条件の概要例を示している。
【００７７】
表１に示す例では、ｂ値が予め定めた設定値ｂ１以上であり、かつ、ｔ１値が予め定めた設定値ｔａと設定値ｔｂで挟まれる範囲内の場合には補正値αａを設定し、ｂ値が予め定めた設定値ｂ１以上であり、かつ、ｔ１値が予め定めた設定値ｔｂ以上である場合には補正値αｂを設定している。
【００７８】
図８は、補正値αａおよび補正値αｂを設定する例を示している。図８中の四角の印および丸印で示される測定点は、ｂ値がｂ１以上の値であって、温度の差分（Ｔ(ｔ)−Ｔ（ｔ−１））が所定温度範囲内となるｔ１値をとる点を表している。四角の印で表される測定点はｔ１値が設定値ｔａと設定値ｔｂで挟まれる範囲内となる場合を示し、この範囲内を通る測定値の予測関数に対しては補正値αａを設定する。
【００７９】
また、図８中の丸印で表される測定点はｔ１値が設定値tｂ以上の範囲内となる場合を示し、この範囲内を通る測定値の予測関数に対しては補正値αｂを設定する。
【００８０】
表２は、測定開始時の測定温度の判定条件の概要例を示している。
【表２】
2010122164.tif 000004
【００８１】
表２に示す例では、測定開始時の測定温度Ｔ０の値が予め定めた設定値Ｔａと設定値Ｔｂで挟まれる範囲内の場合には補正値αＡを設定し、測定開始時の測定温度Ｔ０の値が予め定めた設定値ＴＢと設定値ＴＣで挟まれる範囲内の場合には補正値αＢを設定している。
【００８２】
図９は、補正値αＡおよび補正値αＢを設定する例を示している。図８において、測定開始時（ｔ０）の測定温度Ｔ０がＴＡとＴＢの間に挟まれる測定値の予測関数に対しては補正値αＡを設定し、測定開始時（ｔ０）の測定温度Ｔ０がＴＢとＴＣの間に挟まれる測定値の予測関数に対しては補正値αＢを設定する。
【００８３】
補正値算出部２ｂ２は、判定条件とその判定条件に対応する補正値αを予め定めて記憶しておき、パラメータ決定部２ｂ１１から得られたパラメータｂ、ｔ１値、あるいは測定温度Ｔ０基づいて補正値αを求める。
【００８４】
図７（ａ）は補正値αによる補正を行う前の予測関数を示し、図７（ｂ）はこの段階における予測関数を示している。図７では、破線は温度差がないときに得られる予測関数を示し、実線は温度差があるときの予測関数を示している。また、薄い実線は補正値αで補正した後の状態を示している（S110）。
【００８５】
次に、パラメータａが仮決定のパラメータａ´であるため、前記式（１）で表される予測関数Ｔuに第２のデータ点［ｔ２，Ｔ２］を代入することによってパラメータａを再決定する。ここで、式（１）で表される予測関数Ｔuに代入するデータ点は第１のデータ点［ｔ１，Ｔ１］を用いることもできるが、第２のデータ点［ｔ２，Ｔ２］は平衡温度に近い測定点であるため、この平衡温度に近い測定点のデータを用いることによって予測関数の精度が高まることが期待される（S111）。
【００８６】
これによって、式（１）で表される予測関数Ｔuのパラメータａ，ｂおよび補正値αが定まる。図７（ｃ）中では、予測関数Ｔuは薄い破線で表される。
【００８７】
熱平衡予測温度は、式（１）の予測関数においてｔ＝Ｃ（Ｃは時間定数）とすることで算出することができる（S112）。算出した熱平衡予測温度Ｔuは表示部３ｂに表示する（S113）。
【００８８】
次に、本発明による補正値αの設定例を表３を用いて説明する。
【００８９】
表３は本発明による補正値αの設定例を示している。
【表３】
2010122164.tif 000005
【００９０】
表３は、ｂ値とｔ１値とを判定条件として補正値αを設定する例であり、前記表１で示した例に対応するものである。
【００９１】
例（ａ）では、パラメータｂの値（ｂ値）が０．０１５以上であり、かつｔ１値が４０sec以上の条件を満たす場合に、補正値αとして０．００１０を設定する例を示している。また、例（ｂ）では、パラメータｂの値（ｂ値）が０．０１５以上であり、かつｔ１値が３０secと４０secの範囲内の条件を満たす場合に、補正値αとして０．０００５を設定する例を示している。
【００９２】
この態様によれば、パラメータｂの値（ｂ値）とｔ１値とを判定条件として補正値αを設定することによって、あるいは、測定開始時の測定温度Ｔ０を判定条件として補正値αを設定することによって、測定開始時における測定温度と体温との温度差が大きい場合であっても予測誤差を低減することができる。
【００９３】
以下、本発明の予測による予測誤差と、特許文献４に示す従来の予測による予測誤差とを表４に示す。ここで、測定条件は１０℃の温度環境による測定であり、予測誤差は「（予測値）−（１０分後の実測値）」で表すものとする。
【００９４】
【表４】
2010122164.tif 000006
【００９５】
上記比較結果から、本願発明による予測誤差は０．０℃〜０．１℃であるのに対して、特許文献４に示す予測関数を用いた場合の予測誤差は０．２℃〜０．４℃であり、本願発明による測定誤差は、従来の予測関数を用いた場合の予測誤差よりも減少していることを確認することができる。
【図面の簡単な説明】
【００９６】
【図１】本発明の電子体温計の概略構成を説明するための図である。
【図２】本発明の電子体温計の動作例を説明するためのフローチャートである。
【図３】予測演算用データの一例を説明するための図である。
【図４】本発明の予測手段の概略構成を説明するための図である。
【図５】本発明の電子体温計の予測温度を求める動作を説明するためのフローチャートである。
【図６】本発明の電子体温計の予測温度を求める動作を説明するための温度特性曲線図である。
【図７】本発明の電子体温計の予測温度を求める動作を説明するための温度特性曲線図である。
【図８】本発明の電子体温計の予測温度を求める動作を説明するための温度特性曲線図である。
【図９】本発明の電子体温計の予測温度を求める動作を説明するための温度特性曲線図である。
【図１０】従来のパラメータｂと測定点の測定温度に基づいてパラメータｃを定めるテーブルを説明するための図である。
【図１１】温度差による予測誤差を説明するための図である。
【符号の説明】
【００９７】
１  温度検出手段
１ａ  温度センサ
１ｂ  温度測定部
２  予測手段
２ａ  予測演算用データ保持部
２ｂ  予測関数決定部
２ｂ１  パラメータ算出部
２ｂ１１  パラメータ決定部
２ｂ１２  パラメータ再決定部
２ｂ２  補正値算出部
２ｅ  熱平衡予測温度算出部
２ｃ  温度勾配算出部
２ｄ  計時部
３  表示手段
３ａ  表示切替部
３ｂ  表示部
４  ブザー
５  前回値記憶手段
６  制御手段
【図１】
2010122164.tif 000007
【図２】
2010122164.tif 000008
【図３】
2010122164.tif 000009
【図４】
2010122164.tif 000010
【図５】
2010122164.tif 000011
【図６】
2010122164.tif 000012
【図７】
2010122164.tif 000013
【図８】
2010122164.tif 000014
【図９】
2010122164.tif 000015
【図１０】
2010122164.tif 000016
【図１１】
2010122164.tif 000017