心音信号处理分析（附matlab代码）

在matlab仿真平台上对某一心音信号进行采样、小波分解、消噪处理、信号重构等信号处理。

clc
clear
[ufs]=wavread（‘f:\matlab\work\11.wav‘）;
y=resample（u115）;
fs=fs/15;
figure
plot（u）;
title（‘源心音信号‘）;
figure
plot（y）;
title（‘采样后的心音信号‘）;
u_beg=4000;
u_end=8000;
xn=y（u_beg:u_end）;
figure;
plot（xn）;
title（‘截取一段心音信号‘）;
ls=length（xn）;
w_name=‘db6‘;

%对心音信号进行小波分解
[cl]=wavedec（xn5w_name）;%5层小波分解（层：尺度）
ca5=appcoef（clw_name5）;% 提取一维小波变换低频系数（第五层的逼近系数）
cd5=detcoef（cl5）;%提取一维小波变换高频系数 （细节系数）
cd4=detcoef（cl4）;
cd3=detcoef（cl3）;
cd2=detcoef（cl2）;
cd1=detcoef（cl1）;
figure;
subplot（611）;plot（ca5）;title（‘低频系数‘）;
subplot（612）;plot（cd5）;title（‘第五层高频系数‘）;
subplot（613）;plot（cd4）;title（‘第四层高频系数‘）;
subplot（614）;plot（cd3）;title（‘第三层高频系数‘）;
subplot（615）;plot（cd2）;title（‘第二层高频系数‘）;
subplot（616）;xlabel（‘小波的各层频率系数‘）;plot（cd1）;title（‘第一层高频系数‘）;




%消除噪声处理
[thrsorhkeepapp]=ddencmp（‘den‘‘wv‘xn）;
s1=wdencmp（‘gbl‘clw_name3thrsorhkeepapp）;
figure;
subplot（212）;
plot（s1）;
title（‘小波去噪后的信号‘）;
subplot（211）;
plot（xn）;
title（‘含有噪声的信号‘）;

%去除噪声，重构信号对一维小波系数进行单支重构（重构之后得到的是信号，单支重构的目的在于获得原信号中某频率段的分量信号。分量信号的长度（点数）、采样
%频率均与原信号相同。将单支重构得到的各个分量信号直接按对应点求和，可恢复原信号。）

%重构第5层逼近信号（重构各分解尺度（层）上的信号）
a5=wrcoef（‘a‘clw_name5）;%是对尺度5上的低频部分进行重构，a是低频，d是高频，cl是wavedec所得
%重构第1~5层细节信号
d5=wrcoef（‘d‘clw_name5）;
d4=wrcoef（‘d‘clw_name4）;
d3=wrcoef（‘d‘clw_name3）;
d2=wrcoef（‘d‘clw_name2）;
d1=wrcoef（‘d‘clw_name1）;


figure;
subplot（211）;
plot（xn）;
title（‘噪声污染的信号‘）;
 
subplot（212）;
y=a5+d4;
plot（y）;
title（‘小波分解重构去噪后的信号‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件     275962  2012-10-13 11:30  心音信号处理分析\11.wav

     文件       1848  2012-12-15 14:54  心音信号处理分析\xinhaochuli.m

     文件      89600  2013-05-30 15:31  心音信号处理分析\心音信号处理分析.doc

     目录          0  2013-05-30 15:31  心音信号处理分析

----------- ---------  ---------- -----  ----

               367410                    4