第二章C语言算法实验

DSP实验第二章、C语言算法实验实验一卷积（Convolve）算法实验一、实验目的1、了解卷积算法的原理。2、学会编写卷积算法。二、实验原理若LTI系统对δ(n)的响应为h(n)，h(n)称为系统的单位脉冲响应，则由时不变特性可得出系统对δ(n-k)的响应为h(n-k)；离散时间LTI系统对输入x(n)的响应为，

此式称为卷积和，通常记为：

y(n)=x(n)∗h(n)开始K<=Length1-1i<=pOutput[k]=Output[k]+x[k-1]*y[i],i++P++,P<=Length1-1K++P=Length1-2k>Length1-1k<=Length1-1I<=pOutput[k]=Output[k]+x[Length1-1-i]*y[iLength1-1-p+i],i++p-1END三、实验步骤1、熟悉卷积的基本原理；2、阅读所提供的样例实验程序；3、运行CCS软件，对样例程序进行跟踪，分析结果；4、填写实验报告；5、本实验所提供的样例子程序操作如下；1）源代码查看：用下拉菜单中Project/Open，打开”convolve.pjt”，双击”Source”可查看源程序，并加载”convolve.out”。2）在程序最后“i=0”处，设置断点；3）单击”Run”运行程序，程序运行到断点处停止。4）用下拉菜单中的View/Graph/Time/Frequency打开一个图形观察窗口；（设置该图形观察窗口的变量及参数，采用双踪观察两路输入变量Input及Impulse的波形，波形长度为80，数值类型为32位浮点型；

DisplayType：DualTime StartAddress-upperdisplay：Input StartAddress-lowerdisplay：Impulse AcquisitionBufferSize：80 DisplayDataSize：80

DSPDataType：32-bitfloatingpoint5）再用下拉菜单中的View/Graph/Time/Frequency打开一个图形观察窗口，以观察卷积结果波形，该观察窗口的参数设置为：变量为Output，长度为159，数据类型为32位浮点数；

DisplayType：SingleTime StartAddress：Output AcquisitionBufferSize：159 DisplayDataSize：159

DSPDataType：32-bitfloatingpoint6）调整观察窗口，观察两路输入波形和卷积结果波形；7）关闭“convolve.pjt”工程文件，关闭各窗口，实验结束；实验二相关（Correlation）算法实验一、实验目的1、了解相关算法。2、学习相关算法的实现方法。二、实验设备三、实验原理广义平稳随机信号x(n)和y(n)的相关函数的定义为：如果x(n),y(n)是各态遍历的，则上式的集合平均可以由单一样本序列的时间平均来实现，即：如果观察的点数N为有限值，则求r(m)估计值的一种方法是实际计算中，有偏估计为：无偏估计为：开始k=0K<length?L<length?r[length-1+k]=r[k]+x[j]*y[j+k]j++K++K=0K<length?L<length?r[length-1-k]=r[k]+x[j+k]*y[j]K++END四、实验步骤1、熟悉基本原理，阅读实验提供的程序；2、运行CCS软件，记录相关系数；3、填写实验报告；4、实验程序操作说明1）启动CCS2.0，用Project/open打开“cor01.pjt”工程文件：双击“cor01.pjt”及”Source”可查看各源程序：加载“cor01.out”：2）在程序最后“i=0”处，设置断点；3）单击“Run”运行程序，或按F5运行程序；程序将运行至断点处停止;用下拉菜单中的View/Graph的”Time/Frequency”打开一个图形观察窗口；采用双踪观察输入x_real及y_real的波形，长度为128，数值类型为32位浮点型；DisplayType：DualTime StartAddress-upperdisplay：x_real StartAddress-lowerdisplay：y_real AcquisitionBufferSize：128 DisplayDataSize：128

DSPDataType：32-bitfloatingpoint4）再打开一个图形窗口，以观察相关运算的结果，该观察窗口的参数设置为：变量为r_real，长度255，数据类型为32位浮点型；

DisplayType：SingleTime StartAddress：r_real AcquisitionBufferSize：255 DisplayDataSize：255

