均匀采样信号的鲁棒Savistky-Golay滤波(MATLAB)

S-G滤波器又称S-G卷积平滑器,它是一种特殊的低通滤波器,用来平滑噪声数据。该滤波器被广泛地运用于信号去噪,采用在时域内基于多项式最小二乘法及窗口移动实现最佳拟合的方法。与通常的滤波器要经过时域-频域-时域变换不同,S-G滤波直接处理时域数据进行平滑,其平滑效果随窗口宽度不同而不同。相对于均值平滑滤波,S-G滤波更能保留相对极大值、极小值和宽度等分布特征。该滤波算法的另一优点是其运算量相对较小,对计算机的内存及数据处理能力要求较低。鉴于此,采用一种鲁棒Savistky-Golay滤波器对均匀采样信号进行降噪平滑,运行环境MATLAB 2021。


for i = 1:step:Nind1 = max(1,i-half_window) ;ind2 = min(N,i+half_window) ;y_spl = y_input(ind1:ind2) ;x_spl = [ind1:ind2]' ;weights = ones(size(y_spl)) ;p = inf(order+1,1) ; s = inf ; ds = 1 ; dp = 1 ;X = repmat(x_spl,1,order+1).^[order:-1:0] ;c = 4.685 ; nb_iter = 0 ;while (abs(ds) >= opts.TolFun) || (abs(dp) >= opts.TolX)if nb_iter > opts.MaxIterbreakendnb_iter = nb_iter+1 ;old_p = p ;p = (sqrt(weights).*X)\(sqrt(weights).*y_spl) ;y_spl_calc = X*p ;r = y_spl_calc-y_spl ;tau = median(abs(r-median(r)))/0.6745 ;if tau == 0s = std(r)*sqrt(1+1/numel(y_spl)+(i-mean(x_spl)).^2/sum((x_spl-mean(x_spl)).^2)) ;breakendz = r/tau ;weights = (abs(z)<c).*(1-(z/c).^2).^2 ;weights = weights.*(weights>=0) ; old_s = s ;s = std(r,weights)*sqrt(1+1/numel(y_spl)+(i-sum(weights.*x_spl)./sum(weights)).^2/sum((x_spl-sum(weights.*x_spl)./sum(weights)).^2)) ;dp = max((old_p-p)./p) ;ds = (old_s-s)/s ;endy3(i,:) = p' ;y2(i,1) = s ;y1(i,1) = polyval(p,i) ;endif y_input_type == "row"y1 = y1' ;y2 = y2' ;x3 = y3' ;y_input = y_input' ;x_interp = 1:step:N ;x_query = 1:N ;
elsex_interp = [1:step:N]' ;x_query = [1:N]' ;
endif step == 1% No need for interpolation.return
endy1 = interp1(x_interp,y1(1:step:N),x_query,opts.Interpolation,"extrap") ;
y2 = interp1(x_interp,y2(1:step:N),x_query,opts.Interpolation,"extrap") ;
y3 = interp1(x_interp,y3(1:step:N,:),x_query,opts.Interpolation,"extrap") ;

图片

图片

图片

图片

  • 工学博士,担任《Mechanical System and Signal Processing》《中国电机工程学报》《控制与决策》等期刊审稿专家,擅长领域:现代信号处理,机器学习,深度学习,数字孪生,时间序列分析,设备缺陷检测、设备异常检测、设备智能故障诊断与健康管理PHM等。代码通过知乎学术咨询获得:https://www.zhihu.com/consult/people/792359672131756032?isMe=1

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/bicheng/42100.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

掌握 IPython:灵活运用系统环境变量

掌握 IPython&#xff1a;灵活运用系统环境变量 在数据科学和编程的广阔天地中&#xff0c;IPython 以其强大的交互式功能成为了探索和实验的利器。但你知道吗&#xff1f;IPython 同样能够与系统环境变量进行交互&#xff0c;这使得我们能够更灵活地控制我们的代码和环境。本…

ubuntu22 使用ufw防火墙

专栏总目录 一、安装 sudo apt update sudo apt install ufw 二、启动防火墙 &#xff08;一&#xff09;启动命令 sudo ufw enable &#xff08;二&#xff09;重启命令 sudo ufw reload 三、配置规则 #允许SSH连接 sudo ufw allow ssh #如果sshd服务端口指定到了8888&a…

Oracle PL / SQL 插入insert

INSERT是SQL的数据操作语言的一部分的SQL关键字之一。 DML可以操作关系数据库中的数据。 让我们从INSERT语句的第一种形式开始&#xff0c;INSERT ... VALUES。 INSERT ... VALUES 以下代码显示如何使用DML的INSERT VALUES子语句。 INSERT INTO authors ( id, name, birth_d…

RedHat运维-LinuxSELinux基础2-文件绑定SELinux上下文

1. 所有资源&#xff0c;比如进程、文件和________&#xff0c;都被打上了SELinux上下文&#xff1b; 2. 所有资源&#xff0c;比如进程、文件和________&#xff0c;都被打上了SELinux上下文&#xff1b; 3. 所有资源&#xff0c;比如进程、文件和________&#xff0c;都被打上…

HTML5使用<blockquote>标签:段落缩进

使用<blockquote>标签可以实现页面文字的段落缩进。这一标签也是每使用一次&#xff0c;段落就缩进一次&#xff0c;并且可以嵌套使用&#xff0c;以达到不同的缩进效果。语法如下&#xff1a; <blockquote>文字</blockquote> 【实例】使用<blockquote&…

常用 Android 反编译工具apktooldex2jarenjarifyjd-guijadx

apktool 官网地址 : https://ibotpeaches.github.io/Apktool/ 反编译命令 : java -jar apktool.jar d demo.apk -o demoapktool d app.apk -r --only-main-classes -o destDird 表示反编译 app.apk 是目标apk文件路径 -r 表示不反编译资源文件 –only-main-classes 表示不反…

