实验三FIR滤波器设计与实现.docx

电子科技大学电子工程学院标准实验报告（实验）课程名称 DSP技术电子科技大学教务处制表电子科技大学实验报告学生姓名：学号：指导教师：实验地点：实验时间：一、实验室名称：DSP技术实验室二、实验项目名称：FIR滤波器设计与实现三、实验学时：4四、实验目的：熟悉BF609开发板WL-BF609-EDU硬件平台。熟悉CCES开发软件平台的使用，掌握CCES集成开发环境的基本操作和常用功能，掌握CCES工程的创建、程序编写、编译和调试。掌握DSP中FIR滤波器设计、实现的方法。五、实验内容：了解BF609开发板WL-BF609-EDU。熟悉CCES集成开发环境的基本操作和常用功能。学习实验指导书中的低通滤波器设计与实现，验证滤波效果。。独立设计、实现FIR高通滤波器，并验证滤波效果。六、实验环境：预装开发环境Cross Core Embedded Studio 1.0.2的计算机。BF609开发板一套。ADDS HPUSB-ICE仿真器一套。实验步骤：输入条件：1MHz的点频信号，峰值为1；10M的点频信号，峰值为0.5；采样时钟40MHz。滤波器：17阶低通滤波器。输出：保留10MHz的点频信号，滤除1M的点频信号。1.用MATLAB设计FIR高通滤波器FIR滤波器原理有M个权系数（抽头）的FIR滤波器，如下图所示。滤波器的输入为随机过程，输出为其中，表示横向滤波器的权系数。图M抽头的FIR滤波器定义输入信号向量和权向量分别为则输出可表示为2.FIR滤波器高通滤波器实现(学习实验指导书中的低通滤波实验，独立完成高通滤波的DSP实现)生成Rb1.dat中数据，即滤波器权系数的m代码：close all;clear all;clc;f1=1e6;T1=1/f1;f2=10e6;T2=1/f2;T=1/40e6;t=0:T:1.6e-6-T;s1=cos(2*pi*f1*t);s2=0.5*cos(2*pi*f2*t);Rs1=round(s1*64);Rs2=round(s2*64);fid = fopen(Rs1.dat,w);fprintf(fid,%g\t,,Rs1);fclose(fid)fid = fopen(Rs2.dat,w);fprintf(fid,%g\t,,Rs2);fclose(fid)n=16;wn=0.35;b=fir1(n,wn,high);freqz(b,1);Rb=round(b*1024);fid = fopen(Rb1.dat,w);fprintf(fid,%g\t,,Rb);fclose(fid)结果：Rb1.dat中的数据为-2,-5,-4,17,42,12,-114,-280,666,-280,-114,12,42,17,-4,-5,-2,八、FIR高通滤波器代码输入条件：1MHz的点频信号，峰值为1；10M的点频信号，峰值为0.5；采样时钟40MHz。滤波器：17阶低通滤波器。输出：保留10MHz的点频信号，滤除1M的点频信号。1、MATLAB仿真代码closeall;clearall;clc;f1=1e6;T1=1/f1;f2=10e6;T2=1/f2;T=1/40e6;t=0:T:1.6e-6-T;s1=cos(2*pi*f1*t);s2=0.5*cos(2*pi*f2*t);n=16;wn=0.2;b=fir1(n,wn,high);freqz(b,1);Rs1=round(s1*64);Rs2=round(s2*64);Rb=round(b*1024);subplot(4,3,1); plot(t,s1,b);grid on;holdon;title(s1);subplot(4,3,2); plot(t,s2,k);grid on;holdon;title(s2);subplot(4,3,3); plot(t,s1+s2,r);grid on;holdon;title(s1+s2);subplot(4,3,4); plot(abs(fft(s1)),b);grid on;holdon;title(fft(s1));subplot(4,3,5); plot(abs(fft(s2)),k);grid on;holdon;title(fft(s2));subplot(4,3,6); plot(abs(fft(s1+s2)),r);grid on;holdon;title(fft(s1+s2));subplot(4,3,7); plot(abs(fft(conv(s1,b))),b);grid on;holdon;title(conv(s1,b));subplot(4,3,8); plot(abs(fft(conv(s2,b))),k);grid on;holdon;title(conv(