Fir_design2.zip_明汉窗滤波_汉明窗 matlab_汉明窗滤波_矩形滤波_矩形窗 MATLAB

N=21;wc=pi/4; %理想低通滤波器参数 n=0:N-1;r=(N-1)/2; hdn=sin(wc*(n-r))/pi./(n-r); %计算理想低通单位脉冲响应 if rem(N,2)~=0 hdn(r+1)=wc/pi; end %N为奇数时，处理n=r点的0/0型 %wn1=boxcar(N); %矩形窗 %hn1=hdn.*wn1'; %加窗 w=2*[0:511]/512; hn1=fir1(N-1,wc/pi,boxcar(N)); hw1=fft(hn1,512); figure(1) subplot(1,2,1); %作图 stem(n,hn1,'.');grid; xlabel('n'); ylabel('h(n)'); title('矩形窗设计的h(n) '); subplot(1,2,2); plot(w,20*log10(abs(hw1)));grid; xlabel('w/pi'); ylabel('Magnitude(dB)'); title('幅频特性（dB）'); %wn2=hanning(N); %hanning窗 %hn2=hdn.*wn2'; %加窗 hn2=fir1(N-1,wc/pi,hanning(N)); hw2=fft(hn2,512); figure(2) subplot(1,2,1); %作图 stem(n,hn2,'.');grid; xlabel('n'); ylabel('h(n)'); title('hanning窗设计的h(n) '); subplot(1,2,2); plot(w,20*log10(abs(hw2)));grid; xlabel('w/pi'); ylabel('Magnitude(dB)'); title('幅频特性（dB）'); %wn3=hamming(N); %hamming窗 %hn3=hdn.*wn3'; %加窗 hn3=fir1(N-1,wc/pi,hamming(N)); hw3=fft(hn3,512); figure(3) subplot(1,2,1); %作图 stem(n,hn2,'.');grid; xlabel('n'); ylabel('h(n)'); title('hamming窗设计的h(n) '); subplot(1,2,2); plot(w,20*log10(abs(hw3)));grid; xlabel('w/pi'); ylabel('Magnitude(dB)'); title('幅频特性（dB）'); %wn4=blackman(N); %Blackman窗 %hn4=hdn.*wn4'; %加窗 hn4=fir1(N-1,wc/pi,blackman(N)); hw4=fft(hn4,512); figure(4) subplot(1,2,1); %作图 stem(n,hn2,'.');grid; xlabel('n'); ylabel('h(n)'); title('Blackman窗设计的h(n) '); subplot(1,2,2); plot(w,20*log10(abs(hw4)));grid; xlabel('w/pi'); ylabel('Magnitude(dB)'); title('幅频特性（dB）');