【面向就业的Linux基础】从入门到熟练,探索Linux的秘密(十一)-git(3)

Git是目前最流行的版本控制系统之一&#xff0c;在现代软件开发中扮演着重要的角色。它能够有效地跟踪文件变化、协作开发&#xff0c;并存储项目的历史记录。本文的目的是向读者介绍Git的基本概念和工作原理&#xff0c;帮助初学者快速上手使用Git&#xff0c;并帮助有经验的开…

[目标检测]labelme标注数据转yoloV8需要的.txt格式

1、首先需要知道yoloV8目标检测的标签格式&#xff1a; yolov8标签数据格式说明 After using a tool like Roboflow Annotate to label your images, export your labels to YOLO format, with one *.txt file per image (if no objects in image, no *.txt file is required).…

Java.lang.Thread类和Java的主线程

一.Java.lang.Thread类 支持多线程编程 常用方法 二.主线程 ◆Java程序启动时&#xff0c;一个线程立即随之启动&#xff0c;通常称之为程序的主线程 ◆main()方法即为主线程入口 ◆产生其他子线程的线程 ◆必须最后完成执行&#xff0c;因为它执行各种关闭动作 示例 使用…

0-1背包、完全背包算法模板从递归到记忆化搜索到动态规划

0-1背包、完全背包算法模板从递归到记忆化搜索到动态规划 不管是0-1背包还是完全背包&#xff0c;我们都可以将问题转换成为选择或者不选的问题&#xff0c;这个问题在前面的回溯算法模板是一样的。 举个栗子&#xff1a; 假如有1、4、5这三个数&#xff0c;问组成和为12的数需…

202406 CCF-GESP Python 五级试题及详细答案注释

202406 CCF-GESP Python 五级试题及详细答案注释1 单选题(每题 2 分,共 30 分) 第 1 题 在Python中,print((c for c in "GESP"))的输出是( )。 A. (G, E, S, P) B. [G, E, S, P] C. {G, E, S, P} D. 以上选项均不正确答案:D解析:在Python中,(c for c in&q…

Keras实战之图像分类识别

文章目录 整体流程数据加载与预处理搭建网络模型优化网络模型学习率Drop-out操作权重初始化方法对比正则化加载模型进行测试 实战&#xff1a;利用Keras框架搭建神经网络模型实现基本图像分类识别&#xff0c;使用自己的数据集进行训练测试。 问&#xff1a;为什么选择Keras&am…

这门HCIE认证正式发布!

华为认证AI解决方案架构专家HCIE-AI Solution Architect V1.0&#xff08;中文版&#xff09;自2024年6月28日起&#xff0c;正式在中国区发布。 基于“平台生态”战略&#xff0c;围绕“云-管-端”协同的新ICT技术架构&#xff0c;华为公司打造了覆盖ICT领域的认证体系&#x…

C++ | Leetcode C++题解之第217题存在重复元素

题目&#xff1a; 题解&#xff1a; class Solution { public:bool containsDuplicate(vector<int>& nums) {unordered_set<int> s;for (int x: nums) {if (s.find(x) ! s.end()) {return true;}s.insert(x);}return false;} };

2024年江苏省研究生数学建模科研创新实践大赛C题气象数据高精度融合技术研究论文和代码分析

经过不懈的努力&#xff0c; 2024年江苏省研究生数学建模科研创新实践大赛C题气象数据高精度融合技术研究论文和代码已完成&#xff0c;代码为C题全部问题的代码&#xff0c;论文包括摘要、问题重述、问题分析、模型假设、符号说明、模型的建立和求解&#xff08;问题1模型的建…

Go 语言 UUID 库 google/uuid 源码解析:UUID version1 的实现

google/uuid 库地址 关于 UUID 的总体介绍可以查看这篇文章&#xff0c;其包含阅读此篇文章的前置内容。 UUID version 1 在 RFC 4122 文件中定义&#xff0c;其实现基于节点 ID、时钟序列以及当前时间&#xff08;距离格里历改日【1582年10月15日】 的100纳秒数&#xff0c;具…

服务器BMC基础知识总结

前言 因为对硬件方面不太理解&#xff0c;所以打算先从服务器开始学习&#xff0c;也想和大家一起分享一下&#xff0c;有什么不对的地方可以纠正一下哦&#xff01;谢谢啦&#xff01;互相学习共同成长~ 1.BMC是什么&#xff1f; 官方解释&#xff1a;BMC全名Baseboard Mana…

【深度学习】-WASB-调试说明

要改这么几个地方&#xff1a; 代码仓库&#xff1a;/Desktop/code/python_project/WASB-SBDT-main/ 篮球数据集xx_xx_11.xml只保留最后一个11.xml 并把11下直接放置11 video&#xff1a; 这里的东西被我改了&#xff0c;要以仓库为准

SpringSecurity中文文档(Servlet X.509 Run-As)

X.509 Authentication X. 509证书身份验证的最常见用途是在使用 SSL 时验证服务器的身份&#xff0c;最常见的用法是在浏览器中使用 HTTPS。浏览器会自动检查服务器提供的证书是否已由其维护的受信任证书颁发机构列表中的一个颁发(数字签名)。 您还可以使用带有“相互身份验证…

Java与Selenium配置及常见报错解决方法

Java与Selenium配置及常见报错解决方法 1. 简介 Java与Selenium是自动化测试中常用的工具和技术。正确配置开发环境并处理可能出现的问题&#xff0c;对于顺利开展测试工作至关重要。本文将为您详细介绍Java与Selenium的配置过程&#xff0c;并提供常见报错的解决方法。 2. …