用频率采样法设计FIR高通滤波器.doc

此文档收集于网络，如有侵权，请联系网站删除实验总成绩： 装 订 线报告份数： 西安邮电大学 数字信号处理课内上机报告专业班级： 学生姓名： 学号(班内序号): 实验要求：给定参数N，利用频率采样法设计出相应FIR高通滤波器。实验目的：学习用频率采样法设计FIR高通滤波器的原理及其方法步骤；学习编写数字滤波器设计程序方法，并能进行正确编程；根据设计指标要求，给出设计步骤。实验代码：M=31;%所需频率采样点个数Wp=0.6*pi;%通带截止频率m=0:M/2;%阻频带上的采样点Wm=2*pi*m./(M+1);%阻带截止频率mtr=ceil(Wp*(M+1)/(2*pi);%向正方向舍入ceil(3.5)=4;ceil(-3.2)=-3;Ad=Wm=Wp;Ad(mtr)=0.28;Hd=Ad.*exp(-j*0.5*M*Wm);%构造频域采样向量H(k)Hd=Hd conj(fliplr(Hd(2:M/2+1);%fliplr函数实现矩阵的左右翻转conj是求复数的共轭h=real(ifft(Hd);%h（n）=IDFTH(k)w=linspace(0,pi,1000);%用于产生0,pi之间的1000点行矢量H=freqz(h,1,w);%滤波器的幅频特性图figure(1)plot(w/pi,20*log10(abs(H);%参数分别是归一化频率与幅值xlabel(归一化频率);ylabel(增益/分贝);axis(0 1 -50 5);f1=200;f2=700;f3=800%待滤波正弦信号频率fs=2000;%采样频率figure(2)subplot(211)t=0:1/fs:0.25;%定义时间范围和步长s=sin(2*pi*f1*t)+sin(2*pi*f2*t)+sin(2*pi*f3*t);%滤波前信号plot(t,s);%滤波前的信号图像xlabel(时间/秒);ylabel(幅度);title(信号滤波前时域图);subplot(212)Fs=fft(s,512);%将信号变换到频域AFs=abs(Fs);%信号频域图的幅值f=(0:255)*fs/512;%频率采样plot(f,AFs(1:256);%滤波前的信号频域图xlabel(频率/赫兹);ylabel(幅度);title(信号滤波前频域图);figure(3)sf=filter(h,1,s);%使用filter函数对信号进行滤波%输入的参数分别为滤波器系统函数的分子和分母多项式系数向量和待滤波信号输入subplot(211)plot(t,sf)%滤波后的信号图像xlabel(时间/秒);ylabel(幅度);title(信号滤波后时域图);axis(0.2 0.25 -2 2);%限定图像坐标范围subplot(212)Fsf=fft(sf,512);%滤波后的信号频域图AFsf=abs(Fsf);%信号频域图的幅值f=(0:255)*fs/512;%频率采样plot(f,AFsf(1:256)%滤波后的信号频域图xlabel(频率/赫兹);ylabel(幅度);title(信号滤波后频域图);实验截图:装 订 线指导教师评语： 实 验 成 绩: 指导(辅导)教师 : 说明：1. 模板以生产实习课程为例，因此红色字体部分为可更改部分，请根据所授课程的实际情况填写。2. 本实习报告封皮适用于通信与信息工程学院所有实践、实验类课程（如：生产实习、科研训练、企业实习、下一代网络实验、信息安全基础实验、移动通信基础实验、信息安全专业课程设计、通信工程专业课程设计、认识实习、通信软件设计课内实验、基于JAVA的网络编程、通信网络实验、移动通信实验等）。3. 封面中，时间一栏指的是学生提交报告的时间，如：2011年3月10日。4. 封底中，教师评语要求详尽、具体，能根据学生的实际情况填写，并能指出学生学习本课程中所