a=wavread("jiasiqi.wav"); %將音訊訊號jiasiqi.wav讀入subplot(2,1,1), %分配畫布，一幅圖上共兩個圖，這是第一個plot(a);title("original signal"); %畫出原始訊號，即前面這個音訊訊號的原始波形grid %新增網格線N=256; %設定短時傅立葉變換的長度，同時也是漢明窗的長度h=hamming(N); %設定漢明窗for m=1:N %用漢明窗擷取訊號，長度為N，主要是為了減少截斷引起的柵欄效應等b(m)=a(m)*h(m)endy=20*log(abs(fft(b))) %做傅立葉變換，取其模值，即幅頻特性，然後用分貝（dB）表示subplot(2,1,2) %分配畫布，第二副圖plot(y);title("短時譜"); %畫出短時譜grid %新增網格線
a=wavread("jiasiqi.wav"); %將音訊訊號jiasiqi.wav讀入subplot(2,1,1), %分配畫布，一幅圖上共兩個圖，這是第一個plot(a);title("original signal"); %畫出原始訊號，即前面這個音訊訊號的原始波形grid %新增網格線N=256; %設定短時傅立葉變換的長度，同時也是漢明窗的長度h=hamming(N); %設定漢明窗for m=1:N %用漢明窗擷取訊號，長度為N，主要是為了減少截斷引起的柵欄效應等b(m)=a(m)*h(m)endy=20*log(abs(fft(b))) %做傅立葉變換，取其模值，即幅頻特性，然後用分貝（dB）表示subplot(2,1,2) %分配畫布，第二副圖plot(y);title("短時譜"); %畫出短時譜grid %新增網格線