DSP技术及应用课程设计报告.doc

dsp技术及应用课程设计报告学院：通信与信息工程学院专业： 电子信息工程专业 班级： 电信*班 姓名： 艾依河里的鱼 时间： 2012年1月1日 一、设计任务及要求设计任务：本课程设计要求应用des5402pp-u实验板实现fir滤波器，能够在ccs集成开发环境中使用图形显示工具显示输入和输出波形。设计要求：利用des5402pp评估板的模拟信号输出通道产生一个1khz的方波，然后将产生的方波信号通过模拟信号输入通道输入des5402pp-u实验板，对方波进行低通滤波，得到一个1khz的正弦信号。二、设计原理本设计中我们使用一个38阶的对称结构的fir低通滤波器，其采样频率fs为25khz，通带截止频率1.2khz，阻带截止频率为2.8khz，阻带衰减为40db。我们可以应用matlab中的设计fir滤波器的指令fir()，也可以应用matlab的工具fdatool来进行设计，然后将所设计的滤波器的系数导入dsp程序中。下面对fdatool进行简单的介绍：fdatool(the filter design and analysis tool)是一个设计分析滤波器的工具，它具有很好的用户界面。借助它可以快速的设计数字iir或fir滤波器。你可以设置滤波器的性能指标，也可以直接配置滤波器的系数，还可以从matlab工作空间导入滤波器。在fdatool界面下，你可以方便的看到设计的滤波器的幅频特性、相频特性、零极点分布、单位脉冲响应、阶跃响应、滤波器系数、群时延等。要使用fdatool，只要在matlab命令窗口中键入，按回车即可。fdatool有两种模式：设计模式和输入模式。在设计模式中，你可以从零开始设计，或者对已经设计的滤波器进行改进，或者分析滤波器。在输入模式中，你可以导入先前已经设计好的滤波器或者已经存在于matlab工作空间的滤波器系数，你还可以分析导入的滤波器。本设计中我们主要使用fdatool的设计模式进行fir滤波器的设计和对所设计的滤波器进行分析，然后可以用file菜单下的export将滤波器的系数导入matlab的workspace和将系数写入dsp汇编程序中。在写系数到汇编程序中，要注意小数的表示方法。三、程序代码.title fir5402.asm .ref _c_int00 .mmregs；寄存器地址设置gpiocr .set 3ch ；5402通用i/o控制寄存器gpiosr .set 3dh ；5402通用i/o状态寄存器dxr10 .set 23h ；串口0数据发送寄存器drr10 .set 21h ；串口0数据接收寄存器spsa0 .set 38h ；串口0子地址寄存器spcr10 .set 39h ；多通道串口0寄存器1spcr20 .set 39h ；多通道串口0寄存器2xcr10 .set 39h ；串口0发送控制寄存器pcr0 .set 39h ；串口0引脚控制寄存器dxr11 .set 43h ；串口1数据发送寄存器spsa1 .set 48h ；串口1子地址寄存器spcr11 .set 49h ；多通道串口1寄存器1spcr21 .set 49h ；多通道串口1寄存器2xcr11 .set 49h ；串口1发送控制寄存器pcr1 .set 49h ；串口1引脚控制寄存器srgr1 .set 49h ；采样率发生寄存器1srgr2 .set 49h ；采样率发生寄存器2swcr .set 2bh ；软件等待状态寄存器 ；定义变量timer_count .set 060h ；定时器计数mcbsp_init .set 061h ；mcbsp初始化wave_ptr .set 062his_new_data .set 063h ；新数据到达标志new_ad .set 064hout_wave_ptr .set 065ht_ar2 .set 066ht_ar3 .set 067htemp .set 068hwave_buf .set 0900h out_wave_buf .set 01000hn .set 38 ; fir tapsfir_coef_buf .set 100h ; fir coef bufferfir_data .set 200h ; fir windows data buffer ! .text;-; interrupt vector table !;-rs b _c_int00 nop nopnmi b _ret nop nopsint17 b _ret nop nopsint18 b _ret nop nopsint19 b _ret nop nopsint20 b _ret .word 0,0sint21 b _ret .word 0,0sint22 .word 01000h .word 0,0,0sint23 .word 0ff80h .word 0,0,0sint24 .word 01000h .word 0,0,0sint25 .word 0ff80h .word 0,0,0sint26 .word 01000h .word 0,0,0sint27 .word 0ff80h .word 0,0,0sint28 .word 01000h .word 0,0,0sint29 .word 0ff80h .word 0,0,0sint30 .word 01000h .word 0,0,0int0 b _ret nop nopint1 b _ret nop nopint2 b _ret nop noptint b timer ；内部定时中断 nop nopbrint0 b receive ；缓冲串口接收中断 nop nopbxint0 b transmit ；缓冲串口发送中断 nop noptrint b _ret nop noptxint b _ret nop nopint3 b _ret nop nophpint b _ret nop nopq26 .word 0ff80h .word 0,0,0q27 .word 01000h .word 0,0,0q28 .word 0ff80h .word 0,0,0q29 .word 01000h .word 0,0,0q30 .word 0ff80h .word 0,0,0q31 .word 01000h .word 0,0,0;-; 复位中断服务程序 ;-_c_int00: stm #2020h,pmst ; vector table start: 0x2000，看iptr中的值 ssbx intm ; close all int ! (ssbx intm) ssbx sxm ; extend sign ! ssbx frct ; faction multiply stm #10h,26h ; stop timer0 ! stm #10h,36h ; stop timer1 ! stm #0ffh,sp ; sp= 0x0ff ld #0,dp ; dp=0 stm #0ffffh,ifr ; clear all int ! 中断标志寄存器 stm #20h,imr ; bit5-bxint0, bit4-brint0, enable bxint0 ；中断屏蔽寄存器 stm #02492h,swwsr ; all 2 waits ! stm #0,swcr;/*-; 配置dsp的时钟des5402pp-u板的上电默认配置没有将dsp的时钟设置到100m，而只有20m。下面的程序将pll改为10，以便dsp以100m全速运行。pll的时钟工作方式寄存器clkmd，地址为0058h，见dsp技术及应用教材248页。;-*/ stm #0,58h wait_pll: ldm 58h,a and #1,a bc wait_pll,aneq stm #97ffh,58h ; switch pll*10 - 100m clk rpt #100 nop;*;初始化ram, variable.;* stm #wave_buf,ar6 stm #1,ar0 rpt #24 mvpd wave_data,*ar6+0 ; move wave data to wave_buf - 0x0900h stm #wave_buf,wave_ptr ; save wave bufferptr stm #out_wave_buf,out_wave_ptr ; save out wave bufferptr, because use ; circular address, so when save out_wave ; date,must + ! stm #fir_coef_buf,ar6 rpt #n-1 mvpd fir_coef,*ar6+ ; move fir coef to fir_coef_buf(in data mem) stm #fir_coef_buf,t_ar2 ; stm #fir_data,t_ar3 stm #0aa55h,mcbsp_init ; when mcbsp_init=0xaa55, transmit for init ac01 stm #5,timer_count ; init variable timer_count ! stm #0,is_new_data ; if is_new_data =1 , has new ad data by slave ac01 ! stm #1000h,ar0 stm #1800h,ar5 stm #400h,bk ld *ar0,t ld #0,b;*; the following codes are used to initalize mcbsp0 !;*;- stm #0,spsa0 ; spsa0为串口0的子地址寄存器，choose spcr10 stm #2000h,spcr10 ; receive sign_extend in drr stm #1,spsa0 ; choose spcr20 stm #0c0h,spcr20 ; fs - int ! stm #2,spsa0 ; choose rcr10 stm #40h,39h stm #3,spsa0 ; choose rcr20 stm #0,39h stm #4,spsa0 ; choose xcr10 stm #40h,xcr10 stm #5,spsa0 ; choose xcr20 stm #0,39h stm #0eh,spsa0 ; choose pcr0 stm #0dh,pcr0 ; fs is low active, rise edge of clkx, falling edge ; of clkr stm #7h,spsa0 stm #8000h,39h rpt #0ffffh nop stm #00h,dxr10 ldm 22h,a stm #1,38h stm #0c1h,39h ; start mcbsp0 send ! ld #603h,a rsbx intm ; enable all int !;-; the following codes initalizeac01, all ac01 setup same !;- idle 1 idle 1 idle 1 idle 1 ; 软件复位两个ac01 ld #600h,a idle 1 ; wait for int . ld #3,a idle 1 ; send #3 ld #00105h,a ; fs=10m/2/a/b= 25k ; flp=10m/2/40/a= 25k (all pass) ; fhp=fs/200= 125hz (not use !) ; a=05,b=40(0x28) idle 1 ; 送05h到ac01的数据寄存器1 ld #3,a idle 1 ; send #3 ld #00228h,a idle 1 ; 送28h到ac01的数据寄存器2 ld #3,a idle 1 ; send #3 ld #00300h,a idle 1 ; 送00h到ac01的数据寄存器3 ld #3,a idle 1 ; send #3 ld #00405h,a idle 1 ; 送05h 到ac01的数据寄存器4(ad & da 0 db) ld #3,a idle 1 ; send #3 ld #00505h,a idle 1 ;送05h到ac01的数据寄存器5(highpass filter ;disable) ld #3,a idle 1 ; send #3 ld #00600h,a idle 1 ; 送00h 到ac01的数据寄存器6 ld #3,a idle 1 ; send #3 ld #00700h,a idle 1 ; 送0到ac01的数据寄存器7 ld #3,a idle 1 ; send #3 ld #800h,a idle 1 ; 送0到ac01的数据寄存器8 ld #0h,a stm #0,38h stm #2001h,39h ;开始mcbsp0 接收! stm #30h,imr ; bit5-bxint0, bit4-brint0, 使能 bxint0;-;*; the following codes are used to initalize timer !;*;- stm #10h,tcr ; stop timer ! stm #0ffffh,prd ;装载prd, 周期寄存器 stm #0fh,tcr ; timer start, each about 100hz . orm #8,imr ; enable timer int !;-;* *; the following codes are main function;* *main: bitf is_new_data,#1 ; test is_new_data 1 ? cc fir,tc ;1, then call fir ! new ad data in new_ad(0x64) b main ；检测到新数据到达标志（变量is_new_data）为1时，调；用fir滤波程序。;* *; the following codes are served for mcbsp0 transmit !;* *;-transmit: pshm st1 pshm st0 pshm bk pshm ar7 pshm ah pshm al cmpm mcbsp_init,#0aa55h bc init,tc stm #25,bk ; each ac01 has 25 dots to d/a ! mvdm #wave_ptr,ar7 ; restore wave buffer ptr ! nop ld *ar7+%,a nop stlm a,dxr10 ; send wave data ! master ac01 make source wave ; slave ac01 make filter result wave ! bd wave_make_end mvmd ar7,#wave_ptr ; save wave buffer ptr !init: stlm a,dxr10 bitf dxr10,#800h ; if setup last register of ac01, set mcbsp_init=0 ! bc wave_make_end,ntc ; not equal ,jump ! stm #0,mcbsp_init ; when mcbsp=0 or 1, transmit for make wave !wave_make_end: popm al popm ah popm ar7 popm bk popm st0 popm st1_ret: rete;- ;*; the following codes are serived for mcbsp0 receive ! use ar5,bk!;*;-receive: pshm st1 pshm st0 ld #0,dp mvkd drr10,ar5 ; save ad data pshd *ar5+% ; save ad data in ad_buffer(0x1800-0x1c00) popd new_ad stm #1,is_new_data ; has new ad data popm st0 popm st1 rete-;*; the following codes are serived for vc5402 timer !;*;-timer: pshm ah pshm al ldm timer_count,a sub #1,a bc not_show_led,aneq ; if a != 0 ,jump ! bitf st1,#2000h ; test xf bc show_led,tc ssbx xf ; set xf=1 b show_conshow_led rsbx xf ; set xf=0show_con ld #5,a ; not_show_led stlm a,timer_count popm al popm ah rete;- ;*; the following codes is used to fir( new ad data in new_ad(0x64); fir buffer coff- 0x100(ar2), fir window buffer - 0x200(ar3);*fir: pshm st1 pshm st0 pshm bk pshm ar7 pshm ah pshm al pshm ar0 pshm ar2 pshm ar3mvdm #t_ar2,ar2 ; restore ar2 mvdm #t_ar3,ar3 ; restore ar3 ld new_ad,-2,a ; new ad data - a, remove two lsbs ;- stm#1,ar0stm#n,bk ; set circular addressing sizestla,*ar3+%rptza,#(n-1) ; 0 - a, then repeat 34 timesmac*ar2+0%,*ar3+0%,a ; done fir filter, result in ahmvmdar3,#t_ar3 ; save ar3mvmdar2,#t_ar2 ; save ar2 stm #1,ar0 stm #100h,bk ; each ac01 has 25 dots to da ! mvdm #out_wave_ptr,ar7 ; restore wave buffer ptr ! sth a,-2,temp ; save fir result in temp ld temp,2,a ; shift 2 bit for ac01 da ! stl a,*ar7+0% ; save fir result to wave buffer ! stl a,temp mvmd ar7,#out_wave_ptr ; save wave buffer ptr ! stm #0,is_new_data ; new ad data has been processed ! ld #0,dp cmpm out_wave_ptr, 1000h bc con_fir,ntc ;- ccs_show: ；此处加断点，查看结果 nop ; if out_wave_ptr = 0xc00, nop ; may show out_wave in ccs ! show 128 dots ! ; input buffer - 0x1800 ; output buffer - 0x1000 (show from 0x1020) con_fir: popm ar3 popm ar2 popm ar0 popm al popm ah popm ar7 popm bk popm st0 popm st1 ret;*; the following data is used by make wave ! when using, must copy to; wave buffer - 0x6000;*wave_data: ; freq about 1k hz .word 03ffch ; +2 volt .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 03ffch .word 0c000h ; -2 volt .word 0c000h .word 0c000h .word 0c000h .word 0c000h .word 0c000h .word 0c000h .word 0c000h .