1. 项目概述从零搭建一个ECG信号采集电路如果你对生物医学电子感兴趣或者想亲手做一个能“听见”心跳的电路那么这个基于运算放大器的ECG采集项目会是一个绝佳的实践。心电图信号极其微弱通常在毫伏级别并且淹没在强大的工频干扰和肌电噪声中。直接用一个普通放大器去放大得到的只会是一团噪音。这个项目的核心就是设计一套“信号调理链”像一条精密的流水线把原始的心电信号一步步提纯、放大最终变成我们能在屏幕上清晰看到的波形甚至计算出心率。整个电路由三个关键模块串联而成仪表放大器INA、60Hz陷波滤波器和150Hz低通滤波器。仪表放大器负责以高共模抑制比放大微弱的差分心电信号陷波滤波器专门“挖掉”50Hz或60Hz的工频干扰低通滤波器则限制信号带宽滤除高频噪声。我们将使用LTspice进行仿真设计验证理论计算然后在面包板上用经典的741运算放大器搭建实物电路最后用Arduino进行数据采集和可视化。这个过程会涉及模拟电路设计的核心从理论计算、仿真优化到硬件实现与调试的全流程。无论你是电子工程的学生还是DIY爱好者跟着走一遍你对运放应用和生物信号处理的理解会上一个台阶。2. 核心电路模块设计与原理剖析2.1 仪表放大器INA如何捕捉微伏级差分信号仪表放大器是处理生物电信号的“门户”。心电信号通过体表电极采集本质上是两个测量点之间的电位差差分信号但同时会叠加一个很大的、对两个测量点都相同的干扰共模信号比如工频干扰。普通运放电路难以区分这两者而仪表放大器的价值就在于其极高的共模抑制比。我们采用经典的三运放仪表放大器结构。它的核心优势在于第一级由两个同相放大器组成提供高输入阻抗避免从人体抽取过多电流第二级是一个差分放大器将第一级输出的差分信号进一步放大并抑制共模信号。其闭环增益公式为G (1 2R2/R1) * (R4/R3)。在这个项目中目标增益设为1000倍60dB。我们首先选定R120kΩ然后根据公式计算其他电阻值。这里有一个关键点公式中的R2和R3需要匹配R4和R5也需要匹配在差分放大级任何失配都会直接降低共模抑制比。在仿真和实际搭建时我们应尽量使用精度为1%的金属膜电阻。注意为什么仿真增益约1600会偏离理论值1000这很可能源于仿真模型中运放的非理想特性如有限开环增益、输入偏置电流以及我们计算时忽略了第二级差分放大器的电阻匹配条件。在实际电路中使用高精度电阻和低失调电压的运放如OP07、AD620等专用芯片是减少误差的关键。对于741这类通用运放失调电压较大不适合直接用于放大毫伏级信号本项目主要用于原理验证。2.2 60Hz陷波滤波器精准剔除电源干扰工频干扰50Hz或60Hz取决于地区是生物电采集中最顽固的噪声源其幅度可能比心电信号大几十倍。陷波滤波器也叫带阻滤波器其频率响应曲线在特定频率点出现一个尖锐的“凹陷”能将该频率点的信号大幅衰减。我们采用双T型有源陷波滤波器。它由一个无源双T网络和一个运放构成。无源双T网络本身在中心频率f0处具有零输出特性但Q值品质因数很低凹陷很宽。加入运放构成正反馈后可以显著提高Q值让凹陷变得又深又窄只精准消除60Hz干扰而不影响附近频率的心电信号如QRS波的主要能量在5-20Hz。中心频率f0的计算公式为f0 1 / (2πRC)。为了简化设计通常先选取方便易得的电容值如本项目中的0.1μF和0.2μF然后反推电阻值。计算时需注意双T网络要求两组RC乘积相等。提高Q值可以通过调节运放的正反馈量来实现但反馈过强会导致电路不稳定甚至自激振荡。仿真时观察幅频特性曲线一个理想的陷波滤波器应在60Hz处有超过40dB的衰减且凹陷两侧的过渡尽可能陡峭。2.3 150Hz低通滤波器划定信号带宽的边界心电信号的有效频率成分大致在0.05Hz到150Hz之间。设置一个150Hz截止频率的低通滤波器目的是限制信号带宽滤除高频噪声如肌电干扰通常100Hz、电极接触噪声以及可能混入的射频干扰。我们选择二阶赛伦-凯Sallen-Key低通滤波器拓扑。这种结构使用一个运放和少量RC元件就能实现每十倍频程-40dB的衰减斜率性能优于简单的一阶RC滤波器。其设计关键在于确定滤波器的类型如巴特沃斯、切比雪夫。巴特沃斯型具有最平坦的通带响应是我们追求信号保真度的首选。设计公式涉及截止频率f_c和品质因数Q。对于巴特沃斯响应Q值约为0.707。我们已知f_c150Hz可以先设定一个电容值如C10.066μF然后利用标准设计公式或查表法计算电阻值。仿真时进行AC扫描验证-3dB点是否准确落在150Hz并观察通带是否平坦阻带衰减是否足够快。一个设计良好的低通滤波器在150Hz处增益下降3dB在300Hz两倍截止频率处应能有约-12dB的衰减。3. LTspice仿真从理论到虚拟验证3.1 INA模块的仿真与误差分析在LTspice中搭建三运放仪表放大器电路。使用理想的运放模型如UniversalOpamp2进行初步计算验证是个好习惯它能快速检验电路拓扑和增益公式的正确性。按照计算好的电阻值R120k R210k R31k R4500k放置元件。进行AC分析交流小信号分析设置扫描频率从1Hz到10kHz。观察输出波形图通常以dB表示增益。我们会发现在低频段增益曲线并非完全平坦的60dB1000倍而可能有一个峰值且-3dB带宽可能比预期窄。这揭示了第一个非理想性运放的增益带宽积限制。741运放的增益带宽积典型值为1MHz。当闭环增益设为1000时其理论带宽只有1MHz/1000 1kHz。这意味着对于频率高于1kHz的信号增益会开始下降。而我们仿真看到的约1600的增益64dB很可能是因为在测试频率点如100Hz上运放的开环增益仍然足够高但电阻值的微小偏差仿真中为理想值但公式计算可能未考虑差分级匹配被放大所致。实操心得仿真时建议将电阻值从理想值替换为E96系列的标准值如10k用9.76k或10.2k的最近似值并引入1%的容差进行蒙特卡洛分析。这能让你提前看到元件公差对电路性能尤其是共模抑制比的影响。你会发现即使很小的失配也可能导致CMRR从理论上的无穷大下降到80dB甚至更低。3.2 滤波器模块的仿真与参数优化分别搭建双T陷波滤波器和赛伦-凯低通滤波器的仿真电路。对于陷波滤波器AC扫描范围应围绕60Hz如10Hz到200Hz重点关注陷波深度和-3dB带宽。仿真结果可能显示中心频率偏移如63.655Hz这源于计算时π的取值、电容标称值与实际值的差异。我们可以通过微调电阻值将固定电阻改为电位器模型进行参数扫描来“校准”中心频率。对于低通滤波器扫描范围可以更宽1Hz到1kHz。验证截止频率是否为150Hz并观察在通带内150Hz增益是否基本为0dB无放大衰减在阻带衰减斜率是否接近-40dB/dec。如果发现通带有凸起峰值说明Q值过高接近不稳定边缘需要减小与运放正反馈相关的电阻值。关键仿真技巧使用.step指令对关键元件如陷波器的R低通滤波器的C进行参数扫描直观观察元件值变化对频率响应的影响。这能帮你理解电路对元件精度的敏感度为后续硬件选型提供依据。例如你会发现陷波滤波器的中心频率对电容值非常敏感因此在实际电路中应选用高精度、低漂移的C0G或NP0材质电容。3.3 级联系统仿真与动态范围评估将三个模块按INA - 陷波 - 低通的顺序级联起来进行整体仿真。这里需要引入一个重要的概念动态范围和各模块增益分配。假设心电信号最大幅值为1mV目标最终输出幅值为便于观察的1V则系统总增益需要1000倍。全部增益都放在INA是否合理不一定。如果INA增益过高其输出可能在后级滤波器处理大信号时提前接近运放的电源电压±9V导致削波失真。更合理的分配可能是INA增益设为200-500提供初步放大和共模抑制陷波滤波器增益为1单位增益低通滤波器增益为2-5补偿滤波带来的信号损失并调整总增益至1000。在LTspice中使用瞬态分析Transient Analysis输入一个幅值为1mV、频率为1Hz模拟心率并叠加了60Hz、幅值为10mV干扰的正弦波作为测试信号。观察最终输出波形是否干净、有无失真。同时检查每个运放输出节点的电压瞬时值确保其在任何时刻都远离电源轨如保持在±7V以内为信号留出足够的摆幅空间。4. 硬件实现面包板上的挑战与调试4.1 元件布局、供电与接地艺术在面包板上搭建多级运放电路布局是成功的一半。首要原则是电源去耦在每个运放的电源引脚V和V-附近尽可能靠近引脚的地方跨接一个0.1μF的陶瓷电容到地。这为运放的高速电流需求提供了低阻抗通路能有效抑制通过电源线耦合的高频噪声和防止自激振荡。对于整个电路板还应在大电容如100μF电解电容处为电源总线滤波。接地策略至关重要。混乱的地线会引入严重的噪声。应采用星型接地或单点接地准备一块独立的、干净的“信号地”区域可以用一根粗导线或铜箔所有运放的接地端、滤波电容的接地端、输入信号的参考地都集中连接到这一点。然后用另一根单独的导线将这一点连接到电源的接地端。绝对避免形成地线环路。开始搭建时不要一次性连接整个系统。应分模块搭建、分模块测试。先搭建INA用函数发生器输入一个差分小信号如1mV 1kHz用示波器双通道分别测量输入和输出计算实际增益并与仿真对比。注意此时示波器探头的接地夹必须接在电路的信号地上。4.2 INA模块的实测与高误差排查按照仿真电路搭建INA使用±9V双电源供电。实测中很可能发现增益远高于设计值如本项目的1600倍误差高达60%。这通常由以下原因导致电阻值误差面包板常用的碳膜电阻公差多为5%甚至更差。一个5%误差的电阻在增益公式中可能被放大。用万用表逐一测量并记录实际使用的每个电阻值代入公式重新计算理论增益你会发现问题所在。741运放的输入失调电压通用741的输入失调电压典型值为2mV最大可达6mV。这个电压会被放大1000倍在输出端产生2-6V的直流偏移这很可能使输出饱和顶到电源轨你测量到的“增益”实际上是失真的饱和波形。解决方法是在运放的同相端和反相端之间接入一个由电位器构成的可调调零电路或者在后续级加入隔直电容但需注意心电信号含有低频成分。共模信号处理不当如果输入信号的共模部分如两个输入端对地的平均电压不在运放的共模输入电压范围内对于±9V供电的741大约在-7V到7V运放将无法正常工作。确保你的函数发生器或电极输入共模电平在合理范围内。避坑指南在测试高增益放大器时先将增益调低例如临时减小R4或增大R3确保电路基本功能正常、输出无饱和再逐步调整电阻至目标增益值进行测试。同时始终用示波器监视输出波形确保是干净的正弦波放大而不是畸变或振荡。4.3 滤波器模块的实测与频率响应校准搭建陷波和低通滤波器。测试滤波器最有效的方法是使用函数发生器的频率扫描功能配合示波器或万用表交流电压档测量输出幅度。对于陷波滤波器将函数发生器设置为正弦波输出一个恒定幅度如1Vpp的信号。从10Hz缓慢扫频至200Hz观察输出幅度。当扫到60Hz附近时输出应出现一个明显的凹陷。用示波器可以精确找到幅度下降-3dB的两个频率点f1和f2中心频率f0 √(f1*f2)-3dB带宽BW f2 - f1。实测中心频率可能与60Hz有偏差这主要是由电容和电阻的容差造成。可以通过并联或串联小电阻、小电容进行微调。例如若中心频率偏高可在主电容上并联一个小电容如几皮法到几十皮法。对于低通滤波器从10Hz扫频至1kHz。在150Hz附近输出幅度应下降至输入幅度的70.7%-3dB。同样截止频率可能偏移。赛伦-凯滤波器的截止频率对RC乘积敏感微调其中一个电阻通常能校准。一个关键测试单独测试滤波器时输入一个纯净正弦波。但当接入整个系统后由于前级INA可能已经引入噪声或振荡会导致滤波器性能下降。因此每级测试时最好用函数发生器模拟前级的理想输出进行测试。4.4 系统集成与生物信号采集将三级电路按顺序连接INA的输出接陷波滤波器的输入陷波滤波器的输出接低通滤波器的输入。集成时最大的挑战是级间阻抗匹配和噪声累加。理想情况下前一级的输出阻抗应远小于后一级的输入阻抗以避免信号衰减。运放电压跟随器缓冲器输出阻抗极低是级联的理想选择。在我们的电路中INA和两个滤波器本身都基于运放输出阻抗较低直接级联一般问题不大。但最好在每级之间预留测试点方便隔离排查问题。使用函数发生器模拟一个“干净”的心电信号可以是一个幅值1mV、频率1.2Hz的类正弦波叠加一个更高频的QRS脉冲输入集成电路观察最终输出。此时你可能会发现波形底部有高频毛刺或低频漂移。高频毛刺可能来自电源噪声或布局不当可尝试加强电源去耦。低频漂移基线漂移可能来自电极极化电位、呼吸运动或运放的失调电压温漂在硬件层面较难根除通常需要在软件中进行数字高通滤波处理。最后进行真人ECG采集。使用一次性心电电极片按照导联I的放置方法正极LA置于左踝负极RA置于右腕右踝RL作为驱动右腿电路的参考地本项目简化为接信号地。将LA和RA分别接入INA的两个差分输入端RL接入电路板信号地。务必注意安全确保电路与市电完全隔离使用电池或隔离电源供电所有裸露导线绝缘良好。上电后在示波器上应能看到周期性的心电波形。由于电路相位反转你看到的可能是倒置的波形这可以通过交换INA的输入线或后期软件处理来纠正。5. 基于Arduino的数据采集与心率计算5.1 硬件连接与模拟输入配置Arduino Uno的ADC模数转换器输入范围是0-5V而我们的电路输出是±V的双极性信号如±5V。直接连接会导致负电压部分被削波。因此需要在信号送入Arduino前增加一个电平移位电路将双极性信号抬升到0-5V的单极性范围。一个简单的方案是使用一个运放构成的加法器将ECG信号V_ecg与一个2.5V的直流偏置V_ref相加。当V_ecg -2.5V时输出为0V当V_ecg 0V时输出为2.5V当V_ecg 2.5V时输出为5V。这样ECG波形就以2.5V为中心在0-5V内变化。2.5V参考电压可以由Arduino的5V电源通过两个等值电阻分压获得并经过一个电压跟随器缓冲后使用。将移位后的信号连接到Arduino的A0引脚。在代码中配置ADC参考电压为默认的5V。采样率是关键参数。根据奈奎斯特采样定理要无失真地恢复150Hz的信号采样率至少需要300Hz。考虑到心电波形特征尤其是陡峭的R波通常需要更高的采样率如250-500Hz。Arduino Uno的ADC在保持一定精度下最高采样率约为10kHz完全满足要求。我们可以使用analogRead()函数并通过控制循环速度或使用定时器中断来稳定采样率。5.2 心率检测算法与代码实现心率检测的核心是识别QRS波群中的R波峰值。一个在嵌入式系统上简单有效的算法是阈值检测法。数据缓冲与预处理连续采样一定数量的点如对应2-3秒的数据存入数组。首先进行软件数字滤波以进一步降噪。可以实现一个移动平均滤波器低通来平滑高频噪声以及一个一阶高通滤波器如减去移动平均来消除基线漂移。动态阈值计算心电信号幅度会因人、因电极位置而异。需要动态更新阈值。可以跟踪信号的最大值峰值和最小值谷值。将阈值设置为阈值 零点 灵敏度 * (峰值 - 零点)其中“零点”可以是信号的平均值或中值灵敏度经验值可取0.5-0.7。R波检测遍历滤波后的数据。当发现一个数据点高于阈值且其前后点呈上升-下降趋势时标记为一个R波候选。为了防止对同一个R波重复检测设置一个不应期如对应200-300ms的时间窗口在检测到一个R波后在不应期内忽略新的峰值。心率计算记录连续两个R波之间的时间间隔RR间期单位为毫秒。心率HR 60000 / RR间期(ms)。为了显示稳定可以计算最近几个RR间期的平均值。// 简化伪代码逻辑 const int sampleRate 500; // 500 Hz const int refractoryPeriod sampleRate * 0.25; // 250ms 不应期 float threshold 512.0; // 初始阈值对应2.5V int lastPeakIndex -refractoryPeriod; // 上一个R波位置 long lastRTime 0; float bpm 0; void loop() { int sensorValue analogRead(A0); // 1. 软件滤波 (此处省略具体滤波代码) float filteredValue softwareFilter(sensorValue); // 2. 更新动态阈值 (基于近期数据统计) updateThreshold(filteredValue); // 3. 检测R波 int currentIndex getCurrentSampleIndex(); if (currentIndex - lastPeakIndex refractoryPeriod) { if (filteredValue threshold isPeak(filteredValue)) { // 发现R波 lastPeakIndex currentIndex; long currentTime millis(); if (lastRTime 0) { long rrInterval currentTime - lastRTime; bpm 60000.0 / rrInterval; // 计算瞬时心率 // 可对bpm进行滑动平均滤波 } lastRTime currentTime; } } // 4. 发送数据到串口绘图或显示 Serial.print(filteredValue); Serial.print( ); Serial.println(threshold); // 可选同时输出阈值线 }5.3 串口可视化与常见问题调试利用Arduino IDE的串口绘图器是一个极佳的调试工具。将滤波后的信号值和阈值线实时打印到串口可以直观看到波形和检测过程。常见问题与排查波形倒置如前所述这是电路相位反转造成的。在Arduino代码中可以用1023 - sensorValue来反转或者直接在心电电极连接时交换正负极。基线剧烈漂移软件高通滤波器的截止频率设置过低如0.5Hz可能导致恢复缓慢设置过高如1Hz可能损伤ST段。需要折中调整。硬件上检查电平移位电路的参考电压是否稳定。误检或漏检R波误检通常因噪声尖峰超过阈值引起。加强软件滤波如增加移动平均窗口或提高阈值灵敏度系数。漏检R波幅度太小。检查电极接触是否良好皮肤用酒精清洁涂抹导电膏或降低阈值灵敏度系数。也可能是高通滤波器滤掉了过多的低频成分。心率值跳动剧烈瞬时心率对单个RR间期非常敏感。实现一个长度为4-8的滑动窗口对连续计算出的心率值或RR间期进行平均能有效稳定显示结果。串口绘图卡顿或数据丢失确保串口波特率设置正确如115200并且循环中没有不必要的延迟。采样率不要超过串口传输的承载能力。如果采样率很高可以考虑在Arduino端进行大部分计算只将心率和关键波形点发送给电脑。完成这些步骤后你就能在电脑屏幕上看到实时跳动的心电波形和稳定的心率数值。这个从传感器、模拟电路到嵌入式软件的全链路实践完整地再现了一个简易心电图机的核心功能。尽管它离医疗级设备还有巨大差距但其中涉及的信号放大、噪声抑制、滤波设计和数字处理思想是生物医学电子领域永恒的基础。