数字信号处理论文 数字信号处理实验报告

　实验一 快速傅里叶变换(FFT)及其应用

　一、 实验目的在理论学习的基础上，通过本实验，加深对 FFT的理解，熟悉 MATLAB^的有关函数。

　熟悉应用FFT对典型信号进行频谱分析的方法。

　了解应用FFT进行信号频谱分析过程中可能出现的问题，以便在实际中正确应用 FFT。

　熟悉应用FFT实现两个序列的线性卷积和相关的方法。

　二、 实验类型

　演示型

　三、 实验仪器

　装有MATLA爵言的计算机

　四、 实验原理

　在各种信号序列中，有限长序列信号处理占有很重要地位，对有限长序列，我们可以

　使用离散Fouier变换(DFT)。这一变换不但可以很好的反映序列的频谱特性， 而且易于用快

　速算法在计算机上实现，当序列 x(n)的长度为N时，它的DFT定义为：

　JV-1 $生

　反变换为：

　如-器冃吋

　科—

　有限长序列的 DFT是其Z变换在单位圆上的等距采样， 或者说是序列Fourier变换的等

　距采样，因此可以用于序列的谱分析。

　FFT并不是与DFT不同的另一种变换，而是为了减少 DFT运算次数的一种快速算法。它

　是对变换式进行一次次分解，使其成为若干小点数的组合，从而减少运算量。常用的 FFT

　是以2为基数的，其长度 A - o它的效率高，程序简单，使用非常方便，当要变换的序列长度不等于2的整数次方时，为了使用以2为基数的FFT，可以用末位补零的方法， 使

　其长度延长至2的整数次方。

　(一)在运用DFT进行频谱分析的过程中可能的产生三种误差

　1 .混叠

　序列的频谱是被采样信号频谱的周期延拓，当采样速率不满足 Nyquist定理时，就会

　发生频谱混叠，使得采样后的信号序列频谱不能真实的反映原信号的频谱。 避免混叠现象的唯一方法是保证采样速率足够高， 使频谱混叠现象不致出现， 即在确定采样频率之前， 必须

　对频谱的性质有所了解，在一般情况下，为了保证高于折叠频率的分量不会出现， 在采样前，

　先用低通模拟滤波器对信号进行滤波。

　2?泄漏

　实际中我们往往用截短的序列来近似很长的甚至是无限长的序列，这样可以使用较短 的DFT来对信号进行频谱分析， 这种截短等价于给原信号序列乘以一个矩形窗函数， 也相当

　于在频域将信号的频谱和矩形窗函数的频谱卷积，所得的频谱是原序列频谱的扩展。

　泄漏不能与混叠完全分开， 因为泄漏导致频谱的扩展， 从而造成混叠。为了减少泄漏的影响，可以选择适当的窗函数使频谱的扩散减至最小。

　3 .栅栏效应

　DFT是对单位圆上Z变换的均匀采样，所以它不可能将频谱视为一个连续函数， 就一定

　意义上看，用DFT来观察频谱就好像通过一个栅栏来观看一个图景一样， 只能在离散点上看

　到真实的频谱，这样就有可能发生一些频谱的峰点或谷点被“尖桩的栅栏” 所拦住，不能别

　我们观察到。

　减小栅栏效应的一个方法就是借助于在原序列的末端填补一些零值，从而变动 DFT的点数，这一方法实际上是人为地改变了对真实频谱采样的点数和位置， 相当于搬动了每一根

　“尖桩栅栏”的位置，从而使得频谱的峰点或谷点暴露出来。

　(二)用FFT计算线性卷积

　用FFT可以实现两个序列的圆周卷积。在一定的条件下，可以使圆周卷积等于线性卷 积。一般情况，设两个序列的长度分别为 N和2,要使圆周卷积等于线性卷积的充要条件

　是FFT的长度NA N+ N2-1，对于长度不足 N的两个序列，分别将他们补零延长到 No

　当两个序列中有一个序列比较长的时候，我们可以采用分段卷积的方法。有两种方法：

　.重叠相加法

