JPWO2010134525A1 - 飲酒検知システム及びコンピュータプログラム - Google Patents
飲酒検知システム及びコンピュータプログラム Download PDFInfo
- Publication number
- JPWO2010134525A1 JPWO2010134525A1 JP2011514424A JP2011514424A JPWO2010134525A1 JP WO2010134525 A1 JPWO2010134525 A1 JP WO2010134525A1 JP 2011514424 A JP2011514424 A JP 2011514424A JP 2011514424 A JP2011514424 A JP 2011514424A JP WO2010134525 A1 JPWO2010134525 A1 JP WO2010134525A1
- Authority
- JP
- Japan
- Prior art keywords
- frequency
- time
- drinking
- series waveform
- time series
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000035622 drinking Effects 0.000 title claims abstract description 288
- 238000001514 detection method Methods 0.000 title claims description 38
- 238000004590 computer program Methods 0.000 title claims description 16
- 230000010365 information processing Effects 0.000 claims abstract description 43
- 238000004364 calculation method Methods 0.000 claims description 90
- 238000012731 temporal analysis Methods 0.000 claims description 50
- 238000000700 time series analysis Methods 0.000 claims description 50
- 238000004458 analytical method Methods 0.000 claims description 42
- 230000008859 change Effects 0.000 claims description 41
- 230000033001 locomotion Effects 0.000 claims description 41
- 230000007423 decrease Effects 0.000 claims description 20
- 238000005259 measurement Methods 0.000 claims description 19
- 230000001965 increasing effect Effects 0.000 claims description 8
- 230000003247 decreasing effect Effects 0.000 claims description 3
- 239000004744 fabric Substances 0.000 description 88
- 239000007787 solid Substances 0.000 description 66
- 239000006260 foam Substances 0.000 description 42
- 239000011324 bead Substances 0.000 description 38
- 206010006326 Breath odour Diseases 0.000 description 25
- 229920005989 resin Polymers 0.000 description 25
- 239000011347 resin Substances 0.000 description 25
- 238000000034 method Methods 0.000 description 24
- 238000012360 testing method Methods 0.000 description 22
- 238000010586 diagram Methods 0.000 description 21
- -1 polyethylene terephthalate Polymers 0.000 description 19
- LFQSCWFLJHTTHZ-UHFFFAOYSA-N Ethanol Chemical compound CCO LFQSCWFLJHTTHZ-UHFFFAOYSA-N 0.000 description 13
- 239000000835 fiber Substances 0.000 description 12
- 239000000463 material Substances 0.000 description 10
- 238000002474 experimental method Methods 0.000 description 9
- 239000012528 membrane Substances 0.000 description 9
- 229920000139 polyethylene terephthalate Polymers 0.000 description 9
- 239000005020 polyethylene terephthalate Substances 0.000 description 9
- 210000003205 muscle Anatomy 0.000 description 8
- 230000002889 sympathetic effect Effects 0.000 description 7
- 210000005036 nerve Anatomy 0.000 description 6
- 239000004745 nonwoven fabric Substances 0.000 description 6
- 230000006870 function Effects 0.000 description 5
- 229920002215 polytrimethylene terephthalate Polymers 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 210000004177 elastic tissue Anatomy 0.000 description 4
- 238000009499 grossing Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000035945 sensitivity Effects 0.000 description 4
- 238000001228 spectrum Methods 0.000 description 4
- 238000012795 verification Methods 0.000 description 4
- 239000004698 Polyethylene Substances 0.000 description 3
- 210000000709 aorta Anatomy 0.000 description 3
- 210000002376 aorta thoracic Anatomy 0.000 description 3
- 230000002238 attenuated effect Effects 0.000 description 3
- 230000000747 cardiac effect Effects 0.000 description 3
- 238000005187 foaming Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000010355 oscillation Effects 0.000 description 3
- 229920000573 polyethylene Polymers 0.000 description 3
- 229920003225 polyurethane elastomer Polymers 0.000 description 3
- 238000012545 processing Methods 0.000 description 3
- 230000000644 propagated effect Effects 0.000 description 3
- 239000004743 Polypropylene Substances 0.000 description 2
- 239000004793 Polystyrene Substances 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 2
- 210000001217 buttock Anatomy 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 210000002837 heart atrium Anatomy 0.000 description 2
- 230000008035 nerve activity Effects 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 229920000728 polyester Polymers 0.000 description 2
- 239000011112 polyethylene naphthalate Substances 0.000 description 2
- 229920001155 polypropylene Polymers 0.000 description 2
- 229920002223 polystyrene Polymers 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000000284 resting effect Effects 0.000 description 2
- 239000000523 sample Substances 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 210000001765 aortic valve Anatomy 0.000 description 1
- 210000001367 artery Anatomy 0.000 description 1
- 230000001746 atrial effect Effects 0.000 description 1
- 235000013405 beer Nutrition 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 230000017531 blood circulation Effects 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 230000001447 compensatory effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 235000013305 food Nutrition 0.000 description 1
- 230000036571 hydration Effects 0.000 description 1
- 238000006703 hydration reaction Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 238000010030 laminating Methods 0.000 description 1
- WABPQHHGFIMREM-UHFFFAOYSA-N lead(0) Chemical compound [Pb] WABPQHHGFIMREM-UHFFFAOYSA-N 0.000 description 1
- 239000010985 leather Substances 0.000 description 1
- 239000002649 leather substitute Substances 0.000 description 1
- 210000004705 lumbosacral region Anatomy 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 208000010125 myocardial infarction Diseases 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000001734 parasympathetic effect Effects 0.000 description 1
- 229920003023 plastic Polymers 0.000 description 1
- 239000004033 plastic Substances 0.000 description 1
- 239000002985 plastic film Substances 0.000 description 1
- 229920006255 plastic film Polymers 0.000 description 1
- 229920003207 poly(ethylene-2,6-naphthalate) Polymers 0.000 description 1
- 229920006267 polyester film Polymers 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 239000002994 raw material Substances 0.000 description 1
- 239000012508 resin bead Substances 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 230000035807 sensation Effects 0.000 description 1
- 235000011888 snacks Nutrition 0.000 description 1
- 210000002784 stomach Anatomy 0.000 description 1
- 229920003002 synthetic resin Polymers 0.000 description 1
- 239000000057 synthetic resin Substances 0.000 description 1
- KKEYFWRCBNTPAC-UHFFFAOYSA-L terephthalate(2-) Chemical compound [O-]C(=O)C1=CC=C(C([O-])=O)C=C1 KKEYFWRCBNTPAC-UHFFFAOYSA-L 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 229920001169 thermoplastic Polymers 0.000 description 1
- 229920002725 thermoplastic elastomer Polymers 0.000 description 1
- 239000004416 thermosoftening plastic Substances 0.000 description 1
- 239000002759 woven fabric Substances 0.000 description 1
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6887—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
- A61B5/6892—Mats
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
- A61B5/02444—Details of sensor
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/48—Other medical applications
- A61B5/4845—Toxicology, e.g. by detection of alcohol, drug or toxic products
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/68—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient
- A61B5/6887—Arrangements of detecting, measuring or recording means, e.g. sensors, in relation to patient mounted on external non-worn devices, e.g. non-medical devices
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2503/00—Evaluating a particular growth phase or type of persons or animals
- A61B2503/20—Workers
- A61B2503/22—Motor vehicles operators, e.g. drivers, pilots, captains
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B2562/00—Details of sensors; Constructional details of sensor housings or probes; Accessories for sensors
- A61B2562/16—Details of sensor housings or probes; Details of structural supports for sensors
- A61B2562/168—Fluid filled sensor housings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1102—Ballistocardiography
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W2540/00—Input parameters relating to occupants
- B60W2540/24—Drug level, e.g. alcohol
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Physics & Mathematics (AREA)
- Veterinary Medicine (AREA)
- Biophysics (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- General Health & Medical Sciences (AREA)
- Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Cardiology (AREA)
- Toxicology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Pharmacology & Pharmacy (AREA)
- Physiology (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Auxiliary Drives, Propulsion Controls, And Safety Devices (AREA)
Abstract
飲酒の有無の判定をより正確に行うようにする。エアパックにより検出した背部の脈波の周波数に関する時系列変動の傾向を求める周波数動的情報処理手段６１０と、この周波数動的情報処理手段６１０により得られた周波数に関する時系列変動の傾向が非飲酒時における周波数に関する時系列変動の傾向と乖離している場合に飲酒状態と判定する飲酒状態判定手段６５０を有している。非飲酒時における周波数に関する時系列変動と比較して飲酒状態か否かを判定する構成であって、人の体調によって変化する脈波の周波数を周波数解析するだけでなく、その時系列変動を用いて判定する構成であるため、従来よりも正確に飲酒の有無を判定することができる。
Description
本発明は、呼気ではなく、人の動脈の脈波を用いて飲酒状態を検知する飲酒検知システム及びコンピュータプログラムに関する。
飲酒運転を防止する手段として、近年、アルコール・インターロック装置を自動車に搭載することが試みられている。しかしながら、呼気による飲酒検知装置では、同乗者の呼気によってごまかす懸念もある。
一方、本出願人は、運転中の運転者の生体状態を非侵襲で監視するシステムとして、例えば、特許文献１〜３において、シートクッション部に圧力センサを配置し、臀部脈波を採取して分析し、人の状態を推定する手法を開示している。具体的には、脈波の時系列信号を、それぞれ、ＳａｖｉｔｚｋｙとＧｏｌａｙによる平滑化微分法により、極大値と極小値を求める。そして、５秒ごとに極大値と極小値を切り分け、それぞれの平均値を求める。求めた極大値と極小値のそれぞれの平均値の差の二乗をパワー値とし、このパワー値を５秒ごとにプロットし、パワー値の時系列波形を作る。この時系列波形からパワー値の大域的な変化を読み取るために、ある時間窓Ｔｗ（１８０秒）について最小二乗法でパワー値の傾きを求める。次に、オーバーラップ時間Ｔｌ（１６２秒）で次の時間窓Ｔｗを同様に計算して結果をプロットする。この計算（移動計算）を順次繰り返してパワー値の傾きの時系列波形を得る。一方、脈波の時系列信号をカオス解析して最大リアプノフ指数を求め、上記と同様に、平滑化微分によって極大値を求め、移動計算することにより最大リアプノフ指数の傾きの時系列波形を得る。そして、パワー値の傾きの時系列波形と最大リアプノフ指数の傾きの時系列波形が逆位相となっており、さらには、パワー値の傾きの時系列波形で低周波、大振幅の波形が生じている波形を、入眠予兆を示す特徴的な信号と判定し、その後に振幅が小さくなったポイントを入眠点と判定している。
また、特許文献４として、内部に三次元立体編物を挿入した空気袋（エアパック）を備え、このエアパックを人の腰部に対応する部位に配置し、エアパックの空気圧変動を測定し、得られた空気圧変動の時系列データから人の生体信号を検出し、人の生体の状態を分析するシステムを開示している。また、非特許文献１及び２においても、腰腸肋筋に沿うようにエアパックセンサを配置して人の生体信号を検出する試みが報告されている。腰部付近の脈波は、心拍に伴う下行大動脈を流れる血流の循環変動を示すものであり、特許文献１及び２の臀部脈波を利用する場合よりも、心拍変動に即した人の状態変化をより正確に捉えることができる。
そして本出願人は、非特許文献３として、特許文献１〜４、非特許文献１，２の技術をさらに発展させ、エアパックセンサから得られる脈波の周波数解析を行うことで飲酒の有無の判定が可能であることを報告している。
「非侵襲型センサによって測定された生体ゆらぎ信号の疲労と入眠予知への応用」、落合直輝（外６名）、第３９回日本人間工学会 中国・四国支部大会 講演論文集、平成１８年１１月２５日発行、発行所：日本人間工学会 中国・四国支部事務局
「非侵襲生体信号センシング機能を有する車両用シートの試作」、前田慎一郎（外４名）、第３９回日本人間工学会 中国・四国支部大会 講演論文集、平成１８年１１月２５日発行、発行所：日本人間工学会 中国・四国支部事務局
「飲酒状態の非侵襲センシングシステム」、小島重行（外１０名）、自動車技術春季学術講演会前刷集、Ｎｏ．３７−０７、１５−１８，２００７年５月２３日発行、発行所：社団法人自動車技術会
非特許文献３は、上記のように、エアパックセンサから得られる背部の大動脈の脈波の時系列波形を周波数解析して飲酒の有無の判定が可能であることを示したが、より精度よく飲酒判定を行うことが望ましい。
本発明は上記に鑑みなされたものであり、人が飲酒したか否かの検知をより確実に行うことができる飲酒検知システム及びコンピュータプログラムを提供することを課題とする。
上記課題を解決するため、本発明は、人の胴部に対応して配置されるエアパックの空気圧変動を検出するセンサから得られる生体信号を解析して飲酒の有無を推定する飲酒解析推定部を備えた飲酒検知システムであって、前記飲酒解析推定部が、前記生体信号の時系列波形からその周波数に関する時系列変動を求める周波数動的情報処理手段と、前記周波数動的情報処理手段により得られた周波数に関する時系列変動の傾向が非飲酒時における周波数に関する時系列変動の傾向と乖離している場合に飲酒状態と判定する飲酒状態判定手段とを具備することを特徴とする飲酒検知システムを提供する。
前記周波数動的情報処理手段は、前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、前記飲酒状態判定手段は、前記卓越周波数時系列波形が非飲酒時よりも高いか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、前記飲酒状態判定手段は、所定時間範囲における前記卓越周波数時系列波形の卓越周波数の値のバラツキ度合いが、非飲酒時よりも拡大しているか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段を有し、前記飲酒状態判定手段は、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段と、前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記卓越周波数の傾きの時系列変化を卓越周波数傾き時系列波形として出力する卓越周波数傾き時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記卓越周波数傾き時系列解析演算手段により得られる卓越周波数傾き時系列波形の基線位置に対して、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも乖離しているか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否か、及び、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定する構成とすることができる。
また、本発明は、人の胴部の脈波を検出可能な部位に対応して配置されるエアパックと、前記エアパックの空気圧変動を検出するセンサとを備えた生体信号測定装置から得られる生体信号の時系列信号を分析して飲酒の有無を推定する飲酒検知システムの記憶部に設定される飲酒解析推定部を構成するコンピュータプログラムであって、前記飲酒解析推定部が、前記生体信号の時系列波形からその周波数に関する時系列変動を求める周波数動的情報処理手段と、前記周波数動的情報処理手段により得られた周波数に関する時系列変動の傾向が非飲酒時における周波数に関する時系列変動の傾向と乖離している場合に飲酒状態と判定する飲酒状態判定手段とを具備することを特徴とするコンピュータプログラムを提供する。
前記周波数動的情報処理手段は、前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、前記飲酒状態判定手段は、前記卓越周波数時系列波形が非飲酒時よりも高いか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、前記飲酒状態判定手段は、所定時間範囲における前記卓越周波数時系列波形の卓越周波数の値のバラツキ度合いが、非飲酒時よりも拡大しているか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段を有し、前記飲酒状態判定手段は、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段と、前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記卓越周波数の傾きの時系列変化を卓越周波数傾き時系列波形として出力する卓越周波数傾き時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記卓越周波数傾き時系列解析演算手段により得られる卓越周波数傾き時系列波形の基線位置に対して、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも乖離しているか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより判定する構成とすることができる。
前記周波数動的情報処理手段は、前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段と、前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段とを有し、前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否か、及び、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定する構成とすることができる。
本発明によれば、生体信号測定装置により人の上体より検出した生体信号の周波数に関する時系列変動の傾向を求め、非飲酒時における周波数に関する時系列変動と比較して飲酒状態か否かを判定する構成である。人の体調によって変化する生体信号の周波数を周波数解析するだけでなく、その時系列変動を用いて判定する構成であるため、従来よりも正確に飲酒の有無を判定することができる。
以下、図面に示した本発明の実施形態に基づき、本発明をさらに詳細に説明する。図１は、本実施形態に係る飲酒検知システム６０の分析対象である心部揺動波（人の上体の背部から検出される心房の動き及び大動脈の揺動に伴う生体信号）を採取する生体信号測定装置１を組み込んだ自動車用のシート５００の外観を示した図である。この図に示したように、生体信号測定装置１は、シートバック部５１０に組み込まれて用いられる。ここで、生体信号測定装置１によって採取される信号には、生体信号成分以外のノイズ信号がより少ないことが望ましい。そこで、本実施形態の生体信号測定装置１は、以下に説明するように、自動車の走行中等の振動環境下においても、センサの出力信号自体に含まれるノイズ信号を少なくできるような工夫がなされている。
生体信号測定装置１は、エアパックユニット１００と、第１のビーズ発泡樹脂弾性部材２０と、第２のビーズ発泡樹脂弾性部材３０とを有して構成されている。エアパックユニット１００は、収容体１５と、該収容体１５に収容した２つのエアパック１０を備えて構成される。各エアパック１０は、図３及び図４に示したように、表側エアパック１１と裏側エアパック１２とが積層されて構成され、収容体１５の左右にそれぞれ配置される。表側エアパック１１は、３つの小空気袋１１１が縦方向に連接されている一方、そのそれぞれは空気の流通がないように形成されている。各小空気袋１１１内には、復元力付与部材としての三次元立体編物１１２が配置されている。
裏側エアパック１２は、３つの小空気袋１１１を連接してなる表側エアパック１１の全長と同じ長さの大空気袋１２１とこの大空気袋１２１内に収容される復元力付与部材としての三次元立体編物１２２とを備えて構成される(図４参照）。表側エアパック１１と裏側エアパック１２とは、長手方向に沿った一方の側縁同士が接合され、接合された側縁を中心にして２つ折りにされて、相互に重ね合わせられて用いられる(図３（ｄ）及び図４参照)。
本実施形態では、このように表側エアパック１１と裏側エアパック１２とが相互に重ね合わせられたエアパック１０が左右に配置される。左右に配置することにより、着座者の背への当たりが左右均等になり、違和感を感じにくくなる。また、左右の表側エアパック１１，１１のいずれか一方を構成するいずれかの小空気袋１１１にセンサ取付チューブ１１１ａが設けられ、その内側に空気圧変動を測定するセンサ１１１ｂが固定されている。なお、センサ取付チューブ１１１ａは密閉されている。裏側エアパック１２を構成する大空気袋１２１にセンサを配設することもできるが、容量の大きい空気袋に設けると、生体信号による空気圧変動が吸収されてしまう場合があるため、小空気袋１１１に設けることが好ましい。但し、図４に示したように、予め、大空気袋１２１に取付チューブ１２１ａを設けその部位にセンサを配設しておき、必要に応じて、大空気袋１２１の空気圧変動を測定することで、小空気袋１１１の測定結果の検証に利用できるようにしておいてもよい。小空気袋１１１は、このような生体信号による空気圧変動に敏感に反応させるために、大きさは、幅４０〜１００ｍｍ、長さ１２０〜２００ｍｍの範囲が好ましい。小空気袋１１１の素材は限定されるものではないが、例えば、ポリウレタンエラストマー（例えば、シーダム株式会社製、品番「ＤＵＳ６０５−ＣＤＲ」）からなるシートを用いて形成することができる。センサ１１１ｂとしては、小空気袋１１１内の空気圧を測定できるものであればよく、例えば、コンデンサ型マイクロフォンセンサを用いることができる。
大空気袋１２１の大きさ及び小空気袋１１１を３つ連接した場合の全体の大きさとしては、自動車のシート５００のシートバック部５１０に用いる場合、幅４０〜１００ｍｍ、全長４００〜６００ｍｍの範囲とすることが好ましい。長さが短い場合、シートバック部５１０において、着座者が、腰部付近の一部分のみに異物感を感じるため、４００ｍｍ以上の長さとして、できるだけ、着座者の背全体に対応させることが好ましい。
空気圧変動を検出するセンサ１１１ｂは、本実施形態では、着座者の左側に配置されるエアパック１０を構成する表側エアパック１１の中央の小空気袋１１１に設けている。この小空気袋１１１の位置は、着座者の背部から採取される心房の動き及び大動脈（特に、「下行大動脈」）の揺動に伴う生体信号（心部揺動波）を検知可能な領域に相当する。この心部揺動波を検知可能な領域は、着座者の体格により一律ではないが、身長１５８ｃｍの日本人女性から身長１８５ｃｍの日本人男性までの様々な体格の被験者２０名で測定したところ、該小空気袋１１１（幅６０ｍｍ、長さ１６０ｍｍ）をシートバック部５１０の中心寄りの側縁と下縁の交差部Ｐ(図２及び図３参照）が、シートクッション部５２０の上面からシートバック部５１０の表面に沿った長さＬ：２２０ｍｍ、シートバック部５１０の中心からの距離Ｍ：８０ｍｍとなるように設定したところ、上記全ての被験者において心部揺動波を検知できた。小空気袋１１１の大きさが、幅４０〜１００ｍｍ、長さ１２０〜２００ｍｍの範囲の場合、交差部Ｐの位置を、シートクッション部５２０の上面からシートバック部５１０の表面に沿った長さで１５０〜２８０ｍｍ、シートバック部５１０の中心から６０〜１２０ｍｍの範囲に設定することが好ましい。
上記した２つのエアパック１０をシートバック部５１０において容易に所定の位置に設定できるようにユニット化しておくことが好ましい。従って、図２〜図４に示したような収容体１５にエアパック１０を装填したエアパックユニット１００として構成とすることが好ましい。収容体１５は、両側にエアパック１０を収容する袋状のエアパック収容部１５１を有し、２つのエアパック収容部１５１間に接続部１５２を有している。
２つのエアパック収容部１５１には、それぞれエアパック１０が挿入される。また、エアパック収容部１５１には、エアパック１０とほぼ同じ大きさの三次元立体編物４０を、エアパック１０の裏側エアパック１２の背面側に重ねて挿入することが好ましい（図３（ｄ）参照）。三次元立体編物４０を配置することにより、エアパック１０が該三次元立体編物４０によっていわば浮くように支持されるため、シートバック部５１０からの外部振動が伝わりにくくなる。すなわち、三次元立体編物４０を配置することにより、高周波小振幅の外部振動が入力された場合には、三次元立体編物４０のパイルと空気圧の変動から、エアパック内にバネ定数の低い、バネ・マス・ダンパ系が作られる。そして、それが三次元立体編物４０を内蔵したエアパック１０において、低・高周波入力に対するフィルタ（ローパスフィルタ・ハイパスフィルタ）として作用し、該外部振動を減衰する。
接続部１５２は、２つのエアパック部１５１を所定間隔をおいて支持できるものであればよく、幅６０〜１２０ｍｍ程度で形成される。接続部１５２も、袋状に形成し、その内部に三次元立体編物４５を挿入することが好ましい(図３（ｄ）及び図４参照）。これにより、該接続部１５２を通じて入力される振動も、該三次元立体編物４５を挿入することにより効果的に除振でき、センサ１１１ｂを備えたエアパック１０への外部振動の伝達を抑制できる。
なお、上記したように、小空気袋１１１は、例えば、ポリウレタンエラストマー（例えば、シーダム株式会社製、品番「ＤＵＳ６０５−ＣＤＲ」）からなるシートを用いて形成されるが、裏側エアパック１２を形成する大空気袋１２１及び収容体１５も、同じ素材を用いて形成することが好ましい。また、小空気袋１１１、大空気袋１２１、エアパック収容部１５１及び接続部１５２内に装填される各三次元立体編物は、例えば、特開２００２−３３１６０３号公報に開示されているように、互いに離間して配置された一対のグランド編地と、該一対のグランド編地間を往復して両者を結合する多数の連結糸とを有する立体的な三次元構造となった編地である。
一方のグランド編地は、例えば、単繊維を撚った糸から、ウェール方向及びコース方向のいずれの方向にも連続したフラットな編地組織（細目）によって形成され、他方のグランド編地は、例えば、短繊維を撚った糸から、ハニカム状（六角形）のメッシュを有する編み目構造に形成されている。もちろん、この編地組織は任意であり、細目組織やハニカム状以外の編地組織を採用することもできるし、両者とも細目組織を採用するなど、その組み合わせも任意である。連結糸は、一方のグランド編地と他方のグランド編地とが所定の間隔を保持するように、２つのグランド編地間に編み込んだものである。このような三次元立体編物としては、例えば、以下のようなものを用いることができる。なお、各三次元立体編物は、必要に応じて複数枚積層して用いることもできる。
（１）製品番号：４９０７６Ｄ（住江織物（株）製）
材質：
表側のグランド編地・・・３００デシテックス／２８８ｆのポリエチレンテレフタレート繊維仮撚加工糸と７００デシテックス／１９２ｆのポリエチレンテレフタレート繊維仮撚加工糸との撚り糸
裏側のグランド編地・・・４５０デシテックス／１０８ｆのポリエチレンテレフタレート繊維仮撚加工糸と３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメント
材質：
表側のグランド編地・・・３００デシテックス／２８８ｆのポリエチレンテレフタレート繊維仮撚加工糸と７００デシテックス／１９２ｆのポリエチレンテレフタレート繊維仮撚加工糸との撚り糸
裏側のグランド編地・・・４５０デシテックス／１０８ｆのポリエチレンテレフタレート繊維仮撚加工糸と３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメント
（２）製品番号：４９０１１Ｄ（住江織物（株）製）
材質：
グランド編地（縦糸）・・・６００デシテックス／１９２ｆのポリエチレンテレフタレート繊維仮撚加工糸
グランド編地（横糸）・・・３００デシテックス／７２ｆのポリエチレンテレフタレート繊維仮撚加工糸
連結糸・・・・・・・・・８００デシテックス／１ｆのポリエチレンテレフタレートモノフィラメント
材質：
グランド編地（縦糸）・・・６００デシテックス／１９２ｆのポリエチレンテレフタレート繊維仮撚加工糸
グランド編地（横糸）・・・３００デシテックス／７２ｆのポリエチレンテレフタレート繊維仮撚加工糸
連結糸・・・・・・・・・８００デシテックス／１ｆのポリエチレンテレフタレートモノフィラメント
（３）製品番号：４９０１３Ｄ（住江織物（株）製）
材質：
表側のグランド編地・・・４５０デシテックス／１０８ｆのポリエチレンテレフタレート繊維仮撚加工糸の２本の撚り糸
裏側のグランド編地・・・４５０デシテックス／１０８ｆのポリエチレンテレフタレート繊維仮撚加工糸の２本の撚り糸
連結糸・・・・・・・・・３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメント
材質：
表側のグランド編地・・・４５０デシテックス／１０８ｆのポリエチレンテレフタレート繊維仮撚加工糸の２本の撚り糸
裏側のグランド編地・・・４５０デシテックス／１０８ｆのポリエチレンテレフタレート繊維仮撚加工糸の２本の撚り糸
連結糸・・・・・・・・・３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメント
（４）製品番号：６９０３０Ｄ（住江織物（株）製）
材質：
表側のグランド編地・・・４５０デシテックス／１４４ｆのポリエチレンテレフタレート繊維仮撚加工糸の２本の撚り糸
裏側のグランド編地・・・４５０デシテックス／１４４ｆのポリエチレンテレフタレート繊維仮撚加工糸と３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメント
材質：
表側のグランド編地・・・４５０デシテックス／１４４ｆのポリエチレンテレフタレート繊維仮撚加工糸の２本の撚り糸
裏側のグランド編地・・・４５０デシテックス／１４４ｆのポリエチレンテレフタレート繊維仮撚加工糸と３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメントとの組み合わせ
連結糸・・・・・・・・・３５０デシテックス／１ｆのポリトリメチレンテレフタレートモノフィラメント
（５）旭化成せんい（株）製の製品番号：Ｔ２４０５３ＡＹ５−１Ｓ
第１のビーズ発泡樹脂弾性部材２０と第２のビーズ発泡樹脂弾性部材３０とは、シートバック部５１０の表皮部材とエアパック１０を収容した収容体１５（エアパックユニット１００）との間に配設され、２つのエアパック１０の全長に相当する長さを有し、２つのエアパック１０の頂部間の長さに相当する幅を有している。従って、長さが４００〜６００ｍｍ、幅が２５０〜３５０ｍｍ程度の大きさのものを用いることが好ましい。これにより、２つのエアパック１０が共に覆われるため、２つのエアパック１０の凹凸を感じにくくなる。
第１のビーズ発泡樹脂弾性部材２０は、平板状に形成されたビーズ発泡体と、その外面に貼着される被覆材とから構成されている。ビーズ発泡体としては、ポリスチレン、ポリプロピレン及びポリエチレンのいずれか少なくとも一つを含む樹脂のビーズ法による発泡成形体が用いられる。なお、発泡倍率は任意であり限定されるものではない。被覆材は、ビーズ発泡体の外面に接着により貼着され、高い伸度と回復率を有する素材であり、好ましくは、伸度２００％以上、１００％伸長時の回復率が８０％以上である弾性繊維不織布が用いられる。例えば、特開２００７−９２２１７号公報に開示された熱可塑性エラストマー弾性繊維が相互に溶融接着された不織布を用いることができる。具体的には、ＫＢセーレン（株）製、商品名「エスパンシオーネ」を用いることができる。
第２のビーズ発泡樹脂弾性部材３０は、第１のビーズ発泡樹脂弾性部材２０と同様にビーズ発泡体を備えて構成されるが、その外面を覆う被覆材としては、第１のビーズ発泡樹脂弾性部材２０において用いた弾性繊維不織布よりも伸縮性の小さい素材、例えば、熱可塑性ポリエステルからなる不織布が用いられる。具体的には、帝人（株）製のポリエチレンナフタレート（ＰＥＮ）繊維（１１００ｄｔｅｘ）から形成した２軸織物（縦：２０本／ｉｎｃｈ、横：２０本／ｉｎｃｈ）を用いることができる。
第２のビーズ発泡樹脂弾性部材３０は、第１のビーズ発泡樹脂弾性部材２０と同様にビーズ発泡体を備えて構成されるが、その外面を覆う被覆材としては、第１のビーズ発泡樹脂弾性部材２０において用いた弾性繊維不織布よりも伸縮性の小さい素材、例えば、熱可塑性ポリエステルからなる不織布が用いられる。具体的には、帝人（株）製のポリエチレンナフタレート（ＰＥＮ）繊維（１１００ｄｔｅｘ）から形成した２軸織物（縦：２０本／ｉｎｃｈ、横：２０本／ｉｎｃｈ）を用いることができる。
第１のビーズ発泡樹脂弾性部材２０と第２のビーズ発泡樹脂弾性部材３０とを積層する順序は限定されるものではないが、シートバック部５１０の表皮部材５１１に近い側に、弾性の高い第１のビーズ発泡樹脂弾性部材２０を配設することが好ましい。また、第１及び第２のビーズ発泡樹脂弾性部材２０，３０を構成するビーズ発泡体は、厚さ約５〜６ｍｍ程度とし、その外面に、厚さ約１ｍｍ以下の上記した弾性繊維不織布や熱可塑性ポリエステルからなる不織布を貼着して形成される。なお、本実施形態では、第１のビーズ発泡樹脂弾性部材２０の表皮部材５１１に対向する面、第２のビーズ発泡樹脂弾性部材３０のエアパックユニット１００に対向する面に、それぞれＰＥＮフィルムなどのポリエステルフィルムを貼着している。これにより、生体信号の伝達性が向上する。
本実施形態において人体支持手段を構成するシート５００のシートバック部５１０は、表皮部材５１１と該表皮部材５１１の背面側に配設されるクッション支持部材５１２とを備えてなり、該表皮部材５１１とクッション支持部材５１２との間にエアパック１０を保持した収容体１５（エアパックユニット１００）と第１及び第２のビーズ発泡樹脂弾性部材２０，３０が組み込まれる。この際、クッション支持部材５１２側にまずエアパック１０を保持した収容体１５（エアパックユニット１００）が配置され、その表面側に第２のビーズ発泡樹脂弾性部材３０が、さらにその表面側に第１のビーズ発泡樹脂弾性部材２０が配置された上で、表皮部材５１１により被覆される。なお、クッション支持部材５１２は、例えば、三次元立体編物をシートバック部５１０の左右一対のサイドフレームの後端縁間に張って形成することもできるし、合成樹脂板から形成することもできる。表皮部材５１１は、例えば、三次元立体編物、合成皮革、皮革、あるいはこれらの積層体などを左右一対のサイドフレームの前縁間に張って設けることができる。
このように、本実施形態においては、表皮部材５１１の裏面側に所定の大きさの第１のビーズ発泡樹脂弾性部材２０及び第２のビーズ発泡樹脂弾性部材３０が積層して配置され、さらにその後方に左右一対のエアパック１０を保持した収容体１５（エアパックユニット１００）が配置される構成であるため、着座者が背にエアパック１０の凹凸を感じることなくなり、生体信号を測定するためのエアパック１０を有する構成でありながら、座り心地が向上する。
次に、飲酒検知システム６０の構成について図６に基づいて説明する。飲酒検知システム６０には、生体信号測定装置１により検出された生体信号である心部揺動波（以下、場合により「エアパック脈波」という）の時系列波形を解析して飲酒の有無を推定する飲酒解析推定部６００が組み込まれている。なお、本実施形態で用いた生体信号測定装置１は、上記のようにノイズ対策を施しているため、検出信号へのノイズの混入は少ないが、特に、自動車運転中のような動的環境下では検出信号に心部揺動波以外のノイズが含まれることが多くなる。従って、その場合には、飲酒解析推定部６００で処理する前の前処理として、検出信号を心部揺動波が含まれている所定周波数でフィルタリングするなどして処理し、この前処理した検出信号を心部揺動波（エアパック脈波）の時系列波形として用いることが好ましい。
本実施形態では、飲酒解析推定部６００を、飲酒検知システム６０の記憶部に設定したコンピュータプログラムから構成している。すなわち、飲酒解析推定部６００は、生体信号測定装置１から得られるエアパック脈波の時系列波形からその周波数に関する時系列変動を求める周波数動的情報処理手段６１０と、周波数動的情報処理手段６１０により得られた周波数に関する時系列変動の傾向が非飲酒時における周波数に関する時系列変動の傾向と乖離している場合に飲酒状態と判定する飲酒状態判定手段６５０とを具備する。
周波数動的情報処理手段６１０は、さらに、卓越周波数時系列波形演算手段６２１、卓越周波数変動時系列解析演算手段６２２、卓越周波数傾き時系列解析演算手段６２３、周波数演算手段６３１、周波数変動時系列解析演算手段６３２、周波数傾き時系列解析演算手段６３３を有して構成される。
なお、これらのコンピュータプログラムは、フレキシブルディスク、ハードディスク、ＣＤ−ＲＯＭ、ＭＯ（光磁気ディスク）、ＤＶＤ−ＲＯＭなどの記録媒体へ記憶させて提供することもできるし、通信回線を通じて伝送することも可能である。
卓越周波数時系列波形演算手段６２１は、エアパック脈波の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める。例えば、９０秒、３００秒といった時間窓におけるエアパック脈波の時系列波形を周波数解析して求める。この卓越周波数時系列波形は、図１１に示したように、３００秒の時間窓毎にオーバーラップさせることなく、卓越周波数を求めるようにしてもよいが（すなわち、図１１の場合、２０分間の測定時間において、４点の卓越周波数が求められる）、図１２及び図１３に示したように、例えば、９０秒の時間窓、３００秒の時間窓を設定して、その９０％の時間をオーバーラップさせて、それぞれの時間窓毎に卓越周波数を求めるようにしても良い。図１１の処理よりも、図１２及び図１３の処理の方が卓越周波数の変遷の傾向が反映されるという利点がある。
卓越周波数変動時系列解析演算手段６２２は、卓越周波数時系列波形演算手段６２１により得られる卓越周波数時系列波形に、所定の時間幅の時間窓（好ましくは１８０秒）を設定し、周波数の平均値を求める（図７の[3],[4]のステップ参照）。次に、所定のオーバーラップ時間（好ましくは１６２秒）で設定した所定の時間窓（好ましくは１８０秒）毎にその間に算出された卓越周波数の平均値を求める移動計算を行い、プロットする。そして、時間窓毎にプロットされた周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する（図７の[7]のステップ参照）。
卓越周波数傾き時系列解析演算手段６２３は、卓越周波数時系列波形演算手段６２１により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に卓越周波数の傾きを求める移動計算を行い、時間窓毎に得られる卓越周波数の傾きの時系列変化を卓越周波数傾き時系列波形として出力する。具体的には、まず、ある時間窓Ｔｗ１における周波数の傾きを最小二乗法により求めてプロットする（図７の[3],[5]のステップ参照）。次に、オーバーラップ時間Ｔｌ（図７の[6]のステップ）で次の時間窓Ｔｗ２を設定し、この時間窓Ｔｗ２における卓越周波数の傾きを同様に最小二乗法により求めてプロットする。この計算（移動計算）を順次繰り返し、エアパック脈波の卓越周波数の傾きの時系列変化を卓越周波数傾き時系列波形として出力する（図７の[8]のステップ参照）。なお、時間窓Ｔｗの時間幅は１８０秒に設定することが好ましく、オーバーラップ時間Ｔｌは１６２秒に設定することが好ましい。これは、本出願人による上記特許文献３（ＷＯ２００５／０９２１９３Ａ１公報）において示したように、時間窓Ｔｗの時間幅及びオーバーラップ時間Ｔｌを種々変更して行った睡眠実験から、人の状態を示す特徴的な信号波形が最も感度よく出現する値として選択されたものである。
周波数演算手段６３１は、生体信号測定装置１から得られるエアパック脈波の時系列波形における周波数の時系列波形を求める。具体的には、まず、エアパック脈波の時系列波形を平滑化微分して極大値（ピーク）を求める。例えば、ＳａｖｉｔｚｋｙとＧｏｌａｙによる平滑化微分法により極大値を求める。次に、５秒ごとに極大値を切り分け、その５秒間に含まれる時系列波形の極大値（波形の山側頂部）間の時間間隔の逆数を個別周波数ｆとして求め、その５秒間における個別周波数ｆの平均値を当該５秒間の周波数Ｆの値として採用する（図７の[1]のステップ参照）。そして、この５秒毎に得られる周波数Ｆをプロットすることにより、周波数の時系列波形を求める（図７の[2]のステップ参照）。
周波数変動時系列解析演算手段６３２は、周波数演算手段６３１により得られたエアパック脈波の周波数の時系列波形（図７の[2]のステップ参照）に、所定の時間幅の時間窓（好ましくは１８０秒）を設定し、周波数の平均値を求める（図７の[3],[4]のステップ参照）。次に、所定のオーバーラップ時間（好ましくは１６２秒）で設定した所定の時間窓（好ましくは１８０秒）毎にエアパック脈波の周波数の平均値を求める移動計算を行い、プロットする。そして、時間窓毎にプロットされた周波数の平均値の時系列変化を周波数変動時系列波形として出力する（図７の[7]のステップ参照）。
周波数時系列解析傾き演算手段６３３は、周波数演算手段６３１により得られたエアパック脈波の周波数の時系列波形から、所定の時間幅の時間窓を設定し、時間窓毎に最小二乗法により該エアパック脈波の周波数の傾きを求め、その時系列波形を出力する。具体的には、まず、ある時間窓Ｔｗ１における周波数の傾きを最小二乗法により求めてプロットする（図７の[3],[5]のステップ参照）。次に、オーバーラップ時間Ｔｌ（図７の[6]のステップ参照）で次の時間窓Ｔｗ２を設定し、この時間窓Ｔｗ２における周波数の傾きを同様に最小二乗法により求めてプロットする。この計算（移動計算）を順次繰り返し、エアパック脈波の周波数の傾きの時系列変化を周波数傾き時系列波形として出力する（図７の[8]のステップ参照）。
なお、周波数演算手段６３１は、図７のステップに従った処理（以下、「ピーク検出法」という）ではなく、図８のステップに従った処理（以下、「ゼロクロス法」という）で実行することができる。図７のピーク検出法は、周波数演算手段６３１が、生体信号測定装置１から得られるエアパック脈波の時系列波形を平滑化微分して極大値（ピーク値）を求めているが、図８のゼロクロス法では、エアパック脈波の時系列波形において、正から負に切り替わる地点（以下、「ゼロクロス地点」という）を求める。そして、５秒毎にゼロクロス地点を切り分け、その５秒間に含まれる時系列波形のゼロクロス地点間の時間間隔の逆数を個別周波数ｆとして求め、その５秒間における個別周波数ｆの平均値を当該５秒間の周波数Ｆの値として採用する（図８の[1]のステップ参照）。そして、この５秒毎に得られる周波数Ｆをプロットすることにより、周波数の時系列波形を求める（図８の[2]のステップ参照）。
その後は、図７の場合と同様に、周波数変動時系列解析演算手段６３２によって移動計算を行い、周波数変動時系列波形を求める（図８の[3],[4],[6],[7]のステップ参照）。また、周波数傾き時系列解析演算手段６３３によって移動計算を行い、周波数傾き時系列波形を求める（図８の[3],[5],[6],[8]のステップ参照）。
後述の飲酒状態判定手段６５０において、周波数傾き時系列波形、周波数変動時系列波形及びその基線を用いて人の状態を判定する場合、ピーク検出法、ゼロクロス法のいずれを用いてもよい。いずれがより人の状態を明確に示すかは個人差があるため、人により、適する方法を予め設定できるようにすることが好ましい。
飲酒状態判定手段６５０は、本実施形態では、上記した各種の周波数動的情報処理手段６１０により得られた周波数に関する時系列変動を基に、飲酒状態か否かを判定する以下のような複数のプログラムのうちの少なくとも一つから構成される。
（ａ）卓越周波数時系列波形演算手段６２１により得られる卓越周波数時系列波形が非飲酒時よりも高いか否かにより飲酒状態を判定する手段
（ｂ）所定時間範囲における卓越周波数時系列波形の卓越周波数の値のバラツキ度合いが、非飲酒時よりも拡大しているか否かにより飲酒状態を判定する手段
（ｃ）卓越周波数変動時系列解析演算手段６２２により得られる卓越周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する手段
（ｄ）卓越周波数傾き時系列解析演算手段６２３により得られる卓越周波数傾き時系列波形の基線位置に対して、卓越周波数変動時系列解析演算手段６２２により得られる卓越周波数変動時系列波形の基線位置が非飲酒時よりも乖離しているか否かにより飲酒状態を判定する手段
（ｅ）周波数変動時系列解析演算手段６３２により得られる周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する手段
（ｆ）周波数傾き時系列解析演算手段６３３により得られる周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定する手段
（ｂ）所定時間範囲における卓越周波数時系列波形の卓越周波数の値のバラツキ度合いが、非飲酒時よりも拡大しているか否かにより飲酒状態を判定する手段
（ｃ）卓越周波数変動時系列解析演算手段６２２により得られる卓越周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する手段
（ｄ）卓越周波数傾き時系列解析演算手段６２３により得られる卓越周波数傾き時系列波形の基線位置に対して、卓越周波数変動時系列解析演算手段６２２により得られる卓越周波数変動時系列波形の基線位置が非飲酒時よりも乖離しているか否かにより飲酒状態を判定する手段
（ｅ）周波数変動時系列解析演算手段６３２により得られる周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する手段
（ｆ）周波数傾き時系列解析演算手段６３３により得られる周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定する手段
飲酒状態判定手段６５０は、（ａ）〜（ｆ）の複数の判定手段を備えていることにより、いずれか一つにおいて飲酒状態と判定されたならば、その旨（飲酒状態であること）を出力するように設定すれば、飲酒の検知感度を高めることができる。その一方、（ａ）〜（ｆ）のうち、例えば、２つ以上の手段において飲酒状態と判定されたならば、その旨（飲酒状態であること）を出力するように設定することもできる。後者によれば、誤検知を少なくすることができるという利点がある。これは、本実施形態の飲酒検知システム６０を用いる用途により、例えば、診断用、飲酒運転検知用などにより、適宜の組み合わせで設定することが可能である。また、いずれも、予め、非飲酒時の測定結果（卓越周波数時系列波形、卓越周波数変動時系列波形、卓越周波数傾き時系列波形、周波数変動時系列波形、周波数傾き時系列波形、周波数傾き時系列波形から得られる正の傾きの積分値と負の傾きの積分値等）を個人毎に記憶させておく。そして、ある時点における測定結果が得られたならば、飲酒状態判定手段６５０が、記憶されている非飲酒時の測定結果を読み込み、当該ある時点における測定結果と比較する。それにより、（ａ）〜（ｆ）の各判定手段において、非飲酒時と乖離した傾向を示すと「飲酒状態」と判定できる。但し、どの程度乖離している場合に、「飲酒状態」と判定するかについては、上記（ａ）〜（ｆ）の各判定手段に対応して、個人毎に予め閾値を設定しておくことが好ましい（つまり、非飲酒時のデータと飲酒時のデータとを比較して、どの程度乖離していると飲酒状態となるかについてのデータを予めとり、閾値を定めておく）。
なお、上記実施形態においては、生体信号測定装置１を構成するエアパック１０、第１及び第２のビーズ発泡樹脂弾性部材２０，３０を自動車用のシートに組み込んでいるが、自動車用のシートに限らず、ベッドなどの寝具、病院設備における診断用の椅子等に組み込むこともできる。
（試験例１）
健康な２０歳代から４０歳代の日本人男性４名を被験者としてそれぞれ上記シート５００に着座させて実験を行った。なお、以下において、エアパック脈波は、上記シート５００に着座させて生体信号測定装置１（以下、場合により「エアパックセンサ」という）から得られた心部揺動波である。また、指尖容積脈波は、光学式指尖容積脈波計（（株）アコム製、フィンガークリッププローブ ＳＲ−５Ｃ）により測定して得られたものであり、呼気アルコール濃度は、「東海電子（株）製 ＡＬＣ−ｍｉｎｉ」により測定して得られたものである。被験者には，あらかじめ飲酒実験とは別の日にエタノールパッチテストを行い、活性型（ＮＮ型）であることを確認した。尚、被験者の体重および身長は、被験者Ａは体重７１ｋｇ、身長１６７ｃｍ、被験者Ｂは体重６８ｋｇ、身長１７８ｃｍ、被験者Ｃは体重６５ｋｇ、身長１７１ｃｍ、被験者Ｄは体重５９ｋｇ、身長１６６ｃｍの被験者である。
健康な２０歳代から４０歳代の日本人男性４名を被験者としてそれぞれ上記シート５００に着座させて実験を行った。なお、以下において、エアパック脈波は、上記シート５００に着座させて生体信号測定装置１（以下、場合により「エアパックセンサ」という）から得られた心部揺動波である。また、指尖容積脈波は、光学式指尖容積脈波計（（株）アコム製、フィンガークリッププローブ ＳＲ−５Ｃ）により測定して得られたものであり、呼気アルコール濃度は、「東海電子（株）製 ＡＬＣ−ｍｉｎｉ」により測定して得られたものである。被験者には，あらかじめ飲酒実験とは別の日にエタノールパッチテストを行い、活性型（ＮＮ型）であることを確認した。尚、被験者の体重および身長は、被験者Ａは体重７１ｋｇ、身長１６７ｃｍ、被験者Ｂは体重６８ｋｇ、身長１７８ｃｍ、被験者Ｃは体重６５ｋｇ、身長１７１ｃｍ、被験者Ｄは体重５９ｋｇ、身長１６６ｃｍの被験者である。
非侵襲センサであるエアパックセンサと他のセンサ（光学式指尖容積脈波計、呼気アルコール濃度計）による生体信号の計測は、飲酒前に２０分間の計測を１回行い、その後、飲酒（ビール：５００ｍｌ）をして、飲酒後に最も血液中のアルコール濃度が高くなるとされている２０〜４０分後の２０分間に１回目の計測を行った。その後、被験者Ａについては時間経過による変化を見るために、一定時間を空けて９０〜１１０分後、１６０〜１８０分後の２回の計４回測定を行った。被験者は飲酒開始時間より前、最低３時間以上食事をしていない状態とし、空腹状態に近い状態を再現した。また通常の飲酒状態を模擬するため、飲酒時に適量のおつまみを摂取した。また実験中には、これ以外には水分補給のみ行い、他のものの飲食は行わないようにした。生体信号の計測前後には、呼気アルコール濃度を計測した。その変動の様子を図９に示す。
図１０に被験者Ａの場合の指尖容積脈波とエアパック脈波（心部揺動波）のそれぞれ計測開始５分間分のデータでの原波形と５分間分の周波数解析結果を示す。指尖容積脈波とエアパック脈波の卓越周波数のピーク位置はほぼ一致しており、同等の周波数特性を有した生体信号が得られていることがわかる。
図１１は、卓越周波数時系列波形演算手段６２１により求められた被験者Ａのエアパック脈波の３００秒毎の時間窓をオーバーラップさせずに求めた卓越周波数の変遷を示す。図１１には、指尖容積脈波の卓越周波数の３００秒分ごとの変遷も示している。図９及び図１１から、飲酒による呼気アルコール濃度の上昇に対応してエアパック脈波の周波数の変動が認められ、飲酒により心拍数が上昇した状態が検出できる。そして時間の経過と共に、呼気アルコール濃度の低下が生じ、飲酒前の心拍数に近づいていることもわかる。また、図１１から、飲酒前と比較して、飲酒後２０〜４０分間、飲酒後９０〜１１０分間では卓越周波数のバラツキ度合いが拡大していることがわかる。飲酒後１６０〜１８０分間の場合も、飲酒前より卓越周波数のバラツキ度合いが若干拡大している。従って、この卓越周波数の値のバラツキ度合いが拡大するか否かを判定することにより、飲酒状態を判定することができることがわかる。なお、この点は、次に述べるように、図１２及び図１３の場合でも同様である。
図１２は、時間窓９０秒に設定して９０％のオーバーラップ時間で各時間窓毎に求めた、卓越周波数時系列波形演算手段６２１により得られた卓越周波数の変遷の詳細を示し、図１３は、時間窓３００秒に設定して９０％のオーバーラップ時間で各時間窓毎に求めた、卓越周波数時系列波形演算手段６２１により得られた卓越周波数の変遷の詳細を示す。
図９から、被験者Ａは、飲酒後２０〜４０分間に酒気帯び運転法定基準の境界値である０．１５ｍｇ／ｌ程度の呼気アルコール濃度が検出され、飲酒後９０〜１１０分間も０．０５ｍｇ／ｌからそれ以下の呼気アルコール濃度が検出されている。
図１２では、飲酒前の卓越周波数時系列波形が概ね１．４〜１．５５Ｈｚで推移し、飲酒後１６０〜１８０分間の卓越周波数時系列波形が概ね１．４５〜１．５５Ｈｚの間で推移しているのに対し、飲酒後２０〜４０分間は概ね１．５２〜１．７５Ｈｚ、飲酒後９０〜１１０分間は概ね１．４９〜１．７５Ｈｚと、卓越周波数時系列波形は全体的に高い傾向にある。従って、飲酒状態判定手段６５０は、この卓越周波数時系列波形を飲酒前と比較することにより、飲酒状態を判定することができることがわかる。
また、卓越周波数のバラツキ度合いは、飲酒前が０．１５Ｈｚであるのに対し、飲酒後２０〜４０分間が０．２３Ｈｚ、飲酒後９０〜１１０分間が０．２６Ｈｚと拡大している。飲酒後１６０〜１８０分間は０．１Ｈｚと再度バラツキ度合いが小さくなっている。従って、このバラツキ度合いの閾値を設定することにより飲酒状態を判定することができることがわかる。
図１３の場合には、飲酒前の卓越周波数時系列波形が概ね１．４４〜１．４６Ｈｚ（８００秒における突出値は体動によるものと思われる）（バラツキ度合い０．０２Ｈｚ）で推移しているのに対し、飲酒後２０〜４０分間は概ね１．５５〜１．６４Ｈｚ（バラツキ度合い０．０９Ｈｚ）、飲酒後９０〜１１０分間は概ね１．５５〜１．６２Ｈｚ（バラツキ度合い０．０７Ｈｚ）、飲酒後１６０〜１８０分間は概ね１．４４〜１．５２Ｈｚ（バラツキ度合い０．０８Ｈｚ）となっており、図１２の場合と同様に、飲酒により、卓越周波数時系列波形は全体的に高い傾向になると共に、バラツキ度合いが拡大する傾向にあることがわかる。
図１４（ａ），（ｂ）及び図１５（ａ），（ｂ）は、卓越周波数変動時系列解析演算手段６２２により得られる卓越周波数変動時系列波形と、卓越周波数傾き時系列解析演算手段６２３により得られる卓越周波数傾き時系列波形と出力したグラフであり、図１４は時間窓９０秒に設定して求めたグラフであり、図１５は時間窓３００秒に設定して求めたグラフである。これらのグラフから明らかなように、卓越周波数変動時系列波形は、その基線位置が、飲酒後２０〜４０分間及び飲酒後９０〜１１０分間は、いずれも概ね１．６Ｈｚ前後であり、飲酒前及び飲酒後１６０〜１８０分間の基線位置が概ね１．５Ｈｚ前後であると比較すると、明らかに高くなっていることがわかる。よって、例えば、閾値１．６Ｈｚに設定し、基線位置が１．６Ｈｚ以上に至ったならば、飲酒状態判定手段６５０が飲酒と判定するようにすることができる。なお、閾値を何Ｈｚに設定するかは、個人毎に設定できるようにすることが好ましい。
図１６（ａ），（ｂ）は、周波数変動時系列解析演算手段６３２及び周波数傾き時系列解析演算手段６３３により得られた被験者Ａのエアパック脈波の周波数変動時系列波形と周波数傾き時系列波形を示す。
周波数変動時系列解析演算手段６３２により得られた周波数変動時系列波形の基線位置は、飲酒前が１．３７〜１．４４Ｈｚ程度であるのに対し、飲酒後２０〜４０分間が１．５〜１．５５Ｈｚの範囲、飲酒後９０〜１１０分間が１．５Ｈｚ〜１．６Ｈｚの範囲、飲酒後１６０〜１８０分間が１．４〜１．５Ｈｚ程度であり、非飲酒時よりも飲酒時の方が基線位置が高くなっていることがわかる。従って、この基線位置の閾値を設定すれば、飲酒状態判定手段６５０により飲酒状態を判定することができる。
また、図１６の周波数変動時系列波形を見ると、飲酒前は周波数の値が試験時間の全域に亘り僅かに上下しながらほぼ水平に推移しているのに対し、飲酒後２０〜４０分間、飲酒後９０〜１１０分間、飲酒後１６０〜１８０分間では、試験開始から約１０００秒程度まで右下がりの傾向で推移し、１０００秒経過後、試験終了まで上昇傾向にある。この試験は、シート５００に２０分間着座する試験であるが、時間経過と共に被験者は疲労する。このとき、飲酒前であれば、疲労の進行に対して交感神経の代償作用が機能し、活性度を上げて対応しようとするため、周波数変動時系列波形は右下がりにはならない。しかし、飲酒により生体が休息モードになると、交感神経の代償作用が機能しないために心拍変動が低下気味になる。これにより、周波数変動時系列波形は右下がりに推移していくと考えられる。従って、周波数変動時系列波形の基線位置に加えて、周波数変動時系列波形が右下がりに推移するか否かにより、飲酒状態を判定すれば、判定精度がよりよくなる。なお、試験終了間際において周波数変動時系列波形が上昇しているのは、試験がもうすぐ終了するという刺激が入るためである。
図１７は、図１６（ｂ）に示した周波数傾き時系列波形の積分値を示す。具体的には、図１６（ｂ）の周波数傾き時系列波形を正の傾きと負の傾きで別々に分けて積分を行った。周波数変動の傾き時系列波形の積分値を飲酒前後で比較すると、呼気アルコール濃度の増加に伴い、負の成分は増加傾向にあるが、一方呼気アルコール濃度の低下が生じると飲酒前の状態と差の少ない状態に戻っていくのがわかる。また正の成分は、呼気アルコール濃度の増加に伴い減少していき、呼気アルコール濃度の減少が生じると増加していく。これは周波数変動の傾き時系列波形はエアパック脈波の周波数成分のゆらぎ度合を現しているためであり、飲酒によって、エアパック脈波の周波数変動のゆらぎの仕方が変化していることがわかる。従って、飲酒状態判定手段６５０は、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定することができる。
図１８に図１６（ｂ）の周波数変動の傾き時系列波形の周波数解析結果を示す。飲酒による呼気アルコール濃度の増加に伴い、卓越周波数ピークが低周波側に移っていくのがわかる。その後、呼気アルコール濃度の低下に伴い、高周波側に移ってゆく。これは図１７に示す積分値の変動の様子と一致し、飲酒によりエアパック脈波の周波数の変動のゆらぎからカオス性がなくなり、ゆらぎが単調になっていることを示唆する。
図１９はエアパック脈波のリアプノフ指数を用いて、時間窓１００秒、オーバーラップ時間を９０秒間に変えたエアパック脈波のパワースペクトルを示す。飲酒によるエアパック脈波のリアプノフ指数のゆらぎの強弱はパワースペクトルの強弱で表わされる。このパワースペクトルの低下からも、飲酒による呼気アルコール濃度の増加は、エアパック脈波の周波数変動のゆらぎを小さくすることが分かる。
図２０は被験者Ａの指尖容積脈波から得られた脈拍数変動のウェーブレット解析結果を示す。飲酒により、交感神経の影響を表すＬＦ／ＨＦのバースト波の出る割合が増加し、またＨＦの基線のレベルが低下している。飲酒により交感神経の亢進を示していると考えられる。飲酒後には、飲酒前の安静時に比べて、心臓副交感（迷走）神経活動が抑制され、心臓交換神経活動が亢進した状態にあることが推測される。さらに脈波のカオス解析から、交感神経緊張状態や心筋梗塞にてＩＡＢＰ（大動脈バルンパンピング）中でＩＡＢＰにより機械的に脈波が生成される場合などでは、脈波が単純化してカオス性が低下することが指摘されている。これらの知見より、飲酒による呼気アルコール濃度の増加に伴う指尖容積脈波・エアパック脈波の周波数のゆらぎを小さくしているのは、飲酒に伴う交感神経の亢進とそれに伴う脈波の単純化が関連していることが考えられる。また、時間の経過による呼気アルコール濃度の減少に伴い、この交感神経の亢進は抑えられていくが、飲酒前の状態と異なり、緊張を伴ったもので、リラックス状態には戻りきっていないことも確認された。尚、ほかの３名の試験事例も本試験例と同様の結果になった。
また、生体信号測定装置としては、上記したエアパック１０を用いたものに限らず、図２１に示したものを用いることもできる。図２１に示した生体信号測定装置２００は、三次元立体編物２１０、三次元立体編物支持部材２１５、フィルム２１６、板状発泡体２２１，２２２、振動センサ２３０を有して構成される。
三次元立体編物２１０は、図１等に示した生体信号測定装置１と同様のものを用いることができる。三次元立体編物２１０は、厚み方向の荷重−たわみ特性が、測定板上に載置して直径３０ｍｍ又は直径９８ｍｍの加圧板で加圧した際に、荷重１００Ｎまでの範囲で、人の臀部の筋肉の荷重−たわみ特性に近似したバネ定数を備えることが好ましい。具体的には直径３０ｍｍの加圧板で加圧した際の当該バネ定数が０．１〜５Ｎ／ｍｍの範囲、又は、直径９８ｍｍの加圧板で加圧した際の当該バネ定数が１〜１０Ｎ／ｍｍであるものを用いることが好ましい。人の臀部の筋肉の荷重−たわみ特性に近似していることにより、三次元立体編物と筋肉とが釣り合い、生体信号が伝播されると、三次元立体編物が人の筋肉と同様の振動を生じることになり、生体信号を大きく減衰させることなく伝播できる。
板状発泡体２２１，２２２は、ビーズ発泡体により構成することが好ましい。ビーズ発泡体としては、例えば、ポリスチレン、ポリプロピレン及びポリエチレンのいずれか少なくとも一つを含む樹脂のビーズ法による発泡成形体が用いることができる。ビーズ発泡体からなる板状発泡体２２１，２２２は、個々の微細なビーズを構成している発泡により形成された球状の樹脂膜の特性により、微小な振幅を伴う生体信号を膜振動（横波）として伝播する。この膜振動（横波）が三次元立体編物に弦振動として伝わり、これらの膜振動（横波）と弦振動が重畳され、生体信号は、膜振動（横波）と弦振動が重畳されることによって増幅された機械振動として、後述する振動センサ２３０により検出される。従って、生体信号の検出が容易になる。
板状発泡体２２１，２２２をビーズ発泡体から構成する場合、発泡倍率は２５〜５０倍の範囲で、厚さがビーズの平均直径以下に形成されていることが好ましい。例えば、３０倍発泡のビーズの平均直径が４〜６ｍｍ程度の場合では、板状発泡体２２１，２２２の厚さは３〜５ｍｍ程度にスライスカットする。これにより、板状発泡体２２１，２２２に柔らかな弾性が付与され、振幅の小さな振動に共振し、フィルム上を伝わる横波に減衰が生じにくくなる。なお、板状発泡体２２１，２２２は、本実施形態のように、三次元立体編物２１０を挟んで両側に配置されていても良いが、いずれか片側、好ましくは、シートバック側のみに配置した構成とすることもできる。
ここで、三次元立体編物２１０は、幅４０〜１００ｍｍ、長さ１００〜３００ｍｍの範囲の短冊状のものが用いられる。この大きさのものだと、三次元立体編物２１０に予備圧縮（連結糸に張力が発生する状態）を生じやすくなり、人と三次元立体編物２１０との間で平衡状態が作りやすい。本実施形態では、人が背部が当接した際の違和感軽減のため、脊柱に対応する部位を挟んで対象に２枚配設するようにしている。三次元立体編物２１０を簡単に所定位置に配置するようにするため、図２１に示したように、三次元立体編物２１０は三次元立体編物支持部材２１５に支持させた構成とすることが好ましい。三次元立体編物支持部材２１５は、板状に成形され、脊柱に対応する部位を挟んで対称位置に、縦長の配置用貫通孔２１５ａ，２１５ａが２つ形成されている。三次元立体編物支持部材２１５は、上記板状発泡体２２１，２２２と同様に、板状に形成されたビーズ発泡体から構成することが好ましい。三次元立体編物支持部材２１５をビーズ発泡体から構成する場合の好ましい発泡倍率、厚さの範囲は上記板状発泡体２２１，２２２と同様である。但し、生体信号により膜振動（横波）をより顕著に起こさせるためには、三次元立体編物２１０，２１０の上下に積層される板状発泡体２２１，２２２の厚さが、三次元立体編物支持部材２１５の厚さよりも薄いことが好ましい。
三次元立体編物支持部材２１５に形成した配置用貫通孔２１５ａ，２１５ａに、２つの三次元立体編物２１０，２１０を挿入配置した状態で、三次元立体編物２１０，２１０の表側及び裏側にフィルム２１６，２１６を積層する。なお、配置用貫通孔２１５ａ，２１５ａの形成位置（すなわち、三次元立体編物２１０，２１０の配設位置）は、心房と大動脈（特に、「下行大動脈」）の拍出に伴う動きによって生じる揺れ及び大動脈弁の動き（心部揺動波）を検知可能な領域に相当する位置とすることが好ましい。この結果、三次元立体編物２１０，２１０は、上下面が板状発泡体２２１，２２２によりサンドイッチされ、周縁部が三次元立体編物支持部材２１５によって取り囲まれており、板状発泡体２２１，２２２及び三次元立体編物支持部材２１５が共振箱（共鳴箱）の機能を果たす。
また、三次元立体編物支持部材２１５よりも、三次元立体編物２１０，２１０の方が厚いものを用いることが好ましい。つまり、三次元立体編物２１０，２１０を配置用貫通孔２１５ａ，２１５ａに配置した場合には、三次元立体編物２１０，２１０の表面及び裏面が、該配置用貫通孔２１５ａ，２１５ａよりも突出するような厚さ関係とする。これにより、フィルム２１６，２１６の周縁部を配置用貫通孔２１５ａ，２１５ａの周縁部に貼着すると、三次元立体編物２１０，２１０は厚み方向に押圧されるため、フィルム２１６，２１６の反力による張力が発生し、該フィルム２１６，２１６に固体振動（膜振動（横波））が生じやすくなる。一方、三次元立体編物２１０，２１０にも予備圧縮が生じ、三次元立体編物の厚さ形態を保持する連結糸にも反力による張力が生じて弦振動が生じやすくなる。なお、フィルム２１６，２１６は、三次元立体編物２１０，２１０の表側及び裏側の両側に設けることが好ましいが、いずれか少なくとも一方に設けた構成とすることも可能である。フィルム２１６，２１６としては、例えば、ポリウレタンエラストマーからなるプラスチックフィルム(例えば、シーダム株式会社製、品番「ＤＵＳ６０５−ＣＤＲ」）等を用いることができる。
振動センサ２３０は、上記したフィルム２１６，２１６を積層する前に、いずれか一方の三次元立体編物２１０に固着して配設される。三次元立体編物２１０は一対のグランド編地と連結糸とから構成されるが、各連結糸の弦振動がグランド編地との節点を介してフィルム２１６，２１６及び板状発泡体２２１，２２２に伝達されるため、振動センサ２３０は感知部２３０ａを三次元立体編物２１０の表面（グランド編地の表面）に固着することが好ましい。振動センサ２３０としては、マイクロフォンセンサ、中でも、コンデンサ型マイクロフォンセンサを用いることが好ましい。本実施形態では、マイクロフォンセンサを配置した部位（すなわち、三次元立体編物２１０を配置した配置用貫通孔２１５ａ）の密閉性を考慮する必要がないため、マイクロフォンセンサのリード線の配線は容易に行うことができる。生体信号によって生じる人の筋肉を介した体表面の振動は、三次元立体編物２１０だけでなく、板状発泡体２２１，２２２、フィルム２１６にも伝播され、それらが振動（弦振動、膜振動（横波））して減衰を防止しつつ重畳されて増幅する。よって、振動センサ２３０は、三次元立体編物２１０に限らず、振動伝達経路を構成する板状発泡体２２１，２２２及びフィルム２１６に、その感知部２３０ａを固定することもできる。
生体信号測定装置２００としては、図２１に示したものに限定されず、例えば、図２２（ａ）に示したように、２つの三次元立体編物２１０，２１０の両方を覆うことのできる大きさのフィルム２１７を少なくとも一方に用いても良い。また、図２２（ｂ）に示したように、略長方形の三次元立体編物を両端縁から中央部に向かって折り曲げ、重なり合った部分の中央部を縫い合わせたランバーサポート２１８を配置するようにしてもよい。ランバーサポート２１８は、三次元立体編物支持部材２１５に面ファスナ等を用いて固定される。ランバーサポート２１８をこのようにして設けることにより、狭いスペース中でストローク感を高めるのに寄与する。
上記した生体信号測定装置２００は、例えば、図２３に示したように、自動車用シート１０００のシートバックフレーム１１００に被覆される表皮１２００の内側に配置される。なお、配置作業を容易にするため、生体信号測定装置２００を構成する三次元立体編物２１０、三次元立体編物支持部材２１５、フィルム２１６、板状発泡体２２１，２２２、振動センサ２３０等は予めユニット化しておくことが好ましい。
上記した生体信号測定装置２００によれば、生体信号により、筋肉の荷重−たわみ特性に近似する荷重−たわみ特性をもつ板状発泡体２２１，２２２やフィルム２１６に膜振動（横波）が生じると共に、人の筋肉の荷重−たわみ特性に近似した荷重−たわみ特性を有する三次元立体編物２１０に弦振動が生じる。そして、三次元立体編物２１０の弦振動は再びフィルム２１６等の膜振動（横波）に影響を与え、これらの振動が重畳して作用する。その結果、生体信号に伴って体表面から入力される振動は、減衰することなく弦振動と膜振動（横波）との重畳によって増幅された固体振動として直接振動センサ２３０により検出されることになる。
図１等に示したエアパック１０内の空気圧変動を検出する生体信号測定装置１の場合、体積と圧力が反比例関係にあるため、密閉袋の体積を小さくしないと圧力変動を検出しにくい。これに対し、図２１〜図２２に示した生体信号測定装置２００によれば、空気圧変動ではなく、上記のように、機械的増幅デバイス（三次元立体編物２１０、板状発泡体２２１，２２２、フィルム２１６又はフィルム２１７）に伝播される増幅された固体振動を検出するものであるため、その容積（体積）が検出感度の観点から制限されることはほとんどなく、心部揺動波という振幅の小さな振動を感度良く検出できる。このため、多様な体格を有する人に対応できる。従って、図２１〜図２２に示した生体信号測定装置２００は、乗物用シートのように、多様な体格を有する人が利用し、さらに多様な外部振動が入力される環境下においても感度良く生体信号を検出できる。
（試験例２）
図２２に示した生体信号測定装置２００を、図２３に示したように自動車用シート１０００に装着し、試験例１と同様の飲酒実験を行った。
図２２に示した生体信号測定装置２００を、図２３に示したように自動車用シート１０００に装着し、試験例１と同様の飲酒実験を行った。
被験者は健康な２０歳代から４０歳代の日本人男性４名（Ｇ，Ｈ，Ｉ，Ｊ）としてそれぞれ上記シート１０００に着座させて心部揺動波の計測を行った。心部揺動波の計測と同時に、指尖容積脈波を光学式指尖容積脈波計（（株）アコム製、フィンガークリッププローブ ＳＲ−５Ｃ）により計測し、心部揺動波の計測の前後で呼気アルコール濃度を、「東海電子（株）製 ＡＬＣ−ｍｉｎｉ」により計測した。被験者には，あらかじめ飲酒実験とは別の日にエタノールパッチテストを行い、活性型（ＮＮ型）であることを確認した。尚、被験者の体重および身長は、被験者Ｇは体重７６ｋｇ、身長１７８ｃｍ、被験者Ｈは体重６４ｋｇ、身長１６７ｃｍ、被験者Ｉは体重５１ｋｇ、身長１７３ｃｍ、被験者Ｊは体重６１ｋｇ、身長１７４ｃｍである。その他の試験方法等は試験例１と全く同様である。
結果を図２４〜図２８に示す。なお、図２４は、心部揺動波の計測前後に測定した呼気アルコール濃度の変動の様子を示したものである。また、被験者Ｉ，Ｊについては、飲酒実験の日とは別の日に飲酒無しで、ほぼ同時刻の実験スケジュールにて比較検証実験を行った。図２５は、心部揺動波の原波形の周波数解析結果であり、図２６は、図２５の卓越周波数の時系列変化を示している。図２７は、検証として示した指尖容積脈波の原波形の周波数解析結果であり、図２８は、図２７の卓越周波数の時系列変化を示している。
まず、図２５の心部揺動波の原波形の周波数解析結果と図２７の指尖容積脈波の原波形の周波数解析結果とを比較すると、本試験例において使用した生体信号測定装置２００から得られた生体信号である心部揺動波が、指尖容積脈波とピーク位置がいずれの被験者もほぼ一致しており、指尖容積脈波と同等の周波数特性を有していることがわかる。
そして、図２６の卓越周波数の時系列変化より、被験者Ｇ，Ｈ，Ｉは飲酒により卓越周波数が上昇していることがわかる。これは飲酒による心拍数の上昇によって生じていると考えられる。図２４と比較した場合、呼気アルコール濃度は、飲酒後１回目の計測（１２００−２４００秒）の際にピークを迎えているが、卓越周波数はやや遅れて、飲酒後２回目の計測（５４００−６６００秒）の際にピークを迎えていることがわかる。
被験者Ｊは、飲酒した際には卓越周波数はほとんど変化が見られない。しかし、飲酒無しの場合には、時間経過と共に卓越周波数が低下している。このことから、被験者Ｆの場合も、飲酒無しとの比較において、飲酒した際には相対的な卓越周波数の上昇が見られることがわかる。この点は、被験者Ｈも同様である。つまり、飲酒した場合には、飲酒無しとの比較において、相対的に、卓越周波数が時間経過と共に上昇する現象が捉えられる。
なお、図２６の心部揺動波の卓越周波数の時系列変化は、図２８の指尖容積脈波の卓越周波数の時系列変化とほぼ同様の傾向を示しており、本試験例において使用した生体信号測定装置２００から得られた生体信号を利用することが有効であることがわかる。
図２９は、本試験例において使用した生体信号測定装置２００から得られた生体信号（心部揺動波）を、周波数変動時系列解析演算手段６３２によって処理し、その結果得られた周波数変動時系列波形を用いて呼気アルコール濃度との相関を捉えたグラフである。具体的には、着座疲労の影響の少ない０−６００秒間の周波数変動時系列波形の平均値を求め、その平均値と飲酒前との値との差分値を求め、得られた差分値を呼気アルコール濃度との相関で示したものである。なお、図２９（ａ）はピーク検出法を用いて周波数変動時系列波形を用いた場合の結果であり、（ｂ）はゼロクロス法を用いて周波数変動時系列波形を用いた場合の結果である。
図２９（ａ），（ｂ）から、ピーク検出法、ゼロクロス法共に、飲酒による呼気アルコール濃度の上昇に伴い、周波数変動の差分値がプラス方向に増加し、その後の呼気アルコール濃度の低下に伴い、元の周波数に近づいていることがわかる。一方、飲酒無しの場合には差分値がマイナス方向に増加していくことが分かる。このことから、飲酒無しの場合と飲酒している場合とでは大きな差がある。そこで、このような周波数変動の差分値を求め、所定以上の差分値が生じている場合には、飲酒状態と推定するようにすることが可能である。また、図２９（ａ），（ｂ）に示したような呼気アルコール濃度と差分値との対応テーブルを個人毎にコンピュータに記憶させておけば、周波数変動の差分値を算出することで呼気アルコール濃度を推定することもできる。
１ 生体信号測定装置
１０ エアパック
１１ 表側エアパック
１１１ 小空気袋
１１１ｂ センサ
１１２ 三次元立体編物
１２ 裏側エアパック
１２１ 大空気袋
１２２ 三次元立体編物
１５ 収容体
１００ エアパックユニット
２０ 第１のビーズ発泡樹脂弾性部材
３０ 第２のビーズ発泡樹脂弾性部材
４０，４５ 三次元立体編物
５００ シート
５１０ シートバック部
５１１ 表皮部材
５１２ クッション支持部材
５２０ シートクッション部
６０ 飲酒検知システム
６００ 飲酒解析推定部
６１０ 周波数動的情報処理手段
６２１ 卓越周波数時系列波形演算手段
６２２ 卓越周波数変動時系列解析演算手段
６２３ 卓越周波数傾き時系列解析演算手段
６３１ 周波数演算手段
６３２ 周波数変動時系列解析手段
６３３ 周波数傾き時系列解析手段
６５０ 飲酒判定手段
１０ エアパック
１１ 表側エアパック
１１１ 小空気袋
１１１ｂ センサ
１１２ 三次元立体編物
１２ 裏側エアパック
１２１ 大空気袋
１２２ 三次元立体編物
１５ 収容体
１００ エアパックユニット
２０ 第１のビーズ発泡樹脂弾性部材
３０ 第２のビーズ発泡樹脂弾性部材
４０，４５ 三次元立体編物
５００ シート
５１０ シートバック部
５１１ 表皮部材
５１２ クッション支持部材
５２０ シートクッション部
６０ 飲酒検知システム
６００ 飲酒解析推定部
６１０ 周波数動的情報処理手段
６２１ 卓越周波数時系列波形演算手段
６２２ 卓越周波数変動時系列解析演算手段
６２３ 卓越周波数傾き時系列解析演算手段
６３１ 周波数演算手段
６３２ 周波数変動時系列解析手段
６３３ 周波数傾き時系列解析手段
６５０ 飲酒判定手段
Claims (16)
- 生体信号測定装置により人の上体から採取した生体信号を解析して飲酒の有無を推定する飲酒解析推定部を備えた飲酒検知システムであって、
前記飲酒解析推定部が、
前記生体信号の時系列波形からその周波数に関する時系列変動を求める周波数動的情報処理手段と、
前記周波数動的情報処理手段により得られた周波数に関する時系列変動の傾向が非飲酒時における周波数に関する時系列変動の傾向と乖離している場合に飲酒状態と判定する飲酒状態判定手段と
を具備することを特徴とする飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、
前記飲酒状態判定手段は、前記卓越周波数時系列波形が非飲酒時よりも高いか否かにより飲酒状態を判定する構成である請求項１記載の飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、
前記飲酒状態判定手段は、所定時間範囲における前記卓越周波数時系列波形の卓越周波数の値のバラツキ度合いが、非飲酒時よりも拡大しているか否かにより飲酒状態を判定する構成である請求項１記載の飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段を有し、
前記飲酒状態判定手段は、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成である請求項２記載の飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段と、
前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記卓越周波数の傾きの時系列変化を卓越周波数傾き時系列波形として出力する卓越周波数傾き時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記卓越周波数傾き時系列解析演算手段により得られる卓越周波数傾き時系列波形の基線位置に対して、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも乖離しているか否かにより飲酒状態を判定する構成である請求項２記載の飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成である請求項１記載の飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより判定する請求項１記載の飲酒検知システム。 - 前記周波数動的情報処理手段は、
前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否か、及び、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定する構成である請求項１記載の飲酒検知システム。 - 生体信号測定装置により人の上体から採取した生体信号を分析して飲酒の有無を推定する飲酒検知システムの記憶部に設定される飲酒解析推定部を構成するコンピュータプログラムであって、
前記飲酒解析推定部が、
前記生体信号の時系列波形からその周波数に関する時系列変動を求める周波数動的情報処理手段と、
前記周波数動的情報処理手段により得られた周波数に関する時系列変動の傾向が非飲酒時における周波数に関する時系列変動の傾向と乖離している場合に飲酒状態と判定する飲酒状態判定手段と
を具備することを特徴とするコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、
前記飲酒状態判定手段は、前記卓越周波数時系列波形が非飲酒時よりも高いか否かにより飲酒状態を判定する構成である請求項９記載のコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記生体信号の所定時間範囲における時系列波形を周波数解析して卓越周波数を得て、卓越周波数時系列波形を求める卓越周波数時系列波形演算手段を有し、
前記飲酒状態判定手段は、所定時間範囲における前記卓越周波数時系列波形の卓越周波数の値のバラツキ度合いが、非飲酒時よりも拡大しているか否かにより飲酒状態を判定する構成である請求項９記載のコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段を有し、
前記飲酒状態判定手段は、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成である請求項１０記載のコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記卓越周波数の平均値の時系列変化を卓越周波数変動時系列波形として出力する卓越周波数変動時系列解析演算手段と、
前記卓越周波数時系列波形演算手段により得られる卓越周波数時系列波形を、所定のオーバーラップ時間で設定した所定の時間窓毎に前記卓越周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記卓越周波数の傾きの時系列変化を卓越周波数傾き時系列波形として出力する卓越周波数傾き時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記卓越周波数傾き時系列解析演算手段により得られる卓越周波数傾き時系列波形の基線位置に対して、前記卓越周波数変動時系列波形の基線位置が非飲酒時よりも乖離しているか否かにより飲酒状態を判定する構成である請求項１０記載のコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否かにより飲酒状態を判定する構成である請求項９記載のコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより判定する請求項９記載のコンピュータプログラム。 - 前記周波数動的情報処理手段は、
前記生体信号の時系列波形における周波数の時系列波形を求める周波数演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の平均値を求める移動計算を行い、時間窓毎に得られる前記周波数の平均値の時系列変化を周波数変動時系列波形として出力する周波数変動時系列解析演算手段と、
前記周波数演算手段により得られた前記生体信号の周波数の時系列波形において、所定のオーバーラップ時間で設定した所定の時間窓毎に前記周波数の傾きを求める移動計算を行い、時間窓毎に得られる前記周波数の傾きの時系列変化を周波数傾き時系列波形として出力する周波数傾き時系列解析演算手段と
を有し、
前記飲酒状態判定手段は、前記周波数変動時系列波形の基線位置が非飲酒時よりも高いか否か、及び、前記周波数傾き時系列波形を、正の傾きと負の傾きとに分け、それぞれの傾きを積分し、飲酒前の正常状態の積分値と比較して、正の傾きの積分値が予め設定した差分以上減少し、負の傾きの積分値が予め設定した差分以上増加しているか否かにより飲酒状態を判定する構成である請求項９記載のコンピュータプログラム。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2011514424A JP5669287B2 (ja) | 2009-05-19 | 2010-05-18 | 飲酒検知システム及びコンピュータプログラム |
Applications Claiming Priority (4)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2009120849 | 2009-05-19 | ||
| JP2009120849 | 2009-05-19 | ||
| PCT/JP2010/058363 WO2010134525A1 (ja) | 2009-05-19 | 2010-05-18 | 飲酒検知システム及びコンピュータプログラム |
| JP2011514424A JP5669287B2 (ja) | 2009-05-19 | 2010-05-18 | 飲酒検知システム及びコンピュータプログラム |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPWO2010134525A1 true JPWO2010134525A1 (ja) | 2012-11-12 |
| JP5669287B2 JP5669287B2 (ja) | 2015-02-12 |
Family
ID=43126202
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2011514424A Expired - Fee Related JP5669287B2 (ja) | 2009-05-19 | 2010-05-18 | 飲酒検知システム及びコンピュータプログラム |
Country Status (4)
| Country | Link |
|---|---|
| US (1) | US9149231B2 (ja) |
| EP (1) | EP2433559B1 (ja) |
| JP (1) | JP5669287B2 (ja) |
| WO (1) | WO2010134525A1 (ja) |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP5679541B2 (ja) * | 2010-01-26 | 2015-03-04 | 株式会社デルタツーリング | 生体信号検出装置 |
| US8698639B2 (en) | 2011-02-18 | 2014-04-15 | Honda Motor Co., Ltd. | System and method for responding to driver behavior |
| US9292471B2 (en) | 2011-02-18 | 2016-03-22 | Honda Motor Co., Ltd. | Coordinated vehicle response system and method for driver behavior |
| JP5977938B2 (ja) | 2011-11-25 | 2016-08-24 | 株式会社デルタツーリング | 生体状態推定装置及びコンピュータプログラム |
| US9192334B2 (en) | 2013-01-31 | 2015-11-24 | KHN Solutions, Inc. | Method and system for monitoring intoxication |
| US9788772B2 (en) * | 2013-01-31 | 2017-10-17 | KHN Solutions, Inc. | Wearable system and method for monitoring intoxication |
| US10153796B2 (en) | 2013-04-06 | 2018-12-11 | Honda Motor Co., Ltd. | System and method for capturing and decontaminating photoplethysmopgraphy (PPG) signals in a vehicle |
| US9751534B2 (en) | 2013-03-15 | 2017-09-05 | Honda Motor Co., Ltd. | System and method for responding to driver state |
| US9272689B2 (en) * | 2013-04-06 | 2016-03-01 | Honda Motor Co., Ltd. | System and method for biometric identification in a vehicle |
| US10537288B2 (en) | 2013-04-06 | 2020-01-21 | Honda Motor Co., Ltd. | System and method for biological signal processing with highly auto-correlated carrier sequences |
| US10213162B2 (en) | 2013-04-06 | 2019-02-26 | Honda Motor Co., Ltd. | System and method for capturing and decontaminating photoplethysmopgraphy (PPG) signals in a vehicle |
| US10499856B2 (en) | 2013-04-06 | 2019-12-10 | Honda Motor Co., Ltd. | System and method for biological signal processing with highly auto-correlated carrier sequences |
| CN103538477A (zh) * | 2013-11-01 | 2014-01-29 | 李金泽 | 一种酒驾自动断电装置 |
| US10542961B2 (en) | 2015-06-15 | 2020-01-28 | The Research Foundation For The State University Of New York | System and method for infrasonic cardiac monitoring |
| US20190102963A1 (en) * | 2017-09-29 | 2019-04-04 | Toyota Motor Engineering & Manufacturing North America, Inc. | Systems and methods for securing an object in a vehicle |
| US11006895B2 (en) | 2018-03-22 | 2021-05-18 | KHN Solutions, Inc. | Method and system for transdermal alcohol monitoring |
| US11602306B2 (en) | 2021-01-12 | 2023-03-14 | KHN Solutions, Inc. | Method and system for remote transdermal alcohol monitoring |
Family Cites Families (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4789342B2 (ja) | 2001-05-10 | 2011-10-12 | 株式会社デルタツーリング | クッション材、シート及びパイル糸の植毛方法 |
| JP2004121625A (ja) * | 2002-10-04 | 2004-04-22 | Seiko Instruments Inc | 脈波検出装置及びフーリエ変換処理装置 |
| JP4347621B2 (ja) | 2003-05-21 | 2009-10-21 | 株式会社デルタツーリング | 生体評価システム、コンピュータプログラム及び記録媒体 |
| JP4247055B2 (ja) | 2003-05-21 | 2009-04-02 | 株式会社デルタツーリング | 運転席用座席システム |
| CN100475143C (zh) | 2004-03-25 | 2009-04-08 | 株式会社三角工具加工 | 载荷体状态判定装置、车辆座椅 |
| JP4959178B2 (ja) | 2005-02-28 | 2012-06-20 | 株式会社デルタツーリング | クッション材及び圧力変動検出装置 |
| JP4663387B2 (ja) * | 2005-04-22 | 2011-04-06 | 株式会社日立製作所 | 生体光計測装置 |
| JP4751683B2 (ja) | 2005-09-28 | 2011-08-17 | Ｋｂセーレン株式会社 | 着色弾性繊維不織布 |
| JP5044230B2 (ja) | 2007-02-14 | 2012-10-10 | 株式会社デルタツーリング | 生体信号分析装置、シート及び生体信号分析方法 |
| JP5179798B2 (ja) * | 2007-05-22 | 2013-04-10 | 株式会社デルタツーリング | 座席構造 |
| JP2009219541A (ja) | 2008-03-13 | 2009-10-01 | Toyota Motor Corp | 飲酒状態検出装置 |
| WO2010021229A1 (ja) | 2008-08-19 | 2010-02-25 | 株式会社デルタツーリング | 生体信号測定装置及び生体状態分析システム |
-
2010
- 2010-05-18 EP EP10777758.3A patent/EP2433559B1/en not_active Not-in-force
- 2010-05-18 JP JP2011514424A patent/JP5669287B2/ja not_active Expired - Fee Related
- 2010-05-18 WO PCT/JP2010/058363 patent/WO2010134525A1/ja active Application Filing
- 2010-05-18 US US13/321,367 patent/US9149231B2/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| WO2010134525A1 (ja) | 2010-11-25 |
| EP2433559A1 (en) | 2012-03-28 |
| JP5669287B2 (ja) | 2015-02-12 |
| EP2433559B1 (en) | 2018-04-04 |
| US20120130261A1 (en) | 2012-05-24 |
| EP2433559A4 (en) | 2014-12-31 |
| US9149231B2 (en) | 2015-10-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP5669287B2 (ja) | 飲酒検知システム及びコンピュータプログラム | |
| JP5429946B2 (ja) | 生体状態分析装置及びコンピュータプログラム | |
| US9622708B2 (en) | Biological body state estimation device and computer program | |
| JP5044230B2 (ja) | 生体信号分析装置、シート及び生体信号分析方法 | |
| JP5561552B2 (ja) | 生体信号測定装置及び生体状態分析システム | |
| JP5553303B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
| JP5582478B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
| JP4959178B2 (ja) | クッション材及び圧力変動検出装置 | |
| JP5892678B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
| JP5800357B2 (ja) | 生体信号検出機構 | |
| JP5327584B2 (ja) | 生体状態分析装置、コンピュータプログラム及び記録媒体 | |
| JP5751475B2 (ja) | 生体状態推定装置及びコンピュータプログラム | |
| JP5236424B2 (ja) | 脈波検出装置及び生体状態分析装置 | |
| JP5679556B2 (ja) | 飲酒状態検知装置、コンピュータプログラム及び記録媒体 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130517 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20140704 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140901 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20141114 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20141213 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 5669287 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |