matlab处理音乐数字信号.docVIP

数字信号处理 学院：物理工程学院 姓名：金义飞 读取音乐信号并观察波形与频谱 程序 [w,fs,bit]=wavread(你是我的眼.wav); %音频信号w。采样率fs 精度bit wav=(w(:,1)); %取第一列信号 sound(w,fs) %听取抽样后的wav音乐 sound(w,2*fs) %听取抽样后的wav音乐 sound(w,fs/2) %听取抽样后的wav音乐 figure; %创建一个窗口 subplot(2,1,1);plot(wav); %画出时域图 fwav=fft(wav); %对音频信号做傅里叶变换 lwav=round(length(fwav)/2); %音频信号正半轴长度 nwav=[0:lwav-1]; wwav=nwav/(lwav);f=wwav/2*fs; %x轴坐标 实际频率 subplot(2,1,2); plot(wwav,abs(fwav(1:lwav))); %画出音乐频域图 图像 现象说明： 从音谱即时域谱可以看出音乐信号随时间的幅度变化；从频谱可看到音乐信号主要分布在低频段，高频成分较少，在0.5pi 以后几乎无音乐信号 音乐信号的抽取（减抽样） 三倍减抽样 1，程序 clear all;close all;clc %清除存储空间，关闭所有窗口，清除命令 [w,fs,bit]=wavread(你是我的眼.wav) wav=(w(:,1)); fwav=fft(wav); lwav=round(length(fwav)/2); j=0;d1=3; for i=1:d1:lwav; j=j+1; dwav1(j)=wav(i); end %本段每隔三个点取一个点，相当于三倍减抽样 sound(dwav1,fs/d1); figure; subplot(2,1,1);plot(dwav1); fwav=fft(dwav1); lwav=round(length(fwav)/2); nwav=[0:lwav-1]; wwav=nwav/(lwav);f=wwav/fs; subplot(2,1,2); plot(wwav,abs(fwav(1:lwav))); 2图像 20倍减抽样 1，程序 clear all;close all;clc [w,fs,bit]=wavread(你是我的眼.wav) wav=(w(:,1)); fwav=fft(wav); lwav=round(length(fwav)/2); j=0;d1=20; for i=1:d1:lwav; j=j+1; dwav1(j)=wav(i); end sound(dwav1,fs/d1); figure; subplot(2,1,1);plot(dwav1); fwav=fft(dwav1); lwav=round(length(fwav)/2); nwav=[0:lwav-1]; wwav=nwav/(lwav);f=wwav/fs; subplot(2,1,2); plot(wwav,abs(fwav(1:lwav))); 分别选取抽样间隔为3倍频和20倍频进行减抽样，后者产生混叠，前者不混叠，当选取3倍频减抽样时，频率样间隔变大但未超过最小抽样间隔，尽管频谱有所延伸，但并不产生混叠，即减抽样后频谱上限未超过折叠频率；而选取20倍频减抽样时，由于时域抽样间隔变大超过最