资源描述：
《滤波器设计实例》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、%凯泽窗设计示例1：[n,Wn,beta,ftype]=kaiserord([7131723],[010],[0.010.010.01],100);%得出滤波器的阶数n=38，beta=3.4fc1=10;fc2=20;fs=100;w1=2*fc1/fs;w2=2*fc2/fs;%将模拟滤波器的技术指标转换为数字滤波器的技术指标window=kaiser(n+1,beta);%使用kaiser窗函数b=fir1(n,[w1w2],window);%使用标准频率响应的加窗设计函数fir1freqz(b,1,512);%数字
2、滤波器频率响应t=(0:100)/fs;s=sin(2*pi*t*5)+sin(2*pi*t*15)+sin(2*pi*t*30);%混和正弦波信号sf=filter(b,1,s);%对信号s进行滤波%凯泽窗设计FIR滤波器高通，带通，低通滤波器设计示例2：clearall;wb=0.317*pi;ws=0.272*pi;wa=0.136*pi;As=65;tr_width=wb-ws;%计算过渡带宽N=(ceil((As-7.95)*2*pi/(14.36*tr_width)+1))/2;%计算凯泽窗的长度disp(['
3、窗函数的长度为',num2str(N)]);%在Matlab命令窗口显示N的数值beta=0.1102*(As-8.7);%计算凯泽窗的线性相位斜率值，即凯泽窗系数n=[0:1:N-1];disp(['线性相位斜率为',num2str(beta)]);%在Matlab命令窗口显示线性相位斜率值w_kai=(kaiser(N,beta))';%凯泽窗函数hd_L=ideal_lp(wa,N);%调用理想脉冲函数程序,Lowpass，这里调用函数ideal_lph_L=hd_L.*w_kai;%实际脉冲响应=理想脉冲响应*凯泽
4、窗函数[db_L,mag,pha,grd,w]=freqz_m(h_L,[1]);%求幅频响应的各个参数，db为幅度响应，w为对应的频率值（弧度表示）%绘图程序figure(4)plot(w/pi,db_L);holdonhd_B=ideal_lp(wb,N)-ideal_lp(wa,N);%调用理想脉冲函数程序,Bandpassh_B=hd_B.*w_kai;[db_B,mag,pha,grd,w]=freqz_m(h_B,[1]);plot(w/pi,db_B,'m');holdonwc_H=(wb+ws)/2;hd_
5、H=ideal_lp(pi,N)-ideal_lp(wc_H,N);%调用理想脉冲函数程序,Highpassh_H=hd_H.*w_kai;[db_H,mag,pha,grd,w]=freqz_m(h_H,[1]);plot(w/pi,db_H,':r');legend('lowpass','bandpass','highpass',1);title('幅度响应/dB');axis([01-12010]);xlabel('以π为单位的频率');ylabel('分贝数/dB');%相关函数1function[db,mag,p
6、ha,grd,w]=freqz_m(b,a);%db是0到π内的相对振幅，mag是0到π内的绝对振幅，pha是0到π内的相位响应%grd是0到π内的群延迟，w是0到π内的501个频率样本向量%b是直接型滤波器分子多项式系数，a是直接型滤波器分母多项式系数[H,w]=freqz(b,a,1000,'whole');H=(H(1:1:501))';w=(w(1:1:501))';mag=abs(H);db=20*log10((mag+eps)/max(mag));%num2str(eps)的数值是2.22e-16pha=ang
7、le(H);grd=grpdelay(b,a,w);%相关函数2functionhd=ideal_lp(wc,N);alpha=(N-1)/2;%计算冲击脉冲响应的相位n=[0:1:(N-1)];m=n-alpha+eps;%加入eps，确保分母是非零数hd=sin(wc*m)./(pi*m);%直接传递函数设计低通滤波器示例3：Fs=50;N=25;tao=(N-1)/2;wc=2*15/50;%n是奇数H(N)为偶数对称，采用N为奇数的I型，取样化为归一化频率m=fix(wc/(2/N)+1);%fix向下取整omeg
8、a=[0:N-1]*2*pi/N;ABS_H=[ones(1,m),0.4,zeros(1,N-2*m-1),0.4,ones(1,m-1)];H=ABS_H.*exp(-j*tao.*omega);h=ifft(H);figure(1);subplot(2,1,1);stem(omega,abs(H)
