f7_10.m

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

%样本点数
N=4096;
%正常时心音信号
[x,fxs,bits] = wavread('normal',N);
%肺动脉瓣狭窄时心音信号
[y,fys,bits] = wavread('pulmonarystenosis',N);

%设置参数
%时域窗函数
g=window(9,'Kaiser');
%频域窗函数
h=window(27,'Kaiser');
t=1:N;

%计算Choi-Willams分布
%%正常状态
x=hilbert(x);
x=x-mean(x);
[tfr,t,f]=tfrcw(x,t,128,g,h,10);
%%显示:三维图
figure(1)
mesh(t,f*fxs/2,abs(tfr));
xlabel('样本序号 n/个');
ylabel('频率 f/Hz');
zlabel('幅值 A/V');
%%等高线图
figure(2)
contour(t,f*fxs/2,abs(tfr));
xlabel('样本序号 n/个');
ylabel('频率 f/Hz');

%%异常状态
y=hilbert(y);
y=y-mean(y);
[tfr,t,f]=tfrcw(y,t,128,g,h,10);

%%显示:三维图
figure(3)
mesh(t,f*fys/2,abs(tfr));
xlabel('样本序号 n/个');
ylabel('频率 f/Hz');
zlabel('幅值 A/V');

%%等高线图
figure(4)
contour(t,f*fys/2,abs(tfr));
xlabel('样本序号 n/个');
ylabel('频率 f/Hz');