巴特沃斯滤波器课程设计.docx

1、巴特沃斯滤波器课程设计摘 要本篇论文叙述了数字滤波器的基本原理，IIR数字滤波器的设计方法喝IIR数字高通滤波器设计在MATLABE上的实现与IIR数字滤波器在世纪中应用。无限脉冲响应（IIR）数字滤波器是冲击函数包含无限个抽样值的滤波器，一般是按照预定的模拟滤波器的逼近函数来转换成相应的数字滤波器，现有的逼近函数如巴特沃斯，切比雪夫。其设计过程都是由模拟滤波器的系统函数去变换出相应的数字滤波器的系统函数。关键字：数字滤波器，MATLAB，巴特沃斯，切比雪夫，双线性变换法ABSTRACTThe queue phenomenon in the telecom offices is a norma

2、l issue. To improve the customers satisfaction and to support the company changing, we should solve this problem properly. The basic goal to resolve queue problem is the appropriate tradeoff between the customers wavy demand and the telecom offices limited service capability.This paper is based on t

3、he queuing theory and demand management theory. And based on the data collection and customer survey and interview, the author uses some statistical methods to reflect the actuality. Then the author finds the reason of queuing in telecom office from customers view. Furthermore, the author analyses t

4、he real demand of the customers by sorting them into types of paying and time and price sensitivity. To follow up, three solutions had been brought forward: firstly, distributing the customers; secondly, stopping phone by different number; lastly, promoting the demand during the non- fastidious.Key

5、Words: queuing theory, demand management, telecom offices第1章引言1.1 课程设计目的了解巴特沃斯高通数字滤波器的概念及原理，掌握巴特沃斯高通数字滤波器的设计方法，综合运用专业知识及基础知识解决实际工程技术问题。 1.2 研究思路通过对目前主流各型滤波器的学习与研究，掌握滤波器设计的基本知识，再通过具体学习巴特沃斯型滤波器，正确地设计出一款巴特沃斯型高通滤波器。再运用模拟软件验证设计思路与实验数据。最终完全掌握巴特沃斯型滤波器的设计方法。第2章研究的理论基础2.1 数字滤波器数字滤波器是指输出输入均为数字信号，通过数字运算处理改变输入信

6、号所含频率成分的相对比例，或者滤除某些频率成分的数字期间或程序。因为数字滤波器通过数值运算实现滤波，所以狮子滤波器处理精度高，稳定，体积小，不存在阻抗匹配问题。可以实现模拟滤波器无法实现的特殊滤波功能。按照不同的分类方式，数字滤波器可以有很多种类型，但总的来说可以分为两大类：经典滤波器和现代滤波器。经典滤波器的特点是其输入信号中有用的成分和希望滤除的成分分别占有不同的频带，通过一个合适的选频滤波器滤除干扰，得到纯净信号，打到滤波目的。但是，如果信号和干扰信号的频谱相互重叠，经典滤波器则无法有效滤除干扰。现代滤波器是根据随机信号的一些统一特性，在某种最佳准则下，最大限度的抑制干扰，同时最大限度地

7、恢复信号，达到最佳滤波目的。2.2 模拟滤波器模拟滤波器分为无源滤波器和有源滤波器，其中无源滤波器由R、L、C组成，有源滤波器由集成运放和R、C组成，不需要使用电感。集成运放的开环电压增益和输入阻抗均很高，输出电阻小，构成有源滤波电路后还具有一定的电压放大和缓冲作用。但集成运放带宽有限，所以目前的有源滤波电路的工作频率难以做得很高。典型的模拟滤波器有巴特沃斯滤波器，切比雪夫滤波器，椭圆滤波器，贝塞尔滤波器等。这些典型的滤波器各有各的特点：巴特沃斯滤波器具有单调下降的幅频特性。切比雪夫滤波器的幅频特性在通带或阻带有等波纹特性，可以提高选择性。贝塞尔滤波器通带内有较好的线性特性。椭圆滤波器的特性相

8、对前三者是最好的，但通带和阻带内均呈现等波纹性，相位特性的非线性也稍严重。设计师需要根据具体要求选择滤波器的类型。第3章巴特沃斯滤波器的设计3.1 确定滤波器阶数假设通带频率在0.3内，通带幅度波动小于1Db，阻带频率范围在0.5rad内，阻带衰减大于12Db。数字指数：P=0.3，s=0.5求p，s利用频率畸变公式得：， 确定滤波器阶数： 确定系统函数： 3.2 巴特沃斯带阻滤波器系数计算带阻滤波器(又称陷波器)的作用是抑制某一频率信号而通过其他频率的信号。在Matlab中用b，a=butterord(n，Wn，stop，s)语句来实现带阻滤波器，其中b，a是滤波器传递函数多项式的系数;n为

