数值信号处理课程的设计的报告.docVIP

西安课程设计报告内容 实验一 实验题目： Using the bilinear transformation and a lowpass analog Butterworth prototype, design a highpass digital filter operating at a rate of 20kHz and having passband extending to 5kHz with maximum passband attenuation of 0.5dB,and stopband ending at 4kHz with a minimum stopband attenuation of 10dB. 实验目的： 进一步加深对用双线形变换法设计数字高通滤波器的理解。 学会利用低通滤波器设计高通滤波器。 熟练掌握低通滤波器的设计方法。 进一步熟悉MATLAB提供的各种设计滤波器函数。 实验原理： 双线形变换法是利用s=2*（1-z-1）/T*（1+z-1）将s域转换到z域，从而得到系统函数H（Z）。由于巴特沃斯低通滤波器和高通滤波器只有在归一化频率时和去归一化时关系式互为倒数，因此只要在这两处变换系数，就可以利用沃斯低通滤波器来设计高通滤波器。 实验步骤简述： （1）先将数字频率转换成角频率。 （2）再将数字角频率转换成模拟角频率。 （3）计算模拟低通滤波器传输函数。 （4）再将模拟低通滤波器传输函数转换成模拟高通滤波器传输函数。 （5）利用双线性变法将模拟高通滤波器转换成数字滤波器。 程序流程图： 六 源程序： 程序代码： clear all; fp=5000;fs=4000; Fs=20000; %采样频率 rp=0.5;rs=10; wp=2*pi*fp/Fs; ws=2*pi*fs/Fs; wap=2*Fs*tan(wp/2);was=2*Fs*tan(ws/2); [N,wn]=buttord(wap,was,rp,rs,s); %确定器阶数N和截止频率wn [z,p,k]=buttap(N); %确定归一化传输函数 [bp,ap]=zp2tf(z,p,k); %零极点增益滤波器参数转换为传输函数形式 [bs,as]=lp2hp(bp,ap,wn); %低通滤波器转换为某一截止频率的高通滤波器 [bz,az]=bilinear(bs,as,Fs); %双线性变换法 [h,f]=freqz(bz,az,256,Fs); %求频率响应 subplot(2,1,1); plot(f,abs(h));grid on; xlabel(f);ylabel(|H(w)|); title(|H(w)|的幅度谱图); db=20*log10((abs(h)+eps)/max(abs(h))); subplot(2,1,2); plot(f,db);grid on; xlabel(f);ylabel(|H(w)|); title(|H(w)|以dB为单位的幅度谱图); clear all; fp=5000; fs=4000; F=20000; rp=0.5; rs=10; T=1/F;wp=2*pi*fp/F; ws=2*pi*fs/F; Wp=2*tan(wp/2)/T; Ws=2*tan(ws/2)/T; W1=1/Ws;W2=1/Wp; [N,wn]=buttord(W1,W2,rp,rs,s); %巴特沃思滤波器的阶估计，计算阶数，3db截止频率 Wn=1/wn;fprintf(滤波器的阶数N=%.0f\n,N); [b,a]=butter(N,1,s); [bs,as]=lp2hp(b,a,Wn); [numd,dend]=bilinear(bs,as,F);% disp(分子系数b) fprintf(%.4e ,numd);fprintf(\n); disp(分母系数a); fprintf(%.4e ,dend);fprintf(\n); w=linspace(0,pi,1024); h=freqz(numd,dend,w);%计算频率响应H(z),在[0，2pi]上采样,w决定采样频率点 subplot(211),plot(w/pi,abs(h)); xlabel(归一化频率);ylabel(|H(w)|); subplot(212),plot(w/pi,20*log10(abs(h))); axis([0 1 -50 1]); xlabel(归一化频率);yl