安康学院学年论文(设计)
开始 读入窗口长度N 选择窗函数形式 由相应的频率响应计算hd?n? 调用窗函数子程序求h?n? 加窗计算h?n??hd?n?w?n? 计算频率响应 绘制仿真结果曲线
图2.1 窗函数仿真设计流程图
3.应用凯瑟窗设计高通滤波器
3.1凯瑟窗需计算的参数
元素dev必须包含特定化顶峰波动的值或通带、止带允许的偏移量。
然后dev?10??rp/20??1/10??rp/20??1,10???rs/20?,(10??rp/20??1)/10??rp/20??1。应用?n,Wn,beta,ftype??kaiserord?fc,mags,dev,fs?得到kaiser窗的尺寸并进行尺寸修正。
[5]
????????3.2凯瑟窗设计及仿真
将模拟频率转化为数字频率,设取样时间为T,由已知的指标,可得T?程序如下ωs=0.5π,ωp=0.6π,dp=0.5 dB,ds=50 dB, Wdel=wc-wp; Beta=1; %ds>21
N=ceil((As-8)/2.285/wdel) %N=5 Wn=Kaiser (N,bata); ws=(wp+wc)/2/pi; b=fir1(N-1,ws,’high’,wn); figure(1)
第5页 共7页
1?0.0001s, f安康学院学年论文(设计)
freqz(b,1);
所得仿真结果如下图3.1所示:
图3.1 Kaiser窗滤波器幅频与相频响应
4.结语
总之,数字滤波器的应用广泛,利用MATLAB语言,可以较容易地设汁出滤波器,且滤波器的调整也较方便。利用MATLAB及计数字滤波器有多种方法.我们要灵活运用这蝗方法,从而在工程上有提高设计效率。
滤波器从功能上分类, 可分为低通、 高通、 带通&、带阻滤波器, 这种理想滤波器是不可能实现的, 只能作为逼近标准用。 从单位脉冲响应分类\可分为无限脉冲响应IIR滤波器和有限脉冲响应FIR滤波器。FIR数字滤波器的单位脉冲响应是有限长序列,IIR滤波器由于不具备FIR滤波器所固有的稳定性和设计方法的简单性。所以一般不使用,Matlab所有滤波器没计函数都将返回一个FIR滤波器,FIR数字滤波器具有描述方便、系统总是稳定的,易于实现的。
[参考文献]
[1] Joint Video Team, Draft ITU-T Recommendation and Final Draft International Standard
第6页 共7页
安康学院学年论文(设计)
of JointVideo Specification, ITU-T Rec. H.264 and ISO/IEC 14496-10 AVC[S]. March 2005. [2]王宏.MATLAB6.5及其在信号处理中的应用[M].北京:清华大学出版社.2005:159.160. [3]毛青青,徐分亮.窗函数及其应用[J].中国水运(学术版).2007年第7卷第2期.
[4]张志田,何其文,基于窗函数法的FIR数字滤波器设计[J].湖南科技学院学报,2007年第28卷09
期.
[5]程佩青,数字信号处理教程[M].北京清华大学出版社2001.
第7页 共7页