　将长序列分成与短序列相仿的片段，分别用 FFT对它们作线性卷积，再将分段卷积各

　段重叠的部分相加构成总的卷积输出。

　.重叠保留法

　这种方法在长序列分段时， 段与段之间保留有互相重叠的部分， 在构成总的卷积输出时

　只需将各段线性卷积部分直接连接起来，省掉了输出段的直接相加。

　五、实验内容和要求

　1、一个连续信号含两个频率分量，经采样得

　x(n) sin[2 0.125n] cos[2 (0.125 f)n] n 0,1, ,N 1

　已知N 16, f分别为1/16和1/64 ,观察其频谱；当 N 128时，f不变，其结果有 何不同，为什么？

　代码：

　N=16;

　n=0:N-1;

　Df=1/16;

　Xn=si n( 2*pi*0.125* n)+cos(2*pi*(0.125+Df)* n); Xk=fft(X n,N);

　subplot(245);

　stem( n,abs(Xk));

　xlabel(' n');ylabel('X(k)');title('N=16,Df=1/16, subplot(241);

　stem( n,abs(X n));

　xlabel(' n');ylabel('X( n)');title('N=16,Df=1/16, Df=1/64;

　Xn=si n( 2*pi*0.125* n)+cos(2*pi*(0.125+Df)* n); Xk=fft(X n,N);

　subplot(246);

　stem( n,abs(Xk));

　xlabel(' n');ylabel('X(k)');title('N=16,Df=1/64, subplot(242);

　stem( n,abs(X n));

　xlabel(' n');ylabel('X( n)');title('N=16,Df=1/64,

　N=128;

　n=0:N-1;

　Df=1/16;

　Xn=si n( 2*pi*0.125* n)+cos(2*pi*(0.125+Df)* n);

　Xk=fft(X n,N);

　subplot(247);

　stem( n,abs(Xk));

　xlabel(' n');ylabel('X(k)');title('N=128,Df=1/16, subplot(243);

　stem( n,abs(X n));

　xlabel(' n');ylabel('X( n)');title('N=128,Df=1/16, Df=1/64;

　Xn=si n( 2*pi*0.125* n)+cos(2*pi*(0.125+Df)* n); Xk=fft(X n,N);

　subplot(248);

　stem( n,abs(Xk));

　xlabel(' n');ylabel('X(k)');title('N=128,Df=1/64, subplot(244);

　stem( n,abs(X n));

　频谱图');

　频谱图');

　时序图');

　频谱图');

　时序图');

　频谱图');

　时序图');

　频谱图');

　时序图');

　频谱图:

　2、h=1EDMfi4.fl7*￡、=口丄1=：.飞邓因N=1^QM風新鶴團N= 123,0^1^ 颉語區hW*ib殊駆40nou(n 6)

　2、

　h=1EDMfi4.fl7*￡

　、=口丄1=：.飞邓因

　N=1^QM風新鶴團

　N= 123,0^1^ 颉語區

　hW*ib殊駆

　40

　no

　u(n 6)的线性卷积和

　11 和 h(n)

　u(n)

　用FFT分别实现的x(n) 0.8n,0 n

　10点、20

　10点、20点圆周卷积，记录其波形，并说明他们之间的关系。

　代码：

　fun cti on [ y ] = circ onv( x1,x2,N )

　%UNTITLED3 Summary of this fun ctio n goes here

　% Detailed expla nati on goes here

　if len gth(x1)>N

　error('N should higher tha n or equal to the len gth of x1') end

　if len gth(x2)>N

　error(”)

　end

　x1= [x1,zeros(1,N-le ngth(x1))];

　x2=[x2,zeros(1,N-le ngth(x2))];

　m=0:1:N-1; x2=x2(mod(-m,N)+1); H=zeros(N,N);

　for n=1:1:N

　H(n ,:)=cirshift(x2, n-1,N);

　end

　y=x1*H';

　fun cti on [ y] = cirshift( x,m,N )

　%UNTITLED4 Summary of this fun ctio n goes here

　% Detailed expla nati on goes here

　if len gth(x)>N

　error(”)

　end

　x=[x zeros(1,N-le ngth(x))];

　n=0:1:N-1;

　n=mod( n-m,N);

　y=x( n+1);

　n=0:11;

　x=0.8.A n;

　h=[1,1,1,1,1,1]；

　N=17;

　x=[x,zeros(1,5)];

　h=[h,zeros(1,11)];

　X=fft(x);

　H=fft(h);

　Y=X.*H;

　y=ifft(Y);

　subplot(321);stem(x);xlabel(' n');ylabel('x');

　title('序列 x(n)');

　subplot(323);stem(h);xlabel(' n');ylabel('h');

　title('序列 h(n)');

　subplot(325);stem(y);xlabel(' n');ylabel('y');

　title(' 序列y(n),线性卷积');

　n=0:11;

　N=20;

　x=0.8.a n;

　h=[1,1,1,1,1,1];

　y=circ onv (x,h,N);

　x=[x,zeros(1,5)];

　h=[h,zeros(1,11)]; subplot(322);stem(x);xlabel(' n');ylabel('x'); title('序列 x(n)');

　subplot(324);stem(h);xlabel(' n');ylabel('h');

　title('序列 h(n)');

　subplot(326);stem(y);xlabel(' n');ylabel('y');

　title(' 序列y(n)，圆周卷积');

　频谱图:

　蔚 h]n)呼列h[n)痞列h?房列烛线性卷识库別」『)』區弔卷积

　蔚 h]n)

　呼列h[n)

　痞列h?

　房列烛线性卷识

　库別」『)』區弔卷积

　说明他们之间的关系？

　周期卷积是线性卷积的周期延拓。

　六、思考题

　实验中的信号序列xc( n)和xd( n)，在单位圆上的Z变换频谱Xc(ej )和Xd(ej ) 会相同吗？如果不同，说出哪一个低频分量更多一些，为什么？

　答:不同；其中，Xd(n)的低频分量较多，由图形可以看出，在低频处， Xd(n)的取值较

　多，呈递减趋势。

　对一个有限长序列进行 DFT等价于将该序列周期延拓后进行 DFS展开，因为 DFS

　也只是取其中一个周期来计算，所以FFT在一定条件下也可以用以分析周期信号序列。 如果

　实正弦信号sin(2 fn), f 0.1用16点FFT来做DFS运算，得到的频谱是信号本身的真实谱

　吗？为什么？

　答:可以把有限长非周期序列假设为一无限长周期序列的一个主直周期， 即对有限长非

　周期序列进行周期延拓，延拓后的序列完全可以采用 DFS进行处理，即采用复指数。

　实验二IIR数字滤波器的设计

　一、 实验目的掌握双线性变换法及脉冲响应不变法设计 IIR数字滤波器的具体设计方法及其原

　理，熟悉用双线性变换法及脉冲响应不变法设计低通、高通和带通 IIR数字滤波器的MATLAB 编程。

　观察双线性变换及脉冲响应不变法设计的滤波器的频域特性， 了解双线性变换法及

　脉冲响应不变法的特点。

　熟悉Butterworth滤波器、Chebyshev滤波器和椭圆滤波器的频率特性。

　二、 实验类型

　设计型

　三、 实验仪器

　装有MATLAB^言的计算机

　四、 实验原理

　1 .脉冲响应不变法

　用数字滤波器的单位脉冲响应序列 h(n)模仿模拟滤波器的冲激响应 ha(t),让h(n)正

　好等于ha(t)的采样值，即h(n) ha(nT)，其中T为采样间隔，如果以Ha(s)及H(z)分

　别表示ha(t)的拉式变换及h(n)的Z变换，则

　H(z)z esT

　H(z)

　z esT

　Ha(s j*m)

　2?双线性变换法

　S平面与z平面之间满足以下映射关系:

　jj ; z re1 Ts

　j

　j ; z re

　1 Is

　s平面的虚轴单值地映射于z平面的单位圆上，s平面的左半平面完全映射到 z

　s平面的虚轴单值地映射于

　圆内。双线性变换不存在混叠问题。

　双线性变换是一种非线性变换

　过预畸而得到校正。

　以低通数字滤波器为例，将设计步骤归纳如下：

　确定数字滤波器的性能指标：通带临界频率 fc、阻带临界频率fr、通带波动

　阻带内的最小衰减 At、采样周期T、采样频率fs；

　确定相应的数字角频率， c 2 fcT； r 2 frT；

　计算经过预畸的相应模拟低通原型的频率，

　2 C \ 2 r

　c 〒tg(-2)， r Ttg^2)；