DSPDataType：32-bitfloatingpoint5）调整观察窗口，观察两路输入信号相关运算的结果；程序中，可赋0或1，赋0时，完成相关函数无偏估计的计算，赋1时，完成相关函数有偏估计的计算；x_real和y_real为参与相关运算的两路信号，当x_real=y_real时，完成自相关函数的计算，而当x_real不等于y_real时，完成互相关函数的计算。6）修改以上参数，进行“RebuildAll”,并重新加载程序，运行程序可以看到不同的实验结果。7）关闭“cor01.pjt”工程文件，关闭各窗口，实验结束。相关运算：x_real和y_real波形图Mode=0时，r_real波形图Mode=1时，r_real的波形图实验三FFT的DSP实现方法一、实验目的1、加深对DFT算法原理和基本性质的理解。2、熟悉FFT算法原理和FFT子程序的应用。3、学习用FFT对连续信号和时域信号进行谱分析，了解可能出现的分析误差及其原因。二、实验设备（略）四、实验内容和步骤1、复习DFT的定义，性质和用DFT作频谱分析的有关内容；2、复习FFT算法原理与编程思想，并对照DIT-FFT运算流程图和程序框图，了解本实验提供的FFT子程序；3、阅读本实验所提供的样例子程序；4、运行CCS软件，对样例程序进行跟踪，分析结果；记录必要的参数。5、填写实验报告。6、提供样例程序操作说明1)实验前准备在模拟信号产生单元中，一路信号源产生低频正弦信号（S1置“L”），另一路信号源产生高频正弦波信号（S11置“H”）,检查模拟信号输出端口”A””B”应短接；实验箱上电，用示波器分别观测out1和out2输出的模拟信号，调节电位器SPR1，SPR2（out1输出信号的频率调节和幅值调节）和电位器SPR1，SPR2(out2输出信号的频率调节和幅值调节)，直到满意为止；本样例实验程序建议：低频正弦波信号为1000Hz/1v;

高频正弦波信号为6KHz/1v；实验箱掉电，做以下连接和检查：短接输出端口“A”“B”；短接JAD1的INM,INMF;短接JAD2的INP,INPF;用导线连接out2（模拟信号输出）和JAD3脚（AD_IN）;检查JD是否断开。2）实验a）启动CCS2.0，用Project/open打开“ExpFFT01.pjt”工程文件：双击“ExpFFT01.pjt”及”Source”可查看各源程序：加载“ExpFFT01.out”；b）在主程序中，K++处设置断点；c）单击“Run”运行程序，或按F5运行程序；程序将运行至断点处停止：d）用View/Graph/Time/Frequency打开一个图形观察窗口；设置该观察图形窗口变量及参数；采用双踪观察在启动地址分别为0X3000h和0x3080h，长度为128的单元中数值的变化，数值类型为16位有效整型变量，这两段存储单元中分别存放的是经A/D转换后的输入信号和对该信号进行FFT变换后的结果；DisplayType：DualTime StartAddress-upperdisplay：0x3000 StartAddress-lowerdisplay：0x3080 AcquisitionBufferSize：128 DisplayDataSize：128

DSPDataType：16-bitsignedintegere）单击“Animate”运行程序，或按F10运行；调节观察窗口并观察变换结果；f）单击”Halt”暂停程序运行，关闭窗口，本实验结束；实验结果：在CCS2.0环境，同步观察输入信号波形及其FFT变换结果；初始化输入数组排序计算第一层中间值计算层数计算对应曾步长计算计算各层中间结果计算层数=0？计算结果输出实验四FIR滤波器的DSP实现方法一、实验目的1、掌握用窗函数法设计FIR数字滤波器的原理和方法。2、熟悉线性相位FIR数字滤波器特性。3、了解各种窗函数对滤波特性的影响。二、实验设备计算机，CCS2.0版软件，DSP硬件仿真器，实验箱。291．FIR滤波器基本概念1.FIR滤波器没有反馈回路，因此它是无条件稳定系统，其单位冲激响应h（n）是一个有限长序列。2.FIR滤波算法实际上是一种乘法累加运算。要点302.FIR滤波器中z-1的实现对于N级的FIR滤波器，在数据存储器中开辟一个称之为滑窗的N个单元的缓冲区，存放最新的N个输入样本；从最老的样本开始，每读一个样本后，将此样本向下移位，读完最后一个样本后，输入最新样本至缓冲区的顶部。用线性缓冲区实现z-1的优点是，新老数据在存储器中存放的位置直接明了。特点

(1)用线性缓冲区法实现z-1

31C54x片内没有I/O资源，CPU通过外部译码可以寻址64K的I/O单元。有两条实现输入和输出的指令：PORTRPA,Smem；将为PA的端口内容送；数据存储器SmemPORTWSmem,PA；将地址为Smem的数据；存储器内容送端口PA这两条指令至少需要2个字和2个机器周期。如果I/O设备是慢速器件，则需要插入等待状态。数据的输入/输出

32N=6的线性缓冲区存储器图33在数据存储器中开辟一个称之为滑窗的N个单元的缓冲区，滑窗中存放最新的N个输入样本；每次输入新样本时，以新样本改写滑窗中的最老的数据，而滑窗中的其它数据不作移动；利用片内BK（循环缓冲区长度）寄存器对滑窗进行间接寻址，循环缓冲区地址首尾相邻。利用循环缓冲区实现Z-1的优点是不需要移动数据，不存在一个机器周期中要求能一次读和一次写的数据存储器，因而可以将循环缓冲区定位在数据存储器的任何位置（线性缓冲区要求定位在DARAM，因为需要同时读和同时写）。

特点

