实验三FIR滤波器的DSP实现.docVIP

实验三 FIR滤波器的DSP实现 一、实验目的 1、熟悉数字滤波的基本原理和实现方法； 2、熟悉FIR数字滤波器的设计方法； 3、熟悉FIR滤波器的DSP实现方法； 二、实验设备 硬件：计算机、DSP实验箱 软件：WINDOWS操作系统、CCS软件和MATLAB软件 实验内容 1、录制一段音频信号，用实现FIR实现滤波过程； 2、在CCS中编写程序，实现（1）中设计的滤波器； 3、基本原理 X(t)为自己录制的一段音频信号（语音或者音乐），模数转换过程即为将模拟信号转换成数字信号的过程，数字信号处理器完成FIR滤波器的功能，滤波器的系数由MATLB设计实现。滤波后的数据存下来，通过数模转换播放出来。 本实验通过用DSP实现规定的FIR滤波器，使同学们了解FIR滤波器的设计方法，及DSP上FIR滤波器的实现方法。 设h(n)(n=0,1,2,3…，N-1)为滤波器的单位冲激响应序列，输入信号为x(n),则 滤波器就是要实现下列差分方程 假设滤波器的单位脉冲响应为h[0],h[1],h[N-1],x[k]为待滤波的输入信号，则滤波器的输出为： 进行滤波实际上就是计算序列的线性卷积。 四、实验步骤 1、FIR滤波器的MATALAB的设计 FIR滤波器的设计任务是选择有限长度的h(n)。使传输函数H()满足要求。FIR滤波器的设计方法有多种，如窗函数法、频率采样法及其它各种优化设计方法x=x(:,1); y = abs(fft(x(1:2048))); subplot(211);plot(x);title(原始信号); subplot(212);plot(y);title(频谱信号); s = round(32767* x); fid = fopen(D:\Windows XP 关机.dat,w); fprintf(fid,1651 1 0 0 100\n); for i=1:256 fprintf(fid,%d\n,s(i)); end fclose(fid); 2）启动滤波器设计分析器 在MATLAB的start菜单中选择Toolboxes-Filter Design-Filter Design Analysis Tools(fdatool)，或者在命令行中输入fdatool来启动滤波器设计分析器。启动后界面如图2.2所示： 图2.2： 3）滤波器设计 在选项中选择或输入滤波器参数，然后点击”Design Filter”按钮完成滤波器设计。设计成功后的结果如图2.3所示。 图2.3 4）导出FIR滤波器系数 在fdatool中，选择Targets-Code Composer Studio?IDE。 在出现的对话框中选择输出文件类型为C.header file,输出系数类型为signed 16-bit integer，如图2.4所示。 点击Generate按钮，选择路径，即可输出前一步设计出的FIR滤波器系数表。 图2.4 2、滤波器的DSP实现 FIR是将待滤波的数据序列与滤波系数序列相乘后再相加运算,同时要模仿FIR结构中的延迟线将数据在存储器中滑动 以往的通用处理器或单片机在进行此类乘加运算的时候，都会耗费大量的CPU和总线资源DSP作为专用信号处理芯片，它在结构上针对数字信号处理的特点进行了改进和优化,并且增加了特殊的指令专门用于数字处理。 设计好FIR滤波器后，如何写入到DSP中，使DSP按照设计要求实现低通滤波的效果，本实验的滤波算法是直接从TMS320C54x的DSPDIB库中调用FIR滤波函数。DSPLIB是一个为C语言程序员开发TMS320C54x而建立的经过优化的DSP函数库。在工程中要添加头文件“*.h ”，在C源程序中要包含头文件。 1）数据的载入 头文件的更改 将头文件中的11个点的数据复制到test.dat文件中 载入语音数据 2）信号的输出 观察输入输出信号 生成输入输出图像 观察频域输入图像 观察频域输入图像 3、滤波后语音信号的读取 在MATALAB中将之前滤波后的数据转换为语音信号并试听，比较其与滤波前的语音信号的差别。 代码如下： fid = fopen(D:\Windows XP 关机.dat); y=fscanf(fid,%x); y=y/32768; wavwrite(y,22050,16,D:\Windows XP 关机后.wav); 2