wav信號的波形分析

武漢大學教學實驗報告電子信息學院電子信息科學與技術(shù)專業(yè)201311wav信號的波形分析指導教師卜方玲姓名 年級2010級學號成一、預習部分實驗?zāi)康膶嶒灮驹碇饕獌x器設(shè)備（含必要的元器件、工具）實驗?zāi)康膰L試利用短時傅里葉變換分析非平穩(wěn)信號的頻譜變化。熟悉MATLABwavreadwavrecordwavplayfftifft用。實驗基本原理（，從頻率域可以很清楚地了解它們由哪些二、實驗操作部分實驗數(shù)據(jù)、表格及數(shù)據(jù)處理實驗操作過程（可用圖表示）實驗結(jié)論將Windowsding.wavfftplot觀察主要的正弦分量；參考代碼：[y,fs]=wavread('Vuvuzela.wav');%取其中的一個聲道,譬如說，右聲道yr=y(:,2);%截取前1024個點yr=yr(1:1024);2048、1024FFT，觀察產(chǎn)生的不同頻譜；參考代碼：%求取幅度譜并顯示，首先是2048點的YR2048=fft(yr,2048);figure('numbertitle','off','name','2048點FFT');plot(linspace(-pi,pi,2048),fftshift(abs(YR2048)));%求取1024點的YR1024=fft(yr,1024);figure('numbertitle','off','name','1024點FFT');plot(linspace(-pi,pi,1024),fftshift(abs(YR1024)));運行結(jié)果圖：ifft始信號的波形；參考代碼：%反變換，結(jié)果應(yīng)與原始信號相同yr1024=real(ifft(YR1024));figure('numbertitle','off','name','1024點fft的反變換');subplot(2,1,1);plot(yr1024)title('1024點fft的反變換')subplot(2,1,2);plot(yr);title(運行結(jié)果圖：從頻譜中找到幅度最大的正弦分量，構(gòu)造一個同樣幅度的正弦信號，將其波形與原始信號比較，并且試聽一下。從頻譜中找到幅度最大的正弦分量，構(gòu)造一個同樣幅度的正弦信號，將其波形與原始信號比較，并且試聽一下。參考代碼：%重構(gòu)%尋找幅度最大的正弦分量[maxpeak,peaki]=max(abs(YR1024(1:512)));MAXSIN=zeros(1,1024);MAXSIN(peaki)=maxpeak;MAXSIN(1026-peaki)=maxpeak;maxsin=ifft(MAXSIN);figure('numbertitle','off','name','用最大正弦分量重構(gòu)信號');subplot(2,1,1)plot(maxsin);title(subplot(2,1,2);plot(yr1024);title('原始信號');%試聽wavplay(yr1024);wavplay(maxsin);運行結(jié)果圖：實驗結(jié)論：2048FFT1024FFTIFFT似。三、實驗效