(57)【要約】【課題】予測関数を用いた測定時の熱平衡予測温度の算出を高速かつ簡便に行うことができる温度計及び温度計の温度予測方法を提供する。  【解決手段】被測定部位の温度を検出する温度検出手段と、検出した温度に応じて平衡温度を算出する予測手段とを備え、予測手段は、ｘを測定開始からの測定時間、ｙを平衡予測温度、Ａ，Ｂをあらかじめ定められた関数として定まる値としたときに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂ（式１とする）からなる予測関数に基づいて熱平衡予測温度を算出する。これにより、測定時の熱平衡温度の算出過程が複雑化することなく、短時間で簡便に熱平衡予測温度の算出を行うことが可能となる。【選択図】図１
2015045602.tif 000002
【特許請求の範囲】
【請求項１】
被測定部位の温度を検出する温度検出手段と、
検出した温度に応じて平衡温度を算出する予測手段とを備え、
前記予測手段は、
ｙ＝Ａｌｏｇ（ｘ）＋Ｂ・・・（式１）
（ｘ：測定開始からの測定時間
ｙ：予測温度
Ａ，Ｂはあらかじめ定められた関数により定められる値）
からなる予測関数に基づいて予測温度を算出する
ことを特徴とする温度計。
【請求項２】
パラメータＡは、所定時間ｘ１において実測値の曲線に式１の曲線が沿うように設定したＡの値Ａ１と、
別の所定時間ｘ２において実測値の曲線に式１の曲線が沿うように設定したＡの値Ａ２と、の関係に基づいて決定する
ことを特徴とする請求項１に記載の温度計。
【請求項３】
パラメータＢは、前記所定時間ｘ１において、実測値の曲線に式１の曲線が沿うように設定し、前記所定時間ｘ１における実測値の値をとるように設定したときのＢの値をＢ１、前記所定時間ｘ２において、実測値の曲線に式１の曲線が沿うように設定し、前記所定時間ｘ１における実測値の値をとるようにオフセットさせたときのＢの値をＢ２としたときの、Ｂ２とＢ１との差であるＢ３の値と、
前記所定時間ｘ２において、実測値の曲線に式１の曲線が沿うように設定し、所定時間ｘ１における実測値の値をとるようにオフセットさせた曲線の平衡状態における値と、平衡状態における実測値との差Δｙと、
の関係に基づいて決定する
ことを特徴とする請求項２に記載の温度計。
【請求項４】
Ａ１とＡ２との関係を２次式で規定することを特徴とする請求項２に記載の温度計。
【請求項５】
Ｂ３とΔｙとの関係を１次式で規定することを特徴とする請求項３に記載の温度計。
【請求項６】
被測定部位の温度を検出する温度検出工程と、
検出した温度に応じて平衡温度を算出する予測工程とを備え、
前記予測工程は、
ｙ＝Ａｌｏｇ（ｘ）＋Ｂ・・・（式１）
（ｘ：測定開始からの測定時間
ｙ：予測温度
Ａ，Ｂはあらかじめ定められた関数により定められる値）
からなる予測関数に基づいて予測温度を算出する
ことを特徴とする温度計の温度予測方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、温度計に関し、特に平衡温度を算出する予測手段を備えた温度計、及び温度計の温度予測方法に関する。
【背景技術】
【０００２】
温度計、特に被測定者の体温を測定する体温計においては、測定開始直後には温度センサの温度と体温との間に温度差があるため、温度センサが被測定者の体温との間に熱平衡状態となるまでに時間がかかる。そのため、正確な体温を測定するのに時間を要する。
【０００３】
また、従来の体温計は、口内や腋下に体温計を保持して測定することが一般的であったが、年齢が低い場合や体調によっては、長時間体温計を保持しておくことが容易ではない場合もある。
【０００４】
そのため、近年温度検知素子としてサーモパイルを用い、被験者の額などに対して非接触に対向させて、温度の検出を行う体温計も用いられている。しかし、サーモパイルを用いた体温計は、測定する時間は短いものの、サーミスタを用いたものに比べてコストが高くなってしまう。
【０００５】
そこで、従来のサーミスタを用いた体温計などの温度計では、実測値の変化から熱平衡状態での平衡温度を予測して、この予測値を順次更新して表示することによって、測定時間を短縮することが行われている。
【０００６】
平衡温度の予測は、測定温度とその変化率を監視して、この測定温度と変化率を変数とする予測関数を演算することによって行われるのが一般的である。このような予測関数を用いて平衡温度を予測する体温計が提案されている。
【０００７】
特許文献１には、平衡温度を予測する予測関数を複数用意しておき、所定間隔で平衡温度を予測し、今回の平衡温度予測値と前回の平衡温度予測値との差が所定範囲外の場合には、新たな予測関数を選択する一方、所定範囲内の場合には、当初の予測関数に基づいて平衡温度予測値を表示しつつ、予測演算を続行する体温計の技術が開示されている。
【０００８】
また、特許文献２には、被測定部位の温度を検出する温度検出手段を備え、検出した温度を複数の所定時間においてサンプリングして、予測関数のパラメータを定めて熱平衡予測温度を算出する技術が開示されている。
【先行技術文献】
【特許文献】
【０００９】
【特許文献１】特公平２−５９４１８号公報（第４ページ、図１）
【特許文献２】特許第４７０２７８１号公報（第７ページ、図１）
【発明の概要】
【発明が解決しようとする課題】
【００１０】
特許文献１に示された技術は、複数の予測関数の中から使用する予測関数を測定ごとに選択するため、予測のための制御が複雑となる。また、あらかじめ複数の予測関数を用意しておく必要があるため構成が複雑となり、適切な予測関数が選択されない場合に、予測値を求めることができなくなるおそれがある。
【００１１】
また、特許文献２に示された技術は、被測定部位の温度を検出して、検出した温度を所定時間においてサンプリングしつつ、予測関数のパラメータを定めているため、サンプリング点における検出温度によってパラメータを定めながら、熱平衡予測温度を算出することとなり、熱平衡予測温度の算出過程が複雑となり、時間がかかってしまうおそれがある。
【００１２】
本発明の目的は、予測関数を用いた測定時の熱平衡予測温度の算出を高速かつ簡便に行うことができる温度計及び温度計の温度予測方法を提供することである。
【課題を解決するための手段】
【００１３】
本発明の温度計は、被測定部位の温度を検出する温度検出手段と、検出した温度に応じて平衡温度を算出する予測手段とを備え、予測手段は、ｘを測定開始からの測定時間、ｙを予測温度、Ａ，Ｂをあらかじめ定められた関数により定められる値としたときに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂ（式１とする）からなる予測関数に基づいて予測温度を算出する。
【００１４】
このような構成により、予測関数として対数関数を用いた上で、そのパラメータを、所定時間におけるサンプリング結果によって計算しながら変更するのではなく、あらかじめ定められた関数として定めるようにすることによって、測定時の熱平衡温度の算出過程が複雑化することなく、短時間で簡便に熱平衡予測温度の算出を行うことが可能となる。
【００１５】
また、パラメータＡは、所定時間ｘ１において、実測値の曲線に式１の曲線が沿うように設定したＡの値Ａ１と、別の所定時間ｘ２において、実測値の曲線に式１の曲線が沿うように設定したＡの値Ａ２と、の関係に基づいて決定することができる。
【００１６】
このような構成により、異なる時間における実測値の温度勾配において、対数関数の温度勾配が等しくなるように設定した対数部の係数の関係に基づいて、予測関数の対数部の係数を決定することで、実測値の変化に対応して予測関数のカーブを決定することが可能となる。
【００１７】
さらに、パラメータＢは、所定時間ｘ１において、実測値の曲線に式１の曲線が沿うように設定し、所定時間ｘ１における実測値の値をとるように設定したときのＢの値をＢ１、所定時間ｘ２において、実測値の曲線に式１の曲線が沿うように設定し、前記所定時間ｘ１における実測値の値をとるようにオフセットさせたときのＢの値をＢ２としたときの、Ｂ２とＢ１との差であるＢ３の値と、所定時間ｘ２において、実測値の曲線に式１の曲線が沿うように設定し、所定時間ｘ１における実測値の値をとるようにオフセットさせた曲線の平衡状態における値と、平衡状態における実測値との差Δｙと、の関係に基づいて決定することができる。
【００１８】
このような構成により、所定時間ｘ１における温度勾配に対応させたときの予測関数の係数と、所定時間ｘ２における温度勾配に対応させたあと、所定時間ｘ１において実測値の値に等しくなるようにｙ軸方向に平行移動したときの係数との差異から、実測値の変化に対応した予測関数のオフセット値を決定することができる。
【００１９】
また、Ａ１とＡ２の関係は２次式で規定し、Ｂ３とΔｙとの関係は１次式で規定することがより好ましい。このようにすることにより、対数関数からなる予測関数をより正確に実測値の時間変化に合うようにすることができる。
【００２０】
本発明の温度計の温度予測方法は、被測定部位の温度を検出する温度検出工程と、検出した温度に応じて平衡温度を算出する予測工程とを備え、予測工程は、ｘを測定開始からの測定時間、ｙを予測温度、Ａ，Ｂをあらかじめ定められた関数として定まる値としたと
きに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂ（式１とする）からなる予測関数に基づいて熱平衡予測温度を算出する。
【００２１】
このような構成により、予測関数として対数関数を用いた上で、そのパラメータを、所定時間におけるサンプリング結果によって計算して変更するのではなく、あらかじめ定められた関数として定めるようにすることによって、熱平衡温度の算出過程が複雑化することなく、短時間で簡便に熱平衡予測温度の算出を行うことが可能となる。
【発明の効果】
【００２２】
本発明によれば、予測関数を用いた測定時の熱平衡予測温度の算出を高速かつ簡便に行うことができる温度計及び温度計の温度予測方法を提供することができる。
【図面の簡単な説明】
【００２３】
【図１】本発明の実施形態における体温計の構成の概略図である。
【図２】本発明の実施形態における体温計の動作例を示すフローチャートである。
【図３】本発明の実施形態における実測値と予測値との関係の一例を示したグラフである。
【図４】本発明の実施形態における実測値と、段階的な予測値との関係の一例を示したグラフである。
【図５】本発明の実施形態における予測式のパラメータの補正式を求めるための実験方法を示した図である。
【図６】本発明の実施形態における予測式のパラメータの補正式を求めるための実験方法を示した図である。
【図７】本発明の実施形態における実験条件の一覧表である。
【図８】本発明の実施形態における予測式のパラメータＡの値の補正式を求めるための計算処理を示すフローチャートである。
【図９】本発明の実施形態における予測式のパラメータＡの値の補正式を求めるための、実験結果の関係をプロットしたプロット結果を示した図である。
【図１０】本発明の実施形態における予測式のパラメータＢの値の補正式を求めるための計算処理を示すフローチャートである。
【図１１】本発明の実施形態における予測式のパラメータＢの値の補正式を求めるための、実験結果の関係をプロットしたプロット結果を示した図である。
【図１２】本発明の実施形態における予測手段の内部構成を示した図である。
【図１３】本発明の実施形態における平衡予測温度を算出する方法を示したフローチャートである。
【図１４】本発明の実施形態における実験例から得られる数値の一覧表である。
【図１５】本発明の実施形態における実験例の実験２、結果２の数値を用いて実測予測値と予測値との関係を示したグラフである。
【図１６】本発明の実施形態における実験例の実験２、結果３の数値を用いて実測予測値と予測値との関係を示したグラフである。
【図１７】本発明の実施形態における実験例の実験２、結果４の数値を用いて実測予測値と予測値との関係を示したグラフである。
【図１８】本発明の実施形態における実験例の実験２、結果５の数値を用いて実測予測値と予測値との関係を示したグラフである。
【図１９】本発明の実施形態における実験例の実験２、結果６の数値を用いて実測予測値と予測値との関係を示したグラフである。
【図２０】本発明の実施形態における実験例の実験２、結果７の数値を用いて実測予測値と予測値との関係を示したグラフである。
【図２１】本発明の実施形態における実験例の実験２、結果８の数値を用いて実測予測値と予測値との関係を示したグラフである。
【発明を実施するための形態】
【００２４】
以下、図面に基づいて、本発明の実施形態について説明する。本実施形態においては、温度計として体温計を用いている。
【００２５】
図１は、本実施形態における体温計の構成の概略図である。本実施形態の体温計は、被測定部位の温度を検出する温度検出手段１と、温度検出手段１で検出された実測値を用いて熱平衡時の温度を算出する予測手段２と、予測手段２で予測した予測値や温度検出手段１で検出された実測値を表示する表示手段３と、予測値が表示されたり、表示種別が変更されたときに報知するブザー４と、前回の実測値や予測値を記憶する前回値記憶手段５と、表示手段３、ブザー４、体温計の駆動を制御する制御手段６とを備えている。
【００２６】
温度検出手段１は、被測定部位の温度を測定する温度センサー１１と、温度センサー１１から出力される検出信号を温度信号に変換して実測値として出力する温度測定部１２とを備えている。温度測定部１２の内部には、図示しないＡ／Ｄ変換器を備えており、実測値はデジタル値に変換されて出力される。
【００２７】
表示手段３は、予測値や実測値を表示する表示部３２と、表示部３２の表示を予測値から実測値に切り替えるなどの表示切り替えを制御する表示切替部３１とを備えている。
【００２８】
表示部３２は、予測手段２で予測した予測値を表示したり、前回値記憶手段５に記憶されている前回実測値または前回予測値や、測定値が安定した後に予測値から実測値に切り替えて実測値を表示したりすることができる。また、現在予測中であることや、予測値の表示中であること、あるいは実測値の表示中であること等、現在の表示態様についても、必要に応じて表示することができる。予測手段２の構成については後述する。
【００２９】
ブザー４は、予測値表示時や予測値表示から実測値への表示切替を行った際などに駆動されて報知を行う。
【００３０】
制御手段６は、体温計全体の制御を行うとともに、測定値が安定したか否かの安定検出を行い、表示切替部３１に対して表示切替の制御を行う。また、ブザー４を駆動して報知を行う制御も行う。
【００３１】
本実施形態の体温計の動作例を、図２のフローチャートに基づいて説明する。図２のフローチャートに示されている各工程の制御は、制御手段６内の図示しないＣＰＵや、制御プログラムが記憶されたメモリ等で行われている。
【００３２】
体温計のスイッチがオンされると、表示手段３の表示部３２の表示が正常に動作していることを示すために、表示セグメントを全点灯表示させる等の表示動作が行われる（ステップＳ１０１）。また、ブザー４による報知を行って、ブザー４が正常に動作していることを報知するとともに、使用者にスイッチがオン状態となったことを報知する（ステップＳ１０２）。なお、本実施形態においては、スイッチがオンされたときに表示手段３は表示セグメントを全点灯表示するようにしているが、他の表示形態とすることもできる。表示手段３の表示部３２は、全点灯表示を行った後、前回値記憶手段５に記憶されている前回の測定値を読み出して表示する（ステップＳ１０３）。
【００３３】
続いて、温度検出手段１が温度測定を開始して、検出した実測値を予測手段２及び制御手段６に送信する（ステップＳ１０４）。温度検出手段１による温度測定が開始されると、表示手段３は予測中であることを表示する（ステップＳ１０５）。予測手段２は、予測演算用データの検出が開始されたかどうかを判定し（ステップＳ１０６）、予測値の算出
処理を行う（ステップＳ１０７）。予測手段２における具体的な予測値の算出処理については後述する。
【００３４】
予測手段２は、ステップＳ１０７において算出した予測値を表示手段３に送信して、表示手段３において予測値を表示する（ステップＳ１０８）。また、それと共に、制御手段６からの指令によりブザー４による報知を行い、利用者に予測値が表示されたことを報知する（ステップＳ１０９）。
【００３５】
温度検出において熱平衡状態に近づくと、制御手段６は検出温度の変化が安定状態に達したことを検出したかどうかを判定し（ステップＳ１１０）、安定状態に達すると、表示手段３は実測値を表示する（ステップＳ１１１）。さらに、ブザー４による報知を行うことによって、利用者に安定状態となって表示手段３が実測値を表示するようになったことを報知する（ステップＳ１１２）。
【００３６】
図３は、本実施形態における実測値と予測値との関係の一例を示したグラフである。縦軸は実測値または予測値の温度を示しており、横軸はサンプリング数を示している。グラフのうち、実線は実測値を示しており、破線は予測値を示している。図３においては、サンプリング間隔を０．５秒としているため、例えばサンプリング数６００の点は３００秒（５分）を、サンプリング数１２００の点は６００秒（１０分）を表している。
【００３７】
次に、図４に基づいて、本実施形態における予測値の算出方法について説明する。図４は、本実施形態における実測値と、段階的な予測値との関係の一例を示したグラフである。まず、実測Ｔｈ１は、実測した温度の値の時間変化を示したグラフである。予測１は、実測値をｙ＝Ａｌｏｇ（ｘ）＋Ｂの式で予測値を算出するにあたり、実測Ｔｈ１における第１所定時間（本実施形態においては１０秒時点）において、実測値の曲線に曲線が沿うように、ｙ＝Ａｌｏｇ（ｘ）＋ＢのＡの値を調整したグラフである。また、予測１．５は、実測Ｔｈ１における第２所定時間（本実施形態においては３分時点）において、実測値の曲線に曲線が沿うように、ｙ＝Ａｌｏｇ（ｘ）＋ＢのＡの値を調整した後に、Ｂの値を調整して、第１所定時間における実測値の値をとるようにしたグラフである。
【００３８】
さらに、予測２は、予測１におけるｙ＝Ａｌｏｇ（ｘ）＋ＢのＡの値を、事前に求めておいたＡの補正式に基づいて、Ａの値を補正し、また、Ｂの値を調整して、第１所定時間における実測値の値をとるようにしたグラフである。このＡの補正式の算出方法については後述する。また、予測３は、予測２におけるｙ＝Ａｌｏｇ（ｘ）＋ＢのＢの値を、事前に求めておいたＢの補正式に基づいて、Ｂの値を補正したグラフである。このＢの補正式の算出方法についても後述する。
【００３９】
次に、先ほどの予測式におけるＡ及びＢの補正式の算出方法について説明する。本実施形態においては、事前の実験結果から得られた値を用いて、Ａ及びＢの補正式を算出している。
【００４０】
図５及び図６は、本実施形態における補正式を求めるための実験方法を示した図であり、図５は、サーミスタ１０１を人体の額１１０に近づけた様子を示した図であり、図６は、サーミスタ１０１を水槽内の水１２１に浮かせたゴム１２２上に乗せる様子を示した図である。
【００４１】
図５の例では、サーミスタ１０１を挟むように、人体の額１１０側にＳＵＳ板１０２（厚さ０．１８ｍｍ）を設け、人体の額１１０とは反対側にスポンジ１０３を設けている。そして、スポンジ部を保持して、人体の額１１０に近づけて接触させるように移動させる。
【００４２】
図６の例では、サーミスタ１０１、ＳＵＳ板１０２、及びスポンジ１０３の構成は図５と同じであるが、スポンジ１０３の上に錘１０４を乗せ、水槽内の水１２１に浮かせたゴム１２２の上に乗せるようにする。水槽としては恒温水槽を用いて、中の水の温度を等しく保つようにしている。
【００４３】
図７に、本実施形態における実験条件を一覧で示している。実験１及び実験３は、図５の装置を用いて、人体の額１１０で行った実験例であり、実験２，４及び５は、図６の装置を用いて、スポンジ１０３の上の乗せる錘１０４の重量を変えたり、水槽内の水１２１の温度を変えたりして行った実験例である。
【００４４】
これらの実験例を用いて、温度測定開始時点から測定温度が平衡状態となるまで（概ね１０分程度）の経過時間と測定温度の関係を抽出しつつ、図８及び図１０のフローチャートに示されているような計算処理を行う。
【００４５】
図８は、予測式におけるＡの値の補正式を求めるための計算処理を示すフローチャートである。まず、図７に挙げた実験例において、計測時間ｔ＝０にセットして（ステップＳ２０１）、その後第１の所定時間ｔ１が経過したかどうかを判定する（ステップＳ２０２）。第１の所定時間ｔ１が経過するまではステップＳ２０２の処理を繰り返す。第１の所定時間ｔ１が経過すると、ｔ１における実測値のカーブに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂの式でカーブフィッティング（曲線どうしを沿わせる）したときのＡの値Ａ１を求める（ステップＳ２０３）。
【００４６】
続いて、第２の所定時間ｔ２が経過したかどうかを判定する（ステップＳ２０４）。第２の所定時間ｔ２が経過するまではステップＳ２０４の処理を繰り返す。第２の所定時間ｔ２が経過すると、ｔ２における実測値のカーブに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂの式でカーブフィッティング（曲線どうしを沿わせる）したときのＡの値Ａ２を求める（ステップＳ２０５）。そして、Ａ１、Ａ２の値をＡ１とＡ２の関係を示したグラフ平面上にプロットする（ステップＳ２０６）。
【００４７】
次に、これらの処理がすべての実験例について行われてプロットが終了したかどうかを判定する（ステップＳ２０７）。すべて終了していない場合（ＮＯ）は、ステップＳ２０１に戻って次の実験例についての処理を行い、すべて終了した場合（ＹＥＳ）は、ステップＳ２０８に進む。
【００４８】
Ａ１とＡ２との関係をプロットすると、図９のようなプロット結果を示した図となる。ステップＳ２０８では、これらの実験例のプロットの結果から、Ａ１とＡ２の関係として、２次式Ａ１＝ｐＡ２2＋ｑＡ２で補間できるようなｐ，ｑの値を求める。そして、後述するパラメータ保持部２１にＡ＝ｐｘ2＋ｑｘを保持する（ステップＳ２０９）。
【００４９】
図１０は、予測式におけるＢの値の補正式を求めるための計算処理を示すフローチャートである。まず、図７に挙げた実験例において、計測時間ｔ＝０にセットして（ステップＳ３０１）、その後第１の所定時間ｔ１が経過したかどうかを判定する（ステップＳ３０２）。第１の所定時間ｔ１が経過するまではステップＳ３０２の処理を繰り返す。第１の所定時間ｔ１が経過すると、ｔ１における実測値のカーブに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂの式でカーブフィッティング（曲線どうしを沿わせる）するようにＡの値を調整して、ｔ１における実測値の値を通るようにしたときのＢの値Ｂ１を求める（ステップＳ３０３）。
【００５０】
続いて、第２の所定時間ｔ２が経過したかどうかを判定する（ステップＳ３０４）。第２の所定時間ｔ２が経過するまではステップＳ３０４の処理を繰り返す。第２の所定時間
ｔ２が経過すると、ｔ２における実測値のカーブに、ｙ＝Ａｌｏｇ（ｘ）＋Ｂの式でカーブフィッティング（曲線どうしを沿わせる）するようにＡの値を調整して、ｔ２における実測値の値を通るようにしたときのＢの値Ｂ２’を求める（ステップＳ３０５）。次に、ステップＳ３０５で求めた式が、第１の所定時間ｔ１における実測値の値を通るように、Ｂの値を変更してオフセットさせたときのＢの値Ｂ２を求める（ステップＳ３０６）。そして、Ｂ２とＢ１との差分値Ｂ３を求める（ステップＳ３０７）。
【００５１】
次に、測定温度が平衡状態となっているかどうかを判定する（ステップＳ３０８）。平衡状態となるまでは、ステップＳ３０８の処理を繰り返す。平衡状態となった後に、平衡状態におけるｙ＝Ａｌｏｇ（ｘ）＋Ｂ２の値と実測値との差分Δｙを求める（ステップＳ３０９）。そして、Δｙの値とＢ３の値をΔｙとＢ３との関係を示したグラフ平面上にプロットする（ステップＳ３１０）。
【００５２】
次に、これらの処理がすべての実験例について行われてプロットが終了したかどうかを判定する（ステップＳ３１１）。すべて終了していない場合（ＮＯ）は、ステップＳ３０１に戻って次の実験例についての処理を行い、すべて終了した場合（ＹＥＳ）は、ステップＳ３１２に進む。
【００５３】
ΔｙとＢ３との関係をプロットすると、図１１のようなプロット結果を示した図となる。ステップＳ３１２では、これらの実験例のプロットの結果から、ΔｙとＢ３の関係として、１次式Δｙ＝ｒＢ３＋ｓで補間できるようなｒ，ｓの値を求める。そして、後述するパラメータ保持部２１にＢ＝ｒｘ＋ｓを保持する（ステップＳ３１３）。
【００５４】
図１４は、実験例から得られる数値の一覧表である。図１４において、列Ａは、５．５〜１０秒の間を０．５秒間隔で実測値をサンプリングして、それらのサンプリングした実測値がｙ＝Ａｌｏｇ（ｘ）＋Ｂの曲線にフィット（曲線どうしを沿わせる）するようにカーブフィッティングしたときのＡの値であり、列Ｂは、同様に６０〜１８０秒の間を０．５秒間隔で実測値をサンプリングして、それらのサンプリングした実測値がｙ＝Ａｌｏｇ（ｘ）＋Ｂの曲線にフィット（曲線どうしを沿わせる）するようにカーブフィッティングしたときのＡの値である。先ほどの、図８のフローチャートに基づいてＡの補正式を算出する際の値としては、これらの列Ａと列Ｂの値を用いている。
【００５５】
列Ｃは、５．５〜１０秒の間を０．５秒間隔で実測値をサンプリングして、それらのサンプリングした実測値がｙ＝Ａｌｏｇ（ｘ）＋Ｂの曲線にフィット（曲線どうしを沿わせる）するようにカーブフィッティングしたときのＢの値であり、列Ｄは、６０〜１８０秒の間を０．５秒間隔で実測値をサンプリングして、それらのサンプリングした実測値がｙ＝Ａｌｏｇ（ｘ）＋Ｂの曲線にフィット（曲線どうしを沿わせる）するようにした曲線が、１０秒時点での実測値を通るようにＢの値を変えてオフセットしたときのＢの値である。
【００５６】
そして、列Ｅは、列Ｃと列Ｄとの差の値である。また、列Ｆは、ｙ＝Ａｌｏｇ（ｘ）＋ＢのＡの値として列Ｂの値を、Ｂとして列Ｄの値を採用した式が１０分時点での予測温度として算出する算出値と、１０分時点における実測値との差の値である。列Ｆにおいて、＃ＤＩＶ／０！と表示されている欄は、１０分未満で実験を停止して、値の取得ができなかったことを示している。
【００５７】
図９は、図１４の各結果における列Ａと列Ｂとの値の関係をプロットした図である。これらのプロット結果から、図８のフローチャートのステップＳ２０８で示したように、両者の関係を２次式で求めると、ｙ＝０．０４１６ｘ2＋０．２８１６ｘが求められる。図８のフローチャートのステップＳ２０９で示したように、これらの関係式をパラメータ保
持部に保持するようにする。
【００５８】
同様に、図１１は図１４の各結果における列Ｆと列Ｇとの値の関係をプロットした図である。これらのプロット結果から、図１０のフローチャートのステップＳ３１２で示したように、両者の関係を１次式で求めると、ｙ＝０．２６５ｘ＋０．１９７１が求められる。図１０のフローチャートのステップＳ３１３で示したように、これらの関係式もパラメータ保持部に保持するようにする。
【００５９】
図１２は、図１に示された予測手段２の内部構成を示した図である。予測手段２の内部には、パラメータ保持部２１、カーブフィッティング部２２、及び熱平衡予測温度算出部２３を有し、実測値がカーブフィッティング部２２に入力されている。カーブフィッティング部２２においては、複数の実測値を用いてｙ＝Ａｌｏｇ（ｘ）＋Ｂの式へのカーブフィッティング（曲線どうしを沿わせる）を実行する。パラメータ保持部２１においては、各実験結果から算出されたＡ，Ｂの補正式を保持している。熱平衡予測温度算出部２３においては、後述する算出方法によって、パラメータ保持部２１において保持されているパラメータを用いて、カーブフィッティング部２２が算出したパラメータＡ，Ｂを補正する。
【００６０】
図１３は、事前の実験結果から得られた値を用いて、Ａ及びＢの補正式を算出した後に、実際に平衡予測温度を算出する方法の流れを示したフローチャートであり、図２のステップＳ１０７の予測値算出処理を示したものである。まず測定開始したかどうかを判定し（ステップＳ４０１）、測定開始したら、計測時間ｔ＝０にセットして（ステップＳ４０２）、実測値をサンプリングし、その結果を蓄積する（ステップＳ４０３）、その後第１の所定時間ｔ１が経過したかどうかを判定する（ステップＳ４０４）。第１の所定時間ｔ１が経過するまではステップＳ４０３の処理に戻る。
【００６１】
第１の所定時間ｔ１が経過したら、蓄積された複数の実測値を用いてｙ＝Ａｌｏｇ（ｘ）＋Ｂの式へのカーブフィッティング（曲線どうしを沿わせる）を実行する（ステップＳ４０５）。次に、パラメータ保持部に保持されている補正式を用いてパラメータＡ，Ｂを補正する（ステップＳ４０６）。そして、これらのＡ，Ｂの値を用いて、ｙ＝Ａｌｏｇ（ｘ）＋Ｂの式で平衡予測温度を算出する（ステップＳ４０７）。
【００６２】
図１５〜図２１は、本実施形態において求められたＡ及びＢの補正式を用いて、いくつかの実験例の結果から、実測値と予測値との関係を示した例である。それぞれ、図１５は図１４の表における実験２、結果２の値、図１６は実験２、結果３の値、図１７は実験２、結果４の値、図１８は実験２、結果５の値、図１９は実験２、結果６の値、図２０は実験２、結果７の値、図２１は実験２、結果８の値を用いたものである。すでに説明したとおり、それぞれの図における予測３の値が、パラメータＡ，Ｂの値を補正した最終的な予測値であるが、いずれの例においてもおおむね実測値に沿った変化を示しており、良好な予測値が得られることがわかる。
【００６３】
なお、予測式におけるパラメータＡ及びＢの補正式は、体温計の形状や材質によって異なるものとなることがある。そのため、それぞれの体温計に応じて事前の実験に基づいて、パラメータＡ及びＢの補正式を定めるのがよい。
【００６４】
また、本実施形態においては、人体の額で測る体温計を用いて説明してきたが、本発明は額で測る体温計に限られたものではなく、実験例で水の温度を測定したように、他の物体の温度を測定する温度計にも適用可能である。
【００６５】
また、本実施形態においては、第１の所定時間として１０秒時点の値、第２の所定時間
として３分時点の値を用いて計算を行ったが、それぞれ、１０秒、３分に限ったものではなく、適宜変更することも可能である。さらに、パラメータＡの補正式は２次式による補正、パラメータＢの補正式は１次式による補正を行っているが、これらについても他の形式の補正式を用いて補正式を構成するようにすることも可能である。
【産業上の利用可能性】
【００６６】
本発明の温度計は、予測手段における平衡温度の算出にあたって、事前に決定した補正パラメータを用いて算出を行うため、高速かつ簡便に平衡温度の算出を行うことが可能となる。
【符号の説明】
【００６７】
１    温度検出手段
１１  温度センサー
１２  温度測定部
２    予測手段
２１  パラメータ保持部
２２  カーブフィッティング部
２３  熱平衡予測温度算出部
３    表示手段
３１  表示切替部
３２  表示部
４    ブザー
５    前回値記憶手段
６    制御手段
【図１】
2015045602.tif 000003
【図２】
2015045602.tif 000004
【図３】
2015045602.tif 000005
【図４】
2015045602.tif 000006
【図５】
2015045602.tif 000007
【図６】
2015045602.tif 000008
【図７】
2015045602.tif 000009
【図８】
2015045602.tif 000010
【図９】
2015045602.tif 000011
【図１０】
2015045602.tif 000012
【図１１】
2015045602.tif 000013
【図１２】
2015045602.tif 000014
【図１３】
2015045602.tif 000015
【図１４】
2015045602.tif 000016
【図１５】
2015045602.tif 000017
【図１６】
2015045602.tif 000018
【図１７】
2015045602.tif 000019
【図１８】
2015045602.tif 000020
【図１９】
2015045602.tif 000021
【図２０】
2015045602.tif 000022
【図２１】
2015045602.tif 000023