DSP应用系统设计综合设计报告。fir滤波器.doc

《DSP应用系统设计》综合设计报告 题目： 姓名 ＿ 学号 ＿ ＿ 手机 Email 专 业 班级 ＿ ＿ 开课学期 至 学年＿ 学期 提交时间 年 月 日 摘 要一、设计要求 通带最大衰减 阻带边界频率 阻带最小衰减 二、系统方案 三、算法设计 采用MATLAB设计FIR数字滤波器的方法有多种，本实验采用fir1函数b=fir1(N,Wn)或b = fir1(N,Wn,low) 其中N+1为滤波器阶数，b为返回的长度为N+1的系数向量。Wn为归一化边界频率，0 Wn 1.0，Wn 处增益为-6dB。Wn =1.0对应的非归一化频率为Fs/2，其中Fs为采样频率。 例，设Fs=5000Hz，低通滤波器的边界频率fp=800Hz，则对应的归一化数字角频率为Wn=fp/(Fs/2)=0.32，取N=63（对应的滤波器阶数为64）。于是，该低通fir滤波器的设计指令为b=fir1(63,0.32)。执行该指令得到的系数为： b=[ 0.0002 , -0.0006 , -0.0009 , -0.0004 , 0.0008 , 0.0016 , 0.0010 , -0.0012 , -0.0030 , -0.0021 , 0.0016 , 0.0050 , 0.0042 , -0.0018, -0.0080 -0.0075 , 0.0014 , 0.0119 , 0.0128 , -0.0000 , -0.0171 , -0.0211 , -0.0034 0.0244 , 0.0354 , 0.0111 , -0.0370 , -0.0664, -0.0326 , 0.0739 , 0.2109 , 0.3069 , 0.3069 , 0.2109 , 0.0739 -0.0326 -0.0664 , -0.0370 , 0.0111 , 0.0354 , 0.0244 , -0.0034,-0.0211 , -0.0171 , -0.0000 , 0.0128 , 0.0119 , 0.0014 , -0.0075 ,-0.0080 , -0.0018 , 0.0042 , 0.0050 , 0.0016 , -0.0021 , -0.0030,-0.0012 , 0.0010 , 0.0016 , 0.0008, -0.0004 , -0.0009 , -0.0006 , 0.0002] 注意: 上述滤波器系数为浮点格式，当DSP程序中采用整数格式时。滤波器系数需相应地变为整数，方法为乘以32768并取整。 b=[7,-19,-31,-14,26,53,32,-39,-97,-69,52,165,136,-58,-261,-247,47,389,420,0,-559,-693,-111,799,1162,362,-1212,-2176,-1068,2420,6912,10056,10056,6912,2420,-1068,-2176,-1212,362,1162,799,-111,-693,-559,0,420,389,47,-247,-261,-58,136,165,52,-69,-97,-39,32,53,26,-14,-31,-19,7] （2）设计高通滤波器，格式为：b = fir1(N,Wn,high) 参数意义同上。 （3）设计带通滤波器，格式为：b = fir1(N,Wn)或b = fir1(N,Wn, bandpass) 其中Wn = [W1 W2], W1、W2为边界频率。 （4）设计带通滤波器，格式为：b = fir1(N,Wn)或b = fir1(N,Wn, stop) 其中Wn = [W1 W2], W1、W2为边界频率。 采用函数plot(b)可以绘制滤波器的冲击响应图： 采用函数freqz(b)可以绘制滤波器的频率响应图： 四、软件实现 .c为主程序，用于设置滤波器、测试信号参数，产生测试信号，调用fir滤波程序。 （2）fir.asm为采用汇编语言编写的fir滤波模块。 （3）fir1.cmd为命令文件。 （1）主程序 #include math.h #define Fs 8000 /* 采样频率 */ #define T 1/Fs /* 采样时间 */ #define f1 500 /* 正弦信号1频率 */