【心音信号】基于matlab GUI心音诊断系统【含Matlab源码 762期】

一、心音诊断系统简介

1 心音
心脏收缩舒张时产生的声音，可用耳或听诊器在胸壁听到，亦可用电子仪器记录下来（心音图）。可分为第一心音（S1）第二心音（S2）。（正常情况下均可听到）。第三心音(S3通常仅在儿童及青少年可听到)，第四心音（S4正常情况很少听到）。从心脏产生的心音经过组织的介导传到胸壁表面，其中以骨传导最好。
心音是心脏及心血管系统机械运动状况的反映，其中包含着心脏各个部分本身及相互之间作用的生理和病理信息。心音信号的识别与分类对心血管系统疾病的诊断具有重要的意义，其准确性、可靠性的好坏决定着诊断与治疗心脏病患者的效果。早期的心音识别与分类是医生根据听诊结果来完成的，显然这一过程具有一定的主观性且可靠性不高。随着信号处理与分析技术的不断发展，对心音的研究也逐步由定性分析进入了定量分析的阶段。国内外许多生物医学工程研究人员将传统的模式识别方法，以及神经网络方法用于心音的识别与分类，期望实现心音的自动解释和自动诊断，以便向临床医生提供实用的辅助诊断信息。此外，心音的识别与分类还有助于对心音产生机制的理解。
心音信号研究主要是采用微电子技术，检测技术，现代数字信号处理技术和生物医学工程技术，研究和揭示心音与心脏病之间的关系。
无论是多大单位的记录，人的体重，心率，血压等生理参数，都是时变的，称为心率或者血压的变异性。
心音是在心动周期中，由于心肌收缩和舒张，瓣膜启闭，血流冲击心室壁和大 动脉等因素引起的机械振动，通过周围组织传到胸壁，将耳紧贴胸壁或将听诊器放在胸壁 定部位，听到的声音 通常很容易听到第一和第二心音，有时在某些情况下听到第三或第四心音。
第一心音：S1,发生在心脏收缩期开始，音调低沉，持续时间较长(约0.15秒)。产生的原因包括心室肌的收缩，房室瓣突然关闭以及随后射血入主动脉等引起的振动。第一心音的最佳听诊部位在锁骨中线第五肋间隙或在胸骨右缘。相对于心电图上QRS波开始后0.02_{0.04s，占时0.08}0.15s
第二心音：S2,发生在心脏舒张期的开始，频率较高，持续时间较短(约0.08秒)。产生的原因是半月瓣关闭，瓣膜互相撞击以及大动脉中血液减速和室内压迅速下降引起的振动。第二心音的最佳听诊部位在第二肋间隙右侧的主动脉瓣区和左侧的肺动脉瓣区。相对于T波终末部。
第三心音和第四心音：
第三心音S3发生在第二心音后0.1～0.2秒，频率低，它的产生与血液快速流入心室使心室和瓣膜发生振动有关，通常仅在儿童能听到，因为 较易传导到体表。相当于T波后距第二心音0.12~0.20s。
第四心音S4由心房收缩引起，也称心房音,相当于心电图上P波后0.15~0.18s，振幅低。

心音杂音对正常心音形成了一定的干扰，但心音杂音的出现对心音信号分析具有实际应用价值和临床意义。
根据杂音出现的时间与S1,S2心音的关系，可分为早，中，晚期杂音。杂音的强度一般可视其振幅与S1比较分类。大于S1,强。小于S1,大于1/3S1，中。低于S1的1/3则为低，仅有轻微振动则为极低振幅杂音。
由于心音信号属于强噪声背景下的人体微弱生物信号，由于心音信号是由复杂的生命体发出的不稳定的自然信号。心音的改变和杂音的出现往往是心脏产生器质性病变的早期症状，心脏内部的物理结构发生变化将直接影响和改变心音信号。

目前广泛应用 的心电图检查是心脏变时性和变传导性的最佳检测方法，但不能用来检测心脏的变力性先天心脏瓣膜受损。心脏传导组织病变引起的心脏机械活动障碍不会首先反应在心电图上，却能首先反应在心音信号上。当冠心病的阻塞达到70%以上时才能引起心电图信号的改变，实际上，达到25%就可以改变心音信号。

