从零搭建自动化心电图系统：仪表放大器、双T陷波滤波与LabVIEW心率检测

1. 项目概述与核心价值

心电图（ECG）系统，听起来像是医院里那些庞大、精密的专业设备，离我们很遥远。但你知道吗？其核心原理——捕捉并放大我们身体表面那微乎其微的“生物电信号”——其实完全可以用一些基础的电子元件和软件工具来实现。我自己动手搭建过几套这样的系统，从最初的信号微弱到几乎看不见，到后来能清晰稳定地捕捉到每一次心跳的P波、QRS波群和T波，这个过程充满了挑战，也收获了许多教科书里学不到的实战经验。今天，我就来拆解一个我实际验证过的、结构相对简洁但效果可靠的自动化心电图系统方案。它仅使用一个仪表放大器、两个滤波器（陷波和低通），再配合LabVIEW进行数据采集与分析，就能实现从信号拾取到心率自动计算的全流程。无论你是电子爱好者、生物医学工程的学生，还是对生理信号测量感兴趣的创客，这个项目都能让你亲手触摸到“生命信号”的脉搏，理解医疗电子设备最底层的工作原理。

这个系统的核心价值在于其“麻雀虽小，五脏俱全”的完整性。它没有使用复杂的集成芯片或昂贵的专业模块，而是从最基础的运放、电阻、电容开始，让你透彻理解每一级电路的作用：仪表放大器如何从强噪声背景中提取出微伏级的心电信号；陷波滤波器如何精准“挖掉”恼人的50/60Hz工频干扰；低通滤波器又如何像筛子一样滤除肌肉颤动等高频噪声。最后，通过LabVIEW将调理后的信号可视化，并编写算法自动识别R波、计算心率，完成从模拟电路到数字处理的闭环。整个过程，你会遇到增益计算偏差、噪声抑制不足、软件触发不稳定等实际问题，而解决这些问题的过程，正是从理论走向实践的关键。

2. 系统整体架构与设计思路

2.1 信号链路与模块化设计

一个完整的心电图采集系统，其信号链路可以清晰地划分为三个核心阶段：前端模拟信号调理、模数转换与数据采集、后端数字信号处理与显示。我们这个项目聚焦于前两个阶段，并用LabVIEW串联起整个流程。

第一级：仪表放大器（Instrumentation Amplifier, INA）。这是整个系统的“门户”和“第一道放大器”。心电信号通过体表电极拾取时，其幅度通常在0.5mV到5mV之间（QRS波群峰值），并且淹没在高达数伏特的共模干扰（如工频干扰、电极接触噪声）中。普通运放电路对此无能为力。仪表放大器的设计精髓在于其极高的输入阻抗（避免从人体汲取电流，影响信号）和极高的共模抑制比（CMRR），能像“火眼金睛”一样，只放大两个输入端之间的微小电压差（差模信号），而对两个输入端共有的干扰电压（共模信号）进行强力抑制。我们采用经典的三运放结构来实现它，这是理解生物电放大器的必修课。

第二级：陷波滤波器（Notch Filter）。经过放大，信号中的差模成分被凸显，但一个最强的干扰源——50Hz或60Hz的工频干扰——依然存在。这是因为我们的身体就像一个天线，会感应到环境中无处不在的交流电辐射。陷波滤波器，也叫带阻滤波器，它的任务就是在频率响应曲线上“挖”一个非常窄的“坑”，专门针对这个特定频率进行深度衰减。在我们的设计中，目标是消除60Hz干扰（根据地区电网频率选择50Hz或60Hz），让这个频率点的信号增益急剧下降，从而在波形上几乎看不到工频的“毛刺”。

第三级：低通滤波器（Low-Pass Filter）。心电信号的有效频率成分主要集中在0.05Hz到150Hz之间（诊断级ECG会更高，但我们这个简易系统关注核心波形，设定在350Hz以下已足够）。然而，环境中的射频干扰、肌肉电信号（EMG）等噪声频率远高于此。低通滤波器的作用就是设置一个“关卡”，只允许低于某个截止频率（如335Hz）的信号通过，高于此频率的一律大幅衰减。这能有效平滑波形，让P、QRS、T波的特征更加清晰。

