数字信号处理实验二探究.docVIP

数字信号处理实验二 -频率采样型滤波器 姓名：周翔宇 班级：电信硕31 学号：2140508028 一 实验目的 1.1 学习使用频率采样型结构实现 FIR 滤波器，初步熟悉FIR 滤 波器的线性相位特点。 1.2 直观体会频率采样型滤波器所具有的“滤波器组”特性，即 在并联结构的每条支路上可以分别得到输入信号的各次谐波。 1.3 学习如何使用周期冲激串检测所实现滤波器的频域响应。 二 实验内容 频率采样型滤波器是由一个梳状滤波器和若干路谐振器构成的，可用公式表 述如下： 其中 r 值理论上为1，实际中取非常接近1 的值。 为了使系数为实数，可以将谐振器的共轭复根合并，不失一般性，假设N 为偶数，于是可以得到如图1 所示的结构。 其中，，。 以下实验中假设频率采样型滤波器阶数。 三．实验过程及结果分析 1. 构造滤波器输入信号，其中，。基波频率，，，，，，，，。设时域信号的采样频率，绘制出采样时刻从0到的采样信号波形，其中采样点数为，确认时域信号采样正确。 Matlab源码： A = [0.5 1 0.5 2]; Q = [0 pi/2 pi -pi/2]; N = 16; L = 2*N; f0 = 50; fs = N*f0; T = 1/fs; t = [0:T:(L-1)*T]; s = zeros(1,length(t)); for k=1:1:4 s = s + A(k)*cos(2*pi*(k-1)*f0*t+Q(k)); end k = 0:L-1 stem(k,s(k+1)); 运行结果： 实验结果分析： 信号由直流分量，一次谐波（基波），二次谐波，三次谐波组成；信号基波频率为频率为=f0*N，采样点数为*N；相当于对直流分采样了相等的；对（波采样了两个谐波采样四个周期，；对三次谐波了六个周期每个周期 2. 对采样信号的第二个周期，进行离散傅里叶变换，画出幅频特性和相频特性图，观察并分析其特点。 Matlab源代码： A = [0.5 1 0.5 2]; Q = [0 pi/2 pi -pi/2]; N = 16; L = 2*N; f0 = 50; fs = N*f0; T = 1/fs; t = [0:T:(L-1)*T]; s = zeros(1,length(t)); for k=1:1:4 s = s + A(k)*cos(2*pi*(k-1)*f0*t+Q(k)); end k = N:L-1; s_analys = s(k+1); y = fft(s_analys); show_module = abs(y); show_angle = angle(y); figure(1); stem(k,show_module(k-N+1)); figure(2); stem(k,show_angle(k-N+1)); 结果： 结果分析： 采样信号由直流分量，一次谐波，二次谐波，三次谐波组成；对应频域有四个冲激，依次为直流分量，一次谐波，二次谐波，三次谐波的冲激；这与他们的幅度矩阵[0.5 1 0.5 2]组成，由于一次谐波，二次谐波，三次谐波到频域冲激变为所以频域四种信号对应得幅度比为这与所的是一致的；相频特性其初始相位0 0.5pi pi -0.5pi]得，它与图中 3. ，，，，，，计算滤波器抽头系数画出该滤波器的频谱图，观察并分析其幅频特性和相频特性。 Matlab源代码： N=16; H = zeros(1,16); H(1) = 1; H(2) = exp(-1i*pi*(N-1)/N); H(3) = exp(-2i*pi*(N-1)/N); H(15) = exp(-14i*pi*(N-1)/N); H(16) = exp(-15i*pi*(N-1)/N); h = ifft(H) show_module_h = abs(H); show_angle_h = angle(H); figure(1); k= 0:15; stem(k,show_module_h(k+1)); figure(2); stem(k,show_angle_h(k+1)); w=0:pi/200:2*pi; h1=zeros(1,length(w)); for d=1:length(w) for n=1:N re=cos(w(d)*(n-1)); im=-sin(w(d)*(n-1)); e=complex(re,im); h1(d)=h1(d)+h(n)*e;