JPS5975038A

JPS5975038A - 心電計の波形認識におけるｐ波の検出方法

Info

Publication number: JPS5975038A
Application number: JP57185787A
Authority: JP
Inventors: 修栗原; 石原　博海
Original assignee: HAGIWARA DENKI KK; Suzuken KK
Current assignee: HAGIWARA DENKI KK; Suzuken KK
Priority date: 1982-10-22
Filing date: 1982-10-22
Publication date: 1984-04-27

Abstract

(57)【要約】本公報は電子出願前の出願データであるた
め要約のデータは記録されません。

Description

【発明の詳細な説明】本発明は、心臓病の検査に使用される心電図自動解析装
置のための、心電計の波形認識におけるＰ波の検出方法
に関するものである。

心電図自動解析装置においては、房室ブロック１度・ｗ
ｐｗ及び房拡大等の病名診断のためにＰ波を用いていた
。

従って、Ｐ波を正確に求める事により、房室ブロック１
度・ｗｐｗ及び房拡大の診断を正確に行うことができる
。

従来の心電図自動解析装置について、第１図〜第３図を
参照し、以下説明する。

心電図自動解析装置は、制御するためのプログラムを記
憶するＲＯＭ２と、用途に応じて変更しく仔るデータを
記１意するためのＲＡＭ３と、プログラムメモリとして
の前記ＲＯＭ２の内容に従って順次制御および演算を行
なうＣ’Ｐ　Ｕ　１と、プログラムに応じて一定時間毎
に割込み信号を前記ｃｐＵ１に対して発する割込みタイ
マ４と、前記ＣＰＵ１と外部回路とのインターフェイス
を図る入出力回路５とを有し、前記ＣＰＵＩ・ＲＯＭ２
・ＲＡＭ３・割込みタイマ４および入出力回路５とがそ
れぞれアドレスバス・データバスおよびコントロールハ
ス等から成るハス６により接続されている。また、入出
力回路５には、心電計１０からのアナログ信号をデジタ
ル信号に変換するためのＡ−Ｄ変換器７と、メソセージ
やデータを表示するためのＣＲＴ８と、コマンドやデー
タを入力するためのキーボード９とが接続されている。

第１図の心電図自動解析装置の動作について、第２図の
システムプログラムを参照して以下説明する。

電源を投入することによりＲＯＭ２内のプログラムは自
動的にスタート１０１を行ないステップ１０２に進む。

ステップ１０２では、ＣＲＴ８に質問事項（例えば、検
診者番号：より・年齢・日付及び性別等）を表示し、オ
ペレータが前記質問事項に対しキーホード９からデータ
を入力すると、ステップ１０３に進む。ステップ１０３
では、心電計の計測時間をセットし、キーボード９のス
タートスイッチを操作することにより、心電計１０と割
込みタイマ４が作動し、ステップ１０４に進む。ステッ
プ１０４では、心電計１０の複数の出力が前記割込みタ
イマ４に同期してサンプルボールドされ、その後順次Ａ
／Ｄ変換した値がＲＡＭ３にデータＤＡとして読込まれ
、ステップ１０５に進む。ステップ１０５では、全誘導
のデータ取込みが終了すれば次のステップ１０６に進む
。また、全誘導のデータ取込のが終了してなければステ
ップ１０３に進み、全誘導のデータを取込むまでステッ
プ１（１３，１０４および１０５をループする。

ステップ１０６では、−誘導のデータ＝ＤＡ毎に移動平
均法および分割平均法、荷重平均法等を用い、雑音の除
去を行なったデータ：ＤＢを得、そのデータを前記ＲＡ
Ｍ３に記憶し、ステップ１０７に進む。ステップ１０７
では、前記データ：１）　Ｂを一定レベルにてスライス
して第３図に示したＱＲ３波のおおよその位置を１４認
した後、このＱＲ３波の位置を基準にして予め定められ
た補正可能な範囲を定めておき直線近似で基線の補整を
行なったデータ：　ｒ　（ｘｉを記憶した後、ステップ
１０８へ進む。基線の補正を行なったデータ：ｆ（ｘ）
は、アドレスをＸとして第３図に示すようなグラフとし
て表わすことができる。すなわち、アドレスＸのメモリ
にデータｆ　（ｘｉを集録した集まりとなる。ステップ
１０８では、前記ＱＲ３波の位置から各パターンの大体
の区分を検出した後、詳細なパターンの検出を行ない、
ステップ１０９へ進む。

