心电图的基础知识

一、基础概念1. 心电向量心脏的整体的综合心电向量是由每一个心肌细胞矢量叠加而成的而且心脏电激动不是单点一次性触发从窦房结起步依次激动心房、房室结、左右心室激动区域随时间推移不断移动综合向量的大小、朝向也持续改变。空间立体心电向量环就是心脏每一瞬间算出 1 个综合心电向量有大小、三维朝向把所有时刻向量的端点按时间依次衔接在空中画出一条立体闭合轨迹而心电图就是这条立体轨迹在某一条观测轴线上的一维投影曲线多导联和向量环的关系12 个导联相当于从 12 个不同视线方向分别截取立体向量环的一维投影反向推算还原完整三维向量环形态定位病变。心电向量图直接画出平面上向量端点移动轨迹环本身横轴纵轴都是空间方位常规心电图只取单条轴线投影、按时间展开横轴固定为时间。2. 心电信号在心电图上的表示(ElectrocardiographyECG或者EKG)心肌在每次收缩之前先必产生电激动电流传布全身各处产生不同的电位且强弱方向不断变化通过心电图机把这种变动着的电位连续描记成的曲线这个反映心电活动的曲线就是心电图。心脏心肌激动产生的是空间立体心电向量环三维的电流变化但心电图机只能记录二维平面电位变化。每一组导联本质是设定一组「观测视角」把立体的心电向量投影到这个视角对应的平面上输出波形多导联是从胸腔不同方位布设观测电极分别采集心脏电活动在对应方向上的投影信号所有导联叠加才能还原整颗心脏立体的电激动全貌没有任意一个导联能代表整颗心脏全部位置的信号。心脏每个瞬间只有1 个三维综合心电向量12 导联对应 12 条朝向不同的观测轴线把同一个向量分别投影到各轴线就能独立生成 12 条时序心电曲线最终并排打印在同一张心电图图纸上共用统一时间轴。3. 心电图波形拆解这是标准单导联正常心动周期心电图波形解析图采用通用心电图纸网格标尺横向计时、纵向计电压完整展示一次心跳从心房激动到心室激动、复极全流程的波形分段如下图除极心肌细胞静息时内外存在电位差内负外正接收到电刺激后离子跨膜流动、电位反转内正外负这个电位翻转过程就是除极。复极复极是心肌细胞完成除极、触发收缩之后把细胞膜电位从「内正外负」的兴奋状态重新变回静息时「内负外正」初始电位的电学复位过程属于和除极反向的离子流动变化是下一次心电激动能正常触发的前提。P 波P wave心房除极波代表左右心房心肌除极收缩前电激动产生的电位变化是单次心动周期最早出现的波形。QRS 复合波QRS complex全称心室除极复合波由 Q、R、S 三个子波组成Q 波R 波前向下的负向小波R 波整个波形最高的正向主峰S 波R 波之后向下回落的负向波整体记录心室整体除极的电活动波形陡峭、振幅最大。心房复极波形被QRS波完全覆盖所有没有提及。T 波T wave全称心室复极波代表心室完成收缩后、心肌恢复静息电位复极的电位变化波形宽缓圆润。U 波U wave全称后继复极小波出现在 T 波之后一般和心室浦肯野纤维复极相关振幅很低。PR 间期PR interval从 P 波起点到 QRS 波起点反映心房开始激动、信号经房室结传导抵达心室的总耗时房室传导时间。QRS 复合波间期QRS complex durationQ 波起点至 S 波终点衡量心室整体除极耗费时长。JT 间期JT intervalJ 点开始到 T 波终点聚焦心室除极结束后、全程复极的时长。QT 间期QT intervalQ 波起点至 T 波终点覆盖心室从除极启动到复极完成的完整周期临床评估心室整体电活动时长核心指标。节段两个波段之间的平直基线段不含首尾波主体PR 段PR segmentP 波终点到 QRS 波起点对应电信号在房室结、希氏束等传导通路传递的阶段贴近基线。ST 段ST segmentJ 点到 T 波起点心室除极完毕、复极刚开始的平稳阶段偏移抬高 / 压低是判断心肌缺血、心梗关键依据。J 点J pointQRS 波终点和 ST 段起始的衔接拐点是测量 ST 段偏移量的基准点位。所有导联单次心动周期都完整包含这套电活动时序对于朝向相反的导联只不过是将波形换个方向。二、ECG数据集「ECG信号处理——1认识MIT-BIH数据库」2024年10月25日「ECG信号处理——2MIT-BIH数据库在线分析ATM」2024年11月4日「ECG信号处理——3WFDB工具包与MIT-BIH数据库读取」2024年11月13日从下面这个官方链接页面可以下载到所有48组MIT-BIH心电数据https://www.physionet.org/content/mitdb/1.0.0/数据介绍来源MIT‑BIH Arrhythmia Database Directory.hea头文件ASCII 文本格式记录整条心电记录的基础元信息采样率、导联类型、通道数量、采样总量、信号增益、基线偏移、数据存储格式是解析.dat的索引依据可用记事本直接打开查看。.dat原始心电波形数据二进制格式存储双通道连续心电电压采样值采用 MIT 自定义紧凑二进制编码3 字节存 2 个 12 位采样是绘制 ECG 波形、做信号预处理的原始数据源。.atr专家标注注释文件二进制格式心脏病专家人工标注结果包含每个 R 波峰值坐标、心搏类别正常 / 室早 / 房早等、心律类型、噪声标记用来训练心律失常分类模型、校验算法识别精度。辅助拓展文件.xwsPhysioNet 配套波形查看工具xwaves的工程配置文件保存该条记录的视图缩放、导联显示、标注筛选等界面参数供软件可视化快速调用。.at_标注类临时备份 / 残缺缓存文件属于下载、解压过程产生的不完整注释副本非标准业务文件正式解析时无需使用。1. 数据筛选标准数据库共收录48 条有效记录分为两组第一组共 23 条编号 100~124存在缺号采用纯随机抽样选取每条记录时长均略超过 30 分钟。本组用于模拟临床常规场景下心律失常检测算法会遇到的各类正常波形与噪声伪影抽样采用随机数法筛选磁带及 30 分钟分段仅排除两条心电通道信号质量均不满足人工专家判读标准的片段。第二组共 25 条编号 200~234存在缺号定向筛选收录大量罕见但临床意义重大的心电现象这类病例无法通过纯随机抽样充分覆盖。本组包含复杂室性、交界性、室上性心律失常及各类传导阻滞多条记录的节律形态、QRS 波多变性、信号噪声会对自动心律失常检测器形成极高识别难度是本数据库在算法评测领域最具代表性的测试样本。2. 心电图导联配置绝大多数记录采用双通道采集上方通道改良肢体 Ⅱ 导联MLII电极贴于胸壁正常心搏在该通道 QRS 波清晰突出是算法识别首选通道。下方通道多为改良 V1 导联少数为 V2/V51 例为 V4电极同样贴于胸壁。该导联组合为贝斯以色列心律失常实验室标准采集方案。下方通道导联轴线与心脏平均电轴近乎正交正常窦性心搏波形多呈双向、电位接近等电位线难以分辨而异位搏动在该通道波形特征更突出以记录 106 为例。特例说明记录 114 双通道信号上下倒置该情况临床偶发心律失常检测算法需兼容此种通道颠倒场景。记录 102、104 因患者胸区包扎手术敷料无法采集改良 Ⅱ 导联因此上方通道改用改良 V5 导联。3.模拟磁带录制与回放设备心电信号的采集过程人体心电模拟电压信号 → 便携记录仪采集 → 写入模拟磁带外出随身记录阶段 → 实验室回放磁带 → ADC采样量化生成数字心电MIT-BIH最终数据库文件磁带信号固有缺陷通道时序偏移Skew部分记录仪双通道信号最大时序偏移可达 40ms。偏移分为两类一是录放设备磁头微小角度差带来的固定偏移二是磁带录放时纵向微小抖动造成的可变偏移偏移幅度与固定偏移处于同一量级。该问题同样存在于 AHA 心电数据库会对双通道实时分析算法造成干扰。磁带速度扰动磁带轻微快慢波动对常规心律失常检测器无明显影响低速回放极易出现磁带粘连、打滑打滑片段均已添加注释标记。心率变异性研究需重点考量磁带抖晃失真本批原始磁带无抖晃补偿机制。机械结构引入固定频率噪声伪影下表为录放设备机械部件对应的干扰频率实时基准| 频率 (Hz)| 噪声来源 || ---- | ---- ||0.042 | 记录仪压带轮 ||0.083|2 倍速回放时的设备主轴 ||0.090 | 记录仪主轴 ||0.167 | 实时回放时的设备主轴 ||0.10~0.18 | 收带盘频率随磁带卷径增大持续降低||0.20~0.36 | 供带盘频率随磁带卷径减小持续升高|干扰强度排序实时回放 0.167Hz 噪声最强其次为 0.090Hz 噪声、2 倍速回放 0.083Hz 噪声0.042Hz 噪声幅值最低。其余传动结构衍生频率0.42Hz、1.96Hz、9.1Hz、42Hz无明显干扰。供带盘噪声显著强于收带盘额外在 0.100.18Hz、0.300.54Hz 产生干扰。起搏信号记录说明48 条记录中 102、104、107、217 包含起搏心搏。原始模拟磁带对起搏脉冲信号保真度不足无法仅依靠脉冲幅值 / 上升沿斜率 / 脉宽完成实时识别但数字化数据库完整还原磁带信号可适配基于磁带离线分析的起搏伪影检测算法。4.数字化采集流程回放设备模拟输出经0.1~100Hz 带通滤波抗混叠、防止模数转换器饱和带宽覆盖心电全部有效频率成分。采用麻省理工生物医学工程中心与贝斯以色列医院自研硬件进行模数转换单通道采样率 360Hz该采样频率便于算法内置 60Hz 工频陷波滤波器。记录仪为电池供电数据库中 60Hz 工频噪声大多产生于回放环节2 倍速数字化的记录工频噪声等效为 30Hz 及其倍频。双通道采样时刻近乎同步通道间采样时序偏差仅微秒级远小于前文磁带带来的数十毫秒级偏移。模数转换器为单极性、11 位量化电压量程 ±5mV采样值范围 0~2047数值 1024 对应 0 电位基线。原始数据受存储容量限制采用 8 位一阶差分编码存储结合采样率与量化精度系统最大可记录信号变化率 ±225mV/s。实际场景下仅少数强噪声片段偶尔超出该上限对整体信号质量无实质影响。本 CD-ROM 版本已将差分数据还原为原始幅值采用 3 字节打包存储两组 12 位采样值存储格式详见 signal (5) 文档。每个采样点以24 位二进制补码形式存储幅值字节排布采用小端序低有效字节优先存放。5.心电标注体系标注生成流程采用简单斜率型 QRS 检测器生成初始标注全部识别事件标记为正常心搏每条 30 分钟记录打印两份完整心电图纸标注初稿打印在图纸边缘两份图纸分送两名心脏病专家独立判读补充漏检心搏、删除误检标记、修正全部异常心搏分类同步添加心律标注、信号质量标注与文字备注将两份专家标注转录至数字文件逐心搏自动比对生成标注差异图纸两名专家共同复核所有分歧点协商统一标注修正内容二次录入后审计程序校验标注逻辑一致性并提取每条记录最长、最短 RR 间期排查漏检 / 误检。标注位置校准数据库早期版本心搏标记置于 R 波附近但人工标注无法精准对齐波峰1983 年后全部记录117 除外执行统一校准流程对上方心电通道做数字带通滤波强化 QRS 波修正滤波器相位偏移后将所有心搏标记对齐至 QRS 局部极值点少量强噪声心搏人工微调位置。记录 117 的标注直至 1998 年 3 月完成对齐。校准后所有标注基准点统一落在 R 波波峰精度满足波形平均、心率变异性研究需求但前文磁带抖晃缺陷仍需纳入分析考量。数据库总计约 109000 条心搏标注1980 年数据库正式发布后持续修正标注错误关键修正节点发布首 7 年修正 16 处心搏标注记录 104/108/114/203/207/217/222记录 214 原有左束支传导阻滞心搏全部误标为正常心搏心律标注大幅更新新增二联律、三联律、起搏心律标记。1998 年 10 月修正记录 203 心律标注2001 年 10 月发现并修正第 17 处心搏标注错误记录 2092003 年 4 月人工微调记录 119 中 26 个室性早搏标注位置最大偏移 74ms2003 年 5 月修正第 18 处心搏标注错误记录 2142005 年 4 月记录 222 中原房颤片段部分 / 全部修正为房扑标注2008 年 4 月修正 3 处心搏标注记录 108 两处、记录 215 一处2010 年 6 月发现并修正第 22、23 处心搏标注错误均在记录 203。感谢 Bob Bruce、Pat Hamilton、尹登峰、Roger Mark、Sebastian Vasquez、Mariano Llamedo Soria 反馈标注漏洞。6.波形图符号释义一心搏标记标准层符号 全称与释义· / N 正常窦性心搏L 左束支传导阻滞心搏R 右束支传导阻滞心搏A 房性早搏a 室内差异性传导房早J 房室交界性早搏S 室上性早搏V 室性早搏PVCF 室性 - 正常心搏融合波[ 心室扑动 / 颤动起始! 心室扑动波] 心室扑动 / 颤动终止e 房性逸搏j 房室交界性逸搏E 室性逸搏/ 起搏心搏f 起搏 - 正常心搏融合波x 未下传 P 波阻滞型房早Q 无法分类心搏| 孤立类 QRS 波形伪影二心律标注绘制于心搏标记下方标记 释义(AB 房性二联律(AFIB 心房颤动(AFL 心房扑动(B 室性二联律(BII 二度房室传导阻滞(IVR 自主室性心律(N 正常窦性心律(NOD 房室交界性心律(P 起搏心律(PREX 预激综合征WPW(SBR 窦性心动过缓(SVTA 室上性快速心律失常(T 室性三联律(VFL 心室扑动(VT 室性心动过速三信号质量 / 备注标注绘制于心搏标记上方标记 释义qq 通道信号质量变更首字符 (c/n) 代表上方通道干净 / 噪声第二位代表下方通道U 双通道重度噪声 / 信号丢失心电图无法判读M / MISSB 漏搏P / PSE 长间歇T / TS 磁带打滑失真专业术语注释便于算法 / 心电研究查阅Holter动态心电记录仪MLIIModified Limb Lead II改良胸贴式肢体 Ⅱ 导联QRS complex心室除极复合波PVCPremature Ventricular Contraction室性早搏Skew双通道时序偏移Anti-aliasing filter抗混叠滤波器Notch filter工频陷波滤波器A/D converter (ADC)模数转换器Fusion beat融合心搏Bundle branch block束支传导阻滞Atrial flutter/fibrillation房扑 / 房颤Ventricular tachycardia/flutter室速 / 室扑Baseline isoelectric line等电位基线Annotator心电标注专家Bandpass filter带通滤波器。三、ECG信号中常见噪声与处理方法ECGSYN一个纯净心电信号生成器ECGSYN 是一个用于生成合成心电图ECG信号的 MATLAB 函数。它通过模拟心电图信号的动态系统来生成具有真实特征的合成 ECG 信号。下载地址ECGSYN - A realistic ECG waveform generator v1.0.01.去噪效果评价指标SNRSignal-to-Noise Ratio信噪比是衡量信号质量的常用指标它表示信号的强度与噪声的强度之比。在去噪后信号的信噪比越高说明去噪效果越好。其中xi是原始信号x ^ i \hat{x}_ix^i​是去噪后的信号N是信号的长度。MSEMean Squared Error均方误差是衡量信号重建误差的指标它表示原始信号与去噪后信号之间的平均平方误差。MSE 越小说明去噪效果越好。其中xi是原始信号x ^ i \hat{x}_ix^i​是去噪后的信号N是信号的长度。RMSERoot Mean Squared Error根均方误差是均方误差的平方根它与均方误差类似也是衡量信号重建误差的指标。RMSE 越小说明去噪效果越好。其中xi是原始信号x ^ i \hat{x}_ix^i​是去噪后的信号N是信号的长度。Peak-to-Peak Amplitude峰-峰值是衡量信号幅度变化的指标它表示去噪后信号的最大值与最小值之差。在去噪过程中应避免信号的过度衰减保持信号的幅度特征。其中x ^ i \hat{x}_ix^i​是去噪后的信号。2.常见噪声及去噪算法1基线漂移及其去噪基线漂移介绍https://blog.csdn.net/qq_42936724/article/details/144987265?spm1001.2014.3001.5501一、产生原理按诱因分类基线漂移是叠加在心电波形上的极低频慢变干扰会让整体信号基准线缓慢上下偏移。效果如下图人体生理活动最主要来源呼吸运动吸气 / 呼气时胸腔扩张收缩肺部组织电阻率、心脏空间朝向周期性改变体表采集到的心电综合向量随之波动频率稳定在0.1~0.3Hz对应成人正常呼吸频次漂移波形接近正弦周期波动。体位 / 肢体微动翻身、抬手、肌肉轻微牵拉直接改变电极按压位置、皮肤形变造成电极 - 皮肤接触阻抗突变产生非平稳的慢趋势漂移出汗会软化皮肤角质层进一步加剧阻抗波动。电极接触层面问题电极胶老化、粘贴松动、皮肤角质过厚接触电阻随外力、湿度动态变化等效引入缓慢直流偏移长期监测时该问题在动态 Holter 记录中尤为突出MIT-BIH 数据集原始磁带记录也会叠加这类漂移。采集仪器硬件误差前置放大电路存在固有直流偏置工作时元器件温漂、供电电压小幅波动会叠加缓慢的基线偏移早年模拟磁带录放设备的转速抖动也会衍生低频基线畸变。其他次要因素体表局部体温变化、皮下组织代谢带来的导电性微变形成长时程缓慢漂移趋势。二、基线漂移的信号特点与特征低频特性基线漂移主要表现为 0.05~1 Hz 的低频波动与心电信号的高频特征如 PQRST 波形成鲜明对比。缓慢变化基线漂移的变化速度相对较慢通常在数秒到数分钟内完成一个周期的变化。幅度变化基线漂移的幅度通常较小但足以掩盖心电信号的某些特征波影响心电图的准确分析和诊断。非线性特性基线漂移的产生原因复杂多样其变化过程具有一定的非线性特性这增加了去除基线漂移的难度。方法一高通滤波基于基线漂移低频、变化缓慢、幅度变化小的特征使用高通滤波器进行去噪。低频特性基线漂移主要表现为 0.05~1 Hz 的低频波动。高通滤波器可以有效地阻断低于其截止频率的低频成分从而去除基线漂移。通常高通滤波器的截止频率设置在 0.5 Hz 到 1 Hz 之间以确保心电信号的高频特征波如 PQRST 波不受影响。缓慢变化由于基线漂移的变化速度相对较慢高通滤波器能够将这种缓慢变化的低频成分滤除保留心电信号的快速变化特征。幅度变化虽然基线漂移的幅度较小但高通滤波器能够有效地去除这种幅度变化避免其掩盖心电信号的特征波提高心电图的准确分析和诊断能力。适用场景适用于基线漂移频率较低且与心电信号特征波频率有明显区别的场景。方法二小波变换https://blog.csdn.net/qq_42936724/article/details/145248080?spm1001.2014.3001.5501小波变换将信号分解为逼近信号和细节信号逼近信号反映了信号的整体趋势和主要特征而细节信号则包含了信号的局部变化和噪声成分。入信号和小波函数做卷积将输入信号分解为一系列的信号。如下图小波分解基本思想小波分解的基本过程是通过对信号进行滤波将信号分解成不同频带的成分。每一次分解都包括两个步骤1、低频部分逼近信息通过尺度函数滤波得到。2、高频部分细节信息通过小波函数滤波得到。在每一层分解中低频部分会继续进行下一次分解而高频部分被提取为细节信息。小波分解数学解释设原始信号为 在进行第 层的小波分解时信号首先被分解为逼近信息低频部分和细节信息高频部分。这两个部分通过滤波器对信号进行处理公式如下第j jj层逼近信息低频部分其中h ( n ) h(n)h(n)是尺度函数低通滤波器的系数用于捕捉信号的低频部分平滑信息。尺度函数的作用是对信号进行平滑处理提取信号的主结构。Cj1,n是上一层的逼近系数。第j jj层细节信息高频部分其中g ( n ) g(n)g(n)是小波函数高通滤波器的系数用于捕捉信号的高频部分细节信息。小波函数的作用是捕捉信号的变化或细节尤其是快速变化的部分。Cj1,n仍然是上一层的逼近系数。多层小波分解思想在多层小波分解中分解过程是递归的。第0层原始信号f ( x ) f(x)f(x)通过尺度函数系数g ( n ) g(n)g(n)和小波函数系数h ( n ) h(n)h(n)滤波得到第 1 层的逼近信息和细节信息。第1层从第0层得到的逼近信息低频部分继续进行滤波得到第 2 层的逼近信息和细节信息。以此类推直到分解达到指定的层数N。注意分解的层数 决定了信号被分解为多少个尺度的逼近信息和细节信息。较高的层数可以提供更精细的信号表示但也会增加计算复杂度。每分解一层信号的时间分辨率减半。因此分解的层数N通常受信号长度限制。小波基的选择与特点小波基的选择对于小波变换的效果影响很大。不同的小波基有不同的数学结构和性质因此适用于不同类型的信号处理任务。小波基主要由尺度函数low-pass filter和小波函数high-pass filter组成。常见的小波基包括哈尔小波Haar Wavelet、道尔小波Daubechies Wavelet、Symlet小波等。1哈尔小波Haar Wavelet公式与定义哈尔小波是最简单的小波基属于最基础的小波变换它的尺度函数ϕ ( t ) \phi(t)ϕ(t)和小波函数ψ ( x ) \psi(x)ψ(x)在数学上定义如下尺度函数母小波小波函数哈尔小波在时间-频率平面中提供了简单的离散化方式是具有最简单结构的分解方式。它是一致不连续的小波主要是通过一步跳跃式的变化来进行局部信息的提取。适用场景与优缺点1适用场景哈尔小波主要用于图像处理、信号压缩、数字图像和音频的初步处理等。因为其计算简便所以经常用于一些实时处理系统。2优缺点优点简单高效计算速度快易于实现。缺点对信号变化的表达能力差不能很好地捕捉信号的平滑和细腻结构尤其不适合不连续或者具有显著跳变的信号。2道尔小波Daubechies Wavelet公式与定义道尔小波是一类具有正交性的小波具有更好的频率局部化特性。它们通过多项式形式定义最常用的是道尔小波的第一种Db1 Db1小波实际上是哈尔小波以及其后的几种Db2Db3…。道尔小波的阶数即小波函数的平滑程度越高信号的局部变化就越平滑频率局部化性越好。对于更高阶的小波如 Db2尺度函数和小波函数的定义更为复杂。举例来说Db2小波的尺度函数和小波函数分别为相比于Db1此时的小波有更强的平滑性和更高的频率分辨能力适合于信号中的细节和复杂结构的提取。适用场景与优缺点1适用场景信号处理道尔小波特别适用于去噪、信号压缩和重建。图像处理特别是图像压缩算法中有广泛应用。时频分析对平滑的信号或带有一定平稳性的信号非常有效。2优缺点优点平滑性较好(在数学上平滑性好代表小波 / 尺度函数本身连续、高阶可微曲线过渡顺滑)具有较强的时频局部化特性适合平稳信号的处理。缺点计算相对复杂需要更多的计算资源。3Symlet 小波Symlet Wavelet公式与定义Symlet小波是道尔小波的一种变体它通过进一步优化对称性来克服道尔小波在某些信号上产生的失真。Symlet小波同样是正交小波并且它是对称的这使得它在很多应用场景中更有优势尤其在信号的时域重构中。Symlet 小波可以通过以下关系来定义ϕ ( t ) \phi(t)ϕ(t)是尺度函数决定信号的平滑部分。ψ ( x ) \psi(x)ψ(x)是小波函数捕捉信号的细节部分。hk和gk是滤波器的系数具体系数值与具体的Symlet小波阶数有关。N是滤波器的长度通常是2或4表示滤波器的阶数。Symlet 小波的阶数越高其平滑性和对称性越好。以 Symlet-2 为例Symlet-2也称为 S2的尺度函数和小波函数系数如下Symlet-2 小波通过这些系数生成尺度函数和小波函数进行时移和尺度缩放等操作。适用场景与优缺点1适用场景Symlet-2小波适用于对信号对称性和平滑性要求较高的应用如音频处理、图像分析和生物医学信号分析等。它特别适合处理平稳或复杂的信号能够有效进行信号去噪、特征提取和压缩。2优缺点优点高对称性Symlet 小波比 Daubechies 小波具有更好的对称性这对于一些需要保留信号对称特征的应用非常有用如音频、图像等。平滑性Symlet 小波的平滑度较高适用于平稳信号和复杂信号的处理。正交性和 Daubechies 小波一样Symlet 小波也保持了正交性能够有效地分解和重构信号适用于信号压缩、去噪等任务。缺点计算复杂度较高由于Symlet小波具有较高的对称性它的计算复杂度通常要比其他低阶小波如Haar小波更高因此在实时信号处理中可能需要更多的计算资源。内存占用大、相对较慢的计算速度由于Symlet小波的多项式阶数较高使用较高阶的Symlet小波时处理需要更多的内存和时间尤其对于长时间序列信号。离散小波变换本质就是用一对数字滤波器(h,g)进行卷积下采样这对滤波器是生成尺度函数、小波函数的种子但滤波器本身不等于ϕ ( t ) \phi(t)ϕ(t)和ψ ( t ) \psi(t)ψ(t)。给定小波基名字db4就唯一确定一套h , g h,gh,g给定h , g h,gh,g也唯一确定ϕ \phiϕ,ψ \psiψ。小波基选择总结表格补充正交性是指小波基之间的内积为零即它们是彼此“独立”的。这意味着小波基可以在信号分解和重构时避免冗余信息。正交小波可以将信号能量分布到不同频带且不同频带之间没有相互干扰。适合信号分解、重构、压缩等应用。双正交性是指两个不同的小波基满足正交性一个用于信号分解另一个用于信号重构。双正交小波可以通过不同的小波基实现更灵活的信号分析和重构。它同时保留了正交小波的优点并在某些方面提供更高效的处理能力如更好的边界处理。紧支性是指小波函数的非零值仅在有限的范围内存在即小波函数的“支撑”区域是有限的。紧支性小波在计算时可以快速收敛且在实际应用中计算效率高。紧支性的小波基适合处理具有局部特性的信号。支撑长度是指小波函数在时间域中非零值区域的长度。支撑长度越短时间局部化越强处理短时信号时更加精确。然而支撑长度越短频域局部化的能力可能会下降时间-频率分辨率的权衡。对称性是指小波基函数在时间域的对称程度。对称小波基在处理边界问题时可以减少失真或伪影。常见的小波基中Mexicanhat是完全对称的而Daubechies是近似对称的。小波重构小波重构是小波变换的逆过程即通过逼近信息和细节信息逐层恢复原始信号的过程。重构过程利用多尺度分析的分解结果即逼近系数和细节系数以及相应的小波基函数进行信号的还原。 从第N层开始向上递归进行重构逐层计算每一层的逼近系数和细节系数被用于重构上一层的逼近系数。直到回复原始信号。小波重构特点a完全重构小波变换的正交性或准正交性保证了信号可以完全重构而无损。b局部重构只用某些特定的逼近或细节信息可以实现局部信号的重构常用于压缩或降噪。c滤波器精度重构质量与所使用的小波基、滤波器系数的精度密切相关。小波分解与重构的关系小波去噪原理具体步骤如下1小波分解使用小波变换将信号分解为不同的频率成分逼近系数和细节系数。低频部分逼近系数通常包含了基线漂移。2阈值处理通过设定阈值去除低频分量中的噪声。阈值处理通常应用于小波系数用硬阈值或软阈值来去除噪声。硬阈值是直接将小于阈值的系数置为零而软阈值则是在硬阈值的基础上对大于阈值的系数进行缩减。3小波重构通过逆小波变换将去噪后的小波系数重构回原始信号。这样去除基线漂移后的信号就能恢复为更接近原始的清晰信号。硬阈值方法简单但可能会导致去噪信号出现不连续性和突变尤其在处理低幅信号时容易丢失一些重要的细节信息。 软阈值的优点是能平滑信号减少硬阈值可能产生的“块状效应”但它可能导致信号的一些细节丢失因为它对所有系数进行了缩减。2肌电干扰及其去噪参考文章https://blog.csdn.net/qq_42936724/article/details/145585714?spm1001.2014.3001.5501肌电干扰的产生与特点产生原理肌电干扰EMG interference通常是由肌肉的收缩或放松引起的。当肌肉纤维收缩时会产生动作电位这些电位会传播并产生电场被电极检测到形成肌电信号。然而这些信号也可能被其他设备或电极误检测到从而在非目标信号中产生干扰。例如在心电图ECG信号中肌电干扰会以不规则且快速变化的锯齿波形出现如下图** 信号特点与特征 **肌电干扰的频率范围较宽主要集中在10-150Hz之间。其频谱特性类似于零均值白高斯噪声。在静息状态下肌电干扰通常较小但在一些特殊情况下如患有无法控制的震颤的受试者、需要努力保持肢体或身体姿势的残疾人以及儿童中肌电干扰较为常见。带通滤波器处理肌电干扰带通滤波器处理肌电干扰通带频率范围根据需要去除的肌电干扰频率范围来设定例如若肌电干扰集中在10-120Hz可设置带通滤波器的通带范围为该区间之外的频率。滤波器阶数滤波器的阶数越高其频率选择性越好但可能会引入更多的相位失真。滤波器类型常见的有巴特沃斯滤波器、切比雪夫滤波器等不同类型的滤波器在通带平坦度、阻带衰减等方面存在差异。原始信号的快速震荡的小幅度高频波动可以被滤除使心电信号变得更加平顺与此同时明显的基线漂移也可以被滤除。通过合理设置带通滤波器的截止频率的系数可以做到这一点。3工频干扰及其去噪「ECG信号处理——9工频干扰去噪与实现方法陷波/带阻滤波器」2025年3月4日工频干扰概述工频干扰的原因工频干扰主要是由电力系统中的交流电源50Hz或60Hz产生的电磁波。这些干扰通过电极接触、导线等途径进入到心电信号中表现为周期性、规则的波动。在许多国家和地区电力系统的标准频率为50Hz如欧洲、亚洲等地区而在美国和其他地区为60Hz。工频干扰的频率范围工频干扰的主要频率为50Hz干扰信号会影响心电图的质量导致测量误差。其频率范围通常集中在50Hz及其谐波上影响ECG信号的准确性。下图为带有工频干扰的心电信号带阻滤波器去除工频干扰陷波滤波器的原理陷波滤波器Notch Filter又称为 带阻滤波器Band-Stop Filter, BSF但通常特指一种 用于抑制特定窄带频率分量的滤波器特别是在工频干扰如 50Hz 或 60Hz去除时应用广泛。陷波滤波器是一种专门设计用于抑制特定频率信号的滤波器。在处理工频干扰时通常设计一个以50Hz为中心频率的陷波滤波器。该滤波器将50Hz及其附近频率的信号有效地抑制掉同时几乎不影响其他频率成分。陷波滤波器的设计陷波滤波器具有非常窄的带宽可以在频谱中为50Hz提供一个“陷阱”将其从信号中去除。设计时需要保证其频率响应在50Hz频率附近下降到足够低的水平以实现干扰信号的有效抑制。常见参数包括中心频率50Hz带宽0.5~2Hz具体取决于干扰的特性优缺点优点陷波滤波器可以非常精确地去除50Hz工频干扰对心电信号的影响较小不会削弱其他频率成分。缺点当工频干扰的频率发生变化时例如电源不稳定陷波滤波器的性能可能下降。因此设计时需要考虑滤波器的频率响应范围和适应性。可以通过信噪比和均方误差判断带阻滤波器的效果。