9、滤波器的阶数;wn=fl*2*pi fu*2*pi分别为阻带低端一3 dB和高端一3 dB处的角频率;s表示模拟滤波器的。巴特沃特带阻滤波器返回分子分母多项式系数，还可以作出增益曲线，以及求出相位响应以及群延时响应。带阻滤波器遵循规律： Wp和Ws为二元矢量且Wp>Ws，如Wp=0.1,0.8,Ws=0.2,0.7。 b,a= buttord(ceil（n/2）,W1,W2,stop) % n为buttord（）设计出的带阻滤波器阶数 % 根据MATLAB帮组文件，butter(n,W1,W2，stop)将返回2*n阶滤波器系数 % Wn为带通滤波器截止频率，Wn=W1,W2,是2元向量

10、 % 注意，带通滤波器阶数是2倍关系 % b为H（z）的分子多项式系数 % a为H（z）的分母多项式系数第4章巴特沃斯滤波器设计仿真4.1 巴特沃斯低通滤波器实例仿真采样率为8000Hz，设计一个巴特沃斯低通滤波器，通带截止频率为2KHz，阻带截止频率为3KHz，通带波纹系数为1，阻带衰减系数为20，采样频率为10KHz，则只需在MATLAB的命令窗口下键入程序如下： f_N=8000; f_p=2100;f_s=2500;R_p=3;R_s=25;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2); n,Wn=buttord(Wp,Ws,R_p,R_s); b,a=butter(n,

11、Wn); freqz(b,a,1000,8000) subplot(2,1,1);axis(0 4000 -30 3) 程序运行后所设计出的巴特沃斯低通滤波器的频率响应如图4-1所示。横坐标为：angular frequency(red/s) 纵坐标为：phase(degrees)和gain(dB)。从图中可以得出，仿真后的频率响应满足的设计要求。 图4-1 巴特沃斯低通滤波器的频率响应图4.2 巴特沃斯高通滤波器实例仿真采样率为8000Hz，设计一个巴特沃斯高通滤波器，通带截止频率为2KHz，阻带截止频率为3KHz，通带波纹系数为1，阻带衰减系数为20，采样频率为 10KHz，f_p=100

12、0 Hz;f_s=700 Hz,Rp=3dB,Rs=20dB。程序如下： f_N=8000; f_p=1000;f_s=700;R_p=3;R_s=20;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2); n,Wn=buttord(Wp,Ws,R_p,R_s); b,a=butter(n,Wn，high); freqz(b,a,1000,8000) subplot(2,1,1);axis(0 4000 -30 3) 程序运行后所设计出的巴特沃斯高通滤波器的频率响应如图4-2所示。横坐标为：frequency(Hz) 纵坐标为：phase(degrees)和Magnitude(dB)。

13、从仿真图中得出，仿真后的频率响应满足设计的要求。图4-2 巴特沃斯高通滤波器的频率响应图4.3 巴特沃斯带通滤波器实例仿真采样率为10000Hz，设计一个巴特沃斯带通滤波器，通带截止频率为2KHz，阻带截止频率为3KHz，通带波纹系数为1，阻带衰减系数为20，采样频率为10KHz，f_p=1000 ,1500Hz;f_s=600,1900 Hz,Rp=3dB,Rs=20dB。程序如下： f_N=10000; f_p=1000,1500;f_s=600,1900;R_p=3;R_s=20; Ws=f_s/(f_N/2);Wp=f_p/(f_N/2); n,Wn=buttord(Wp,Ws,R_p

14、,R_s); b,a=butter(n,Wn); freqz(b,a,1000,10000) subplot(2,1,1);axis(0 5000 -30 3) 程序运行后所设计出的巴特沃斯带通滤波器的频率响应如图4-3所示。横坐标为：frequency(Hz) 纵坐标为：phase(degrees)和Magnitude(dB)。从仿真图中得出，仿真后频率响应图满足的设计要求。图4-3 巴特沃斯带通滤波器的频率响应图4.4 巴特沃斯带阻滤波器实例仿真与低通滤波器的幅频、相频特性类似，带阻滤波器的幅频、相频特性如图四所示。在该设计中，电阻R1，R2，R3，电容C1的取值需要确定，它们的计算公式可