　根据Q c和Q r计算模拟低通原型滤波器的阶数 N,并求得低通原型的传递函数

　Ha(S)；

　用上面的双线性变换公式代入 H a(s)，求出所设计的传递函数 H(Z)；

　分析滤波器特性，检查其指标是否满足要求。

　五、实验内容和要求

　1 ?用脉冲响应不变法设计一个巴特沃斯数字低通滤波器 ，要求在0-0.2 n内衰耗不

　大于3dB,在0.6 n -n内衰耗不小于 60dB,采样频率 Fs=500 Hz.

　代码：

　Omeg ap=0.2*pi*500;

　0megas=0.6*pi*500;

　rp=3;

　as=60;

　fs=500;

　wp=Omegap/fs;

　ws=Omegas/fs;

　[n ,Omegac]=buttord(Omegap,Omegas,rp,as,'s');

　[b,a]=butter( n,Omegac,'s');

　[bz, az]=impi nv ar(b,a,fs);

　w0=[wp,ws];

　hx=freqz(b z,az, w0);

　[h,w]=freqz(b z, az);

　plot(w*fs/(2*pi),abs(h));

　grid

　xlabel('频率 hz')

　ylabel('频率响应幅度')

　频谱图:

　100 150领率hz250

　100 150

　领率hz

