DSP原理与应用实验报告

1、DSP原理与应用实验报告于矿班级：学院：指导教师:实验一代数汇编指令基础实验验目的：1. 通过调试目标代码，掌握指令的功能，熟悉指令2. 通过指令的熟悉，能够指令应用丁实际项目中。二、实验原理：Ti公司的代数汇编指令。三、实验程序：.title"算术指令综合实验".mmregs.sect".vect".copy"vectors.asm".textStart:；AR7=#767;A=#38CAH;DP=#08AH;RSA=#0123H;DP=#188H;ASM=#0AHAR7=AMMR(=HI(A)<<ASM|B=*AR2-

2、<<16AR2+A=DBL(*AR5-)-AAR7+)=#1234HDP=#04HA=#9876HAR6=#230H*AR6+=#9ACDHARP=#6J38H=ANOPNOPAR0=#003AHA=*AR6+0T=#08HSXM=1B=*AR6-<<TSAR7=#0220H*AR7=#0E5D8HA=*AR7+0B<<-123AH=#665AHAR2=#023AHBK=#0A=*AR2-0%<<12BK=#002AHAR3=#0256H-*AR3-=HI(B)<<8ASM=#-8*AR3+=A<<ASMAR4=#0216

3、H*AR4=#8775HB=#025AHA=B<<ASMMMR(*AR4+)=BA=MMR(13H)AR2=#0245H*AR2=A<<73AH=A<<ASM*AR2+0%=B<<4;*AR2=#1CHASM=*AR2T=*AR4+LTD(*AR2+)B=RND(*AR2+)*AR3=#0F57AHA=UNS(*AR3)B=*AR3+*AR3+0B=TTRN=#12ACHA=DBL(*AR4+)B=DUAL(*AR2-)DBL(*AR3-)=ADUAL(*AR4+)=B*AR3+=HI(B)<<ASM|T=*AR3+A=#3456HI

4、F(AGT)*AR4+=HI(A)<<ASMB=#0F679HIF(BLEQ)*AR3-=HI(B)<<ASMA=#0F98DHIF(AGT)*AR2+=BRCIF(ALT)*AR3-=TB=#0125CHIF(BGEQ)*AR2+=BRCCMPS(A,*AR4-)CMPS(B,*AR2+);B=20DP=#40A=#1234H22=AA=A+9AH25=AAR3=#0236H*AR3=#0F775HAR5=#024AH*AR5=#09ACDHNEXT:NOPA=#9ABCH*AR5+=AA=A+#1000HA=A+#08ADEHSXM=0A=#07AB8H<&l

5、t;16A=A+#04ADEHA=A-#08ADEH<<16SXM=1B=#0FF7CHA=#0889AHA=A-#09ACDH<<16B=A-*AR5-C=1B=B-A<<ASM*AR5=AC16=0B=B-*AR3+A=B+*AR5+<<16A=A-B<<ASMB=*AR3+<<16-*AR4-<<16A=A+*AR4<<12B=A+*AR5<<-12A=B-#06789H<<16B=B+*AR7+0B+CARRYA=A-*AR2-BORROWSUBC(25,A)A=A-

6、UNS(*AR7+)T=#9ACDHB=DADST(*AR4,T)A=DADST(*AR7,T)C16=1A=A+DBL(*AR5+)A=A-DBL(*AR5-)T=#7654HC16=0A=DADST(*AR5+,T)C16=1A=DSADT(*AR5-,T)A=DBL(*AR5+)-A*AR3+=HI(B)|B=A+*AR5+0%<<16*AR4-=HI(A)|A=*AR3-<<16-BGOTONEXT;SXM=0A=#89ABHA=A+#4567H<<16*AR3=#9999HA=A&*AR3-B=#8897HB=B+#079ADH<&l

