DSP课程设计-基于MATLAB的音乐信号处理和分析.docVIP

提供各专业全套毕业设计，球球 153893706 《数字信号处理》 课程设计报告 题目: 基于MATLAB的音乐信号处理和分析 专业：测控技术与仪器 年级：2011级 姓名： 学号： 基于MATLAB的音乐信号处理和分析 一、读取音乐信号并观察波形与频谱 [w,fs,bit]=wavread(爸爸去哪儿.wav); wav=(w(:,1)); sound(w,fs) 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; subplot(2,1,2); plot(wwav,abs(fwav(1:lwav))); 以二倍的抽样率听声音信号时，音乐播放的特别快，像被压缩了，播放的时间比原信号短。 以二分之一的抽样率听声音信号时，音乐播放的特别慢，像被拉长了，播放的时间比原信号长。 原信号频谱截止频率为0.5*pi 。 二、（1）音乐信号的5倍减抽样 [w,fs,bit]=wavread(爸爸去哪儿.wav); wav=(w(:,1)); fwav=fft(wav); lwav=round(length(fwav)/2); j=0;d1=5; for i=1:d1:lwav j=j+1; dwav1(j)=wav(i); end figure; subplot(2,1,1); plot(dwav1); fdwav1=fft(dwav1); sound(dwav1,fs/d1) lwav1=round(length(fdwav1)/2); nwav=[0:lwav1-1]; wwav=nwav/(lwav1);f=wwav/2*fs; subplot(2,1,2); plot(wwav,abs(fdwav1(1:lwav1)));%减抽样音乐频域 （2）音乐信号的10倍减抽样（有混频） [w,fs,bit]=wavread(爸爸去哪儿.wav); wav=(w(:,1)); fwav=fft(wav); lwav=round(length(fwav)/2); j=0;d1=10; for i=1:d1:lwav j=j+1; dwav1(j)=wav(i); end figure; subplot(2,1,1); plot(dwav1); fdwav1=fft(dwav1); sound(dwav1,fs/d1) lwav1=round(length(fdwav1)/2); nwav=[0:lwav1-1]; wwav=nwav/(lwav1);f=wwav/2*fs; subplot(2,1,2); plot(wwav,abs(fdwav1(1:lwav1))); 原信号频谱截止频率为0.5*pi。 减抽样后的音乐信号听起来变得尖锐，有失真。 抽样率随着抽样间隔的增大而逐渐变小，声音越来越失真，音调变得急促，而尖锐，信号产生混叠 音乐信号的AM调制 （1）调制频率为pi/2 [w,fs,bit]=wavread(爸爸去哪儿.wav); wav=(w(:,1)); lam=length(wav); amw=pi/2; cwav=cos(amw*[0:lam-1]); amwav=wav.*cwav; sound(amwav,fs) figure; subplot(2,1,1);plot(amwav); title(爸爸去哪儿调制后) famwav=fft(amwav); lwav=round(length(famwav)/2); nwav=[0:lwav-1]; wwav=nwav/(lwav); subplot(2,1,2); plot(wwav,abs(famwav(1:lwav))); (2)调制频率为2*pi/3 [w,fs,bit]=wavread(爸爸去哪儿.wav); wav=(w(:,1)); lam=length(wav); amw=2*pi/3; cwav=cos(amw*[0:lam-1]); amwav=wav.*cwav; sound(amwav,fs) figure; subplot(2,1,1);plot(amwav); title(爸爸去哪儿调制后) famwav=fft(amwav); lwav=round(length(famwav)/2); nwav=[0:lwav-1]; wwav=nwav/(lwav); subplot(2,1,2); plot(wwav,abs(famwav(1:lwav))); 原信号的调制相当于频谱搬移，左移一个右移一个，当调制频率（余弦频率）小于