　250

　分别用脉冲响应不变法和双线性变换法分别设计一个巴特沃斯滤波器数字低通滤波

　器,使其特性逼近一个巴特沃斯模拟滤波器的性能指标如下： 通带截止频率为2nX

　2000 (rad/s ),阻带截止频率为 2 nX 3000 (rad/s )，通带衰耗不大于 3dB,阻带衰 耗不小于15dB,采样频率Fs=100000 Hz,观察记录所设计数字滤波器的幅频特性曲 线，记录带宽和衰减量，检查是否满足要求。比较这两种方法的优缺点。

　代码(1):

　omeg ap=2*pi*2000;

　omegas=2*pi*3000;

　rp=3;

　as=15;

　fs=10000;

　wp=omegap/fs;

　ws=omegas/fs;

　[n,o megac]=buttord(omegap,omegas,rp,as,'s');

　[b,a]=butter( n,omegac,'s');

　[bz, az]=impi nv ar(b,a,fs);

　w0=[wp,ws];

　hx=freqz(bz,a z, w0);

　[h,w]=freqz(b z, az);

　plot(w*fs/(2*pi),abs(h));

　grid

　xlabel（'频率 hz'） ylabel（'频率响应幅度'）

　屏幕图（1）:

　In-JII5H E o.a fetB拦寻餅te

　In-

　JII

　5

　H E o.a fetB拦寻餅te

　Q

　50

　1

　o DO

　4D

　3 o

　30D

　N □ h 旺車 D co

