基于matlab的乐信号处理和分析.docVIP

基于matlab的音乐信号处理和分析 1音乐信号的音谱和频谱的观察 使用windows下的录音机录制一段音乐信号或者采用其他软件截取一段音乐信号（要求：时间不超过5s，文件格式为WAV） 使用wavread语句读取音乐信号获取抽样率；（注意：读取的信号是双声道信号，即为双列向量，需要分列处理）； 输出音乐信号的波形和频谱，观察现象； 使用sound语句播放音乐信号，注意不同抽样率下的音调变化，解释现象 程序： clear all;close all;clc [y,fs,bit]=wavread(E:\music\11); %读取音乐信号 size(y) %看音乐信号是双列还是单列 y1=y(:,1); %取单列 fs=fs %获取音乐信号的抽样率 N=length(y1); %音乐信号的长度 Fy1=fft(y1,N); %对信号做快速傅里叶变换 w=2/N*[0:N-1]; figure %画音乐信号的波形和频谱 subplot(2,1,1);plot(y1); grid on title(音乐信号的波形) xlabel(time/s) ylabel(Magnitude) subplot(2,1,2);plot(w,abs(Fy1)); grid on title(音乐信号的频谱) xlabel(Frequency/pi) ylabel(Magnitude) sound(y1,fs) %以抽样率fs播放音乐信号 sound(y1,2*fs) %以抽样率2fs播放音乐信号 sound(y1,1/2*fs) %以抽样率1/2fs播放音乐信号 程序运行结果如下： 音乐信号的抽样率fs=44100 音乐信号的波形和频谱如下图： 抽样率较高时，听到的音乐音调较高，速度较快；抽样率较低时，听到的音乐音调低，速度慢。由此可知，音乐信号的音调和抽样率有关。 2音乐信号的抽取（减抽样） 观察音乐信号频率的上限，选择适当的抽取间隔对信号进行减抽样（给出两种抽取间隔，代表混叠与非混叠）； 输出减抽样音乐信号的波形和频谱，观察现象，给出理论解释； 播放减抽样音乐信号，注意抽样率的改变，比较不同抽取间隔下的声音，解释现象。 音乐信号的频率上限为0.2pi，这里分别采用3和60作为非混叠与混叠的抽样间隔进行抽样，程序如下： clear all;close all;clc [y,fs,bit]=wavread(E:\music\22);%读取音乐信号 y1=y(:,1); %取单列 x1=y1(1:180000); %截取音乐信号 d=3;j=0; for i=1:d:length(x1) j=j+1; x2(j)=x1(i); end %对音乐信号以3为间隔做减抽样 N1=length(x1); %截断的音乐信号的长度 N2=length(x2); %以3为间隔减抽样后音乐信号的长度 Fx1=fft(x1,N1); %对截断信号做快速傅里叶变换 Fx2=fft(x2,N2); %对减抽样后的信号做快速傅里叶变换 w1=2/N1*[0:N1-1]; w2=2/N2*[0:N2-1]; figure %画出截断后的信号的波形和频谱 subplot(2,1,1);plot(x1); grid on title(截断后的波形) xlabel(time/s) ylabel(Magnitude) subplot(2,1,2);plot(w1,abs(Fx1)); grid on title(截断信号的频谱) xlabel(Frequency/pi) ylabel(Magnitude) figure %画出以3为间隔减抽样后的信号的波形和频谱 subplot(2,1,1);plot(x2); grid on title(对截断信号以3为间隔减抽样后的波形) xlabel(time/s) ylabel(Magnitude) subplot(2,1,2);plot(w2,abs(Fx2)); grid on title(对截断信号以3为间隔减抽样后的频谱) xlabel(Frequency/pi) ylab