f7_9.m

matlab时频分析配套源码第6章的程序

%样本点数
N=5000;
[normal,Fs,bits] = wavread('normal',N);
[mitralstenosis,Fs,bits] = wavread('mitralstenosis',N);
[aorticinsufficiency,Fs,bits] = wavread('aorticinsufficiency',N);
[pulmonarystenosis,Fs,bits] = wavread('pulmonarystenosis',N);

figure(1);
%正常时心音波形
subplot(2,1,1);
plot(normal);
xlabel('样本序号 n/个');
ylabel('幅值 A/V');
title('正常时心音波形');
axis ([0 N -0.1 0.1]);

figure(2);
%大动脉狭窄时心音波形
subplot(3,1,1);
plot(aorticinsufficiency);
% xlabel('样本序号 n/个');
ylabel('幅值 A/V');
title('大动脉狭窄时心音波形');
axis ([0 N -0.2 0.2]);
%二尖瓣狭窄时心音波形
subplot(3,1,2);
plot(mitralstenosis);
% xlabel('样本序号 n/个');
ylabel('幅值 A/V');
title('二尖瓣狭窄时心音波形');
axis ([0 N -0.1 0.1]);

%肺动脉瓣狭窄时心音波形
subplot(3,1,3);
plot(pulmonarystenosis);
xlabel('样本序号 n/个');
ylabel('幅值 A/V');
title('肺动脉瓣狭窄时心音波形');
axis ([0 N -0.2 0.2]);