fb=10;  fc=13;                   % 设置滤波器截止频率  fs=100;   wb=2*pi*fb/fs;     ws=2*pi*fc/fs;    wc=0.5*(wb+ws); tr_width=ws-wb;               % 过渡带宽 M=ceil(1.8*pi/tr_width);  hd=ideal_lp(wc,M);             % 产生一理想低通滤波器  w_box=(boxcar(M))';               % 矩形窗窗函数  h=hd.*w_box; h=hd.*w_box;h=hd.*w_box;  [db,mag,pha,grd,w]=freqz_m(h,[1]);  n=[0:1:M-1];  subplot(2,2,1); stem(n,hd);                   % 理想脉冲响应  xlabel('n'); ylabel('hd(n)');  title('Ideal Impulse Response');  subplot(2,2,2); stem(n,w_box);                  % 矩形窗  xlabel ('n'); ylabel('w(n)');  title('Boxcar Window'); subplot(2,2,3); stem(n,h);                  % 实际脉冲响应  xlabel('n'); ylabel('h(n)');  title('Actual Impulse Response');  subplot(2,2,4); plot(w*fs/(2*pi),db);      % 幅度响应（dB）  axis([0 40 -50 0]);  xlabel('Frequency (Hz)'); ylabel('Decibels'); title('Magnitude Response in dB');  grid;