2 异常心音
包括S2、S2的异常及收缩期，舒张期的附加音（或额外音）。
① 第一心音异常。指S1增强，减弱或分裂。估计S1的响度，最好是用听诊，心音图判断能力有限。S1增强、减弱或强弱不等的临床情况见表1。S1分裂指M1与T1相距>0.04sec，可见于正常儿童、青年及体瘦者，无重要意义。S1异常分裂时M1与T1相距可大于0.06秒，见于电激动延迟(如右束支传导阻滞等)、机械活动延迟（如房间隔缺损、严重二尖瓣狭窄等。听诊S1分裂在二尖瓣及三尖瓣区最清楚坐位及呼气时更清楚。
② 第二心音异常。包括S2增强，减弱或分裂。
S2增强又分P2增强及A2增强，P2增强见于肺血流量增多（如间隔缺损），肺血管阻力增加；肺静脉压力增高(如二尖瓣狭窄)。P2亢进一般在肺动脉瓣区听到。A2增强见于体循环阻力增高或血流量增多，向肺动脉瓣及心尖区传导，见于 高血压等。
S2减弱又分P2减弱及A2减弱。P2减弱见于肺动脉压低、肺血流量减少或 肺动脉瓣狭窄等。A2减弱见于体循环阻力低、血流减少、血压低、 主动脉瓣狭窄或严重关闭不全。
S2分裂可为生理性。右室射血结束稍晚于左室，吸气时P2延迟出现，此时可听到S2分裂。但呼气时A2与P2接近或重叠，分裂消失。这见于青少年，在肺动脉瓣区听诊明显，坐位时可消失。
S2异常分裂包括宽分裂、固定分裂、逆分裂、分裂减窄。宽分裂是呼气时S2分裂不消失，见于右室射血时间延长或左室射血时间缩短等情况。S2固定分裂指呼吸时A2 -P2间期无明显改变或变动<0.02sec，见于分流量较大的房间隔缺损等。S2逆分裂指A2在P2之后，吸气时A2-P2分裂不显，呼气时P2提早出现，分裂增宽。见于主动脉瓣关闭延迟。S2分裂减窄常由于严重肺动脉高压P2较早出现所致。

3 心音信号分析方法
传统的谱分析方法通过快速傅立叶变换将时，频域关联起来。但FFT时频域分离，并以信号的频率特性时不变，或统计特性稳定为前提。
传统的稳态分析方法反映的是信号的静态频谱特征，对于包括人体心音信号在内的生物学生理信号，由于环境的影响而表现为非平稳时变特性。因此采用经典谱分析方法难以准确反映心音信号的动态变化。
在传统心音分析的基础上，提出了很多方法：
1.短时傅立叶变换（STFT）
2.小波变换
3.其他时频分析方法

二、源代码

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

% Last Modified by GUIDE v2.5 03-May-2011 02:33:20

