数字信号处理课程设计 基于 MATLAB 的音乐信号处理和分析

H,W]=freqz(b,a);Рsubplot(3,1,3);Рplot(W,abs(H));Рtitle('IIR滤波器频谱');Рxlabel('w/pi');Рylabel('幅度');Рgrid on;Рm=filter(b,a,y4);Рsound(m,fs);Рlm=length(m);%滤波后信号长度Рtm=(0:lm-1)/fs;%设置横坐标Рwm=2/lm*[0:lm-1];РM=fft(m);Рfigure(3)Рsubplot(2,1,1);plot(tm,m);Рgrid on;Рtitle('IIR滤波器滤波后波形');Рxlabel('t');Рylabel('幅度');Рsubplot(2,1,2);plot(wm,abs(M));Рtitle('IIR滤波器滤波后频谱');Рxlabel('w/pi');Рylabel('幅度');Р%矩形窗和布莱克曼窗РN=33;wc=0.3*pi;%基于经验的指标,其中N为理想低通滤波器阶数,wc为截止频率Рhd=ideal(N,wc);%调用理想低通滤波器函数Рw1=boxcar(N);%产生各种窗函数Рw2=blackman(N);Рh1=hd.*w1';%加窗设计各种FIR滤波器Рh2=hd.*w2';Рth1=(0:32)/fs;Рth2=(0:32)/fs;РM=21184;Рfh1=fft(h1,M);%矩形窗频谱函数Рw=2/M*[0:M-1];Рfh2=fft(h2,M);%布莱克曼窗频谱函数Рfigure(4)Рsubplot(2,2,1);plot(th1,h1)Рtitle('矩形窗时域');Рsubplot(2,2,2);plot(w,abs(fh1));Рtitle('矩形窗频域');Рsubplot(2,2,3);plot(th2,h2);Рtitle('布莱克曼窗时域');