【信号处理】基于HMM的睡眠状态检测matlab源码

Matlab科研辅导帮

已于 2024-09-23 22:35:18 修改

阅读量8

点赞数

文章标签： matlab 人工智能深度学习计算机视觉数据挖掘

于 2021-06-29 23:06:56 首次发布

本文链接：https://blog.csdn.net/m0_60703264/article/details/120135632

版权

在每个睡眠阶段采集检测目标的睡眠生理信号，睡眠生理信号包括心率信号和呼吸信号；对呼吸信号进行预处理得到呼吸轮廓信号；对心率信号进行预处理得到心率周期的均匀时间序列；分别提取呼吸轮廓信号和均匀时间序列的特征信息，并根据预设的训练样本集获取该睡眠生理信号对应的先验知识；根据特征信息和先验知识以及检测目标的目标特征生成特征集合，并通过主成分分析法从特征集合中确定最优特征子集合得到检测样本；将检测样本与样本集训练得到的检测模型进行匹配，并根据匹配结果输出睡眠呼吸障碍检测信息。

``` function varargout = SleepStatus(varargin) % SLEEPSTATUS MATLAB code for SleepStatus.fig % SLEEPSTATUS, by itself, creates a new SLEEPSTATUS or raises the existing % singleton. % % H = SLEEPSTATUS returns the handle to a new SLEEPSTATUS or the handle to % the existing singleton. % % SLEEPSTATUS('CALLBACK',hObject,eventData,handles,...) calls the local % function named CALLBACK in SLEEPSTATUS.M with the given input arguments. % % SLEEPSTATUS('Property','Value',...) creates a new SLEEPSTATUS or raises the % existing singleton*. Starting from the left, property value pairs are % applied to the GUI before SleepStatusOpeningFcn gets called. An % unrecognized property name or invalid value makes property application % stop. All inputs are passed to SleepStatusOpeningFcn 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 SleepStatus

% Last Modified by GUIDE v2.5 21-Oct-2020 17:16:55

% Begin initialization code - DO NOT EDIT guiSingleton = 1; guiState = struct('guiName', mfilename, ... 'guiSingleton', guiSingleton, ... 'guiOpeningFcn', @SleepStatusOpeningFcn, ... 'guiOutputFcn', @SleepStatusOutputFcn, ... 'guiLayoutFcn', [] , ... 'guiCallback', []); if nargin && ischar(varargin{1}) guiState.gui_Callback = str2func(varargin{1}); end

if nargout [varargout{1:nargout}] = guimainfcn(guiState, varargin{:}); else guimainfcn(guiState, varargin{:}); end % End initialization code - DO NOT EDIT

% --- Executes just before SleepStatus is made visible. function SleepStatus_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 SleepStatus (see VARARGIN)

% Choose default command line output for SleepStatus handles.output = hObject;

set(handles.text2,'string','Sleep status detector','foregroundcolor','b','fontsize',20);

set(handles.pushbutton1,'String','打开文件'); set(handles.pushbutton2,'String','播放且判断'); set(handles.pushbutton3,'String','停止');

set(handles.edit1,'String','实时判断结果'); set(handles.edit2,'String','过往判断结果'); set(handles.edit2,'Max',10); set(handles.edit2,'Min',0);

cla(handles.axes1); cla(handles.axes2);

global player timercounter accuratetime nlast nlast1 allresult num player=[]; timercounter =0; accuratetime = 0; nlast =2; nlast1 =2; num =1; allresult={}; % Update handles structure guidata(hObject, handles);

% UIWAIT makes SleepStatus wait for user response (see UIRESUME) % uiwait(handles.figure1);

% --- Outputs from this function are returned to the command line. function varargout = SleepStatus_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 pushbutton1. function pushbutton1Callback(hObject, eventdata, handles) % hObject handle to pushbutton1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global wavfilename timercounter allresult num [FileName,PathName] = uigetfile('*.wav','Select the wav-file'); if FileName==0 wavfilename =[]; return end wavfilename = [PathName,FileName]; [x,fs] = audioread(wavfilename); x = mean(x,2); axes(handles.axes1); plot(x); ylim([-1 1]); timercounter = 0; accuratetime =0; allresult={}; set(handles.edit2,'string',''); num = 1;

% --- Executes on button press in pushbutton2. function pushbutton2Callback(hObject, eventdata, handles) % hObject handle to pushbutton2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global wavfilename tAxis player myStruct startminute start_seconds if isempty(wavfilename) return end %default

durT = 3060; %seconds fs = 32000; durS = fsdurT; %samples

[x,fs] = audioread(wavfilename); x = mean(x,2); durT = size(x,1)/fs; durS = size(x,1);

dt = 1/fs; tAxis = dt:dt:durT;

frameRate = 25; %fps frameT = 1/frameRate;

mag = 1;

axes(handles.axes1); cla(handles.axes1); plot(tAxis, x); ylim([-mag mag]) xlim([0 durT]) xlabel('Time [s]')

playHeadLoc = 0; % hold on; ax = plot([0 playHeadLoc], [-mag mag], 'r', 'LineWidth', 2);

player = audioplayer(x, fs); myStruct.playHeadLoc = playHeadLoc; myStruct.frameT = frameT; myStruct.ax1 = handles.axes1; myStruct.ax2 = handles.axes2; myStruct.edit1 = handles.edit1; myStruct.edit2 = handles.edit2; tmr1 = timer('ExecutionMode', 'FixedRate', ... 'Period', 0.5, ... 'TimerFcn', {@timerCallback}); % tmr2 = timer('ExecutionMode', 'FixedRate', ... % 'Period', 0.1, ... % 'TimerFcn', {@timerCallback2}); handles.timer1 = tmr1; % handles.timer2 = tmr2;

guidata(hObject, handles);

starttime = clock(); startseconds = starttime(6); startminute = start_time(5);

start(tmr1); % start(tmr2);

playblocking(player); % play(player); if ~isplaying(player)
stop(handles.timer1); % stop(handles.timer2); end

% function timerCallback2(hObj, eventdata) % global accuratetime % accuratetime = accurate_time + 0.1;

function timerCallback(hObj, eventdata) global wavfilename myStruct tAxis timercounter accuratetime startminute startseconds nlast nlast1 allresult num [x1,fs] = audioread(wavfilename); x1= mean(x1,2); lengthmax = length(x1);

% timercounter = accuratetime;

starttime = clock(); seconds = starttime(6); minute = start_time(5);

timercounter = (minute*60+seconds)-(startminute*60+startseconds); % disp(timercounter) axes(myStruct.ax1); plot(tAxis,x1); line([timercounter timercounter],[-1 1],'color','r');

axis([tAxis(1) tAxis(end) -1 1]);

time1 = timer_counter; time2 = time1 + 1;%1s

idx1 = fix(fstime1); idx2 = fix(fstime2); if (idx1<1) idx1 = 1; end if(idx1>=lengthmax) idx1 = lengthmax; end if(idx2>=lengthmax) idx2 = lengthmax; end if(idx1>=idx2) return end x1 = x1([idx1:idx2]); axes(myStruct.ax2); plot(x1,'color',[1, 0.6732, 0]); axis([1 1*fs+1 -.5 .5]);

%判断结果 load hmm.mat addpath (genpath('HMMall'),'-end');

warning('off'); classfolder={'NormalBreathe','FastBreathe','Snore'}; ResultStr = classfolder; k = x1; datacell = MFCC(k,fs); % size(datacell) datacell = datacell'; for j=1:size(hmm,2) HMM = hmm{j};
loglik(j) = mhmm_logprob(datacell,HMM.prior, HMM.transmat, HMM.mu, HMM.Sigma, HMM.mixmat);%计算概率 end

[d,n] = max(loglik);

% fprintf('文件名:%s，识别结果:%s\n', [classfolder{i},'\',fname],char(ResultStr{n})); njudge = nlast1;%延迟1秒 if njudge ==1 set(myStruct.edit1,'string',[num2str(time1),' s: ',char(ResultStr{njudge})]); set(myStruct.edit1,'string',[num2str(time1),' s: ',char(ResultStr{njudge})]); set(myStruct.edit1,'foregroundcolor','r'); set(myStruct.edit1,'fontsize',20); allresult{num} = num2str(time1); num = num+1; set(myStruct.edit2,'string',allresult); set(myStruct.edit2,'foregroundcolor','r'); elseif(njudge==2) set(myStruct.edit1,'string',[num2str(time1),' s: ',char(ResultStr{njudge})]); set(myStruct.edit1,'foregroundcolor','g'); set(myStruct.edit1,'fontsize',20); else set(myStruct.edit1,'string',[num2str(time1),' s: ',char(ResultStr{njudge})]); set(myStruct.edit1,'foregroundcolor','b'); set(myStruct.edit1,'fontsize',20); end nlast1 = nlast; nlast = n; %结果留着下下次显示

rmpath (genpath('HMMall'));

warning('on'); % timercounter = accuratetime; % timercounter = timercounter + 0.5;

function edit1_Callback(hObject, eventdata, handles) % hObject handle to edit1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of edit1 as text % str2double(get(hObject,'String')) returns contents of edit1 as a double

% --- Executes during object creation, after setting all properties. function edit1_CreateFcn(hObject, eventdata, handles) % hObject handle to edit1 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called

% Hint: edit controls usually have a white background on Windows. % See ISPC and COMPUTER. if ispc && isequal(get(hObject,'BackgroundColor'), get(0,'defaultUicontrolBackgroundColor')) set(hObject,'BackgroundColor','white'); end

% --- Executes on button press in pushbutton3. function pushbutton3_Callback(hObject, eventdata, handles) % hObject handle to pushbutton3 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA) global player if ~isempty(player) stop(player); stop(handles.timer1); % stop(handles.timer2); end

% --- Executes when user attempts to close figure1. function figure1_CloseRequestFcn(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)

% Hint: delete(hObject) closes the figure delete(hObject); global player if ~isempty(player) stop(player); stop(handles.timer1); % stop(handles.timer2); end

function edit2_Callback(hObject, eventdata, handles) % hObject handle to edit2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles structure with handles and user data (see GUIDATA)

% Hints: get(hObject,'String') returns contents of edit2 as text % str2double(get(hObject,'String')) returns contents of edit2 as a double

% --- Executes during object creation, after setting all properties. function edit2_CreateFcn(hObject, eventdata, handles) % hObject handle to edit2 (see GCBO) % eventdata reserved - to be defined in a future version of MATLAB % handles empty - handles not created until after all CreateFcns called