ステップ１０９では、ステップ１０８にて検出した波形
パターンの区分を利用し、これらの波形について、波高
値および波形幅、波形の形、時間間隔などを計測し、ス
テップ１１０に進む。ステップｌｌＯでは、全誘導の計
測が終了していればステップ１１１に進み、また、全誘
導の計測が終了するまでステップ１０６〜１１０をルー
プし゛、全誘導の計測が終了するとステップ１１１に進
む。

ステップ１１１では、前記計」り結果を診断論理に用い
て診断結果を得、ステップ１１２に進む。ステップ１１
２では、前記診断結果をＣＲＴ８およびプリンタ１１に
出力し、ステップ１１３に進みプログラムは終了する。

従来Ｐ波の検出方法としては、コロナ社より昭和５１年
５月１５Ｅ］７版（新版）発行の「新版電子計算機ハン
ドブック」のページ４８４〜４８６に記載の「心電図自
動診断システム」の項で詳細に説明されているような、 ■）心電計のアナログ出力を比較器を使用してコンバレ
ー１−（スライス）してＰ帯を得、そのＰ帯から波形パ
ターンの区分を正しく識別する方法、２）心電計の出力
をＡ−Ｄ変換し、そのＡ−Ｄ変換により得られたサンプ
リングデータの基線を補正しノイズ除去した後、データ
ｆ　（Ｘ）をコンパレートしてＰ帯を得、そのＰ帯から
波形パターンの区分を正しく識別する方法等があった。

しかし、従来の心電計の波形をコンパレート（スライス
）してＰ帯を検出する方法に於ては、心電波の基線がド
リフ１−や電極の取付等の影響を受は変化し、誤ったＰ
帯を検出し正確なＰ波の頂点を検出することができなか
った。また、基線補正及びノイズ除去を行った後のデー
タｆ　（ＸｌをコンパレートとしてＰ帯を得る方法に於
ては、小さなＰ波の場合にはノイズとみなし検出するこ
とができなかった。

本発明はかかる欠点を無くすために、心電計より取り込
んだデータを雑音の除去及び基線の補正をしてデータｆ
（×）を得、該データｆ　ｆＸｌ中に含まれる原点Ｎ及
びＴ波の終点Ｍを予め検出し、その後前記原点Ｎから検
索始点Ｉと前記Ｔ波の終点Ｍから検索終点Ｊとをそれぞ
れ定め、さらに必要に応じて前記原点Ｎから第２の検索
終点Ｋを定め、前記検索始点■と検索終点ＪまたばＫと
の間を所定間隔にて検索点ｉを順次変え、これら各検索
点ｉを中心とする前記データｆ　（Ｘ）のサンプリング
数２ｙ（ＩＩの和からその前後の前記データｆ（×）の
サンプリング数７個の和をそれぞれ減算して得られる演
算結果Ｆ　（Ｘ）を前記各検索点ｉの移行に伴って順次
得て、これら各演算結果Ｆ　（Ｘ）と前記演算結果の最
大値に応じて定まる所望の値とを比較してＰ帯を順次検
出し、ついで各Ｐ帯における前記各演算結果Ｆ　（Ｘ）
の最大値を検出し、この最大値に対応した検索点のデー
タｆ　（ｘｉを基準点として設定し、その後前記基準点
と前記各データｆ（×）との比較に基づきそれぞれの変
曲点を求め、これら変曲点の中心をＰ波の頂点として検
出するようにした心電計の波形認識におけるＰ波の検出
方法を提供することにある。

本発明によれば、Ｐ波の検出方法にあたり前記な変化を
も拡大して検出でき、小さなＰ波の検出も可能となり、
また、ノイズの多いデータにおいても多点の積算を行う
ことによりノイズの影響を少なくでき、正確なＰ波の検
出を行うことができる。

以下、本発明の実施例について第３図〜第６図を使用し
説明する。

商用周波数の２ｎ倍（但し、ｎ＝１・２・３・・・）の
サンプリング周波数：Ｓとして２００１１ｚ（５ｍｓｅ
ｃ）を用い、心電計の出力をＡ−Ｄ変換して得られたデ
ータは、雑音の除去及び基線の補正をし得られるデータ
ｆ（に）に変換される。そのデータｆ（×）の一部が第
３図に示されている。