15、以通过Matlab编程来得到，最终根据计算结果选择相应的标称值。采样率为10000Hz，要求设计一个带阻滤波器，f_p=1000 ,1500Hz;f_s=1200,1300Hz,Rp=3dB,Rs=30dB。 程序如下： f_N=10000; f_p=1000,1500;f_s=1200,1300;R_p=3;R_s=30;Ws=f_s/(f_N/2);Wp=f_p/(f_N/2); n,Wn=buttord(Wp,Ws,R_p,R_s); b,a=butter(n,Wn，stop); freqz(b,a,1000,10000) subplot(2,1,1);axis(0 5000 -35 3

16、) 程序运行后所设计出的巴特沃斯带阻滤波器的频率响应如图4-4所示。从图中得出，频率响应满足的设计要求。图4-4 巴特沃斯带阻滤波器的频率响应图第5章基于需求管理的CD电信营业厅排队优化措施5.1 顾客的需求分析需求管理的首要任务便是对顾客的总需求进行分析，研究总需求是否可以被分为不同的子需求。（略）5.1.1 两类缴费顾客的对比本文主要研究的是去营业厅缴纳话费的那部分顾客。提取2006年1月至10月办理费用缴纳业务的用户数据平均值，得到各种缴费方式的用户比例情况，如图5-1所示（红旗连锁代收方式数据不全，不计入其中）。数据来源：CD电信财务数据库图5-1各种缴费方式用户比例图中，（略）5.1

17、.2 顾客对时间和价格的敏感度分析我们用一个时间和价格敏感度的二维坐标矩阵对缴费顾客进行需求划分。图5-2TS矩阵通过图5-2所示，该二维坐标矩阵，我们将顾客分为四类16：（略）5.2 基于需求管理的优化措施5.2.1 分流顾客5.2.1.1分流顾客的基本思路 分流顾客是（略）5.2.1.2分流顾客的理论依据分流顾客的措施的理论依据在于（略）5.2.1.3分流顾客的具体措施根据（略）5.2.2 分号停机5.2.2.1分号停机的基本思路 分号停机策略是（略）5.2.2.2分号停机的理论依据通常情况下，（略） 5.2.2.3分号停机的模型建立（略）由此，依据分号停机的理论基础，我们在建立分号停机模

18、型的时候，可以借用信号与系统中的关于线性时不变系统平稳冲激响应的方法。时不变性的含义是，如果系统的输入在时间上有一个平移，则由其产生的响应也产生一个同样的平移20，即： （5-1）处理该问题时，我们把停机缴费顾客具有的三个特性，抽象为一个系统，该系统具有以下三个特性：（略）设表示具有这三个特性的系统，并且在模型中始终保持不变，表示没有采用策略前第N天的停机顾客数，表示第N天的缴费顾客数 21，则： （表示卷积） （5-2）（略）据此，令表示分号停机时第N天的停机顾客数，表示分号停机时第N天的缴费顾客数，表示被停机的尾号为单号的数量，表示被停机的尾号为双号的数量，那么，则： （5-3）由于系统具

19、有时不变性，采用分号停机策略前的一个月内，准确的说是6号，尾号为双号的停机用户的数量应该和策略后，即当月21号时，尾号为双号的停机用户数相等，也就是说延迟15天停机不改变尾号为双号的停机顾客总数，则： （5-4）另外，由于系统具有线性的特点，我们假设每月的停机单双号数量相等，有： （5-5） （5-6）将公式5-4、5-5、5-6，代入公式5-3，可获得以下近似恒等变形： （5-7）上式表明，（略）通过MATLAB建模，利用06年CD电信10、11月新华营业厅日缴费顾客数的数据，仿真分号停机策略实施后的日缴费户数情况如图5-3、图5-4所示（MATLAB程序命令见附录三）。图5-3分号停机策略

20、实施后停机缴费顾客数图5-4分号停机策略实施后日缴费顾客数通过以上两图（略）5.2.3 促进非高峰期需求5.2.3.1促进非高峰期需求的基本思路促进非高峰期需求策略是（略）5.2.3.2促进非高峰期需求的理论依据促进非高峰期需求策略是（略）5.2.3.3促进非高峰期需求的具体措施促进非高峰期需求的基本措施有两种（略）5.3 优化实施后的排队模型为了（略）参考文献1克里斯托弗.H.洛夫洛克. 服务营销(第三版).北京:中国人民大学出版社,2001.78802胡运权等.运筹学教程(第2版).北京:清华大学出版社,2003.1281323 Richard B.Chase, Nicholas J. A

