音乐信号的处理matlab.docxVIP

1、音乐信号的音谱和频谱观察 [y,fs,bit]=wavread(2); sound(y,fs); plot(y); y1=y(:,1); M=512; fh1=fft(y1,M); w=2/M*[0:M-1]; figure;plot(w,abs(fh1)); 0.25 0.2 0.1S 0.1 0 05 Q . I I I I 」 1 . 5 0 0.2 0.4 C 6 0 8 1 1 2 1 4 1 S 1.8 2 2、音乐信号的抽取(减抽样) [y,fs,bit]=wavread(2); sound(y,fs); figure; subplot(1,2,1); plot(y); title(原音乐信号波形); y1=y(:,1); M=512; fh1=fft(y1,M); w=2/M*[0:M-1]; subplot(1,2,2); plot(w,fh1); title(原音乐信号频谱); grid on; j=0; d=10; %减抽样间隔 取3和10代表混叠与非混叠 y3=y1(1:100000); for i=1:d:length(y3); j=j+1; yd(j)=y1(i); end sound(yd,fs/d); figure; ■。卫00 5 ■。卫0 0 5 subplot(1,2,1); plot(yd); title(减抽样时域); fh2=fft(yd); l1=length(fh2); w2=2/l1*(0:l1-1); subplot(1,2,2); plot(w2,abs(fh2)); title(减抽样频域); 原音乐信号频谱 0.15 0.1 0 05 0 -0 05 -0J -0.16 2500 r - 2000 - 1600 - 1000 - 500 - L 3、 音乐信号的AM调制和巴特沃斯IIR滤波器滤波 [w,fs,bit]=wavread(2); w1=w(:,1); l=length(w1); ww=2/l*(0:l-1); n=0:length(w1)-1; t1=w1.*cos(0.6*pi*n); T1=fft(t1); plot(t1); title(AM低/高频调制后波形 figure; plot(ww,abs(T1)); title(AM低/高频调制后频谱 sound(t1,fs); t2=t1.*cos(0.6*pi*n); T2=fft(t2); figure; plot(t2); title(AM低/高频解调后波形 figure; plot(ww,abs(T2)); title(AM低/高频解调后频谱 sound(t2,fs); [N,wc]=buttord(0.2,0.3,2,15); )； )； )； )； [b1,a]=butter(N,wc); [h,wi]=freqz(b1,a); figure; plot(wi/pi,abs(h)); title(iir滤波器频谱); %用iir对信号滤波 y=filter(b1,a,t2); Y=fft(y); figure; plot（y）; title（低/高频解调滤波后时域）; figure; plot(ww,abs(Y)); title（低/高频解调滤波后的频域）; sound(y,fs); AM低/高频调制后拔节 10 80.604 1 0 8 0.6 04 0.2 0 -0.2 -0.4 -0.6 -0.8 0 0 5 1 2 2.6 3 x 10  用窗函数法设计 FIR滤波器(矩形窗和布莱克曼窗) [w,fs,bit]=wavread(2); w1=w(:,1); n=0:length(w1)-1; 1.4 1.2 1 0.8 0 6 04 0.2 0 t1=w1.*cos(0.6*pi*n); j1=t1.*cos(0.6*pi*n); %用FIR对信号同步解调 N=55;wc=pi/4;% boxcar N 取 9hd=ideal(N,wc); w=blackman(N);%blackman 和 boxcarh=hd.*w; M=512; fh=fft(h,M); ww=2/M*(0:M-1); plot(ww,abs(fh)); title(矩形窗fir滤波器频率响应); y=conv(j1,h); Y=fft(y,M); figureplot(y); title(矩形窗fir滤波后时域); figureplot(ww,abs(Y)); title(矩形窗fir滤波后频域,); sound(y,fs); 矩形商fi褫波器频率响应 三余弦和白噪声的实验 [y,fs,bit]=wavread(2); yi=y(：,i); f=fft(y1); l=length(f); w