7.3 实验二——用SIMULATOR模拟实现数字信号处理.docVIP

7.3 实验二——用SIMULATOR模拟实现数字信号处理 本实验主要通过四个练习来了解如何利用SIMULATOR实现基本的信号处理方法。实验包括以下几个练习： 练习一：时域卷积运算 练习二：DFT运算 练习三：时域相关运算 练习四：利用相关函数计算信号的功率谱 7.3.1 练习一——时域卷积运算 练习一的目的是使学生掌握VisualDSP++的使用，同时加深对数字信号处理理论的理解。这是一个较简单的计算两个输入序列的卷积程序。两个输入序列的长度相同，分别写入到几个文件中。程序编译运行时将数据读入存储器中，并对其进行卷积运算，计算结果存放于另一个缓冲区中。 1． 程序说明 设输入序列X的长度为M，Y的长度为N，则程序按照因果系统计算卷积： , 程序中一个输入序列X指定到数组inputx中，另一个输入序列Y指定到数组inputy中，卷积的输出指定到数组output中。根据卷积原理，两个输入序列长度为M=N=64，输出序列长度为M+N-1=127。 程序主要由主程序和卷积子程序（位于pm_code段中）二个模块组成。在主程序中，先将输入数组X复制到输出数组的后半部分，然后调用卷积子程序conv计算卷积。 卷积计算子程序中，将计算按照输出长度分成两段完成，在程序中利用left和right两个二重循环实现。Left循环计算m=0~N-1，right循环计算m=N~N+M-1。在计算卷积时，序列X采用正常序寻址，利用指令F5=DM(I9,M9)读取X(k)，其中地址增量寄存器M9＝1，开始计算时寄存器I9指向序列X的第一个元素。而序列Y是逆序寻址，利用指令F0=DM(I1,M2)读取Y(m-k)，其中地址增量寄存器M2＝－1，开始计算时寄存器I1指向序列Y的最后一个元素。 实验测试程序主要模块的源程序如下： /******************************************************************** conv.asm ADSP-2106x 程序功能：计算X[]与Y[]的卷积 ******************************************************************/ #include def21060.h /* 存储器映射的IOP 寄存器定义 */ #define LENGTH_X 64 /* 符号常数 */ #define LENGTH_Y 64 .SECTION/DM dm_data; /* 定义数据存储区中的变量 */ .VAR inputx[LENGTH_X]= pulse1.dat; .VAR inputy[LENGTH_Y]= sin64.dat; .VAR output[LENGTH_X+LENGTH_Y-1]= zeros.dat; .SECTION/PM pm_rsti; /* 复位向量段 */ NOP; USTAT2= 0x108421; /* 复位后程序从这里开始执行 */ DM(WAIT)=USTAT2; /* 设置外部存储器访问为0等待 */ JUMP start; .SECTION/PM pm_code; start: M1=1; I0=inputx; L0=0; I3=output+LENGTH_Y-1; L3=0; /* 把数组inputx复制到output的后半段（位于LENGTH_Y－1开始） */ LCNTR=LENGTH_X, do copydata until lce; F0=DM(I0,M1); copydata:DM(I3,M1)=F0; CALL conv; /* 调用子程序,进行卷积计算 */ end: IDLE; /* 等待 */ /********* 计算时域相关性的子程序　****************************/ conv: I2=output; M1=1; M2=-1; M9=1; F15=0; /* 用于给累加器清零 */ /* 计算前 LENGTH_X 个数据 */ LCNTR=LENGTH_X, DO left UNTIL LCE; I9=I2; I1=inputy+LENGTH_Y-1; /* 指向inputy的结尾 */ F9=PASS F15, F0=DM(I1,M2), F5=PM(I9,M9); LCNTR=LENGTH_Y, DO in