　2(

　?rr

　代码(2):

　omeg ap=2*pi*2000;

　omegas=2*pi*3000;

　rp=3;

　as=15;

　fs=10000;

　wp=omegap/fs;

　ws=omegas/fs;

　omegap1=2*fs*ta n(wp/2);

　omegas 1=2*fs*ta n(ws/2);

　[n,o megac]=buttord(omegap1,omegas1,rp,as,'s');

　[b,a]=butter( n,omegac,'s');

　[bz, az]=bili near(b,a,fs);

　w0=[wp,ws];

　hx=freqz(bz,a z, w0);

　[h,w]=freqz(b z,az);

　plot(w*fs/(2*pi),abs(h));

　grid

　xlabel('频率 hz')

　ylabel('频率响应幅度')

　频谱图(2):

　￡QU 1

　￡QU 1UUU 1SLIU 2C100 250D 3UUU 3S0Q 4Q0D 45UU *kJUO

　-^^10642

　—□.□□.0.

　W鯉関昏褂虧

　比较优缺点：

　脉冲响应不变法的优点：1，模拟频率到数字频率的转换时线性的。

　2，数字滤波器单位脉冲响应的数字表示近似原型的模拟滤

　波器单位脉冲响应，因此时域特性逼近好

　缺点： 会产生频谱混叠现象，只适合带限滤波器

　双线性变换法优点： 克服多值映射得关系，可以消除频率的混叠

　缺点：是非线性的，在高频处有较大的失真。

　编写滤波器仿真程序，完成对实际采集的心电图信号序列 x(n)(具体数据见下面)

　的总响应序列y(n)，可直接调用 MATLAB filter 函数实现仿真。

　附：人体心电图采样信号在测量过程中往往受到工业高频干扰，所以，必须经过低通滤 波处理后，才能作为判断心脏功能的有用信息。下面的序列就是一个实际心电图信号采 样序列样本x(n)，其中存在高频干扰，实验时，将其作为输入信号，滤除其中的干扰 成分。

　4, 2, 0, 4, 6, 4, 2, 4, 6, 6, 4, 4, 6, 6'

　2, 6,12,8, 0, 16, 38, 60, 84, 90, 66, 32, 4, 2, x( n)

　4,8,12, 12, 10, 6, 6, 6, 4, 0, 0, 0, 0, 0, 2, 4, 0, 0, 0,

　2, 2, 0, 0, 2, 2, 2, 2,0

　代码:

　clear all;

　wp=0.2*pi;

　ws=0.3*pi;

　rp=1;

　rs=15;

　[n,wn]=buttord(wp/pi,ws/pi,rp,rs);

　[b,a]=butter( n,w n);

　N=0.5/(0.02);

　figure(1);

　freqz(b,a,N);

　grid;

　xlabel('频率 hz');

　ylabel('频率响应幅度');

　xn=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38, -60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2, -4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];

　figure (2);

　subplot(2,1,1);

　stem(x n,'.');

