用频率采样法设计IR高通滤波器.docVIP

西安邮电大学 数字信号处理课内上机报告 专业班级： 学生姓名： 学号(班内序号): 实验要求： 给定参数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）=IDFT[H(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(信号滤波后频域图); 实验截图:  指导教师评语： 实 验 成 绩: 指导(辅导)教师 : 说明： 模板以《生产实习》课程为例，因此红色字体部分为可更改部分，请根据所授课程的实际情况填写。 本实习报告封皮适用于通信与信息工程学院所有实践、实验类课程（如：《生产实习》、《科研训练》、《企业实习》、《下一代网络实验