用凯泽窗设计线性相位带通FIR滤波器_第1页
用凯泽窗设计线性相位带通FIR滤波器_第2页
用凯泽窗设计线性相位带通FIR滤波器_第3页
用凯泽窗设计线性相位带通FIR滤波器_第4页
用凯泽窗设计线性相位带通FIR滤波器_第5页
全文预览已结束

下载本文档

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

文档简介

1、 用凯泽窗设计线性相位带通FIR滤波器 一 源程序： 1.计算理想带通滤波器单位抽样响应ha(n)的M函数文件dbpfr.mfunction hd=dbpfr(omega,N)dw=omega(2)-omega(1);alfas=40;M=ceil(alfas-7.95)/(2.286*dw);if mod(M,2)=0; N=M-1;else N=M;endn=0:N-1;wc1=1/2*(omega(1)+omega(2);wc2=1/2*(omega(3)+omega(4);alpha=(N-1)/2;hd=sin(wc2*(n-alpha)./(pi.*(n-alpha)-sin(wc

2、1*(n-alpha)./(pi.*(n-alpha);jj=find(n=alpha);hd(jj)=(wc2-wc1)/pi;2.计算凯泽窗函数wk(n)的M函数文件kwnd.mfunctionwk,N=kwnd(omega,alfas)dw=omega(2)-omega(1);alfas=40;M=ceil(alfas-7.95)./(2.286*dw);if mod(M,2)=0 N=M-1;else N=M;endif alfas=50 beta=0.1102*(alfas-8.7)else if alfas(eps*s) T=x/(2*n); n=n+1; D=D.*T.2; s=

3、s+D;end通带允许起伏:Rp=-1dB wp1=0.3*pi wp2=0.5*pi阻带衰减: As=40dB ws1=0.15*pi ws2=0.65*piclc;clear;Rp =-1;alfas=40;omega=0.15*pi 0.3*pi 0.5*pi 0.65*pi;wk,N=kwnd(omega,alfas)hd=dbpfr(omega,N)figure(1),clf;stem(hd);grid;title( hd(n);figure(2),clf;stem(wk);grid;title(w(n);hn=hd.*wkfigure(3),clfstem(hn);grid;tit

4、le(hd(n);Hn=fft(hn,128);figure(4),clfsubplot(1,2,1);plot(abs(Hn(1:64);grid;title(幅频响应曲线);subplot(1,2,2);plot(20.*log(abs(Hn(1:32);grid;title(对数幅频响应曲线);Hd=fft(hd,128);figure(5),clfplot(abs(Hd(1:65);grid;title(理想幅频响应曲线);二 截图三 思考题：1 窗函数的指标包括：旁瓣峰值幅度、主瓣宽度、过渡带宽度、阻带最小衰减2 根据阻带衰减可以选择窗函数的类型，根据过渡带宽计算窗函数的长度。3 N取奇数和偶数时，有

人人文库> 全部分类> 教育资料 > 辅导培训

温馨提示

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

最新文档

评论

0/150

提交评论