%讀入語音信號
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');