　title(' 心电图信号采样序列’)；

　yn=filter(b,a,x n);

　subplot(2,1,2)

　stem(y n,'.');

　title(' 滤波后的心电图信号');

　频谱图：

　40

　20

　-20

　-40

　20

　30

　六、思考题

　双线性变换法中 Q和3之间的关系是非线性的， 在实验中你注意到这种非线性关系了

　吗？从那几种数字滤波器的幅频特性曲线中可以观察到这种非线性关系？

　答：在双线性变换法中，模拟频率与数字频率不再是线性关系，所以一个线性相位模拟 器经过双线性变换后得到的数字滤波器不再保持原有的线性相位了。 如以上实验过

　程中，采用双线性变化法设计的 butter和chebyl数字滤波器，从图中可以看到这

　种非线性关系。

　2.能否利用公式

　完成脉冲响应不变法的数字滤波器设计？为什么？

　答：IIR数字滤波器的设计实际上是求解滤波器的系数和， 它是数学上的一种逼近问题，

　即在规定意义上（通常采用最小均方误差准则）去逼近系统的特性。如果在 S平面

　上去逼近，就得到模拟滤波器；如果在 z平面上去逼近，就得到数字滤波器。但是

　它的缺点是，存在频率混迭效应，故只适用于阻带的模拟滤波器。

　实验三FIR数字滤波器的设计

　、实验目的掌握用窗函数法，频率采样法及优化设计法设计 FIR滤波器的原理及方法，熟悉

　相应的MATLAB编程。

　熟悉线性相位FIR滤波器的幅频特性和相频特性。

　了解各种不同窗函数对滤波器性能的影响。

　二、 实验类型

　设计型

　三、 实验仪器

　装有MATLAB^言的计算机

　四、 实验原理

　线性相位实系数 FIR滤波器按其N值奇偶和h(n)的奇偶对称性分为四种： 1 . h(n)为偶对称，N为奇数

　丹J二皿(口) + ^空"功(兰丄4打)5雅丁

　2 歴=i 2

　H(ej 3的幅值关于 3 =0 n 2 n成偶对称。

　h(n)为偶对称，N为偶数

　H(ej 3的幅值关于 3 =成奇对称，不适合作高通。

　h(n)为奇对称，N为奇数

　\r_i丹牡)二辽胡(斗

　\r_i

　丹牡)二辽胡(斗

　H(ej 3的幅值关于 3 =0 n 2 n成奇对称，不适合作高通和低通。

　h(n)为奇对称，N为偶数

　冷)二［壬

　冷)二［壬2扯纠-1 +町血仙“

　泸1 2

　H(ej 3 3 =0 2n= 0，不适合作低通。

　FIR滤波器的常用设计方法:

　(一)窗口法

　窗函数法设计线性相位 FIR滤波器步骤

　确定数字滤波器的性能要求：临界频率 p、 st，滤波器单位脉冲响应长度 N。

　根据性能要求，合理选择单位脉冲响应 h(n)的奇偶对称性，从而确定理想频率响

　应H d (ej )的幅频特性和相频特性。

　求理想单位脉冲响应 hd(n)，在实际计算中，可对 Hd(ej )按M(M远大于N)点

　等距离采样，并对其求IDFT得hM (n)，用hM (n)代替hd(n)。

　选择适当的窗函数 w(n)，根据h(n) hd(n)w(n)求所需设计的fir滤波器单位脉 冲响应。

　求H(ej )，分析其幅频特性，若不满足要求，可适当改变窗函数形式或长度 N，

　重复上述设计过程，以得到满意的结果。

　窗函数的傅氏变换 W(ej )的主瓣决定了 H(ej )的过渡带宽。W(ej )的旁瓣大小和 多少决定了 H(ej )在通带和阻带范围内波动的幅度。

　常用的几种窗函数有：

　矩形窗 w(n)=R n(n)；

　2)373

　 TOC \o "1-5" \h \z Hanning 窗 b 匕 - i「丨- i ：; ; ；

　N」1

　Hamming 窗 - i'： -■ i - i ■ -1 ' —

　- 丁一— ?

　Blackmen 窗 ’

　Kaiser 窗

　h(0) 。

　式中Io(x)为零阶贝塞尔函数。

　(二)频率采样法

　频率采样法是从频域出发，将给定的理想频率响应 Hd(ej )加以等间隔采样

　瓦/严)1力厂臥⑹

　旳=—Jr

　N

　然后以此H d (k)作为实际FIR数字滤波器的频率特性的采样值 H(k)，即令

　丹@)=瓦(閒=兄(严)|加&7l"N—1

　q=——k

　N

　由H(k)通过IDFT可得有限长序列h(n)

　1 NT

　h(n) = —V丘閃总㈣的异=0入…川_1

　将上式代入到z变换中去可得

　Hg申已(片)

　Hg

　申已(片)

　其中①(3是内插函数五、实验内容及步骤27-1￡(斧)二工曰仗)0⑷-

　其中①(3是内插函数

　五、实验内容及步骤

　27-1

　￡(斧)二工曰仗)0⑷-

　i=0

　1 sin(^/2)

　M(<y)=

　N sin(d?/2)

　N=50,编程并画出矩形窗、汉宁窗、海明窗和布莱克曼窗的时域波形和归一化的幅

　度谱，并比较各自的主要特点。

　代码：

　wvtool(boxcar(50),ha nnin g(50),hammi ng(50),blackma n(50)));grid;

　频谱图：

　TIrrte dlornfiin Fr&qu^niCySamplesNoiirealized Fiequenc^ (xi: rad/^arnple)Factor Q %side? nhe

　TIrrte dlornfiin Fr&qu^niCy

　Samples

　Noiirealized Fiequenc^ (xi: rad/^arnple)

　Factor Q %

