基于matlab的海明窗带阻滤波器的设计

已经答辩通过的程序。保证可用。 分别是两个调用函数与一个主函数。 含验证程序。

function [dbmagphagrdw] = freqz_m（ba）;
% Modified version of freqz subroutine
% ------------------------------------
% [dbmagphagrdw] = freqz_m（ba）;
% db = Relative magnitude in dB computed over 0 to pi radians
% mag = absolute magnitude computed over 0 to pi radians 
% pha = Phase response in radians over 0 to pi radians
% grd = Group delay over 0 to pi radians
% w = 501 frequency samples between 0 to pi radians
% b = numerator polynomial of H（z） （for FIR: b=h）
% a = denominator polynomial of H（z） （for FIR: a=[1]）
%
[Hw] = freqz（ba1000‘whole‘）;
H = （H（1:1:501））‘; w = （w（1:1:501））‘;
mag = abs（H）;
db = 20*log10（（mag+eps）/max（mag））;
pha = angle（H）;
% pha = unwrap（angle（H））;
grd = grpdelay（baw）;
% grd = diff（pha）;
% grd = [grd（1） grd];
% grd = [0 grd（1:1:500）; grd; grd（2:1:501） 0];
% grd = median（grd）*500/pi;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        855  2010-11-30 22:45  freqz_m2.m

     文件        129  2010-11-30 15:37  ideal_bs.m

     文件        335  2010-11-30 17:40  jet.asv

     文件        600  2010-11-30 18:07  jet.m

     文件       1486  2010-11-30 23:04  jet1.m

     文件        882  2010-12-01 09:44  jet2.m

----------- ---------  ---------- -----  ----

                 4287                    6