数字信号处理滤波器的设计.doc

 《数字信号处理》课程设计报告 设计课题 滤波器设计与实现 专业班级 电信1101班 姓 名 学 号 201105 报告日期 2013年12月 《数字信号处理》 课程设计任务书 题 目 滤波器设计与实现 学生姓名 甘源滢 学号 201105020103 专业班级 电信1101班 设 计 内 容 与 要 求 一、设计内容： 设计一个模拟低通巴特沃斯滤波器，通带最大衰减；阻带起始频率，阻带最小衰减，画出其幅度谱和相位谱。 二、设计要求 1 设计报告一律按照规定的格式，使用A4纸，格式、封面统一给出模版。报告内容 设计题目（）设计原理 （）设计（）收获和体会（）参考文献 √2， Ωc是3dB截止频率。在Ω=Ωc附近，随Ω加大，复制迅速下降。复制特性与Ω和N的关，N越大，通带越平坦，过渡带越窄，过渡带与阻带幅值下降的速度越快，总是频响特性与理想低通滤波器的误差越小。 以S代替jΩ,将幅度平方函数|Ha(jΩ)|^2写成s的函数： Ha(S)Ha(-S)=1/(1+(S/jΩ)^2N) 复变量S=δ+jΩ,此公式表明幅度平方函数由2N各极点，极点Sk用下公式表示： Sk=(-1)^1/2N(jΩc)=Ωce^jπ(1/2+(2k+1)/2N) 公式中K=0,1,2……,2N-1。2N个极点等间隔分布在半径为Ωc的圆上，间隔是 Π/Nrad。 为形成因果稳定的滤波器，2N个极点中只取S平面左半平面的N个极点构成Ha(S),而右半平面的N个极点构成Ha(-S)。Ha(S)的表达式为 Ha(S)=Ωc^N/Π(S-Sk) 设N=3，极点由六个，它们分别是 S0=Ωce^j2π/3 S1=-Ωc S2=Ωce^-j2π/3 S3=Ωce^jπ/3 S4=Ωc S5=Ωce^-jπ/3 取S平面的左半平面的极点S0,S1,S2组成系统函数Ha(S),即 Ha（S)=Ωc^3/(S+Ωc)(S-Ωc^j2π/3)(S-Ωc^-j2π/3) 由于不同的技术指标对应的边界频率和滤波器幅频特性不同，为使设计公式和图表统一将频率归一化。巴特沃斯滤波器采用对3dB截止频率Ωc归一化，归一化后的系统函数为 Ga(S/Ωc)=1/Σ(S/Ωc-Sk/Ωc) 令ρ=η+jλ=s/Ωc,λ=Ω/Ωc,λ称为归一化频率，ρ称为归一化复变量，这样的巴特沃斯滤波器的归一化低通原型系统函数为 Ga(ρ)=1/Σ(ρ-ρk) 公式中，ρk=sk/Ωc称为归一化极点，用下公式表示， ρk=e^jπ（1/2+(2k+1)/2N) k=0,1,…,N-1 显然，sk=Ωcρk 这样，只要根据技术指标求出阶数N,按照ρk=e^jπ（1/2+(2k+1)/2N)公式求出N个极点，再按Ga(ρ)=1/Σ(ρ-ρk)得到归一化低通原型系统函数Ga(ρ)，如果给定Ωc，再去归一化，即将ρ=S/Ωc，代入Ga(ρ)中，便可得到期望设计的系统函数Ha(S)。 2.2低通巴特沃斯滤波器的设计步骤如下 根据技术指标Ωρ，αρ，Ωs,αs,用λsp=Ωs/Ωρ ksp=√(10^αs/10-1)/(10^α/10-1) N=lgksp/lgλsp求出滤波器的阶数N。 按照ρk=e^jπ（1/2+(2k+1)/2N) k=0,1,…,N-1公式求出归一化极点ρk，将ρk代入Ga(ρ)=1/Σ(ρ-ρk)，得到归一化低通原型系统函数GA(ρ)。也可以根据阶数N直接查表得到ρk，Ga(ρ)。 将Ga(ρ)去归一化。将ρ=S/Ωc代入Ga(ρ)，得到实际的滤波器系统函数Ha（S)=G(ρ)|ρ=S/Ωc代入Ga(ρ)，得到实际的滤波器系统函数Ha（S)=G(ρ)|ρ=S/Ωc 这里Ωc为3dB截止频率，如果技术指标没有给出Ωc，可以由Ωc=Ωρ（10^αs/10-1）^（-1/2N) Ωc=Ωs（10^αs/10-1）^（-1/2N) 求出。 2.3函数说明buttord函数 [N,wc]=buttord(wp,ws,ap,as) 用于计算巴特沃斯数字滤波器的阶数N和3dB截止频率wc。其中，调用参数wp,ws分别为数字滤波器的通带、阻带截止频率的归一化值，要求：0≤wp≤1,0≤ws≤1（1表示数字频率pi）。当ws≤wp时，为高通滤波器；当wp和ws为二元矢量时，为带通或