基于matlab的音乐信号处理.doc

最新范本,供参考！最新范本,供参考！最新范本,供参考！Matlab音乐信号处理——程序******************1********************************************closeall;clearall;clc;%信号音谱[w,fs,b]=wavread('宽恕.wav');w1=w(:,2);subplot(2,2,1);plot(w1);title('原信号音谱');sound(w1,fs);%分别以2fs,fs,fs/2的频率播放声音%信号频谱f=fft(w1);l=length(f);ww=2/l*(0:l-1);subplot(2,2,2);plot(ww,abs(f));title('原信号频谱');gridon;%减抽样时域j=0;d=3;%减抽样间隔，分别取3和10，代表混叠与非混叠w3=w1(1:100000);%信号截断，提高执行速度fori=1:d:length(w3);j=j+1;dw(j)=w1(i);%抽样出的点组成新的抽样序列end%sound(dw,fs/d);%以d为间隔抽样，所以抽样率为fs/dsubplot(2,2,3);plot(dw);title('减抽样时域');%减抽样频域f1=fft(dw);l1=length(f1);w2=2/l1*(0:l1-1);subplot(2,2,4);plot(w2,abs(f1));title('减抽样频域');**********************2*****************************************************注：该程序分为3个部分，以“%%%%%%%%%”为标记，在执行第一部分时，将后面内容全部作为注释屏蔽掉（按ctrl+r），在分别执行IIR滤波及FIR滤波相关的内容时，只将第一部分绘图语句作为注释，其余保留。%调制[w,fs,b]=wavread('宽恕.wav');w1=w(:,1);l=length(w1);ww=2/l*(0:l-1);n=0:length(w1)-1;t1=w1.*cos(0.6*pi*n)';%低频0.4，高频0.6T1=fft(t1);subplot(2,2,1);plot(t1);title('AM低/高频调制后时域');subplot(2,2,2);plot(ww,abs(T1));title('AM低/高频调制后频域');sound(t1,fs);%解调j1=t1.*cos(0.6*pi*n)';%0.4/0.6J1=fft(j1);subplot(2,2,3);plot(j1);title('AM低/高频解调后时域');subplot(2,2,4);plot(ww,abs(J1));title('AM低/高频解调后频域');sound(j1,fs);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%用IIR滤波器[N,wc]=buttord(0.25,0.65,2,20);%wc=0.25,wst=0.5/0.65[b1,a]=butter(N,wc);[h,wi]=freqz(b1,a);subplot(2,2,1);plot(wi/pi,abs(h));title('iir滤波器频谱');%用iir对信号滤波y=filter(b1,a,j1);Y=fft(y);subplot(2,2,2);plot(y);title('低/高频解调滤波后时域');subplot(2,2,3);plot(ww,abs(Y));title('低/高频解调滤波后的频域');sound(y,fs);%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%fir滤波Wc=pi/4;nn=33;%滤波器阶数wfr=blackman(nn);%用boxcar、blackman窗替换hd=ideal_filter(Wc,nn);h=hd.*wfr';fh=fft(h,l);fh_db=20*log10((abs(fh)+eps)/max(abs(fh)));subplot(2,2,1);plot(ww,fh_db);title('fir滤波器频率响应');%用fir滤波器对信号滤波y1=conv(j1,h);Y1=fft(y1,l);subplot(2,2,2);plot(y1);title('滤波后时域');subplot(2,2,3);plot(ww,abs(Y1));title('滤波后频域');sound(y1,fs);*****************3************************************