　side? nhe 乳怙rrnMcrr 阳 1 cJB M^m nhf? wicltl- i 3dH'i (:

　黑屋-dul<

　并比较各自的王要特点 ？

　矩形窗优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频干 扰和泄漏，甚至出现负谱现象； 汉宁窗主瓣加宽并降低， 旁瓣则显著减小，从减小泄漏观点

　出发，汉宁窗优于矩形窗?但汉宁窗主瓣加宽，相当于分析带宽加宽，频率分辨力下降；海 明窗的第一旁瓣衰减为一 42dB ?海明窗的频谱也是由 3个矩形时窗的频谱合成，但其旁瓣 衰减速度为20dB/ (10oct),这比汉宁窗衰减速度慢；布莱克曼窗的幅度函数主要由五部分 组成，他们的位移都不同，其幅度也是不同的 WRg (w)使旁瓣再进一步抵消。旁瓣峰值

　幅度进一步增加，其幅度谱主瓣宽度是矩形窗的 3倍。设计程序时用 backman函数调用。

　N=15，带通滤波器的两个通带边界分别是 1 0.3 , 2 0.5 。用汉宁(Hanning )

　窗设计此线性相位带通滤波器，观察它的实际 3dB和20dB带宽。N=45，重复这一设计，

　观察幅频和相位特性的变化，注意长度 N变化的影响。

　N=15汉宁代码：

　clc;clear all

　n=15;

　w1=0.3;w2=0.5;

　wn=[w1,w2];

　b2=fir1( n,wn ,ha nning(n+1));

　freqz(b2,1);

　title('汉宁窗，N=15');

　N=15汉宁频谱：

　-100*200 2o.a 0 4 as o.e o.FNonnallzed Frequency ( n raeWsample)OJl CLW无負p-*i

　-100

　*20

　0 2

　o.a 0 4 as o.e o.F

　Nonnallzed Frequency ( n raeWsample)

　OJl CLW

　无負p-*i常f

