通过Matlab实现Hermite基函数进行信号拟合,可应用于信号降噪

2024-05-06 18:28

本文主要是介绍通过Matlab实现Hermite基函数进行信号拟合,可应用于信号降噪,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

利用Hermite基函数的Hermite近似,在不牺牲精度的情况下,实现对时序信号的降噪,文中图片以心电信号QRS波群滤除高频干扰为例。

1.知识背景

        Hermite正交多项式是一类重要的正交多项式,它们起源于数学中的Hermite函数和特殊函数理论,它们在物理学、工程学和数学中都有广泛的应用。

        Hermite正交多项式满足一组特殊的正交性质,即在一定的权函数下,它们彼此正交。在通常的情况下,这个权函数是e^(-x^2),即高斯权函数。这种正交性质使得Hermite多项式在傅立叶变换、量子力学、概率论等领域有着重要的应用。

        Hermite正交多项式的一大特点是其递归关系。设Hn(x)表示Hermite多项式的第n阶多项式,则有如下递归关系:

H_n(x) = 2xH_{n-1}(x)-2(n-1)H_{n-2}(x)

        Hermite正交多项式是通过递归关系定义的,下面是前六个Hermite正交多项式的表达式:

1. ~H_0(x)=1\\ 2. ~H_1(x)=2x\\ 3. ~ H_2(x)=4x^2-2 \\ 4. ~H_3(x)=8x^3-12x\\ 5. ~ H_4(x)=16x^4-48x^2+12\\ 6. ~H_5(x)=35x^5-160x^3+120x

        Hermite基函数是以Hermite多项式为基础构建的一组函数集合,通常用于函数逼近和信号处理等领域,Hermite基函数的形式如下:

\Phi _n(t,\sigma ) = \frac{1}{\sqrt{\sigma 2^nn!\sqrt{\pi}}}e^{-\frac{t^2}{2\sigma ^2}}H_n(t/\sigma )

其中,\Phi _n​表示第n个Hermite基函数,H_n是对应的第n个Hermite多项式。Hermite基函数的特点是它们在实数轴上有限且彼此正交,因此在信号处理领域常用于构建正交基函数系列,用于信号分解滤波和重建等任务。

        σ为30时,Hermite基函数\Phi _n图像如下所示(本文取前六个基函数对信号进行降噪):

2.最终拟合

        通过Hermite基函数对信号进行拟合实现降噪,拟合公式如下:

 \hat{x}(t) = \sum_{n=0}^{N-1}c_n(\sigma)\Phi _n(t,\sigma )

其中,\hat{x}(t)表示拟合后的信号,c_n(\sigma)表示Hermite基函数的系数权重,不同的σ对应不同的权重,对于给定的σ,系数c_n(\sigma)计算公式如下:

c_n(\sigma)=\sum_{t}^{}x(t)\Phi _n(t,\sigma )

        最佳权重c_n(\sigma)的选择根据均方误差(Mean Squared Error,MSE)的最小值进行选取,MSE公式如下:

MSE=\frac{1}{n}\sum_{i=1}^{n}(\hat{x}_i-x_i)^2

其中x_i表示拟合前的信号。

拟合效果如下图所示:

3.应用

该方法可应用于心电信号QRS波群滤除高频干扰,同理可应用于其他时序信号的滤波。

如需代码,请加下方微信:

这篇关于通过Matlab实现Hermite基函数进行信号拟合,可应用于信号降噪的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/965036

相关文章

中文分词jieba库的使用与实景应用(一)

知识星球:https://articles.zsxq.com/id_fxvgc803qmr2.html 目录 一.定义: 精确模式(默认模式): 全模式: 搜索引擎模式: paddle 模式(基于深度学习的分词模式): 二 自定义词典 三.文本解析   调整词出现的频率 四. 关键词提取 A. 基于TF-IDF算法的关键词提取 B. 基于TextRank算法的关键词提取

水位雨量在线监测系统概述及应用介绍