パターン認識のプログラムにより、原点（ＱＲ８Ｒ８軸
点始点びＴ波の終点を、それぞれＮ及びＭのデータ番！
　（ＲＡＭのアドレス）として検出した後、第５図のＰ
帯検出プログラムを実行する。

Ｐ帯検出プログラムは、ステップ２０１で開始しステッ
プ２０２に進む。ステップ２０２では、サンプリングデ
ータｆ　（Ｘｉの積算時間Ｔ１を設定し、そのデータの
サンプリンタ冒固数ｙをザンプリングデータ周波数Ｓと
積算時間Ｔ１との積として求めると共に、最初の積算点
の中心とする検索始点Ｉを原点アドレスＮから前記デー
タのサンプリング個数ｙの２倍を引算して得、その後ス
テップ２０３に進む。

積算時間Ｔ１は、積算点を多くしノイズの影響を少なく
するため、２０〜３Ｑｍｓｅｃが適当である。本実施例
においては、前記積算時間Ｔ１が２５ｍ５ｅｃまた前記
データのサンプリング個数ｙが５として以下説明する。

ステップ２０３では、原点ＮからＰ波の最大確認範囲に
応じた検索終点Ｊを前記検索始点Ｉから最大確認時間Ｔ
２として、Ｊ＝Ｉ−Ｔ２ｘＳの関係から求め、ステップ
２０４に進む。但し、前記最大確認時間Ｔ２は良く知ら
れた時間で原点Ｎから４００　ｍ、　ｓｅｃ前後である
。この場合には、Ｔ２−３３５ｍ５ｅｃとして以下説明
する。従って、前記検索終点Ｊは、■−６７である。ス
テップ２０４では、Ｔ波の終点ＭからＴ３時間たったま
でをＰ波確認領域とし、その検索終点Ｋをに＝Ｍ＋Ｔ３
ＸＳとして求め、ステップ２０５に進む。この場合には
、Ｔ３＝７５ｍｓｅｃとして、前記検索終点ＫがＭ＋１
５である。ステップ２０５では、積算値の演算結果の最
大値を格納するためのアドレス又はレジスフをＢとして
、前記アドレスＢに０をセントすると共に、検索点ｉに
検索始点Ｉをセントし、さらに前記演算結果を格納する
ＲＡＭのアドレスＸを前記検索点ｉからＡだけオフセッ
トされ検索点ｉに対応したアドレスＡ＋ｉにセットし、
ステップ２０６に進む。

ステップ２０６では、Ｐ波の検索最終中心点ＫまたはＪ
よりも検索点ｉが大きいか等しければステップ２０７に
進む。ステップ２０７では、検索点ｉを中心としてＴ１
時間の２倍の各データｆ（×）の積算和からその両側の
それぞれＴ１時間の各データｆ　（Ｘｌの積算和を減算
して得られるデータＦ（Ｘ）を得るため、２０８に進む・。ステ・ンブ２０８では、アドレスＸ（
すなわちＡ＋　ｉ）に前記各演算結果Ｆ　（Ｘ）を収納
し、ステップ２０９に進む。ステップ２０９では、前記
演算結果Ｆ　（Ｘ）の最大値を検出するため、アドレス
ＢのデータＦ　（Ｂ）よりも演算結果Ｆ　（Ｘ）の絶対
値の方が大きげればステップ２１０に進み、アドレスＢ
に前記演算結果Ｆ　（Ｘ）の絶対値を収納し、ステップ
２１１に進む。また、ステップ２０９では、前記データ
Ｆ　（Ｂ）よりも前記結果Ｆ　［Ｘｌの絶対値の方が小
さければステップ２１１に進む。ステップ２１１では、
検索点のアドレスｉ及び演算結果のアドレスを更新する
ために、１＝ｉ−１を実行した後、ステップ２０６にも
どる。ステップ２０６では、まだＰ波の検索終点Ｋまた
はＪよりも検索点ｉが大きいが等しげれば、前述の様に
ステップ２０７，２０８，２０９゜（２１０）、２１１
へ進み、前記検索点ｉが前記検索終点ＪまたばＫよりも
小さくなるまでループし、小さくなるとステップ２１２
に進む。その演算結果Ｆ　（Ｘ）がデータｆ　ｆｘ）と
対応して、第４図（ｂ）に示しである。