21、quilano, F. Robert Jacobs. 运营管理(原书第9版).北京:机械工业出版社,2003.1451634Madeleine E Pullman, Gary Thompson. Strategies for integrating capacity with demand in service networks, Journal of Service Research, 2003, Vol.5,No.6:16205Mark M. Davis等. 服务管理(第1版).北京:人民邮电出版社, 2006.1321486Madeleine E Pullman, Gary M Thom

22、pson, Evaluating capacity- and demand-management decisions at a ski resort, Cornell Hotel and Restaurant Administration Quarterly, 2002, Vol.43,No.6:25377陈敏. 基于需求管理的成都地区供电企业服务营销策略研究:工商管理硕士（MBA）学位论文. 成都:四川大学,20048Kenneth J Klassen, Thomas R Rohleder, Combining operations and marketing to manage capac

23、ity and demand in services, The Service Industries Journal, 2001, Vol.21,No.2:25309Rhyne, David M, the Impact of Demand Management on Service System Performa, the Service Industries Journal, 1988, Vol.8,No.4:263110Douglas R Moodie，Demand management: The evaluation of price and due date negotiation s

24、trategies using simulation，Production and Operations Management, 1999, Vol.8, No.2:152311Eitan Altman, Uri Yechiali, Analysis of customers impatience in queues with server vacations，Queuing Systems, 2006, Vol.52,No.4:616312盛骤, 谢式千等. 概率论与数理统计(第3版). 浙江大学出版社, 2001.636513James A Fitzimmons等. 服务管理(第5版).

25、北京:机械工业出版社, 2007.18520714Ward Whitt, Sensitivity of Performance in the Erlang-A Queuing Model to Changes in the Model Parameters, Operations Research, 2006, Vol.54,No.2:141815James A Fitzimmons等. 服务管理(第5版). 北京:机械工业出版社, 2007.29830416蔡滟. 浅谈消费者的需求管理. 商业时代学术评论, 2006, Vol.7,No.5:444817罗鹏. 服务营销中的排队管理. 营销策

26、略, 2006,Vol.6,No.7:505118李平英. 排队现象及其管理研究. 山东农业大学学报, 2003,Vol.1,No.9:687219姜兵. 银行营业网点客户排队问题研究:高级工商管理硕士（EMBA）学位论文. 南京: 南京大学, 200520乐友正. 信号与系统(第1版). 北京:清华大学出版社, 2004.374221Charles L. Phillips等.Signals, Systems, and Transfor(第3版).北京:机械工业出版社, 2006.586322James A Fitzimmons等. 服务管理(第5版). 北京:机械工业出版社, 2007.12

27、913723Kenneth J Klassen, Thomas R Rohleder, Combining operations and marketing to manage capacity and demand in services, The Service Industries Journal, 2001, Vol.21,No.2:272924张准. 论服务业中的顾客需求管理. 内蒙古财经学院学报, 2004, Vol.4,No.6:444625刘洋. 中国移动自有营业厅服务管理研究:硕士学位论文.北京:北京邮电大学, 2006致谢我要首先感谢我的课题设计老师，谢谢她一直帮我完善论文

28、，其次是我的同学我们一起查找资料，一起用心完成论文，我最开始不知道怎么做，就找老师，上网搜资料，去图书馆查找等等，终于皇天不负有心人，在不断地资料搜寻当中我们渐渐的明白了我们将要做的设计原理和相关知识。 在我们一起查找资料的过程中渐渐学会了很多东西，让我知道了滤波器这个东西，还有让我不再懒惰，这篇论文完成的真的不容易。 通过整个设计过程，让我感到最深刻的就是MATLAB功能的强大性掌握变成各种函数和语句的重要性。当然首先要了解所要编程运行的对象的原理。在课程设计的过程中，我深深的感受到我们所学的东西太少了，需要学习的东西太多了，真的是学无止境。但同时也是快乐的，让我们大家朝着自己的目标努力奋斗

29、吧，快乐生活，快乐学习，终将有一天我们会成功的 。附录附录一：06年11月CD电信新华营业厅日缴费顾客半小时到达数据（略）附录二：顾客调查问卷（略）附录三：仿真分号停机策略的MATLAB程序命令（略）附录四：分号停机策略实施前后日缴费顾客数量对比表（略）外文资料原文Combining operations and marketing to manage capacity and demand in servicesKenneth J Klassen, Thomas R RohlederINTRODUCTION: DEMAND IN SERVICESServices have long struggled with the challenge of managing capacity and demand.（略）.译文在服务中结合运作与营销的需求与能力管理Kenneth J Klassen, Thomas R Rohleder介绍：服务中的需求服务长久以来与能力和需求管理的挑战作斗争。（略） （注：可编辑下载，若有不当之处，请指正，谢谢!）