提示：这个三级级联的顺序是经过深思熟虑的。放大器必须先工作，因为原始信号太微弱，直接滤波会导致信号被噪声淹没而无法有效提取。先放大，再逐级滤除特定噪声，是生物电信号处理的黄金流程。

2.2 关键设计参数与目标

在动手之前，我们必须明确整个系统的性能目标，这决定了后续每一个元器件的选型和参数计算：

1. 总电压增益：目标是将约1mV的心电信号（QRS波）放大到1V左右，以便于后续的模数转换器（ADC）采集和软件显示。因此，需要的总电压增益约为1000倍（60dB）。
2. 带宽：系统需要能通过0.5Hz（用于观察基线漂移和低频成分）到350Hz（包含心电信号的主要谐波）的信号。这由各级电路的频率响应共同决定。
3. 噪声抑制：
   • 对60Hz工频干扰的抑制能力（陷波滤波器深度）至少达到-40dB，即干扰幅度被衰减到原来的1%。
   • 对350Hz以上频率的抑制（低通滤波器滚降）应足够陡峭，以滤除大部分高频噪声。
4. 输出阻抗：最后一级电路（低通滤波器输出）的输出阻抗应足够低，以驱动后续的数据采集（DAQ）设备，通常要求小于1kΩ。

明确了这些目标，我们就可以开始逐一设计每个电路模块了。

3. 核心电路模块详解与设计

3.1 仪表放大器：微弱信号的“金耳朵”

我们采用经典的三运放仪表放大器结构。它由两级组成：第一级是两个同相放大器，提供高输入阻抗和差分增益；第二级是一个差分放大器，将差分信号转换为单端输出，并进一步提供增益，同时抑制共模信号。

元件选型与参数计算：

• 运放：选择通用、低噪声的精密运放，如TI的TL07x系列或Analog Devices的OP07系列。它们具有较高的输入阻抗和较低的偏置电流，适合生物电应用。我们需要三个这样的运放。

• 增益电阻：增益主要由第一级两个运放外围的电阻R1、R2和第二级的电阻R3、R4决定。整个电路的差分增益公式为：Gain = (1 + 2*R2/R1) * (R4/R3)这个公式是设计的核心。为了获得约1000倍的总增益，我们需要合理分配两级增益。通常，第一级增益不宜过高，以避免共模信号被放大到超出运放输出范围。一个常见的策略是让第一级增益在10-100之间，第二级增益在10-100之间，二者相乘得到总增益。

  在参考设计中，作者使用了R1=3.3kΩ, R2=33kΩ, R3=1kΩ, R4=48kΩ（原文中“48 ohms”应为笔误，结合增益计算，应为48kΩ）。代入公式： 第一级增益 = 1 + 2*(33k/3.3k) = 1 + 20 = 21 第二级增益 = 48k / 1k = 48 总增益 = 21 * 48 = 1008 这与目标增益1000倍非常接近。

• 仿真验证：在LTspice中搭建电路，进行AC交流扫描分析（AC Sweep）。输入一个1mV的差分小信号，扫描频率从1Hz到1kHz。观察输出信号的幅度。在关键频段（如10Hz），增益应为1000倍（即60dBV）。同时，可以施加一个共模信号（如1V， 60Hz），观察输出端信号的大小，来评估共模抑制比（CMRR）。一个设计良好的仪表放大器，其共模输出应远小于差模输出。

实操心得：电阻的精度直接影响增益的准确性和CMRR。尽量使用1%精度的金属膜电阻。对于R1和R2，如果使用两个独立的电阻，它们的匹配度至关重要，不匹配会严重降低CMRR。在实际制作中，可以考虑使用高精度的电阻网络或通过筛选匹配电阻对。

3.2 双T型陷波滤波器：精准狙击工频干扰

工频干扰是生物电测量中最顽固的敌人。我们选择有源双T型陷波滤波器，因为它能提供非常深的阻带和较高的Q值（选择性），且电路相对简单。

工作原理：双T型网络由电阻和电容构成一个桥式电路，在中心频率处产生相消干涉，使得该频率信号被极大衰减。加入运放构成有源形式，可以缓冲输出并调整Q值。

