wav信号的波形分析

武汉大学教学实验报告电子信息学院电子信息科学与技术专业201311wav信号的波形分析指导教师卜方玲姓名 年级2010级学号成一、预习部分实验目的实验基本原理主要仪器设备（含必要的元器件、工具）实验目的尝试利用短时傅里叶变换分析非平稳信号的频谱变化。熟悉MATLABwavreadwavrecordwavplayfftifft用。实验基本原理（，从频率域可以很清楚地了解它们由哪些二、实验操作部分实验数据、表格及数据处理实验操作过程（可用图表示）实验结论将Windowsding.wavfftplot观察主要的正弦分量；参考代码：[y,fs]=wavread('Vuvuzela.wav');%取其中的一个声道,譬如说，右声道yr=y(:,2);%截取前1024个点yr=yr(1:1024);2048、1024FFT，观察产生的不同频谱；参考代码：%求取幅度谱并显示，首先是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)));运行结果图：ifft始信号的波形；参考代码：%反变换，结果应与原始信号相同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(运行结果图：从频谱中找到幅度最大的正弦分量，构造一个同样幅度的正弦信号，将其波形与原始信号比较，并且试听一下。从频谱中找到幅度最大的正弦分量，构造一个同样幅度的正弦信号，将其波形与原始信号比较，并且试听一下。参考代码：%重构%寻找幅度最大的正弦分量[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','用最大正弦分量重构信号');subplot(2,1,1)plot(maxsin);title(subplot(2,1,2);plot(yr1024);title('原始信号');%试听wavplay(yr1024);wavplay(maxsin);运行结果图：实验结论：2048FFT1024FFTIFFT似。三、实验效