fir滤波器线性相位 实验三用凯塞窗设计线性相位带通FIR滤波器-副本.docVIP

fir滤波器线性相位 实验三用凯塞窗设计线性相位带通FIR滤波器-副本 中国地质大学(北京) 实 验 报 告 课程名称：数字信号处理 实验名称：用凯塞窗设计线性相位带通FIR滤波器 姓 名： 张淑坤 学 号: 1004133105 班 级指导教师：陈玉东 评 分： 实验时间：2015-12-31 用… 行为金融学对投资模式的影响及对我国投资现象解释 作者：杨培垌 来源：《商业时代》2013年第14期 内容摘要：本文介绍了行为金融学，对理性投资者在有效市场追求预期效用最大化的挑战展开研究，并结合我国投资市场的有效性、资产定价、证券市场黑庄频频，崇拜股… 课程名称：战略管理 论文题目：丰田汽车公司的核心竞争力 姓名：崔毅 班级：工管121 学号：201210944109 目录 1 引言 .........................................................… 中国地质大学(北京) 实 验 报 告 课程名称：数字信号处理 实验名称：用凯塞窗设计线性相位带通FIR滤波器 姓 名： 张淑坤 学 号: 1004133105 班 级指导教师：陈玉东 评 分： 实验时间：2015-12-31 用凯塞窗设计线性相位带通FIR滤波器 实验目的： 基于MATLAB环境，熟悉利用窗函数法设计线性相位FIR滤波器的原理和方法。 实验环境： 硬件环境:计算机,软件环境：MATLAB平台。 实验环境原理： 凯塞窗是对于给定的阻带衰减，使主瓣具有最大能量意义下的最佳窗函数，因此具有最陡的过渡带。 窗函数的形式为 I0[ ??(1? wk(n)? 2n2 ) ]N?1， 0?n?N?1 I0(?) ? 其中?为形状参数，I0[x]是第一类修正零阶贝塞尔函数，其幂函数展开式为 1x I0[x]?1??[()k]2 k?1k!2 给定滤波器的过渡带宽度?? (rad)和阻带衰减As(dB)，则滤波器的长度和形状参数?可由下列经验公式给出： N? As?7.95 2.286Δ? As?50dB?0.1102(As?8.7)， ? β??0.5842(As?21)0.4?0.07886(As?21)，21?As?50dB ?0，As?21dB? 设计流程图 实验内容要求： 1. 编写计算理想带通滤波器单位抽样响应hd(n)的M函数文件dbpfilter_FIR.m，各变量定义 如下： dbpfilter_FIR(d_omega, N, beta) % d_omega ---输入数字频率数组（向量）； % d_omega(1)---阻带下边缘截止频率 % d_omega(2)---通带下边缘截止频率 % d_omega(3)---通带上边缘截止频率 % d_omega(4)---阻带上边缘截止频率 % N ---数字带通滤波器的长度 % hd ---理想数字带通滤波器单位冲激响应 % h ---实际数字带通滤波器单位冲激响应（所设计的） %wk ---凯泽窗序列 % beta ---凯泽窗参数 2. 编写计算凯泽窗函数wk(n)的M函数文件kaiser_WF.m，各变量定义如下： function wk= kaiser_WF(N, beta) %iser_WF.m文件，设计下列带通FIR滤波器： 通带允许起伏 ≤l dB，?p1=0.3?，?p2=0.5? 阻带衰减 ≤40 dB，?s1=0.15?，?s2=0.65 实验内容及结果 1、MATLAB程序 function [hd]=ideal_LP(wc,N) %Ideal Lowpass filter computation %------------------------------------------ %[hd]=ideal_LP(wc,M) %hd ---理想低通滤波器单位冲激响应（0alpha=(N-1)/2; n=[0:(N-1)]; m=n-alpha+eps; hd=sin(wc*m)./(pi*m); function wk= kaiser_WF(N, beta) for n=1:1:N wk(n)=bessel_IM (beta*sqrt(1-(1-2*(n-1)/(N-1))))/bessel_IM(beta); end function s=bessel_IM (x) N=,0.65]*pi delta_w=d_omega(2)-d_omega(1); N=(As-7.95)/(2.286*delta_w) N=ceil(