等波纹最佳逼近法设计滤波器.docVIP

程序 %等波纹最佳逼近法设计滤波器 xt=signal_add_noise; Fs=1000; fp=110;fs=130;Rp=0.2;Rs=48; fb=[fp,fs];m=[1,0]; % 确定remezord函数所需参数f,m,dev dev=[(10^(Rp/20)-1)/(10^(Rp/20)+1),10^(-Rs/20)]; [n,fo,mo,W]=remezord(fb,m,dev,Fs); % 确定remez函数所需参数 hn=remez(n,fo,mo,W); % 调用remez函数进行设计低通fir滤波器的h(n) [H,f]=freqz(hn,1,1024,Fs); %求滤波器的频率响应 subplot(411);stem(hn);title((a) h(n)波形) axis([0,N-1,min(hn),max(hn)]);xlabel(n);ylabel(h(n)); subplot(412);plot(f,20*log10(abs(H)/max(abs(H))));grid title((b)幅度响应);xlabel(f(hz));ylabel(分贝数); %利用设计的滤波器对信号去噪，画出去噪后的波形和频谱 M=1000;Fs=1000;T=