基于matlab的音乐信号处理

-Matlab音乐信号处理程序*1*close all;clear all;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(原信号频谱);grid on; %减抽样时域j=0;d=3; %减抽样间隔，分别取3和10，代表混叠与非混叠w3=w1(1:100000); %信号截断，提高执行速度for i=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.65b1,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*clear all;close all;clc;w,fs,b=wavread(宽恕.wav);w1=w(:,1);l=length(w1);ww=2/l*(0:l-1); %三余弦混合噪声n=0:length(w1)-1;x=0.05*(cos(2*pi*3000*n/fs)+cos(2*pi*5000*n/fs)+cos(2*pi*8000*n/fs);fx=fft(x);subplot(3,2,1);plot(ww,abs(fx);title(噪声频谱);r=w1+x;fr=fft(r);subplot(3,2,2);plot(r);title(加噪后音乐信号的音谱);subplot(3,2,3);plot(ww,abs(fr);title(加噪后音乐信号的频谱);%sound(r,fs);%滤波N,wc=buttord(0.05,0.11,1,30);b1,a=butter(N,wc);h,wi=freqz(b1,a); %用iir对信号滤波R=filter(b1,a,r);fr=fft(R);subplot(3,2,4);plot(R);title(滤波去噪后信号的音谱);subplot(3,2,5);plot(ww,abs(fr);title(滤波去噪后信号频谱);sound(R,fs);*4*clear all;close all;clc;w,fs,b=wavread(宽恕.wav);w1=w(:,1);l=length(w1);ww=2/l*(0:l-1); %白噪声d=rand(size(w1)-0.5;df=fft(d);subplot(3,2,1);plot(ww,abs(df);title(白噪声频谱); g=w1+d;fg=fft(g);subplot(3,2,2);plot(g);title(加白噪声后音乐信号的音谱);subplot(3,2,3);plot(ww,abs(fg);title(加白噪声后音乐信号的频谱);sound(g,fs);%滤波N,wc=buttord(0.08,0.15,2,20);b1,a=butter(N,wc);h,wi=freqz(b1,a); %用iir对信号滤波G=filter(b1,a,g);fg1=fft(G);subplot(3,2,4);plot(G);title(滤波去噪后信号的音谱);subplot(3,2,5);plot(ww,abs(fg1);title(滤波去噪后信号频谱);%sound(G,fs);*另存