DSP系统设计与实现实验二报告.docxVIP

实验二应用程序开发实验一. 实验目的本实验通过五个练习来了解如何利用ADSP-21161实现基本的信号处理方法，这些常用算法包括：时域卷积运算、DFT运算、时域相关运算、利用相关函数计算信号的功率谱和利用SIMD模式实现DFT运算。二．实验内容实验二包括以下五个基本练习：练习一：时域卷积运算；练习二：DFT运算；练习三：时域相关运算；练习四：利用相关函数计算信号的功率谱；练习五：利用SIMD 模式实现DFT 运算。三．实验步骤练习一：时域卷积运算1．实验步骤1）调入程序启动运行VisualDSP++，运行菜单命令\Project\New 来新建一个工程，工程名称自定义，正确设置工程的各个选项，将DSP_exp\unit_2\Conv目录下的源文件（conv.ASM 和conv.ldf）加入到工程中，同时也将备用的数据文件pulse1.dat、pulse4.dat、square1.dat、square2.dat、sin64.dat、trig64.dat 以及初始化文件zeros.dat 加入工程中。2）选择输入数据文件编译程序前，在下述程序行中将输入数据文件改为pulse1.dat 和sin64.dat：.VAR inputx[LENGTH_X]= pulse1.dat;.VAR inputy[LENGTH_Y]= sin64.dat;其中pulse1.dat 为单个冲击脉冲信号，sin64.dat 为正弦信号。3）编译程序分析、理解源程序，在此基础上再编译运行程序。可通过点击按钮命令或选择菜单命令\Project\Build Project 来完成编译。本程序的sessions 为ADSP-21161 ADSP-2116x Simulator。4）观察结果按F5 运行程序，然后执行菜单命令\View\Debug Windows\Plot 命令，分别作出输入数据inputx[]和inputy[] 与输出数据output[]的图形。执行菜单命令\Debug\Step Into 单步运行程序，打开Register File 和对应的存储区，观察寄存器和存储器的变化。5）改变数据再运行?数组inputx保持为pulse1.dat，将数组inputy分别改变为数据文件square1.dat(方波)、trig64.dat（三角波）重复3),4)，来考察程序的执行结果。?数组inputx改变为pulse2.dat（4 个冲击脉冲），将数组inputy分别改变为数据文件sin64.dat，square1.dat(方波)、trig64.dat（三角波）重复3),4)，来考察程序的执行结果。?数组inputx改变为square1.dat，将数组inputy改变为数据文件square2.dat(方波)，重复3),4)，来考察程序的执行结果。2. 实验结果Inputx为64个点的脉冲序列，Inputy为64个点的正弦序列输入波形inputx为：图1 输入序列inputx(四个冲击脉冲的序列)输入波形Inputy为图2 输入序列inputy(正弦序列)脉冲序列与正弦序列的卷积结果为:图3 卷积结果3.结果分析：本次实验中分别进行了单个冲击脉冲与方波，三角波的卷积，四个冲击脉冲组成的序列与正弦序列，方波，三角波的卷积以及两个不同形状的方波的卷积。任何信号与冲激函数的卷积相当于以冲击脉冲所在的位置为中心对原信号进行移位，原信号形状保持不变，而四个冲击脉冲组成的序列与其他函数的卷积的效果就是将与其卷积的函数进行移位，然后再进行相加，最后形成输出叠加的波形。方波序列的卷积会形成梯形波，实验的结果与理论推出结果是一致的。下面讨论一下关于卷积的长度问题：设输入序列X的长度为M，Y的长度为N，则程序按照因果系统计算卷积：程序中一个输入序列X指定到数组inputx中，另一个输入序列Y指定到数组inputy中，卷积的输出指定到数组output中。根据卷积原理，两个输入序列长度为M=N=64，输出序列长度为M+N-1=127。练习二：DFT运算通过本练习加深对各种基本波形的DFT的理解，加深对数字信号处理理论的理解。本练习采用64点离散傅里叶变换。1. 实验步骤1) 调入程序启动运行VisualDSP++，建立新工程，将位于将DSP_exp\unit_2\DFT_MOD 目录下的源文件（dft_mod.asm 和dft_mod.ldf）添加到工程中，同时也将备用的数据文件square64.dat、sin64.dat、trig64.dat、noise64.dat 以及初始化文件zeros64.dat 加入工程中。2）选择输入数据文件编译程序前，在下述程序行中将输入数据文件改为square64.dat：.VAR input[N]= square64.dat;其中square6