JPH0654815A - Method and device for rr interval spectral analysis - Google Patents
Method and device for rr interval spectral analysisInfo
- Publication number
- JPH0654815A JPH0654815A JP4232772A JP23277292A JPH0654815A JP H0654815 A JPH0654815 A JP H0654815A JP 4232772 A JP4232772 A JP 4232772A JP 23277292 A JP23277292 A JP 23277292A JP H0654815 A JPH0654815 A JP H0654815A
- Authority
- JP
- Japan
- Prior art keywords
- interval
- spectrum analysis
- spectrum
- analysis
- range
- 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.)
- Pending
Links
- 238000010183 spectrum analysis Methods 0.000 title claims abstract description 96
- 238000000034 method Methods 0.000 title claims description 24
- 238000004458 analytical method Methods 0.000 claims abstract description 30
- 238000001228 spectrum Methods 0.000 claims description 26
- 230000007774 longterm Effects 0.000 claims description 10
- 238000001514 detection method Methods 0.000 abstract description 3
- 230000002567 autonomic effect Effects 0.000 abstract 1
- 206010061666 Autonomic neuropathy Diseases 0.000 description 9
- 210000000467 autonomic pathway Anatomy 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000008035 nerve activity Effects 0.000 description 6
- 210000005037 parasympathetic nerve Anatomy 0.000 description 4
- 238000007781 pre-processing Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 3
- 230000002889 sympathetic effect Effects 0.000 description 3
- 208000009729 Ventricular Premature Complexes Diseases 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 2
- 206010003119 arrhythmia Diseases 0.000 description 2
- 230000006793 arrhythmia Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000000747 cardiac effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 238000012217 deletion Methods 0.000 description 1
- 230000037430 deletion Effects 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 230000001771 impaired effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 210000005036 nerve Anatomy 0.000 description 1
- 238000012805 post-processing Methods 0.000 description 1
- 238000003825 pressing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
【０００１】[0001]
【産業上の利用分野】本発明はＲＲ間隔スペクトル分析
方法及びその装置、特に長時間心電図解析装置から入力
したＲＲ間隔の変動をスペクトル分析するＲＲ間隔スペ
クトル分析方法及びその装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to an RR interval spectrum analysis method and apparatus, and more particularly to an RR interval spectrum analysis method and apparatus for spectrally analyzing fluctuations in the RR interval input from a long-term electrocardiogram analysis apparatus.
【０００２】[0002]
【従来の技術】一般に、心電図を構成するＲ波とＲ波と
の間隔を、ＲＲ間隔と称することは、よく知られてい
る。近年、心機能障害と自律神経障害との関連性から、
このＲＲ間隔をスペクトル分析する方法、即ち、ＲＲ間
隔変動のパワースペクトルを高速フーリエ変換を用いて
周波数分析する方法が、自律神経に障害があるか否かを
発見する一手法として、注目を浴びて来ている。Ｒ−Ｒ
間隔の変動は、心電図のＲ波のみのトレンド波形である
Ｒ−Ｒタコグラムに表現される。そして、Ｒ−Ｒタコグ
ラムに現れたＲＲ間隔の変動は、自律神経が正常の場合
は、一定のゆらぎがあり、一方、自律神経に障害のある
場合は、このゆらぎがないか、又は非常に小さい、とい
うことが知られている。このような意味を持つＲ−Ｒタ
コグラムにおいて、スペクトル分析を行う範囲を指定
し、この範囲内でスペクトル分析を行うと、自律神経が
正常な場合と異常な場合とでは、分析結果の波形上で、
異なる結果が表れる。即ち、自律神経が正常な場合は、
一般に、自律神経を構成する交感神経と副交感神経を示
す２つのピークが現れる。しかし、自律神経が異常な場
合には、この２つのピークが現れないか、又は非常に小
さい。このようにして、ＲＲ間隔の変動をスペクトル分
析した結果の波形上において、現れた２つのピークの状
態を見ることにより、自律神経障害を発見することが可
能である。従来は、この自律神経障害を発見する一手法
であるＲＲ間隔スペクトル分析方法を、次のように行っ
ていた。先ず、ＲＲ間隔スペクトル分析装置の画面上
に、上述したＲＲ間隔の変動を表すＲ−Ｒタコグラムが
表示される。同時に、上記Ｒ−Ｒタコグラムの波形上
に、２本のカーソルが表示される。次に、この２本のカ
ーソルを、Ｒ−Ｒタコグラムの波形上において、同時に
移動させることにより、スペクトル分析対象範囲を指定
する。このようにして、２本のカーソルで挟まれた領域
が、スペクトル分析対象範囲となる。最後に、この指定
された範囲内で、ＲＲ間隔スペクトル分析を行い、得ら
れた結果を見ながら、上述した２つのピークの有無やパ
ワー値の大きさ等により、自律神経障害を発見すること
ができる。2. Description of the Related Art Generally, it is well known that an interval between R waves forming an electrocardiogram is called an RR interval. In recent years, from the relationship between cardiac dysfunction and autonomic neuropathy,
The method of spectrally analyzing the RR interval, that is, the method of frequency-analyzing the power spectrum of the RR interval fluctuation by using the fast Fourier transform has attracted attention as one method for discovering whether or not there is a disorder in the autonomic nerve. It is coming. RR
The fluctuation of the interval is expressed in the RR tachogram which is a trend waveform of only the R wave of the electrocardiogram. The fluctuation of the RR interval appearing in the RR tachogram has a certain fluctuation when the autonomic nerve is normal, while there is no fluctuation or a very small fluctuation when the autonomic nerve is impaired. It is known that. In the RR tachogram having such a meaning, when a range for spectral analysis is specified and spectral analysis is performed within this range, the autonomic nerve is normal and abnormal. ,
Different results appear. That is, when the autonomic nerve is normal,
In general, two peaks showing sympathetic nerves and parasympathetic nerves forming the autonomic nerve appear. However, when the autonomic nerve is abnormal, these two peaks do not appear or are very small. In this way, it is possible to detect the autonomic neuropathy by observing the states of the two peaks that appear on the waveform obtained as a result of the spectrum analysis of the fluctuation of the RR interval. Conventionally, the RR interval spectrum analysis method, which is one method for discovering this autonomic neuropathy, has been performed as follows. First, on the screen of the RR interval spectrum analyzer, the RR tachogram representing the variation of the RR interval described above is displayed. At the same time, two cursors are displayed on the waveform of the RR tachogram. Next, the two cursors are simultaneously moved on the waveform of the RR tachogram to specify the spectrum analysis target range. In this way, the region sandwiched by the two cursors becomes the spectrum analysis target range. Finally, within this designated range, RR interval spectrum analysis is performed, and while looking at the obtained results, autonomic neuropathy can be detected by the presence or absence of the above-mentioned two peaks and the magnitude of the power value. it can.
【０００３】[0003]
【発明が解決しようとする課題】しかしながら、以上の
従来例においては、以下のような課題があった。即ち、
ＲＲ間隔スペクトル分析方法において使用されるＲＲ間
隔の変動データは、通常は、長時間心電図解析装置が取
り扱う２４時間の長時間心電図に基づいて解析されたも
のが、多い。従って、ＲＲ間隔の変動状態を表すＲ−Ｒ
タコグラムも、長時間のものとなる。しかし、従来のよ
うに、２本のカーソルを、Ｒ−Ｒタコグラムの波形上に
おいて、データの頭から順次移動させて行くことによ
り、スペクトル分析対象範囲を指定することは、Ｒ−Ｒ
タコグラムが長時間になればなる程、時間がかかること
になる。スペクトル分析対象範囲を指定するまでに時間
がかかれば、それだけ、ＲＲ間隔スペクトル分析自体も
遅れることになり、自律神経障害も早期に発見できない
場合がある。本発明の目的は、スペクトル分析対象範囲
の指定を迅速に行うことにより、ＲＲ間隔スペクトル分
析を早く行い、以て自律神経障害の早期発見に資するこ
とにある。However, the above-mentioned conventional examples have the following problems. That is,
Most of the variation data of the RR interval used in the RR interval spectrum analysis method is usually analyzed on the basis of a 24-hour long electrocardiogram handled by a long-term electrocardiogram analyzer. Therefore, RR that represents the fluctuation state of the RR interval
The tachogram will also be long. However, as in the prior art, it is possible to specify the spectrum analysis target range by sequentially moving two cursors from the head of the data on the waveform of the RR tachogram.
The longer the tachogram, the longer it will take. If it takes time to specify the spectrum analysis target range, the RR interval spectrum analysis itself will be delayed accordingly, and autonomic neuropathy may not be detected early. An object of the present invention is to quickly specify the range to be subjected to spectrum analysis, thereby performing RR interval spectrum analysis quickly, thereby contributing to early detection of autonomic neuropathy.
【０００４】[0004]
【課題を解決するための手段】上記課題は、第１ステッ
プＱ１において、ＲＲ間隔データを入力し、第２ステッ
プＱ２において、該入力したＲＲ間隔データのうちでス
ペクトル分析を行うべき範囲を、数値を入力することに
より、指定し、第３ステップＱ３において、該指定した
範囲内でスペクトル分析を行い、第４ステップＱ４にお
いて、その分析結果を出力するようにしたことを特徴と
するＲＲ間隔スペクトル分析方法である第１発明（図
１）、及びＲＲ間隔データを入力する入力部１Ａと、該
入力したＲＲ間隔データのうちでスペクトル分析を行う
べき範囲を数値を入力することにより指定する指定部１
Ｅと、該指定部１Ｅにより指定した範囲内でスペクトル
分析を行なうスペクトル分析部１Ｃと、該スペクトル分
析部１Ｃにより分析された結果を出力する出力部１Ｄと
から成ることを特徴とするＲＲ間隔スペクトル分析装置
である第２発明（図２）により、解決される。The above problems are solved by inputting RR interval data in the first step Q1, and by setting a numerical range of the range in which the spectrum analysis is to be performed in the input RR interval data in the second step Q2. Specified by inputting, the spectrum analysis within the designated range is performed in the third step Q3, and the analysis result is output in the fourth step Q4. 1st invention (FIG. 1) which is a method, and an input unit 1A for inputting RR interval data, and a specifying unit 1 for specifying a range for spectral analysis in the input RR interval data by inputting numerical values.
E, an spectrum analysis unit 1C that performs spectrum analysis within a range designated by the designation unit 1E, and an output unit 1D that outputs a result analyzed by the spectrum analysis unit 1C. This is solved by the second invention (FIG. 2) which is an analyzer.
【０００５】[0005]
【作用】上記のとおり、本発明によれば、上記第１発明
（図１）と第２発明（図２）とにより、入力したＲＲ間
隔データのうちで、スペクトル分析を行うべき範囲を、
数値を入力することにより、指定することができによう
なった。上記構成によれば、従来のように、２本のカー
ソルを移動させることにより、スペクトル分析対象範囲
を指定していたのに比較すれば、指定がより一層迅速に
行えることになった。このため、ＲＲ間隔スペクトル分
析自体も早く行うことができ、その分、自律神経障害も
早期に発見されることとなった。As described above, according to the present invention, according to the first invention (FIG. 1) and the second invention (FIG. 2), the range in which the spectrum analysis is to be performed, among the input RR interval data,
You can specify it by entering a numerical value. According to the above configuration, by moving two cursors as in the conventional case, the spectrum analysis target range is designated, but in comparison, the designation can be performed more quickly. Therefore, the RR interval spectrum analysis itself can be performed earlier, and autonomic neuropathy can be detected earlier accordingly.
【０００６】[0006]
【実施例】以下、本発明を実施例により添付図面を参照
して説明する。本発明は、上述したように、第１発明
（図１）と第２発明（図２）とから構成されている。 Ａ．第１発明について 第１発明は、請求項１に記載したように、第１ステップ
Ｑ１において、ＲＲ間隔データを入力し、第２ステップ
Ｑ２において、該入力したＲＲ間隔データのうちでスペ
クトル分析を行うべき範囲を、数値を入力することによ
り、指定し、第３ステップＱ３において、該指定した範
囲内でスペクトル分析を行い、第４ステップＱ４におい
て、その分析結果を出力するようにしたことを特徴とす
るＲＲ間隔スペクトル分析方法である。以下、第１発明
に係る方法を、その各手順に従って、詳述する。先ず、
第１発明を実施するための装置であるＲＲ間隔スペクト
ル分析装置１（図２）の電源を投入すると、図３に示す
ように、「ＲＲ間隔スペクトル分析」という文字が表示
された画面Ｇ１が現れる。次いで、上記画面Ｇ１が消え
て、画面Ｇ２が表れる。この画面Ｇ２は、メニュー画面
といわれ、第１発明は、この画面Ｇ２の中の各項目Ｇ２
１、Ｇ２２、Ｇ２３、Ｇ２４をカーソル（図示省略）で
選択することにより、以下のように、実施される。 （１）ＲＲ間隔データを入力する（図１の第１ステップ
Ｑ１）。 図１の第１ステップＱ１において、ＲＲ間隔データを入
力するが、その前処理として、装置どうしの接続、後処
理として、データファイルの作成が行われる。 装置どうしの接続 第１発明に係る方法は、例えば、後述する第２発明であ
るＲＲ間隔スペクトル分析装置（図２）とこれに接続さ
れる長時間心電図解析装置２により、実施されるので、
両装置どうしを接続する。 ＲＲ間隔データの入力 上記画面Ｇ２おいて、項目Ｇ２１の「データ受信」を選
択すると、長時間心電図解析装置２（図２）からＲＲ間
隔データが入力される。このＲＲ間隔データは、既述し
たように、心電図を構成するＲ波とＲ波との間隔である
ＲＲ間隔の変動を表すデータであって、長時間心電図解
析装置２（図２）が扱う長時間心電図に基づいて解析さ
れたものである。従って、この入力段階において、どの
時間の範囲のＲＲ間隔データを入力するかを、指定す
る。その後は、この指定された範囲のＲＲ間隔データが
入力される。これら一連の操作は、全て、キー入力によ
り、行われる。尚、画面Ｇ２において、項目Ｇ２２の
「スペクトル分析」は、実際にスペクトル分析を行う場
合に、また、項目Ｇ２３の「環境設定」は、ＲＲ間隔ス
ペクトル分析装置１の種類に応じて、システム、データ
ドライブ等の環境設定を行う場合に、更に、項目Ｇ２４
の「終了」は、動作を終了させる場合に、それぞれ選択
するものである。 データファイルの作成 上記で入力されたＲＲ間隔データのデータファイルを
作成する。以降は、このデータファイルを使用して、動
作が行われる。DESCRIPTION OF THE PREFERRED EMBODIMENTS The present invention will be described below with reference to the accompanying drawings by way of embodiments. The present invention is composed of the first invention (FIG. 1) and the second invention (FIG. 2) as described above. A. About 1st invention As for 1st invention, as described in Claim 1, in 1st step Q1, RR interval data is input, and in 2nd step Q2, the spectrum analysis is performed among this input RR interval data. The power range is specified by inputting a numerical value, the spectrum analysis is performed within the specified range in the third step Q3, and the analysis result is output in the fourth step Q4. RR interval spectrum analysis method. Hereinafter, the method according to the first invention will be described in detail according to each procedure thereof. First,
When the power of the RR interval spectrum analyzer 1 (FIG. 2), which is the device for carrying out the first invention, is turned on, a screen G1 displaying the characters “RR interval spectrum analysis” appears as shown in FIG. . Next, the screen G1 disappears and the screen G2 appears. This screen G2 is called a menu screen, and the first invention is that each item G2 in this screen G2 is
By selecting 1, G22, G23, and G24 with a cursor (not shown), the operation is performed as follows. (1) Input RR interval data (first step Q1 in FIG. 1). In the first step Q1 of FIG. 1, the RR interval data is input, and as a preprocessing thereof, the devices are connected to each other, and as a post processing, a data file is created. Connection between Devices Since the method according to the first invention is carried out, for example, by the RR interval spectrum analyzer (FIG. 2) and the long-term electrocardiogram analyzer 2 connected thereto, which is the second invention described later,
Connect both devices. Input of RR Interval Data When "Data reception" in item G21 is selected on the screen G2, the RR interval data is input from the long-term electrocardiogram analyzer 2 (Fig. 2). As described above, the RR interval data is data representing fluctuations in the RR interval, which is an interval between R waves forming an electrocardiogram, and is a length that the long-term electrocardiogram analyzer 2 (FIG. 2) handles. It was analyzed based on the time electrocardiogram. Therefore, in this input step, the range of time to which the RR interval data is input is specified. After that, the RR interval data in the designated range is input. All of these series of operations are performed by key input. In the screen G2, the item "G22""spectralanalysis" indicates that the actual spectrum analysis is to be performed, and the item "G23""environmentsetting" indicates the system and data depending on the type of the RR interval spectrum analyzer 1. If you want to set the environment for the drive, etc.
“End” is selected when ending the operation. Creation of data file Create a data file of the RR interval data input above. After that, the operation is performed using this data file.
【０００７】（２）上記入力したＲＲ間隔データのうち
でスペクトル分析を行うべき範囲を、数値を入力するこ
とにより、指定する（図１の第２ステップＱ２）。次
に、図１の第２ステップＱ２において、スペクトル分析
対象範囲の指定が行われるが、これは、既述したよう
に、ＲＲ間隔スペクトル分析の前提となるものである。
この場合、図３に示す画面Ｇ２において、項目Ｇ２２の
「スペクトル分析」を選択する。 Ｒ−Ｒタコグラムの表示 項目Ｇ２２の「スペクトル分析」を選択すると、図４の
画面Ｇ３に切り替わり、ＲＲ間隔の変動を表すＲ−Ｒタ
コグラムＧ３１が画面Ｇ３の上段に表示される。 スペクトル分析対象範囲の指定 上記画面Ｇ３に表示されたＲ−ＲタコグラムＧ３１を見
ながら、スペクトル分析を行うべき範囲を、数値を入力
することにより、指定する。具体的な動作は、次のとお
りである。 ー１ 即ち、画面Ｇ３において、設定キー（図示省略）を押す
ことにより、図５の画面Ｇ４に切り替わる。この画面Ｇ
４は、後述するスペクトル分析を行う場合の種々の条件
を設定するための画面である。分析条件の設定は、図５
の画面Ｇ４を構成する各項目Ｇ４１、Ｇ４２、Ｇ４３、
Ｇ４４、Ｇ４５、Ｇ４６の内容を選択することにより、
行う。このうち、第２ステップＱ２のスペクトル分析対
象範囲の指定に最も関連がある項目は、項目Ｇ４２の分
析サンプル数である。この分析サンプル数とは、スペク
トル分析の対象となるＲＲ間隔データの長さ、即ち、Ｒ
波の数をいい、６４、１２８、２５６、５１２、１０２
４ビートの５種類のうちからいずれか１つを選択するよ
うにしたものである。その他、項目Ｇ４１の分析アルゴ
リズムとは、スペクトル分析をする場合のアルゴリズム
をＦＥＴ、ＭＥＭの２種類に定め、いずれか１つを選択
するようにしたものである。項目Ｇ４３の不整脈処理方
法とは、上記（２）の第２ステップＱ２において指定し
たスペクトル分析対象範囲Ｌ（図４）内に、不整脈の一
種である心室性期外収縮が存在すると、スペクトル分析
が正確にできないため、その心室性期外収縮の取扱を削
除、補間、無処理の３種類に定め、いずれか１つを選択
するようにしたものである。項目Ｇ４４のデータ前処理
方法とは、スペクトル分析では等間隔データ列を扱うの
で、ＲＲ間隔の不等間隔データ列から等間隔データ列を
作る場合の作り方をスプライン補間、無処理の２種類に
定め、いずれか１つを選択するようにしたものである。
項目Ｇ４５のＤＣ成分の除去とは、ＲＲ間隔のゆらぎの
基準となるＤＣ成分（平均ＲＲ間隔）のパワー値がゆら
ぎ成分のパワー値に比較して極めて大きくなるので、こ
のＤＣ成分の除去の仕方を除去、無処理の２種類に定
め、いずれか１つを選択するようにしたものである。項
目Ｇ４６の窓関数とは、本発明の対象であるＲＲ間隔デ
ータは有限区間データであるために、一般に無限関数を
取り扱うフーリエ変換に整合させるように、切り出した
有限区間データを繰り返して等価無限関数を作っている
が、有限区間の境界部の歪みによる高調波成分の発生を
回避する目的で等価無限関数全体に掛ける重み関数であ
る窓関数を方形、ハミング、Ｈａｎｎの３種類に定め、
いずれか１つを選択するようにしたものである。 ー２ 次に、上記項目４２の分析サンプル数を含む分析条件を
設定後、リターンキー（図示省略）を押すと、再び、図
４の画面Ｇ３に戻るので、スペクトル分析の開始時刻Ｔ
１をキーボード等から入力する。スペクトル分析の対象
となるＲＲ間隔データの長さは、上記ー１で述べた分
析サンプル数を（図５）、例えば、５１２と設定するこ
とにより、固定されているので、終了時刻をＴ２とすれ
ば、図４に示すように、時刻Ｔ１とＴ２で挟まれた範囲
Ｌが、スペクトル分析対象範囲となる。 Ｒ−Ｒタコグラムの拡大表示 上記のとおり、開始時刻を表す数値「Ｔ１」を入力する
ことにより、スペクトル分析対象範囲Ｌを指定すると、
この範囲Ｌを拡大したＲ−ＲタコグラムＧ３２が、画面
Ｇ３の下段に表示される。(2) Of the input RR interval data, the range for spectral analysis is designated by inputting a numerical value (second step Q2 in FIG. 1). Next, in the second step Q2 of FIG. 1, the spectrum analysis target range is specified, which is a prerequisite for the RR interval spectrum analysis as described above.
In this case, on the screen G2 shown in FIG. 3, the item "G22" of the spectrum analysis is selected. Display of RR tachogram When the item "Spectrum analysis" of the item G22 is selected, the screen is switched to the screen G3 in FIG. 4, and the RR tachogram G31 showing the variation of the RR interval is displayed in the upper part of the screen G3. Designation of Spectrum Analysis Target Range While observing the RR tachogram G31 displayed on the screen G3, the range for spectrum analysis is designated by inputting a numerical value. The specific operation is as follows. -1, that is, by pressing a setting key (not shown) on the screen G3, the screen is switched to the screen G4 in FIG. This screen G
Reference numeral 4 is a screen for setting various conditions when performing a spectrum analysis described later. The analysis conditions are set as shown in Fig. 5.
Items G41, G42, G43, which form the screen G4 of
By selecting the contents of G44, G45, G46,
To do. Of these, the item most relevant to the designation of the spectrum analysis target range in the second step Q2 is the number of analysis samples of item G42. The number of analysis samples is the length of RR interval data to be subjected to spectrum analysis, that is, R
64, 128, 256, 512, 102
One of five types of four beats is selected. In addition, the analysis algorithm of the item G41 is one in which an algorithm for spectrum analysis is set to two types, FET and MEM, and either one is selected. The arrhythmia processing method of the item G43 is that when a ventricular extrasystole, which is a kind of arrhythmia, exists in the spectrum analysis target range L (FIG. 4) designated in the second step Q2 of (2) above, the spectrum analysis is performed. Since it cannot be accurately performed, the handling of the ventricular extrasystole is defined as three types of deletion, interpolation, and no processing, and any one of them is selected. With the data preprocessing method of item G44, since the equidistant data string is handled in the spectrum analysis, the method of making the equidistant data string from the non-equidistant data string of the RR interval is set to two types: spline interpolation and no processing. , Any one of them is selected.
The removal of the DC component of item G45 means that the power value of the DC component (average RR interval), which is the reference for fluctuations in the RR interval, becomes extremely large compared to the power value of the fluctuation component. Are removed and unprocessed, and one of them is selected. The window function of the item G46 is an equivalent infinite function obtained by repeatedly cutting out the finite interval data so that the RR interval data, which is the object of the present invention, is finite interval data so as to match the Fourier transform that generally handles an infinite function. The window function, which is a weighting function to be applied to the entire equivalent infinite function for the purpose of avoiding the generation of the harmonic component due to the distortion of the boundary part of the finite section, is defined as three types of square, Hamming, and Hann.
One of them is selected. -2 Next, when the return key (not shown) is pressed after setting the analysis conditions including the number of analysis samples of the item 42, the screen G3 of FIG. 4 is displayed again.
Input 1 from the keyboard. The length of the RR interval data that is the target of the spectrum analysis is fixed by setting the number of analysis samples described in -1 above (FIG. 5), for example, 512, so the end time is set to T2. For example, as shown in FIG. 4, the range L sandwiched between times T1 and T2 is the spectrum analysis range. Enlarged display of RR tachogram As described above, when the spectrum analysis target range L is specified by inputting the numerical value "T1" representing the start time,
An RR tachogram G32 obtained by enlarging the range L is displayed in the lower part of the screen G3.
【０００８】（３）上記指定した範囲内でスペクトル分
析を行う（図１の第３ステップＱ３）。 図１の第３ステップＱ３において、上記指定した範囲Ｌ
内でスペクトル分析を行う。このスペクトル分析は、既
述したように、ＲＲ間隔変動のパワースペクトルを高速
フーリエ変換を用いて周波数分析することにより、行わ
れる。スペクトル分析中は、図５の画面Ｇ４において選
択したものが、画面上に表示されると共に（図示省
略）、「分析中」の文字が表示され、スペクトル分析が
終了すると、「分析終了」の文字が表示される（図示省
略）。(3) Spectral analysis is performed within the specified range (third step Q3 in FIG. 1). In the third step Q3 of FIG. 1, the specified range L
Perform spectral analysis within. As described above, this spectrum analysis is performed by frequency-analyzing the power spectrum of the RR interval fluctuation using the fast Fourier transform. During the spectrum analysis, the one selected on the screen G4 of FIG. 5 is displayed on the screen (not shown), and the character "in analysis" is displayed. When the spectrum analysis is completed, the character "analysis completed" is displayed. Is displayed (not shown).
【０００９】（４）上記スペクトル分析の分析結果を出
力する（図１の第４ステップＱ４）。 図１の第４ステップＱ４において、その分析結果を出力
する。 この分析結果は、例えば、図６の画面Ｇ５のように、
表示される。画面Ｇ５においては、図１の第２ステップ
Ｑ２において指定したスペクトル分析対象範囲Ｌ内にお
ける（図４）拡大Ｒ−ＲタコグラムＧ５１と、該スペク
トル分析対象範囲ＬのＲ−ＲヒストグラムＧ５２と、該
スペクトル分析対象範囲Ｌ内におけるスペクトル分析の
結果としてのパワースペクトルＧ５３とが、それぞれ表
示されている。上記表示されたパワースペクトルＧ５３
においては、２つのピークＧ５３１とＧ５３２が表れて
いる。これは、自律神経活動を構成する交感神経活動と
副交感神経活動を反映し、ピークＧ５３１が交感神経活
動及び副交感神経活動を示す成分であり、ピークＧ５３
２が副交感神経活動を示す成分である。上記拡大Ｒ−Ｒ
タコグラムＧ５１の縦軸はＲＲ間隔であるが、横軸は、
図５の画面Ｇ４における項目Ｇ４１２のデータ前処理方
法で、いずれかを選択したかにより時間又はビートとな
る。上記Ｒ−ＲヒストグラムＧ５２の縦軸は頻度であ
り、横軸はＲＲ間隔である。上記パワースペクトルＧ５
３の縦軸はパワー値であるが、横軸は、図５の画面Ｇ４
における項目Ｇ４１２のデータ前処理方法で、いずれか
を選択したかによりＨＺ又はｃ／ｂ若しくはＨＺｅｑで
ある。 また、この分析結果である画面Ｇ５のハードコピー
は、例えば、プリンタにより印刷することも可能であ
る。以上、（１）、（２）、（３）、（４）の各手順を
有する第１発明の動作は、図３の画面Ｇ２の項目Ｇ２４
「終了」を選択することにより、終了する。(4) The analysis result of the above spectrum analysis is output (fourth step Q4 in FIG. 1). In the fourth step Q4 of FIG. 1, the analysis result is output. The analysis result is, for example, as shown on the screen G5 of FIG.
Is displayed. In the screen G5, the expanded RR tachogram G51 within the spectrum analysis target range L specified in the second step Q2 of FIG. 1 (FIG. 4), the RR histogram G52 of the spectrum analysis target range L, and the spectrum A power spectrum G53 as a result of spectrum analysis in the analysis target range L is displayed. Power spectrum G53 displayed above
In, two peaks G531 and G532 are shown. This reflects the sympathetic nerve activity and the parasympathetic nerve activity that constitute the autonomic nerve activity, and the peak G531 is a component showing the sympathetic nerve activity and the parasympathetic nerve activity, and the peak G53.
2 is a component showing parasympathetic nerve activity. Expanded R-R
The vertical axis of the tachogram G51 is the RR interval, but the horizontal axis is
Depending on which one is selected in the data preprocessing method of the item G412 on the screen G4 of FIG. 5, it becomes time or beat. The vertical axis of the RR histogram G52 is frequency, and the horizontal axis is RR interval. The power spectrum G5
The vertical axis of 3 is the power value, but the horizontal axis is the screen G4 in FIG.
It is HZ or c / b or HZeq depending on which is selected in the data preprocessing method of the item G412 in. Further, the hard copy of the screen G5, which is the analysis result, can be printed by a printer, for example. As described above, the operation of the first invention having the procedures of (1), (2), (3), and (4) is performed by the item G24 on the screen G2 of FIG.
Select "End" to end.
【００１０】Ｂ．第２発明について 第２発明は、請求項２に記載したように、ＲＲ間隔デー
タを入力する入力部１Ａと、該入力したＲＲ間隔データ
のうちでスペクトル分析を行うべき範囲を数値を入力す
ることにより指定する指定部１Ｅと、該指定部１Ｅによ
り指定した範囲内でスペクトル分析を行なうスペクトル
分析部１Ｃと、該スペクトル分析部１Ｃにより分析され
た結果を出力する出力部１Ｄとから成ることを特徴とす
るＲＲ間隔スペクトル分析装置である。図２は、第２発
明の実施例を示す図である。同図において、参照符号１
はＲＲ間隔スペクトル分析装置、２は長時間心電図解析
装置である。上記ＲＲ間隔スペクトル分析装置１は、Ｒ
Ｒ間隔データを入力してスペクトル分析し、その結果を
出力する装置であり、上記長時間心電図解析装置２は、
長時間心電図を入力してそれを解析し、解析結果である
ＲＲ間隔データ等を出力する装置である。上記ＲＲ間隔
スペクトル分析装置１は、入力部１Ａと、指定部１Ｅ
と、スペクトル分析部１Ｃと、出力部１Ｄとから構成さ
れている。上記入力部１Ａは、ＲＲ間隔データを入力す
る装置である。上記指定部１Ｅは、ＲＲ間隔データのう
ちでスペクトル分析を行うべき範囲を数値を入力するこ
とにより指定する装置である。上記スペクトル分析部１
Ｃは、指定部１Ｅにより指定した範囲内でスペクトル分
析を行なう装置である。上記出力部１Ｄは、スペクトル
分析部１Ｃにより分析された結果を出力する装置であ
る。その他、制御部１Ｂは、スペクトル分析の制御を行
う等ＲＲ間隔スペクトル分析装置１全体の制御を行う装
置である。B. About 2nd invention As for 2nd invention, as described in Claim 2, the input part 1A which inputs RR interval data, and the range which should perform spectrum analysis among this input RR interval data are input numerically. It is characterized by comprising a designating section 1E designated by, a spectrum analyzing section 1C for performing spectrum analysis within a range designated by the designating section 1E, and an output section 1D for outputting a result analyzed by the spectrum analyzing section 1C. And an RR interval spectrum analyzer. FIG. 2 is a diagram showing an embodiment of the second invention. In the figure, reference numeral 1
Is an RR interval spectrum analyzer, and 2 is a long-term electrocardiogram analyzer. The RR interval spectrum analyzer 1 has R
A device for inputting R-interval data, performing spectrum analysis, and outputting the result. The long-term electrocardiogram analysis device 2 is
It is a device for inputting an electrocardiogram for a long time, analyzing it, and outputting RR interval data and the like, which are analysis results. The RR interval spectrum analyzer 1 includes an input unit 1A and a designation unit 1E.
And a spectrum analysis unit 1C and an output unit 1D. The input unit 1A is a device for inputting RR interval data. The designation unit 1E is a device that designates a range of the RR interval data for which spectrum analysis is to be performed by inputting a numerical value. The spectrum analysis unit 1
C is a device for performing spectrum analysis within the range designated by the designation unit 1E. The output unit 1D is a device that outputs the result analyzed by the spectrum analysis unit 1C. In addition, the control unit 1B is a device that controls the entire RR-interval spectrum analyzer 1, such as controlling spectrum analysis.
【００１１】以下、上記構成を有する第２発明の動作を
説明する。先ず、長時間心電図解析装置２からのＲＲ間
隔データ信号Ｓ１が、入力部１Ａに入力する。このＲＲ
間隔データ信号Ｓ１は、長時間のＲＲ間隔データを含む
ので、指定部１Ｅからの指定信号Ｓ５により、所定の時
間を指定すると、制御部１Ｂを介して、所定時間内のＲ
Ｒ間隔データ信号Ｓ２が、入力部１Ａから制御部１Ｂ
に、入力する。制御部１Ｂにおいては、この所定時間内
のＲＲ間隔データ信号Ｓ２に基づいて、その内部記憶装
置又は外部記憶装置（図示省略）に、データファイルが
作成される。指定部１Ｅからの指定信号Ｓ５により、種
々の分析条件を設定すると共に（図５）、スペクトル分
析を行うべき範囲Ｌを数値Ｔ１（図４）を入力すること
により指定すると、制御部１Ｂは、分析対象範囲のＲＲ
間隔データ信号Ｓ３をスペクトル分析部１Ｃに送信す
る。スペクトル分析部１Ｃにおいては、この分析対象範
囲のＲＲ間隔データ信号Ｓ３に基づいて、スペクトル分
析が行われ、該スペクトル分析部１Ｃは、分析結果であ
るパワースペクトルＧ５３（図６）等を含む分析結果信
号Ｓ４を出力する。上記分析結果信号Ｓ４を受信した制
御部１Ｂは、所定の処理をした後の信号Ｓ６を出力部１
Ｄに送信する。出力部１Ｄにおいては、入力信号Ｓ６に
基づいて、拡大Ｒ−ＲタコグラムＧ５１、Ｒ−Ｒヒスト
グラムＧ５２、パワースペクトルＧ５３（図６）がその
表示部１Ｄ１により画面Ｇ５上に表示され、又はそのプ
リンタ１Ｄ２により印刷される。The operation of the second invention having the above structure will be described below. First, the RR interval data signal S1 from the long-term electrocardiogram analyzer 2 is input to the input unit 1A. This RR
Since the interval data signal S1 includes long-time RR interval data, when a predetermined time is designated by the designation signal S5 from the designation unit 1E, the R within the predetermined time is set via the control unit 1B.
The R interval data signal S2 is transmitted from the input section 1A to the control section 1B.
To enter. In the control unit 1B, a data file is created in the internal storage device or an external storage device (not shown) based on the RR interval data signal S2 within the predetermined time. When various analysis conditions are set by the designation signal S5 from the designation unit 1E (FIG. 5) and the range L to be subjected to spectrum analysis is designated by inputting a numerical value T1 (FIG. 4), the control unit 1B RR of analysis target range
The interval data signal S3 is transmitted to the spectrum analyzer 1C. In the spectrum analysis unit 1C, spectrum analysis is performed based on the RR interval data signal S3 in the analysis target range, and the spectrum analysis unit 1C analyzes the analysis result including the power spectrum G53 (FIG. 6) as the analysis result. The signal S4 is output. Upon receiving the analysis result signal S4, the control unit 1B outputs the signal S6 after the predetermined processing to the output unit 1B.
Send to D. In the output section 1D, the enlarged RR tachogram G51, the RR histogram G52, and the power spectrum G53 (FIG. 6) are displayed on the screen G5 by the display section 1D1 or the printer 1D2 thereof based on the input signal S6. Printed by.
【００１２】[0012]
【発明の効果】上記のとおり、本発明によれば、第１ス
テップＱ１において、ＲＲ間隔データを入力し、第２ス
テップＱ２において、該入力したＲＲ間隔データのうち
でスペクトル分析を行うべき範囲を、数値を入力するこ
とにより、指定し、第３ステップＱ３において、該指定
した範囲内でスペクトル分析を行い、第４ステップＱ４
において、その分析結果を出力するようにしたことを特
徴とするＲＲ間隔スペクトル分析方法である第１発明
（図１）、及びＲＲ間隔データを入力する入力部１Ａ
と、該入力したＲＲ間隔データのうちでスペクトル分析
を行うべき範囲を数値を入力することにより指定する指
定部１Ｅと、該指定部１Ｅにより指定した範囲内でスペ
クトル分析を行なうスペクトル分析部１Ｃと、該スペク
トル分析部１Ｃにより分析された結果を出力する出力部
１Ｄとから成ることを特徴とするＲＲ間隔スペクトル分
析装置である第２発明（図２）という技術的手段が講じ
られた。上記構成によれば、従来のように、２本のカー
ソルを移動させることにより、スペクトル分析対象範囲
を指定していたのに比較すれば、指定がより一層迅速に
行えることになった。このため、ＲＲ間隔スペクトル分
析自体も早く行うことができ、その分、自律神経障害も
早期に発見されることとなった。即ち、スペクトル分析
対象範囲の指定を迅速に行うことにより、ＲＲ間隔スペ
クトル分析を早く行い、以て自律神経障害の早期発見に
資するという技術的効果を奏することとなった。As described above, according to the present invention, in the first step Q1, the RR interval data is input, and in the second step Q2, the range in which the spectrum analysis is to be performed within the input RR interval data is selected. , By inputting numerical values, in the third step Q3, the spectrum analysis is performed within the specified range, and the fourth step Q4
In the first invention (FIG. 1), which is an RR interval spectrum analysis method characterized by outputting the analysis result, and an input unit 1A for inputting RR interval data.
And a designating section 1E for designating by inputting a numerical value a range in which spectrum analysis is to be performed in the inputted RR interval data, and a spectrum analyzing section 1C for performing spectrum analysis within the range designated by the designating section 1E. The technical means of the second invention (FIG. 2), which is an RR interval spectrum analyzer, is provided which comprises an output unit 1D for outputting the result analyzed by the spectrum analysis unit 1C. According to the above configuration, by moving two cursors as in the conventional case, the spectrum analysis target range is designated, but in comparison, the designation can be performed more quickly. Therefore, the RR interval spectrum analysis itself can be performed earlier, and autonomic neuropathy can be detected earlier accordingly. That is, by rapidly designating the spectrum analysis target range, the RR interval spectrum analysis can be quickly performed, thereby contributing to the early detection of autonomic neuropathy.
【図１】第１発明の実施例を示す図である。FIG. 1 is a diagram showing an embodiment of a first invention.
【図２】第２発明の実施例を示す図である。FIG. 2 is a diagram showing an embodiment of a second invention.
【図３】本発明の動作説明図である。FIG. 3 is an operation explanatory diagram of the present invention.
【図４】本発明の動作説明図である。FIG. 4 is an explanatory diagram of the operation of the present invention.
【図５】本発明の動作説明図である。FIG. 5 is an operation explanatory diagram of the present invention.
【図６】本発明の動作説明図である。FIG. 6 is an operation explanatory diagram of the present invention.
１Ａ 入力部 １Ｅ 指定部 １Ｃ スペクトル分析部 １Ｄ 出力部 1A Input section 1E Designating section 1C Spectrum analysis section 1D Output section
Claims (9)
ータを入力し、第２ステップＱ２において、該入力した
ＲＲ間隔データのうちでスペクトル分析を行うべき範囲
を、数値を入力することにより、指定し、第３ステップ
Ｑ３において、該指定した範囲内でスペクトル分析を行
い、第４ステップＱ４において、その分析結果を出力す
るようにしたことを特徴とするＲＲ間隔スペクトル分析
方法。1. In a first step Q1, RR interval data is input, and in a second step Q2, a range for spectral analysis in the input RR interval data is specified by inputting a numerical value. In the third step Q3, the spectrum analysis is performed within the designated range, and in the fourth step Q4, the analysis result is output, and the RR interval spectrum analysis method.
心電図に基づいて解析されたものである請求項１記載の
ＲＲ間隔スペクトル分析方法。2. The RR interval spectrum analysis method according to claim 1, wherein the RR interval data is analyzed based on a long-term electrocardiogram of 24 hours.
Ｒ−Ｒタコグラムが表示されている画面上において、数
値を入力することにより、指定する請求項１記載のＲＲ
間隔スペクトル分析方法。3. The range in which the spectral analysis is to be performed
The RR according to claim 1, wherein the RR is designated by inputting a numerical value on the screen on which the RR tachogram is displayed.
Interval spectrum analysis method.
項１記載のＲＲ間隔スペクトル分析方法。4. The RR interval spectrum analysis method according to claim 1, wherein the analysis result is displayed.
項１記載のＲＲ間隔スペクトル分析方法。5. The RR interval spectrum analysis method according to claim 1, wherein the analysis result is printed.
と、該入力したＲＲ間隔データのうちでスペクトル分析
を行うべき範囲を数値を入力することにより指定する指
定部１Ｅと、該指定部１Ｅにより指定した範囲内でスペ
クトル分析を行なうスペクトル分析部１Ｃと、該スペク
トル分析部１Ｃにより分析された結果を出力する出力部
１Ｄとから成ることを特徴とするＲＲ間隔スペクトル分
析装置。6. An input section 1A for inputting RR interval data.
And a designating section 1E for designating by inputting a numerical value a range in which spectrum analysis is to be performed in the inputted RR interval data, and a spectrum analyzing section 1C for performing spectrum analysis within the range designated by the designating section 1E. And an output unit 1D for outputting the result analyzed by the spectrum analysis unit 1C.
装置２から入力されたものである請求項６記載のＲＲ間
隔スペクトル分析装置。7. The RR interval spectrum analyzer according to claim 6, wherein the RR interval data is input from a long-term electrocardiogram analyzer 2.
求項６記載のＲＲ間隔スペクトル分析装置。8. The RR-interval spectrum analyzer according to claim 6, wherein the output unit 1D is a display unit 1D1.
請求項６記載のＲＲ間隔スペクトル分析装置。9. The RR interval spectrum analyzer according to claim 6, wherein the output unit 1D is a printer 1D2.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4232772A JPH0654815A (en) | 1992-08-07 | 1992-08-07 | Method and device for rr interval spectral analysis |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP4232772A JPH0654815A (en) | 1992-08-07 | 1992-08-07 | Method and device for rr interval spectral analysis |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| JPH0654815A true JPH0654815A (en) | 1994-03-01 |
Family
ID=16944502
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP4232772A Pending JPH0654815A (en) | 1992-08-07 | 1992-08-07 | Method and device for rr interval spectral analysis |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JPH0654815A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1323376A2 (en) | 2001-12-27 | 2003-07-02 | Kiyoshi Takizawa | Method and apparatus for ECG diagnosis |
| JP2007517553A (en) * | 2004-01-16 | 2007-07-05 | コンプメディクス リミテッド | Apparatus and method for monitoring, detecting and classifying sleep breathing disorders derived from electrocardiograms |
| US8036733B2 (en) | 2004-05-20 | 2011-10-11 | Kiyoshi Takizawa | Diagnostic parameter calculation method, system for diagnosis and diagnostic program |
-
1992
- 1992-08-07 JP JP4232772A patent/JPH0654815A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP1323376A2 (en) | 2001-12-27 | 2003-07-02 | Kiyoshi Takizawa | Method and apparatus for ECG diagnosis |
| JP2007517553A (en) * | 2004-01-16 | 2007-07-05 | コンプメディクス リミテッド | Apparatus and method for monitoring, detecting and classifying sleep breathing disorders derived from electrocardiograms |
| JP4753881B2 (en) * | 2004-01-16 | 2011-08-24 | コンプメディクス リミテッド | Apparatus and signal processing method for monitoring, detecting and classifying sleep disordered breathing derived from electrocardiogram |
| US8036733B2 (en) | 2004-05-20 | 2011-10-11 | Kiyoshi Takizawa | Diagnostic parameter calculation method, system for diagnosis and diagnostic program |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US4736199A (en) | Method of selectively displaying data | |
| JP4493311B2 (en) | Method and apparatus for performing interactive annotation and measurement functions of time series data with automatic marker sequence creation | |
| Ocorr et al. | Semi-automated optical heartbeat analysis of small hearts | |
| CN112842351B (en) | Dynamic electrocardiogram analysis method and system | |
| JPH0654815A (en) | Method and device for rr interval spectral analysis | |
| CN115708685A (en) | Electrocardiogram analysis method based on three-dimensional electrocardio-stack diagram | |
| JPH06296597A (en) | Analyzing method of brain waves and apparatus therefor | |
| JP3129544B2 (en) | RR interval spectrum analysis method and apparatus | |
| JP3378311B2 (en) | RR interval spectrum analysis method and apparatus | |
| JP2746239B2 (en) | Biological information recording device | |
| JPH06105818A (en) | Method and apparatus for preserving result of rr interval spectral analysis | |
| JP2795867B2 (en) | Plant monitoring equipment | |
| JP3510917B2 (en) | Spectrophotometer | |
| JP3024159B2 (en) | Grid display method | |
| JPH0437464B2 (en) | ||
| JPH06351048A (en) | Video signal analysis display device | |
| JPH05131A (en) | Electrocardiogram diagnosis analysis system for user individually | |
| KR102668103B1 (en) | Cardiac ischemia detection technique and apparatus using electrocardiogram signal characteristics | |
| JPS59145970A (en) | Displaying method of spectrum | |
| JPH0312128A (en) | Electrocardiogram analyzing device | |
| JPH07141311A (en) | Image processor | |
| JP3619823B2 (en) | Spectrophotometer | |
| CN114668406A (en) | Electrocardiogram analysis method and device, electronic equipment and computer-readable storage medium | |
| Deneault | Comparing Geometric Object and Conventional Displays in Patient Monitoring | |
| JPH0555125B2 (en) |