ms24_2.m

《数字信号处理实验(MATLAB版)》这

%矩形窗设计FIR低通滤波器
wc=0.4*pi;            %输入设计指标
N=64;n=0:N-1;
hd=ideal_lp(wc,N);    %建立理想低通滤波器
windows=(boxcar(N))'; %使用矩形窗，并将列向量变为行向量
b=hd.*windows;        %求FIR系统函数系数
[H,w]=freqz(b,1);     %求解频率特性
dbH=20*log10((abs(H)+eps)/max(abs(H))); %化为分贝值
subplot(2,2,1),stem(n,hd,'filled','k');
axis([0,N,1.1*min(hd),1.1*max(hd)]);title('理想脉冲响应');
xlabel('n');ylabel('hd(n)');
subplot(2,2,2),stem(n,windows,'filled','k');
axis([0,N,0,1.1]);title('窗函数特性');
xlabel('n');ylabel('wd(n)');
subplot(2,2,3),stem(n,b,'filled','k');
axis([0,N,1.1*min(b),1.1*max(b)]);title('实际脉冲响应');
xlabel('n');ylabel('h(n)');
subplot(2,2,4),plot(w/pi,dbH,'k');
axis([0,1,-80,10]);title('幅度频率响应');
xlabel('频率(单位:\pi)');ylabel('H(e^{j\omega})');
set(gca,'XTickMode','manual','XTick',[0,wc/pi,1]);
set(gca,'YTickMode','manual','YTick',[-50,-20,-3,0]);grid