生物医学信号分析器 返回
生物医学工程

生物医学信号分析器(心电图/脑电图)

调节心率、噪声水平和低通滤波截止频率,实时可视化ECG波形、RR间期和频谱分析。

ECG波形主要成分
P波:心房除极(0.1s)
QRS综合波:心室除极(0.08s)
T波:心室复极(0.2s)

心率:HR = 60 / RR[s]
心率 (HR)
72 BPM
RR 间期
833 ms
QRS 宽度
80 ms
信噪比(估算)

什么是生物医学信号分析器

🧑‍🎓
心电图里那些上下起伏的波浪线,到底代表什么意思啊?
🎓
简单来说,那条波浪线就是你的心脏在“发电”的轨迹。每一次心跳,心脏肌肉都会产生微弱的电流,传到皮肤表面被电极捕捉到。你看模拟器上那个典型的波形,它就像一个字母组合:P-QRS-T。P波是心房收缩,那个又高又尖的QRS波是心室强力收缩,最后的T波是心室放松恢复。你可以试着拖动“心率”滑块,从40调到200,看看波形是怎么随着心跳加快而“挤在一起”的。
🧑‍🎓
诶,真的吗?那为什么我们测到的心电图总感觉毛茸茸的,不像模拟器里这么光滑?
🎓
问得好!在实际测量中,信号里会混入各种“噪声”。比如你呼吸造成的基线上下漂移,肌肉颤抖产生的高频毛刺,还有无处不在的50Hz工频干扰。你可以在模拟器里把“噪声水平”滑块拉高,立刻就能看到干净信号变成一团糟。这时候就需要“滤波”来救场了,试着调节“滤波截止频率”滑块,看看高频噪声是怎么被一点点滤掉的,只留下我们关心的心电波形。
🧑‍🎓
原来噪声这么麻烦!那下面那个显示“RR间期”和“频谱”的图又是干嘛用的?
🎓
那可是诊断的宝藏!RR间期就是两个R波峰之间的时间,直接决定了心率。但更重要的是它的变化,叫做心率变异性(HRV),能反映你神经系统的紧张程度。频谱图则是把波形从“时间域”转换到“频率域”,让你一眼看出信号里哪些频率成分是主导的。比如,一个健康人的心电频谱在低频段(对应HRV)会有一定能量。你同时调节心率和滤波,观察频谱图上的峰值是怎么移动和变化的,就能直观理解信号处理的核心了。

物理模型与关键公式

心电图信号可以建模为一系列特征波形(P, QRS, T)的合成,并叠加噪声。一个常用的简化数学模型是使用高斯函数及其导数来模拟这些波形成分。

$$ECG(t) = \sum_{i}a_i \cdot \exp\left(-\frac{(t - p_i)^2}{2b_i^2}\right) + \sum_{j}c_j \cdot \frac{d}{dt}\left[\exp\left(-\frac{(t - q_j)^2}{2d_j^2}\right)\right] + Noise(t)$$

其中,$a_i, b_i, p_i$ 分别控制第i个高斯波形的幅度、宽度和中心位置(用于模拟P波和T波);$c_j, d_j, q_j$ 控制其导数波形的参数(用于模拟陡峭的QRS波);$Noise(t)$ 代表叠加的各类噪声。

心率(HR)与RR间期(两个R波峰值之间的时间)直接相关,是信号分析中最基本的计算。

$$HR \, (bpm) = \frac{60}{RR \, (seconds)}$$

$HR$ 是心率,单位是每分钟心跳次数(bpm)。$RR$ 是RR间期,单位是秒(s)。例如,RR间期为1.0秒时,心率为60 bpm。

现实世界中的应用

心律失常诊断:这是心电图最经典的应用。通过分析P波缺失、QRS波形态异常或RR间期不规则,医生可以诊断房颤、室性早搏等多种心律失常。模拟器中改变心率,观察RR间期是否规律,就是在模拟这一诊断基础。

