• 专利附录
  • 专利说明
  • 权利要求
  • 类似技术
  • 同族专利
  • 引用文献
  • 法律状态
  • 优先权
    • 专利摘要:
    A system for measuring heart rate variability (HRV) of a patent comprises: recording means for obtaining and regulating heartbeat-to-heartbeat intervals for a predetermined period of time; processing means for digitizing said intervals, forming a recurrence plot, and assigning a unit mass to each point on the plot representing a measured interval, and calculating the determinant by the expression Q det =Q xx Q yy wherein: Qxx is the quadrupole moment relative to the X axis of the principal coordinate, Qyy is the quadrupole moment relative to the Y axis of the principal coordinate; and Qdet is the product of Qxx and Qyy.
    公开号:US20010008954A1
    申请号:US09/760,594
    申请日:2001-01-16
    公开日:2001-07-19
    发明作者:Jacob Levitan；Meir Lewkowicz
    申请人:LEV-EL DIAGNOSTICS OF HEART DISEASE Ltd；
    IPC主号:A61B5-0245
    专利说明:
    [0001] The present invention relates to a method and a system for measuring heart rate variability (HRV). [0001] BACKGROUND OF THE INVENTION
    [0002] The study of HRV has been in use for many years as part of clinical, prognostic work; there are international guidelines for evaluating conventional HRV parameters. The conventional parameters are partly frequency domain parameters (power spectra), and partly time domain parameters (various RMS estimates). These methods, though in general successful, are not always conclusive. [0002]
    [0003] Over the last few years, new methods of analyzing R-R intervals of PQRS plots representative of the human heartbeat wave have appeared, all of them showing improved diagnostic and prognostic performance. It has been shown that so-called scale dependent methods outperform scale-independent measure with respect to separating healthy subjects from patients suffering from certain cardiac dysfunctions. But in clinical practice, it is of interest to examine whether, within a group of heart patients, one can extract a subgroup of patents who are at risk, e.g., with respect to sudden cardiac death, rather than to verify the presumably known fact that they do not belong to a group of healthy subjects. It has been shown that, while scale-dependent methods worked in the former case, one had to use scale-independent measures in the latter case. [0003]
    [0004] In practical medicine, recurrence plots of the R-R intervals are used for diagnostic purposes by visual inspection. Since the density of the points is ignored in a visual presentation of the plot, sometimes similar patterns are found for recordings with different HRV's. Although some attempts have been made to include the density of points, this procedure is performed manually and is thus dependent on the performer's skins. Hence, crucial information about the topology of the recurrence plot might be lost. [0004]
    [0005] FIG. 1 illustrates an ECO signal wave. The electro-physiological features of the heart are generally measured by an electrocardiograph, and the electro-physiological recording of the heart function is known as the ECG or EKG. The six features P, Q, R, S, T, U (FIG. 1) describe the sequence of two cycles wherein the R potential is the highest peak. It is therefore easy to distinguish the other five features and the R-peak of the next sequence. The R-R distance is measured in milliseconds and represents the inverse heart rate (HR). The HR is normally not constant, but continually oscillates around its mean level. These short-term cyclic changes are primarily caused by cardiac autonomic modulation. [0005]
    [0006] The calculation of HR and its variability can be used to estimate autonomic activity as and in particular, to evaluate autonomic nervous system influences on heart functions. The autonomic nervous system (ANS) comprises all of the efferent nerves through the visceral organs, including the cardiovascular system, the and the peripheral involuntary muscles. The ANS is generally described as a combination of two main systems that balances and interacts; the sympathetic, regulated by adrenergic activities, and the parasympathetic, cholinergically regulated. One of the main nerves controlling the activity of the heart is the fast-acting, parasympathetic Vagus nerve. [0006]
    [0007] It is generally accepted today that HRV measurement is also a valuable took for the determination of the status of the ANS. Changes in vagal activity cause immediate large changes in instantaneous HR, whereas changes in sympathetic activities are associated with more gradual, slow changes. [0007]
    [0008] The measurement of HR and its rhythmicity, HRV, are only used as a diagnostic tool in cardiology. A stable heart rate is a sign that the heart does not respond to external influences, which responses are mainly regulated by the ANS. Such a situation is dangerous for the individual and is considered to be a pathological symptom. Research has indicated that a quantification of HRV, the discrete beat-to-beat variability in the heart, plays an important prognostic role as an indicator of risk associated with a large variety of diseases, behavioral disorders, mortality and also aging, independent of other risk factors. [0008]
    [0009] Depressed, low HRV has been shown to be a powerful predictor of cardiac events after myocardial infarct. It is therefore crucial to establish a measure of HRV and to quantitatively classify the HRVs of different pathological cases, in order to discriminate between healthy HR profiles and those of patients at risk. [0009]
    [0010] The commercially available medical device for detection of HVR is the Holter 24-hour recording and analysis instrument. A Holter instrument monitor continuously records heart patterns from electrodes attached to the patient for a 24-hour period. The Holter recording technique records the ECG on analog magnetic tape, and a Holter scanner analyzes the tape 60 or 120 times to produce a final report. A Holter scanner report may contain statistical calculations of the heart activity and a detailed report of abnormal cardiological events, such as sinus pauses and propped beats. A limited factor of this technology is the long, 24-hour measurement time and the lack of a graphic print-out summarizing the entire examination period as well as a material evaluation of the total measurement. [0010]
    [0011] U.S. Pat. No. 5,682,901 (Kamen) discloses a method and apparatus for measuring autonomic activity of a patient during a short duration. The method utilizes a visual description of the recurrence plot and separates between direct pathological states, according to different patterns. The method, however, suffers from the fact that the figures of the recurrence plots do not allow the inspection of the density of the data points, which varies over the whole contour, but only renders information of the general shape of the plot. The method and apparatus of Kamen include performing a calculation of the correlation dimension in order to quantify the degree of heart failure, but such a calculation necessitates a dimension that is partly bases on a visual, subjective inspection. Such a visual inspection is known to be unreliable. [0011] SUMMARY AND OBJECTS OF THE INVENTION
    [0012] It is therefore a broad object of the present invention to provide an accurate method and system for measuring HRV over a relatively short period of time of 60 minutes or less. [0012]
    [0013] It is a further object of the present invention to provide a method and a system for measuring HRV, giving results which are more precise and easier to evaluate. [0013]
    [0014] It is a still object of tie present invention to provide a method and system enabling a prognosis of the status of a patient with a history of heart failure or impaired heart function, by quantification of the degree of heart failure or heart function impairment. [0014]
    [0015] It is a yet further object of the present invention to provide a method and a system allowing for the classification of patients with a history of heart failure or impaired heart function into the following three groups, ordered according to the risk of death due to heart failure: (1) patients with a minimal risk of sudden death, comparable to that of healthy individuals; (2) patients with an increased risk of sudden death, and (3) patients with a high risk of sudden death. [0015]
    [0016] In accordance with the present invention, there is therefore provided a system for measuring heart rate variability (HtV) of a patient, comprising recording means for obtaining and recording heartbeat-to-heartbeat intervals for a predetermined period of time; process means for digitizing said intervals, forming a recurrence plot, and assigning a unit mass to each point on the plot representing a measured interval, and calculating the determinant by the expression[0016]
    [0017] wherein: [0017]
    [0018] Q[0018] xx is the quadrople moment relative to the X axis of the principal coordinate,
    [0019] Q[0019] yy is the quadrople moment relative to the Y axis of the principal coordinate; and
    [0020] Q[0020] det is the product to Qxx and Qyy.
    [0021] The invention further provides a method for measuring the heart rate variability (HRV) of a patient, comprising collecting data of heatbeat-to-heartbeat intervals; determining the intervals during a predetermined period of time; generating a recurrence plot from said determined intervals, and calculating the determinant by the expression[0021]
    [0022] wherein: [0022]
    [0023] Q[0023] xx is the quadrupole moment relative to the X axis of the principal coordinate,
    [0024] Q[0024] yy is the quadrupole moment relative to the X axis of the principal coordinate; and
    [0025] Q[0025] det is the product of Qxx and Qyy, BRIEF DESCRIPTION OF THE DRAWINGS
    [0026] The invention will now be described in connection with certain preferred embodiments with reference to the following illustrative figures so that it may be more fully understood. [0026]
    [0027] With specific reference now to the figures in detail, it is stressed that the particulars shown are by way of example and for purposes of illustrative discussion of like preferred embodiments of the present invention only, and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspect of the invention. In this regard, no attempt is made to show structural details of the invention in more detail than is necessary for a fundamental understanding of the invention, the description taken with the drawings making apparent to those skilled in the art how the several forms of the invention may be embodied in practice. [0027]
    [0028] In the drawings: [0028]
    [0029] FIG. 1 is an illustration of a conventional electrocardiograph measurement; [0029]
    [0030] FIG. 2 is a block diagram of the system accord to the present invention for measuring HRV; [0030]
    [0031] FIG. 3 illustrates a recurrence or Poincaré plot; [0031]
    [0032] FIG. 4 is a block diagram of the method of the present invention for measuring HRV; and [0032]
    [0033] FIGS. [0033] 5A-5C are recurrence plots indicative of patients with various degrees of risk due to heart failure. DETAILED DESCRIPTION OF THE EXEMPLARY EMBODIMENTS
    [0034] The system of measuring HRV according to the present invention will now be described with reference to FIG. 2. [0034]
    [0035] The patient is placed in a resting position on a table and electrodes attached to the patient are connected to a recorder [0035] 2, either by wires or by means of wireless technology. Any instruments that can measure and record the R-R distance with the sensitivity and accuracy of 1 msec can be applied, such as, for example, a conventional ECG apparatus, a Holter 24-hour recorder, or a specifically designed instrument. The recording is made for a predetermined period of time, from 15 minutes to 24 hours.
    [0036] The recorded R-R intervals can then be obtained from an ECG by a PQRS detector [0036] 4, or from a specific dedicated recorder, and are transferred to a processing system 6, which may be a personal computer equipped with a software program which detects the R peaks, calculates the R-R intervals between adjacent peek in msec, and processes the data obtained. Optionally, detector 4 may also include terminal 14 for inputting R-R intervals from other sources, and a patient data bank file 16 for fixture use, The plots and calculations, together with clinical and personal data regarding the patient, will be presented on a screen 8; a printer 10 can provide a printout.
    [0037] The R-R interval data is plotted in a two- or three-dimensional recurrence or Poincaré plot [0037] 12 (FIG. 3), as an illustration of the HRV. By using a specific mathematical procedure called the Quadrupole Moment Method (QMM) as will be described below, a certain index or value Qdet, a quantitative measure for the HRV, is calculated.
    [0038] Referring now to FIG. 4, the method according to the present invention will be described. [0038]
    [0039] The patient's ECG is collected, recorded by the recorder [0039] 2 and digitized at 18. From the digitized ECG, the R-R intervals are measured and arranged in a two-dimensional array, e.g., in columns, where one column is the number of the interval, represented by the integers 1, 2, 3 . . . , and the second column is the R-R interval expressed in msec, typically between 500-1500 msec for each R-R interval. The detrending is achieved by calculating the local average value of the R-R time series in a running window, typically of the length 2 m, where m is a positive interger. This average value is subtracted from the R-R series at the center position of the running window; resulting in the locally determined R-R interval time series.
    [0040] The fist R-R interval, together with the second R-R interval is marked as a point on a two-dimensional coordinate system [0040] 20, wherein the X coordinate is the first R-R interval and the Y coordinate is the second interval. The second interval is hereafter marked as the next X coordinate and the consecutive interval as the next Y-coordinate. This procedure is continued, the Y coordinate being partly switched to an X-coordinate in the following point. The procedure starts from the first number in the file and obviously ends with the last. This procedure is termed “detrending,“ meaning that every deviation in the heart rate variability is measured relative to an overall average and there is no due consideration given to short term local fluctuations that are not caused by internal variables of the heart rhythm, but are enforced by external conditions. Detrending becomes especially crucial in measurements over an extended period of time, e.g., 24 hours, and in every measurement where the patients measured are left in uncontrolled conditions. The importance of detrending resides in the accuracy of both prognosis and diagnosis, which depend thereon.
    [0041] The evolving recurrence plot [0041] 22 of scattered points in FIG. 4 exhibits the density of the points. Each data point in the recurrence plot 12 (FIG. 3) is assigned a unit mass. The two-dimensional body is then analyzed in terms of a gravitational, multi-pole expansion Φ:
    [0042] wherein: [0042] φ = ∫p  ( r _ ′ ) r _ - r _ ′    2  r _ ′ ≈ M r + R _ cm · r _r 3 + 1 2  ∑ Q ij  x i  x j r 5 + …   ( 1 )
    [0043] ρ(r′) is the density of mass (here equal to 1) at the point (r′); [0043]
    [0044] {overscore (r)} is the position vector of the observation point; [0044]
    [0045] {overscore (r)}′ is the position vector of the mass at r′; [0045]
    [0046] d[0046] 2{overscore (r)}′ is an infinitesimal area around r′;
    [0047] M is a number of masses; [0047]
    [0048] {overscore (R)}[0048] cm=∫{overscore (r)}′d2{overscore (r)}′=0
    [0049] is the dipole moment; [0049]
    [0050] Q[0050] ij is the normalized quadrupole moment: Q ij = 1 M  ∑ masses  ( 3  x i ′  x j ′ -r ′ 2  δ ij)( 2 )
    [0051] x[0051] i is the X or Y coordinate of the vector {overscore (r)}i; and
    [0052] x[0052] ixj is the product of two such coordinates.
    [0053] For the principal axis: [0053] Q ij = ( Q xx0 0Q yy ) ⇒
    [0054] geometrical description by Q[0054] det≡Qxx·Qyy.
    [0055] The first term in equation (1) is the monopole, and represents the number of data points, which is irrelevant to the analysis of the HRV. The second term is the gravitational dipole moment and vanishes by choosing the origin of the coordinate system in the center of the mass. The first significant term is hence the third term in the equation (1), the quadrupole moment. This moment is expressed as a 2×2 matrix and is rendered independent of the number of data points by dividing it by the number of data points. By diagonalization, the off-diagonal elements of this matrix are set equal to zero. This procedure defines, in fact, a new coordinate axis, the Principal Axis. The two non-vanishing terms are symbolized by Q[0055] xx and Qyy. The expressions are given in equation (2).
    [0056] Finally, the determinant of the matrix Q[0056] det=Qxx·Qyy is calculated at 24. This determinant is one of the parameters which can be used in the further investigation of the state of the heart of the patient whose ECG has been measured.
    [0057] In certain cases, as with different-shaped reoccurrence plots for different detrending scales, the multipole expansion has to be continued. [0057]
    [0058] More comparative parameters are extracted from a higher dimesional recurrence plot, i.e., for n's dimension, the recurrence plot is constructed from the n's dimensional vector, R-R[0058] i, R-Ri1, R-Ri+2, . . . R-Ri+n. Here, the multipole expansion is performed on the solution of the Laplace Equation in n dimensions.
    [0059] Eventually, diagnostic and prognostic evaluation is performed at [0059] 26, on the basis of the measurement results. For example, as seen in FIGS. 5A to 5C, a patient with a minimal risk of sudden death, comparable to that of a healthy individual, will have a recurrence plot as shown in FIG. 5A; a patient with an increased risk will have a recurrence plot as shown in FIG. 5B, and a patient with a high risk will have a recurrence plot as shown in FIG. 5C.
    [0060] It will be evident to those skilled in the art that the invention is not limited to the details of the foregoing illustrated embodiments and that the present invention may be embodied in other specific forms without departing from the spirit or essential attributes thereof. The present embodiments are therefore to be considered in all respects as illustrative and not restrictive, the scope of the invention being indicated by the appended claims rather than by the foregoing description, and all changes which come within the meaning and range of equivalency of the claims are therefore intended to be embraced therein. [0060]
    权利要求:
    Claims (11)
    [1" id="US-20010008954-A1-CLM-00001] 1. A system for measuring heart rate variability (HRV) of a patient, said system comprising:
    recording means for obtaining and recording heartbeat-to-heartbeat intervals for a predetermined period of time;
    processing means for digitizing said intervals, forming a recurrence plot, and assigning a unit mass to each point on the plot representing a measured interval, and calculating the determinant by the expression
    Ti Qdet =Q xx Q yy
    wherein:
    Qxx is the quaduole moment relative to the X axis of the principal coordinate,
    Qyy is the quadrupole moment relative to the Y axis of the principal coordinate; and
    Qdet is the product of Qxx and Qyy.
    [2" id="US-20010008954-A1-CLM-00002] 2. The system as claimed in
    claim 1 , wherein said recording means further comprises an input terminal for inputting data concerning heartbeat-heartbeat intervals obtained from other sources.
    [3" id="US-20010008954-A1-CLM-00003] 3. The system as claimed in
    claim 1 , wherein said recurrence plot is formed on a screen.
    [4" id="US-20010008954-A1-CLM-00004] 4. The system as claimed in
    claim 1 , wherein said recurrence plot is produced by a printer.
    [5" id="US-20010008954-A1-CLM-00005] 5. The system as claimed in
    claim 1 , further comprising a patient data bank file for future use.
    [6" id="US-20010008954-A1-CLM-00006] 6. A method for measuring the heat rate variability (HRV) of a patient, said method comprising:
    collecting data of heartbeat-to-heartbeat intervals;
    determining the intervals during a predetermined period of time;
    generating a recurrence plot from said determined intervals, and
    calculating the determinant by the expression
    Q det =Q xx Q yy
    wherein:
    Qxx is the quadrupole moment relative to the X axis of the principal coordinate,
    Qyy is the quadrupole moment relative to the Y axis of the principal coordinate; and
    Qdet is the product of Qxx and Qyy.
    [7" id="US-20010008954-A1-CLM-00007] 7. The method as claimed in
    claim 6 , further comprising digitizing said intervals prior to generating a reoccurrence plot.
    [8" id="US-20010008954-A1-CLM-00008] 8. The method as claimed in
    claim 6 , further comprising exhibiting said recurrence plot on a screen.
    [9" id="US-20010008954-A1-CLM-00009] 9. The method as claimed in
    claim 6 , further comprising printing said recurrence plot.
    [10" id="US-20010008954-A1-CLM-00010] 10. The method as claimed in
    claim 6 , evaluating the state of the palient's sympathetic and parasympathetic nervous system, based on the HRV measurements obtained.
    [11" id="US-20010008954-A1-CLM-00011] 11. The method as claimed in
    claim 6 , further comprising transferring the HRV measurements and calculated determinant to a patient data bank file.
    类似技术:
    公开号 | 公开日 | 专利标题
    US7412283B2|2008-08-12|Method and system for processing electrocardial signals
    US6389308B1|2002-05-14|System and device for multi-scale analysis and representation of electrocardiographic data
    US5159932A|1992-11-03|Myocardial ischemia detection system
    US5609158A|1997-03-11|Apparatus and method for predicting cardiac arrhythmia by detection of micropotentials and analysis of all ECG segments and intervals
    US6390986B1|2002-05-21|Classification of heart rate variability patterns in diabetics using cepstral analysis
    EP1639497B1|2007-07-18|Method and apparatus for extracting causal information from a chaotic time series
    EP0489010B1|1998-11-18|Cardiac biopotential analysis system and method
    US6731974B2|2004-05-04|Method and system for measuring heart rate variability
    US6280395B1|2001-08-28|System and method for determining muscle dysfunction
    US20030093002A1|2003-05-15|Function indicator for autonomic nervous system based on phonocardiogram
    Gitter et al.1995|Fractal analysis of the electromyographic interference pattern
    US7221976B2|2007-05-22|Analysis of the alternans cycle to cycle and/or the variability of the ventricular repolarization wave in an ECG signal
    JP2002523163A|2002-07-30|Systems and methods to aid clinical judgment
    EP0468999A1|1992-02-05|Cerebral biopotential analysis system and method.
    US20090076402A1|2009-03-19|Method and apparatus for determining susceptibility for arrhythmias using wedensky modulated electrocardiography and alternans analysis
    US6607480B1|2003-08-19|Evaluation system for obtaining diagnostic information from the signals and data of medical sensor systems
    US6491629B1|2002-12-10|Method for determining at least one diagnostic piece of information from signal patterns of medical sensor systems
    CN104321013A|2015-01-28|Method and apparatus for providing a visual representation of sleep quality based on ECG signals
    US7367950B1|2008-05-06|Distributed apparatus and method for non-invasive measurement of functional state and adaptive response
    US20020156392A1|2002-10-24|Method and apparatus for inspecting biological rhythms
    US6223073B1|2001-04-24|Noninvasive method for identifying coronary disfunction utilizing electrocardiography derived data
    JP2004350952A|2004-12-16|Autonomic nerve activity evaluation device, its method, and autonomic nerve activity evaluation program
    EP1146433A1|2001-10-17|Method and apparatus for detecting sleep disorders using the variability of the RR interval in an ECG waveform
    Borun et al.1966|Computer analysis of Frank lead electrocardiographic data recorded in an epidemiologic study
    EP0896282A1|1999-02-10|A noninvasive method for identifying coronary disfunction utilizing electrocardiography derived data
    同族专利:
    公开号 | 公开日
    US6731974B2|2004-05-04|
    IL134123D0|2001-04-30|
    JP4668505B2|2011-04-13|
    JP2004523250A|2004-08-05|
    CA2397086A1|2001-07-26|
    AU2701801A|2001-07-31|
    CA2397086C|2005-10-18|
    WO2001052715A2|2001-07-26|
    IL134123A|2005-06-19|
    WO2001052715A3|2002-01-17|
    EP1118945A1|2001-07-25|
    引用文献:
    公开号 | 申请日 | 公开日 | 申请人 | 专利标题
    US20060178588A1|2005-01-03|2006-08-10|Lee Brody|System and method for isolating effects of basal autonomic nervous system activity on heart rate variability|
    US20080004539A1|2002-03-01|2008-01-03|Christine Ross|Analysis of heart rate variability data in animals for health conditions assessment|
    US20080125666A1|2004-09-16|2008-05-29|Stuart Crozier|Medical Monitoring System|
    WO2012017432A1|2010-08-05|2012-02-09|Lev-El Diagnostics Of Heart Disease Ltd.|Method for measuring heart rate variability|
    US8620418B1|2013-01-04|2013-12-31|Infobionic, Inc.|Systems and methods for processing and displaying patient electrocardiograph data|
    US20190110711A1|2017-10-12|2019-04-18|Nihon Kohden Corporation|Apparatus and method for processing physiological information|US4018219A|1975-11-03|1977-04-19|American Home Products Corporation|Heart rate variability measurement|
    US5769793A|1989-09-08|1998-06-23|Steven M. Pincus|System to determine a relative amount of patternness|
    US5201321A|1991-02-11|1993-04-13|Fulton Keith W|Method and apparatus for diagnosing vulnerability to lethal cardiac arrhythmias|
    US5333615A|1992-06-22|1994-08-02|William Craelius|Apparatus for digitally recording and analyzing electrocardial and other bioelectric signals|
    DE4320519A1|1993-06-22|1995-01-05|Max Delbrueck Centrum|Method of diagnosing cardiac dysrhythmias|
    WO1995003739A1|1993-08-03|1995-02-09|Peter Walter Kamen|A method of measuring autonomic activity of a patient|
    JP3196871B2|1994-09-30|2001-08-06|アトムメディカル株式会社|Diagnostic device|
    US5560368A|1994-11-15|1996-10-01|Berger; Ronald D.|Methodology for automated QT variability measurement|
    US5609158A|1995-05-01|1997-03-11|Arrhythmia Research Technology, Inc.|Apparatus and method for predicting cardiac arrhythmia by detection of micropotentials and analysis of all ECG segments and intervals|
    US5755671A|1995-10-05|1998-05-26|Massachusetts Institute Of Technology|Method and apparatus for assessing cardiovascular risk|
    US5749900A|1995-12-11|1998-05-12|Sulzer Intermedics Inc.|Implantable medical device responsive to heart rate variability analysis|
    US6035233A|1995-12-11|2000-03-07|Intermedics Inc.|Implantable medical device responsive to heart rate variability analysis|
    US5995868A|1996-01-23|1999-11-30|University Of Kansas|System for the prediction, rapid detection, warning, prevention, or control of changes in activity states in the brain of a subject|
    US6026320A|1998-06-08|2000-02-15|Cardiac Pacemakers, Inc.|Heart rate variability as an indicator of exercise capacity|US7031765B2|2002-11-11|2006-04-18|Medtronic, Inc|Algorithms for detecting atrial arrhythmias from discriminatory signatures of ventricular cycle lengths|
    PT1477114E|2003-05-13|2009-04-06|Gme Rechte Und Beteiligungen G|Device and method for detecting atrial fibrillation|
    DE10351728A1|2003-10-31|2005-06-30|W.O.M. World Of Medicine Ag|Method and system for the diagnosis and / or monitoring of the cardiovascular system of a living being|
    US8002711B2|2005-03-18|2011-08-23|Respironics, Inc.|Methods and devices for relieving stress|
    BRPI0508896A|2004-03-18|2007-09-11|Helicor Inc|stress relieving methods and devices|
    WO2005109322A2|2004-05-04|2005-11-17|Rush University Medical Center|Methods using recurrence quantification analysis to analyze and generate images|
    WO2010070640A1|2008-12-15|2010-06-24|Ariel Medical K.L. Ltd.|Method and system for measuring heart rate variability|
    US9901266B2|2015-03-06|2018-02-27|Greatbatch Ltd.|Apparatus and method for personalized cardiac signal processing|
    法律状态:
    2001-01-16| AS| Assignment|Owner name: LEV-EL DIAGNOSTICS OF HEART DISEASE LTD., ISRAEL Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:LEVITAN, JACOB;LEWKOWICZ, MEIR;REEL/FRAME:011494/0814 Effective date: 20010114 |
    2004-04-15| STCF| Information on status: patent grant|Free format text: PATENTED CASE |
    2007-10-11| FPAY| Fee payment|Year of fee payment: 4 |
    2011-12-19| REMI| Maintenance fee reminder mailed|
    2012-04-30| FPAY| Fee payment|Year of fee payment: 8 |
    2012-04-30| SULP| Surcharge for late payment|Year of fee payment: 7 |
    2015-11-03| FPAY| Fee payment|Year of fee payment: 12 |
    优先权:
    申请号 | 申请日 | 专利标题
    IL13412300A|IL134123A|2000-01-19|2000-01-19|Method and system for measuring heart rate variability|
    IL134,123||2000-01-19||
    IL134123||2000-01-19||
    [返回顶部]