　N=45汉宁代码： clc;clear all n=45;w1=0.3;w2=0.5; wn=[w1,w2]; b2=fir1( n,wn ,ha nning(n+1)); freqz(b2,1);title('汉宁窗，

　N=45汉宁代码： clc;clear all n=45;

　w1=0.3;w2=0.5; wn=[w1,w2]; b2=fir1( n,wn ,ha nning(n+1)); freqz(b2,1);

　title('汉宁窗，N=45'); N=45汉宁频谱：

　50

　500

　-1AV4>」0VP 一II奶皿 Ud

　■50

　■100

　■150

　0 0.1 0.2 03 0.4 0.5 0.6 0 7 C.8 0.9 1

　Normalized Frequency ( l rad/samplei

　-1500

　0 01 02 0.3 04 0 5 0.6 07 05 09 1

　romalized FrequeiKy (-.n rad/samplei

　分别改用矩形窗和 Blackman窗，设计2中的带通滤波器，观察并记录窗函数对滤波 器幅频特性的影响，比较三种窗的特点。

　N=15矩形代码：

　clc;clear all

　n=15;

　w1=0.3;w2=0.5;

　wn=[w1,w2];

　b2=fir1( n,wn ,boxcar( n+1));

　freqz(b2,1);

　title('矩形窗，N=15');

　N=15矩形频谱:

　Lo oo10.1 n.2 0 3 0.4 0 5 0.6 0.7 n.B 0.9 1IMonrialized Frequency (xti

　Lo oo

　1

　0.1 n.2 0 3 0.4 0 5 0.6 0.7 n.B 0.9 1

　IMonrialized Frequency (xti rad/sarriple)

　oo

　so

　(5 養 Qp) uJSn

　N=45矩形代码：

　clc;clear all n=45;

　w1=0.3;w2=0.5;

　wn=[w1,w2];

　b2=fir1( n,wn ,boxcar( n+1)); freqz(b2,1);

　title('矩形窗，N=45');

　N=45矩形频谱：

　s-p)箸至ugw

　GOO

　c

　-Soo

　■iaoo

　T 1 II

　 i

　1 ?

　1 H

　1

　i *

　1 p

　1

　l

　■

　i

　i

　1

　J

　; f、

　1 II

　■ *■■■■』■ ■ ■■ ■ 0 ?iw ■ ■ ■ ■ ■ ■，

　i fc

　1 ?

　1 ■

　n p

　1 I*

　1 1

　…亠订寸护尺

　i ii

　1 H

　I b

　1 II

　i

　i

　i

　i

　i

　i

　i

　I

　j

　□ 0.1 0.2 0.3 D.4 U.S G E U一丿 0.3 0 y

　Moirn^hz Frtiquemcy (x-ji tad/Hanr-|ple)

　N=15Blackman 代码：

　clc;clear all

　n=15;

　w1=0.3;w2=0.5;

　wn=[w1,w2];

　b2=fir1( n,wn ,Blackma n(n+1)); freqz(b2,1);

　title('Blackman 窗，N=15');

　N=15Blackman 频谱：

　_蚩芋.-d厂护>.-叫a-50Blsjokmarim" P M=1

　_蚩芋.-d厂护>.

　-叫

　a

　-50

　Blsjokmarim" P M=1 5

　-1UU

　□ 1 Q.2 0.3 U 4 U.S D E □ 7 0.3 D 9

　Nuifri dlizud FruqtJLmuy (>.il iy J/burnplLi)

　0

　-GOU

　WOO

　□ n i o ? o 3 n 4 n 5 o R o 7 □ g □ s

　Nnrm FfpquFtn「Y (yn rfl flmplf1)

　N=45Blackman 代码:

　clc;clear all n=45;

　w1=0.3;w2=0.5;

　wn=[w1,w2];

　b2=fir1( n,wn ,Blackma n(n+1));

　freqz(b2,1);

　title('Blackman 窗，N=45');

　BlackmanW'p N=45-150^ rt

　BlackmanW'p N=45

　-150^ rt … rte 厂時

　□ 口/I U.2 U 3 口虫 U 5 0.6 0.7 0 3 □ 9 1

　Norrnahzed Frequency (xti rad/sample)

　N=45Blackman 频谱:

　D

　-50

　「1曲

　600

　-2000

　0

　-500

　*1000

　d50O

　n 1

　o 2

　o m

　n a

　o ￡

　OF

　a 7

　o s

　i

　□

　Norsmalized F requency (xti rad/^^mple)

　比较三种窗的特点？

　矩形窗优点是主瓣比较集中，缺点是旁瓣较高，并有负旁瓣，导致变换中带进了高频干

　扰和泄漏，甚至出现负谱现象；汉宁窗主瓣加宽并降低，旁瓣则显著减小，从减小泄漏 观点出发，汉宁窗优于矩形窗. 但汉宁窗主瓣加宽，相当于分析带宽加宽， 频率分辨力

　下降；布莱克曼窗的幅度函数主要由五部分组成， 他们的位移都不同， 其幅度也是不同

　的WRg ( w)使旁瓣再进一步抵消。旁瓣峰值幅度进一步增加，其幅度谱主瓣宽度是 矩形窗的3倍。设计程序时用 backman函数调用

　六、思考题

　1.定性地说明用本实验程序设计的 FIR滤波器的3dB截止频率在什么位置？它等于理

　想频率响应H d (ej )的截止频率吗？

　答：3dB截止频率在最高幅度下降一半是所对应的频率。因为设计的时候总会有误差，

　所以并不等于理想 Hd(ejw)的截止频率。

　2.如果没有给定h(n)的长度N,而是给定了通带边缘截止频率 p和阻带临界频率 戎

　以及相应的衰减，能根据这些条件用窗函数法设计线性相位 FIR低通滤波器吗？

　答：可以，可以用雷米兹交替法来设计。

相关关键词： 操作系统文件管理实验 新晃思源实验学校：誓师大会，备战中考 新晃思源实验学校图片 新晃思源实验学校简介 “国培计划（2017）”英语科送教活动在白水县田家炳实验中学进行