（2）用循环缓冲区法实现z-134N=6的循环缓冲区存储器图初始化输入低通滤波器各频率参数构建滤波器频域特性波形调用FIR滤波器参数计算子程序计算各参数构建FIR滤波器进行滤波END四、实验内容和步骤1、复习如何设计FIR数字滤波。阅读本实验，掌握设计步骤；2、阅读本实验所提供的样例子程序；3、运行CCS软件，对样例子进行跟踪，分析结果；4、填写实验报告；1）实验前准备在模拟信号产生单元中，一路信号源产生低频正弦信号（S1置“L”），另一路信号源产生高频正弦波信号（S11置“H”）,检查模拟信号输出端口”A””B”应断开；实验箱上电，用示波器分别观测out1和out2输出的模拟信号，调节电位器SPR1，SPR2（out1输出信号的频率调节和幅值调节）和电位器SPR1，SPR2(out2输出信号的频率调节和幅值调节)，直到满意为止；本样例实验程序建议：低频正弦波信号为100Hz/1v;

高频正弦波信号为6KHz/1v；实验箱掉电，做以下连接和检查：短接输出端口“A”“B”；短接JAD1的INM，INMF;短接JAD2的INP，INPF;用导线连接out2（模拟信号输出）和JAD3脚（AD_IN）;检查JD是否断开。2）实验a）启动CCS2.0，用Project/open打开“FIRtestN01.pjt”工程文件：双击“FIRtestN01.pjt”及”Source”可查看各源程序：加载“FIRtestN01.out”：b）在主程序中，K++处设置断点：单击“Run”运行程序，或按F5运行程序；程序将运行至断点处停止：c）用View/Graph/Time/Frequency打开一个图形观察窗口；设置该观察图形窗口变量及参数；采用双踪观察在启动地址分别为0X3000h和0x3100h,长度为256的单元中数值的变化，数值类型为16位有效整型变量，这两段存储单元中分别存放的是经A/D转换后的输入信号和对该信号进行Fir变换后的结果；DisplayType：DualTime StartAddress-upperdisplay：0x3000 StartAddress-lowerdisplay：0x3100 AcquisitionBufferSize：256 DisplayDataSize：256

DSPDataType：16-bitsignedintegerd）单击“Animate”运行程序，或按F10运行；调节观察窗口并观察变换结果；e）单击”Halt”暂停程序运行，激活“FIRtestN01.c的编辑窗口；实验程序说明：该程序为51阶FIR低通滤波器算法程序，采用矩形函数实现，数组h和xmid长度均为51，fs为采样频率，fstop为滤波器截止频率，可以修改以上参数来改变滤波器性能。

f）重新“RebuildAll”后，并加载“load”，单击”Animate”.可得到不同的实验结果；实验五IIR滤波器的DSP实现方法一、实验目的1、熟悉设计IIR数字滤波器的原理与方法。2、掌握数字滤波器的计算机仿真方法。3、观察对实际信号的滤波作用，获得对数字滤波器的感性认识。二、实验设备计算机，CCS2.0版软件，DSP硬件仿真器，实验箱。三、实验原理1、IIR滤波器的设计原理，可以参见相关书籍。2、使用Matlab的FilterDeslgn&AnlysisTool（FDATool）工具设计IIR滤波器的参数，可以参见相关书籍。3、IIR算法的设计，可以参见相关书籍。4、IIR与FIR数字滤波器的比较1)在相同技术指标下，IIR滤波器由于存在着输出对输入的反馈，因而可用比FIR滤波器较少的阶数来满足指标的要求。2)FIR滤波器可得到严格的线性相位，而IIR滤波器选择性愈好，相位的非线性愈严重。3)FIR滤波器主要采用非递归结构，因而从理论上以及从实际的有限精度的运算中，都是稳定的。IIR滤波器必须采用递归结构，极点必须在z平面单位圆内，才能稳定，这种结构，运算中的四舍五入处理，有时会引起寄生振荡。四、实验内容和步骤1、复习有关巴特沃斯滤波器设计和用双线性变换法设计FIR数字滤波器的知识；2、阅读本实验所提供的样例子程序；3、运行CCS软件，对样例子程序进行跟踪，分析结果；4、填写实验报告；5、对样例程序实验操作说明1）实验前准备在模拟信号产生单元中，一路信号源产生低频正弦信号（S1置“L”），另一路信号源产生高频正弦波信号（S11置“H”）,检查模拟信号输出端口”A””B”应断开；实验箱上电，用示波器分别观测out1和out2输出的模拟信号，调节电位器SPR1，SPR2（out1输出信号的频率调节和幅值调节）和电位器SPR1,SPR2(out2输出信号的频率调节和幅值调节)，直到满意为止；本样例实验程序建议： 低频正弦波信号为100Hz/1v;高频正弦波信号为6KHz/1v；实验箱掉电，做以下连接和检查：短接输出端口“A”“B”；短接JAD1的INM，INMF；短接JAD2的INP，INPF；用导线连接out2（模拟信号输出）和JAD3脚（AD_IN）;检查JD是否断开。注：有关以上连接的说明，可参见第八章中语音接口跳线接输出接口配置使用说明以及信号产生单元配置正确完成计算机，DSP仿真器和实验箱的连接后，系统上电。2）实验a）启动CCS2.0，用Proje