参数计算（以消除60Hz为例）： 中心频率公式为：f0 = 1 / (2 * π * R * C)为了计算方便，通常先选定电容C的值。选择0.1uF（即0.1×10^-6 F）是一个常见且电容值容易获取的选择。 则计算电阻R：R = 1 / (2 * π * 60 * 0.1e-6) ≈ 26.5 kΩ这是理论值。在经典的双T有源电路中，为了调节Q值（决定阻带宽度），电阻网络会稍有变化。参考设计中的计算更详细：它通过设定Q=8，计算出了两个关键电阻R1和R2。

• R1 = 1/(2*Q*2*π*f0*C) ≈ 1.66 kΩ，选用1.8kΩ标准值。
• R2 = 2*Q/(2*π*f0*C) ≈ 424 kΩ，通过串联390kΩ和33kΩ得到423kΩ。
• 反馈电阻Rf用于设定通带增益为1（单位增益），计算约为1.79kΩ，选用1.8kΩ。

仿真验证：在LTspice中，输入一个0.1V、频率从1Hz扫到1kHz的正弦波。观察波特图。你会在60Hz处看到一个非常尖锐的下陷（陷波）。测量该点增益，理想情况下应低于-40dB（即0.01倍）。这意味着一个0.1V的60Hz干扰输入，输出只有1mV，干扰被有效抑制。

注意事项：双T网络的性能极度依赖于R和C的匹配。如果用于产生两个“T”的电阻电容对不匹配，陷波频率会偏移，深度会变浅。因此，务必使用高精度、低温漂的元件（如1%电阻，C0G/NP0材质的电容）。在实际焊接时，可以考虑使用可调电阻进行微调，以精准对准50Hz或60Hz。

3.3 一阶低通滤波器：滤除高频“杂音”

心电信号的主要能量集中在低频。我们采用简单有效的一阶有源低通滤波器（巴特沃斯响应），其特点是通带平坦，在截止频率之后以-20dB/十倍频程的斜率衰减。

参数计算： 截止频率公式：f_c = 1 / (2 * π * R * C)我们的目标是让350Hz以下的信号通过。为了留有余量，将截止频率设定在335Hz。 同样，先选定电容C=0.1uF。 计算电阻R：R = 1 / (2 * π * 335 * 0.1e-6) ≈ 4.75 kΩ选用最接近的标准值4.7kΩ。

电路实现：这是一个简单的反相或同相结构。为了保持单位增益并简化设计，常采用同相电压跟随器加RC网络的结构，或者使用Sallen-Key拓扑。一阶滤波器虽然滚降慢，但对于滤除远高于350Hz的噪声（如广播射频）已经足够，且电路简单，相位失真小。

仿真验证：在LTspice中进行AC扫描。在335Hz处，增益应下降至-3dB（约0.707倍）。在1kHz时，增益应显著降低，观察是否达到-10dB以下。

实操心得：如果你发现心电波形上仍有较多高频“毛刺”，可以考虑将截止频率降低到250Hz甚至150Hz。但要注意，过低的截止频率会衰减心电信号本身的高频成分（如QRS波的一些细节），使波形变得圆滑。这是一个需要权衡的地方。在实际调试中，可以通过更换不同阻值的电阻来快速调整截止频率。

4. 从仿真到实物：电路搭建与调试实录

4.1 元器件准备与焊接要点

在完成仿真并确认电路性能符合预期后，就可以着手准备实物了。

物料清单：

• 集成电路：UA741或TL081运放至少3个（建议多备几个）。注意运放的电源引脚（V+， V-）和输出引脚。
• 电阻：根据计算，需要准备多种阻值：3.3kΩ， 33kΩ， 1kΩ， 48kΩ， 1.8kΩ， 390kΩ， 33kΩ（用于串联）， 4.7kΩ。所有电阻建议采用1/4瓦、1%精度的金属膜电阻。
• 电容：0.1uF陶瓷电容或C0G电容若干（用于陷波和低通滤波器）。注意电容的电压等级和材质，对于滤波电路，稳定性好的电容是关键。
• 电源：运放需要双电源供电（如+9V和-9V）。可以使用两个9V电池串联构成中心接地，或者使用专用的双路输出直流稳压电源。
• 其他：面包板、跳线、连接线、电极片（至少3个）、用于连接DAQ设备的BNC或夹子线。