ステップ２１２では、検索点ｉに検索始点Ｉ　（但し、
Ｉ＝Ｎ−２ｙ）をセットし、その後前記演算結果のデー
タＦ　（Ｘ）の最初のアドレスＸにセントするため、Ｘ
＝Ａ＋ｉを実行すると共に、次のデータエントリの最初
のアドレスを前記検索点ｉよりＣだけオフセットされ検
索点に対応したアドレスＣ＋＋にセットした後、ステッ
プ２１３に進む。ステップ２１３では、前記演算結果Ｆ
　（Ｘ）の最大値であるＦ　（Ｂ）からＰ帯の位置を検
出するため、Ｆ　（Ｂ）　Ｘｋ＝Ｆ　（Ｂ’）　　（但
し、ｋは０゜７５〜０．９５本実施例では０．８５を使
用している。）を実行し前記演算結果Ｆ　（Ｘ）のスラ
イス点Ｆ（Ｂ’）をもとめ、ステ、プ２１４に進む。ス
テップ２１４では、前記Ｐ波の検索終点ＫまたはＪより
も検索点ｉが大きければ、ステップ２１５へ進む。

ステップ２１５では、前記結果Ｆ　（Ｘ）の絶対値がス
ライス点Ｆ（Ｂ’）よりも大きければ、ステップ２１６
へ進みアドレスＣにデータａを収納しステップ２１８に
進む。また、前記演算結果Ｆ　（Ｘ）の絶対値がスライ
ス点Ｆ（Ｂ’）よりも大きくなければ、ステップ２１７
へ進みアドレスＣにデータ０を収納しステップ２１８に
進む。ステップ２１８では、１＝ｉ−１を実行すること
によりアドレスＣを更新し、ステップ２１４にもどる。

ステップ２１８では、まだＰ波の検索終点ＫまたはＪよ
りも検索点ｉが大きいか等しければ、前述の様にステッ
プ２１５，２１．６．または’２１７，２１８へ進み、
前記検索点ｉが検索終点ＪまたはＫよりも小さくなるま
でループし、その後ステップ２１９に進み、本プログラ
ムは終了する。

次に、第６図（ａ）に示されているＰ被検出プログラム
に進む。

Ｐ被検出プログラムは、ステップ３０１によりスタート
し、ステップ３０２に進む。ステンプ３０２では、本プ
ログラムの初期化を行うべく、Ｐ化生の最大値を格納す
るためのアドレスＤの内容Ｆ　（Ｄ）をＯにセットし、
検索点ｉに検索始点Ｉ（但し、ＩはＮ−２ｙである。）
をセントし、さらに前記検索点に対応したアドレスにす
るためアドレスＣをＣ＝Ｃ＋　ｉにセットすると共にＰ
波の頂点を格納するアドレスＰをセットし、ステップ３
０３に進む。ステップ３０３では、Ｐ波の検索終点Ｋま
たはＪよりも検索点ｉが大きいか等しければ、ステップ
３０４に進む。ステップ３０４では、比較結果を格納し
たアドレスＣの内容Ｆ　（Ｃ）がａであれは、Ｐ波検出
ルーチン４００　（後で詳細に説明する。）に進み実行
した後、ステップ３０５に進む。また、ステップ３０４
では、前記アドレスＣの内容Ｆ　（Ｃ）がａでなければ
、直ちにステップ３０５に進む。ステップ３０５では、
アドレスを更新するためｊ＝ｉ−１を実行した後ステッ
プ３０３にもどる。前記ステップ３０３では、まだ前記
検索終点Ｊまたは■（よりも検索点ｉが大きいか等しけ
れば、前述の様にステップ３０４゜（４００）、３０５
−１進み、前記検索点ｉが前記検索終点ＪまたはＫより
も大きいか等しくなるまでループし、その後ステップ３
０６に進み、本プログラムは終了する。

第６図（ｂ）に示されている、Ｐ波検出ルーチンについ
て説明する。