在当今社会,随着科技的飞速发展,各种智能监测系统已成为保障公共安全、促进资源管理和环境保护的重要工具。其中,水位雨量在线监测系统作为自然灾害预警、水资源管理及水利工程运行的关键技术,其重要性不言而喻。 一、水位雨量在线监测系统的基本原理 水位雨量在线监测系统主要由数据采集单元、数据传输网络、数据处理中心及用户终端四大部分构成,形成了一个完整的闭环系统。 数据采集单元:这是系统的“眼睛”,

hdu1043(八数码问题,广搜 + hash(实现状态压缩) )

利用康拓展开将一个排列映射成一个自然数,然后就变成了普通的广搜题。 #include<iostream>#include<algorithm>#include<string>#include<stack>#include<queue>#include<map>#include<stdio.h>#include<stdlib.h>#include<ctype.h>#inclu

csu 1446 Problem J Modified LCS (扩展欧几里得算法的简单应用)

这是一道扩展欧几里得算法的简单应用题,这题是在湖南多校训练赛中队友ac的一道题,在比赛之后请教了队友,然后自己把它a掉 这也是自己独自做扩展欧几里得算法的题目 题意:把题意转变下就变成了:求d1*x - d2*y = f2 - f1的解,很明显用exgcd来解 下面介绍一下exgcd的一些知识点:求ax + by = c的解 一、首先求ax + by = gcd(a,b)的解 这个

hdu1394(线段树点更新的应用)

题意:求一个序列经过一定的操作得到的序列的最小逆序数 这题会用到逆序数的一个性质,在0到n-1这些数字组成的乱序排列,将第一个数字A移到最后一位,得到的逆序数为res-a+(n-a-1) 知道上面的知识点后,可以用暴力来解 代码如下: #include<iostream>#include<algorithm>#include<cstring>#include<stack>#in

hdu1171(母函数或多重背包)

题意:把物品分成两份,使得价值最接近 可以用背包,或者是母函数来解,母函数(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v)(1 + x^v+x^2v+.....+x^num*v) 其中指数为价值,每一项的数目为(该物品数+1)个 代码如下: #include<iostream>#include<algorithm>

【C++】_list常用方法解析及模拟实现

相信自己的力量,只要对自己始终保持信心,尽自己最大努力去完成任何事,就算事情最终结果是失败了,努力了也不留遗憾。💓💓💓 目录   ✨说在前面 🍋知识点一:什么是list? •🌰1.list的定义 •🌰2.list的基本特性 •🌰3.常用接口介绍 🍋知识点二:list常用接口 •🌰1.默认成员函数 🔥构造函数(⭐) 🔥析构函数 •🌰2.list对象

【Prometheus】PromQL向量匹配实现不同标签的向量数据进行运算

✨✨ 欢迎大家来到景天科技苑✨✨ 🎈🎈 养成好习惯,先赞后看哦~🎈🎈 🏆 作者简介:景天科技苑 🏆《头衔》:大厂架构师,华为云开发者社区专家博主,阿里云开发者社区专家博主,CSDN全栈领域优质创作者,掘金优秀博主,51CTO博客专家等。 🏆《博客》:Python全栈,前后端开发,小程序开发,人工智能,js逆向,App逆向,网络系统安全,数据分析,Django,fastapi

zoj3820(树的直径的应用)

题意:在一颗树上找两个点,使得所有点到选择与其更近的一个点的距离的最大值最小。 思路:如果是选择一个点的话,那么点就是直径的中点。现在考虑两个点的情况,先求树的直径,再把直径最中间的边去掉,再求剩下的两个子树中直径的中点。 代码如下: #include <stdio.h>#include <string.h>#include <algorithm>#include <map>#

让树莓派智能语音助手实现定时提醒功能

最初的时候是想直接在rasa 的chatbot上实现,因为rasa本身是带有remindschedule模块的。不过经过一番折腾后,忽然发现,chatbot上实现的定时,语音助手不一定会有响应。因为,我目前语音助手的代码设置了长时间无应答会结束对话,这样一来,chatbot定时提醒的触发就不会被语音助手获悉。那怎么让语音助手也具有定时提醒功能呢? 我最后选择的方法是用threading.Time