搭建顺序与技巧：

1. 分模块搭建：不要在面包板上一次性连完所有电路。先搭建仪表放大器，单独测试。用函数发生器产生一个几毫伏、几赫兹的小差分信号输入，用示波器观察输出是否被正确放大1000倍左右，且没有自激振荡。
2. 电源去耦：这是保证电路稳定工作的重中之重。在每个运放的电源引脚（V+和V-）附近，紧贴着引脚，跨接一个0.1uF的陶瓷电容到地（电源地）。这能为运放提供瞬间的电流补偿，抑制电源线上的噪声。
3. 连接陷波滤波器：将仪表放大器的输出接入陷波滤波器的输入。用函数发生器产生一个60Hz、幅度较大的正弦波（模拟干扰），叠加在有用信号上，观察输出端60Hz成分是否被大幅抑制。
4. 连接低通滤波器：最后接入低通滤波器。用函数发生器产生一个高频信号（如1kHz），观察其是否被有效衰减。
5. 整体联调：将所有模块连接起来。此时，可以用函数发生器模拟一个标准的心电波形（频率1-2Hz，幅度1-5mV的类QRS波），观察最终输出波形是否清晰、放大倍数是否正确、基线是否平稳。

4.2 常见故障与排查指南

在实物搭建中，几乎一定会遇到问题。以下是我踩过的一些“坑”及解决方法：

5. LabVIEW自动化采集与分析实现

硬件电路产生了干净的模拟心电信号，接下来就需要LabVIEW将其数字化、可视化，并自动计算心率。

5.1 数据采集（DAQ）配置

你需要一块数据采集卡（DAQ），如National Instruments的USB-6000系列，或者一些兼容的声卡（通过音频输入口，但性能有限）。在LabVIEW中，通常使用“DAQ Assistant”或“DAQmx”函数来配置采集任务。

关键参数设置：

• 物理通道：选择你连接电路输出信号的模拟输入通道（如ai0）。
• 测量类型：选择“电压”。
• 接线端配置：根据你的电路，选择“差分”（Differential）或“参考单端”（RSE）。如果电路输出是相对于地的单端信号，通常选RSE。
• 采样率：根据奈奎斯特定理，至少是信号最高频率的两倍。心电信号我们关心350Hz以下，采样率设为1000 Hz（1 kS/s）是安全且常见的。
• 采样数/持续时间：例如，每次读取3000个样本，在1000 Hz采样率下，就对应3秒的波形。这样可以在前面板上显示一段连续的波形。

5.2 心率检测算法设计

这是LabVIEW程序的核心。心电图中，R波是幅度最高、斜率最大的波，最容易被检测。这里实现一个简单的幅度阈值检测法。

程序框图逻辑：

1. 数据获取：DAQ Assistant读取电压数组。
2. 预处理（可选但推荐）：对原始数据数组进行一个软件低通滤波（例如截止频率40Hz的Butterworth滤波器），以进一步平滑波形，便于检测。
3. 寻找峰值：
   • 计算当前数据段的最大值（Max）和最小值（Min）。
   • 设定一个动态阈值，例如：阈值 = Min + 0.7 * (Max - Min)。这个系数（0.7）需要根据实际波形调整，要高于T波和P波，但低于R波峰值。
   • 遍历电压数组，当一个数据点V[i]的值高于阈值，并且V[i]大于前一个点V[i-1]且大于后一个点V[i+1]时，认为检测到一个R波峰值。
4. 记录与计算：
   • 记录下每个峰值点对应的数组索引（或时间戳）。
   • 计算相邻峰值之间的时间间隔（ΔT）。ΔT = (索引差) / 采样率。
   • 将时间间隔转换为心率：心率 (BPM) = 60 / ΔT。
