频率采样法设计带通滤波器matlab.doc_第1页
频率采样法设计带通滤波器matlab.doc_第2页
频率采样法设计带通滤波器matlab.doc_第3页
频率采样法设计带通滤波器matlab.doc_第4页
全文预览已结束

下载本文档

版权说明：本文档由用户提供并上传，收益归属内容提供方，若内容存在侵权，请进行举报或认领

文档简介

以MATLAB 用频率采样法设计一带通数字滤波器。低通带边缘：w1p=0.2*pi,低阻带边缘：w1s=0.35*pi,高通带边缘w2p=0.65*pi,高阻带边缘w2s=0.8*pi;设计过渡带中的频率样本值为T1=0.109021，T2=0.59417456 % wp1=0.45*pi;wp2=0.65*pi;ws1=0.3*pi;ws2=0.75*pi;Rp=1;Rs=40% =双线型变换法=wp1=0.45*pi; wp2=0.65*pi; ws1=0.3*pi; ws2=0.75*pi;Rp=1; Rs=40;Wp1=tan(wp1/2); Wp2=tan(wp2/2);Ws1=tan(ws1/2); Ws2=tan(ws2/2);BW=Wp2-Wp1; W0=Wp1*Wp2; W00=sqrt(W0);WP=1; WS=WP*(W02-Ws12)/(Ws1*BW);N,Wn=cheb1ord(WP,WS,Rp,Rs,s);B,A=cheby1(N,Rp,Wn,s);BT,AT=lp2bp(B,A,W00,BW);num,den=bilinear(BT,AT,0.5);h,omega=freqz(num,den,64);subplot(2,2,1);stem(omega/pi,abs(h);xlabel(omega/pi);ylabel(|H(z)|);subplot(2,2,2);stem(omega/pi,20*log10(abs(h);xlabel(omega/pi);ylabel(增益.dB);% =直接法=wp1=0.45*pi; wp2=0.65*pi; ws1=0.3*pi; ws2=0.75*pi;Rp=1; Rs=40;Wp=wp1/pi,wp2/pi; Ws=ws1/pi,ws2/pi;N,Wn=cheb1ord(Wp,Ws,Rp,Rs);B,A=cheby1(N,Rp,Wn);h,omega=freqz(B,A,64);subplot(2,2,3);stem(omega/pi,abs(h);xlabel(omega/pi);ylabel(|H(z)|);subplot(2,2,4);stem(omega/pi,20*log10(abs(h);xlabel(omega/pi);ylabel(增益.dB);%cheby2% wp1=0.45*pi;wp2=0.65*pi;ws1=0.3*pi;ws2=0.75*pi;Rp=1;Rs=40% =双线型变换法=wp1=0.45*pi; wp2=0.65*pi; ws1=0.3*pi; ws2=0.75*pi;Rp=1; Rs=40;Wp1=tan(wp1/2); Wp2=tan(wp2/2);Ws1=tan(ws1/2); Ws2=tan(ws2/2);BW=Wp2-Wp1; W0=Wp1*Wp2; W00=sqrt(W0);WP=1; WS=WP*(W02-Ws12)/(Ws1*BW);N,Wn=cheb2ord(WP,WS,Rp,Rs,s);B,A=cheby2(N,Rs,Wn,s);BT,AT=lp2bp(B,A,W00,BW);num,den=bilinear(BT,AT,0.5);h,omega=freqz(num,den,64);subplot(2,2,1);stem(omega/pi,abs(h);xlabel(omega/pi);ylabel(|H(z)|);subplot(2,2,2);stem(omega/pi,20*log10(abs(h);axis(0 1 -100 0);xlabel(omega/pi);ylabel(增益.dB);% =直接法=wp1=0.45*pi; wp2=0.65*pi; ws1=0.3*pi; ws2=0.75*pi;Rp=1; Rs=40;Wp=wp1/pi,wp2/pi; Ws=ws1/pi,ws2/pi;N,Wn=cheb2ord(Wp,Ws,Rp,Rs);B,A=cheby2(N,Rs,Wn);h,omega=freqz(B,A,64);subplot(2,2,3);stem(omega/pi,abs(h);xlabel(omega/pi);ylabel(|H(z)|);subplot(2,2,4);stem(omega/pi,20*log10(abs(h);axis(0 1 -100 0);xlabel(omega/pi);ylabel(增益.dB);%椭圆% wp1=0.45*pi;wp2=0.65*pi;ws1=0.3*pi;ws2=0.75*pi;Rp=1;Rs=40% =双线型变换法=wp1=0.45*pi; wp2=0.65*pi; ws1=0.3*pi; ws2=0.75*pi;Rp=1; Rs=40;Wp1=tan(wp1/2); Wp2=tan(wp2/2);Ws1=tan(ws1/2); Ws2=tan(ws2/2);BW=Wp2-Wp1; W0=Wp1*Wp2; W00=sqrt(W0);WP=1; WS=WP*(W02-Ws12)/(Ws1*BW);N,Wn=ellipord(WP,WS,Rp,Rs,s);B,A=ellip(N,Rp,Rs,Wn,s);BT,AT=lp2bp(B,A,W00,BW);num,den=bilinear(BT,AT,0.5);h,omega=freqz(num,den,64);subplot(2,2,1);stem(omega/pi,abs(h);grid;xlabel(omega/pi);ylabel(|H(z)|);subplot(2,2,2);stem(omega/pi,20*log10(abs(h);grid;xlabel(omega/pi);ylabel(增益.dB);% =直接法=wp1=0.45*pi; wp2=0.65*pi; ws1=0.3*pi; ws2=0.75*pi;Rp=1; Rs=40;Wp=wp1/pi,wp2/pi; Ws=ws1/pi,ws2/pi;N,Wn=ellipord(Wp,Ws,Rp,Rs);B,A=ellip(N,Rp,Rs,Wn);h,omega=freqz(B,A,64);subplot(2,2,3);stem(omega/pi,abs(h);grid;xlabel(omega/pi);ylabel(|H(z)|);subplot(2,2,4);stem(omega/pi,20*log10(abs(h);grid;xlabel(omega/pi);ylabel(增益.dB);wls=0.2*pi;wlp=0.35*pi;whp=0.65*pi;wc=wlp/pi,whp/pi;B=wlp-wls;N=ceil(8/0.15);n=0:N-1;window=hanning(N);h1,w=freqz(window,1);figure(1);stem(window);axis(06001.2);grid;xlabel(n);title(Hanning窗函数);figure(2);plot(w/pi,20*log(abs(h1)/abs(h1(1);axis(01-3500);grid;xlabel(w/pi);ylabel(幅度(dB);title(Hanning窗函数的频谱);hn=fir1(N-1,wc,hanning(N);h2,w=freqz(hn,1,512);figure(3);st
人人文库> 全部分类> 生活休闲 > 科普知识

温馨提示

  • 1. 本站所有资源如无特殊说明，都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
  • 2. 本站的文档不包含任何第三方提供的附件图纸等，如果需要附件，请联系上传者。文件的所有权益归上传用户所有。
  • 3. 本站RAR压缩包中若带图纸，网页内容里面会有图纸预览，若没有图纸预览就没有图纸。
  • 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
  • 5. 人人文库网仅提供信息存储空间，仅对用户上传内容的表现方式做保护处理，对用户上传分享的文档内容本身不做任何修改或编辑，并不能对任何下载内容负责。
  • 6. 下载文件中如有侵权或不适当内容，请与我们联系，我们立即纠正。
  • 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。

最新文档

评论

0/150

提交评论