睡眠与压力研究:心率变异性(HRV)分析广泛应用于评估自主神经功能。HRV降低通常与压力过大、睡眠质量差或某些疾病相关。研究通过长时间监测ECG,分析RR间期的微小波动频谱来实现。

运动负荷试验:在跑步机或动感单车上进行运动时连续监测心电图,观察心率上升是否与运动负荷匹配,以及ST段是否出现压低或抬高,这是诊断冠心病(如心绞痛)的重要手段。

可穿戴设备健康监测:智能手表和手环通过光电或电极持续采集类似ECG的信号,实现心率、心律不齐(如房颤)的初步筛查。这些设备面临更强的噪声干扰,其信号滤波算法正是本模拟器所演示技术的实际应用。

常见误解与注意事项

首先,需要注意“认为完全去除噪声就能达到完美”的常见倾向。在模拟器中过度加强“滤波器截止频率”确实能使波形变得平滑。然而在实际临床场景中,波形的尖锐拐角(高频成分)可能包含重要信息。例如,室性早搏的波形会变宽,但若使用过强的滤波器使其过度圆滑,则存在漏诊此特征的风险。请记住:滤波是一项“权衡取舍”的工作。

其次,关于“心率恒定”的前提假设。模拟器的心率滑块可设置固定值,但真实人体的心率始终存在波动(心率变异性:HRV)。例如即使在静息状态下,RR间期也会以0.85秒、0.92秒、0.88秒…的幅度微妙起伏。这种波动本身即是自主神经功能的指标。分析时不仅需关注短期平均心率,更应注意这种波动幅度。

最后,“将所有显著波峰都视为QRS波”的误判。实际心电图中常会出现表现为大幅伪影(伪信号)的波峰。例如电极脱落导致的基线大幅漂移,或强烈的肌电噪声。将模拟器的噪声水平调至最高时,可重现QRS波上叠加其他波峰的现象。通过算法自动检测QRS波时,需要建立能综合振幅、波形形态及前后文脉进行判断的“逻辑”。

相关工程领域

本模拟器核心的信号处理技术与振动·声学工程完全一脉相承。观察ECG频谱的行为,本质上与从机械振动数据中识别异常声响的“FFT(快速傅里叶变换)分析”完全相同。例如旋转机械的轴不平衡会在特定频率形成峰值,这与识别心电图中电源噪声(50/60Hz峰值)的方法具有共通性。

同时,控制工程的滤波器设计理论可直接应用。生物信号使用的带通滤波器由低通与高通滤波器组合而成,这与反馈控制系统中去除无用频带信号的“陷波滤波器”或“低通滤波器”设计思想完全一致。在模拟器中调整截止频率的感知,直接关联着调整控制系统带宽的感知。

此外还与材料·结构工程的无损检测领域相通。超声波探伤检测中,材料内部缺陷的反射波(信号)常叠加多种噪声。将“真实信号”与“噪声”分离并估算缺陷位置与尺寸的信号处理技术,可以说与检测淹没于噪声中的微小P波、T波的心电图解析技术属于同源分支。

进阶学习指引

第一步可转向建立“作为时间序列数据的生物信号”意识的解析。模拟器主要观察单次心搏,而实际需处理长时间数据。例如读取5分钟ECG数据作为数组,检测所有R波位置生成RR间期时间序列,进而计算其标准差(SDNN:心率变异性指标之一)。在编程环境(如Python的NumPy/SciPy)中实际编写代码最为有效。

数学基础方面,线性代数与统计学知识至关重要,因为它们是现代信号处理与机器学习解析的根基。例如主成分分析(PCA)可用于从多维生物信号(如多导联心电图)中分离噪声或提取特征,其理解必须建立在协方差矩阵与特征值·特征向量的知识体系上。

推荐后续具体研习“小波变换”。与FFT提供信号整体频率成分不同,小波变换能同时提供“时间与频率的双重信息”。这将成为解析非平稳信号——如突发心律失常、随睡眠阶段变化的脑电波(EEG)——的强有力工具。建议在掌握模拟器中的傅里叶变换后,将其作为进阶方向接触学习。