5. 输出与显示：
   • 将计算出的瞬时心率或平均心率显示在前面板的数值控件中。
   • 将原始或滤波后的波形数组送入“波形图”（Waveform Graph）控件显示。

实操心得：单纯的幅度阈值法在信号有波动或噪声时容易误检或漏检。一个更稳健的方法是结合斜率（差分）信息。计算信号的差分（V[i] - V[i-1]），R波上升沿的斜率极大。可以设定一个斜率阈值，只有当幅度超过阈值且斜率超过另一个阈值时，才判定为R波。这在LabVIEW中可以通过“差值”函数和比较逻辑轻松实现，能显著提高检测鲁棒性。

5.3 前面板布局与用户体验

一个清晰的前面板能让测试过程直观高效。

• 主波形图：占据主要区域，实时显示心电波形。可以设置自动缩放，或固定Y轴范围（如-1.5V 到 +1.5V）。
• 心率显示：用大号字体显示实时计算出的心率（BPM）。
• 控制按钮：添加“开始采集”、“停止”、“保存数据”等按钮。
• 参数调节（高级）：可以加入旋钮或输入框，让用户实时调整心率检测的幅度阈值和斜率阈值，以适应不同人的信号强度。
• 状态指示：用一个LED指示灯闪烁来对应每次检测到的心跳，增加直观性。

6. 系统测试与性能评估

6.1 使用函数发生器进行模拟测试

在连接真人测试前，务必先用函数发生器进行系统验证，这能排除人体变量，聚焦于电路和软件本身的性能。

1. 模拟静息心率：将函数发生器设置为输出频率1Hz（对应60 BPM）、峰峰值20mV的正弦波或更接近真实QRS波的三角波。将其输出连接到电路输入端（模拟差分信号）。观察LabVIEW显示的波形是否清晰，幅度是否被放大到约20V（增益1000倍），显示的心率是否在60 BPM左右。
2. 模拟运动后心率：将频率改为2Hz（对应120 BPM）。观察心率显示是否相应变化。
3. 抗干扰测试：在有用信号（1Hz， 20mV）上，叠加一个60Hz、100mV的干扰信号。观察最终输出波形中，60Hz的纹波是否被有效抑制。对比打开和关闭陷波滤波器（或将其旁路）时的波形差异，你会直观感受到滤波器的威力。

6.2 人体实测与结果分析

这是最激动人心的环节。使用一次性心电电极片，按照标准导联I的位置放置：正极（红色）接右腕，负极（白色）接左腕，地线（黑色）接右踝（或左踝）。确保皮肤清洁，电极粘贴牢固。

1. 静息状态测量：让受试者静坐放松。启动LabVIEW采集。你应该能看到清晰、稳定的心电波形，包含P波、QRS波群和T波。记录下稳定的心率值。在我的测试中，静息心率大约在70-80 BPM，与设计案例中的75 BPM相符。
2. 运动后测量：让受试者进行几分钟的原地高抬腿或开合跳。运动后立即连接电极测量。此时波形幅度可能变小（因出汗导致电极接触阻抗变化），噪声可能增大，心率会显著加快。记录心率值。案例中从75 BPM升至137 BPM是合理的。
3. 误差分析：
   • 幅度误差：预期QRS波幅度约1V（1mV * 1000倍），实测可能为0.9V。这可能是电阻精度、运放增益带宽积限制、电极接触阻抗等因素共同导致的。10%左右的误差在分立元件搭建的系统中是可以接受的。
   • 心率误差：与标准心率计对比。误差可能来源于心率检测算法的灵敏度（阈值设置）、波形噪声导致的误判等。通过优化LabVIEW中的检测算法（如加入斜率判断、去抖动逻辑），可以显著提高精度。

通过这个从零搭建的简易自动化心电图系统，你不仅获得了一个能工作的设备，更重要的是，你深入理解了生物电信号放大、滤波、采集和处理的完整链条。每一个环节的参数选择、每一个问题的排查，都是宝贵的工程经验。你可以在此基础上继续探索：增加导联数量、设计更高阶的滤波器、实现更复杂的 arrhythmia（心律失常）检测算法，或者将系统移植到嵌入式平台（如STM32+OLED）做成真正的便携设备。