【肌电信号】基于matlab GUI脉搏信号分析(去噪+特征提取)【含Matlab源码 862期】
一、脉搏信号分析简介
人体脉搏信号是一种微弱信号, 信噪比较低。在检测和采集时, 由于受仪器、人体等方面的影响, 所采集的信号中常存在如下3种噪声:1) 基线漂移、人体呼吸等低频干扰, 频率小于1 Hz;2) 由于肢体抖动、肌肉紧张而引起的干扰, 它的频率范围较大;3) 工频干扰, 是固定频率的干扰, 频率为50 Hz。
这些干扰信号会极大的影响对脉搏信号进行进一步的识别与分析;因此, 在对脉搏信号进行进一步处理之前必须对噪声进行处理, 过去常常采用快速傅里叶变换或者数字滤波器等方法进行去噪, 去噪结果容易产生相位失真, 对脉搏波信号的去噪效果并不好。由于脉搏波信号属于非平稳信号, 而小波变换具有多分辨率的特点, 因此, 采用小波分析能够较好的处理脉搏波等非平稳信号。基于小波变换的信号去噪方法, 一般有阈值法、平移不变量法和模极大值法。
1 小波变换阈值法去噪原理
小波阈值降噪的基本思想是:用如下模型代表一个含噪声的一维信号:
其中s (t) 为原始信号, n (t) 为服从N (0, σ2) 的高斯白噪声。根据式 (1) , 对f (t) 作离散小波变换, 因为正交小波变换具有很强的数据相关性, 进行小波变化时, 它能够把信号的能量集中在一些大的有限的系数中, 而噪声的能量却分布于整个小波域内, 变换后f (t) 的小波系数一部分为信号对应的小波系数, 另一部分为噪声对应的小波系数;并且幅值比较大的小波系数一般以信号为主, 而幅值较小的系数在很大程度上是噪声。总体上来讲, 对f (t) 进行小波分解后, 信号的系数要大于噪声的系数, 这样就可以选择一个合适的临界阈值λ, 如果分解后得到的系数大于这个临界阈值λ时, 就认为此时的分解系数主要是由信号引起的, 就保留这个系数 (硬阈值方法) 或者按照某一固定量向零收缩 (软阈值方法) 这个系数, 如果分解系数小于这个临界阈值λ, 就认为此时的分解系数主要是由噪声引起的, 直接舍弃分解系数;经过这一步骤后, 用得到的小波系数进行小波重构, 就能去除噪声信号n (t) 。
小波阈值法的主要步骤如下:
-
计算含噪声信号的正交小波变换。选择合适的小波和小波分解层数j, 将含噪信号进行小波分解至j层, 得到相应的小波分解系数。
-
对分解得到的小波系数进行阈值处理, 其阈值的处理方法有2种:
其中s表示阈值处理后的信号, t表示阈值。 -
进行小波逆变换。将经阈值处理过的小波系数进行重构, 得到去噪后的信号。分解后的信号如果不做阈值处理, 则重构后的信号仍为原信号。用阈值函数进行处理并重构后, 将得到去噪后的信号, 这个信号经过小波处理后, 去除了噪声, 逼近于原始信号。
在小波阈值方法中, 阈值的选取至关重要, 选用不同的阈值, 滤波信噪比将有明显的差别。若阈值太大, 则信号细节损失太多, 有可能会造成信号失真;而阈值太小, 则达不到预期的去噪效果[4]。虽然阈值选取方法不少, 如常用的固定门限阈值、Stein无偏似然法估计阈值、启发式阈值和最小最大准则门限阈值等[5], 但是阈值的选取不是一个简单问题, 它需考虑抑制无用噪声与保留信号细小变化之间的权衡问题。
图1为一段含噪的脉搏波信号, 采用阈值法对其进行降噪处理, 在matlab7.0平台上, 经过反复仿真, 选择db9小波进行八层分解。
图1 一段含噪脉搏波
图2 小波阈值法去噪结果
经过小波阈值处理后得到的信号波形如图2所示, 直观上可以看出噪声得到了比较好的抑制;但是, 在脉搏波交替的地方出现较大的上、下峰值, 这些峰值并不是原始脉搏波本身所包含的, 而是在去噪的过程中产生的人为干扰, 这种现象类似于Fourier去噪时所产生的伪吉布斯现象。
要避免伪吉布斯现象, 可以采用小波平移不变量去噪法。
2 小波平移不变量去噪法
小波平移不变量去噪法即在一定的平移量范围内对所获得的信号按某一个平移量进行平移, 获得一个在时域上与原始信号具有一定相位差的信号, 然后对这个信号进行去噪处理, 得到去噪后的信号, 对该信号做反向的平移, 则得到与原始信号同相位的信号, 然后按另一个平移量进行平移——去噪——反平移, 重复这一过程, 对所获得结果求平均, 就得到带噪信号去噪后的估计信号。这就是基于平移不变的小波阈值去噪的思想。
采用平移不变量法对图1的原始脉搏波进行去噪处理, 得到波形图如图3所示, 从图中可以看出伪吉布斯现象得到了很好的抑制。
图3 平移不变量法去噪结果
3 小波变换模极大值滤波法
对检测的脉搏信号中混入的随机噪音, 由于随机噪音的频谱很不规则, 与有效信号的频谱差别不大, 可以利用随机噪声的奇异性加以滤波。奇异性的大小用Lipshczti指数来度量。随机噪音的Lipshczti指数与有效信号本身的奇异点的Lipshczti指数大小不一样, 从而它们的小波变换模的极大值在不同尺度下的传播行为也不一样[7], 随尺度变大而幅值减小的, 认为是噪声产生的模极大值, 予以去除;反之, 则认为是由信号奇异点产生的, 予以保留。对所有保留的模极大值进行重建, 即得到去噪后的信号[8]。利用这一特征可将有效信号从随机噪音中提取出来。
具体算法如下:
- 根据信号的特征, 确定小波基、分解层数, 并对含噪声的信号进行小波分解;
- 在最大分解尺度2 j上搜索极值点, 并设定阈值去除小的模极值点, 得到{wj};
- 在j层极大值点的位置, 为2j-1尺度的极值点构造一个搜索领域, 通常为2 j尺度模极值点位置±3, ±4个数据点;
- 根据搜索范围, 把2j-1尺度极大值点落在该区域的点保留, 其他的值置0;
- 令j=j-1, 重复操作, 直到j=2;
- 在j=2层存在极值点的位置查询j=1时相应极值点, 其余位置将极值点置为0;
- 重构。
采用模极大值法去噪后的脉搏波信号如图4所示, 从图中可以看出处理后的信号较为光滑。
4 结论
文中分别采用了3种小波去噪方法对脉搏波信号进行了去噪处理, 采用小波变换阈值去噪法去噪后, 信号主要噪声得到抑制, 且能反映原始信号的特征的尖峰点得到很好的保留。但是在信号的不连续点处, 去噪后会出现伪吉布斯现象, 且阈值的选择对去噪效果有着很重要的影响。为此, 采用平移不变量法进行处理, 处理后, 伪吉布斯现象得到了很好的抑制。模极大值法能有效地保留信号地奇异点信息, 去噪后的信号没有多余振荡, 是原始信号的一个非常好的估计, 但计算量很大, 速度较慢。
二、部分源代码
%% Program Start
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%% 名称：信号处理课程设计
%% 功能：脉搏检测系统GUI主界面
%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%
function varargout = mainjiemian(varargin)
% MAINJIEMIAN MATLAB code for mainjiemian.fig
% MAINJIEMIAN, by itself, creates a new MAINJIEMIAN or raises the existing
% singleton*.
%
% H = MAINJIEMIAN returns the handle to a new MAINJIEMIAN or the handle to
% the existing singleton*.
%
% MAINJIEMIAN('CALLBACK',hObject,eventData,handles,...) calls the local
% function named CALLBACK in MAINJIEMIAN.M with the given input arguments.
%
% MAINJIEMIAN('Property','Value',...) creates a new MAINJIEMIAN or raises the
% existing singleton*. Starting from the left, property value pairs are
% applied to the GUI before mainjiemian_OpeningFcn gets called. An
% unrecognized property name or invalid value makes property application
% stop. All inputs are passed to mainjiemian_OpeningFcn via varargin.
%
% *See GUI Options on GUIDE's Tools menu. Choose "GUI allows only one
% instance to run (singleton)".
%
% See also: GUIDE, GUIDATA, GUIHANDLES
% Edit the above text to modify the response to help mainjiemian
% Last Modified by GUIDE v2.5 29-Apr-2021 12:18:07
% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name', mfilename, ...
'gui_Singleton', gui_Singleton, ...
'gui_OpeningFcn', @mainjiemian_OpeningFcn, ...
'gui_OutputFcn', @mainjiemian_OutputFcn, ...
'gui_LayoutFcn', [] , ...
'gui_Callback', []);
if nargin && ischar(varargin{1})
gui_State.gui_Callback = str2func(varargin{1});
end
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before mainjiemian is made visible.
function mainjiemian_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to mainjiemian (see VARARGIN)
% Choose default command line output for mainjiemian
handles.output = hObject;
ha=axes('units','normalized','position',[0 0 1 1]);
uistack(ha,'down')
II=imread('1.jpg');
image(II)
colormap gray
set(ha,'handlevisibility','off','visible','off');
% Update handles structure
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes mainjiemian wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = mainjiemian_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes on button press in togglebutton1.
function togglebutton1_Callback(hObject, eventdata, handles)
untitled;
% hObject handle to togglebutton1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of togglebutton1
% --- Executes on button press in togglebutton2.
function togglebutton2_Callback(hObject, eventdata, handles)
yuchuli;
% hObject handle to togglebutton2 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of togglebutton2
% --- Executes on button press in togglebutton3.
function togglebutton3_Callback(hObject, eventdata, handles)
tezheng;
% hObject handle to togglebutton3 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Hint: get(hObject,'Value') returns toggle state of togglebutton3
% --- Executes when figure1 is resized.
function figure1_SizeChangedFcn(hObject, eventdata, handles)
% hObject handle to figure1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
if nargout
[varargout{1:nargout}] = gui_mainfcn(gui_State, varargin{:});
else
gui_mainfcn(gui_State, varargin{:});
end
% End initialization code - DO NOT EDIT
% --- Executes just before jidian is made visible.
function jidian_OpeningFcn(hObject, eventdata, handles, varargin)
% This function has no output args, see OutputFcn.
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% varargin command line arguments to jidian (see VARARGIN)
% Choose default command line output for jidian
handles.output = hObject;
ha=axes('units','normalized','position',[0 0 1 1]);
uistack(ha,'down')
II=imread('666.jpg');
image(II)
colormap gray
set(ha,'handlevisibility','off','visible','off');
axes(handles.axes1);
axes(handles.axes2);
axes(handles.axes3);
axes(handles.axes4);
axes(handles.axes5);
axes(handles.axes6);
% Update handles structure
guidata(hObject, handles);
% UIWAIT makes jidian wait for user response (see UIRESUME)
% uiwait(handles.figure1);
% --- Outputs from this function are returned to the command line.
function varargout = jidian_OutputFcn(hObject, eventdata, handles)
% varargout cell array for returning output args (see VARARGOUT);
% hObject handle to figure
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% Get default command line output from handles structure
varargout{1} = handles.output;
% --- Executes when figure1 is resized.
function figure1_SizeChangedFcn(hObject, eventdata, handles)
% hObject handle to figure1 (see GCBO)
% eventdata reserved - to be defined in a future version of MATLAB
% handles structure with handles and user data (see GUIDATA)
% --- Executes on button press in pushbutton1.
function pushbutton1_Callback(hObject, eventdata, handles)
%%%%%把代码放到这，运行时自动加载
[filename,filepath]=uigetfile('C:\Users\dell\Desktop\信号处理课设-----杨硕\初始脉搏数据\MaiBobefore.txt');
filename=[filepath,filename];
[t,Pluse_pre]=textread(filename,'%f%f','headerlines',1);%读入2个浮点值，并跳过文档的第1行
[m,n]=size(Pluse_pre);
n=3;
s3=Pluse_pre;
%%%%%%%%%%%%—————提取2000个点进行数据处理——————%%%%%%%%%%%%%%%%%%%%
fs=360;%采样率
x0=s3(1:2000);%取1到2000共2000个点
t=1:length(x0);%length(x0)指x0数组元素的个数
axes(handles.axes1);
plot(t,x0)
xlabel('采样点');
ylabel('magtitude');
title('标准脉搏信号')
box1=msgbox('正在加载请稍候','提示');
%%%%%%%%%%%%%%%%%%%%%%%%%%%%初步去除基线漂移%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%中值滤波%%%%%%%%%%%%%%%%%%%%%%%%%%%%
L1=medfilt1(x0,330); %一维中值滤波，x0为数组，即要处理原始波形，n是中值滤波器的参数，L1是滤波以后的结果（数组）
L2=x0-L1;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%频谱%%%%%%%%%%%%%%%%%%%%%%%%
N=length(x0);%样点个数
df=fs/(N-1);%分辨率
f=(0:N-1)*df;%其中每点的频率
Y=fft(L2(1:N))/N*2;%真实的幅值
axes(handles.axes2);
plot(f(1:N/2),abs(Y(1:N/2)));%傅里叶变换后的频谱图是对称的，这里需要一半就可以了
xlabel('频率/Hz');
ylabel('振幅');
axis ( [0 100 0 0.4] );
title('中值滤波后脉搏信号频率谱')
drawnow;
三、运行结果
四、matlab版本及参考文献
1 matlab版本
2014a
2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.
[4]徐洁.基于小波分析的脉搏波信号处理[J].电子设计工程. 2013,21(11)
文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。
原文链接：qq912100926.blog.csdn.net/article/details/116266676
- 点赞
- 收藏
- 关注作者
评论（0）