第7章有限脉冲响应数字滤波器的设计.ppt

第七章 习题 16 （第三版 习题19） 解： 技术指标： 频率：通带截止频率 15kHz, 阻带截止频率 20kHz, 采样频率 幅度：通带 阻带 FIR滤波器：采样等波纹最佳逼近法设计。 [M, fo, mo, w]=remezord(f, m, rip, Fs) hn=remez(M, fo, mo, w) 需要用的MATLAB命令： 参数： close all,clear all Fs=60000; f=[15000, 20000]; m=[1, 0];rip=[0.02, 0.01]; [M, fo, mo, w]=remezord(f, m, rip, Fs) M=M+2;hn=remez(M, fo, mo, w); subplot(3,2,1),stem([0:length(hn)-1],hn,.),title(单位脉冲响应) N=1024;Hw=fft(hn,N); H=20*log10(abs(Hw)/max(Hw)); subplot(3,2,2),plot(0.5*Fs*[0:N/2-1]/(0.5*N),H(1:N/2)) grid on,axis([0,Fs/2,-100,5]),title(滤波器的损耗函数) T=1/Fs;n=0:60;t=n*T; x=cos(2*pi*10000*t)+cos(2*pi*17000*t)+cos(2*pi*25000*t)+cos(2*pi*15000*t); X=fft(x,1024); subplot(3,2,3),plot(t*1000,x),title(输入信号),xlabel(ms) subplot(3,2,5),plot(0.5*Fs*[0:511]/512000,abs(X(1:512))),xlabel(kHz) y=filter(hn,1,x);Y=fft(y,1024) subplot(3,2,4),plot(1000*t,y),title(输出信号),xlabel(ms) subplot(3,2,6),plot(0.5*Fs*[0:511]/512000,abs(Y(1:512))),xlabel(kHz) 滤波器的阶数：M = 20 边界频率：fo = 0 0.5000 0.6667 1.0000 希望逼近的幅度特性：mo = 1 1 0 0 误差加权向量：w = 1 2 运行结果： 调用remez函数设计的低通FIR数字滤波器(20阶)的h(n)及损耗函数曲线 调用remez函数设计的低通FIR数字滤波器(22阶)的h(n)及损耗函数曲线 IIR滤波器：采用模拟椭圆滤波器作为原型，用双线性变换法将模拟滤波器转换为数字滤波器。 需要用的MATLAB命令： [N, wpo]=ellipord(wp, ws, Rp, As); [B, A]=ellip(N, Rp, As, wpo, ); wp=2*atan(pi*15000/Fs)/pi ws=2*atan(pi*20000/Fs)/pi Rp=0.35,As=40 close all,clear all Rp=0.35;As=40;Fs=60000; ws=2*atan(pi*20000/Fs)/pi;wp=2*atan(pi*15000/Fs)/pi; [N, wpo]=ellipord(wp, ws, Rp, As) [B, A]=ellip(N, Rp, As, wpo),[Hk,w]=freqz(B, A); subplot(3,1,1),plot(0.5*Fs*w/pi ,20*log10(abs(Hk)));