f7_6.m

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

%读入语音信号
load gabor;
t=0:337;
%采样频率
fs=1000;
%绘制时域波形
figure(1);
plot(t,gabor,'LineWidth',2);
xlabel('时间 t');
ylabel('幅值 A');
%计算功率谱
nfft=10240;
dsp=abs(fft(gabor,nfft));
%绘制功率谱图
figure(2)
plot((0:nfft/2-1)/nfft*fs,dsp(1:nfft/2),'LineWidth',2);
xlabel('频率 f');
ylabel('功率谱 PSD');
%计算短时傅立叶变换
%%窗函数
h=window(85,'hamming');
[tfr,t,f]=tfrstft(gabor,1:338,256,h);
%%绘制等高线图
figure(3)
contour(t,(0:127)/256*fs,abs(tfr(1:128,:)).^2);
xlabel('时间 t');
ylabel('频率 f');
%%三维图
figure(4)
mesh(t,(0:127)/256*fs,abs(tfr(1:128,:)).^2);
xlabel('时间 t');
ylabel('频率 f');
zlabel('幅值 A');