% Begin initialization code - DO NOT EDIT
gui_Singleton = 1;
gui_State = struct('gui_Name',       mfilename, ...
                   'gui_Singleton',  gui_Singleton, ...
                   'gui_OpeningFcn', @xingyingzhengduan_OpeningFcn, ...
                   'gui_OutputFcn',  @xingyingzhengduan_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 xingyingzhengduan is made visible.
function xingyingzhengduan_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 xingyingzhengduan (see VARARGIN)

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

% Update handles structure
guidata(hObject, handles);

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


% --- Outputs from this function are returned to the command line.
function varargout = xingyingzhengduan_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;


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


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


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


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


% --------------------------------------------------------------------
function New_Callback(hObject, eventdata, handles)
% hObject    handle to New (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
[FileName PathName]=uiputfile('*.*','文件另存为');%创建一个保存文件的对话框，并返回路径。
[s ,mess ,messid] = mkdir(PathName,FileName);%创建用户文件夹

while s~=1
    errordlg( mess , '新建文件出错');
    pause(2);
    helpdlg('请重新选择路径或更改文件名','提示');
    pause(2);
    [FileName PathName]=uiputfile('*.*','文件另存为');
    [s ,mess ,messid] = mkdir(PathName,FileName);
end
pathName = strcat(PathName , '\');
pathName = strcat(pathName , FileName);

handles.pathName = pathName;%定义PathName为全局变量

guidata(hObject , handles);%更新数据


% --------------------------------------------------------------------
function Open_Callback(hObject, eventdata, handles)
% hObject    handle to Open (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)
[fileName pathName]=uiputfile({'*.dat' , '*.txt'},'打开文件路径');

handles.fileName = fileName;%定义FileName为全局变量
handles.PathName = PathName;

guidata(hObject , handles);


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



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


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

clc
clear all
close(gcf);%u关闭当前GUI


% --- Executes on key press with focus on figure1 and none of its controls.
function figure1_KeyPressFcn(hObject, eventdata, handles)
% hObject    handle to figure1 (see GCBO)
% eventdata  structure with the following fields (see FIGURE)
%	Key: name of the key that was pressed, in lower case
%	Character: character interpretation of the key(s) that was pressed
%	Modifier: name(s) of the modifier key(s) (i.e., control, shift) pressed
% handles    structure with handles and user data (see GUIDATA)


% --- Executes on selection change in popupmenu1.
function popupmenu1_Callback(hObject, eventdata, handles)
% hObject    handle to popupmenu1 (see GCBO)
% eventdata  reserved - to be defined in a future version of MATLAB
% handles    structure with handles and user data (see GUIDATA)

% Hints: contents = get(hObject,'String') returns popupmenu1 contents as cell array
%        contents{get(hObject,'Value')} returns selected item from popupmenu1


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

% Hint: popupmenu 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


  

 

  
 
 

  
1
  
2
  
3
  
4
  
5
  
6
  
7
  
8
  
9
  
10
  
11
  
12
  
13
  
14
  
15
  
16
  
17
  
18
  
19
  
20
  
21
  
22
  
23
  
24
  
25
  
26
  
27
  
28
  
29
  
30
  
31
  
32
  
33
  
34
  
35
  
36
  
37
  
38
  
39
  
40
  
41
  
42
  
43
  
44
  
45
  
46
  
47
  
48
  
49
  
50
  
51
  
52
  
53
  
54
  
55
  
56
  
57
  
58
  
59
  
60
  
61
  
62
  
63
  
64
  
65
  
66
  
67
  
68
  
69
  
70
  
71
  
72
  
73
  
74
  
75
  
76
  
77
  
78
  
79
  
80
  
81
  
82
  
83
  
84
  
85
  
86
  
87
  
88
  
89
  
90
  
91
  
92
  
93
  
94
  
95
  
96
  
97
  
98
  
99
  
100
  
101
  
102
  
103
  
104
  
105
  
106
  
107
  
108
  
109
  
110
  
111
  
112
  
113
  
114
  
115
  
116
  
117
  
118
  
119
  
120
  
121
  
122
  
123
  
124
  
125
  
126
  
127
  
128
  
129
  
130
  
131
  
132
  
133
  
134
  
135
  
136
  
137
  
138
  
139
  
140
  
141
  
142
  
143
  
144
  
145
  
146
  
147
  
148
  
149
  
150
  
151
  
152
  
153
  
154
  
155
  
156
  
157
  
158
  
159
  
160
  
161
  
162
  
163
  
164
  
165
  
166
  
167
  
168
  
169
  
170
  
171
  
172
  
173
  
174
  
175
  
176
  
177
  
178
  
179
  
180
  
181
  
182
  
183
  
184
  
185
  
186
  
187
  
188
  
189
  
190
  
191
  
192
  
193
  
194
  
195
  
196
  
197
  
198
  
199
 

三、运行结果

四、matlab版本及参考文献

1 matlab版本
2014a

2 参考文献
[1] 沈再阳.精通MATLAB信号处理[M].清华大学出版社，2015.
[2]高宝建,彭进业,王琳,潘建寿.信号与系统——使用MATLAB分析与实现[M].清华大学出版社，2020.
[3]王文光,魏少明,任欣.信号处理与系统分析的MATLAB实现[M].电子工业出版社，2018.
[4]全雪峰,黄文海,常梦星.可视化心音诊断教学系统设计[J].工业控制计算机. 2012,25(12)

文章来源: qq912100926.blog.csdn.net，作者：海神之光，版权归原作者所有，如需转载，请联系作者。

原文链接：qq912100926.blog.csdn.net/article/details/115610936