dsp高通滤波器课程设计报告

《dsp高通滤波器课程设计报告》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。

1、DSP课程设计报告题目FIR高通滤波器设计姓名学号教学院系专业年级指导教师目录一、设计题目1二、设计目标1三、算法研究与参数计算11、FIR的原理和参数生成公式12、利用MATLAB计算滤波系数13、输入信号参数计算2U!五、调试过程41、调试前准备52、MATLAB的使用53、编写及编译程序54、设置断点和探针65、打开观察窗U6六、实验结果及分析61、输入信号的时域波形和频域波形62、输出信号的时域波形和频域波形7七、设计心得8一、设计题目FIR高通滤波器设计二、设计目标设计一个FIR高通滤波器，通带边界频率为6000Hz，采样频率为20000HzoFIR滤波器的设计

2、用MATLAB窗函数法进行。三、算法研究与参数计算1、FIR的原理和参数生成公式(1)N阶有限冲激响应滤波器(FIR)公式：NZ2-1v(w)=A(fc)[x(n-fr)+x(”一(?/一1+介))]公式1-1*-0N=0X2…"FIR设计原理：咕•门.'根据系数h是偶对称为T简化运算产生如卜计算方法如果一个FIR滤波有一个冲激响应.h(0).h(l)，….h(Nd).和x(n)描绘输入的时常滤波n*输出滤波y(n)的n给山以下方积式：y(n);h(0)x(n)+h(1)x(n-l)+h(2)x(n-2)+".+h(N-1〉x[n，(n-1)]公式1—2图3-12、利用M

3、ATLAB计算滤波系数在MATLAB界面输入图3-2所示程序，可得到滤波系数并生成INC文件y=round(firl(64,6000/20000*2/high')*32768);fid=fopen(’0126.inc’w’);fprintf(fid/’.word%d?,y);fclose(fid)图3-2输入freqz（y，1，512），MATLAB中显示高通滤波器的滤波特性曲线。如图3-3所示。图3-33、输入信号参数计算MATLAB中输入图3-4中所示程序，包含两种频率成分的正弦信号，一种信号频率1000Hz，一种信号6000Hz。CommandWindow»i=

4、0:1:255;y=round((sin(2*pi*[i]*1000/20000)+sin(2*pi*[i]*6000/20000))*32768/2);fid=fopen(,1207030126.dat’，’w’)；fprintf(fid/16512010);fprintf(fid,%d’，y);fclose(fid)图3-4其频谱特性曲线如0.1020.3040.50.607NormalizedFrequency(xnrod/somplc)080.9■WW200080S•4—4—(gprpfuceesx1O•2NormaluedFrequuncy(xnrad/

5、sarnple)图3-5(SSJ6-S)S2H编写源程序参考资料，编写汇编语言源程序：HIGHPASS.set1;ifyouwanttouse,pleasesetthevalueto1.globalstart,fir.mmregsCOFF_FIR_START:.sect”coff_firn.copy’’IO126.incnK_FIR_BFFR.set64d_data_buffer.usect’’fir_bfrn，64FIR_DP.usect”fir_vars’’，0d_filin.usect’’fir_vars’’,loutput.usect”fir_vars’’，lin

6、put.usect’’fir_vars’'，ld_filout.usectnfir_vars’’，100hstacksize.set256stack.usectnfir_vars”，stacksize•asgAR4，FIR_DATA_P•asgAR6，INBUF_P•asgAR7,OUTBUF_P•asgAR3,OUTBUF•asgAR2，INBUF.sect’’fir—prog”nopstart:stm#stack+stacksize，SPLD#FIR_DP，DPSTM#ddatabuffer,FIRDATAPRPTZA，#K一FIR一BFFR-1STLA，*FIR_D

7、ATA_P+STM#dfilinJNBUFPSTM#d_filout,OUTBUF_PSTM#output,OUTBUFSTM#input，INBUFSTM#100h，BKfir一loop:NOP;AddBreakpoint&porbepointLD*INBUF一P,ASTLA,*INBUFCALLfirSTHA，*OUTBUF一P+%STHA，*OUTBUFmainend:bfir一loopfir:;SSBXSXMSSBXFRCTSTMSTLSTM#d_data_buffer，FIR_DATA_PA，*FIRDATAP#(d