自-dsp应用系统设计综合设计报告 .docVIP

《DSP应用系统设计》综合设计报告 题目:切比雪夫FIR滤波器滤波 专 　 业 　_ 通信工程_＿　 班级 _ 　　 _ 　 开课学期　 　 提交时间 　 摘要 数字信号处理（Dｉｇital Sigｎaｌ Pｒｏcesｓｉｎg,简称DSP)是一门涉及许多学科而又广泛应用于许多领域的新兴学科。 DSP（Ｄigｉｔａl　Sｉgnal　Processor）是一种独特的微处理器，是以数字信号来处理大量信息的器件。其工作原理是接收模拟信号,转换为0或1的数字信号。再对数字信号进行修改、删除、强化，并在其他系统芯片中把数字数据解译回模拟数据或实际环境格式。它不仅具有可编程性，而且其实时运行速度可达每秒数以千万条复杂指令程序,远远超过通用微处理器,是数字化电子世界中日益重要的电脑芯片。它的强大数据处理能力和高运行速度,是最值得称道的两大特色。 一、设计要求 设计一个切比雪夫Ⅱ型带阻数字滤波器。参数：Fs＝16KＨz,阻带边界频率、,阻带最小衰减；通带边界频率、,通带最大衰减。 二、系统方案 １、在CＣS4集成开发环境下编写ＦIR滤波器程序，并保存。 2、根据设计要求的参数利用MATLＡB计算出切比雪夫I型低通数字滤波器的系数,放大取整后代入上一步写好的程序中，进行仿真调试。 三、算法设计 １、滤波器设计 ?数字滤波器是将输入信号序列，按规定的算法进行处理，从而得到所期望的输出序列。一个线性位移不变系统的输出序列y（n）和输入序列x（n)之间的关系，应满足常系数线性差分方程: y(ｎ)=i=0 x(n)为输入序列，y（ｎ)为输出序列,ai，b y（ｎ）＝i=0 对（2)式进行z变换，可得ＦIＲ滤波器的传递函数： 　　 　　　 H（z）＝ Y（z ２、FIR滤波器函数： 四、软件实现 (一)、滤波器系数设计 1、在ＭATLＡB中调用契比雪夫滤波器的函数得到符合设计要求的数字滤波器系数,并将系数放大取整,备后面的程序使用。 其中用到两个函数:cｈｅb1oｒd(wｐ,ws,Rp,Ａs)、cheby2(N,Rs,wｃ,’ftype’) ｃｈeb１ord(wｐ,ｗs，Ｒｐ,Ａｓ)函数用来得到切比雪夫Ｉ型数字滤波器的阶数Ｎ和阻带截止频率wc。 Chebｙ１(N,Ｒs,wｃ,’ftｙpe’）函数用来得到线性相位的切比雪夫Ｉ型数字滤波器的滤波器系数。 (二）、FIR滤波 ?根据FIR滤波原理,用ｃ语言编写程序。程序中调用了TI提供的库函数firｓ（);故需要在程序开始添加头文件Ｄｓplib.h。 五、实验结果与分析 ?1、matｌab中滤波器的频率响应: ２、ccｓ中dｅbug模式下观察输入输出波形： ?（１）输入波形 （２)输出波形 六、结论 根据输入输出的图形可以看出,输入信号是５00Ｈｚ、1800Hz，3300Hｚ三种频率信号的叠加，经过滤波处理后只有500Hz的信号输出,说明滤波成功。输出波形有轻微的失真,可见滤波对原始信号会产生一定的损伤。 七、附件:源程序清单 ?Ｍaｔｌab中计算滤波器系数： ?Ｆs=1600０; ａp＝1; as=３0； fp=８00; fs=１600； wp＝２*fp／Fs; wｓ=２*ｆs/Fs; [N,wso］=cheb1ｏrd(ｗp,ｗs，ap，aｓ)； ｈ=cheby１(N,as,wso）; ccs环境下FIR滤波程序： #incｌｕdｅ ｍａth.h＂ #inｃｌｕde Ｄspｌib.h＂ ＃dｅｆｉｎe　uinｔ unｓiｇned inｔ ＃definｅ PI 3.1415926 ＃ｄefiｎｅ L 16 #dｅfine　L2 L/2 ＃defiｎe Ｆs 8０00 #define T　1/Fs #defｉnｅ f1 300 #deｆiｎｅ　f2　1８00 ＃dｅｆiｎe　f3 330０ #ｄefiｎe ｗ1　(2*PＩ＊ｆ1＊T) ＃ｄeｆine　w２ （２*PＩ*f2*T） #ｄeｆinｅ　ｗ3 (2*PI*f3*T) #ｄefinｅ a1 0.３3 ＃define ａ2 0．30 #ｄｅfine ａ3 ０.3０ DATA coefｆ[L2]={０,0,2,9,27，6１，１0１,13０}; //DATA coeff[L2]=｛15,0,-9２,0,230,0,-307}; DATＡ　iｎ[L]； DATA out[Ｌ]; ＤATＡ iｎｄｅｘ[L+2］； main() 　{ ?uint i,n=０； float sｉｇnal; ?ｆor（i=０;iL；ｉ+＋) ? ｛ ?in［i]　＝ 0; ?oｕt［i] = 0； ｝ ?for（i=0；i＜Ｌ+2;i＋+） ? 　{ ??indｅx[i］ = 0； ?　 ｝ while(1） 　 { ? signa