1. Learning Objectives
By the end of this study guide, the student should be able to relate the physiology of the cardiovascular system to the engineering principles used to measure, model, and monitor it.
Specific Objectives
- Describe the anatomy of the heart and the structure of the circulatory system (systemic and pulmonary circuits).
- Explain the electrical conduction system of the heart and the ionic basis of the cardiac action potential.
- Interpret a standard 12-lead electrocardiogram (ECG) waveform — P wave, QRS complex, T wave, and intervals.
- Explain the cardiac cycle in terms of pressure–volume relationships, valve action, and heart sounds.
- Apply engineering models (pump–resistor–capacitor analogies, fluid dynamics) to describe blood flow and blood pressure.
- Identify the principles of operation of cardiovascular instrumentation: ECG electrodes, pressure transducers, pulse oximetry, and ultrasound Doppler flow measurement.
- Apply signal-processing techniques — filtering, amplification, noise removal (50/60 Hz mains interference), and heart-rate-variability (HRV) analysis — to biosignals.
- Relate common cardiovascular disorders (arrhythmias, hypertension, myocardial infarction) to their ECG signatures and diagnostic measurements.
Table of Contents
- 2. Overview of the Cardiovascular System
- 3. Anatomy of the Heart and Blood Vessels
- 4. The Heart's Electrical Conduction System
- 5. The Cardiac Cycle
- 6. The Electrocardiogram (ECG)
- 7. Hemodynamics — an Engineering View
- 8. Cardiovascular Instrumentation
- 9. Signal Processing of Cardiac Signals
- 10. Common Cardiovascular Disorders
- 11. Interactive Simulations
- 12. Self-Assessment Quiz
- 13. Summary & Study Tips
- 14. References and Further Reading
2. Overview of the Cardiovascular System
The cardiovascular system is the body's transport network. It consists of the heart (a muscular pump), blood vessels (arteries, veins, capillaries), and blood (the transport medium). Its primary functions are:
- Transport: delivering oxygen and nutrients to tissues and removing carbon dioxide and metabolic wastes.
- Regulation: distributing hormones, regulating body temperature, and maintaining fluid and pH balance.
- Protection: carrying immune cells and clotting factors that defend against injury and blood loss.
2.1 The Two Circuits
Right side of the heart → pumps deoxygenated blood to the lungs via the pulmonary arteries → gas exchange in the alveolar capillaries → oxygenated blood returns to the left atrium via the pulmonary veins. Low-pressure circuit (~25/8 mmHg).
Left side of the heart → pumps oxygenated blood through the aorta to all body tissues → oxygen and nutrients diffuse out; CO₂ and wastes diffuse in → deoxygenated blood returns via the superior and inferior vena cava to the right atrium. High-pressure circuit (~120/80 mmHg).
2.2 Why This Matters to an Electrical & Communication Engineer
3. Anatomy of the Heart and Blood Vessels
3.1 Gross Anatomy of the Heart
The heart is a fist-sized muscular organ (~300 g) located in the mediastinum between the lungs, slightly tilted to the left. It has four chambers:
| Chamber | Function | Walls / Characteristics |
|---|---|---|
| Right Atrium (RA) | Receives deoxygenated blood from the body via vena cava | Thin-walled; contains the sinoatrial (SA) node |
| Right Ventricle (RV) | Pumps blood to the lungs (pulmonary circulation) | Moderate wall thickness; crescent-shaped |
| Left Atrium (LA) | Receives oxygenated blood from the lungs via pulmonary veins | Thin-walled |
| Left Ventricle (LV) | Pumps blood to the entire body (systemic circulation) | Thickest wall (~10–15 mm); conical |
3.2 Valves of the Heart
- Atrioventricular (AV) valves: Tricuspid valve (RA→RV) and mitral (bicuspid) valve (LA→LV). Prevent backflow into the atria during ventricular contraction. Anchored by chordae tendineae to papillary muscles.
- Semilunar valves: Pulmonary valve (RV→pulmonary artery) and aortic valve (LV→aorta). Prevent backflow from the arteries into the ventricles.
3.3 Layers of the Heart Wall
- Epicardium — outer protective layer (visceral pericardium).
- Myocardium — the cardiac muscle layer; the thickest in the left ventricle; the contractile "engine" of the pump.
- Endocardium — inner smooth lining in contact with blood.
3.4 Blood Vessels
| Vessel Type | Structure & Function | Engineering Analogy |
|---|---|---|
| Arteries | Thick elastic/muscular walls carrying blood away from the heart; withstand pulsatile pressure; damp pressure oscillations | Pressure-rated pipelines with surge tanks |
| Arterioles | Small arteries whose smooth muscle regulates flow — the main site of peripheral resistance | Control valves / variable resistors |
| Capillaries | One-cell-thick walls; site of O₂, CO₂, nutrient and waste exchange | Mass-transfer membranes / heat exchangers |
| Venules & Veins | Thin-walled, low-pressure return lines; contain one-way valves; aided by muscle-pump action; ~64% of blood volume resides here | Low-pressure return lines with check valves |
4. The Heart's Electrical Conduction System
The heart is myogenic — it generates its own electrical impulses without nervous stimulation. The conduction system is a specialized network of excitable cells:
- Sinoatrial (SA) node — the natural pacemaker (~60–100 impulses/min), located at the junction of the right atrium and superior vena cava. Its cells spontaneously depolarize (funny current, If).
- Atrial conduction pathways — impulses spread across both atria, causing atrial contraction.
- Atrioventricular (AV) node — the only electrical connection between atria and ventricles; introduces a deliberate delay (~0.1 s) allowing ventricular filling.
- Bundle of His — conducts the impulse from the AV node into the ventricular septum.
- Right and Left Bundle Branches — the left branch splits into anterior and posterior fascicles.
- Purkinje fibres — rapidly distribute the impulse (~4 m/s) throughout the ventricular myocardium, triggering coordinated ventricular contraction from apex upward.
4.1 The Cardiac Action Potential
Unlike the short neuronal action potential, the ventricular myocyte action potential lasts ~250–300 ms and has five phases:
| Phase | Ionic Events | Membrane Behaviour |
|---|---|---|
| Phase 0 — Rapid depolarization | Fast Na⁺ channels open; Na⁺ rushes in | Sharp upstroke (to ~ +20 mV) |
| Phase 1 — Early repolarization | Na⁺ channels inactivate; transient K⁺ efflux (Ito) | Small notch |
| Phase 2 — Plateau | Ca²⁺ influx through L-type Ca²⁺ channels balanced by K⁺ efflux | Prolonged plateau (~200 ms); triggers Ca²⁺-induced Ca²⁺ release → contraction |
| Phase 3 — Repolarization | Ca²⁺ channels close; K⁺ efflux continues (IKr, IKs) | Return toward resting potential |
| Phase 4 — Resting potential | K⁺ equilibrium (~ −90 mV) maintained by Na⁺/K⁺ ATPase and inward rectifier K⁺ current (IK1) | Stable resting potential |
4.2 The Electrophysiology of Automaticity
Pacemaker cells of the SA and AV nodes have an unstable Phase 4 — a slow spontaneous depolarization (the pacemaker potential) caused by the inward "funny" Na⁺ current (If) and T-type Ca²⁺ currents, followed by L-type Ca²⁺-mediated upstrokes. The SA node's intrinsic rate (~100/min) is normally slowed to ~70/min by vagal (parasympathetic) tone.
5. The Cardiac Cycle
The cardiac cycle is the sequence of mechanical events — contraction (systole) and relaxation (diastole) — driven by the electrical events of Section 4. One cycle lasts about 0.8 s at a heart rate of 75 bpm.
5.1 Phases of the Cycle
| Phase | Events | Valve Status (AV / SL) |
|---|---|---|
| 1. Atrial systole | Atrial contraction tops up ventricular filling ("atrial kick", ~10–20% of filling) | Open / Closed |
| 2. Isovolumetric contraction | Ventricles contract; pressure rises sharply; all valves closed; volume constant | Closed / Closed |
| 3. Rapid ejection | Ventricular pressure exceeds aortic/pulmonary pressure; semilunar valves open; blood ejected rapidly | Closed / Open |
| 4. Reduced ejection | Ejection slows as pressures equilibrate | Closed / Open |
| 5. Isovolumetric relaxation | Ventricles relax; pressure falls; all valves closed | Closed / Closed |
| 6. Rapid filling | Ventricular pressure falls below atrial; AV valves open; passive filling (~70–80%) | Open / Closed |
| 7. Diastasis | Slow final filling as pressures equalize | Open / Closed |
5.2 Pressure–Volume Loop
The left ventricular pressure–volume loop is the standard graphical description of the cycle:
- Bottom edge: ventricular filling (pressure low, volume rising from ~50 to ~120 mL).
- Right edge: isovolumetric contraction (vertical line up at constant volume ~120 mL).
- Top edge: ejection (volume falls, pressure peaks ~120 mmHg).
- Left edge: isovolumetric relaxation (vertical line down at ~50 mL).
The area enclosed by the loop equals the stroke work per beat. Key volumes:
| Parameter | Typical Value (adult at rest) |
|---|---|
| End-diastolic volume (EDV) | ~120 mL |
| End-systolic volume (ESV) | ~50 mL |
| Stroke volume (SV = EDV − ESV) | ~70 mL |
| Cardiac output (CO = HR × SV) | ~5 L/min |
| Ejection fraction (EF = SV/EDV) | 55–70% (key clinical indicator) |
5.3 Heart Sounds
- First heart sound (S1, "lub"): closure of the AV valves at the start of systole — lower pitch, longer.
- Second heart sound (S2, "dub"): closure of the semilunar valves at the start of diastole — higher pitch, shorter.
- Third (S3) and Fourth (S4) sounds: may be pathological (ventricular filling gallop, atrial contraction against stiff ventricle).
- Murmurs: turbulent flow through defective (stenotic or regurgitant) valves.
6. The Electrocardiogram (ECG)
The ECG is a recording of the summed electrical activity of the heart as detected by electrodes on the body surface. It is the single most important non-invasive cardiac diagnostic tool.
6.1 The Standard ECG Waveform
| Feature | Represents | Normal Duration / Amplitude |
|---|---|---|
| P wave | Atrial depolarization (initiated by SA node) | < 0.12 s; < 0.25 mV |
| PR interval | Atrial depolarization + AV nodal delay | 0.12 – 0.20 s |
| QRS complex | Ventricular depolarization (septum → apex → walls) | < 0.12 s; largest amplitude |
| ST segment | Period when ventricles are fully depolarized (plateau phase) | Isoelectric (baseline) |
| T wave | Ventricular repolarization | Same polarity as QRS, broader |
| QT interval | Total ventricular depolarization + repolarization (rate-dependent) | < ~0.44 s (corrected: QTc < 0.44 s) |
6.2 The 12-Lead ECG
A standard ECG records 12 "leads" — different electrical views of the heart (each lead is a differential voltage between electrode pairs):
- Limb leads (I, II, III, aVR, aVL, aVF): bipolar and augmented unipolar leads in the frontal plane. Lead II is the classical rhythm strip.
- Precordial (chest) leads V1–V6: placed across the chest; view the heart in the horizontal plane.
By examining which leads show changes, clinicians localize the affected region of the heart (e.g., inferior wall MI shows changes in II, III, aVF).
6.3 ECG Paper and Rate Calculation
- Standard paper speed: 25 mm/s; each small square (1 mm) = 0.04 s; each large square (5 mm) = 0.2 s.
- Heart rate (bpm) = 300 ÷ number of large squares between R waves (e.g., 4 large squares → 75 bpm), or 1500 ÷ small squares.
- Standard calibration: 1 mV deflection = 10 mm (10 mm/mV).
7. Hemodynamics — an Engineering View
7.1 Blood Pressure
Arterial blood pressure is reported as systolic/diastolic (e.g., 120/80 mmHg). Pulse pressure = systolic − diastolic; mean arterial pressure (MAP) ≈ diastolic + (pulse pressure/3).
| Category (adult) | Systolic (mmHg) | Diastolic (mmHg) |
|---|---|---|
| Normal | < 120 | < 80 |
| Pre-hypertension | 120–139 | 80–89 |
| Hypertension (Stage 1) | 140–159 | 90–99 |
| Hypertension (Stage 2) | ≥ 160 | ≥ 100 |
7.2 Ohm's Law Analogy
- Pressure (ΔP) ↔ Voltage (V)
- Flow (Q) ↔ Current (I)
- Resistance (R) ↔ Resistance (R)
7.3 Poiseuille's Law
For laminar flow in a vessel, the resistance is:
R = 8ηL / (πr⁴) ⇒ Q = πΔP r⁴ / (8ηL)
where η = blood viscosity, L = vessel length, r = vessel radius. Note the fourth-power dependence on radius — a 50% decrease in radius increases resistance 16-fold. Blood is non-Newtonian and flow in small vessels is pulsatile and partly turbulent, so Poiseuille's law is an idealization.
7.4 Laplace's Law
Wall tension in a spherical/ellipsoidal vessel or chamber: T = P × r / (2h) (for a sphere) where P = internal pressure, r = radius, h = wall thickness. This explains why a dilated failing heart (large r) or a dilated aneurysm (high P × r) is at risk of wall rupture.
7.5 Cardiac Output and Its Regulation
CO = HR × SV (≈ 5 L/min at rest, up to 25 L/min in athletes). Stroke volume is governed by:
- Preload (end-diastolic volume — ventricular filling) — the Frank–Starling law: greater stretch → stronger contraction (within limits).
- Afterload (aortic pressure the ventricle must overcome to eject).
- Contractility (intrinsic strength, increased by sympathetic stimulation, decreased by hypoxia/acidosis).
Feedback control: Baroreceptors in the aortic arch and carotid sinus detect pressure changes and, via the medulla and autonomic nerves, adjust heart rate (vagus/sympathetic) and vascular resistance — a classic closed-loop feedback system with negative feedback maintaining blood pressure homeostasis.
8. Cardiovascular Instrumentation
This section links each physiological quantity to the sensor and electronic front-end used to measure it — the core of bio-instrumentation for ECE students.
8.1 ECG Systems
| Block | Function & Typical Specifications |
|---|---|
| Electrodes (Ag/AgCl) | Convert ionic current in tissue to electronic current; require low offset potential and stable contact resistance |
| Instrumentation amplifier | Gain ~1000; CMRR > 100 dB; input impedance > 10 MΩ; bias currents < 100 nA |
| Right-leg drive | Inverts and feeds back common-mode voltage to the body, cancelling 50 Hz interference |
| Filters | High-pass ~0.05 Hz (removes baseline drift); low-pass ~150 Hz (anti-alias, removes EMG); notch 50/60 Hz (mains) |
| Isolation | Patient safety — withstand defibrillator surges and leakage limits (IEC 60601) |
| ADC / DSP | Sampling ≥ 500 Hz per lead; QRS detection, ST analysis, arrhythmia classification |
8.2 Blood Pressure Measurement
- Auscultatory (sphygmomanometer): cuff inflated above systolic pressure, deflated slowly; Korotkoff sounds heard via stethoscope — first sound = systolic, disappearance = diastolic.
- Oscillometric (automated digital): cuff pressure oscillations — mean pressure at maximum oscillation; systolic/diastolic estimated by empirical ratios. Used in most modern monitors.
- Intra-arterial (invasive): catheter + strain-gauge pressure transducer (Wheatstone bridge) — the gold standard, continuous waveform.
8.3 Pulse Oximetry
A non-invasive optical technique: red (660 nm) and infrared (940 nm) LEDs shine through a fingertip, and a photodiode measures transmitted intensity. Because oxyhemoglobin and deoxyhemoglobin absorb the two wavelengths differently, the ratio of pulsatile (AC) to non-pulsatile (DC) absorbance at the two wavelengths yields SpO₂. The AC component is caused by arterial blood pulsation — an elegant lock-in style demodulation of a physiological signal from strong DC ambient and venous components.
8.4 Doppler Ultrasound
Blood velocity is measured via the Doppler shift: fd = 2f₀ v cosθ / c, where f₀ = transmitted frequency, v = blood velocity, θ = beam–flow angle, c = speed of sound. Continuous-wave (for high-velocity jets) and pulsed-wave Doppler (depth-resolved) are used; cardiac output can be estimated by integrating velocity across the aortic valve cross-section.
8.5 Pacemakers and Defibrillators
- Cardiac pacemaker: an implanted pulse generator delivering timed voltage pulses (typically ~5 V, 0.5–1 ms, via bipolar leads) to maintain an adequate heart rate. Sensing circuits detect intrinsic R waves to inhibit or trigger pacing on demand — a real-time control system.
- Implantable Cardioverter-Defibrillator (ICD): monitors rhythm and delivers high-energy (~30 J) shocks to terminate lethal ventricular arrhythmias.
9. Signal Processing of Cardiac Signals
9.1 ECG Signal Characteristics and Noise
| Signal Component | Approximate Frequency Range |
|---|---|
| QRS complex (diagnostic content) | 5 – 50 Hz |
| P and T waves | 0.5 – 10 Hz |
| Full diagnostic bandwidth | 0.05 – 150 Hz |
| Noise / Artifact | Frequency / Character | Mitigation |
|---|---|---|
| Mains (50/60 Hz) interference | 50 Hz + harmonics, common-mode | High CMRR amp, right-leg drive, notch filter |
| Baseline wander | < 0.5 Hz (respiration, electrode motion) | High-pass filter ≥ 0.05 Hz |
| Electromyogram (EMG) muscle noise | 20 – 500 Hz, random | Low-pass filter, patient relaxation |
| Electrode contact motion | Sudden spikes/offsets | Good skin preparation, Ag/AgCl gel electrodes |
9.2 QRS Detection
The basis of real-time heart-rate monitoring. Classical algorithm (Pan–Tompkins, 1985):
- Band-pass filtering (5–15 Hz) to emphasize QRS energy.
- Differentiation — highlights the steep QRS slope.
- Squaring — makes all values positive, emphasizes high frequencies.
- Moving-window integration — smooths, giving a pulse-like feature at each beat.
- Adaptive thresholding — two-level (signal/ noise) thresholds with refractory period (~200 ms) reject false detections.
9.3 Heart Rate Variability (HRV)
Beat-to-beat (R–R) intervals are not constant; HRV quantifies this variation, typically from 5-min or 24-h ECG records. Frequency-domain measures mirror spectral estimation:
- HF band (0.15–0.4 Hz): modulated by parasympathetic (vagal) activity and respiration.
- LF band (0.04–0.15 Hz): both sympathetic and parasympathetic influences.
- LF/HF ratio: used as an index of sympatho-vagal balance.
Time-domain measures include SDNN (standard deviation of NN intervals) and RMSSD. Low HRV is associated with stress, diabetes, and poor outcomes after MI.
9.4 Biotelemetry and mHealth
Modern cardiovascular monitoring uses Bluetooth Low Energy (BLE) to stream ECG/PPG from wearables to smartphones and cloud servers — an application of packet communications, compression (lossless ECG compression algorithms), and network QoS (reliable, low-latency delivery of alarm-class data). Kenya's growing telemedicine infrastructure makes this a highly relevant career area.
10. Common Cardiovascular Disorders
| Condition | Physiological Basis | Electrical / Measurement Signature |
|---|---|---|
| Myocardial infarction (MI) | Coronary artery occlusion → ischemia → necrosis of myocardium | ST-segment elevation (STEMI) or depression + T-wave inversion; pathological Q waves develop; troponin biomarkers rise |
| Atrial fibrillation (AF) | Chaotic atrial activity (~350–600/min); irregular AV conduction | No P waves; irregularly irregular R–R intervals; fibrillatory baseline |
| Heart block (AV block) | Impaired AV conduction | Prolonged PR (1st degree); dropped beats (2nd degree); AV dissociation (3rd degree) |
| Ventricular fibrillation (VF) | Disorganized ventricular activity; no effective pumping | Chaotic waveform of varying shape/amplitude; cardiac arrest — requires immediate defibrillation |
| Hypertension | Chronically elevated SVR and/or CO | BP ≥ 140/90 mmHg repeatedly; LV hypertrophy on ECG/echo |
| Valvular stenosis / regurgitation | Narrowed or leaky valve → murmur; pressure overload or volume overload | Murmurs on auscultation/phonocardiogram; characteristic ECG and echo changes |
| Long-QT syndrome | Delayed ventricular repolarization (ion-channel mutations) | Prolonged QT/QTc > 0.44 s; risk of torsades de pointes |
11. Interactive Simulations
11.1 ECG Waveform Simulator (Synthetic Lead II)
This simulator synthesizes a normal sinus-beat ECG using the piecewise Gaussian model (McSharry et al.). Adjust the heart rate and observation window; estimate the rate using the 300 method and check your answer.
11.2 R–R Interval & HRV Calculator
Enter a series of R–R intervals (in seconds, comma-separated) from a rhythm strip. The tool computes mean rate, SDNN (overall variability) and RMSSD (short-term beat-to-beat variability).
11.3 Cardiac Output Calculator (CO = HR × SV)
Compare: a trained athlete at rest may have SV ~ 100 mL at HR 50 → 5 L/min with far less metabolic cost.
12. Self-Assessment Quiz
Attempt each question before revealing the answer.
Q1. Which structure is the normal pacemaker of the heart, and what is its intrinsic rate range?
Q2. What electrical event does each of the following represent: P wave, QRS complex, T wave?
Q3. Why is the PR interval clinically important, and what is its normal range?
Q4. An ECG strip (25 mm/s) shows an R–R interval of 4 large squares. What is the heart rate, and what formula did you use?
Q5. Using the hydraulic–electric analogy, state the relationship between mean arterial pressure, cardiac output, and systemic vascular resistance.
Q6. Why does a small reduction in arteriole radius cause such a large increase in vascular resistance?
Q7. List four sources of noise/artifact in a surface ECG and one mitigation technique for each.
Q8. What is the Frank–Starling law of the heart?
Q9. Explain how a pulse oximeter estimates arterial oxygen saturation using two wavelengths.
Q10. Give the ECG hallmark of (a) STEMI, (b) atrial fibrillation, (c) third-degree AV block.
Q11. A patient's echocardiogram gives EDV = 120 mL and ESV = 60 mL at HR = 70 bpm. Calculate SV, CO, and ejection fraction.
Q12. Why does an implantable pacemaker need both sensing and pacing functions? What engineering constraint dominates its power budget?
13. Summary & Study Tips
13.1 Key Takeaways
- The heart is a myogenic dual pump: SA node → AV node (delay) → His–Purkinje → synchronized ventricular contraction.
- Cardiac action potentials have a Ca²⁺-mediated plateau (Phase 2) linking excitation to contraction and creating the long refractory period.
- The cardiac cycle couples electrical timing to mechanical pumping; the pressure–volume loop quantifies pump work; CO = HR × SV.
- The ECG is a differential millivolt signal whose waves map to depolarization/repolarization events; standard speed 25 mm/s, 300-method for rate.
- Hemodynamics obeys engineering laws: Ohm's-law analogue (MAP = CO × SVR), Poiseuille's r⁴ law, and Laplace's law.
- Cardiovascular instruments — ECG front-ends, oscillometric BP monitors, pulse oximeters, Doppler flow meters, pacemakers — apply instrumentation, signal processing, and communications principles directly.
13.2 Study Tips
- Draw the ECG on grid paper repeatedly until you can sketch P–QRS–T with correct timings from memory.
- Trace one impulse through the conduction system while narrating the ECG waves it produces — link electrophysiology to the waveform.
- Practice calculations: heart rate from R–R, SV/CO/EF from volumes, MAP from BP, resistance changes from Poiseuille's law.
- Build the analogy table (pressure↔voltage, flow↔current, compliance↔capacitance) and solve simple circuit equivalents of circulation.
- Use the simulations in Section 11: estimate rates visually before reading the value — this trains the clinical speed you need in labs.
- Read real ECGs: use open ECG databases (e.g., PhysioNet MIT-BIH) and apply the Pan–Tompkins pipeline you can implement in MATLAB/Python.
14. References and Further Reading
- Guyton, A. C. & Hall, J. E., Textbook of Medical Physiology, 14th ed., Elsevier — Chapters on the heart as a pump and cardiac arrhythmias.
- Klabunde, R. E., Cardiovascular Physiology Concepts, 3rd ed., Wolters Kluwer.
- Webster, J. G. (ed.), Medical Instrumentation: Application and Design, 4th ed., Wiley — chapters on biopotential electrodes, ECG, and blood-pressure measurement.
- Malmivuo, J. & Plonsey, R., Bioelectromagnetism: Principles and Applications of Bioelectric and Biomagnetic Fields, Oxford University Press (online edition).
- Pan, J. & Tompkins, W. J., "A Real-Time QRS Detection Algorithm," IEEE Trans. Biomed. Eng., BME-32(3):230–236, 1985.
- McSharry, P. E. et al., "A dynamical model for generating synthetic electrocardiogram signals," IEEE Trans. Biomed. Eng., 50(3):289–294, 2003.
- Task Force of the European Society of Cardiology, "Heart rate variability: standards of measurement, physiological interpretation and clinical use," Circulation, 93:1043–1065, 1996.
- PhysioNet MIT-BIH Arrhythmia Database — open ECG data for practice: https://physionet.org.
- Semmlow, J. L., Biosignal and Medical Image Processing, 3rd ed., CRC Press.