频率采样法设计数字滤波器・实验程序及运行结果实验内容 1:直接频率采样设计法自定义一个能够产生已绝对值的幅度响应、相对 dB 标尺的幅度响应、相位响应和群时延响应。函数名称定义为 freqz_m。函数代码：freqz_m%功能扩展的频响特性函数 freqz_m.mfunction[db,mag,pha,grd,w]=freqz_m(b,a)%H：频响特性的样值向量%w：频响特性的位置向量%db：幅度衰减(dB)%mag：H 的模%pha：H 的相位%grd：H 的群时延%[H,w]=freqz(b,a,1000,'whole');%H=(H(1:1:501))';w=(w(1:1:501))';[H,w]=freqz(b,a,1024,'whole');H=(H(1:1:512))';w=(w(1:1:512))';mag=abs(H);db=20*log10((mag+eps)/max(mag));pha=angle(H);grd=grpdelay(b,a,w);%grpdelay 群时延函数设计分析：k 取 k—匕 N 的最大整数，因此 k—cN=2，故频率米样值 H(k)为cc2 兀 c2 兀 gH(k)—1,1,1,0,,0,1,1g-15 个零-源程序：shzxhchlshiyan8_l.mV'%shzxhchlshiyan8_1%频率采样法的直接设计方法(FIR 低通滤波器)clearall;closeall;clc;clf;M=20;tao=(M-l)/2;%M:采样点数(20 点),tao：群时延(9.5)k=0:M-1;wl=(2*pi/M)*k;%w1 描述各采样点频率值Hrs=[l,l,l,zeros(l,15),l,l];%频域采样值(频率采样样本向量)%Hg(0)=Hg(1)=Hg(2)=Hg(18)=Hg(19)=1%kc=(wc/2n)N,floor(kc)=2Hdr=[l,l,0,0];wdl=[0,0.25,0.25,l];%Hdr 与 wdl 共同确定 1 和 0 出现的频率位置边界%当设置 wd 上[0,0.25,0.25,1]时，是锐截止，无过渡带%当设置 wdl=[0,0.2,0.3,1]时，是设置频率在 0.2n 到 0.3n 之间的过渡带k1=0:floor((M-1)/2);k2=floor((M-1)/2)+1:M-1;%(频率采样位置向量，0-9、10-19)%k2=ceil((M-l)/2):M-l%floor:向小的整数方向取整，ceil：向大的整数方向取整%fix：取小数的整数部分，round:四舍五入取整 angH=[-tao*(2*pi)/M*kl,tao*(2*pi)/M*(M-k2)];%计算-T(2n/N)kc(相位计算)H=Hrs.*exp(j*angH);%求传输函数h=real(ifft(H,M));%求脉冲响应(单位冲激响应)，即系数值%h=abs(ifft(H,M));%%%不能使用绝对值函数[db,mag,pha,grd,w]=freqz_m(h,1);%分子系数即 h,分母系数恒等于 1(FIR 滤波器)%w：频率的位置向量%db:幅度衰减(dB)%mag：H 的模%pha：H 的相位%grd：H 的群时延[Hr,ww,a,L]=Hr_type2(h);%Hr:振幅响应的样值向量%ww:在[0,6n]区间内计算 Hr 的 512 个频率点%a：第 II 类线性相位 FIR 低通滤波器的系数%L：Hr 的阶次%%subplot(2,2,1);plot(wl(1:11)/pi,Hrs(1:11)...
