数字信号处理上机实验2.docVIP

数字信号处理实验报告 班 级：信息25 姓 名：赵梦然 学 号：2120502123 频率采样型滤波器 一．实验目的 1. 学习使用频率采样型结构实现 FIR 滤波器，初步熟悉FIR 滤 波器的线性相位特点。 2. 直观体会频率采样型滤波器所具有的“滤波器组”特性，即 在并联结构的每条支路上可以分别得到输入信号的各次谐波。 3.学习如何使用周期冲激串检测所实现滤波器的频域响应。 实验内容 频率采样型滤波器是由一个梳状滤波器和若干路谐振器构成的，可用公式表 述如下： 其中 r 值理论上为1，实际中取非常接近1 的值。 为了使系数为实数，可以将谐振器的共轭复根合并，不失一般性，假设N 为 偶数，于是可以得到如图1 所示的结构。 其中，，。 以下实验中假设频率采样型滤波器阶数。 1.构造滤波器输入信号，其中，。基波频率，，，，，，，，。 设时域信号的采样频率，绘制出采样时刻从0到的采样信号波形，其中采样点数为，确认时域信号采样正确。 2.对采样信号的第二个周期，进行离散傅里叶变换，画出幅频特性和相频特性图，观察并分析其特点。 3.设，，，，，，计算滤波器抽头系数画出该滤波器的频谱图，观察并分析其幅频特性和相频特性。 4.编程实现图1所示的频率采样型滤波器结构，其中,取第3步中的值。为了简化编程，梳状滤波器可以调用CombFilter.m，谐振器可以调用Resonator2.m，使用help CombFilter和help Resonatoe2查看如何配置参数。将第1步生成的采样信号通过该滤波器，画出输出信号第二个周期的时域波形和频谱，并与第2步的频谱图进行对比，观察并分析二者的区别。 5.分别画出图1中前4路谐振器的输出信号第二个周期的时域波形，观察并分析输出信号的特点。 6.将输入信号换成周期为的冲激串，画出输出信号第二个周期的幅频特性，并与第3步的滤波器幅频特性进行对比，观察并分析二者的关系。 7.思考并回答下列问题 （1）在第2步的幅频特性中，各次谐波的幅度与相应的时域信号幅度有什么关系？ （2）实验中为什么要观察第二个周期，如果直接观察第一个周期会怎么样？ （3）如果取，观察会出现什么情况。 实验报告要求 1. 按照实验内容编写 Matlab 程序，给出运行结果（图），并逐项进行分析讨 论。 2. 提交完整的 Matlab 源程序。 3. 回答思考题，撰写实验报告。 四．实验过程及结果分析 第一问： MATLAB代码： %%第一问 f0=50; N=16; fs=N*f0; L=2*N; T=1/fs; t=0:T:(L-1)*T; s=0.5*cos(2*pi*0*f0*t+0)+1*cos(2*pi*1*f0*t+pi/2)+0.5*cos(2*pi*2*f0*t+pi)+2*cos(2*pi*3*f0*t-pi/2); %%作图 i=0:1:L-1; stem(i,s(i+1),.); title(时域信号采样); 运行结果： 第二问： MATLAB代码： %%第二问 f0=50; N=16; fs=N*f0; L=2*N;%%为了取第二个周期 T=1/fs; t=N*T:T:(L-1)*T; %%取采样信号的第二个周期 s=0.5*cos(2*pi*0*f0*t+0)+1*cos(2*pi*1*f0*t+pi/2)+0.5*cos(2*pi*2*f0*t+pi)+2*cos(2*pi*3*f0*t-pi/2); a=fft(s); %%作FFT y=abs(a); %%求其幅频特性 ang=angle(a); %%求其相频特性 i=0:1:N-1; subplot(2,1,1); %%subplot（m,n,p）,m表示是图排成m行，n表示图排成n列，p表示图所在的位置.是将多个图画到一个平面上的工具 stem(i,y(i+1),.); title(幅频特性图); subplot(2,1,2); stem(i,ang(i+1),.); title(相频特性图); 运行结果： 第三问： MATLAB代码： %%第三问 N=16; H=[1 exp(-j*pi*(N-1)/N) exp((-j*2*pi*(N-1)/N)) 0 0 0 0 0 0 0 0 0 0 0 -exp((-j*14*pi*(N-1)/N)) -exp((-j*15*pi*(N-1)/N))]; %%已知的H（n），表示成一个矩阵。 h=ifft(H); %%对H(n)求IFFT h %%输出抽头系数h %%画频谱图 y=abs(H); ang=angle(H); i=0:1:N-1