基于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为间