7、t;16A=A|B<<-12A=#8897HA=A+#079ADH<<16A=BA#0567DH<<12DP=#047AH=7AH&#0ACD6H.end四、实验步骤：1、输入以上程序，并进行编译；2、打开codeExplore,并把编译好的程序装载，并进行调试3、逐步调试，并观察各种特殊寄存器的值和预期值是否对应，着重了解各种寻址的特点，及相对应的指令。五、实验心得第一次实验并没有掌握程序的操作方法以及对丁程序的不理解，导致操作不熟练，而且程序错误较多，花费了大量的时间在程序的纠错和调试上。实验二DSPFIR低通滤波器一、实验目的1、通过实验程序理

8、解DSPS序的编写。2、掌握CIC-500系统的A/D转化，申口通信的编程。3、掌握FIR低通滤波器在DSP中的具体实现，进一步了解各种寻址方式的实际应用。二、实验原理FIR滤波器的主要操作控制为设计一个低通滤波器，条件如下：1、采样频率为9.26KHZ,系统工作频率为10MHZ2、截至频率为Flpc为1KHZ3、输入信号为14位带符号数值。4、共取80阶点作运算。滤波器设计好后的系数放在"lp_coeff.asm”文件内其他文件应拷贝到同一文件夹下。三、FIR设计方法1、数字滤波器设计步骤数字滤波器的实质是一个离散的多项式表达式，通过这个表达式来满足一定的需求。因此数字滤波器的设计

9、步骤如下：A、了解滤波器的技术要求B、实现系数计算C、实现，即把传函转换为合适的滤波器网络或结构»有限字长效应分析，主要考虑滤波器系数、输入数据量化和固定字长对滤波器性能的影响E、工程实现:2、滤波器的技术指标p:峰值通带偏差;:阻带偏差p:通带边缘频率S:阻带边缘频率3、用窗口法设计滤波器的步骤第一步：指定理想的毗望的滤波器频率响应HD;第一步：通过傅立叶反变换求期望的滤波器的冲击响应，第三步：选择一个满足通带或衰减指标的窗函数，然后利用滤波器长度与过渡带宽之间的关系确定滤波器的系数数目。第四步：对于选取的窗函数求系数h(n)。w(n)的值，并且将hDn的值与w(n)相乘求得实际的

10、FIR4、标准频率选择性滤波器的理想冲击响应总结滤波器类型hDnhD0低通2fc罕室2fc高通sinnc-2fc-1-2fc带通sinn2wsinn2fc-2f1n2n12f2-2f1带阻sinn1sinn22f-2f2n1n21-2f2+2f15、常用窗函数及重要特征总结窗函数名归一化过渡带宽HZ通带波纹dB相对于旁瓣的主瓣dB阻带衰减dB窗函数表达式矩形窗0.9/N0.761413211Hanning窗3.1/N0.05463144一一2n0.50.5cosNHamming®3.3/N0.01944153一一一2n0.540.46cosN布莱克曼5.5/N2.93/N0.0017

11、0.0274577550一_一2n4n0.420.5cos0.08cosN1N1凯塞4.32/N5.71/N0.002750.000275709020.5I。12n/N1四、随机数产生方法实验程序范例.title"lowpassfilter".mmregs.width80.length55;.setsect”.text”,0x1800,0;theseassemblerdirectivesspecify;.setsect”.data”,0x0200,1;theabsoluteaddressesofdifferent;.setsect"vectors”,0x0180,

12、0;sectionsofcode.sect"vectors"interruptvectortableresidesat.copy"lp_vecs.asm"location0x0180.dataseed.word07e6dh;seedforrandomvariabletemp.word0XN.word0,0,0,0,0,0,0,0,0,0;80datalocationsfor80XN1.word0,0,0,0,0,0,0,0,0,0;stagedelayline.XN2.word0,0,0,0,0,0,0,0,0,0;XN3.word0,0,0,0,0,0

13、,0,0,0,0;XN4.word0,0,0,0,0,0,0,0,0,0;XN5.word0,0,0,0,0,0,0,0,0,0;XN6.word0,0,0,0,0,0,0,0,0,0;XN7.word0,0,0,0,0,0,0,0,0;XNLAST.word0;OUTPUT.word0;extrawordforthebitbucket.text.copy"lp_coeff.asm.copy"lp_ac01.asmdcallAC01INIT;initializeAnaloginterface.start:intm=1;disableallinterruptsDP=#0nop