Ｐ波検出ルーチンは、ステップ３０４で比較結果を格納
したアドレスＣの内容Ｆ　（Ｃ）がａであれば起動され
、ステップ４０１によりスタートし、ステップ４０２に
進む。ステップ４０２では、前記演算結果Ｆ　（Ｘ）の
最大値の検索範囲ρを設定するため、ｌ１＝Ｔ４×Ｓ　
（但し、Ｔ４は検索時間で４０〜８０ｍ　ｓｅｃである
）のＴ　４　＝　６０ｍ　ｓｅｃとして実行し、！−１
２を得た後、ステップ４０３に進む。ステップ４０３ば
、Ｐ帯における前記演算結、ＩＦ　（Ｘ）の最大値を検
出するため、アドレスＤのデータＦ　（Ｄ）と前記各演
算結果ｌＦくＸ）１とを比較し、Ｆ　（Ｄ）　＜　ＩＦ
　（Ｘ）　ｌ　　であれば、ステップ４０４へ進み、マ
ｔ、−Ｆ（Ｄ）　＜　ｌ　Ｆ　（Ｘ）１でなければ、ス
テップ４０５に進む。ステ・７プ４０４では、アドレス
ＤにアドレスＸの内容Ｆ　（Ｘ）を収容するとともに、
アドレスＤ＋］にアドレスＸの１の部分だＬＪを収納し
、ステップ４０５に進む。ステップ４０５では、アドレ
スＸを更新するために１＝ｉ−１を実行するとともに、
前記検索範囲βを更新するためにρ−β−１を実行し、
ステップ４０６に進む。ステップ４０６では、前記検索
範囲βがＯでなければ、ステップ４０３にもどり、前記
検索範囲ｐがＯになるまでステップ４０３、　　（４０
４）、４０５をループし、前記検索範囲βがＯになると
ステップ４０７に進む。ステップ４０７ては、アドレス
ＸにアドレスＤ＋１の内容ｆ（Ｄ−１４）すなわちｉを
収納するとともに、アドレスＥにアドレスＸの内容ＮＸ
Ｉ（すなわちＰ帯における最大値）を収納し、ステップ
４０８に進む。ステップ４０８では、アドレスＸを更新
するためにｉ　＝　ｉ　＋　ｌを実効し、ステップ４０
９に進む。ステップ４０９では、原点ｎに近い側のＰ波
の変曲点を求むべく、アドレスＥのデータｆ　（Ｅ）と
アドレスＸのデータｆ　（Ｘ）とを比較し、前記Ｐ帯に
おける演算結果Ｆ　（Ｘ）の最大値Ｆ　（Ｄ）が正の値
の時には（正のＰ波の時には）、ｆ（Ｅ）＜　ｆ　（ｘ
ｉであればステップ４０８にもどり順次アドレスＸを更
新しｆ　　（Ｅ）≧ｆ　（ｘ）になるまで捜しｆ（Ｅ）
≧ｆ　（Ｘ）になるとステップ４１０に進む。また、前
記Ｐ帯における演算結果Ｆ　（Ｘ）の最大値Ｆ　（Ｄ）
が負の値の時にはく負のＰ波の時には）、ｆ　　（Ｅ）
≧ｆ　ｆＸ）であればステップ４０Ｂにもどり順次アド
レスＸを更新しｆ　　（Ｅ）　＜　ｆ　（ｘｉになるま
で捜しｆ　　（Ｅ）　＜　ｆ　（ｘ）になるとステップ
４１０に進む。ステップ４１０では、変曲点のアドレス
ＸをデータとしてアドレスＤ＋２に収納するとともに、
次の変曲点すなわちＴ波に近い側のＰ波の変曲点すなわ
ちＴ波に近い側のＰ波の変曲点を求むべく、演算結果の
最大値のアドレスのｉの部分のデータｆ（×）をアドレ
スＸにセ・ノドし、ステップ４１１に進む。ステップ４
１１では、アドレスを更新するために１＝ｉ−１を実行
し、ステップ４１２に進む。ステップ４１２では、アド
レスＥのデータｆ　　（Ｅ）とアドレスＸのデータｆ　
（ｘ）とを比較し、前記Ｐ帯における演算結果Ｆ、（Ｘ
）の最大値Ｆ　（Ｄ）が正の値の時にはく正のＰ波の時
には）、ｆ　　（Ｅ）　＜　ｆ　（ｘ）であればステッ
プ４１１にもどり順次アドレスを更新し、ｆ　　（Ｅ）
≧ｆ　（ｘｉになるとステップ４１３に進む。また、前
記Ｐ帯仝こおける演算結果Ｆ　（Ｘ）の最大値Ｆ　（Ｄ
）が負の値の時には（負のＰ波の時には）、ｆ（Ｅ）≧
ｆ（×）であればステップ４１１にもどり順次アドレス
Ｘを更新しｆ　　（Ｅ）　＜　ｒ（ｘ＋になるとステッ
プ４１３に進む。ステップ４１３では、変曲点のアドレ
スＸをデータとしてアドレスＤ＋３に収納し、ステップ
４１４に進む。ステップ４１４では、真のＰ波の原点を
求めるために、Ｆ　（Ｐ）　−（Ｆ　（Ｄ＋２）＋Ｆ（
Ｄ＋３））／２の演算を実行するとともに、アドレスＰ
に前記演算結果Ｆ　（Ｐ）を収納し、次の演算結果の収
納アドレスを更新すべく　Ｐ＝Ｐ＋１を実行し、さらに
次のＰ帯の最大値を検出すべくアドレスＤのデータＦ　
（Ｄ）を０にセントし、ステップ４１５に進む。ステッ
プ４１５では、真のＰ波の頂点から時間Ｔ５（但し、Ｔ
５は３０〜６０　ｍ５ｅｃ）が４５　ｍ５ｅｃのところ
から再度次のＰ波を検索するために、アドレスｉをセッ
トすべくｉ＝Ｆ　（Ｐ）−Ｔ５ＸＳを実行し、ステップ
４１６へ進みＰ波検出ルーチンを終了する。

