KR910005476B1 - 생체의 온도측정방법 및 그 장치 - Google Patents

생체의 온도측정방법 및 그 장치 Download PDF

Info

Publication number
KR910005476B1
KR910005476B1 KR1019880701033A KR880701033A KR910005476B1 KR 910005476 B1 KR910005476 B1 KR 910005476B1 KR 1019880701033 A KR1019880701033 A KR 1019880701033A KR 880701033 A KR880701033 A KR 880701033A KR 910005476 B1 KR910005476 B1 KR 910005476B1
Authority
KR
South Korea
Prior art keywords
temperature
value
data
shape parameter
time
Prior art date
Application number
KR1019880701033A
Other languages
English (en)
Other versions
KR890700222A (ko
Inventor
게이지 야마구찌
Original Assignee
데루모 가부시끼가이샤
도자와 미쯔오
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 데루모 가부시끼가이샤, 도자와 미쯔오 filed Critical 데루모 가부시끼가이샤
Publication of KR890700222A publication Critical patent/KR890700222A/ko
Application granted granted Critical
Publication of KR910005476B1 publication Critical patent/KR910005476B1/ko

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K13/00Thermometers specially adapted for specific purposes
    • G01K13/20Clinical contact thermometers for use with humans or animals
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/42Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Measuring Temperature Or Quantity Of Heat (AREA)
  • Measuring And Recording Apparatus For Diagnosis (AREA)

Abstract

내용 없음.

Description

[발명의 명칭]
생체의 온도측정방법 및 그 장치
[도면의 간단한 설명]
제1도는 본 발명에 의한 실시예의 전자체온계의 기본 구성을 표시하는 블록도.
제2a도 및 제2b도는 제1실시예의 전자체온계의 구성을 표시하는 블록도.
제3a도∼e도는 제1실시예의 전자체온계에 의한 검온프로세스를 표시하는 플로우 챠트.
제4a도는 제1실시예의 형상인식처리를 표시하는 플로우챠트.
제4b도는 제1실시예의 곡선해석처리를 표시하는 플로우챠트.
제4c도는 제1실시예의 예측연산처리를 표시하는 플로우챠트.
제5도는 제1실시예의 전자체온계에 의한 검온 프로세스를 표시하는 타이밍챠트.
제6도는 실측하고, 또한 통계적처리를 하여 선택한 복수의 대표적인 온도상승커브를 표시하는 그래프도.
제7도는 체온측정계의 열전도모델의 일예를 표시하는 개념도.
제8도는 제1실시예의 전자체온계에 의한 겨드랑이 밑 검온의 평균적인 승온커브의 경우를 표시하는 그래프도.
제9도는 제1실시예의 전자체온계에 의한 겨드랑이 밑 검온의 체동(體動)에 의한 승온커브에 변동이 생겼을 경우를 표시하는 그래프.
제10도는 제1실시예의 전자체온계에 의한 겨드랑이 밑 검온의 승온커브가 극히 완만하게 상승하고 있는 경우를 표시하는 그래프도.
제11a도 및 제11b도는 제2실시예의 전자체온계의 구체적인 구성을 나타내는 블록도.
제12a도∼c도는 제2실시예에서의 형상인식처리, 곡선해석처리 및 예측연산처리를 표시하는 플로우챠트.
제13도는 제2실시예에서의 검온데이타 TO,피크데이타 TOP및 그 이동평균치 TOA의 추이관계를 표시하는 그래프도.
제14도는 제2실시예에서의 형상파라미터 a의 최적치와 그 형상파라미터 a를 결정하기 위한 변수(S1/S2)와의 관계를 표시하는 그래프.
제15a도 및 제15b도는 제3실시예의 전자체온계의 구체적인 구성을 나타내는 블록도.
제16a도∼e도는 제3실시예의 전자체온계에 의한 검온 프로세스를 표시하는 플로우챠트.
제17도는 제3실시예에서의 본 계측개시후의 각 소정시각에서의 복수의 검출온도데이타를 플롯트한 관계를 표시하는 그래프도.
제18도 및 제19도는 제2실시예의 전자체온계에 의한 겨드랑이 밑 검온의 측정경과를 표시하는 그래프도.
제20도 내지 제22도는 제3실시예의 전자체온계에 의한 구강검온의 측정경과를 표시하는 그래프도이다.
[발명의 상세한 설명]
[기술분야]
본 발명은 생체의 온도측정방법 및 그 장치에 관한 것이며, 특히 장래시각에 있어서의 검출온도예측기능을 갖춘 생체의 온도측정방법 및 그 장치에 관한 것이다.
[배경기술]
종래에는 이러한 종류의 온도측정장치(예를 들면 전자체온계)는 미리 온도상승커브를 완전히 규정한 예측식을 구비하고 있으며, 실체의 검출온도에 상기 예측식에서 구한 첨가량을 가산함으로써 평형온도의 선행표시를 하고 있었다. 이 때문에 예측식중의 각 정수(파라미터)는 가령 각 전자체온계의 제조공정에서 실제 측정에 사용하는 온도탐침에 의한 실측치의 통계처리등을 행함으로써 통계적으로 가장 예측오차가 적어지는 값으로 설정할 필요가 있었다.
그런데 온도상승커브에는 개인차가 있고, 또한 동일인인 경우에도 겨드랑이검온과 구강검온사이에는 상당한 상이가 있는 것으로 알려져 있다. 이와 같은 경우에는 탐침의 열특성불균일을 보정한 단일 예측식을 구비하고 있어도, 실제로는 정확한 평형온도의 선행표시가 불가하다.
특개소 58-225326호의 전자체온계는 복수의 예측식을 구비함으로써 이 문제를 해결하고 있다. 즉 대량의 측정결과에 의거하여 통계적 처리로 규정한 복수의 예측식을 미리 구비하여 두고, 측정시에는 피측정온도의 상승커브와 선택한 하나의 예측식을 비교하는 등으로 조건 설정(예측식의 파라미터)를 시행착오적으로 변경하여 문제를 해결하고 있다. 그러나 미리 복수의 예측식을 규정하여 두어야 하므로, 현실적인 문제로서 양산시에 생기는 온도탐침의 열특성불균일에 따르는 조정등의 번잡을 회피할 수 없다. 또, 예측정밀도를 높이기 위해서는 상승커브가 다른 다수의 예측식을 구비하지 않으면 안된다. 또한 온도상승부근에서 예측식의 선택이 적당하지 않을 때에는 예측치의 추이가 오버슈우트하기도 한다.
특개소 59-187233호의 전자체온계는 현실의 피측정온도의 상승커브에 부합하는 예측식을 조립함으로써 상기 문제를 해결하고 있다. 즉 피측정체온의 시간미분의 대수치 T과 샘플링시간 t와의 사이에는 직선적인 관계(TL=A-τ't)가 있는 것에 착안하여, 정수 A 및 τ'를 회귀법에서 구하고 있다. 그러나 대수치 T은 피측정체온데이타 그 자체는 아니므로, 체온데이타에 미분연산 및 대수연산을 행함으로써 오차가 야기되며, 그 오차는 정수 A 및 τ'의 설정에 큰 비율로 영향을 미친다. 더욱이 피측정체온데이타가 잡음성분을 내포하고 있을때에는 결과적으로 잡음성분이 예측결과에 지수함수적으로 영향을 주게 되어, 대단히 불안정한 예측추이를 초래하였다.
[발명의 개시]
본 발명은 상술한 종래기술의 결점을 제거하는 것이며, 그 목적하는 바는 탐침의 열특성불균일 또는 개인이나 피측정부위의 차이에 의한 온도상승커브에 차이가 있어도 항상 정확한 선행표시가 가능한 생체의 온도 측정방법 및 그 장치를 제공하는데 있다.
본 발명의 타목적은 검출온도가 잡음성분을 내포하고 있어도 안정된 예측추이를 얻을 수 있는 생체의 온도측정방법 및 그 장치를 제공하는데 있다.
본 발명의 또 하나의 목적은 장래의 임의의 시각에 있어서의 검온치를 용이하고도 정확하게 예측하는 생체의 온도측정방법 및 그 장치를 제공하는데 있다.
본 발명의 또다른 목적은 장래의 매우 긴 시간을 경과한 후의 열평형온도차를 정확히 예측하는 생체의 온도측정방법과 그 장치를 제공하는데 있다.
본 발명의 또 하나의 목적은 생체의 운동등에 의하여 발생하는 온도변화의 유동의 영향을 경감시킴으로써 예측온도의 신뢰성을 한층 향상시킨 생체의 온도측정방법 및 그 장치를 제공하는데 있다.
본 발명은 타목적은 측정대상을 생체의 특정부위, 예를 들면 인체의 겨드랑이, 구강, 직장등 부위에 한정함으로써 예측정밀도를 한층 향상시킨 생체의 온도측정방법 및 그 장치를 제공하는데 있다.
본 발명의 타목적은 측정의 조기단계에서 검온커브의 형상 또는 그 상승형상을 정확히 인식함으로써 비교적 빠른 단계에서 유효한 예측표시가 가능한 생체의 온도측정방법 및 그 장치를 제공하는데 있다.
본 발명의 타목적은 다수의 온도데이타를 축적 혹은 동시에 처리하는 일이 없이 간단한 구성과 데이터처리방법으로 상기 목적을 달성하는 생체의 온도측정방법 및 그 장치를 제공하는데 있다.
본 발명에 의한 생체의 온도측정방법은 상기의 목적을 달성하기 위하여 검온커브의 형상을 반영하기 위한 형상파라미터 및 상기 검온커브에 중복시키기 위한 계수파라미터의 각 수치를 미정으로 하는 소정의 예측함수식을 갖추는 공정과, 온도를 검출하여 그 온도를 표시하는 온도데이타를 발생하는 온도검출공정과, 측정개시후의 경과시간을 계시하여 그 경과시간을 표시하는 시간데이타를 발생하는 시간검출공정과, 상기 검출된 소정의 온도데이타에 의거해서 상기형상파라미터의 값을 설정하는 형상파라미터의 설정공정과, 상기 설정된 형상파라미터를 포함하며, 또한 상이한 복수 시점의 온도데이터를 목적변수로 하고, 또한 상기 복수시점의 시간데이타함수를 설명변수로 하는 복수의 상기 예측함수식으로서 구성되는 연립방정식을 풀므로써 상기 계수파라미터의 값을 설정하는 계수파라미터설정공정과, 상기 설정된 형상파라미터 및 계수파라미터의 각 값으로 특정하는 상기 예측함수식에 의해 장래시각에서의 검출온도를 예측연산하는 예측연산공정등을 갖추는 것을 그 개요로 한다.
또, 바람직하게는 형상파리미터설정공정은 복수의 온도데이타에 의거하는 소정 온도상승구배정보에 근거하여 형상파라미터치를 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 형상파라미터설정공정은 우선 피검온커브가 제1의 소정구배를 표시하는 점을 검출하고, 또한 그 점에 관해서 전후하는 제2의 구배 및 제3의 구배를 검출하여 이들을 비교함으로서 형상파라미터의 값을 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 형상파라미터설정공정은 제2의 구배 S1,및 제3의 구배 S2에 의거해서 형상파라미터 a의 값을 다음식,
Figure kpo00001
단,
n〈1 : 정수
∼Q: 정수
에 따라서 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 형상파라미터설정공정은 측정개시후의 초기단계에서의 복수의 온도데이타에 의거해서 형상파라미터의 값을 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 형상파라미터 설정공정은 각 소정시각에서 복수온도데이타 T에 의거하는 정보 X에 따라서 형상파라미터 a의 값을 다음식,
Figure kpo00002
단,
0∼D5: 정수
0∼X3: T0∼T3
4=(X3-X0)/(X1-X0
에 따라서 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 계수파라미터 설정공정은 상이하는 2점에서의 온도 데이터 T(t1),T(t2) 및 각 검온시점의 시간데이타 t,t에 의거해서 다음 식의 2원연립방정식,
Figure kpo00003
Figure kpo00004
를 풀므로써 계수파라미터 A및 A의 각 수치를 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 계수파라미터 설정공정은 상이하는 2점의 온도데이타로서 측정개시점 부근의 온도데이타와 현시점의 온도데이타를 사용함을 그 하나의 양태로 한다.
또, 바람직하게는 예측연산공정은 형상파라미터 a 및 계수파라미터 A0및 A1의 각 수치로 특정하는 예측함수식에 의하여 장래시각 t에서의 검출온도 T(t)를 다음식,
Figure kpo00005
에 따라서 에측연산함을 그 하나의 양태로 한다.
또, 본 발명에 의한 생체의 온도측정장치는 상기의 목적을 달성하기 위하여 검온커브의 형상을 반영하기 위한 형상파라미터 및 상기 검온커브에 중복시키기 위한 계수파라미터의 각 수치를 미정으로 하는 소정의 예측함수식을 기억하고 있는 기억수단과, 온도를 검출하여 그 온도를 표시하는 온도 데이터를 발생하는 온도 검출수단과, 측정개시후의 경사시간을 계시하여 그 경과시간을 표시하는 시간데이타를 발생하는 시간검출수단과, 상기 검출한 소정의 온도 데이터에 의거해서 상기 형상파라미터의 값을 설정하는 형상파라미터 설정수단과, 상기 설정한 형상파라미터의 값을 포함하며, 또한 상이하는 복수시점의 시간데이타의 함수를 설명변수로 하는 복수의 상기 예측함수식으로 구성되는 연립방정식을 풀므로써 상기 계수파라미터의 값을 설정하는 계수파라미터 설정수단과, 상기 설정한 형상파라미터 및 계수파라미터의 각 수치로 특정하는 상기 예측함수식에 의하여 장래 시각에서의 검출온도를 예측연산하는 예측연산수단을 갖추는 것을 그 개요로 한다.
또, 바람직하게는 온도 검출수단은 검출한 온도중의 최대레벨의 것을 순차검출하여 유지출력하는 피크홀드 수단을 갖추는 것을 그 하나의 양태로 한다.
또, 바람직하게는 온도검출수단은 소정온도로 검출한 온도중의 최대레벨의 것을 순차검출하여 유지하는 피크홀드수단과, 상기 피크홀드 수단이 유지한 복수의 레벨에 관한 이동평균치를 구하여 출력하는 평균화수단을 갖추는 것을 그 하나의 양태로 한다.
또, 바람직하게는 형상파라미터설정수단은 복수의 온도데이타에 의거하는 소정의 온도상승구배정보에 따라서 형상파라미터의 값을 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 형상파라미터설정수단은 측정개시후의 초기단계에서의 복수의 온도데이타에 의거해서 형상파라미터의 값을 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 계수파라미터설정수단은 상이하는 2점에서의 온도 데이터 T(t1),T(t2) 및 각 검온시점의 시간데이타 t1,t2에 의거해서 다음식의 2원연립방정식,
Figure kpo00006
Figure kpo00007
를 풀므로써 계수파라미터 A및 A의 각 수치를 설정함을 그 하나의 양태로 한다.
또, 바람직하게는 계수파라미터설정수단은 상이하는 2점의 온도데이타로서 측정개시점 부근의 온도데이타와 현시점의 온도 데이터를 사용함을 그 하나의 양태로 한다.
또, 바람직하게는 예측연산수단은 형상파라미터 a 및 계수파라미터 A0및 A1의 각 수치로 특정하는 예측함수식에 의하여 장래시각 t에서의 검출온도 T(t)를 다음식,
Figure kpo00008
에 따라서 예측연산함을 그 하나의 양태로 한다.
[발명의 실시를 위한 최선의 형태]
이하 첨부도면을 참조하여 본 발명의 실시예를 상세히 설명한다.
[제1실시예]
[구성]
제1도는 본 발명에 의한 실시예의 전자체온계의 기본구성을 표시하는 블록도이다. 이 전자체온계는 기본적으로는 온도 측정부(1)와 예측연산부(2)와 표시부(3)로써 구성된다.
온도측정부(1)는 피측정부위의 온도를 소정주기로 검출하여 그 온도를 표시하는 온도데이타 T를 라인(103)에 출력하는 부분이다.
또, 예측연산부(2)는 미리 검온커브의 형상을 반영하기 위한 형상파라미터 및 상기 검온커브에 중복시키기 위한 계수파라미터의 각 수치를 미정으로 하는 소정의 예측함수식을 갖추고 있으며, 측정개시전은 온도측정부(1)로부터의 검온데이타 T를 감시하여 소정의 측정개시조건을 판정하고 측정개시후는 온도측정부(1)로부터의 검온데이터 T 및 내장하는 측정경과시간 계측기능으로부터의 시간데이타를 감시함으로써, 우선 검온데이타가 제1의 구배를 나타내는 점을 검출하고, 또한 그 점에 대해서 전후하는 제2 및 제3의 구배를 검출하여 이들을 비교함으로써 상기 형상파라미터의 값을 설정하고, 그리고 상기 설정된 형상파라미터의 값을 포함하며, 또한 2개점의 검온데이타를 목적변수로 함과 동시에, 그 검출시점의 시간데이타의 함수를 설명변수로 하는 2개의 상기 예측함수식으로 구성되는 2원연립방정식을 풀므로써 상기 계수파라미터의 값을 설정하고, 다시 상기 설정된 각 파라미터의 값으로 특정하는 상기 예측함수식에 의하여 장래시각 t에서의 검출온도를 예측연산하며, 연산결과의 예측온도 T(t)를 라인(121)에 출력하는 부분이다.
또 표시부(3)는 측정경과와 더불어 축차적으로 구한 예측온도 T(t)를 수치표시하는 부분이다.
제2a도 및 제2b도는 제1실시예의 전자체온계의 구체적인 구성을 표시하는 블록도이다.
도면에서, 온도측정부(1)는 써어미스터등의 감온소자(4) 및 온도계측회로(5)를 포함하며, 그 온도계측회로(5)는 예측연산부(2)가 라인 (102)을 통하여 보내는 소정주기의 데이터샘플링신호 C1에 따라서 감온소자(4)의 검출온도에 따른 어낼로그전기신호(101)를 샘플링하여 이를 디지털의 온도데이타 T로 변환하고 라인(103 및 104)에 출력하는 부분이다.
예측연산부(2)는 데이터판독기입수단(6), 시간측정수단(7), 계측제어수단(8), 온도데이타의 기억수단(9), 검온커브의 형상인식수단(10), 검온커브의 곡선해석수단(11), 예측연산수단(12) 및 데이터선택부(13)를 포함하며, 이들 각 블록구성은 예컨대 마이크로컴퓨터(CPU)가 도시하지 않는 메모리(ROM 또는 RAM)에 격납되어 있는 제3a도∼e도 및 제4a도∼ c도의 프로그램을 실행함으로써 실현가능하다. 또 예측연산부(2)는 사용자에게 예측유효를 알리기 위한 부저(14)를 포함한다.
계측제어수단(8)은 본 전자체온계의 전체의 동작을 총괄제어하는 수단이며, 측정개시전에는 전력절약을 위하여 예로서 5초에 1회의 비율로 온도계측회로(5)에 온도 데이터 T를 발생시켜, 라인(104)을 통해서 이 온도데이타를 상시 모니터하고, 소정의 측정개시조건을 만족시켰다고 판정하였을 때, 즉 온도 데이터 T가 어떤 온도이상을 표시하며, 또한 어떤 값 이상의 온도변화를 수반하는 것이라고 판정하였을 때에는 라인(105)을 통해서 제어신호 C2를 보내고, 데이터판독 기입수단(6), 시간측정수단(7) 및 기억수단(9)등의 각 기능블록을 작동시켜서 측정을 개시시킨다. 또 측정개시후는 라인(100)을 통해서 CPU가 발생하는 1초주기의 클록신호 CLOCK을 받아, 이에 응동하도록 설치된 후술하는 타이머개입의 처리에 있어서 이하와 같은 각 블록 구성의 동작을 진행시킨다.
또, 계측제어수단(8)은 계측제어의 진행에 필요한 판단정수를 기억하고 있는 각종 레지스터군을 갖추고 있으며, 이 중 레지스터 t은 실제의 체온측정상에서는 이미 예측연산을 계속하는 것이 무의미하게 하는 등의 측정경과시간치(예를 들면 600초)를 기억하고 있고, 레지스터 t는 예측함수의 파라미터설정에 필요한 온도 데이터의 하나를 샘플링하는 시간으로서 측정개시시점부근의 시간(예로서 1초)을 기억하고 있고, 레지스터K는 검온커브가 제1의 구배를 표시하는 점을 검출하기 위한 온도 증가분(예를 들면 8초간에 0.15℃의 증가분)을 기억하고 있고, 레지스터 q는 소정시간간격(예를 들면 8초간격)으로 연산한 각 예측온도의 차의 절대치에 대해서 예측유효로 간주하기 위한 허용치(예를 들면 0.02℃)를 기억하고 있고, 레지스터 DT는 제1의 구배검출시점으로부터 제3의 구배를 검출할때까지의 시간폭정보(예를 들면 16초)를 기억하고 있다.
데이터 판독기입수단(6)은 계측제어수단(8)이 라인(102)을 통해서 보내는 1초 주기의 샘플링신호 C1과 더불어 라인(103)상에 출력된 검온데이타 T를 예측연산부(2)내에 판독기입하는 수단이다.
또, 데이터 판독기입수단(6)은 예컨대 복수의 연속된 검온데이터 T를 신구갱신하면서 (FIFO와 같이 하면서)기억할 수 있는 것으로, 그 출력단자에는 상기 복수의 검온데이타 T로부터 이동평균치 Tav를 출력하는 것도 무방하다. 이와같이 하면 각 검온데이타 T가 평균화되고, 그 결과 예측온도의 추이가 안정화하는 효과가 있다.
시간측정수단(7)은 측정개시후의 경과시간을 계시하여 그 측정 경과시간을 표시하는 측정경과시간데이터t1를 라인(107)에 출력하는 수단이며, 측정개시후에 계측제어수단(8)이 라인(105)을 통해서 보내는 1초주기의 신호를 카운트업함으로써 측정개시시점으로부터의 측정경과시간을 계시한다.
기억수단(9)은 데이터수단(6)이 판독기입한 검온데이타 T0를 순차로 기억수단의 레지스터 T8로부터 T0를 향해서 시프트인/시프트아웃트하면서, 항상 최신의 것으로부터 합계 9개분의 검온데이타를 기억하고 있는 수단이다.
형상인식수단(10)은 체온측정의 초기단계에서 관측되는 검온커브가 어떠한 상승형상으로서 파악되는 가를 간단하고도 유효한 방법으로 인식하는 수단이며, 이 수단의 목적은 체온측정의 초기단계에서 예측함수식의 형상파라미터 a의 값을 설정하는데 있다.
이하, 형상파라미터 a의 값의 설정개념을 제6도를 참조하면서 설명한다.
제6도는 실측하고, 또 통계적 처리를 하여 선택한 복수의 대표적인 온도상승커브를 표시하는 그래프도이다. 도면에서 그래프의 횡축은 측정경과시간 t(초)를 표시하며, 종축은 검온데이타 T(℃)를 표시한다. 이 그래프도에서 명백한 바와 같이, 이들 온도상승커브는 모두 측정개시점부근에서 상당히 급격한 상승특성을 나타내고 있으나, 그 후는 제1의 구배(도면중에서는 측정경과시간 tB1∼tB6까지 사이에 발생하고 있다)를 표시하며, 또한 그 후는 완만한 상승특성을 나타내고 있다.
그런데 이와 같은 온도상승커브의 형상은 후술하는 본 발명의 예측식, 즉,
Figure kpo00009
의 우측변 제2항의 요소(1/t )로서 표시할 수가 있다. 즉, 형상파라미터 a를 적당하게 선택하면 온도상승 커브중 어느 하나의 형상을 표현할 수가 있다. 그리고 도면에는 한 예로서, 각 온도상승커브의 형상에 대응하는 형상파라미터 a의 값을 표기하고 있다. 이에 의하면 측정개시부근에서는 a의 값이 클수록 온도상승커브는 급격한 상승특성을 표시하며, a의 값이 작을수록 온도상승커브는 완만한 상승특성을 표시하고 있음을 알 수 있다. 그러나 측정개시부근의 상승특성이 급격한가 완만한가 하는 것만으로는 복수의 온도상승커브를 그렇게 쉽게 분별할 수 없다.
그래서 형상인식수단(10)은 우선 온도상승커브가 제1의 상승구배를 표시하는 점(도면의 어깨부분)을 검출하고, 다음에 이 점을 경계로 하여, 이에 의해서 앞서의 상대적으로 급격한 제2의 상승구배와, 그 보다도 뒤의 상대적으로 완만한 제3의 상승구배를 검출하여 이들을 비교함으로써, 복수의 온도상승커브를 분별하고 있다. 이 비교는 예로서 도면의 a=0.6의 온도상승커브에서 말하면, 제2의 상승구배는 매우 급격하며, 제3의 상승구배는 매우 완만하다. 이에 비하여 a=0.1의 온도상승커브에서는 제2의 상승구배는 비교적 완만하며, 제3의 상승구배도 비교적 완만하다. 따라서 어깨부분에 전후하는 제2 및 제3의 상승구배를 조사하면, 복수의 온도상승커브를 분별하기 위한 보다 유용한 정보가 얻어진다.
그래서 양자의 비(제2의 상승구배)/(제3의 상승구배)를 채용함으로써 형상상이가 한층 더 명확해지고, 또한 그 비의 값은 형상파라미터 a의 값에 비례하는 것이 기대된다.
또, 형상인식수단(10)에 있어서, 형상파라미터 a의 값의 설정방법은 그외의 것도 생략할 수가 있다. 더욱 간단한 방법으로서 예컨대 소정시각에 있어서의 온도의 구배나, 측정개시후의 소정시각에서의 복수의 온도데이타등에 의해서 형상파라미터 a를 표현하는 것도 가능하다. 요컨대 형상파라미터 a에 대해서 상관성이 높은 변수 혹은 변수의 조에 의해서 형상파라미터 a의 값을 설정하기 위한 수단이 형상인식수단(10)이다.
이상과 같은 개념에 의거하여, 실시예의 형상인식수단(10)은 체온측정의 초기단계에서의 검온커브에 근거해서 형상파라미터 a의 값을 설정한다.
또, 형상인식수단(10)은 형상파라미터 a의 값의 설정에 필요한 정보를 기억하는 레지스터군을 갖추고 있으며, 이 중에 레지스터 T는 측정개시후의 비교적 빠른 시점 t(예컨대 1초경과시점)에서의 검온데이터 T를 기억하고, 레지스터 T는 계측제어수단(8)이 제1의 구배(어깨부)를 검출한 시점의 검온데이타 T를 기억하고, 레지스터 S1은 검온데이타의 차분수치 S1=T-T를 기억한다. 또한 레지스터 S2는 어깨부의 검출로부터 소정시간 DT(예를 들면 16초)경과후의 검온데이터를 T0로 할 때에, 검온데이터의 차분수치 S2=T8-T를 기억하고, 레지스터 a는 다음식,
Figure kpo00010
의 연산에 의해서 구한 형상파라미터 a의 값을 기억한다. 또, 이때에 레지스터 Q1은 정수 Q1(예컨대 0.024)을 기억하고 있고, 레지스터 Q2는 정수 Q2(예컨대 -0.128)을 기억하고 있다.
곡선해석수단(11)은 상기 형상인식수단(10)이 설정한 형상파라미터 a의 값을 포함하는 예측함수식을 사용하여 그 예측함수식의 계수파라미터 A0및 A1의 각 수치를 설정하는 수단이며, 구체적으로는 측정개시점 부근의 온도데이터 T0(t1)와 현시점의 온도 데이터 T0(t2) 및 각 대응하는 검출시점의 시간데이터 t1,t2에 의거해서, 다음식의 2차연립방정식,
Figure kpo00011
Figure kpo00012
을 풀므로서, 예측함수식의 계수파라미터 A0및 A1의 각 수치를 설정하는 수단이다. 이리하여 예측함수식의 계수파라미터 A0및 A1의 각 수치가 설정되어, 이미 설정되어 있는 형상파라미터 a의 값과 함께 예측식이 일시적으로 특정된다.
또, 이 곡선해석수단(11)은 수시로 특정한 예측식의 유효성을 조사하기 위하여, 현시점의 온도데이터 T0(t2)로서, 예컨대 현시점으로부터 8초전의 검온데이터 T0를 사용함으로써 8초전의 계수파라미터 A0'및 A1'의 각 수치를 설정하기도 한다. 이 때문에 곡선해석수단(11)은 현 시점의 검온데이터 T0(t2)로서 현 시점의 검온데이터 T을 사용해서 구한 계수파라미터 A0및 A1의 각 수치를 기억하는 레지스터 A0,A1과, 현 시점의 검온데이터 T(t)로서 현 시점보다 8초전의 검온데이터 T0를 사용해서 구한 계수파라미터 A0',A1'의 각 수치를 기억하는 레지스터 A0',A1'를 갖춘다.
예측연산수단(12)은 곡선해석수단(11)이 특정한 현시점의 예측식을 사용하여 장래의 바람직하게는 임의 시각에서의 검온수치를 예측연산하는 수단이며, 구체적으로는 설정된 형상파라미터 a 및 계수파라미터 A0및 A1의 각 수치로 특정하는 예비함수식에 의하여, 장래시각 t0에서의 검출온도의 예측치 T(t)를 다음식,
Figure kpo00013
에 따라서 예측연산하고, 라인(120)에 예측온도데이터 T(t)를 출력한다.
또, 예측연산수단(12)은 8초전의 계수파라미터 A',A'의 각 수치로 특정하는 예측함수식에 의하여 마찬가지로 장래시각 t0에서의 검출온도의 예측식 T(t)'를 다음식,
Figure kpo00014
에 따라서 예측연산하고, 라인(119)을 통해서 그 예측온도데이터 T(t)'를 상기 예측온도데이터 T(t)와 함께 계측제어수단(8)으로 보낸다. 이에 의하여 계측제어수단(8)에서는 2개의 예측온도데이터 T(t),T(t)'를 비교하여 예측치의 유효성(일정성)을 판단한다.
또, 예측연산수단(12)은 예측온도의 제공에 필요한 정보를 기억하는 레지스터군을 가지고 있으며, 이 중 레지스터 T는 현시점에서 예측한 예측온도 T를 기억하고, 레지스터 T'를 기억하고, 레지스터 t는 장래시각 t(예컨대 10분경과시)를 기억하고 있다.
데이터선택부(13)는 통상적인 체온측정에서, 예측연산을 계속하는 것이 무의미할 정도로 충분히 시간이 경과한 후에는, 예측온도 T(t)에 의한 선행표시를 중지하고 그대신 검온데이터 T0에 의한 직시표시로 변환하는 수단이며, 그 데이터선택부(13)는 측정개시후의 소정의 예측정지조건을 판정할때까지는 예측연산수단(12)측에 접속하고 있으므로, 표시부(3)는 예측온도 T(t)를 표시하고 있다.
[예측원리]
그와 같은 구성에 있어서, 본 발명에 의한 온도예측의 동작원리는 다음과 같다.
본 발명자는 체온측정계에서의 열전도이론해석을 함으로써, 체온측정시의 온도탐침의 온도상승커브가 어떠한 형상이 되는가를 예측하여 보았다. 즉, 본 해석에서는 예로서 제7도와 같은 체온측정계의 모델을 사용하고, 그 측정계를 탐침, 피부, 피하조직의 3개 영역으로 나누고, 각 영역의 온도분포를 체온측정의 과정에서 같다고 가정하고 있다. 즉, 각 영역을 미소체적의 개념으로 취급하게 되는데, 피하조직에 대해서는 그 열용량을 무한대로 가정한다. 단, 피부, 피하조직의 명칭은 인체측을 2층 모델로 가정하기 위하여 편의적으로 붙인 것으로, 실제의 인체구조와 엄밀하게는 대응하지 않는다. 또, 필요하다면 금후의 전개에 따라서 영역을 다분할함으로써 보다 좋은 근사모델로 개량할 수도 있다.
제7도의 측정모델에 있어서, 탐침과 피부간의 열전달률을 h1, 그의 경계면적을 A1, 및 피부와 피부조직간의 열전달률을 h2, 그의 경계면적은 A2로 한다. 또 피하조직의 온도는 열용량이 무한대라고 하는 가정에 의해, 시간에 대하여 일정치 Tsat를 취하는 것으로 한다. 이렇게 해서 탐침을 피부측정부위에 장착한 후에 탐침이 피부로부터 흡수하는 열량은 탐침의 내부에너지의 증가분과 같으므로 다음의 (1)식이 성립한다.
Figure kpo00015
동일하게 해서 피부가 피하조직 및 탐침으로부터 흡수하는 열량은 피부의 내부에너지의 증가분과 같으므로 다음의 (2)식이 성립한다.
Figure kpo00016
여기서,
T
Figure kpo00017
P, C,V: 탐침의 온도, 밀도, 비열, 체적
Figure kpo00018
,C,V: 피부의 온도, 밀도, 비열, 체적
sat[=T(∞)] : 피하조직의 온도=평형온도이다.
다음에 상기(1)식 및 (2)식으로 이루어진 연립선형미분방정식을 풀면, 다음의 (3)식을 얻을 수 있다.
Figure kpo00019
그러나,
Figure kpo00020
이다.
상기(3)식은 고계선형미분방정식이므로 라플라스변환을 사용하여 풀수가 있다. 즉,
Figure kpo00021
라고 각 항을 계산하면,
Figure kpo00022
로 되고, 여기서,
Figure kpo00023
이다.
다시 이것을 T에 대하여 풀면,
Figure kpo00024
가 얻어진다.
여기서 s2+(K1-K2+K3)s+K1·K3=O의 해를 m1,m2라고 해두면,
Figure kpo00025
로 된다.
지금, m1
Figure kpo00026
2일때에는 (4)식으로부터,
Figure kpo00027
로 된다.
여기서 £`1kx=1/(s-k)가 알려져 있으므로, T(t)의 식은 다음과 같이 얻어진다.
Figure kpo00028
그러나,
Figure kpo00029
이다. 또, m1=m2일때는 (4)식으로부터,
Figure kpo00030
이다.
또 여기서 £ekx=1/(s-k), £ekx=1/(s-k)2가 알려져 있으므로 T(t)의 식은 다음과 같이 얻어진다.
Figure kpo00031
그러나,
Figure kpo00032
이다.
이렇게해서 탐침의 온도상승곡석의 이론식은 (5)식 및 (6)식과 같이 주어진다. 그런데 상기한 (5)식 및(6)식에 있어서, m1,m2및 M1∼M4는 탐침 및 피부의 각 물성치(밀도, 비열, 체적등)을 비롯한 체온측정계에 포함되어 있는 각종 물리량의 함수로서 주어져 있으며 이들의 값은 체온계의 매 측정마다 변화한다. 따라서 측정시에는 탐침으로 검출되는 온도데이타를 기초로 하여 m1,m2및 M1∼M4를 결정할 필요가 있다.
또, 전자체온계중에는 탐침을 측정부위에 장착한 후, 어떤 방법, 예컨대 탐침이 소정온도를 검지하기까지는 온도데이터의 판독을 개시하지 않는 것도 있으므로, 이와 같은 전자체온계에 대해서는 다시(5)식 및 (6)식을 이하 식과 같이 변형하면 편리하다.
Figure kpo00033
Figure kpo00034
그러나,
Figure kpo00035
이며, 여기서,
Δt : 탐침장착후, 측정개시까지의 시간
t : 측정개시시각을 t=0으로 한 시간
이다.
상기 (7)식에 있어서는 가령 m1,m2를 고정치로 할 수 있으면, 측정할때에 시계열적으로 검출되는 온도데이터를 사용하여 회귀분석법 또는 연립방정식을 풀므로써, Tsat,P,Q를 비교적 용이하게 구할 수가 있다. 그러나, m1,m2는 피측정자 또는 측정조건의 상위에 따라 측정마다 변화하는 값이며, 더욱이 측정마다 이와 같은 모든 가변요소를 도입하여 최적의 예측함수를 발견하고, 이로써 보편성이 높은 온도예측을 행하는 것이 본 발명의 목적이다. 이 경우에 측정할때에 검출되는 온도데이터를 사용하여 상기(7)식의 연립방정식을 푸는 것에 의하여 m1,m2및 Tsat,P,Q를 구하는 방법은 수학적으로는 있으나 검출한 온도데이터에는 잡음성분이 포함되어 있다는 것과(7)식중에는 지수부를 포함하고 있다는 상승효과에 의하여 매우 불안정한 결과를 준다.
여기서 상기(7)식을 데일러전개하여 다음의 (9)식을 얻는다.
Figure kpo00036
그리고 본 실시예에서는 예컨대 4차이후에는 생략하고 다음의 (10)식을 얻었다.
Figure kpo00037
또한 이상의 것은 상기 (8)식에 관해서도 동일하다.
그런데 본 발명자는 상기(10)식을 사용하여, 또한 이산적으로 추출한 4개의 온도데이터 T(t) 및 그 시점의 데이터 t에 의거하여 다음식의 4원연립방정식,
Figure kpo00038
을 풀므로써, 예측함수의 계수퍼레미터 A0∼A3을 설정하는 방식의 전자체온계를 이미 제안하고 있다. 이와 같이 하면, 선행하는 온도데이타에는 모든 물리적 조건이 포함되어 있으므로, 이들 온도데이터 및 시간데이터의 상관관계로부터 예측함수의 계수파라미터의 값을 설정할 수 있고, 당면한 최적의 예측함수를 특정지울 수 있다. 그리고 이 특정된 예측함수를 사용하여 장래시각 t에서의 온도를 다음식,
Figure kpo00039
에 따라서 예측연산할 수가 있으며, 이 방법은 대단히 보편성이 높고, 안정된 예측결과를 제공한다.
그러나 상기의 경우는 사용하는 온도데이터는 4개이지만, 매회의 당면 최적예측함수를 특정하기 위해서는 항상 최신의 검온커브의 풀스케일을 커버하도록 4개의 온도데이터를 이산적으로 추출할 필요가 있으므로, 이 때문에 오래된 검온데이터를 버릴 수가 없어, 결국 방대한 용량의 검온데이터메모리를 필요로 하였다.
더욱이, 이 점은(10)식의 우측변의 항수를 적게하면 다소나마 메모리 용량을 감소시킬 수 있는 기대를 할 수 있으나, 지나치게 항수를 적게 하면 예측함수의 상승특성도 완만하여져서 유효한 선행표시가 불가하다.
그래서 본 발명에서는 메모리용량의 감소와 유효한 선행표시가 불가하다.
그래서 본 발명에서는 메모리용량의 감소와 유효한 선행표시등을 일거에 달성하기 위해서, 상기 (10)식의 예측함수를 변형하여 다음의 (11)식을 채용하기로 하였다. 즉,
Figure kpo00040
이다.
[동작]
제3a도∼e도 및 제4a도∼c도는 제1실시예의 전자체온계에 의한 검온프로세스를 나타내는 플로우챠트이다. 제5도는 제1실시예의 전자체온계에 의한 검온프로세스를 표시하는 타이밍챠트이다. 제3a도에 있어서, 전원을 투입하면, 개시공정(S100)에 입력하고, 먼저 온도측정부(1) 및 계측제어수단(8)이 작동하여 비교적 거치른 온도계측공정(S101)에 입력한다. 온도계측공정(S101)에서는 계측제어수단(8)이 전력의 절약을 위해서 예컨대 5초에 1회의 비율로 온도계측회로(5)에 온도를 검출시켜, 이 검온데이터(T)를 모니터한다. 판단공정(S102 및 S103)은 예측방식에 의한 체온측정을 개시할 것인가 아닌가를 판단하는 부분이며, 공정(S102)에서는 소정온도(T), 예컨대 30℃를 초과하고 있는지 아닌지를 판단하고, 공정(S103)에서는 예컨대 1초간에 환산하여 0.1℃이상의 온도상승이 인정되느냐 아니냐의 판단을 하고 있다. 이 판단은 실제상으로는 체온계를 겨드랑이 밑 또는 구강내와 같은 피측정부위에 삽입하였는가 여부의 판단이 된다. 그리고 상기 어느 조건도 충족되면, 공정(S104)로 나아가고, 계측제어를 위한 각종 프로그램스위치SW1∼SW4를 클리어 한다. 공정(S105)에서는 라인(105)을 통하여 시간측정수단(7)을 클리어스타트한다. 즉, 시간측정수단(7)의 시간측정용 카운터를 리세트(제5도의 측정경과시간 t0에 대응)함과 동시에, 측정경과시간의 계측을 개시시킨다. 또한 공정(S106)에서는 라인(105)을 통하여 데이터판독기입수단(6)의 데이터읽어 넣는 기능을 작동시키고, 공정(S107)에서는 기억수단(9)의 데이터기억기능을 작동시킨다. 그리고 이상의 초기설정제어를 실행한 후에는 공정(S108)을 1초에 1회의 비율의 타이머중단기능을 작용시키고, 공정(S109)에서는 CPU가 아이돌루틴(idle routine)을 실행하여 타이머중단의 발생을 기다린다.
제3a도의 공정(S109)에서 타이머중단이 발생하면, 제3b도의 공정(S200)에 입력한다. 공정(S201)에서는 후의 일련의 처리를 실행하기 위하여 타이머중단기능을 소멸시킨다. 공정(S202)에서는 계측제어수단(8)이 t>t(가령 600초)인가 여부를 조사한다. 측정개시시점에서는 t>t을 만족시키지 못하므로 공정(S203)으로 진행하여 SW=1인가 여부를 조사한다. SW1은 측정개시부근의 검온데이터를 기억하는 처리를 위한 스위치이다. SW1=1이 아니면 공정(S220)에서 t=t인가 여부를 조사한다. 레지스터 t에는 측정개시부근의 경과시간정수 t(예컨대 1초)가 기억되어 있으며, 이것은 제5도의 t의 타이밍에 상당한다. 만일 t=t가 아니면 제3a도의 공정(S108)으로 되돌아 간다. 그러나, 본 실시예에서는 제1초째에 t=T의 조건을 만족시키므로, 공정(S221)으로 진행하여 레지스터 T에 현 시점의 검온데이터 T8을 기억시키고, 공정(S222)에서 SW1을 논리 1로 세트한다. 이후는 SW1의 처리를 하지 않는다.
또, 공정(S203)에서 SW1=1이면, 공정(S204)에서 SW2=1인가 여부를 조사한다.
SW2는 제1의 상승구배를 검출하는 처리용의 스위치이다. 공정(S204)에서 SW2=1이 아니면 공정(S230)에서 |T8-T0|<K인가 여부를 조사한다. 그리고 만일 |T8-T0|<K가 아니면 제1의 상승구배를 아직 검출할 수 없으므로 플로우는 제3a도의 공정(S108)으로 되돌아 간다. 그러나 곧, |T8-T0|<K를 만족시키면 공정(S231)에서 레지스터 S1에 검온데이터의 차분수치
1=T8-T를 기억시킨다. 제5도와 같이 차분수치
1=T8-T는 간접적으로 제2의 상승구배를 표기하고 있다.
또한 공정(S232)에서는 레지스터 T에 현 시점의 검온데이터 T8을 기억시키고, 공정(S232)에서는 레지스터 t에 현 시점의 검온데이터 T를 기억시키며, 공정(S234)에서는 SW2를 논리 1로 세트한다. 이후는 SW2의 처리를 하지 않는다.
또, 공정(S204)에서 SW2=1이면 공정(S205)로 진행하고, SW3=1인가 여부를 조사한다. SW3은 상기 제1의 구배의 검출후, 다시 소정시간 DT(예로서 16초)의 경과를 기다려 검온커브의 형상인식처리를 하기위한 스위치이다. 만일 SW3=1이 아닐때에는 공정(S240)으로 진행하여 t=t+DT인가 여부를 조사한다. 그리고 t=t+DT가 아닌 경우에는 조건을 만족시킬 때까지 제3a도의 공정(S108)으로 되돌아 간다. 그러나, 곧 t=t+DT의 조건을 만족시키면, 다음에 설명하는 공정(S300)의 형상인식처리를 실행한다. 이것은 제5도의 t의 시점이다.
[형상인식처리]
제4a도는 제1실시예에서의 형상인식처리를 표시하는 플로우챠트이다. 공정(S31)에서는 레지스터 S2에 검온데이터의 차분수치 S2=T8-T를 기억시킨다. 이 차분수치 S2=T8-T는 제5도의 제3의 상승구배를 간접적으로 나타내고 있다. 공정(S302)에서는 형상파라미터 a의 값을,
Figure kpo00041
에 의해서 구하고, 공정(S303)에서 제3b도의 공정(S241)로 되돌아 간다. 여기에서 정수 Q1(예컨대 0.042) 및 Q2(예컨대 -0.128)의 각 수치는 통계적 처리에 의해서 설정되고 있다. 그리고 공정(S241)에서는 SW3을 논리 1로 세트하고, 이후는 SW3의 처리를 하지 않는다. 즉, 이후는 체온측정의 초기단계에서 설정한 형상파라미터 a의 값을 사용한다.
또, 공정(S205)에서 SW3=1일때에는 다음에 설명하는 공정(S400)의 곡선해석처리를 실행한다.
[곡선해석처리]
제4b도는 제1실시예에서의 곡선해석처리를 나타내는 플로우챠트이다. 공정(S401)에서는 현 시점의 검온데이터 T8을 사용해서 계수파라미터 A0의 값을 구하고, 이것을 레지스터 A0에 기억시킨다. 여기에서 계수파라미터 A0의 값을 다음식,
Figure kpo00042
에 따라서 구한다. 공정(S402)에서는 상기 구해진 계수파라미터 A0의 값을 사용하여 계수파라미터 A1의 값을 구하고, 이것을 레지스터 A1에 기억시킨다. 여기에서 계수파라미터 A1의 값을 다음식,
Figure kpo00043
에 의해서 구한다. 공정(S403)에서는 8초전의 검온데이터 T0를 사용하여 계수파라미터 A0'의 값을 구하고, 이것을 레지스터 A'에 기억시킨다. 여기에서, 계수파라미터 A'의 값을 다음식,
Figure kpo00044
에 의해서 구한다. 공정(S404)에서는 상기 구하여진 계수파라미터 A0'의 값을 사용하여 계수파라미터 A1'의 값을 구하고, 이것을 레지스터 A1'에 기억시킨다. 여기에서 계수파라미터 A1'의 값은 다음식,
Figure kpo00045
에 의하여 구한다. 공정(S405)에서는 제3b도의 공정(S500)으로 돌아가서, 다음과 같은 예측연산처리를 실행한다.
[예측연산처리]
제4c도는 제1실시예에서의 예측연산처리를 표시하는 플로우챠트이다. 공정(S501)에서는 형상파라미터 a 및 현 시점의 계수파라미터 A0,A1의 각 수치를 사용해서 장래시각 t(예를들면 600초)에서의 예측온도 T(t)를 구하고, 이것을 레지스터 T에 기억시킨다. 여기에서 예측온도 T(t)는 다음식,
Figure kpo00046
에 의하여 구한다. 공정(S502)에서는 형상파라미터 a 및 8초전의 계수파라미터 A0',A1'의 각 수치를 사용해서 장래시각 t에서의 예측온도 T(t)'를 구하고, 이것을 레지스터 T'에 기억시킨다. 여기에서 예측온도 T(t)'는 다음식,
Figure kpo00047
에 의하여 구한다. 공정(S503)에서는 제3b도의 공정(S206)으로 복귀해서, 구하여진 예측온도 T(t)를 표시한다.
또한 플로우는 제3c도의 공정(S207)로 진행하여, SW4=1인가 여부를 조사한다. SW4는 예측온도의 유효성(일정성)을 조사하는 처리의 스위치이다. 공정(S207)에서 SW4=1이 아니면 공정(S250)으로 진행하여, 8초간격으로 구한 예측온도의 차분의 절대치가
Figure kpo00048
인가 여부를 조사한다. 만일
Figure kpo00049
가 아니면 현시점의 T와 8초전의 예측치 T'가 대략 동일하다고 볼 수 없으므로(예측치 T는 유효하지 않다.) 제3a도의 공정(S108)으로 복귀한다. 그러나
Figure kpo00050
이면 예측은 안정되었다고 간주되므로, 공정(S251)에서 부저(14)를 울려, 사용자에 그 취지를 알린다. 공정(S252)에서는 SW4에 논리1을 세트하고, 이후는 SW4의 처리(부저14의 울림)를 행하지 않는다. 그 후에는 어제 체온계를 떼어도, 사용자에게는 유효한 선행온도표시가 얻어진다. 만일, 검온을 빨리 끝내고자 하면 이 시점의 표시치를 검온치로 인식할 수 있다. 그러나 측정중에 몸이 움직여서, 정상적인 측정이 되지 못하였다는 불안이 있는 경우에는, 혹은 미열의 의심이 있고 고정밀도의 검온이 필요한 경우등에는, 다시 검온을 계속한다. 이 경우, 플로우는 제3a도의 공정(S108)으로 돌아가서 다음의 타이머중단발생에 의하여 다시 제3b도의 곡선처리공정(S400)으로부터 예측온도의 표시공정(S202)에서 t>t(가령 600초)를 만족시키게 되면, 검온데이터 T0자체가 평형온도에 가까워지므로, 그 이후는 통상적인 체온측정으로는 이미 예측연산을 계속하는 것이 무의미하게 되며, 따라서 공정(S210)으로 진행하고, 데이터 선택부(13)를 단자 B측에 접속한 이후는 검온데이터 T0그 자체를 표시한다.
제8도 내지 제10도는 제1실시예의 전자체온계에 의한 겨드랑이 밑 검온의 측정경과를 표시하는 그래프도에 관한 것이며, 도면의 횡축은 측정경과시간 t(분)를 표시하고 있고, 종축은 온도 T℃를 표시하고 있다. 그리고 그래프는 검온데이터 T0의 추이와, 예측10분치 T(600)의 추이와, 실측 10분치 Tmax(600)을 표시하고 있으며, 또 "ERROR"는 부저울림(예측유효) 시점의 예측 10분치 T(600)과 실측10분치 Tmax(600)과의 차분 "JE" 는 부저울림(예측유효)의 판단(공정(S250))에 사용한 허용치
Figure kpo00051
의 값을 표시하고 있다.
제8도는 평균적인 승온커브의 경우를 표시하고 있다. 도면에 의하면 평균적인 승온커브는 안정되고 있으며, 형상파라미터 a의 값의 자동설정도 정확하다. 따라서, 예측 10분치T(600)는 부저울림의 당초로부터 실측 10분치 Tmax(600)을 정확하게 표시하고 있으며, 그 후의 추이도 안정(대략 일정)되고 있다.
제9도는 체동에 의하여 승온커브에 변동이 생기고 있는 결과를 표시하고 있다. 또, 형상파라미터 a의 값은 0.154에 설정되고, 제8도의 승온커브보다도 완만하다. 본 발명에 의하면, 이와 같은 경우에도, 당초의 예측 10분치 T(600)는 그후의 승온커브의 변동에 따라 약간 변동할 뿐이다.
이것은 측정의 초기단계이며, 또한 확실한 방법으로 형상파라미터 a의 값을 설정하고 있으므로 그 후에 승온커브 T0가 변동되어도, 이 변동에 따라서 계수파라미터 A,A1의 각 수치가 약간 병동할 뿐이며, 그러므로 예측 10분치T(600)는 안정되게 추이한다. 종래 방식에서는 이와같은 경우에 예측치가 크게 변동하고 있었다.
제10도는 승온커브가 매우 완만하게 상승하고 있는 경우를 표시하고 있다. 이와같은 경우는 종래의 예측 방식으로는 상당히 낮은 경향의 예측을 하고 있었으나, 본 실시예에 의하면 형상파라미터a의 값이 승온커브의 형상보다도 다소 작게 설정되어 있으며, 이에 의하여 측정의 초기단계에서 예측 10분치T(600)는 다소의 상승경향을 나타내고, 또한 이에 의하여 부저의 울림시각이 자동적으로 지연되어, 따라서 부저울림의 당초부터 양호한 예측정밀도를 얻을 수 있게 되어 있다. 이것은 본 발명에 의하여 승온커브의 개시형상을 형상파라미터 a의 값에 반영시키고 있다는 것과, 예측식의 계수파라미터 A0,A1의 각 수치를 연립방정식에 의하여 실측데이터에 의거한 값으로 설정하고 있다는 사실때문임이 틀림없다.
[제2실시예]
체온측정에 있어서는, 체온계의 열특성과 피부측정부위의 상태 및 부위 그 자체의 특성에 의하여, 측정개시부터 열평형에 이르기까지 사이에 관측되는 온도변화의 모양은 다종다양하다. 그러나 체온계의 열특성을 한정하면, 이것과의 관계에서는 상기 온도변화의 모양을 몇 개의 범주로 분류할 수 있다. 그리고, 극히 큰 분류는 겨드랑이 밑, 입속, 직장등의 측정으로 분류하는 것이다.
그런데, 본 발명에 의한 상기 제1실시예의 전자체온계는 그 보편적인 예측방식의 특징에 의하여 그대로도 겨드랑이 밑, 입속, 직장등의 측정에 적합한 것이었다. 그러나 본 발명의 예측방식을 다시 겨드랑이 밑, 입속, 직장등의 측정중의 어느 하나로 한정하여 적용하는 것을 고려하면, 예측정밀도가 더 향상되는 것을 기대할 수 있다.
제2실시예의 특징은 체온측정을 말하자면, 겨드랑이밑(입속, 직장등으로도 가하다)의 측정으로 한정한점이며, 이에 의하여 제1실시예의 것에 비하여 형상파라미터a의 값의 결정방법을 보다 정밀하게 구체화하여 이로써 예측정밀도를 한층향상시키는 점에 있다.
또, 제2실시예의 특징은, 현실의 체온측정에서 발생하는 온도변화의 변동영향을 경감함으로써 예측온도의 신뢰성을 한층 향상시키는 점에 있다.
또, 제2실시예의 특징은, 예측연산에 관한 처리공정을 단순화하고 혹은 감소시킴으로써, CPU2의 처리 부담을 대폭경감시키는 데 있다.
[원리]
체온측정을 가령, 겨드랑이밑 측정으로 한정하면, 그 온도 변화의 곡선군과 형상파라미터 a의 값과의 상관관계를 한층 정밀하게 규정할 수 있다. 이에 의하여 현상인식수단은 온도변화의 형상을 보다 정밀하게 인식할 수 있어서, 따라서 예측연산의 정밀도가 향상한다.
제14도는 제2실시예에서의 형상파라미터 a의 최적치와 그 형상파라미터 a의 값을 결정하기 위한 변수(S/S)와의 관계를 표시하는 그래프이다.
이 그래프도는 예로서 이하의 방법으로 구할 수 있다. 우선 다수인에 대해서 실측한 다수의(바람직하게는 대표적인) 겨드랑이 밑 검온곡선군을 준비한다. 다음에 어떤 겨드랑이 밑 검온곡선군을 선택해서 전자체온계에 입력한다. 그 때에 미리 형상파라미터 a의 임의의 값을 설정하고 고정한다. 또한 상이한 2점에서의 온도데이터 T(t1),T(t2) 및 각 검온시점에서의 시간데이터 t1·t2에 의거해서 다음식의 2원연립방정식,
Figure kpo00052
를 푸는 것에 의하여, 계수파라미터 A0,A1의 각 수치를 구함과 동시에, 상기 설정한 a의 값과, 축차적으로 구하여지는 계수파라미터 A0,A1의 각 수치를 사용해서 예측연산식,
Figure kpo00053
에 의하여 예측 10분치 T(600)을 구한다. 이리하여 전자체온계가 소정의 예측종료조건(예를들면 예측 10분치가 안정하는 상태)을 만족시켰을때에는 말하자면 외부에서 그 시점의 예측 10분치T(600)와 기지의 실측 10분치TOA(600)를 비교함으로써, 예측오차 γ=|T(600)-TOA(600)|가 소정범위 내에 있는지 여부를 판정한다. 다음에 형상파라미터 a의 모든 가능한 설정치에 대해서 상기의 조작을 행한다. 이에 의하여 어떤 겨드랑이밑 검온곡선에 대한 예측오차 γ가 소정범위내로 되는 등의 복수개의 a의 값(a의 그룹)이 얻어진다.
다음에 별도의 겨드랑이밑 검온곡선을 선택하여 전자체온계에 입력하고, 상기의 조작을 행한다. 이에 의하여 별도의 겨드랑이 밑 검온곡선에 대한 예측오차 γ가 소정범위내의 것이 되는 그러한 다른 a의 그룹이 얻어진다. 또한 모든 겨드랑이 밑 검온곡선을 선택해서 상기의 조작을 행한다.
이리하여 모든 겨드랑이 밑 검온곡선에 대한 예측오차γ가 소정범위내의 것이 되는 그러한 모든 a의 그룹이 얻어진다.
또 한편에 있어서, 상기의 모든 겨드랑이 밑 검온곡선에 대한 형상변수(S/S2)를 검출하고, 겨드랑이 밑 검온곡선을 통해서, 형상변수(S/S)와 a의 그룹을 대응시킨다. 이것이 제14도의 그래프이다.
다음에 제14도의 그래프도에 의거하여 제2실시예에서의 형상파라미터 a와 형상변수(S/S)와의 관계를 규정한다. 제14도에서 명백한 바와 같이, 그래프의 형상은 단순한 직선의 형상으로는 되지 않고, 위의 볼록한 모양의 돌기 경향을 나타내고 있다. 그와 같은 관계는,
Figure kpo00054
또는,
Figure kpo00055
단, n<1 : 정수
1∼C3: 정수
의 형태로 표시된다.
그래서, 제2실시예에서는 형상파라미터 a와 형상변수(S/S)와의 관계를 다음식,
Figure kpo00056
와 같이 규정하였다. 여기에서 제2실시예에서 채용하는 각 정수는 제14도의 그래프특성에 대한 회귀적 통계처리에 의하여 결정되고 있으며, 예를들어 n=0.3, Q=0.04467, Q2=-0.330749, Q=0.393626이다.
[구성]
제11a도 및 b도는 제2실시예의 전자체온계의 구체적인 구성을 나타내는 블록도이다. 이하의 설명에 있어서, 제2a도 및 b도의 제1실시예의 전자체온계와 동등한 구성에는 동일번호를 붙여서 그 설명을 생략한다.
또, 동일번호를 붙인 경우에도 제1실시예와 다소 다른 기능등이 있는 부분에 대해서는 그 설명을 가한다. 제11a도 및 b도에서, 예측연산부(2)는 새롭게 피크홀드수단(61) 및 평균화수단(62)을 포함하며, 또한 제2a도 및 b도의 구성과는 다른 방식으로 실현한 형상인식수단(20), 곡선해석수단(21) 및 예측연산수단(22)을 갖춘다. 그와 같은 구성에서 계측제어수단(8)은 본 계측개시전(예비계측시)에는 4초에 1회의 비율로 온도계측회로(5)에 온도데이터 T를 발생시킨다. 이에 의하여 제1실시예의 경우보다도 예비계측주기가 단축하고, 예비측정정밀도가 향상한다 또, 계측제어수단(8)의 레지스터 t4는 정수로서 2초를 기억하고 있다. 이것은 후술하는 평균화수단(62)을 가졌다는 사실에 의하여, 최초에 사용되는 검온데이터 TOA가 본 계측개시로부터 2초경과후에 유효하게 되기 때문이다.
피크홀드수단(61)은 데이터판독기입수단(6)의 판독기입한 검온데이터 T0중 항상 가장 높은 온도치를 검출하여 기억하는 수단이다. 이 때문에 피크홀드수단(61)은 도시하지 않았지만, 적어도 하나의 피크치 TOP를 기억하는 피크치기억수단과, 그 피크치기억수단이 기억하고 있는 피크치 TOP와 새로 입력되는 검온데이터 T0와의 대소를 비교하는 피크치비교수단을 갖춘다. 그와 같은 구성에 있어서, 피크치기억수단은 최초에는 본 계측개시시점의 검온데이터 T0(또는 본 계측개시조건을 만족시킨 것으로 보고 설정하는 설정치 T=30.0℃)를 피크치 TOP로서 기억한다. 다음에 새로운 검온데이터 T0가 입력될때에는 피크치비교수단에 의하여, 피크치기억수단이 기억하고 있는 피크치 TOP와 새로 입력되는 검온데이터 T0와의 대소를 비교함으로써, 만일 T0
Figure kpo00057
OP의 조건을 만족시킬때에는 피크치기억수단은 새로 입력되는 검온데이터 T0를 기억한다.
평균화수단(62)은 본 계측의 진행과 더불어 피크홀드수단(61)이 출력하는 피크치 TOP에 관한 이동평균치 TOA를 구하는 수단이다. 이 때문에 평균화수단(62)은 도시하지는 않으나, 연속적으로 소정수의 피크치 TOP∼T′OP를 기억하는 피크데이터기억수단과, 그 피크데이타기억수단이 기억하는 소정수의 피크치 TOP∼T′OP를 가산하는 피크치가산수단과, 그 피크치가산수단의 가산결과를 상기 소정수로 나누는 나눗셈수단을 갖는다. 그와 같은 구성에 있어서, 피크데이터기억수단은 각 샘플링시점에서의 최신의 피크치 TOP를 기억함과 동시에, 이미 기억하고 있는 가장 오래된 피크치 T'OP를 지워버린다. 따라서, 피크데이터기억수단은 즉 최초의 샘플링시점(t=0초)에서는 피크치 TOP0를 기억하고, 다음 샘플링시점(t=1초)에서는 피크치 TOP0및 TOP1을 기억하며, 다음 샘플링시점(t=2초)에서는 TOP0, TOP1및 TOP2를 기억한다. 이 시점(t=2초)에서, 피크치가산수단은 최초의 가산결과 TOS1=(TOP0+TOP1+TOP2)를 출력하고, 나눗셈수단은 최초의 이동 평균치 TOA1=TOS1/3을 구하여 출력한다. 또한, 다음의 샘플링시점(t=3초)에서는 피크데이터의 기억수단은 피크치 TOP1, TOP2및 TOP3를 기억한다. 이에 의하여, 피크치가산수단은 다음의 가산결과 TOS2=(TOP1+TOP2+TOP3)을 출력하고, 나눗셈수단은 다음의 이동평균치 TOA2=TOS2/3을 구하여 출력한다.
제13도는 제2실시예에서의 검온데이터 T0, 피크데이터 TOP및 이동평균치 TOA의 추이관계를 표시하는 그래프도이다. 도면에서, 검온데이터 T0는 본 계측개시시점(t=0)으로서, 30.0℃를 통과하고, 그후는 단조로운 상승경향을 표시하고 있으나, t=5초에서 일단 하강으로 바뀌고, t=7초에서 다시 상승으로 바뀌고 있다. 이것은 실제의 체온측정중의 체동등 기타의 원인에 의해서 생기는 검온커브의 변동이다. 일반적으로 미소한 변동은 평균화처리하면 제거된다. 그러나 도면과 같은 현저한 하강부분의 검온데이터는 제1실시예 및 제2실시예와 같은 예측처리에 유용하지 않을뿐만 아니라, 예측결과에 악영향을 준다. 그와 같은 경우에도 피크홀드수단(61)은 예컨대 t=5초에서의 피크치 30.5℃를 기억하여 이 값을 t=8초까지 홀드하므로, 예측처리에서는 하강부분의 검온데이터에 의한 악영향이 제거된다. 또 제2실시예에서는 피크홀드처리를 데이터의 평균화처리보다도 먼저 행하기 때문에, 현실의 검온커브 T0와 피크치커브 TOP가 측정의 전영역에서 잘 일치한다. 따라서, 미리 평균화오차가 개입하지 않는다. 이 평균화오차가 개입하지 않는다는 것은 오히려 측정의 초기단계에서 검온커브가 급격한 상승을 나타낼때에 중요하며, 필요할때에는 검온커브의 상승현상에 대한 이동평균치의 상승현상의 지연정도를 정량적으로 평가할 수 있는 잇점을 가지고 있다. 다음에 평균화수단(62)은 피크치커브 TOP의 이동평균을 취함으로써, 피크치커브 TOP의 변동을 평균화하며 예측연산부(2)의 각 처리에 원활한 검온곡선을 제공한다.
기억수단(9)은 평균화수단(62)에서 구한 이동평균치의 검온데이터 TOA를 순차적으로 레지스터 T8으로부터 T0를 향해서 시프트인/시프트아웃트하면서 항상 최신의 것으로부터 합계 9개의 검온데이터 T8∼T0를 기억한다.
형상인식수단(20)은 겨드랑이밑의 측정으로 한정함으로써, 형상파라미터 a의 값과 변수(S1/S2)와의 사이에 보다 정밀한 상관관계를 얻을 수가 있었다. 즉 형상인식수단(20)의 레지스터a는 다음식,
Figure kpo00058
의 연산에 의해서 구한 형상파라미터 a의 값을 기억한다. 또, 형상인식수단(20)의 레지스터 Q3는 새로 정수 Q3를 기억하기 위하여 마련되어 있다.
곡선해석수단(21)에 있어서는, 8초전의 계수파라미터 A0',A1'의 각 수치를 기억하는 레지스터 A0',A1'와, 그 계수 파라미터 A0',A1'의 각 수치를 구하기 위한 연산처리가 생략되어 있다. CPU2의 연산부담을 경감하기 위해서이다.
예측연산수단(22)은 연속되는 9개분의 예측연산치 T(t)∼T(t)'를 기억하기 위한 예측기억수단을 갖추고 있으며, 그 예측기억수단은 각 샘플링시점에서 연산한 최신의 예측연산치 T(t)를 기억함과 동시에, 이미 기억하고 있는 가장 오래된 예측연산치 사용해서 장래시각 t에서의 예측온도 T(t)'를 구하고, T(t)'를 지워버리는 등을 하며, 연속적으로 9개분의 예측연산치 T(t)∼T(t)'를 기억하도록 구성되어 있다. 이에 의하여 예측수단(22)에 있어서는, 8초전의 온도데이터에 의거한 예측치 T(t)'를 구하는 연산처리를 생략할 수 있고, 따라서 CPU2의 연산부담이 대폭 경감된다. 이점에 대해서, 상기 제1실시예의 전자체온계는 기억수단(9)의 온도데이터 T8∼T0를 유효하게 이용가능하다는 사실에 감하여, 제2실시예와 같은 예측기억수단을 생략한 것이다.
[동작]
제12a도∼c도는 제2실시예에서의 형상인식처리, 곡선해석처리 및 예측연산처리를 표시하는 플로우챠트이다. 또 제2실시예에서의 메인 플로우는 제3a도∼e도의 메인플로우와 동등하게 구성되므로, 이하의 설명에서도 제3a도∼e도를 채용한다. 따라서, 제2실시예의 전체의 동작은 제3b도에서의 형상인식 처리 S300, 곡선해석처리 S400 및 예측연산처리 S500을 각각 제12a도∼c도의 형상인식처리 S300', 곡선해석처리 S400', 및 예측연산처리 S500'으로 바꾸어 놓은 것으로 하고 설명한다. 그리고 제1실시예와는 다른 기능이 있는 부분에 대해서는 그 설명을 가한다.
제3a도의 예비측정공정 S101에서, 계측제어수단(8)은 4초에 1회의 비율로 온도계측회로(5)에 온도를 검출시킨다. 예비계측주기를 단축해서 예비계측정밀도를 향상시킨 것이다. 또, 공정 S102에서는 소정온도 30℃를 넘었는가 여부를 판단하고, 공정 S103에서는 4초간에 0.32℃이상의 온도상승이 인지되는지 여부를 판단한다. 제2실시예에서는 예비계측주기를 4초로 단축한 사실에 대응시켜서 0.32℃이상의 온도상승이 있는지의 여부를 검증하기로 하였다. 또, 제3b도의 공정 S210에서는 피크홀드수단(61) 및 평균화수단(62)을 거친후의 검온데이터 TOA를 표시한다. 또, 공정 S220에서는 t=t로 하여 2초인가 여부를 조사한다. 제2실시예에서는 평균화수단(62)를 갖출 수 있으므로, 레지스터 t의 내용은 최초의 검온데이터 TOA가 이용되는 2초로 하였다.
[형상인식처리]
제12a도에서 제4a도와 동등한 공정에는 동일번호를 붙여서 설명을 생략한다. 공정S310에서는 형상파라미터 a의 값을,
Figure kpo00059
에 따라 구한다. 이에 의하여 겨드랑이 밑 검온시의 형상인식정밀도가 크게 향상하였다.
[곡선해석처리]
제12b도에서는 제4b도의 8초전의 검온데이터 T0를 사용하는 연산공정 S403 및 S404가 생략되고 있다. 제2실시예에서는 예측연산수단(22)내에 예측치 9개분을 기억할 수 있는 예측기억수단을 갖추었으므로, 현시점으로부터 8초전의 검온데이터 T0를 사용하여 계수파라미터 A0',A1'를 구할 필요가 없다. 이에 의하여 CPU2의 처리부담이 대폭 경감된다.
[예측연산처리]
제12c도에서 제4c도와 동등한 공정에는 동일공정번호를 붙여서 설명을 생략한다. 제12c도에서는 제4c도의 8초전에 대한 계수파라미터 A0',A1'의 각 수치를 사용하는 연산공정 S502가 생략되고 있으며, 그 대신에 공정 S510을 삽입하여, 최신의 예측연산온도 T(t)를 예측치기억수단에 기억시킴과 동시에, 8초전에 연산하여 기억시킨 예측연산온도 T(t)'를 지워버리는 처리를 하고 있다.
[제3실시예]
상기 제1실시예의 설명에서 서술한 바와같이, 형상파라미터 a의 값 설정방법에는 그외에도 여러 가지 고려할 수 있다. 제3실시예의 특징은 제1실시예의 설명에서 말한 측정개시후의 각 소정시각에서의 온도구배나, 각 소정시각에서의 복수의 온도데이터등에 의해서 형상파라미터 a의 값을 표현하는 방법을 보다 상세히 구체화한데 있다.
또, 제3실시예의 특징은 체온측정을 구강(그외에 겨드랑이 밑, 직장도 가하다)의 측정으로 한정한 사실이며, 여기에 형상파라미터 a의 값에 대한 다른 결정방법을 적용한 점에 있다.
또, 제3실시예의 특징은 상기 제2실시예와 같이, 온도변화의 변동영향을 경감시킨 점에 있으며, CPU2의 처리부담을 대폭 경감시킨데 있다.
[원리]
제17도는 제3실시예에서의 본 계측개시후의 각 소정시각에서의 복수의 검출온도데이터를 플롯트한 관계를 표시하는 그래프도이다. 상기 제1실시예의 설명에서 설명한 바와 같이, 복수의 온도상승곡선은 그 측정 개시점으로부터의 검출온도의 상승형상의 예리한 정도를 인식하는 것으로서 판별된다. 그래서 제3실시예에서는 제17도와 같은 각 소정시각에서의 복수의 온도 데이터에 의거하여 형상파라미터 a의 값을 표현하기로 하였다.
어떤 검출온도의 상승형상의 예리한 정도를 인식하기 위하여는 예로서 실측데이터의 조를 그대로 채용해서 이들을 정량화하거나 혹은 이들을 기준의 것과 비교하는 방법이 있으며, 혹은 실측데이터의 조를 시간축과의 관계에서 보아(1차미분하여) 얻은 온도변화율의 정보를 정량화 혹은 이들을 기준의 것과 비교하는 방법이 있으며, 또한 실측데이터의 복수의 온도변화율의 조를 시간축과 관계에서 보아(2차미분하여) 얻은 온도변화의 가속도정보를 정량화하거나 혹은 이들을 기준의 것과 비교하는 방법등을 고려할 수 있다. 또, 그 경우에 금일 입수 가능한 1칩 CPU의 처리능력등을 생각한다면, 형상파라미터 a의 변수간의 관계를 규정하는 함수는 선형인 것이 바람직하다. 다행히 미분항을 포함할때에는 현실적으로는 간단한 처리로 실행이 가능하다. 즉, 제17도에 있어서, 예로서 구간(t8∼t16)에 있어서의 1차미분은 각 시점의 검출온도데이터의 차(T16∼T8)를 채택하여 할 수 있고, 또, 구간(t8∼t24)에 있어서의 2차미분은 {T24-T16-(T16-T)}의 형태로 행할 수가 있다. 그래서, 제3실시예에서는 형상파라미터 a와 변수 X와의 관계를 다음식,
Figure kpo00060
단,
0∼D5: 정수
0=T81=T16
2=T243=T32
4=(X3-X0)/(X1-X0
와 같이 규정하였다. 여기에서, T8∼T32는, 즉 t=8초, 16초, 24초, 및 32초의 각 시점에서의 검온데이터 TOA이며, 그 복수의 검온데이터 T8∼T32는 바람직하게는 측정개시후의 검출온도의 상승형상의 예리한 정도를 유효하게 인식하고, 판별하는데 충분한 범위를 커버하는 범위에서 추출한 것이다. 또, 이 범위는 바람직하게는 최초의 예측연산온도 T를 보다 빠른 시점에서 표시하기 위하여 가급적 단축하는 것이다. 또한, 변수 X0∼X4의 수, 샘플링간격, 및 형은 바람직하게는 CPU2의 처리에 과도한 부담을 주지 않는 범위내에서 결정하는 것이다. 또한 제3실시예에서의 각 정수 D0∼D5는 이상을 고려한 결과의 통계적 처리에 의거해서 결정된 것이며, 예컨대, D=-0.2566, D1=0.01601, D2=0.03003,D3=0.35019,D4=0.08913,D5=-12.9657이다.
[구성]
재15a도 및 b도는 제3실시예의 전자체온계의 구체적인 구성을 표시하는 블록도이다. 이하의 설명에서, 제11a도 및 b도의 제2실시예의 전자체온계와 동등한 구성에는 동일번호를 붙여서 그 설명을 생략한다.
제15a도 및 b도에 있어서, 예측연산부(2)에서는 제11a도 및 b도의 구성과는 다른 방식으로 실현한 계측제어 수단(23) 및 형상인식수단(24)을 갖추고 있다.
계측제어수단(23)에 있어서는 레지스터t8∼t32는 복수의 온도데이터를 추출할 각 소정시각의 정수, 예를 들면 t8=8초, t16=16초, t32=32초를 기억하고 있다. 그러나, 각 소정시각의 정수는 이 내용에만 한정되는 것은 아니다. 또한 레지스터 I는 인덱스레지스터이며, 그 인덱스레지스터(I)의 내용은 레지스터 t8∼t32의 내용을 인덱스하여 참조하거나, 후술하는 타목적에서 사용된다.
형상인식수단(24)은 즉, 구강측정으로 한정한 사실에 의하여, 형상파라미터 a의 값과, 변수 X와의 사이에 보다 정밀한 상관을 얻을 수가 있었다. 즉, 형상인식수단(24)의 레지스터a는 다음식,
Figure kpo00061
의 연산에 의해서 구한 형상파라미터 a의 값을 기억한다. 또한 형상인식수단(24)에서, 레지스터 X0∼X4는 변수 X0∼X4를 기억하기 위하여 마련되고, 레지스터 D0∼D5는 정수 D0∼D5를 기억하기 위하여 마련되고 있다.
[동작]
제16a도∼e도는 제3실시예의 전자체온계에 의한 검온프로세스를 나타내는 플로우챠트이다. 제16a도에서 제3a도와 동등한 공정에는 동일번호를 붙여서 설명을 생략한다. 또, 동일공정번호를 붙인 경우에도, 제1실시예와 다소 다른 기능이 있는 부분에 대해서는 그 설명을 가한다. 제16a도의 온도계측공정 S101에서는 계측제어수단(8)이 4초에 1회비율로 온도계측회로(5)에 온도를 검출시킨다. 또, 공정 S102에서는 소정온도 30℃를 초과하였는가 여부를 판단하고, 공정 S103에서는 4초간에 0.32℃이상의 온도상승이 인지되었는가 여부를 판단한다. 그리고, 상기 모든 조건을 만족시키고 있을때에는 공정 S110으로 진행하고, 계측제어를 위한 각종 프로그램스위치 SW1∼SW4 및 인덱스레지스터 I의 내용을 클리어한다.
제16a도으로 공정 S109에서 타이어인터럽트가 발생하면 제16b도의 공정 S200에 입력한다. 제16b도의 공정 S210에서는 검온데이터 TOA를 표시한다.
또 공정 S220에서는 t=t로 하여 2초인가 여부를 조사한다. 또, 공정 S221에서는 레지스터 T에 검온 데이터 TOA를 기억시킨다. 제3실시예에서는 기억수단(9)이 생략되어 있기 때문이다.
공정 S205에서는 SW3=1인가 여부를 조사한다. 제3실시예에서는 SW3은 본 계측개시후의 각 소정시각에서의 복수의 온도데이터를 기억하고, 또한 측정개시후의 비교적 빠른 시점에서 검온곡선의 상승형상의 예리한 정도를 인식(형상파라미터 a의 값을 결정)하기 위한 스위치이다. SW3=1이 아닐때에는 공정 S260으로 진행하여, t1=t8, t16,t24및 t32중 어느 하나인가 아닌가를 조사한다. 이 판정은 인덱스레지스터 I의 내용을 인덱스하여 참조하면서 행한다. 즉 최초는 I=0이며, t=8초를 기다린다. t=8초가 아닌 경우에는 그때까지의 데이터는 형식인식에 유용하지 않으므로, 그대로 공정 S108로 복귀한다. 머지않아 t=8초가 되면, 공정 S260의 판별이 YES가 된다. 이에 의하여 공정S261에서는 인덱스레지스터 I의 내용에 따라서 검온데이터 TOA(최초는 T)를 형상인식수단(24)의 레지스터 X에 기억시킨다. 즉 최초는 I=0이므로 레지스터 X0에 T8을 기억시킨다. 공정 S262에서는 마찬가지로 인덱스레지스터 I의 내용에 따라서 형상인식수단(24)이 형상파라미터 결정을 위한 부분연산, a=f(I,X,D)를 실행한다. 즉, 최초는 I=0이므로 a=(D0, X0)의 부분연산을 행한다. 공정 S263에서는 인덱스레지스터 I의 내용에 1을 가한다. 공정S264에서는 I=3인가 여부를 판별하여, I=3이 아니면, 공정 S108로 돌아간다. 이와 같이 하여 순서대로 소정시각이 되면 상기와 같은 요령의 처리를 한다. 즉, 다음t1=16초에서는 I=1이므로 공정S261에서는 레지스터 X1에 T16을 기억시키고, 공정S262에서는 레지스터 a의 내용에 (D2,X2)의 내용을 가하며, 공정S263에서는 I의 내용에 1을 가한다. 다음의 t1=2초에서는 I=3이므로 I=3이므로, 공정 S261에서는 레지스터 X3에 T32을 기억시키고, 공정 S262에서는 레지스터 a의 내용에 (D3, X3)의 내용을 가하며, 공정 S263에서는 I의 내용에 I을 가한다. 이 시점에서는 공정S264의 판별은 I=3을 만족시키므로 다음에 설명하는 공정S600의 형상인식처리를 행한다. 이것은 제3실시예에서는 항상 본 계측개시후의 t=32초를 경과한 시점이다.
[형상인식처리]
제16d도는 제3실시예에서의 형상인식처리를 나타내는 플로우챠트이다. 공정S601에서는 레지스터 X4에 변수 X4=(X3-X0)/(X1-X0)를 기억시킨다. 공정S602에서는 레지스터a의 내용에 나머지(D4·X4+D5)의 내용을 가한다. 이에 의하여 형상파라미터 a의 값이 계산되었다. 다시 공정S602에서는 a>1인가 여부를 판별하여, 만일 a>1이면, 공정 S606에서, a의 값을 1에 클램프한다. 또, a>1이 아니면 공정 S604에서 a〈0.01인가 여부를 판별하고, 만일 a<0.01이면 a의 값을 0.01에 클램프한다. 또, a<0.01도 아닌때에는 공정 S602의 연산으로 구한 a의 값 그 자체를 사용한다. 공정S607에서는 메인플로우(main flow)로 돌아가고, 공정S265에서는 SW을 논리 1에 셋트한다. 이후는SW3의 처리를 행하지 않는다. 즉, 이후는 체온측정의 초기단계에서 설정한 형상파라미터 a의 값을 사용해서 예측연산을 계속한다.
제18도 및 제19도는 제2실시예의 전자체온계에 의한 겨드랑이 밑 검온의 측정결과를 표시하는 그래프에 관한 것이며, 제18도는 승온커브가 매우 완만하게 상승하고 있는 경우를 표시하고, 제19도는 평균적인 승온커브의 경우를 표시하고 있다. 제2실시예의 전자체온계에서는 형상파라미터 a의 값은 승온커브의 소정어깨부근의 형상을 검출한 시점에서 결정되므로, 부저의 울림시점도 승온커브의 상승형상변화에 따라 변화하고 있다.
제20도 내지 제22도는 제3실시예의 전자체온계에 의한 구강내 검온의 측정경과를 표시하는 그래프도에 관한 것이며, 제20도는 평균적인 승온커브의 경우를 표시하고, 제21도는 제20도 표시의 경우보다도 높은 온도를 측정한 경우를 표시하며, 제22도는 승온커브가 매우 완만하게 상승하고 있는 경우를 표시하고 있다. 제3실시예의 전자체온계에서는 형상파라미터a의 값은 항상 t=32초의 시점에서 결정되므로, 부저의 울림 시점에서도 별로 변화가 없다.
그리고, 상술한 각 실시예는 인체에 사용하는 전자체온계에 관해서 설명하였으나, 기타의 생체에 사용하는 경우에도, 혹은 생체이외의 것의 검온에 사용하는 경우에도 본 발명, 혹은 본 발명의 사상을 용이하게 적용할 수 있다.
또, 상술한 제2실시예 및 제3실시예의 설명에서는 각각 겨드랑이 밑 전용 및 구강전용으로서 실현한 경우의 전자체온계에 관해서 설명하였다. 그러나, 어느 계측알고리즘(특히 형상파라미터a의 값의 결정알고리즘)도 겨드랑이 밑용, 구강용등에 한정되는 것은 아니다. 또한 하나의 전자체온계에 쌍방의 계측알고리즘을 적용해서, 선택적으로 기능발휘시키는 것도 가능하다.
또, 제1실시예 또는 제2실시예에서의 제2의 구배S1과 제3의 구배S2와의 비교는 이들의 비를 찾는 것만에 한하지 않는다. 즉, 차를 찾아서 택하여도 좋다.
이상 설명한 것과 같이, 본 발명에 의하면, 예측식중의 파라미터의 값을 전부 측정시의 실시간온도데이터를 사용해서 산출하고 있으므로, 탐침의 열특성불규칙 또는 개인이나 측정부위의 차에 의한 검온커브에 차가 있어도, 항상 정확한 선행표시를 할 수 있다.
또, 본 발명에 의하면, 실시간온도데이터 그 자체를 목적변수로 하고 있으므로, 연산오차에 의한 영향이 없고, 파라미터의 값설정이 안정되며, 실측온도커브에 중첩된 잡음에 의해서도 예측치가 크게 좌우되거나 하지 않는다.
또, 본 발명에 의하면, 측정경과와 더불어 승온커브를 전부 커버하도록 온도데이터를 끌어내므로, 승온커브가 매우 완만하게 상승하는 경우에도 예측치의 추이는 자연스러운 상승커브를 그리며, 온도의 상승부근에서 오우버슈우트하는 일이 없고, 사용자에게 예측하고 있다는 것을 의식시키지 않고 측정할 수가 있다.
또, 본 발명에 의하면, 측정의 초기단계에서 승온커브의 형상을 정확히 판단함으로써, 예측온도치는 당초부터 평형온도를 나타낸다.
또, 본 발명에 의하면, 예측식에 대해서 임의의 장래시각을 직접적으로 설정할 수 있으므로, 당해 측정계에서 임의 측정시간 경과후의 검온치를 용이하게 제공할 수 있다.
또, 본 발명에 의하면, 곡선형상a의 값을 설정한 후에 2원연립방정식을 풀게되므로, 1점의 온도데이터는 측정계시부근의 온도데이터로 고저할 수 있고, 또 다른 1점의 온도데이터는 현시점의 온도데이터를 사용하면 된다. 따라서, 종래와 같은 온도데이터메모리는 불필요하게 되고, 연산도 간단하고 또한 고속으로 행할 수 있다.
따라서, 보편성이 높고, 염가이며 고정밀도인 전자체온계를 제공할 수 있다.
또, 본 발명에 의하면, 생체의 체동따위에 의하여 발생하는 온도변화의 변동을 유효하게 방지 또는 경감시키므로, 예측온도의 신뢰성이 한층 향상된다.
또, 본 발명에 의하면, 측정대상을 생체의 특정부위로 한정하였으므로, 예측정밀도가 한층 향상된다.

Claims (17)

  1. 장래시각에서의 검출온도예측기능을 갖추는 생체의 온도측정방법에 있어서, 검온커브의 형상을 반영하기 위한 형상파라미터 및 상기 검온커브에 중복시키기 위한 계수파라미터의 각 수치를 미정으로 하는 소정의 예측함수식을 갖추는 공정과, 온도를 검출하여 그 온도를 표시하는 온도데이터를 발생하는 온도검출공정과, 측정개시후의 경과시간을 제시하여 그 경과시간을 표시하는 시간데이터를 발생하는 시간검출공정과, 상기 검출된 소정의 온도데이터에 의거해서 상기 형상파라미터의 값을 설정하는 형상파라미터의 설정공정과, 상기 설정된 형상파라미터치를 포함하여, 또한 상이한 복수시점의 온도데이터를 목적변수로 하고, 또한 상기 복수시점의 시간데이터함수를 설명변수로 하는 복수의 상기 예측함수식으로서 구성되는 연립방정식을 풀므로써 상기 계수파라미터의 값을 설정하는 계수파라미터설정공정과, 상기 설정된 형상파라미터 및 계수파라미터의 각 값으로 특정하는 상기 예측함수식에 의해 장래시각에서의 검출온도를 예측연산하는 예측연산공정을 갖추는 것을 특징으로 하는 생체의 온도측정방법.
  2. 제1항에 있어서, 형상파라미터설정공정은 복수의 온도데이터에 의거하는 소정의 온도상승구배정보에 의해서 형상파라미터의 값을 설정하는 것을 특징으로 하는 생체의 온도측정방법.
  3. 제2항에 있어서, 형상파라미터설정공정은 우선 피검온커브가 제1의 소정구배를 표시하는 점을 검출하고, 또한 그 점에 관해서 전후하는 제2의 구배 및 제3의 구배를 검출하여 이것을 비교함으로써, 형상파라미터의 값을 설정함을 특징으로 하는 생체의 온도측정방법.
  4. 제3항에 있어서, 형상파라미터설정공정은 제2의 구배, S1및 제3의 구배 S2에 의거해서 형상파라미터 a의 값을 다음식,
    Figure kpo00062
    단,
    n<1 : 정수
    1∼Q3: 정수
    에 따라 설정함을 특징으로 하는 생체의 온도측정방법.
  5. 제1항에 있어서, 형상파라미터설정공정은 측정개시후의 초기단계에서의 복수의 온도데이터에 의거해서 형상파라미터의 값을 설정하는 것을 특징으로 하는 생체의 온도측정방법.
  6. 제5항에 있어서, 형상파라미터설정공정은 각 소정시각에서의 복수의 온도데이터 T에 의거하는 정보 X에 근거해서 형상파라미터 a의 값을 다음식,
    Figure kpo00063
    단,
    0∼D5: 정수
    0=T81=T16
    2=T243=T32
    4=(X3-X0)/(X1-X0
    에 따라서 설정함을 특징으로 하는 생체의 온도측정방법.
  7. 제1항에 있어서, 계수파라미터설정공정은 상이하는 2점에서의 온도데이터 T(t1),T(t2) 및 각 검온시점에서의 시간데이터 t2,t2에 의거해서 다음식의 2원연립방정식,
    Figure kpo00064
    를 풀므로써 계수파라미터 A0및 A1의 각 수치를 설정함을 특징으로 하는 생체의 온도측정방법.
  8. 제7항에 있어서, 계수파라미터설정공정은 상이하는 2점에서의 온도데이터로서 측정개시점부근의 온도 데이터와 현시점의 온도데이터를 사용함을 특징으로 하는 생체의 온도측정방법.
  9. 제1항에 있어서, 예측연산공정은 형상파라미터 a 및 계수파라미터 A0및 A1의 각 수치로 특정하는 예측함수식에 의하여 장래시각 t에서의 검출온도 T(t)를 다음식,
    Figure kpo00065
    에 따라서 예측연산함을 특징으로 하는 생체의 온도측정방법.
  10. 장래시각에서의 검출온도 예측기능을 갖추는 생체의 온도측정장치에 있어서, 검온커브의 형상을 반영하기 위한 형상파라미터 및 상기 검온커브에 중복시킥시 위한 계수파라미터의 각 수치르 미정으로 하는 소정의 예측함수식을 기억하고 있는 기억수단과, 온도를 검출하여 그 온도를 표시하는 온도데이터를 발생하는 온도검출수단과, 측정개시후의 경과시간을 계시하여 그 경과시간을 표시하는 시간데이터를 발생하는 시간검출수단과, 상기 검출된 소정의 온도데이터에 의거해서 상기 형상파라미터의 값을 설정하는 형상파라미터의 설정수단과, 상기 설정된 형상파라미터치를 포함하여, 또한 상이한 복수시점의 온도데이터를 목적변수로 하고, 또한 상기 복수시점의 시간데이터함수를 설명변수로 하는 복수의 상기 예측함수식으로서 구성되는 연립방정식을 풀므로써 상기 계수파라미터의 값을 설정하는 계수파라미터설정수단과, 상기 설정된 형상파라미터 및 계수파라미터의 각 값으로 특정하는 상기 예측함수식에 의해 장래시각에서의 검출온도를 예측연산하는 예측연산수단을 갖추는 것을 특징으로 하는 생체의 온도측정장치.
  11. 제10항에 있어서, 온도검출수단은 소정주기에서 검출한 온도중 최대레벨의 것을 순차검출하여 유지출력하는 피크홀드수단을 갖추는 것을 특징으로 하는 생체의 온도측정장치.
  12. 제10항에 있어서, 온도검출수단은 소정주기에서 검출한 온도중 최대레벨의 것을 순차검출하여 유지하는 피크홀드수단과, 상기 피크홀드수단이 유지한 복수의 레벨에 관한 이동평균치를 구하여 출력하는 평균화 수단을 갖추는 것을 특징으로 하는 생체의 온도측정장치.
  13. 제10항에 있어서, 형상파라미터설정수단은 복수의 온도데이터에 의거하는 소정의 온도상승구배정보에 의해서 형상파라미터의 값을 설정함을 특징으로 하는 생체의 온도측정장치.
  14. 제10항에 있어서, 형상파라미터설정수단은 측정개시후의 초기단계에서의 복수의 온도데이터에 의거해서 형상파라미터의 값을 설정하는 것을 특징으로 하는 생체의 온도측정장치.
  15. 제10항에 있어서, 계수파라미터설정수단은 상이하는 2점에서의 온도데이터T(t1),T(t2) 및 각 검온시점에서의 시간데이터 t1,t2에 의거해서 다음식의 2원연립방정식,
    Figure kpo00066
    를 풀므로써 계수파라미터 A0및 A1의 각 수치를 설정함을 특징으로 하는 생체의 온도측정장치.
  16. 제15항에 있어서, 계수파라미터설정공정은 상이하는 2점에서의 온도데이터로서 측정개시점부근의 온도데이터와 현시점의 온도데이터를 사용함을 특징으로 하는 생체의 온도측정장치.
  17. 제10항에 있어서, 예측연산공정은 형상파라미터 a및 계수파라미터 A0및 A1의 각 수치로 특정하는 예측함수식에 의하여 장래시각 t에서의 검출온도 T(t)를 다음식,
    Figure kpo00067
    에 따라서 예측연산함을 특징으로 하는 생체의 온도측정장치.
KR1019880701033A 1986-12-24 1987-10-23 생체의 온도측정방법 및 그 장치 KR910005476B1 (ko)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
JP61-306357 1986-12-24
JP30635786 1986-12-24
JP62104483A JPH0795004B2 (ja) 1986-12-24 1987-04-30 生体の温度測定装置
JP62-104483 1987-04-30
PCT/JP1987/000816 WO1988005161A1 (en) 1986-12-24 1987-10-23 Thermometric method for living organism and apparatus therefor

Publications (2)

Publication Number Publication Date
KR890700222A KR890700222A (ko) 1989-03-10
KR910005476B1 true KR910005476B1 (ko) 1991-07-31

Family

ID=26444945

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1019880701033A KR910005476B1 (ko) 1986-12-24 1987-10-23 생체의 온도측정방법 및 그 장치

Country Status (9)

Country Link
US (1) US5015102A (ko)
EP (1) EP0338076B1 (ko)
JP (1) JPH0795004B2 (ko)
KR (1) KR910005476B1 (ko)
AU (1) AU608549B2 (ko)
CA (1) CA1293567C (ko)
DE (1) DE3789874T2 (ko)
ES (1) ES2007434A6 (ko)
WO (1) WO1988005161A1 (ko)

Families Citing this family (45)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5150969A (en) * 1990-03-12 1992-09-29 Ivac Corporation System and method for temperature determination and calibration in a biomedical probe
JP2769225B2 (ja) * 1990-03-27 1998-06-25 日新製鋼株式会社 ヒートストリーク発生の予測方法及び温度センサーロール
US5259389A (en) * 1990-10-24 1993-11-09 Terumo Kabushiki Kaisha Electronic clincal thermometer
US5158087A (en) * 1992-01-31 1992-10-27 Hewlett-Packard Company Differential temperature measurement for ultrasound transducer thermal control
JPH0910313A (ja) * 1995-06-26 1997-01-14 Sanyo Electric Co Ltd 感覚制御装置
US5730134A (en) * 1996-09-09 1998-03-24 General Electric Company System to monitor temperature near an invasive device during magnetic resonance procedures
WO1999056629A1 (fr) 1998-05-06 1999-11-11 Matsushita Electric Industrial Co., Ltd. Thermometre auriculaire a usage feminin
ATE348565T1 (de) * 1998-10-20 2007-01-15 Omron Healthcare Co Ltd Infrarotthermomether
US6270252B1 (en) 1999-05-18 2001-08-07 Alaris Medical Systems, Inc. Predictive temperature measurement system
WO2002071934A1 (fr) * 2001-03-11 2002-09-19 Lemark Mikhailovich Klukin Procede de diagnostic thermique de la pathologie d'un objet biologique et dispositif pour mettre en oeuvre ce procede
EP1249691A1 (en) * 2001-04-11 2002-10-16 Omron Corporation Electronic clinical thermometer
US6886978B2 (en) * 2001-06-18 2005-05-03 Omron Corporation Electronic clinical thermometer
US7484887B2 (en) * 2003-02-20 2009-02-03 Ysis Incorporated Digitally modified resistive output for a temperature sensor
US7318004B2 (en) * 2005-04-01 2008-01-08 Cardinal Health 303, Inc. Temperature prediction system and method
JP4949649B2 (ja) * 2005-07-12 2012-06-13 テルモ株式会社 電子体温計及びその制御方法
US7275866B2 (en) * 2005-11-09 2007-10-02 Chao-Man Tseng Digital thermometer
US7507019B2 (en) * 2006-05-19 2009-03-24 Covidien Ag Thermometer calibration
US20070268954A1 (en) * 2006-05-19 2007-11-22 Sherwood Services Ag Portable test apparatus for radiation-sensing thermometer
US20070268952A1 (en) * 2006-05-19 2007-11-22 Sherwood Services Ag Thermometer calibration by immersion in non-electrically conductive liquid
US7549792B2 (en) 2006-10-06 2009-06-23 Covidien Ag Electronic thermometer with selectable modes
EP1988375A1 (en) * 2007-05-03 2008-11-05 Renishaw plc Method and apparatus for rapid temperature measurement
CN101435727B (zh) * 2007-11-12 2011-01-26 深圳迈瑞生物医疗电子股份有限公司 温度预测方法及装置
EP3591796B1 (en) 2008-03-07 2020-09-30 Milwaukee Electric Tool Corporation Battery pack for use with a power tool and a non-motorized sensing tool
US20090257469A1 (en) * 2008-04-09 2009-10-15 Jones Mike N Infrared thermometer
US8274278B1 (en) * 2008-12-10 2012-09-25 University Of South Florida Noncontact anhysteric curve plotter and static field to time-varying hysteresisgraph with integrated temperature chamber
US8136983B2 (en) * 2009-01-01 2012-03-20 Mahmoud Razzaghi Sensor and control system
CN102802506A (zh) * 2009-04-15 2012-11-28 亚利桑特保健公司 深部组织温度探测器结构
WO2010120360A1 (en) * 2009-04-15 2010-10-21 Arizant Healthcare Inc. Deep tissue temperature probe constructions
US8226294B2 (en) * 2009-08-31 2012-07-24 Arizant Healthcare Inc. Flexible deep tissue temperature measurement devices
US8794829B2 (en) * 2009-12-31 2014-08-05 Welch Allyn, Inc. Temperature-measurement probe
US8292502B2 (en) 2010-04-07 2012-10-23 Arizant Healthcare Inc. Constructions for zero-heat-flux, deep tissue temperature measurement devices
US8292495B2 (en) 2010-04-07 2012-10-23 Arizant Healthcare Inc. Zero-heat-flux, deep tissue temperature measurement devices with thermal sensor calibration
US9354122B2 (en) 2011-05-10 2016-05-31 3M Innovative Properties Company Zero-heat-flux, deep tissue temperature measurement system
US9271672B2 (en) 2012-04-02 2016-03-01 Podimetrics, Inc. Method and apparatus for indicating the emergence of an ulcer
US9377368B2 (en) * 2012-11-16 2016-06-28 Cochlear Limited Temperature sensor and temperature probe
CA2943381C (en) 2014-03-21 2023-10-03 Podimetrics, Inc. Method and apparatus of monitoring foot inflammation
US20150300888A1 (en) * 2014-04-21 2015-10-22 National Taiwan University Temperature prediction system and method thereof
US9561393B2 (en) * 2014-05-16 2017-02-07 Stephen D. Shoap Fluid hose apparatus with integrated nozzle and related systems and methods
US10401235B2 (en) 2015-09-11 2019-09-03 Qualcomm Incorporated Thermal sensor placement for hotspot interpolation
CN108471946B (zh) 2015-11-06 2023-03-17 珀迪迈垂克斯公司 用于溃疡或前溃疡检测的鞋类系统
RU2018126877A (ru) * 2015-12-21 2020-01-23 Конинклейке Филипс Н.В. Способ прогнозирования температуры стабилизации датчика теплового потока
CN106679848B (zh) * 2017-03-27 2018-12-25 成都天衡电科科技有限公司 一种利用普通温度传感器快速测温的方法
EP3866675A4 (en) 2018-10-15 2022-07-06 Podimetrics, Inc. IPSILATERAL ULCER AND METHODS AND DEVICE FOR DETECTING PRE-ULCERATIONS
US11203023B2 (en) 2019-12-20 2021-12-21 Stephen D. Shoap Modular fluid spray nozzles and related systems and methods
WO2023107483A1 (en) 2021-12-06 2023-06-15 Podimetrics, Inc. Apparatus and method of measuring blood flow in the foot

Family Cites Families (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3111032A (en) * 1961-12-07 1963-11-19 Gen Electric Temperature measurement system
US3280312A (en) * 1962-12-11 1966-10-18 Ibm System for measuring kiln conditions using thermocouples in combination with a mathematical model
US3702076A (en) * 1970-06-15 1972-11-07 Ivac Corp Electronic thermometer
US3878724A (en) * 1972-11-16 1975-04-22 Gerald F Allen Physically realizable optimum predictor filter technique for transducers, sensors, and other massive devices, based upon inverse transfers function synthesis
US3978325A (en) * 1973-09-05 1976-08-31 Control Electronics Co., Inc. Electronic thermometer
US3878728A (en) * 1973-10-18 1975-04-22 Us Health Thermesthesiometer
US4068526A (en) * 1975-07-29 1978-01-17 Control Electronics Co., Inc. Electronic thermometer
JPS5327084A (en) * 1976-08-25 1978-03-13 Hitachi Ltd Apparatus for early estimation of converged te mperature
US4092863A (en) * 1977-01-12 1978-06-06 Johnson & Johnson Electronic thermometer circuit
US4232682A (en) * 1978-08-21 1980-11-11 United States Surgical Corporation Physiological parameter measuring system
JPS5571919A (en) * 1978-11-24 1980-05-30 Toshiba Corp Measurement of temperature
JPS55125423A (en) * 1979-03-22 1980-09-27 Canon Inc Electronic apparatus with thermometer fitted
US4541734A (en) * 1982-06-24 1985-09-17 Terumo Kabushiki Kaisha Electronic clinical thermometer, and method of measuring body temperature
JPS58225326A (ja) * 1982-06-24 1983-12-27 Terumo Corp 電子体温計
JPS58225325A (ja) * 1982-06-24 1983-12-27 Terumo Corp 電子体温計
DE3346285A1 (de) * 1982-12-21 1984-10-11 Terumo K.K., Tokio/Tokyo Elektronisches klinisches thermometer und verfahren zur koerpertemperaturmessung
US4574359A (en) * 1982-12-21 1986-03-04 Terumo Kabushiki Kaisha Electronic clinical thermometer, and method of measuring body temperature
JPS59187233A (ja) * 1983-04-06 1984-10-24 Omron Tateisi Electronics Co 電子体温計
JPS60157031A (ja) * 1984-01-27 1985-08-17 Terumo Corp 電子体温計
US4727500A (en) * 1985-05-01 1988-02-23 Sherwood Medical Company Electronic thermometer with fixed response time
JPS6170429A (ja) * 1985-08-13 1986-04-11 Terumo Corp 熱平衡温度の予測方法
JPS62165132A (ja) * 1986-01-16 1987-07-21 Omron Tateisi Electronics Co 電子体温計
US4878184A (en) * 1986-02-10 1989-10-31 Omron Tateisi Electronics Co. Electronic thermometer with predicting means
JPH0625700B2 (ja) * 1986-03-04 1994-04-06 テルモ株式会社 電子体温計
JPH0792405B2 (ja) * 1986-05-13 1995-10-09 オムロン株式会社 電子体温計
GB2197724B (en) * 1986-11-05 1990-12-19 Citizen Watch Co Ltd Predictive operation type electronic clinical thermometer
US4986669A (en) * 1986-11-19 1991-01-22 Terumo Kabushiki Kaisha Electronic clinical thermometer
JP2797872B2 (ja) * 1992-12-08 1998-09-17 住友金属工業株式会社 タンデムミルの板厚制御方法

Also Published As

Publication number Publication date
ES2007434A6 (es) 1989-06-16
US5015102A (en) 1991-05-14
AU8104787A (en) 1988-07-27
WO1988005161A1 (en) 1988-07-14
EP0338076A4 (en) 1991-08-07
KR890700222A (ko) 1989-03-10
JPS63265126A (ja) 1988-11-01
DE3789874T2 (de) 1994-09-01
EP0338076A1 (en) 1989-10-25
JPH0795004B2 (ja) 1995-10-11
EP0338076B1 (en) 1994-05-18
DE3789874D1 (de) 1994-06-23
CA1293567C (en) 1991-12-24
AU608549B2 (en) 1991-04-11

Similar Documents

Publication Publication Date Title
KR910005476B1 (ko) 생체의 온도측정방법 및 그 장치
EP1183508B1 (en) Predictive temperature measurement system
US5066141A (en) Electronic clinical thermometer
US5738441A (en) Electronic clinical predictive thermometer using logarithm for temperature prediction
EP0332700A1 (en) Electronic thermometer
US4691713A (en) Electronic clinical thermometer
US4866621A (en) Predictive operation type electronic clinical thermometer
US5259389A (en) Electronic clincal thermometer
GB2133563A (en) Electronic clinical thermometer and method of measuring temperature
JPS62165132A (ja) 電子体温計
EP1906162B1 (en) Electronic clinical thermometer, control method and control program for electronic clinical thermometer
JPH0259417B2 (ko)
EP0482562B1 (en) Electronic clinical thermometer
JP2000111415A (ja) 婦人体温計
KR910003827B1 (ko) 전자체온계
JP3389033B2 (ja) 婦人体温計
JPS6170429A (ja) 熱平衡温度の予測方法
JPH0656333B2 (ja) 電子体温計
JPS60147621A (ja) 電子体温計
JPS62192625A (ja) 電子体温計
JPH03567B2 (ko)

Legal Events

Date Code Title Description
A201 Request for examination
G160 Decision to publish patent application
E701 Decision to grant or registration of patent right
GRNT Written decision to grant
FPAY Annual fee payment

Payment date: 20060725

Year of fee payment: 16

LAPS Lapse due to unpaid annual fee