14、pmst=#01a0h;Interruptpointermapsvectorstopage#3(3*128=180h)sp=#0ffah;stackpointerlocatedinCommunicationsKernalimr=#240h;unmaskTDMRINTandHPIINT(hostportinterface)intm=0;enableallinterruptsWAIT:gotoWAITreceive:DP=#seed;waitforreceiveinterrupt.;ThissetsDataMemoryPagePointertopageXN,whichisdefinedearlie

15、rintheprogram.-randomnoiseGenerator(P-5csModulator)a=seed<<1a=seedAatemp=a<<2a=tempaaa=#8000h&aa=a+seed<<16seed=hi(a)<<1a=seed<<11a=a&#0fffch<<15repeat(#12)a=a<<C-1getsampleandrunthroughlowpassfilter-b=DRR1;LOADACCUMULATORWITHWORD;RECEIVEDFROMAIC!XN=

16、A<<0;STORETHEVALUEOFRECEIVED;WORDTOVARIABLEXN!AR0=#XNLAST;LOADAR0WITHADDRESSOFLAST;DELAYELEMENT!A=#0;ZEROACCUMULATORA!repeat(#79);Repeatnextinstructions80times.macd(*AR0-,h0,A);ComputeFIRoutput.OUTPUT=hi(A)<<0;Storethefilteredinputinto;variableOUTPUT.A=OUTPUT<<0;OUTPUT=>Accumula

17、torAA=#0FFFCh&A;TWOLSB'sMUSTBEZEROFORAIC!DXR1=A;SENDTOTRANSMITREGISTER!return_enable;Enableinterruptsandreturn;frominterrupt.transmit:return_enable;Enableinterruptsandreturn;frominterrupt.end四、实验步骤1、连接CIC-500DSP发展实验系统"2、载入程序，从J2端输入音乐，运行一段时间后暂停，观察实验结果。3、观察不同频宽的音乐，把结果进行比较。_五、实验心得通过实验了解到了D

18、SP程序的编写和A/D转化，以及串口通信的编程，进一步了解了各种寻址方式的实际应用，依据上次实验的经验，这次的操作过程更为熟练。实验三FFT算法的实现一、实验目的1、掌握FFT算法。-2、进一步掌握FFT的DSP编程。3、指出FFT算法中的各种寻址方式，结合算法掌握寻址方式的使用技巧。4、认真分析各功能模块的具体功能，和课堂所讲的算法步骤关联起来。二、实验原理实验原理依据课堂所授的蝶形变换设计1024个点的FFT变换，输入数据放在in.dat文件里面。在输入实验的时候把bitrev.asm、fft.asmunpack.asm、power.asm拷贝至同一个目录下-实验主程序：.width80.

19、length55.title"RealFastFourierTransfrom".mmregs.copy"vectors.asm”|.copy"initrfft.asm".textstart:SP=#BOSDP=#0OVM=#0FRCT=#1ASM=#0CALLbit_revCALLfftCALLunpackCALLpowerhangloose:NOPgotohangloose.copy"bitrev.asm”.copy"fft.asm".copy"unpack.asm".copy"power.asm".end三、实验步骤：1、认真阅读四个子程序，理解四个子程序的作用2、认真进行调试，通过时域频域图比照算法的精度3、认真思考，改进算法。四、实验心得掌握FFT算法，以及相关的FFT的DSP编程。认真分析各功能模块的具体功能，和课堂所讲的算法步骤关联起来，通过实验能够了解到FFT算法的一些简单的应用以及算法的运算准则。实验七设计PID控制器或设计矩阵乘法程序一实验目的1、根据以前做实验的所培养的经验，写出PID控制器或矩阵乘法程序二、实验原理1、PID控制原理2、矩阵乘法原理3、矩阵乘法程序main.as