基于小波域优化Savitzky–Golay滤波器的脑电图信号的运动伪影去除方法(MATLAB R2018A)

在获取或采集数据的过程中,不可避免地将噪声引入到数据中,噪声的存在使得原始数据发生变异,对数据的处理及分析产生严重地影响。常用的去噪模型有平滑去噪、均值去噪。其中,平滑去噪又包括移动平均平滑法和Savitzky-Golay卷积平滑法。

Savitzky-Golay是由Savitzky和Golay于1964年提出,之后被广泛地运用于数据平滑除噪,是一种在时域内基于局域多项式最小二乘法拟合的滤波方法。即对样本中的每个数据点邻域内的数据,用一个n阶多项式拟合,此多项式的系数是由最小二乘法准则在拟合误差最小的条件下确定,从而得到该数据点的最佳拟合值,也就是经过SG去噪后的数据值。

鉴于此,采用小波域优化Savitzky--Golay滤波器对脑电图信号的运动伪影进行去除。

Matlab 复制代码
function MAE_D=MAE_PSD(rec, x2,fs)
[Pxxf,freq]=psd_fft(rec, fs);
 Pxxdbf=10*log10(Pxxf);
 [Pxxn,freq]=psd_fft(x2, fs);
 Pxxdbn=10*log10(Pxxn);
%%PSD calculation for delta(0-4Hz)%%
for i=1:46
MAE_d(i)=abs(Pxxdbf(i)-Pxxdbn(i))./(46-0); %%%%delta information is present upto 46 bins, this will change for different sampling frequency
end
MAE_D=(sum(MAE_d)/(fs/2)); %%%%%formulae given in the paper
完整代码:https://mbd.pub/o/bread/ZJ6alZly
end

工学博士,担任《Mechanical System and Signal Processing》《中国电机工程学报》《控制与决策》等期刊审稿专家,擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。

相关推荐
飞哥数智坊1 分钟前
AI编程实战:数据大屏生成初探
人工智能
蚝油菜花3 分钟前
Cua:Mac用户狂喜!这个开源框架让AI直接接管你的电脑,快速实现AI自动化办公
人工智能·开源
蚝油菜花3 分钟前
AutoAgent:无需编程!接入DeepSeek用自然语言创建和部署AI智能体!港大开源框架让AI智能体开发变成填空题
人工智能·开源
每次的天空5 分钟前
Kotlin 作用域函数:apply、let、run、with、also
android·开发语言·kotlin
nuise_5 分钟前
李宏毅机器学习笔记06 | 鱼和熊掌可以兼得的机器学习 - 内容接宝可梦
人工智能·笔记·机器学习
林泽毅6 分钟前
SwanLab x EasyR1:多模态LLM强化学习后训练组合拳,让模型进化更高效
算法·llm·强化学习
小林熬夜学编程8 分钟前
【高并发内存池】第八弹---脱离new的定长内存池与多线程malloc测试
c语言·开发语言·数据结构·c++·算法·哈希算法
刚入门的大一新生14 分钟前
归并排序延伸-非递归版本
算法·排序算法
声网18 分钟前
MiniMax 发布新 TTS 模型 Speech-02,轻松制作长篇有声内容;Meta 高端眼镜年底推出:售价上千美元丨日报
人工智能
独好紫罗兰19 分钟前
洛谷题单3-P1980 [NOIP 2013 普及组] 计数问题-python-流程图重构
开发语言·python·算法