wav信號(hào)的波形分析

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