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

✅博主简介：热爱科研的Matlab仿真开发者，修心和技术同步精进，Matlab项目合作可私信。
🍎个人主页：海神之光
🏆代码获取方式：
海神之光Matlab王者学习之路—代码获取方式
⛳️座右铭：行百里者，半于九十。

更多Matlab仿真内容点击👇
Matlab图像处理（进阶版）
路径规划（Matlab）
神经网络预测与分类（Matlab）
优化求解（Matlab）
语音处理（Matlab）
信号处理（Matlab）
车间调度（Matlab）

⛄一、心音诊断系统简介

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

⛄三、运行结果

⛄四、matlab版本及参考文献

1 matlab版本
2014a

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

3 备注
简介此部分摘自互联网，仅供参考，若侵权，联系删除

🍅 仿真咨询
1 各类智能优化算法改进及应用
生产调度、经济调度、装配线调度、充电优化、车间调度、发车优化、水库调度、三维装箱、物流选址、货位优化、公交排班优化、充电桩布局优化、车间布局优化、集装箱船配载优化、水泵组合优化、解医疗资源分配优化、设施布局优化、可视域基站和无人机选址优化

2 机器学习和深度学习方面
卷积神经网络（CNN）、LSTM、支持向量机（SVM）、最小二乘支持向量机（LSSVM）、极限学习机（ELM）、核极限学习机（KELM）、BP、RBF、宽度学习、DBN、RF、RBF、DELM、XGBOOST、TCN实现风电预测、光伏预测、电池寿命预测、辐射源识别、交通流预测、负荷预测、股价预测、PM2.5浓度预测、电池健康状态预测、水体光学参数反演、NLOS信号识别、地铁停车精准预测、变压器故障诊断

3 图像处理方面
图像识别、图像分割、图像检测、图像隐藏、图像配准、图像拼接、图像融合、图像增强、图像压缩感知

4 路径规划方面
旅行商问题（TSP）、车辆路径问题（VRP、MVRP、CVRP、VRPTW等）、无人机三维路径规划、无人机协同、无人机编队、机器人路径规划、栅格地图路径规划、多式联运运输问题、车辆协同无人机路径规划、天线线性阵列分布优化、车间布局优化

5 无人机应用方面
无人机路径规划、无人机控制、无人机编队、无人机协同、无人机任务分配

6 无线传感器定位及布局方面
传感器部署优化、通信协议优化、路由优化、目标定位优化、Dv-Hop定位优化、Leach协议优化、WSN覆盖优化、组播优化、RSSI定位优化

7 信号处理方面
信号识别、信号加密、信号去噪、信号增强、雷达信号处理、信号水印嵌入提取、肌电信号、脑电信号、信号配时优化

8 电力系统方面
微电网优化、无功优化、配电网重构、储能配置

9 元胞自动机方面
交通流 人群疏散 病毒扩散 晶体生长

10 雷达方面
卡尔曼滤波跟踪、航迹关联、航迹融合