本発明は、Ｐ波のおおよその位置を（１）式を利用して
Ｐ帯として求め、前記Ｐ帯に対応した原波形ｆ　（ｘｉ
より真のＰ波の原点を求める方法を用いることにより、（１）基線のドリフトや電極の取り付けにより生じるノ
イズの影響を無くす事ができ、（２）小さなＰ帯をも検出することができる。

従って、Ｐ波を本方法を使用して正確に検出することが
でき、診断理論を用い不整脈および房拡大等を正確に見
付ける事ができる。

【図面の簡単な説明】

第１図は、心電図自動解析装置のブロック図である。第
２図は、システムプログラムを示す図である。第３図は
、基線補正を行なったデータｆ　（ｘｉをアドレスＸに
より表わした図である。第４図（ａ）は、第３図のＰ波
近傍を拡大して示した図である。第４図（ｂｌは、Ｐ帯検出プログラム中のステップ２０
７の演算結果を示す図である。第４図（Ｃ）は、Ｐ帯検
出プログラムにより得られたＰ帯を示す図である。第５
図は、Ｐ帯検出プログラムを示すフローチャートである
。第６図は、Ｐ波の最大点を検出するプログラム示すフ
ローチャートである。符号の説明　　　１０・・・心電図、ｆ　（Ｘ）・・・データ、Ｎ・・・原
点、Ｍ・・・終点、■・・・検索始点、Ｊ、Ｋ・・・検
索終点、ｉ・・・検索点、Ｆ（Ｘ）　　・・・演算結果
。出願人　　萩原電気株式会社（ばか１名）代理人　　弁
理士　長　谷　照　− （ばか１名）

Claims

【特許請求の範囲】

（１）心電計より取り込んだデータを雑音の除去及び基
線の補正をしてデータｆ　（Ｘｌを得、該データｆ（×
）中に含まれる原点Ｎ及びＴ波の終点Ｍを予め検出し、
その後前記原点Ｎから検索始点■と前記Ｔ波の終点Ｍか
ら検索終点Ｊとをそれぞれ定め、さらに必要に応じて前
記原点Ｎから第２の検索終点Ｋを定め、前記検索始点■
と検索終点ＪまたばＫとの間を所定間隔にて検索点ｉを
順次変え、これら各検索点ｉを中心とする前記データｆ
　（ｘｌのサンプリング数２ｙ（Ｉｌｉｌの和からその
前後の前記データｆ　（Ｘ！のサンプリング数７個の和
をそれぞれ減算して得られる演算結果Ｆ　（Ｘ）を前記
各検索点ｉの移行に伴って順次得て、これら各演算結果
Ｆ（Ｘ）と前記演算結果の最大値に応じて定まる所望の
値とを比較してＰ帯を順次検出し、ついで各Ｐ帯におけ
る前記各演算結果Ｆ　（Ｘ）の最大値を検出し、この最
大値に対応した検索点ｉの前記データｆ　（Ｘ）を基準
点として設定し、その後前記基準点と前記各データｆ　
（ｘ＋との比較に基づきそれぞれの変曲点を求め、これ
ら変曲点の中心をＰ波の頂点として検出するようにした
心電計の波形認識におけるＰ波の検出方法。
（２）前記演算結果Ｆ　（Ｘ）が下記式によって得られ
るようにした特許請求の範囲第１項記載の心