知財求人 - 知財ポータルサイト「IP Force」
特開2024-165466四肢動物の生体信号処理方法、生体信号処理装置、生体信号処理システム、プログラム、および記憶媒体
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】公開特許公報(A)
(11)【公開番号】P2024165466
(43)【公開日】2024-11-28
(54)【発明の名称】四肢動物の生体信号処理方法、生体信号処理装置、生体信号処理システム、プログラム、および記憶媒体
(51)【国際特許分類】
A01K 29/00 20060101AFI20241121BHJP
【ＦＩ】
A01K29/00 A
【審査請求】未請求
【請求項の数】14
【出願形態】ＯＬ
(21)【出願番号】P 2023081691
(22)【出願日】2023-05-17
(71)【出願人】
【識別番号】000003160
【氏名又は名称】東洋紡株式会社
(72)【発明者】
【氏名】小松 陽子
(57)【要約】
【課題】四肢動物が臥位であると推測される時間期間を特定し、当該時間期間における当該四肢動物の生体信号の処理を行うことを可能とする技術を提供することを目的とする。
【解決手段】この発明の一態様の四肢動物の生体信号処理方法は、四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得することと、前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定することと、前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行うこととを含む。
【選択図】図４
【解決手段】この発明の一態様の四肢動物の生体信号処理方法は、四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得することと、前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定することと、前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行うこととを含む。
【選択図】図４
【特許請求の範囲】
【請求項1】
コンピュータおよび記憶媒体を備える装置が実行する四肢動物の生体信号処理方法であって、
四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得することと、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定することと、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行うことと
を含む、生体信号処理方法。
四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得することと、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定することと、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行うことと
を含む、生体信号処理方法。
【請求項2】
前記取得される加速度の値は、互いに交わるＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度ｘ、ｙ、およびｚのベクトル和の大きさに対応する、請求項１に記載の生体信号処理方法。
【請求項3】
前記加速度の値のばらつきの度合いを示す値は、標準偏差、変動係数、および分散のいずれかである、請求項１に記載の生体信号処理方法。
【請求項4】
前記第１の時間期間は、前記ばらつきの度合いを示す値が前記第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間である、請求項１に記載の生体信号処理方法。
【請求項5】
前記第１の時間期間より後の、前記ばらつきの度合いを示す値が第２の閾値以下の第２の時間期間を特定することをさらに含み、
前記出力は、前記生体信号のうちの前記第２の時間期間における生体信号にさらに基づいて行われる、
請求項１に記載の生体信号処理方法。
前記出力は、前記生体信号のうちの前記第２の時間期間における生体信号にさらに基づいて行われる、
請求項１に記載の生体信号処理方法。
【請求項6】
前記第１の時間期間は、前記ばらつきの度合いを示す値が前記第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間であり、
前記第２の時間期間は、前記ばらつきの度合いを示す値が前記第２の閾値以下となる時点から第２の時間経過した時点以降に開始される時間期間である、
請求項５に記載の四肢動物の生体信号処理方法。
前記第２の時間期間は、前記ばらつきの度合いを示す値が前記第２の閾値以下となる時点から第２の時間経過した時点以降に開始される時間期間である、
請求項５に記載の四肢動物の生体信号処理方法。
【請求項7】
請求項５に記載の生体信号処理方法であって、
前記生体信号が、拍動を示す信号を含み、
前記第１の時間期間における生体信号と、前記第２の時間期間における生体信号との各々について、拍動間隔の平均値、単位時間当たりの拍数、拍動間隔のばらつき、単位時間当たりの拍数のばらつき、ＨＦ、およびＬＦ／ＨＦのうちの少なくとも一つを算出し、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されるか否かを判定することをさらに含み、
前記出力は、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されると判定された場合に行われる、
生体信号処理方法。
前記第１の時間期間における拍動間隔の平均値＞前記第２の時間期間における拍動間隔の平均値・・・（ｃ１）
前記第１の時間期間での単位時間当たりの拍数＜前記第２の時間期間での単位時間当たりの拍数・・・（ｃ２）
前記第１の時間期間における拍動間隔のばらつき＞前記第２の時間期間における拍動間隔のばらつき・・・（ｃ３）
第１の時間期間での単位時間当たりの拍数のばらつき＞第２の時間期間での単位時間当たりの拍数のばらつき・・・（ｃ４）
前記第１の時間期間のＨＦ＞前記第２の時間期間のＨＦ・・・（ｃ５）
前記第１の時間期間のＬＦ／ＨＦ＜前記第２の時間期間のＬＦ／ＨＦ・・・（ｃ６）
但し、ＬＦは、拍動間隔を周波数スペクトル変換することにより取得されるパワースペクトルを周波数Ｌｆ１からＬｆ２まで定積分した値であり、ＨＦは、前記パワースペクトルを周波数Ｈｆ１からＨｆ２まで定積分した値であり、Ｈｆ１＞Ｌｆ１、Ｈｆ２＞Ｌｆ２である。
前記生体信号が、拍動を示す信号を含み、
前記第１の時間期間における生体信号と、前記第２の時間期間における生体信号との各々について、拍動間隔の平均値、単位時間当たりの拍数、拍動間隔のばらつき、単位時間当たりの拍数のばらつき、ＨＦ、およびＬＦ／ＨＦのうちの少なくとも一つを算出し、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されるか否かを判定することをさらに含み、
前記出力は、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されると判定された場合に行われる、
生体信号処理方法。
前記第１の時間期間における拍動間隔の平均値＞前記第２の時間期間における拍動間隔の平均値・・・（ｃ１）
前記第１の時間期間での単位時間当たりの拍数＜前記第２の時間期間での単位時間当たりの拍数・・・（ｃ２）
前記第１の時間期間における拍動間隔のばらつき＞前記第２の時間期間における拍動間隔のばらつき・・・（ｃ３）
第１の時間期間での単位時間当たりの拍数のばらつき＞第２の時間期間での単位時間当たりの拍数のばらつき・・・（ｃ４）
前記第１の時間期間のＨＦ＞前記第２の時間期間のＨＦ・・・（ｃ５）
前記第１の時間期間のＬＦ／ＨＦ＜前記第２の時間期間のＬＦ／ＨＦ・・・（ｃ６）
但し、ＬＦは、拍動間隔を周波数スペクトル変換することにより取得されるパワースペクトルを周波数Ｌｆ１からＬｆ２まで定積分した値であり、ＨＦは、前記パワースペクトルを周波数Ｈｆ１からＨｆ２まで定積分した値であり、Ｈｆ１＞Ｌｆ１、Ｈｆ２＞Ｌｆ２である。
【請求項8】
前記拍動間隔の平均値は、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの平均値である、請求項７に記載の生体信号処理方法。
【請求項9】
前記拍動間隔のばらつきは、ＳＤＮＮ、ＲＭＳＳＤ、ＣＶＲＲ、ＮＮ５０、およびｐＮＮ５０のいずれかであり、
ＳＤＮＮは、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの標準偏差であり、
ＲＭＳＳＤは、連続して隣接するＲＲＩの差の２乗の平均値の平方根であり、
ＣＶＲＲは、ＳＤＮＮの値をＲＲＩの平均値で割って１００を掛けた値であり、
ＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える総数を示す値であり、
ｐＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える心拍の割合を示す値である、
請求項７に記載の生体信号処理方法。
ＳＤＮＮは、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの標準偏差であり、
ＲＭＳＳＤは、連続して隣接するＲＲＩの差の２乗の平均値の平方根であり、
ＣＶＲＲは、ＳＤＮＮの値をＲＲＩの平均値で割って１００を掛けた値であり、
ＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える総数を示す値であり、
ｐＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える心拍の割合を示す値である、
請求項７に記載の生体信号処理方法。
【請求項10】
前記生体信号の取得は生体信号計測用衣類を用いて行われる、請求項１から９のいずれかに記載の生体信号処理方法。
【請求項11】
四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得する取得部と、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する第１時間期間特定部と、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行う出力部と
を備える、四肢動物の生体信号処理装置。
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する第１時間期間特定部と、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行う出力部と
を備える、四肢動物の生体信号処理装置。
【請求項12】
請求項１１に記載の生体信号処理装置と、
前記生体信号および／または前記加速度の取得に用いられる計測装置、ならびに、前記生体信号処理装置による出力に基づいて出力を行う出力装置、の少なくとも一方と
を具備する、生体信号処理システム。
前記生体信号および／または前記加速度の取得に用いられる計測装置、ならびに、前記生体信号処理装置による出力に基づいて出力を行う出力装置、の少なくとも一方と
を具備する、生体信号処理システム。
【請求項13】
請求項１１に記載の生体信号処理装置の各部による処理をコンピュータに実行させるためのプログラム。
【請求項14】
請求項１１に記載の生体信号処理装置の各部による処理をコンピュータに実行させるためのプログラムを記憶する、非一時的なコンピュータ読み取り可能な記憶媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、四肢動物の生体信号処理方法、生体信号処理装置、生体信号処理システム、プログラム、および記憶媒体に関する。
【背景技術】
【0002】
近年、動物の生体信号を測定するための技術が注目されている。例えば特許文献１では、分娩を監視すべき動物母体の心拍数を心拍計により測定し、測定された心拍数に基づいて、分娩以前の平時での該心拍数の変動域を超えた一過性の心拍数の減少を弁別して分娩直前信号を通知する技術が開示されている。
【先行技術文献】
【特許文献】
【0003】
【特許文献1】特許第４４８７０７５号公報
【発明の概要】
【発明が解決しようとする課題】
【0004】
しかしながら、生理反応として、心拍数は動作が少ないと少なくなり、動作が多くなると多くなる。したがって、監視のために心拍数を用いる場合、動物の動作を考慮しなければ誤った通知が行われてしまう。さらに、動物が立位である場合には、臥位である場合と比較して当該動物の心拍数が増加する。他の生体信号においても、動物の動作および姿勢の影響は同様に見られるものと考えられる。そこで、本発明は、四肢動物の動きに伴う加速度から、当該四肢動物が臥位であると推測される時間期間を特定し、当該時間期間における生体信号の処理を行うことを可能とする、四肢動物の生体信号処理方法、生体信号処理装置、生体信号処理システム、プログラム、および記憶媒体を提供することを目的とする。
【課題を解決するための手段】
【0005】
すなわち本発明は、以下の発明を含み得る。
［１］ コンピュータおよび記憶媒体を備える装置が実行する四肢動物の生体信号処理方法であって、
四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得することと、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定することと、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行うことと
を含む、生体信号処理方法。
［２］ 前記取得される加速度の値は、互いに交わるＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度ｘ、ｙ、およびｚのベクトル和の大きさに対応する、［１］に記載の生体信号処理方法。
［３］ 前記加速度の値のばらつきの度合いを示す値は、標準偏差、変動係数、および分散のいずれかである、［１］または［２］に記載の生体信号処理方法。
［４］ 前記第１の時間期間は、前記ばらつきの度合いを示す値が前記第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間である、［１］から［３］のいずれかに記載の生体信号処理方法。
［５］ 前記第１の時間期間より後の、前記ばらつきの度合いを示す値が第２の閾値以下の第２の時間期間を特定することをさらに含み、
前記出力は、前記生体信号のうちの前記第２の時間期間における生体信号にさらに基づいて行われる、
［１］から［４］のいずれかに記載の生体信号処理方法。
［６］ 前記第１の時間期間は、前記ばらつきの度合いを示す値が前記第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間であり、
前記第２の時間期間は、前記ばらつきの度合いを示す値が前記第２の閾値以下となる時点から第２の時間経過した時点以降に開始される時間期間である、
［５］に記載の四肢動物の生体信号処理方法。
［７］ ［５］または［６］に記載の生体信号処理方法であって、
前記生体信号が、拍動を示す信号を含み、
前記第１の時間期間における生体信号と、前記第２の時間期間における生体信号との各々について、拍動間隔の平均値、単位時間当たりの拍数、拍動間隔のばらつき、単位時間当たりの拍数のばらつき、ＨＦ、およびＬＦ／ＨＦのうちの少なくとも一つを算出し、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されるか否かを判定することをさらに含み、
前記出力は、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されると判定された場合に行われる、
生体信号処理方法。
前記第１の時間期間における拍動間隔の平均値＞前記第２の時間期間における拍動間隔の平均値・・・（ｃ１）
前記第１の時間期間での単位時間当たりの拍数＜前記第２の時間期間での単位時間当たりの拍数・・・（ｃ２）
前記第１の時間期間における拍動間隔のばらつき＞前記第２の時間期間における拍動間隔のばらつき・・・（ｃ３）
第１の時間期間での単位時間当たりの拍数のばらつき＞第２の時間期間での単位時間当たりの拍数のばらつき・・・（ｃ４）
前記第１の時間期間のＨＦ＞前記第２の時間期間のＨＦ・・・（ｃ５）
前記第１の時間期間のＬＦ／ＨＦ＜前記第２の時間期間のＬＦ／ＨＦ・・・（ｃ６）
但し、ＬＦは、拍動間隔を周波数スペクトル変換することにより取得されるパワースペクトルを周波数Ｌｆ１からＬｆ２まで定積分した値であり、ＨＦは、前記パワースペクトルを周波数Ｈｆ１からＨｆ２まで定積分した値であり、Ｈｆ１＞Ｌｆ１、Ｈｆ２＞Ｌｆ２である。
［８］ 前記拍動間隔の平均値は、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの平均値である、［７］に記載の生体信号処理方法。
［９］ 前記拍動間隔のばらつきは、ＳＤＮＮ、ＲＭＳＳＤ、ＣＶＲＲ、ＮＮ５０、およびｐＮＮ５０のいずれかであり、
ＳＤＮＮは、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの標準偏差であり、
ＲＭＳＳＤは、連続して隣接するＲＲＩの差の２乗の平均値の平方根であり、
ＣＶＲＲは、ＳＤＮＮの値をＲＲＩの平均値で割って１００を掛けた値であり、
ＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える総数を示す値であり、
ｐＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える心拍の割合を示す値である、
［７］に記載の生体信号処理方法。
［１０］ 前記生体信号の取得は生体信号計測用衣類を用いて行われる、［１］から［９］のいずれかに記載の生体信号処理方法。
［１１］ 四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得する取得部と、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する第１時間期間特定部と、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行う出力部と
を備える、四肢動物の生体信号処理装置。
［１２］ ［１１］に記載の生体信号処理装置と、
前記生体信号および／または前記加速度の取得に用いられる計測装置、ならびに、前記生体信号処理装置による出力に基づいて出力を行う出力装置、の少なくとも一方と
を具備する、生体信号処理システム。
［１３］ ［１１］に記載の生体信号処理装置の各部による処理をコンピュータに実行させるためのプログラム。
［１４］ ［１１］に記載の生体信号処理装置の各部による処理をコンピュータに実行させるためのプログラムを記憶する、非一時的なコンピュータ読み取り可能な記憶媒体。
［１］ コンピュータおよび記憶媒体を備える装置が実行する四肢動物の生体信号処理方法であって、
四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得することと、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定することと、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行うことと
を含む、生体信号処理方法。
［２］ 前記取得される加速度の値は、互いに交わるＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度ｘ、ｙ、およびｚのベクトル和の大きさに対応する、［１］に記載の生体信号処理方法。
［３］ 前記加速度の値のばらつきの度合いを示す値は、標準偏差、変動係数、および分散のいずれかである、［１］または［２］に記載の生体信号処理方法。
［４］ 前記第１の時間期間は、前記ばらつきの度合いを示す値が前記第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間である、［１］から［３］のいずれかに記載の生体信号処理方法。
［５］ 前記第１の時間期間より後の、前記ばらつきの度合いを示す値が第２の閾値以下の第２の時間期間を特定することをさらに含み、
前記出力は、前記生体信号のうちの前記第２の時間期間における生体信号にさらに基づいて行われる、
［１］から［４］のいずれかに記載の生体信号処理方法。
［６］ 前記第１の時間期間は、前記ばらつきの度合いを示す値が前記第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間であり、
前記第２の時間期間は、前記ばらつきの度合いを示す値が前記第２の閾値以下となる時点から第２の時間経過した時点以降に開始される時間期間である、
［５］に記載の四肢動物の生体信号処理方法。
［７］ ［５］または［６］に記載の生体信号処理方法であって、
前記生体信号が、拍動を示す信号を含み、
前記第１の時間期間における生体信号と、前記第２の時間期間における生体信号との各々について、拍動間隔の平均値、単位時間当たりの拍数、拍動間隔のばらつき、単位時間当たりの拍数のばらつき、ＨＦ、およびＬＦ／ＨＦのうちの少なくとも一つを算出し、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されるか否かを判定することをさらに含み、
前記出力は、下記式（ｃ１）～式（ｃ６）のうち少なくとも一つの式の大小関係が満足されると判定された場合に行われる、
生体信号処理方法。
前記第１の時間期間における拍動間隔の平均値＞前記第２の時間期間における拍動間隔の平均値・・・（ｃ１）
前記第１の時間期間での単位時間当たりの拍数＜前記第２の時間期間での単位時間当たりの拍数・・・（ｃ２）
前記第１の時間期間における拍動間隔のばらつき＞前記第２の時間期間における拍動間隔のばらつき・・・（ｃ３）
第１の時間期間での単位時間当たりの拍数のばらつき＞第２の時間期間での単位時間当たりの拍数のばらつき・・・（ｃ４）
前記第１の時間期間のＨＦ＞前記第２の時間期間のＨＦ・・・（ｃ５）
前記第１の時間期間のＬＦ／ＨＦ＜前記第２の時間期間のＬＦ／ＨＦ・・・（ｃ６）
但し、ＬＦは、拍動間隔を周波数スペクトル変換することにより取得されるパワースペクトルを周波数Ｌｆ１からＬｆ２まで定積分した値であり、ＨＦは、前記パワースペクトルを周波数Ｈｆ１からＨｆ２まで定積分した値であり、Ｈｆ１＞Ｌｆ１、Ｈｆ２＞Ｌｆ２である。
［８］ 前記拍動間隔の平均値は、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの平均値である、［７］に記載の生体信号処理方法。
［９］ 前記拍動間隔のばらつきは、ＳＤＮＮ、ＲＭＳＳＤ、ＣＶＲＲ、ＮＮ５０、およびｐＮＮ５０のいずれかであり、
ＳＤＮＮは、心電信号において隣り合うＲ波とＲ波との間隔であるＲＲＩの標準偏差であり、
ＲＭＳＳＤは、連続して隣接するＲＲＩの差の２乗の平均値の平方根であり、
ＣＶＲＲは、ＳＤＮＮの値をＲＲＩの平均値で割って１００を掛けた値であり、
ＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える総数を示す値であり、
ｐＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える心拍の割合を示す値である、
［７］に記載の生体信号処理方法。
［１０］ 前記生体信号の取得は生体信号計測用衣類を用いて行われる、［１］から［９］のいずれかに記載の生体信号処理方法。
［１１］ 四肢動物の生体信号と、前記四肢動物の動きに伴う加速度とを取得する取得部と、
前記加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する第１時間期間特定部と、
前記生体信号のうちの前記第１の時間期間における生体信号に基づいて出力を行う出力部と
を備える、四肢動物の生体信号処理装置。
［１２］ ［１１］に記載の生体信号処理装置と、
前記生体信号および／または前記加速度の取得に用いられる計測装置、ならびに、前記生体信号処理装置による出力に基づいて出力を行う出力装置、の少なくとも一方と
を具備する、生体信号処理システム。
［１３］ ［１１］に記載の生体信号処理装置の各部による処理をコンピュータに実行させるためのプログラム。
［１４］ ［１１］に記載の生体信号処理装置の各部による処理をコンピュータに実行させるためのプログラムを記憶する、非一時的なコンピュータ読み取り可能な記憶媒体。
【発明の効果】
【0006】
本発明によれば、四肢動物の動きに伴う加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間が特定され、当該第１の時間期間における生体信号に基づく出力が可能となる。特定される第１の時間期間では、当該四肢動物が臥位であると推測される。このため、四肢動物の動作および姿勢が当該四肢動物の生体信号に影響を及ぼしてしまうような時間期間の生体信号を排除して、当該四肢動物の生体信号に基づく出力を行うことが可能となる。
【図面の簡単な説明】
【0007】
【発明を実施するための形態】
【0008】
以下、図面を参照してこの発明に係る実施形態を説明する。
以下では、下記実施形態に基づき本発明をより具体的に説明するが、本発明はもとより下記実施形態によって制限を受けるものではなく、前・後記の趣旨に適合し得る範囲で適当に変更を加えて実施することも勿論可能であり、それらはいずれも本発明の技術的範囲に包含される。
以下では、下記実施形態に基づき本発明をより具体的に説明するが、本発明はもとより下記実施形態によって制限を受けるものではなく、前・後記の趣旨に適合し得る範囲で適当に変更を加えて実施することも勿論可能であり、それらはいずれも本発明の技術的範囲に包含される。
【0009】
［第１実施形態］
（構成例）
（１）システム構成
図１は、第１実施形態に係る生体信号処理装置１を含む生体信号処理システムＳＹＳの構成の一例を示す。当該システムＳＹＳは、生体信号処理装置１に加えて、例えば、計測装置２および出力装置３を含む。
（構成例）
（１）システム構成
図１は、第１実施形態に係る生体信号処理装置１を含む生体信号処理システムＳＹＳの構成の一例を示す。当該システムＳＹＳは、生体信号処理装置１に加えて、例えば、計測装置２および出力装置３を含む。
【0010】
計測装置２は、四肢動物に取り付けられる。計測装置２としては、加速度計に加え、生体信号の検出に用いられる心電計、心拍センサ、脈波センサ、呼吸数計測器、温度センサ、ならびに、電磁誘導、赤外線、または超音波の受信器、発振器等が挙げられる。呼吸数計測器は、例えば、呼吸により周長が変化する部位の長さの変化を検出する変位センサである。変位センサは、例えば、当該長さの変位により変化する抵抗または静電容量を検出するものである。図１では、計測装置２を便宜的に単一のブロックで表しているが、計測装置２として、例えば、加速度計と、上述したような生体信号の検出に用いられる各種機器との組み合わせが用いられる。当該生体信号の検出に用いられる各種機器としては、１種または２種以上が用いられてもよい。当該生体信号の検出に用いられる各種機器は、例えば、検出される電圧信号等の生体信号をアナログデジタル変換し、当該アナログデジタル変換後の信号を、通信ネットワークＮＷを介して外部装置に送信可能である。
【0011】
生体信号処理装置１は、例えば計測機器およびパーソナルコンピュータ等である。出力装置３は、例えば、スマートフォン、移動端末（例えばタブレット端末）、スマートウォッチ、ヘッドマウントディスプレイ、スマートグラス、スマートコンタクトレンズ、およびパーソナルコンピュータ等のような、音声出力、振動、および／または画面表示等の出力が可能な装置である。図１では、生体信号処理装置１と出力装置３とが通信ネットワークＮＷを介して接続される別個の装置であるものとして示されているが、本実施形態はこれに限定されない。生体信号処理装置１と出力装置３とを組み合わせたものが、単一の生体信号処理装置を構成していてもよい。例えば、生体信号処理装置１と出力装置３との組み合わせが、計測機器、パーソナルコンピュータ、スマートフォン、移動端末（例えばタブレット端末）、スマートウォッチ、およびヘッドマウントディスプレイ等を構成していてもよい。同様に、計測装置２、生体信号処理装置１、および出力装置３のうちの２つ以上が、単一の生体信号処理装置を構成していてもよい。
【0012】
計測装置２は、通信ネットワークＮＷを介して、計測装置２が取り付けられた四肢動物の生体信号および加速度データを生体信号処理装置１に送信する。生体信号処理装置１は、当該生体信号および加速度データを受信する。当該加速度データの受信により、生体信号処理装置１は、上記四肢動物の動きに伴う加速度を取得する。これにより、生体信号処理装置１は、当該加速度に基づいて当該生体信号の処理を行う。
【0013】
生体信号処理装置１による加速度に基づく生体信号の処理について説明する。生体信号処理装置１は、後述するように、時間期間特定部１１２による制御の下、四肢動物の動きに伴う加速度の値のばらつきの度合いを示す値が閾値以下の時間期間を特定する。当該時間期間では、四肢動物が臥位であると推測される。生体信号処理装置１は、生体信号処理部１１３による制御の下、例えば、当該時間期間における生体信号を処理する。生体信号処理装置１は、当該処理の結果に基づくデータを、通信ネットワークＮＷを介して出力装置３に送信可能である。
【0014】
出力装置３は、上記生体信号の処理の結果に基づくデータを受信し、当該データに基づいて音声出力および／または画面表示等の出力を行う。これにより、例えば、四肢動物の分娩接近、および分娩開始等をユーザが知覚することが可能となる。図１では、出力装置３を便宜的に単一のブロックで表しているが、出力装置３として、上述したような装置のうちの１種のみが用いられてもよく、２種以上が用いられてもよい。
【0015】
（２）計測装置の構成例
計測装置２の一例には、生体信号計測用衣類がある。
当該生体信号計測用衣類は、ベルト等の帯状物や衣服等であることが好ましい。当該衣類は、胸部、腹部、脚部、頸部、および頭部から選ばれる１以上の部位を覆うものであることが好ましく、胸部および／または腹部を覆うものであることがより好ましく、胸部を覆うものであることがさらに好ましい。当該衣類は、洗濯可能であり、複数回使用する観点から、着脱可能な電子ユニットを備えていることが好ましい。
計測装置２の一例には、生体信号計測用衣類がある。
当該生体信号計測用衣類は、ベルト等の帯状物や衣服等であることが好ましい。当該衣類は、胸部、腹部、脚部、頸部、および頭部から選ばれる１以上の部位を覆うものであることが好ましく、胸部および／または腹部を覆うものであることがより好ましく、胸部を覆うものであることがさらに好ましい。当該衣類は、洗濯可能であり、複数回使用する観点から、着脱可能な電子ユニットを備えていることが好ましい。
【0016】
生体信号計測用衣類は、生地と、当該生地の体表面側に設けられた電極と、配線と、電子機器接続部と、生体信号を計測する機器とを少なくとも備え、当該生地は、体周方向の一端部側に第１帯状部と、体周方向の他端部側に第２帯状部とを備えることが好ましい。
【0017】
生体信号計測用衣類は、少なくとも電極が形成されている領域において帯状本体部１ｃｍ２当たりに対して該帯状本体部の厚み方向に５ｋｇの荷重をかけたときに、厚みが１０ｍｍ以上小さくなる四肢動物用の衣類であることが好ましい。電極が形成されている領域にクッション性があることで、四肢動物に着用させた場合に、電極の位置ずれが生じ難くなるため生体信号の測定精度が向上する。
【0018】
第１帯状部と第２帯状部は、四肢動物に生体信号計測用衣類を着用させる際に直接または間接に、互いに固定等できるようになっていればよい。例えば、第１帯状部と第２帯状部は、それぞれ第１連結部材、第２連結部材とを有することが好ましい。第１連結部材、第２連結部材として、マジックテープ（登録商標）、フリーマジックテープ（登録商標）等の面ファスナーや、バックル、フック、およびループ等が挙げられる。また、第１連結部材と第２連結部材は、例えば取り外し式ベルトを介して連結されていてもよい。これにより体周方向の長さを調節し易くすることができるため、生体信号計測用衣類の装着を容易とすることができる。
【0019】
生地は、第１帯状部および第２帯状部に加えて本体部を有しており、本体部では、第１帯状部と第２帯状部よりも前後方向の長さが大きい。本体部は、例えば第１帯状部に向かって前後方向の長さが短くなる第１テーパ部を有していてもよい。本体部は、第２帯状部に向かって前後方向の長さが短くなる第２テーパ部を有していてもよい。本体部が第１テーパ部および／または第２テーパ部を有することにより、生体信号計測用衣類の締め付けが容易となる。第１テーパ部および第２テーパ部の各々の形状として、外縁が後方向に向かって湾曲している湾曲形状、および、一定の角度で先細りになっているいわゆる線形テーパ形状等が挙げられる。本体部は、正方形、長方形等の矩形であってもよい。
【0020】
生地の体表面側には、例えば、第一絶縁層、導電層、および第二絶縁層が、当該生地から見て第一絶縁層、導電層、および第二絶縁層の順に設けられている。第二絶縁層は、導電層の一部が露出するように設けられている。
【0021】
電極は、第一絶縁層、導電層、および第二導電層を含む積層体のうち、導電層が露出している第一絶縁層および導電層を含む積層構造として定義できる。電極は、１つに限らず２つ以上、生地の体表面側に設けられていてもよい。
【0022】
生体信号計測用衣類がベルト状の場合、当該ベルト状の生体信号計測用衣類の体表面側は平均摩擦係数ＭＩＵが０．４０以上である高摩擦部を有していることが好ましい。これによりベルト状の生体信号計測用衣類がずれ難くなる。当該衣類の生地の体表面側が高摩擦部を有していてもよく、第１帯状部、第２帯状部の体表面側が高摩擦部を有していてもよい。
【0023】
高摩擦部は、細繊維が露出している部分であることが好ましい。細繊維としては、マイクロファイバー、ナノファイバーが好ましく挙げられる。細繊維の単繊維径は、好ましくは８μｍ以下、より好ましくは１０００ｎｍ以下、さらに好ましくは１０００ｎｍ未満、さらにより好ましくは８００ｎｍ以下、特に好ましくは７５０ｎｍ以下である。一方、単繊維径は、１００ｎｍ以上であってもよく、３００ｎｍ以上であってもよく、５００ｎｍ以上であってもよい。
【0024】
ベルト状の生体信号計測用衣類の体表面側はさらに、平均摩擦係数ＭＩＵが０．４０未満である低摩擦部を有していることが好ましい。これによりベルト状の生体信号計測用衣類の位置の調節が容易となる。当該衣類の生地の体表面側が低摩擦部を有していてもよく、第１帯状部、第２帯状部の体表面側が低摩擦部を有していてもよい。
【0025】
ベルト状の生体信号計測用衣類は、５％伸長時の引張強さが２．８Ｎ／ｃｍ以下であり、１０％伸長時の引張強さが４．０Ｎ／ｃｍ以下であることが好ましい。引張強さが上記範囲であることにより、ベルトが伸長することができる。引張強さは、例えば以下の手順で測定することができる。先ず、引張試験機のチャック間の中心に生地の中心が位置し、且つチャック間距離が４５ｃｍとなるようにベルト状の生体信号計測用衣類をチャックで挟み、次いでロードセル１ｋＮ、引張速度１００ｍｍ／分の条件で、体周方向に伸長させて、次いで２．２５ｃｍ伸長時（５％伸長時）と４．５ｃｍ伸長時（１０％伸長時）における荷重（Ｎ）をそれぞれ求める。さらに上記チャック間におけるベルト状の生体信号計測用衣類の平均幅で、荷重（Ｎ）を割って、５％伸長時の引張強さ（Ｎ／ｃｍ）、１０％伸長時の引張強さ（Ｎ／ｃｍ）をそれぞれ求めることができる。
【0026】
生地の体表面側には、配線が設けられていることが好ましい。配線は、第一絶縁層、導電層、および第二絶縁層を含む積層体のうち、導電層が第二絶縁層に覆われている部分として定義できる。電極で取得された生体信号を、配線を介して伝達することができる。生体信号を伝達するための配線として、上記で定義したものに加え／上記で定義したものの代わりに、リード線が用いられるようにしてもよい。
【0027】
体表面側と第一絶縁層の間に接着層等の他の層が存在していてもよい。第一絶縁層と導電層との間にホットメルト層等の他の層が存在していてもよく、導電層と第二絶縁層との間にホットメルト層等の他の層が存在していてもよい。
【0028】
電子機器接続部は、第一絶縁層、導電層、および第二絶縁層を含む積層体のうち、導電層が第二絶縁層に覆われていない第一絶縁層および第一絶縁層を含む積層体として定義できる。電子機器接続部において、例えばスナップファスナー等のコネクタを介して、生地の体表面側とは反対側の表側面に電子機器を取り付けることができる。また、電子機器接続部の体表面側を覆うように保護層が形成されていてもよい。
【0029】
これら電極、第一絶縁層、導電層、第二絶縁層、配線等の詳細については、特開２０２０－１００９０３号公報を参照することができる。
【0030】
（３）生体信号処理装置に係るハードウェア構成
図２は、第１実施形態に係る生体信号処理装置１のハードウェア構成の一例を示すブロック図である。
図２は、第１実施形態に係る生体信号処理装置１のハードウェア構成の一例を示すブロック図である。
【0031】
生体信号処理装置１は、制御部１１、プログラム記憶部１２、データ記憶部１３、入出力インタフェース（入出力Ｉ／Ｆ）１４、およびバスＢＵＳを含む。プログラム記憶部１２、データ記憶部１３、および入出力インタフェース１４の各々は、バスＢＵＳを介して制御部１１に接続される。
【0032】
制御部１１は、例えば、中央処理ユニット（Central Processing Unit：ＣＰＵ）等のハードウェアプロセッサを有する。
【0033】
プログラム記憶部１２は、記憶媒体として、例えば、ＨＤＤ（Hard Disk Drive）またはＳＳＤ（Solid State Drive）等の随時書込みおよび読出しが可能な不揮発性メモリと、ＲＯＭ（Read Only Memory）等の不揮発性メモリとを組み合わせたものである。プログラム記憶部１２は、ＯＳ（Operating System）等のミドルウェアに加えて、本実施形態に係る各種制御処理の実行のために用いられるプログラムを格納する。プログラム記憶部１２は、例えば、フラッシュメモリ等の補助記憶装置、ＣＤ－ＲＯＭ、ＤＶＤディスク、およびＵＳＢメモリ等の可搬型の記録媒体により実現されてもよい。
【0034】
データ記憶部１３は、記憶媒体として、例えば、ＨＤＤまたはＳＳＤ等の随時書込みおよび読出しが可能な不揮発性メモリと、ＲＡＭ（Random Access Memory）等の揮発性メモリとを組み合わせたものである。データ記憶部１３は、制御部１１が有するハードウェアプロセッサの作業領域として使用され、データを一時的に保持し、バッファおよびキャッシュとして機能する。データ記憶部１３は、例えば、フラッシュメモリ等の補助記憶装置、ＣＤ－ＲＯＭ、ＤＶＤディスク、およびＵＳＢメモリ等の可搬型の記録媒体により実現されてもよい。
【0035】
入出力インタフェース１４は、制御部１による制御の下、通信ネットワークＮＷにより定義される通信プロトコルを使用して、外部装置との間で伝送されるデータの送受信を行う。入出力インタフェース１４は、例えば有線ＬＡＮまたは無線ＬＡＮに対応するインタフェースにより構成される。
【0036】
（４）生体信号処理装置に係るソフトウェア構成
図３は、第１実施形態に係る生体信号処理装置１の制御部１１のソフトウェア構成の一例を示すブロック図である。制御部１１中に示される各種機能部は、一例にすぎず、図３に示されるように区別されている必要はない。単一のブロックとして示されている或る機能部が複数の単位に分割されて実現されてもよいし、異なる複数のブロックとして示されている機能部が或るまとまった単位で実現されてもよい。
図３は、第１実施形態に係る生体信号処理装置１の制御部１１のソフトウェア構成の一例を示すブロック図である。制御部１１中に示される各種機能部は、一例にすぎず、図３に示されるように区別されている必要はない。単一のブロックとして示されている或る機能部が複数の単位に分割されて実現されてもよいし、異なる複数のブロックとして示されている機能部が或るまとまった単位で実現されてもよい。
【0037】
制御部１１は、例えば、取得部１１１、時間期間特定部１１２、生体信号処理部１１３、および出力部１１４を含む。制御部１１が含む各機能部の処理機能は、制御部１１が、プログラム記憶部１２に格納されるプログラムを、制御部１１のハードウェアプロセッサに実行させることにより実現される。プログラム記憶部１２に格納されるプログラムが用いられる場合の説明を行ったが、用いられるプログラムは、通信ネットワークＮＷを通して提供されるものであってもよい。
【0038】
入出力インタフェース１４は、計測装置２から通信ネットワークＮＷを介して送信される生体信号および加速度データを受け取り、当該生体信号および加速度データを制御部１１に入力する。入出力インタフェース１４はまた、制御部１１から出力されるデータを受け取り、制御部１１からの指示にしたがって、当該データを出力装置３に通信ネットワークＮＷを介して送信する。
【0039】
データ記憶部１３は、例えば、加速度データ記憶部１３１、生体信号記憶部１３２、時間期間記憶部１３３、および生体信号処理結果記憶部１３４を含む。これら各種記憶部は、生体信号処理装置１に含まれているものとして説明するが、これら記憶部のうち１つ以上が、例えばクラウドコンピューティング上等、生体信号処理装置１の外部に設けられていてもよい。
【0040】
加速度データ記憶部１３１は、加速度データを記憶する。生体信号記憶部１３２は、生体信号を記憶する。時間期間記憶部１３３は、時間期間特定部１１２による処理で特定される時間期間を記憶する。生体信号処理結果記憶部１３４は、生体信号処理部１１３による処理の結果のデータを記憶する。
【0041】
取得部１１１は、計測装置２により計測される生体信号および加速度データを入出力インタフェース１４を介して取得し、当該加速度データを加速度データ記憶部１３１に記憶させ、当該生体信号を生体信号記憶部１３２に記憶させる処理を実行する。加速度データ記憶部１３１では、各加速度データアイテムがタイミング情報に対応付けられて記憶される。同様に、生体信号記憶部１３２では、生体信号がタイミング情報に対応付けられて記憶される。本明細書では、後述するように加速度データも生体信号の一例として説明される場合がある。この場合には、加速度データ記憶部１３１も生体信号記憶部１３２に含まれるものとみなすことができる。
【0042】
時間期間特定部１１２は、例えば、第１時間期間特定部１１２１および第２時間期間特定部１１２２を含む。
第１時間期間特定部１１２１は、加速度データ記憶部１３１から加速度データを読み出し、計測装置２が取り付けられた四肢動物の運動に伴う加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定し、当該第１の時間期間を示すデータを時間期間記憶部１３３に記憶させる処理を実行する。第１の時間期間では、例えば、当該四肢動物が臥位であると推測される。
第１時間期間特定部１１２１は、加速度データ記憶部１３１から加速度データを読み出し、計測装置２が取り付けられた四肢動物の運動に伴う加速度の値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定し、当該第１の時間期間を示すデータを時間期間記憶部１３３に記憶させる処理を実行する。第１の時間期間では、例えば、当該四肢動物が臥位であると推測される。
【0043】
第１の時間期間は、上記ばらつきの度合いを示す値が第１の閾値以下となる時点から第１の閾値以下である時間期間のうちの一部であっても全体であってもよい。例えば、上記ばらつきの度合いを示す値が第１の閾値以下となる時点から第１の閾値以下である時間期間のうちの、当該第１の閾値以下となる時点から第１の時間経過した時点以降に開始される時間期間が、第１の時間期間であってもよい。
【0044】
第２時間期間特定部１１２２は、上記ばらつきの度合いを示す値が第２の閾値以下の第２の時間期間を特定し、当該第２の時間期間を示すデータを時間期間記憶部１３３に記憶させる処理を実行する。第２の時間期間の始点は、例えば第１の時間期間の終点以降である。第２の時間期間においても、例えば、上記四肢動物が臥位であると推測される。第２の閾値は、第１の閾値と同一であっても異なっていてもよい。
【0045】
第２の時間期間は、上記ばらつきの度合いを示す値が第２の閾値以下となる時点から第２の閾値以下である時間期間のうちの一部であっても全体であってもよい。例えば、上記ばらつきの度合いを示す値が第２の閾値以下となる時点から第２の閾値以下である時間期間のうちの、当該第２の閾値以下となる時点から第２の時間経過した時点以降に開始される時間期間が、第２の時間期間であってもよい。第２の時間は、上記第１の時間と同一であっても異なっていてもよい。
【0046】
生体信号処理部１１３は、例えば、時間期間記憶部１３３に記憶される第１時間期間を示すデータを読み出し、さらに、第１時間期間のタイミング情報に対応付けられた生体信号のデータを生体信号記憶部１３２から読み出す処理を実行する。生体信号処理部１１３は、第１時間期間における生体信号を処理し、当該処理の結果のデータを生体信号処理結果記憶部１３４に記憶させる処理を実行してもよい。当該生体信号の処理では、例えば、後述する（ｉ１）拍動間隔の平均値、（ｉ２）単位時間当たりの拍数、（ｉ３）拍動間隔のばらつき、（ｉ４）単位時間当たりの拍数のばらつき、（ｉ５）ＨＦ、（ｉ６）ＬＦ／ＨＦ、（ｉ７）加速度、（ｉ８）角速度、（ｉ９）単位時間当たりの呼吸数、（ｉ１０）外陰部間の距離、（ｉ１１）体表温度、および（ｉ１２）膣内温度（以下、便宜的に、それぞれ指標（ｉ１）～（ｉ１２）とも称され得る。）のうちの少なくとも１つが算出され、生体信号処理結果記憶部１３４に記憶される。当該生体信号の処理では、指標（ｉ１）～（ｉ６）のうちの少なくとも１つが算出されてもよい。生体信号処理部１１３による処理では、適宜ノイズカット処理等が行われてもよい。ノイズカット処理は生体信号処理部１１３により行われるものに限定されない。
【0047】
生体信号処理部１１３は、時間期間記憶部１３３に記憶される第２時間期間を示すデータを読み出し、さらに、第２時間期間のタイミング情報に対応付けられた生体信号のデータを生体信号記憶部１３２から読み出す処理を実行してもよい。生体信号処理部１１３は、第２時間期間における生体信号を処理し、当該処理の結果のデータを生体信号処理結果記憶部１３４に記憶させる処理を実行してもよい。当該生体信号の処理では、上記で第１の時間期間における生体信号の処理について説明したのと同様の算出処理が行われてもよい。
【0048】
生体信号処理部１１３は、例えば判定部１１３１を含む。
判定部１１３１は、第１の時間期間および／または第２の時間期間における生体信号の処理の結果の上記データに基づいて例えば或る条件式が満足されるか否かを判定し、当該判定の結果のデータを生体信号処理結果記憶部１３４に記憶させる処理を実行する。判定部１１３１は、例えば、第１時間期間における生体信号の処理の結果のデータと、第２時間期間における生体信号の処理の結果のデータとが、１または２以上の条件式を満足するか否かを判定し、当該判定の結果のデータを生体信号処理結果記憶部１３４に記憶させる処理を実行してもよい。
判定部１１３１は、第１の時間期間および／または第２の時間期間における生体信号の処理の結果の上記データに基づいて例えば或る条件式が満足されるか否かを判定し、当該判定の結果のデータを生体信号処理結果記憶部１３４に記憶させる処理を実行する。判定部１１３１は、例えば、第１時間期間における生体信号の処理の結果のデータと、第２時間期間における生体信号の処理の結果のデータとが、１または２以上の条件式を満足するか否かを判定し、当該判定の結果のデータを生体信号処理結果記憶部１３４に記憶させる処理を実行してもよい。
【0049】
出力部１１４は、例えば、生体信号処理結果記憶部１３４に記憶される生体信号の処理の結果のデータを読み出し、当該処理の結果に基づくデータを、入出力インタフェース１４を介して出力装置３に出力する処理を実行する。出力部１１４は、例えば、時間期間記憶部１３３に記憶される第１時間期間および／または第２時間期間について、第１時間期間および／または第２時間期間のタイミング情報に対応付けられた生体信号のデータを生体信号記憶部１３２から読み出し、当該読み出したデータを、入出力インタフェース１４を介して出力装置３に出力する処理を実行してもよい。
【0050】
出力部１１４は、例えば、生体信号処理結果記憶部１３４に記憶される上記判定の結果のデータを読み出し、当該データに応じて通知信号を生成し、当該通知信号を、入出力インタフェース１４を介して出力装置３に出力する処理を実行する。当該出力により、出力装置３による例えば通知信号に基づく、警告等に係る音声出力、振動、および／または画面表示等の出力が可能となる。
【0051】
このようにして、出力部１１４は、生体信号のうちの第１の時間期間および／または第２の時間期間における生体信号に基づいて出力を行う処理を実行することが可能である。
【0052】
（動作例）
以上のように構成された生体信号処理装置１の動作例を説明する。
（１）全体動作フロー
図４は、生体信号処理装置１により実行される動作の一例のフローチャートを示す。以下で説明する動作は一例に過ぎず、本実施形態に係る動作はこれに限定されるものではない。
以上のように構成された生体信号処理装置１の動作例を説明する。
（１）全体動作フロー
図４は、生体信号処理装置１により実行される動作の一例のフローチャートを示す。以下で説明する動作は一例に過ぎず、本実施形態に係る動作はこれに限定されるものではない。
【0053】
当該動作に先立ち、計測装置２が、計測装置２が取り付けられた四肢動物の生体信号および加速度データを計測する。計測装置２は、（図示しない）出力部による制御の下、当該生体信号および加速度データを、通信ネットワークＮＷを介して生体信号処理装置１に送信する。当該送信に応じて、図４のフローチャートに示される動作が開始される。
【0054】
生体信号処理装置１の制御部１１は、取得部１１１による制御の下、四肢動物の加速度データおよび生体信号を取得する（ＳＴ０１）。ＳＴ０１の動作は、以降で説明する各ＳＴの動作が実行されている間に、その少なくとも一部が並行して実行されていてもよい。
【0055】
四肢動物は、家畜、または、愛玩動物等のペットであることが好ましい。より具体的には、四肢動物は、犬、猫、牛、馬、羊、豚、または山羊等であることがより好ましく、犬または牛であることがさらに好ましい。
【0056】
加速度データの計測には、四肢動物の体に取り付けられた加速度計が用いられる。当該加速度計は、例えば胴に取り付けられ、好ましくは胸腹部に取り付けられ、より好ましくは左側の胸腹部に取り付けられる。
【0057】
加速度計としては、一軸加速度計、二軸加速度計、および三軸加速度計等が用いられる。なかでも、二軸加速度計、三軸加速度計が用いられることが好ましく、三軸加速度計が用いられることがより好ましい。以下では、三軸加速度計が用いられる場合を例に挙げて説明するが、本実施形態はこれに限定されない。
加速度計として、圧電型加速度センサ、サーボ型加速度センサ、ひずみゲージ型加速度センサ、および半導体型加速度センサ等が使用されてもよく、圧電型加速度センサが使用されることが好ましい。加速度計としては、ユニオンツール株式会社製のｍｙＢｅａｔ（登録商標）等が挙げられる。
加速度計として、圧電型加速度センサ、サーボ型加速度センサ、ひずみゲージ型加速度センサ、および半導体型加速度センサ等が使用されてもよく、圧電型加速度センサが使用されることが好ましい。加速度計としては、ユニオンツール株式会社製のｍｙＢｅａｔ（登録商標）等が挙げられる。
【0058】
加速度計で取得される加速度は、例えば、四肢動物の動きに伴う加速度と四肢動物に作用する重力加速度との合成値である。当該合成値は、重力加速度の大きさｇ（＝９．８ｍ／ｓ２）に対する比で表されるものであってもよい（単位：無次元量）。
【0059】
四肢動物の動きに伴う加速度（以下、加速度Ａとも称され得る。）の値は、例えば、三軸加速度計を用いて計測されるＸ軸、Ｙ軸、およびＺ軸方向それぞれの加速度ｘ、ｙ、およびｚから算出される。Ｘ軸、Ｙ軸、およびＺ軸は互いに交わり、より好ましくは互いに直交する。加速度Ａは、例えば、次に説明するように加速度ｘ、ｙ、およびｚのベクトル和の大きさに対応する。
【0060】
より具体的には、加速度Ａは、次の式（１）で表されるように、加速度ｘ、ｙ、およびｚの大きさ（式（１）ではそれぞれｘ、ｙ、およびｚと示している。）それぞれを２乗した値の和（二乗和）の平方根から重力加速度の大きさｇ（＝９．８ｍ／ｓ2）分として（ｇ／ｇ）＝１を減じた値である。このような加速度Ａでは、四肢動物の動きに伴う加速度のみが反映されるようになり得る。
式（１）において、四肢動物の動きがない場合、加速度ｘ、ｙ、およびｚのベクトル和は重力加速度とほぼ同一であり、ゆえに加速度ｘ、ｙ、およびｚの二乗和の値がほぼ１となることから、加速度Ａは０に殆ど等しい値となる。他方、式（１）において、四肢動物の動きがある場合、加速度ｘ、ｙ、およびｚのベクトル和の大きさが重力加速度の大きさよりも大きくなり、ゆえに加速度ｘ、ｙ、およびｚの大きさの二乗和の値が１よりも大きくなることから、加速度Ａは０よりも大きくなる。
式（１）において、四肢動物の動きがない場合、加速度ｘ、ｙ、およびｚのベクトル和は重力加速度とほぼ同一であり、ゆえに加速度ｘ、ｙ、およびｚの二乗和の値がほぼ１となることから、加速度Ａは０に殆ど等しい値となる。他方、式（１）において、四肢動物の動きがある場合、加速度ｘ、ｙ、およびｚのベクトル和の大きさが重力加速度の大きさよりも大きくなり、ゆえに加速度ｘ、ｙ、およびｚの大きさの二乗和の値が１よりも大きくなることから、加速度Ａは０よりも大きくなる。
【0061】
【数1】
【0062】
加速度Ａは、例えば、三軸加速度計を用いて計測される加速度ｘ、ｙ、およびｚの大きさそれぞれを２乗した値の和（二乗和）の平方根であってもよい。
【0063】
加速度計または計測装置２により、加速度Ａが算出されてもよく、この場合、計測装置２は、加速度Ａを生体信号処理装置１に送信してもよい。あるいは、計測装置２が加速度ｘ、ｙ、およびｚを生体信号処理装置１に送信し、生体信号処理装置１の制御部１１により加速度Ａが算出されることにより加速度Ａが取得されてもよい。このように生体信号処理装置１に送信される加速度に係るデータを、本明細書では総称して加速度データと称する。
【0064】
生体信号の例には、計測装置２が取り付けられた四肢動物の（ｉ１）拍動間隔の平均値、（ｉ２）単位時間当たりの拍数、（ｉ３）拍動間隔のばらつき、（ｉ４）単位時間当たりの拍数のばらつき、（ｉ５）ＨＦ、（ｉ６）ＬＦ／ＨＦ、（ｉ７）加速度、（ｉ８）角速度、（ｉ９）単位時間当たりの呼吸数、（ｉ１０）外陰部間の距離、（ｉ１１）体表温度、および（ｉ１２）膣内温度（上述したように、それぞれ指標（ｉ１）～（ｉ１２）とも称され得る。）、あるいは、これらの算出に用いられる、例えば拍動を示す信号等の何らかの他の信号がある。加速度、角速度等、および、加速度、角速度等の算出に用いられる信号も、生体に関する情報が示されるという点で、本明細書では生体信号の一例として説明される場合がある。
【0065】
続いて、制御部１１は、時間期間特定部１１２による制御の下、加速度Ａの値のばらつきの度合いを示す値が閾値以下の時間期間を特定する（ＳＴ０２）。
【0066】
特定される時間期間は、例えば、当該ばらつきの度合いを示す値が当該閾値以下となる時点から当該閾値以下であり続ける時間が５分以上、好ましくは１０分以上、より好ましくは２０分以上、さらに好ましくは３０分以上、さらに好ましくは１時間以上である時間期間から選択される。このように或る程度の長さの時間期間が特定されるようにすることにより、四肢動物が瞬間的に静止している時間期間が特定されないようにすることが可能である。上記特定される時間期間は、例えば１分以上、２時間以下、好ましくは２分以上、１時間以下、より好ましくは３分以上、５０分以下、さらに好ましくは５分以上、３０分以下である。
【0067】
或るタイミングに対応する上記ばらつきの度合いを示す値としては、例えば、当該タイミングを含む或る特定の期間における加速度Ａの値から算出される（ｉ）分散、（ｉｉ）標準偏差、および（ｉｉｉ）変動係数のいずれかが用いられる。当該特定の期間は、例えば１分以上、６０分以下、好ましくは２分以上、３０分以下、より好ましくは３分以上、１５分以下、さらに好ましくは５分以上、１０分以下である。例えば、このような特定の期間毎に、加速度Ａの複数の値から、上記ばらつきの度合いを示す値が算出される。ばらつきの度合いを示す値の算出に用いられる特定の期間の長さは必ずしも一定でなくてもよい。
【0068】
上記特定される時間期間では、例えば、四肢動物が臥位であると推測される。例えば、上記ばらつきの度合いを示す値が上記閾値以下となる時点から当該閾値以下である時間期間のうちの、上記閾値以下となる時点から或る時間（以下、待機時間とも称され得る。）経過した時点以降に開始される時間期間が、特定される時間期間であってもよい。当該待機時間は、例えば１分以上、６０分以下、好ましくは２分以上、３０分以下、より好ましくは３分以上、１５分以下、さらに好ましくは５分以上、１０分以下である。これは、四肢動物が臥位となっても、しばらくは臥位となる以前の動作の影響が生体信号にみられることがあり、このような影響が見られる時間帯を、特定される時間期間から排除することを目的とする。
【0069】
（ｉ）分散
分散としては、例えば、次に説明する分散ｓ２と分散Ｖａｒ[Ｘ]のいずれかが用いられる。
分散ｓ２は、上記特定の期間における加速度の値のばらつきの度合いを示し、以下の式（２）で表される。式（２）中、ｎは当該特定の期間内で求められている加速度Ａの数を表し、ｘｉは当該特定の期間内で求められているｉ番目の加速度Ａの値である。すなわち、分散ｓ２は、当該特定の期間内で求められている或る加速度Ａの値と当該特定の期間内で求められている加速度Ａの値の平均との差を２乗した値を各加速度Ａについて算出して和をとり、ｎで割った値となる。
分散Ｖａｒ[Ｘ]は平均との差の二乗の期待値であるとした場合は、確率変数Ｘと平均μを用いて式（３）で表してもかまわない。
分散としては、例えば、次に説明する分散ｓ２と分散Ｖａｒ[Ｘ]のいずれかが用いられる。
分散ｓ２は、上記特定の期間における加速度の値のばらつきの度合いを示し、以下の式（２）で表される。式（２）中、ｎは当該特定の期間内で求められている加速度Ａの数を表し、ｘｉは当該特定の期間内で求められているｉ番目の加速度Ａの値である。すなわち、分散ｓ２は、当該特定の期間内で求められている或る加速度Ａの値と当該特定の期間内で求められている加速度Ａの値の平均との差を２乗した値を各加速度Ａについて算出して和をとり、ｎで割った値となる。
分散Ｖａｒ[Ｘ]は平均との差の二乗の期待値であるとした場合は、確率変数Ｘと平均μを用いて式（３）で表してもかまわない。
【0070】
【数2】
【0071】
【数3】
【0072】
（ｉｉ）標準偏差
標準偏差としては、例えば、次に説明する標準偏差ｓが用いられる。
標準偏差ｓは、上記特定の期間における加速度の値のばらつきの度合いを示し、具体的には以下の式（４）で表される。式（４）中、ｓ２は上述した分散ｓ２である。
標準偏差としては、例えば、次に説明する標準偏差ｓが用いられる。
標準偏差ｓは、上記特定の期間における加速度の値のばらつきの度合いを示し、具体的には以下の式（４）で表される。式（４）中、ｓ２は上述した分散ｓ２である。
【0073】
【数4】
【0074】
（ｉｉｉ）変動係数
変動係数としては、例えば、次に説明する変動係数ＣＶが用いられる。
変動係数ＣＶは、上記特定の期間における加速度の値のばらつきの度合いを示し、具体的には以下の式（５）で表される。式（５）中、ｓは上述した標準偏差ｓである。すなわち、変動係数ＣＶは、標準偏差ｓを、上記特定の期間内で求められている加速度Ａの値の平均で割った値となる。
変動係数としては、例えば、次に説明する変動係数ＣＶが用いられる。
変動係数ＣＶは、上記特定の期間における加速度の値のばらつきの度合いを示し、具体的には以下の式（５）で表される。式（５）中、ｓは上述した標準偏差ｓである。すなわち、変動係数ＣＶは、標準偏差ｓを、上記特定の期間内で求められている加速度Ａの値の平均で割った値となる。
【0075】
【数5】
【0076】
加速度Ａの値のばらつきの度合いを示す値として、標準偏差ｓが用いられる場合、上記閾値は、例えば０．１以下、好ましくは０．０６以下、より好ましくは０．０３以下である。当該閾値は、好ましくは、０．００１以上である。これらの数値範囲は、上述したように加速度Ａを無次元量とした場合のものである。
【0077】
加速度Ａの値のばらつきの度合いを示す値として、分散ｓ２または分散Ｖａｒ[Ｘ]が用いられる場合、上記閾値は、例えば０．０１以下、好ましくは０．００３６以下、より好ましくは０．０００９以下である。当該閾値は、好ましくは、０．０００００１以上である。これらの数値範囲は、上述したように加速度Ａを無次元量とした場合のものである。
【0078】
加速度Ａの値のばらつきの度合いを示す値として、変動係数ＣＶが用いられる場合、上記閾値は、例えば０．１以下、好ましくは０．０６以下、より好ましくは０．０３以下である。当該閾値は、好ましくは、０．００１以上である。これらの数値範囲は、上述したように加速度Ａを無次元量とした場合のものである。
【0079】
続いて、制御部１１は、生体信号処理部１１３および／または出力部１１４による制御の下、上記特定された時間期間における生体信号に基づいて出力を行う（ＳＴ０３）。より具体的には、制御部１１は、生体信号処理部１１３による制御の下、上記特定された時間期間における生体信号を処理して、当該処理の結果のデータを出力してもよい。制御部１１は、当該処理の結果に基づいて或る条件式が満足されるか否かを判定し、当該判定の結果のデータに基づく通知信号を生成し、当該通知信号を出力してもよい。
【0080】
図５は、ＳＴ０２の時間期間を特定する動作を説明するための例示的なグラフを示す。当該グラフにおいて、縦軸は、上記ばらつきの度合いを示す値として用いられる標準偏差ｓの大きさを示し、横軸は時刻を示す。
加速度Ａは、三軸加速度計を用いて計測される加速度ｘ、ｙ、およびｚの大きさそれぞれを２乗した値の和（二乗和）の平方根として無次元量で算出されており、標準偏差ｓは上記特定の期間としての５分間毎に算出されている。加速度ｘ、ｙ、およびｚは、互いに直交するＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度である。当該グラフに示されるように、標準偏差ｓは５分毎にプロットされている。
加速度Ａは、三軸加速度計を用いて計測される加速度ｘ、ｙ、およびｚの大きさそれぞれを２乗した値の和（二乗和）の平方根として無次元量で算出されており、標準偏差ｓは上記特定の期間としての５分間毎に算出されている。加速度ｘ、ｙ、およびｚは、互いに直交するＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度である。当該グラフに示されるように、標準偏差ｓは５分毎にプロットされている。
【0081】
図５では、上記閾値（Ｔｈとして示されている。）として０．０３が用いられる場合の例が示される。図５の例では、上記ばらつきの度合いを示す値が閾値としての０．０３以下となる時点から当該閾値以下であり続ける時間期間として、時間期間ＴＰ１、時間期間ＴＰ２、時間期間ＴＰ３、時間期間ＴＰ４、時間期間ＴＰ５、時間期間ＴＰ６、および時間期間ＴＰ７が示される。例えば、ＳＴ０２の動作において特定される時間期間が、上記ばらつきの度合いを示す値が上記閾値以下となる時点から当該閾値以下であり続ける時間が１０分以上である時間期間から選択される場合、特定される時間期間は、時間期間ＴＰ１、時間期間ＴＰ３、時間期間ＴＰ４、時間期間ＴＰ５、時間期間ＴＰ６、および時間期間ＴＰ７から選択され得る一方で、時間期間ＴＰ２からは選択されない。ＳＴ０２の動作では、例えば、時間期間ＴＰ１に含まれる１０分間の時間期間が選択される。同様に、ＳＴ０２の動作では、例えば、時間期間ＴＰ４に含まれる１０分間の時間期間が選択されてもよい。
【0082】
次に、ＳＴ０２の動作のより詳細を説明した後に、ＳＴ０３の動作のより詳細を説明する。
図６は、生体信号処理装置１により実行されるＳＴ０２の時間期間を特定する動作のより詳細な一例のフローチャートを示す。例えば、ＳＴ０１の動作に応じて、図６のフローチャートに示される動作が開始される。
図６は、生体信号処理装置１により実行されるＳＴ０２の時間期間を特定する動作のより詳細な一例のフローチャートを示す。例えば、ＳＴ０１の動作に応じて、図６のフローチャートに示される動作が開始される。
【0083】
制御部１１は、第１時間期間特定部１１２１による制御の下、加速度Ａの値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する（ＳＴ２１）。第１の時間期間の特定の説明は、ＳＴ０２の動作について説明した通りである。ＳＴ０２の説明において、特定される時間期間を第１時間期間に置き換え、閾値を第１の閾値に置き換え、待機時間を第１の時間に置き換えた説明が成り立つ。図５の例では、例えば、第１の時間期間として、時間期間ＴＰ１に含まれる１０分間の時間期間が選択される。
【0084】
続いて、制御部１１は、第２時間期間特定部１１２２による制御の下、加速度Ａの値のばらつきの度合いを示す値が第２の閾値以下の第２の時間期間を特定する（ＳＴ２２）。第２の時間期間の始点は、例えば第１の時間期間の終点以降である。第２の時間期間の特定の説明は、ＳＴ０２の動作について説明した通りである。ＳＴ０２の説明において、特定される時間期間を第２時間期間に置き換え、閾値を第２の閾値に置き換え、待機時間を第２の時間に置き換えた説明が成り立つ。第２の閾値は、第１の閾値と同一であっても異なっていてもよい。第２の時間は、第１の時間と同一であっても異なっていてもよい。図５の例では、例えば、第２の時間期間として、時間期間ＴＰ４に含まれる１０分間の時間期間が選択される。
【0085】
【0086】
制御部１１は、生体信号処理部１１３による制御の下、第１の時間期間における生体信号を処理し、第２の時間期間における生体信号を処理する（ＳＴ３１）。当該生体信号の処理では、第１の時間期間と第２の時間期間の各々について、例えば、指標（ｉ１）～（ｉ１２）のうちの少なくとも１つが算出される。当該生体信号の処理では、指標（ｉ１）～（ｉ６）のうちの少なくとも１つが算出されてもよい。計測装置２において指標（ｉ１）～（ｉ１２）のうちの少なくとも１つが計測または算出されている場合、計測装置２は、生体信号として、計測または算出されている指標（ｉ１）～（ｉ１２）のうちの少なくとも１つを生体信号処理装置１に送信していてもよい。この場合、ＳＴ３１の動作は適宜省略可能である。本明細書では、便宜的にＳＴ３１の動作が行われたものとして、以下の説明を行う。
【0087】
制御部１１は、判定部１１３１による制御の下、第１の時間期間における生体信号の処理の結果と、第２の時間期間における生体信号の処理の結果とに基づいて、或る条件式が満足されるか否かを判定する（ＳＴ３２）。ＳＴ３２の動作では、１または２以上の条件式が満足されるか否かが判定されてもよい。より詳細については以降で説明する。以降で説明するように、条件式が満足されるか否かには、例えば、四肢動物が健康異常であること、四肢動物に分娩の予兆があること、または四肢動物に投薬の影響があることが反映され得る。
【0088】
上記条件式が満足されると判定された場合、制御部１１は、出力部１１４による制御の下、当該判定結果に基づいて通知信号を生成し、当該通知信号を出力装置３に出力する（ＳＴ３３）。これにより、出力装置３において、例えば、四肢動物が健康異常であること、四肢動物に分娩の予兆があること、または四肢動物に投薬の影響があること等に係る、音声出力、振動、および／または画面表示等の出力が可能となる。
【0089】
上記条件式が満足されると判定されなかった場合、例えば、ＳＴ２２の動作が再度実行され、続いてＳＴ３１の動作、ＳＴ３２の動作が再度実行される。
【0090】
以下、ＳＴ３１の動作で算出される各指標について、ＳＴ３２の動作で用いられ得る条件式の詳細を説明する。
（ｉ１）拍動間隔の平均値
拍動間隔は、心臓の拍動に伴う電気信号を時系列で記録した心電信号において隣り合うＲ波とＲ波の間隔（以下、ＲＲＩまたはＲＲ間隔とも称され得る。）であること、または、心拍信号としての脈波と脈波の間隔であることが好ましく、ＲＲ間隔であることがより好ましい。拍動間隔が、脈波と脈波の間隔である場合、心臓の脈動に伴って変化する血流量や血管の容積変化を時系列で記録した脈波のパルスのピークであるＰ波と次のＰ波までの間隔であるＰＰ間隔であることがより好ましい。
（ｉ１）拍動間隔の平均値
拍動間隔は、心臓の拍動に伴う電気信号を時系列で記録した心電信号において隣り合うＲ波とＲ波の間隔（以下、ＲＲＩまたはＲＲ間隔とも称され得る。）であること、または、心拍信号としての脈波と脈波の間隔であることが好ましく、ＲＲ間隔であることがより好ましい。拍動間隔が、脈波と脈波の間隔である場合、心臓の脈動に伴って変化する血流量や血管の容積変化を時系列で記録した脈波のパルスのピークであるＰ波と次のＰ波までの間隔であるＰＰ間隔であることがより好ましい。
【0091】
拍動間隔は、心拍を計測する心電計、心拍センサ等により計測することができる。心電計および心拍センサは、四肢動物の胸部および／または腹部に配置されることが好ましく、少なくとも胸部に配置されることが好ましい。拍動間隔は、脈波センサを使用して計測してもよい。脈波センサを使用して脈拍を計測する場合、四肢動物の耳等に、波長７００～１２００ｎｍの近赤外線を照射してその反射量を接触または非接触で計測すればよい。
【0092】
第１の時間期間における拍動間隔の平均値＞第２の時間期間における拍動間隔の平均値・・・（ｃ１）（以下、条件式（ｃ１）とも称され得る。）
が満足される場合、妊娠した四肢動物が、例えば分娩に伴う活動量を増やしている、陣痛を伴っていること等が推測され、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１）が満足されているか否かが判定される。
が満足される場合、妊娠した四肢動物が、例えば分娩に伴う活動量を増やしている、陣痛を伴っていること等が推測され、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１）が満足されているか否かが判定される。
【0093】
四肢動物が牛であり、拍動間隔としてＲＲ間隔が用いられる場合、条件式（ｃ１）の関係に加えて次の関係が満たされることにより、条件式（ｃ１）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における拍動間隔の平均値が１２００ｍｓ超、第２の時間期間における拍動間隔の平均値が１２００ｍｓ以下を満たすことが好ましく、第１の時間期間における拍動間隔の平均値が１０００ｍｓ超、第２の時間期間における拍動間隔の平均値が１０００ｍｓ以下を満たすことがより好ましく、第１の時間期間における拍動間隔の平均値が８００ｍｓ超、第２の時間期間における拍動間隔の平均値が８００ｍｓ以下を満たすことがさらに好ましい。
第１の時間期間における拍動間隔の平均値が１２００ｍｓ超、第２の時間期間における拍動間隔の平均値が１２００ｍｓ以下を満たすことが好ましく、第１の時間期間における拍動間隔の平均値が１０００ｍｓ超、第２の時間期間における拍動間隔の平均値が１０００ｍｓ以下を満たすことがより好ましく、第１の時間期間における拍動間隔の平均値が８００ｍｓ超、第２の時間期間における拍動間隔の平均値が８００ｍｓ以下を満たすことがさらに好ましい。
【0094】
（ｉ２）単位時間当たりの拍数
単位時間当たりの拍数は、単位時間当たりの心拍数または脈拍数であることが好ましく、単位時間当たりの心拍数であることがより好ましい。単位時間当たりの心拍数は、例えば上記ＲＲＩから算出することもできる。単位時間当たりの心拍数は、例えば、特定される時間期間内のＲＲＩ毎に６０／（ＲＲＩ×０．００１）を算出して平均をとったものにより、または、６０／｛（特定される時間期間におけるＲＲＩの平均値）×０．００１｝により算出される。あるいは、単位時間当たりの心拍数は、例えば、特定の期間毎に、上述したいずれかの式で単位時間当たりの心拍数を算出し、さらに、特定される時間期間内について、このように算出された単位時間当たりの心拍数の平均をとったものにより、算出される。単位時間当たりの心拍数の算出に用いられる特定の期間の長さは必ずしも一定でなくてもよい。
単位時間当たりの拍数は、単位時間当たりの心拍数または脈拍数であることが好ましく、単位時間当たりの心拍数であることがより好ましい。単位時間当たりの心拍数は、例えば上記ＲＲＩから算出することもできる。単位時間当たりの心拍数は、例えば、特定される時間期間内のＲＲＩ毎に６０／（ＲＲＩ×０．００１）を算出して平均をとったものにより、または、６０／｛（特定される時間期間におけるＲＲＩの平均値）×０．００１｝により算出される。あるいは、単位時間当たりの心拍数は、例えば、特定の期間毎に、上述したいずれかの式で単位時間当たりの心拍数を算出し、さらに、特定される時間期間内について、このように算出された単位時間当たりの心拍数の平均をとったものにより、算出される。単位時間当たりの心拍数の算出に用いられる特定の期間の長さは必ずしも一定でなくてもよい。
【0095】
単位時間当たりの心拍数は、心電信号の他の何かしらのピークに基づいて算出されるものであってもよい。例えば、四肢動物が牛である場合、（ｉ１）で説明したＲ波の特定が難しいことがあるため、指標（ｉ２）として、心電信号の他のピークに基づいて算出された単位時間当たりの心拍数が用いられることが好ましい。当該単位時間の心拍数は、例えば、ＲＲＩを用いる場合について説明したのと同様に算出される。
【0096】
第１の時間期間での単位時間当たりの拍数＜第２の時間期間での単位時間当たりの拍数・・・（ｃ２）（以下、条件式（ｃ２）とも称され得る。）
が満足される場合、妊娠した四肢動物が例えば陣痛等を伴っていることが推測され、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ２）が満足されているか否かが判定される。
が満足される場合、妊娠した四肢動物が例えば陣痛等を伴っていることが推測され、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ２）が満足されているか否かが判定される。
【0097】
四肢動物が牛であり、単位時間当たりの心拍数が用いられる場合、条件式（ｃ２）の関係に加えて次の関係が満たされることにより、条件式（ｃ２）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
上記単位時間が１分間である場合、第１の時間期間での単位時間当たりの拍数が７０以下、第２の時間期間での単位時間当たりの拍数が７０超を満たすことが好ましく、第１の時間期間での単位時間当たりの拍数が８０以下、第２の時間期間での単位時間当たりの拍数が８０超を満たすことがより好ましく、第１の時間期間での単位時間当たりの拍数が９０以下、第２の時間期間での単位時間当たりの拍数が９０超を満たすことがさらに好ましい。
上記単位時間が１分間である場合、第１の時間期間での単位時間当たりの拍数が７０以下、第２の時間期間での単位時間当たりの拍数が７０超を満たすことが好ましく、第１の時間期間での単位時間当たりの拍数が８０以下、第２の時間期間での単位時間当たりの拍数が８０超を満たすことがより好ましく、第１の時間期間での単位時間当たりの拍数が９０以下、第２の時間期間での単位時間当たりの拍数が９０超を満たすことがさらに好ましい。
【0098】
図９は、ＳＴ３２の或る条件式が満足されるか否かを判定する動作を説明するための例示的なグラフを示す。図９の例は、図５の例に対応している。図９の（Ａ）のグラフは、図５に示したグラフを拡大したものに相当する。
図９の例では、例えば、時間期間ＴＰ１から第１の時間期間が選択され、時間期間ＴＰ４から第２の時間期間が選択されたものとする。
図９の（Ｂ）のグラフでは、縦軸は、５分間についての単位時間としての１分間当たりの心拍数の大きさを５分毎のタイミングで示し、横軸は時刻を示す。
図９の例では、例えば、時間期間ＴＰ１から第１の時間期間が選択され、時間期間ＴＰ４から第２の時間期間が選択されたものとする。
図９の（Ｂ）のグラフでは、縦軸は、５分間についての単位時間としての１分間当たりの心拍数の大きさを５分毎のタイミングで示し、横軸は時刻を示す。
【0099】
例えば、第１の時間期間として、時間期間ＴＰ１のうち、時間期間ＴＰ１の始点から第１の時間としての５分経過した時点からの時間期間が特定され、条件式（ｃ２）で用いられる第１の時間期間での単位時間当たりの拍数として、第１の時間期間に含まれる４つのプロットの心拍数の大きさの平均が算出される。さらに、例えば、第２の時間期間として、時間期間ＴＰ４のうち、時間期間ＴＰ４の始点から第２の時間としての５分経過した時点からの時間期間が特定され、条件式（ｃ２）で用いられる第２の時間期間での単位時間当たりの拍数として、時間期間ＴＰ４に含まれる６つのプロットの心拍数の大きさの平均が算出される。図９の例では、このように算出された第１の時間期間および第２の時間期間それぞれでの単位時間当たりの拍数では条件式（ｃ２）の大小関係は満足されない。その結果、例えば、ＳＴ３３の動作に進むことなく、ＳＴ２２の動作が再度実行される。
【0100】
（ｉ３）拍動間隔のばらつき
拍動間隔のばらつきは、ＲＲＩ等の時間間隔に係る情報を時間解析することにより求められる。
第１の時間期間における拍動間隔のばらつき＞第２の時間期間における拍動間隔のばらつき・・・（ｃ３）（以下、条件式（ｃ３）とも称され得る。）
が満足される場合、副交感神経活動が低下し、交感神経活動が優勢となっており、妊娠した四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ３）が満足されているか否かが判定される。
拍動間隔のばらつきは、ＲＲＩ等の時間間隔に係る情報を時間解析することにより求められる。
第１の時間期間における拍動間隔のばらつき＞第２の時間期間における拍動間隔のばらつき・・・（ｃ３）（以下、条件式（ｃ３）とも称され得る。）
が満足される場合、副交感神経活動が低下し、交感神経活動が優勢となっており、妊娠した四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ３）が満足されているか否かが判定される。
【0101】
拍動間隔のばらつきとして、ＳＤＮＮ、ＲＭＳＳＤ、ＣＶＲＲ、ＮＮ５０、およびｐＮＮ５０のいずれかが用いられることが好ましい。
ＳＤＮＮはＲＲＩの標準偏差であり、ＲＭＳＳＤは、連続して隣接するＲＲＩの差の２乗の平均値の平方根であり、ＣＶＲＲは、ＳＤＮＮの値をＲＲＩの平均値で割って１００を掛けた値であり、ＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える総数を示す値であり、ｐＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える心拍の割合を示す値である。ＳＤＮＮおよびＣＶＲＲは、自律神経活動の指標として見なされ得る。ＲＭＳＳＤ、ＮＮ５０、およびｐＮＮ５０は、副交感神経活動の指標として見なされ得る。
ＳＤＮＮはＲＲＩの標準偏差であり、ＲＭＳＳＤは、連続して隣接するＲＲＩの差の２乗の平均値の平方根であり、ＣＶＲＲは、ＳＤＮＮの値をＲＲＩの平均値で割って１００を掛けた値であり、ＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える総数を示す値であり、ｐＮＮ５０は、連続した隣接するＲＲＩの差が５０ｍｓを超える心拍の割合を示す値である。ＳＤＮＮおよびＣＶＲＲは、自律神経活動の指標として見なされ得る。ＲＭＳＳＤ、ＮＮ５０、およびｐＮＮ５０は、副交感神経活動の指標として見なされ得る。
【0102】
四肢動物が牛であり、拍動間隔のばらつきとしてＳＤＮＮが用いられる場合、条件式（ｃ３）の関係に加えて次の関係が満たされることにより、条件式（ｃ３）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における拍動間隔のばらつきが４０超、第２の時間期間における拍動間隔のばらつきが４０以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが１４０超、第２の時間期間における拍動間隔のばらつきが１４０以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが２４０超、第２の時間期間における拍動間隔のばらつきが２４０以下であることがさらに好ましい。これらの数値範囲は、ＲＲＩの単位としてｍｓを用いた場合のものである。次に説明するＲＭＳＳＤおよびＣＶＲＲについての数値範囲についても同様である。
第１の時間期間における拍動間隔のばらつきが４０超、第２の時間期間における拍動間隔のばらつきが４０以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが１４０超、第２の時間期間における拍動間隔のばらつきが１４０以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが２４０超、第２の時間期間における拍動間隔のばらつきが２４０以下であることがさらに好ましい。これらの数値範囲は、ＲＲＩの単位としてｍｓを用いた場合のものである。次に説明するＲＭＳＳＤおよびＣＶＲＲについての数値範囲についても同様である。
【0103】
四肢動物が牛であり、拍動間隔のばらつきとしてＲＭＳＳＤが用いられる場合、条件式（ｃ３）の関係に加えて次の関係が満たされることにより、条件式（ｃ３）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における拍動間隔のばらつきが１２超、第２の時間期間における拍動間隔のばらつきが１２以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが１２０超、第２の時間期間における拍動間隔のばらつきが１２０以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが１４０超、第２の時間期間における拍動間隔のばらつきが１４０以下であることがさらに好ましい。
第１の時間期間における拍動間隔のばらつきが１２超、第２の時間期間における拍動間隔のばらつきが１２以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが１２０超、第２の時間期間における拍動間隔のばらつきが１２０以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが１４０超、第２の時間期間における拍動間隔のばらつきが１４０以下であることがさらに好ましい。
【0104】
四肢動物が牛であり、拍動間隔のばらつきとしてＣＶＲＲが用いられる場合、条件式（ｃ３）の関係に加えて次の関係が満たされることにより、条件式（ｃ３）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における拍動間隔のばらつきが２５超、第２の時間期間における拍動間隔のばらつきが２５以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが３５超、第２の時間期間における拍動間隔のばらつきが３５以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが４５超、第２の時間期間における拍動間隔のばらつきが４５以下であることがさらに好ましい。
第１の時間期間における拍動間隔のばらつきが２５超、第２の時間期間における拍動間隔のばらつきが２５以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが３５超、第２の時間期間における拍動間隔のばらつきが３５以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが４５超、第２の時間期間における拍動間隔のばらつきが４５以下であることがさらに好ましい。
【0105】
四肢動物が牛であり、拍動間隔のばらつきとしてＮＮ５０が用いられる場合、条件式（ｃ３）の関係に加えて次の関係が満たされることにより、条件式（ｃ３）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における拍動間隔のばらつきが３０超、第２の時間期間における拍動間隔のばらつきが３０以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが４５超、第２の時間期間における拍動間隔のばらつきが４５以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが６０超、第２の時間期間における拍動間隔のばらつきが６０以下であることがさらに好ましい。
第１の時間期間における拍動間隔のばらつきが３０超、第２の時間期間における拍動間隔のばらつきが３０以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが４５超、第２の時間期間における拍動間隔のばらつきが４５以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが６０超、第２の時間期間における拍動間隔のばらつきが６０以下であることがさらに好ましい。
【0106】
四肢動物が牛であり、拍動間隔のばらつきとしてｐＮＮ５０が用いられる場合、条件式（ｃ３）の関係に加えて次の関係が満たされることにより、条件式（ｃ３）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における拍動間隔のばらつきが８０超、第２の時間期間における拍動間隔のばらつきが８０以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが９０超、第２の時間期間における拍動間隔のばらつきが９０以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが１００超、第２の時間期間における拍動間隔のばらつきが１００以下であることがさらに好ましい。
第１の時間期間における拍動間隔のばらつきが８０超、第２の時間期間における拍動間隔のばらつきが８０以下であることが好ましく、第１の時間期間における拍動間隔のばらつきが９０超、第２の時間期間における拍動間隔のばらつきが９０以下であることがより好ましく、第１の時間期間における拍動間隔のばらつきが１００超、第２の時間期間における拍動間隔のばらつきが１００以下であることがさらに好ましい。
【0107】
（ｉ４）単位時間当たりの拍数のばらつき
例えば、指標（ｉ２）に関して上記で説明したように、特定される時間期間内の複数のタイミングの各々について、当該タイミングを含む或る特定の期間について単位時間当たりの拍数が取得されている場合、これら単位時間当たりの拍数のばらつきの度合いを示す値（以下、単に、単位時間当たりの拍数のばらつきとも称され得る。）が算出可能である。当該ばらつきとしては、上記で加速度Ａについて説明したのと同様に、例えば、分散、標準偏差、および変動係数のいずれかが用いられる。
第１の時間期間での単位時間当たりの拍数のばらつき＞第２の時間期間での単位時間当たりの拍数のばらつき・・・（ｃ４）（以下、条件式（ｃ４）とも称され得る。）
が満足される場合、妊娠した四肢動物が例えば陣痛等を伴っていることが推測され、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ４）が満足されているか否かが判定されてもよい。
例えば、指標（ｉ２）に関して上記で説明したように、特定される時間期間内の複数のタイミングの各々について、当該タイミングを含む或る特定の期間について単位時間当たりの拍数が取得されている場合、これら単位時間当たりの拍数のばらつきの度合いを示す値（以下、単に、単位時間当たりの拍数のばらつきとも称され得る。）が算出可能である。当該ばらつきとしては、上記で加速度Ａについて説明したのと同様に、例えば、分散、標準偏差、および変動係数のいずれかが用いられる。
第１の時間期間での単位時間当たりの拍数のばらつき＞第２の時間期間での単位時間当たりの拍数のばらつき・・・（ｃ４）（以下、条件式（ｃ４）とも称され得る。）
が満足される場合、妊娠した四肢動物が例えば陣痛等を伴っていることが推測され、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ４）が満足されているか否かが判定されてもよい。
【0108】
（ｉ５）ＨＦ、（ｉ６）ＬＦ／ＨＦ
第１の時間期間のＨＦ＞第２の時間期間のＨＦ・・・（ｃ５）（以下、条件式（ｃ５）とも称され得る。）、
第１の時間期間のＬＦ／ＨＦ＜第２の時間期間のＬＦ／ＨＦ・・・（ｃ６）（以下、条件式（ｃ６）とも称され得る。）
が満足される場合、副交感神経活動が低下していること、および／または、交感神経活動が優勢になっていることにより、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ５）と条件式（ｃ６）との少なくとも一方が満足されているか否かが判定されるが、条件式（ｃ５）と条件式（ｃ６）の一方だけではなく両方が満足されているか否かが判定されることが好ましい。
第１の時間期間のＨＦ＞第２の時間期間のＨＦ・・・（ｃ５）（以下、条件式（ｃ５）とも称され得る。）、
第１の時間期間のＬＦ／ＨＦ＜第２の時間期間のＬＦ／ＨＦ・・・（ｃ６）（以下、条件式（ｃ６）とも称され得る。）
が満足される場合、副交感神経活動が低下していること、および／または、交感神経活動が優勢になっていることにより、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ５）と条件式（ｃ６）との少なくとも一方が満足されているか否かが判定されるが、条件式（ｃ５）と条件式（ｃ６）の一方だけではなく両方が満足されているか否かが判定されることが好ましい。
【0109】
ＨＦは、呼吸からもたらされる副交感神経活動の指標とみなされてもよく、ＬＦは、交感神経活動および副交感神経活動の指標とみなされてもよく、ＬＦ／ＨＦは、交感神経活動の指標とみなされてもよい。ＨＦ、および／または、ＬＦ／ＨＦを指標とすると、四肢動物がリラックスしているか、外因（例えば四肢動物が活動している環境、妊娠した四肢動物以外の生物（例えば人間））または内因（例えば分娩に関連する活動）の影響を受けて四肢動物がリラックスしていないかを推測することができる。
【0110】
ＨＦおよびＬＦの詳細を次に説明する。
ＨＦおよびＬＦは各々、例えば、ＲＲＩを周波数解析してパワースペクトルを計算し、所定の周波数領域のパワーを積分して求めることができる。ＬＦは、拍動間隔を周波数スペクトル変換することにより取得されるパワースペクトルを周波数Ｌｆ１からＬｆ２まで定積分した値であり、ＨＦは、当該パワースペクトルを周波数Ｈｆ１からＨｆ２まで定積分した値であり、Ｈｆ１＞Ｌｆ１、Ｈｆ２＞Ｌｆ２である。
ＨＦおよびＬＦは各々、例えば、ＲＲＩを周波数解析してパワースペクトルを計算し、所定の周波数領域のパワーを積分して求めることができる。ＬＦは、拍動間隔を周波数スペクトル変換することにより取得されるパワースペクトルを周波数Ｌｆ１からＬｆ２まで定積分した値であり、ＨＦは、当該パワースペクトルを周波数Ｈｆ１からＨｆ２まで定積分した値であり、Ｈｆ１＞Ｌｆ１、Ｈｆ２＞Ｌｆ２である。
【0111】
図８は、ＬＦおよびＨＦの算出方法の一例を説明するためのグラフを示す。当該グラフにおいて、縦軸はパワースペクトル密度（単位：ｍｓ２／Ｈｚ）の大きさを示し、横軸は周波数（単位：Ｈｚ）の大きさを示す。ＬＦは、パワースペクトルＦ2を例えば０．０４Ｈｚ（Ｌｆ１）から０．１５Ｈｚ（Ｌｆ２）まで定積分した値であり、図８のグラフにおいて斜線によりハッチングがされている部分の面積である。ＨＦは、パワースペクトルＦ2を例えば０．１５Ｈｚ（Ｈｆ１）から０．４０Ｈｚ（Ｈｆ２）まで定積分した値であり、図８のグラフにおいて縦線によりハッチングがされている部分の面積である。図８では、Ｌｆ２とＨｆ１がいずれも０．１５Ｈｚと等しくなるように積分範囲を設定する場合の例を示したが、積分範囲は上記範囲に限定されない。
【0112】
Ｌｆ１＜Ｈｆ１およびＬｆ２＜Ｈｆ２の関係を満たしていることが好ましく、さらにＨｆ１≧Ｌｆ２の関係を満たしていることがより好ましい。またＬＦの積分範囲が０．０１３３～０．３０Ｈｚであり、且つＨＦの積分範囲が０．３０～０．８０Ｈｚであってもよい。
ＬＦの積分範囲は、少なくとも０．１０Ｈｚを含み、Ｌｆ１＜０．１０＜Ｌｆ２であることが好ましい。Ｌｆ１は０．０１～０．０９Ｈｚであることが好ましく、０．０１～０．０５Ｈｚであることがより好ましい。Ｌｆ２は０．１１～０．３４Ｈｚであることがより好ましい。ＨＦの積分範囲は、少なくとも０．３５Ｈｚを含み、Ｈｆ１＜０．３５＜Ｈｆ２であることが好ましい。Ｈｆ１は０．１０～０．３４Ｈｚであることがより好ましい。Ｈｆ２は０．３６～０．９０Ｈｚであることが好ましく、０．３８～０．９０Ｈｚであることがより好ましい。
以上のＬＦ、ＨＦの算出に当たっては、国際公開第２０１６／０３１６５０号等を参照することができる。
ＬＦの積分範囲は、少なくとも０．１０Ｈｚを含み、Ｌｆ１＜０．１０＜Ｌｆ２であることが好ましい。Ｌｆ１は０．０１～０．０９Ｈｚであることが好ましく、０．０１～０．０５Ｈｚであることがより好ましい。Ｌｆ２は０．１１～０．３４Ｈｚであることがより好ましい。ＨＦの積分範囲は、少なくとも０．３５Ｈｚを含み、Ｈｆ１＜０．３５＜Ｈｆ２であることが好ましい。Ｈｆ１は０．１０～０．３４Ｈｚであることがより好ましい。Ｈｆ２は０．３６～０．９０Ｈｚであることが好ましく、０．３８～０．９０Ｈｚであることがより好ましい。
以上のＬＦ、ＨＦの算出に当たっては、国際公開第２０１６／０３１６５０号等を参照することができる。
【0113】
四肢動物が牛である場合、条件式（ｃ５）の関係に加えて次の関係が満たされることにより、条件式（ｃ５）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間のＨＦが１０００超、第２の時間期間のＨＦが１０００以下であることが好ましく、第１の時間期間のＨＦが３０００超、第２の時間期間のＨＦが３０００以下であることがより好ましく、第１の時間期間のＨＦが５０００超、第２の時間期間のＨＦが５０００以下であることがさらに好ましい。これらの数値範囲でのＨＦの単位はｍｓ２である。
第１の時間期間のＨＦが１０００超、第２の時間期間のＨＦが１０００以下であることが好ましく、第１の時間期間のＨＦが３０００超、第２の時間期間のＨＦが３０００以下であることがより好ましく、第１の時間期間のＨＦが５０００超、第２の時間期間のＨＦが５０００以下であることがさらに好ましい。これらの数値範囲でのＨＦの単位はｍｓ２である。
【0114】
四肢動物が牛である場合、条件式（ｃ６）の関係に加えて次の関係が満たされることにより、条件式（ｃ６）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間のＬＦ／ＨＦが７５以下、第２の時間期間のＬＦ／ＨＦが７５超であることが好ましく、第１の時間期間のＬＦ／ＨＦが５０以下、第２の時間期間のＬＦ／ＨＦが５０超であることがより好ましく、第１の時間期間のＬＦ／ＨＦが２５以下、第２の時間期間のＬＦ／ＨＦが２５超であることがさらに好ましい。
第１の時間期間のＬＦ／ＨＦが７５以下、第２の時間期間のＬＦ／ＨＦが７５超であることが好ましく、第１の時間期間のＬＦ／ＨＦが５０以下、第２の時間期間のＬＦ／ＨＦが５０超であることがより好ましく、第１の時間期間のＬＦ／ＨＦが２５以下、第２の時間期間のＬＦ／ＨＦが２５超であることがさらに好ましい。
【0115】
（ｉ７）加速度、（ｉ８）角速度
加速度は、上述した加速度Ａであることが好ましい。加速度は、上述した加速度Ａの平均であることがより好ましい。
角速度は、加速度データについて説明したのと同様に、例えば、加速度計等を含む計測装置２により計測できるか、計測装置２によって計測されるデータに基づいて算出可能な、四肢動物に係る角速度である。特定される時間期間内で角速度の値が複数計測または算出される場合、上記角速度は、計測または算出される角速度の値の平均であることがより好ましい。
第１の時間期間における加速度＜第２の時間期間における加速度・・・（ｃ７）（以下、条件式（ｃ７）とも称され得る。）、
第１の時間期間における角速度＜第２の時間期間における角速度・・・（ｃ８）（以下、条件式（ｃ８）とも称され得る。）
が満足される場合、四肢動物の活動量が増えており、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ７）と条件式（ｃ８）との少なくとも一方が満足されているか否かが判定されるが、条件式（ｃ７）と条件式（ｃ８）の一方だけではなく両方が満足されているか否かが判定されることが好ましい。
加速度は、上述した加速度Ａであることが好ましい。加速度は、上述した加速度Ａの平均であることがより好ましい。
角速度は、加速度データについて説明したのと同様に、例えば、加速度計等を含む計測装置２により計測できるか、計測装置２によって計測されるデータに基づいて算出可能な、四肢動物に係る角速度である。特定される時間期間内で角速度の値が複数計測または算出される場合、上記角速度は、計測または算出される角速度の値の平均であることがより好ましい。
第１の時間期間における加速度＜第２の時間期間における加速度・・・（ｃ７）（以下、条件式（ｃ７）とも称され得る。）、
第１の時間期間における角速度＜第２の時間期間における角速度・・・（ｃ８）（以下、条件式（ｃ８）とも称され得る。）
が満足される場合、四肢動物の活動量が増えており、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ７）と条件式（ｃ８）との少なくとも一方が満足されているか否かが判定されるが、条件式（ｃ７）と条件式（ｃ８）の一方だけではなく両方が満足されているか否かが判定されることが好ましい。
【0116】
（ｉ９）単位時間当たりの呼吸数
単位時間当たりの呼吸数は、１分当たりの呼吸回数として表されてもよい。呼吸数は、心電図およびパルスオキシメータ等を用いて計測すればよく、呼吸に伴う胸腹部の変位を計測するセンサを用いて計測されてもよい。
第１の時間期間での単位時間当たりの呼吸数＜第２の時間期間での単位時間当たりの呼吸数・・・（ｃ９）（以下、条件式（ｃ９）とも称され得る。）
が満足される場合、四肢動物の呼吸が荒くなっており、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ９）が満足されているか否かが判定される。
単位時間当たりの呼吸数は、１分当たりの呼吸回数として表されてもよい。呼吸数は、心電図およびパルスオキシメータ等を用いて計測すればよく、呼吸に伴う胸腹部の変位を計測するセンサを用いて計測されてもよい。
第１の時間期間での単位時間当たりの呼吸数＜第２の時間期間での単位時間当たりの呼吸数・・・（ｃ９）（以下、条件式（ｃ９）とも称され得る。）
が満足される場合、四肢動物の呼吸が荒くなっており、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ９）が満足されているか否かが判定される。
【0117】
四肢動物が牛である場合、条件式（ｃ９）の関係に加えて次の関係が満たされることにより、条件式（ｃ９）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
上記単位時間が１分間である場合、第１の時間期間での単位時間当たりの呼吸数が３５以下、第２の時間期間での単位時間当たりの呼吸数が３５超であることが好ましく、第１の時間期間での単位時間当たりの呼吸数が３８以下、第２の時間期間での単位時間当たりの呼吸数が３８超であることがより好ましく、第１の時間期間での単位時間当たりの呼吸数が４１以下、第２の時間期間での単位時間当たりの呼吸数が４１超であることがさらに好ましい。
上記単位時間が１分間である場合、第１の時間期間での単位時間当たりの呼吸数が３５以下、第２の時間期間での単位時間当たりの呼吸数が３５超であることが好ましく、第１の時間期間での単位時間当たりの呼吸数が３８以下、第２の時間期間での単位時間当たりの呼吸数が３８超であることがより好ましく、第１の時間期間での単位時間当たりの呼吸数が４１以下、第２の時間期間での単位時間当たりの呼吸数が４１超であることがさらに好ましい。
【0118】
（ｉ１０）外陰部間の距離
外陰部間の距離は、電磁誘導、赤外線、超音波等を用いて測定できるものであればよく、外陰部の一方および他方に各種の受信器および発振器を配置することにより、外陰部が開いたことを検知する方法を採用すればよい。受信器および発振器はそれぞれ、２つ以上を使用してもよい。なかでも、電磁誘導を用いた外陰部間の距離測定は、分娩直前に取り付けることが望ましく、継続的に取り付けて測定する場合には四肢動物の活動に伴い誤作動する虞があり、外陰部が開口すると作動しなくなることから、赤外線、超音波を使用して、外陰部間の距離を測定することが好ましい。特定される時間期間内で外陰部間の距離の値が複数測定される場合、上記外陰部間の距離は、測定される外陰部間の距離の値の平均であることがより好ましい。
第１の時間期間における外陰部間の距離＜第２の時間期間における外陰部間の距離・・・（ｃ１０）（以下、条件式（ｃ１０）とも称され得る。）
が満足される場合、産道から胎児が出てきて外陰部が開口したことを意味し、四肢動物の分娩が開始したことが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１０）が満足されているか否かが判定される。
外陰部間の距離は、電磁誘導、赤外線、超音波等を用いて測定できるものであればよく、外陰部の一方および他方に各種の受信器および発振器を配置することにより、外陰部が開いたことを検知する方法を採用すればよい。受信器および発振器はそれぞれ、２つ以上を使用してもよい。なかでも、電磁誘導を用いた外陰部間の距離測定は、分娩直前に取り付けることが望ましく、継続的に取り付けて測定する場合には四肢動物の活動に伴い誤作動する虞があり、外陰部が開口すると作動しなくなることから、赤外線、超音波を使用して、外陰部間の距離を測定することが好ましい。特定される時間期間内で外陰部間の距離の値が複数測定される場合、上記外陰部間の距離は、測定される外陰部間の距離の値の平均であることがより好ましい。
第１の時間期間における外陰部間の距離＜第２の時間期間における外陰部間の距離・・・（ｃ１０）（以下、条件式（ｃ１０）とも称され得る。）
が満足される場合、産道から胎児が出てきて外陰部が開口したことを意味し、四肢動物の分娩が開始したことが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１０）が満足されているか否かが判定される。
【0119】
四肢動物が牛である場合、条件式（ｃ１０）の関係に加えて次の関係が満たされることにより、条件式（ｃ１０）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第１の時間期間における外陰部間の距離が１０ｃｍ以下であり、第２の時間期間における外陰部間の距離が１０ｃｍ超であることが好ましく、第１の時間期間における外陰部間の距離が１１ｃｍ以下であり、第２の時間期間における外陰部間の距離が１１ｃｍ超であることよりが好ましく、第１の時間期間における外陰部間の距離が１２ｃｍ以下であり、第２の時間期間における外陰部間の距離が１２ｃｍ超であることがさらに好ましい。
第１の時間期間における外陰部間の距離が１０ｃｍ以下であり、第２の時間期間における外陰部間の距離が１０ｃｍ超であることが好ましく、第１の時間期間における外陰部間の距離が１１ｃｍ以下であり、第２の時間期間における外陰部間の距離が１１ｃｍ超であることよりが好ましく、第１の時間期間における外陰部間の距離が１２ｃｍ以下であり、第２の時間期間における外陰部間の距離が１２ｃｍ超であることがさらに好ましい。
【0120】
（ｉ１１）体表温度
体表温度は、鼓膜温度であることが好ましい。体表温度の測定に用いられる温度計は、金属酸化物、半導体のサーミスタからなる熱型赤外センサであってもよく、例えば、温度に応じて電気抵抗が変化することを利用するものである。特定される時間期間内で体表温度の値が複数測定される場合、上記体表温度は、測定される体表温度の値の平均であることがより好ましい。
第１の時間期間における体表温度＞第２の時間期間における体表温度・・・（ｃ１１）（以下、条件式（ｃ１１）とも称され得る。）
が満足される場合、四肢動物の体表温度が低下して、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１１）が満足されているか否かが判定される。
体表温度は、鼓膜温度であることが好ましい。体表温度の測定に用いられる温度計は、金属酸化物、半導体のサーミスタからなる熱型赤外センサであってもよく、例えば、温度に応じて電気抵抗が変化することを利用するものである。特定される時間期間内で体表温度の値が複数測定される場合、上記体表温度は、測定される体表温度の値の平均であることがより好ましい。
第１の時間期間における体表温度＞第２の時間期間における体表温度・・・（ｃ１１）（以下、条件式（ｃ１１）とも称され得る。）
が満足される場合、四肢動物の体表温度が低下して、四肢動物の分娩が近づいていることが推測される。あるいは、四肢動物の体調が悪化している、四肢動物に投薬の影響が強く出ていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１１）が満足されているか否かが判定される。
【0121】
四肢動物が牛である場合、条件式（ｃ１１）の関係に加えて次の関係が満たされることにより、条件式（ｃ１１）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第２の時間期間における体表温度が、第１の時間期間における体表温度よりも０．３～１．０℃低いことが好ましく、０．５～１．０℃低いことがより好ましい。
第２の時間期間における体表温度が、第１の時間期間における体表温度よりも０．３～１．０℃低いことが好ましく、０．５～１．０℃低いことがより好ましい。
【0122】
（ｉ１２）膣内温度
特定される時間期間内で膣内温度の値が複数測定される場合、上記膣内温度は、測定される膣内温度の値の平均であることがより好ましい。
第１の時間期間における膣内温度＞第２の時間期間における膣内温度・・・（ｃ１２）（以下、条件式（ｃ１２）とも称され得る。）
が満足される場合、四肢動物の膣内温度が低下して、四肢動物の分娩が近づいていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１２）が満足されているか否かが判定される。
特定される時間期間内で膣内温度の値が複数測定される場合、上記膣内温度は、測定される膣内温度の値の平均であることがより好ましい。
第１の時間期間における膣内温度＞第２の時間期間における膣内温度・・・（ｃ１２）（以下、条件式（ｃ１２）とも称され得る。）
が満足される場合、四肢動物の膣内温度が低下して、四肢動物の分娩が近づいていることが推測される。ＳＴ３２の動作では、例えば、条件式（ｃ１２）が満足されているか否かが判定される。
【0123】
四肢動物が牛である場合、条件式（ｃ１２）の関係に加えて次の関係が満たされることにより、条件式（ｃ１２）を用いて四肢動物の状態を上述のように推測する精度が向上され得る。
第２の時間期間における膣内温度が、第１の時間期間における膣内温度よりも０．３～１．０℃低いことが好ましく、０．４～１．０℃低いことがより好ましい。
第２の時間期間における膣内温度が、第１の時間期間における膣内温度よりも０．３～１．０℃低いことが好ましく、０．４～１．０℃低いことがより好ましい。
【0124】
条件式（ｃ１）～（ｃ６）のうちの少なくとも１つが満足されることが好ましく、条件式（ｃ１）～（ｃ６）の２以上が満足されることがより好ましく、条件式（ｃ１）～（ｃ６）の４以上が満足されることがさらに好ましく、条件式（ｃ１）～（ｃ６）の全てが満足されることが最も好ましい。
【0125】
また、条件式（ｃ１）～（ｃ１０）のうちの少なくとも１つが満足されることが好ましく、条件式（ｃ１）～（ｃ１０）の２以上が満足されることがより好ましく、条件式（ｃ１）～（ｃ１０）の４以上が満足されることがより好ましく、条件式（ｃ１）～（ｃ１０）の６以上が満足されることがさらに好ましく、条件式（ｃ１）～（ｃ１０）の８以上が満足されることがさらにより好ましく、条件式（ｃ１）～（ｃ１０）の全てが満足されることが最も好ましい。
【0126】
また、条件式（ｃ１）～（ｃ１２）のうちの少なくとも１つが満足されることが好ましく、条件式（ｃ１）～（ｃ１２）の２以上が満足されることがより好ましく、条件式（ｃ１）～（ｃ１２）の４以上が満足されることがより好ましく、条件式（ｃ１）～（ｃ１２）の６以上が満足されることがさらに好ましく、条件式（ｃ１）～（ｃ１２）の８以上が満足されることがさらにより好ましく、条件式（ｃ１）～（ｃ１２）の１０以上が満足されることが特に好ましく、条件式（ｃ１）～（ｃ１２）の全てが満足されることが最も好ましい。
【0127】
四肢動物が牛または犬である場合、ＳＴ３２の動作では、条件式（ｃ２）と（ｃ９）との少なくとも一方が用いられることが好ましく、条件式（ｃ２）が用いられることがより好ましい。
【0128】
条件式（ｃ１）、（ｃ３）～（ｃ８）のうちの少なくとも１つが満足されることが好ましく、条件式（ｃ１）、（ｃ３）～（ｃ８）の２以上が満足されることがより好ましく、条件式（ｃ１）、（ｃ３）～（ｃ８）の４以上が満足されることがさらに好ましく、条件式（ｃ１）、（ｃ３）～（ｃ８）の全てが満足されることがさらにより好ましい。
【0129】
条件式（ｃ２）、（ｃ３）、および（ｃ１２）のうちの少なくとも１つが満足されることが好ましく、条件式（ｃ２）、（ｃ３）、および（ｃ１２）の２以上が満足されることがより好ましく、条件式（ｃ２）および（ｃ３）が満足されること、条件式（ｃ２）および（ｃ１２）が満足されること、条件式（ｃ３）および（ｃ１２）が満足されることがより好ましく、条件式（ｃ２）、（ｃ３）および（ｃ１２）が満足されることがさらに好ましい。
【0130】
上述したＳＴ３２の動作では、例えば、条件式（ｃ１）～（ｃ１２）のうちの少なくとも１つが満足されるか否かが判定され、当該少なくとも１つの条件式が満足されると判定された場合に、ＳＴ３３の動作において、通知信号が生成され、当該通知信号が出力装置３に出力されてもよい。一方、上述したＳＴ３２の動作において、上記で条件式（ｃ１）～（ｃ１２）のうち満足されていることが好ましいと説明した組み合わせの任意のものに係る条件式がすべて満足されているか否かが判定され、満足されていると判定された場合に、ＳＴ３３の動作において、通知信号が生成され、当該通知信号が出力装置３に出力されてもよい。このようにして、出力装置３において、例えば、妊娠している四肢動物の分娩接近、分娩開始等を知らせる警告等に係る音声出力、振動、および／または画面表示等の出力が可能となる、あるいは、四肢動物が健康異常であること、四肢動物に投薬の影響があることに係る、上述したような出力が可能となる。
【0131】
以上、生体信号処理装置１の動作例を詳細に説明したが、上記で説明した動作フローは一例に過ぎない。
例えば、生体信号処理装置１の制御部１１は、時間期間特定部１１２による制御の下、加速度Ａの値のばらつきの度合いを示す値が或る閾値以下の第３の時間期間を特定してもよい。第３の時間期間の始点は、例えば第２の時間期間の終点以降である。第２の時間期間と第３の時間期間は、連続していても連続していなくてもよい。第３の時間期間は、例えば５分以上、５時間以下、好ましくは１０分以上、４時間以下、より好ましくは２０分以上、３時間以下、さらに好ましくは３０分以上、２時間以下である。第３の時間期間は、難産となる場合に設けられることが好ましく、陣痛開始後～胎児摘出までの期間であってもよい。第３の時間期間の特定に用いる閾値については、ＳＴ０２の動作の説明において行った閾値についての説明が成り立つ。当該閾値は、第１の閾値および第２の閾値のいずれとも異なっていてもよく、第１の閾値および第２の閾値の少なくとも一方と同一であってもよい。
続いて、ＳＴ３１の動作について説明したのと同様に、制御部１１は、生体信号処理部１１３による制御の下、第３の時間期間における生体信号を処理してもよい。当該生体信号の処理では、例えば、前述した指標（ｉ１）～（ｉ１２）のうちの少なくとも１つが算出される。指標（ｉ１）～（ｉ６）のうちの少なくとも１つが算出されてもよい。
第３時間期間についてのかかる指標をモニタリングしながら、異常な数値が発見されれば、四肢動物の母体および胎児に適切な処置を施すことが可能となる。
例えば、生体信号処理装置１の制御部１１は、時間期間特定部１１２による制御の下、加速度Ａの値のばらつきの度合いを示す値が或る閾値以下の第３の時間期間を特定してもよい。第３の時間期間の始点は、例えば第２の時間期間の終点以降である。第２の時間期間と第３の時間期間は、連続していても連続していなくてもよい。第３の時間期間は、例えば５分以上、５時間以下、好ましくは１０分以上、４時間以下、より好ましくは２０分以上、３時間以下、さらに好ましくは３０分以上、２時間以下である。第３の時間期間は、難産となる場合に設けられることが好ましく、陣痛開始後～胎児摘出までの期間であってもよい。第３の時間期間の特定に用いる閾値については、ＳＴ０２の動作の説明において行った閾値についての説明が成り立つ。当該閾値は、第１の閾値および第２の閾値のいずれとも異なっていてもよく、第１の閾値および第２の閾値の少なくとも一方と同一であってもよい。
続いて、ＳＴ３１の動作について説明したのと同様に、制御部１１は、生体信号処理部１１３による制御の下、第３の時間期間における生体信号を処理してもよい。当該生体信号の処理では、例えば、前述した指標（ｉ１）～（ｉ１２）のうちの少なくとも１つが算出される。指標（ｉ１）～（ｉ６）のうちの少なくとも１つが算出されてもよい。
第３時間期間についてのかかる指標をモニタリングしながら、異常な数値が発見されれば、四肢動物の母体および胎児に適切な処置を施すことが可能となる。
【0132】
（効果）
第１実施形態に係る生体信号処理装置１は、四肢動物の生体信号と、当該四肢動物の動きに伴う加速度Ａとを取得し、加速度Ａの値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する。加速度Ａの値は、例えば、互いに交わるＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度ｘ、ｙ、およびｚのベクトル和の大きさに対応する。このように特定される第１の時間期間では、例えば、当該四肢動物が臥位であると推測される。
第１実施形態に係る生体信号処理装置１は、四肢動物の生体信号と、当該四肢動物の動きに伴う加速度Ａとを取得し、加速度Ａの値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間を特定する。加速度Ａの値は、例えば、互いに交わるＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度ｘ、ｙ、およびｚのベクトル和の大きさに対応する。このように特定される第１の時間期間では、例えば、当該四肢動物が臥位であると推測される。
【0133】
例えば、人が立位の時に上下である方向の加速度を計測するようにすると、人が立位および座位の場合と比較して、人が臥位の場合に当該加速度の大きさは小さい。このため、当該加速度の大きさのみに基づいて人が臥位であることを推測することが可能である。一方、四肢動物が立位の場合と臥位の場合での姿勢の変化は人の場合と大きく異なるため、同様の方法で四肢動物が臥位であることを推測することは困難である。このため、四肢動物の姿勢の推測には、加速度として、例えば、互いに交わるＸ軸、Ｙ軸、およびＺ軸それぞれの方向の加速度ｘ、ｙ、およびｚのベクトル和の大きさを用いることが好ましい。
【0134】
ここで、四肢動物が臥位となるときには、身体の左側を下にして臥位になる場合と、身体の右側を下にして臥位になる場合とがある。このような２種類の臥位の間では、上記加速度が異なる大きさで安定することがある。安定する大きさのうちの少なくとも一方が、四肢動物が立位の場合の加速度の大きさより大きいことがある。このため、上記加速度の大きさのみに基づいて四肢動物が臥位であることを推測することは困難である。一方、生体信号処理装置１によれば、加速度Ａの値のばらつきの度合いを示す値が第１の閾値以下の第１の時間期間が特定される。四肢動物が上述した２種類の臥位のいずれであっても、加速度Ａの値のばらつきの度合いを示す値は、四肢動物が立位の場合の加速度Ａの値のばらつきの度合いを示す値より小さいことが考えられる。このため、このように特定される第１の時間期間では、四肢動物が臥位であると推測される精度が向上される。
上記では、横臥位（側臥位）を例に挙げて説明を行ったが、四肢動物の種類によっては、四肢動物は、例えば、背中を下とする仰臥位（背臥位）、および、腹を下とする伏臥位（腹臥位）の姿勢も取り得る。これらの臥位の場合も、上記加速度が他の臥位の場合とは異なる大きさで安定することがあるが、上述したのと同様に、加速度Ａの値のばらつきの度合いを示す値は、立位の場合より小さいことが考えられる。このため、四肢動物が取り得る臥位の種類によらず、上述した効果が奏され得る。
上記では、横臥位（側臥位）を例に挙げて説明を行ったが、四肢動物の種類によっては、四肢動物は、例えば、背中を下とする仰臥位（背臥位）、および、腹を下とする伏臥位（腹臥位）の姿勢も取り得る。これらの臥位の場合も、上記加速度が他の臥位の場合とは異なる大きさで安定することがあるが、上述したのと同様に、加速度Ａの値のばらつきの度合いを示す値は、立位の場合より小さいことが考えられる。このため、四肢動物が取り得る臥位の種類によらず、上述した効果が奏され得る。
【0135】
生体信号処理装置１は、上記生体信号のうちの上記第１の時間期間における生体信号に基づいて出力を行う。これにより、四肢動物の動作および姿勢が当該四肢動物の生体信号に影響を及ぼしてしまうような時間期間の生体信号を排除して、四肢動物が臥位であると推測される時間期間での当該四肢動物の生体信号に基づく出力が行われることになる。
【0136】
生体信号処理装置１は、例えば、上記第１の時間期間より後の、上記ばらつきの度合いを示す値が第２の閾値以下の第２の時間期間を特定する。生体信号処理装置１は、例えば、第１の時間期間における生体信号の処理の結果と、第２の時間期間における生体信号の処理の結果とに基づいて、或る条件式が満足されるか否かを判定する。条件式が満足されるか否かには、例えば、四肢動物が健康異常であること、四肢動物に分娩の予兆があること、四肢動物に投薬の影響があることが反映され得る。生体信号処理装置１は、上記判定結果に基づいて通知信号を生成し、上記出力として、当該通知信号を出力装置３に出力する。これにより、出力装置３において、例えば、四肢動物が健康異常であること、四肢動物に分娩の予兆があること、四肢動物に投薬の影響があること等に係る、音声出力、振動、および／または画面表示等の出力が可能となる。
【実施例0137】
以下、実施例を挙げて本発明をより具体的に説明するが、本発明はもとより下記実施例によって制限を受けるものではなく、前・後記の趣旨に適合し得る範囲で適当に変更を加えて実施することも勿論可能であり、それらはいずれも本発明の技術的範囲に包含される。
【0138】
先ず、以下の四肢動物の生体信号計測用衣類（ベルト）を作製した。
当該衣類は、四肢動物の胴体周りに配置される帯状本体部と、帯状本体部の肌側の面に形成された電極とを有し、帯状本体部は、電極形成領域において帯状本体部１ｃｍ2当たりに対して帯状本体部の厚み方向に５ｋｇの荷重をかけたときに、厚みが１０ｍｍ以上小さくなる発泡体（荷重をかけない時の厚み１１～５０ｍｍ、平均摩擦係数ＭＩＵ０．４０以上の摩擦部を有する）を有していた。
当該衣類は、四肢動物の胴体周りに配置される帯状本体部と、帯状本体部の肌側の面に形成された電極とを有し、帯状本体部は、電極形成領域において帯状本体部１ｃｍ2当たりに対して帯状本体部の厚み方向に５ｋｇの荷重をかけたときに、厚みが１０ｍｍ以上小さくなる発泡体（荷重をかけない時の厚み１１～５０ｍｍ、平均摩擦係数ＭＩＵ０．４０以上の摩擦部を有する）を有していた。
【0139】
次に、当該ベルト状の生体信号計測用衣類を、妊娠している黒毛和種の牛に装着させた。上記ベルトには、小型心拍計ｍｙＢｅａｔ（登録商標）（ユニオンツール株式会社製）が装着されており、分娩予定牛の平均心拍数（上述した単位時間当たりの拍数に相当）と加速度を計測した。
【0140】
具体的には、分娩予定日の３０日前の加速度の値のばらつきの標準偏差が０．０３以下となった時点を第１時間期間の始点、当該始点から１５分後を第１時間期間の終点とする、加速度の値のばらつきの標準偏差が０．０３以下である第１の時間期間を特定した。ここでも加速度として上述した無次元量の値を用いた。第１の時間期間では、上記生体信号計測用衣類を装着させた牛が、身体の左側を下とする臥位だったことが確認できた。第１の時間期間における平均心拍数を算出したところ、平均心拍数が８２ｂｐｍであった。
【0141】
上記第１の時間期間の終点の３０日後、加速度の値のばらつきの標準偏差が０．０３以下となった時点を第２の時間期間の始点、当該始点から１０分後を第２の時間期間の終点とする、加速度の値のばらつきの標準偏差が０．０３以下である第２の時間期間を特定した。第２の時間期間では、上記牛が、身体の右側を下側とする臥位だったことが確認できた。第２の時間期間における平均心拍数を算出したところ、平均心拍数が９１ｂｐｍであった。この時、分娩の兆候がみられると推測し、監視を始めた結果、２４時間以内に分娩が始まり安全な出産に至った。これにより、条件式（ｃ２）の有用性が確認された。
【0142】
［第１変形例］
上記第１の実施形態では、加速度のばらつきの度合いを示す値が閾値以下の時間期間が特定されるものとして説明した。しかしながら、本実施形態はこれに限定されない。
上記で説明した条件に加え、加速度および／または角速度の値が別の閾値以下である条件を満たす時間期間が特定されるようにしてもよい。当該加速度としては、上述した加速度Ａが用いられることが好ましい。これにより、四肢動物が臥位である時間期間が特定される精度が向上され得る。
上記第１の実施形態では、加速度のばらつきの度合いを示す値が閾値以下の時間期間が特定されるものとして説明した。しかしながら、本実施形態はこれに限定されない。
上記で説明した条件に加え、加速度および／または角速度の値が別の閾値以下である条件を満たす時間期間が特定されるようにしてもよい。当該加速度としては、上述した加速度Ａが用いられることが好ましい。これにより、四肢動物が臥位である時間期間が特定される精度が向上され得る。
【0143】
［第２変形例］
上述したように、身体の左側を下とする臥位と、身体の右側を下とする臥位とのような、異なる姿勢での臥位の間では、加速度が異なる値で安定することがある。例えば、第２の時間期間として、第１の時間期間の間と同じ姿勢での臥位に四肢動物があると推測される時間期間が特定されてもよい。これは、上記で説明した条件に加え、例えば、第２の時間期間の候補の時間期間での加速度と第１の時間期間での加速度との比較に基づいて実現され得る。例えば、第１の時間期間での加速度の平均値と、上記候補の時間期間での加速度の平均値との差が、或る範囲内にある場合に、当該候補の時間期間を第２の時間期間として特定してもよい。当該範囲は、例えば、複数種の臥位それぞれで安定する加速度に基づいて予め設定され得る。これにより、四肢動物が同じ姿勢での臥位にあると推測される、近い条件の２つの時間期間について、条件式が満足されるか否かが判定されることとなる。このため、条件式を用いることによる、四肢動物が健康異常であること、四肢動物に分娩の予兆があること、または四肢動物に投薬の影響があることの推測の精度が向上され得る。
上述したように、身体の左側を下とする臥位と、身体の右側を下とする臥位とのような、異なる姿勢での臥位の間では、加速度が異なる値で安定することがある。例えば、第２の時間期間として、第１の時間期間の間と同じ姿勢での臥位に四肢動物があると推測される時間期間が特定されてもよい。これは、上記で説明した条件に加え、例えば、第２の時間期間の候補の時間期間での加速度と第１の時間期間での加速度との比較に基づいて実現され得る。例えば、第１の時間期間での加速度の平均値と、上記候補の時間期間での加速度の平均値との差が、或る範囲内にある場合に、当該候補の時間期間を第２の時間期間として特定してもよい。当該範囲は、例えば、複数種の臥位それぞれで安定する加速度に基づいて予め設定され得る。これにより、四肢動物が同じ姿勢での臥位にあると推測される、近い条件の２つの時間期間について、条件式が満足されるか否かが判定されることとなる。このため、条件式を用いることによる、四肢動物が健康異常であること、四肢動物に分娩の予兆があること、または四肢動物に投薬の影響があることの推測の精度が向上され得る。
【0144】
以上、この発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる
実施形態に亘る構成要素を適宜組み合せてもよい。
実施形態に亘る構成要素を適宜組み合せてもよい。
【符号の説明】
【0145】
１…生体信号処理装置、１１…制御部、１１１…取得部、１１２…時間期間特定部、１１２１…第１時間期間特定部、１１２２…第２時間期間特定部、１１３…生体信号処理部、１１３１…判定部、１１４…出力部、１２…プログラム記憶部、１３…データ記憶部、１３１…加速度データ記憶部、１３２…生体信号記憶部、１３３…時間期間記憶部、１３４…生体信号処理結果記憶部、１４…入出力インタフェース、２…計測装置、３…出力装置、ＳＹＳ…生体信号処理システム、ＮＷ…通信ネットワーク、ＢＵＳ…バス。
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
