北京理工大学信与系统实验报告信的频域分析

1、实验3信号的频域分析(综合型实验)一、实验目的1) 深入理解信号频谱的概念，掌握信号的频域分析方法。2) 观察典型周期信号和非周期信号的频谱，掌握其频谱特性。二、实验原理与方法1连续周期信号的频谱分析如果周期信号满足Dirichlet条件，就可展开为傅里叶级数的形式，即x(t)ckejk 0t ( 1)kck x(t)e jk 0dt ( 2)T0 To其中To表示基波周期,o 2 /To为基波频率,(.)表示任一个基波周期内的积分。ToCk成为x(t)的傅里叶系数。周期信上面两式为周期信号复指数形式的傅里叶级数，系数号的傅里叶级数还可由三角函数的线性组合来表示，即x(t) aoak cosk

2、 otk 1kgsink ot ( 3)1其中ao丄x(t)dt,ak 刍x(t)coskotdt, bkx(t)sink otdt ( 4)ToTTo TrnTo T(3)式中同频率的正弦、余弦项合并可以得到三角函数形式的傅里叶级数，即x(t) AoAk cos(k ot k)(5)k 1其中 Aoao, Ak .a2 bk, k arcta(6)ak任何满足Dirichlet条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号，但在实际应用中常采用有限项级数代替，所选级数项越多就越接近原信号。2连续非周期信号的频谱分析

3、对于非周期连续时间信号，信号的傅里叶变换和傅里叶逆变换定义为.t 1 . tX( ) x(t)e j tdt (7)x(t)X( )ej td (8)以上两式把信号的时频特性联系起来，确立了非周期信号x(t)和频谱X()之间的关系。利用MATLAB可以方便地求出非周期连续时间信号的傅里叶变换，几种常见方法如下：1) 符号运算法fourierMATLAB的符号数学工具箱提供了直接求解傅里叶变换和反变换的函数, 函数和ifourier函数，基本调用格式为X fourier (x)x ifourier (X)默认的时域变量为t，频域变量为。例：求x(t) e2|t|的傅里叶变换，代码及运行结果如下：

4、>> syms t>> x=exp(-2*abs(t);>> X=fourier(x)4/(wA2 + 4)所以傅里叶变换结果为X()也可利用int函数直接根据式(7)求傅里叶变换。2) 数值积分法除了采用符号运算的方法外，还可以采用 MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析。quad函数是一个用来计算数值积分的函数。利用quad函数可以计算非周期连续时间信号的频谱。Quad函数的一般调用格式为：y=quad(fu n, a,b)y=quad(fu n,a,b,TOL,TRACE,p1,p2,.)其中fun指定被积函数，可以用i

5、nline命令来创建，也可通过传递函数句柄的形式 来指定，a、b表示定积分的上下限，TOL表示允许的相对或绝对积分误差，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值，如果TOL和TRACE为空矩阵，则使用缺省值，“p1，p2,.”表示被积函数除时间t之外所需的其他额外输 入参数。3) 数值近似法还可以利用 MATLAB的数值计算的方法近似计算连续时间傅里叶变换。傅里叶变换X()可以由(9)式近似计算X()x(t)e j t dtlim0kx(k )e j k(9)当x(t)为时限信号，且足够小，则(9)可以演变成bX( ) x(k )e jk( 10)k a(10)中求和部分又可以

6、表示成一个行向量和一个列向量的乘积jab . e x(k )e jk x(a ),x(a 1) ),.,x(b )k a.j(a 1)(11)jb上式可以很方便地利用MATLAB实现。3离散周期信号的分频域分析基波周期为N的周期序列x(n)可以用N个成谐波关系的复指数序列的加权和表示，即x(n)Ckejk(2 /N)n( 12)k N这里k=<N>表示求和仅需包括一个周期内的N项，周期序列在一个周期内的求和与起点无关。将周期序列表示成式(12)的形式，称为离散傅里叶级数，而系数ck则称为傅里叶系数。离散傅里叶系数ck可由(13)式确定。1/、 jk(2 /N)n一 c、Ckx(n)

7、e j ()( 13)N k N傅里叶系数ck也称为x(n)的频谱系数，而且可以证明Ck是以N为周期的离散频率序列。 这说明了周期的离散时间函数对应于频域为周期的离散频率。这里，我们用周期 N与傅里叶系数ck的乘积来表示周期离散时间信号的频谱，即X(k) N c,x(n)ejk(2/N)n ( 14)k NX(k)可以利用MATLAB提供的函数fft用来计算，调用格式为X=fft(x)该函数返回X(k)一个周期内的值，其中x表示x(n)个周期内的样本值。4.离散非周期信号的频域分析非周期序列x(n)可以表示成一组复指数序列的连续和x(n)乂 X(ej 0 "d(15)其中 X(ej

8、) x(n)e j n ( 16)n式(16)成为x(n)的离散时间傅里叶变换，式(15)和(16)确立了非周期离散时间信号x(n)及其离散时间傅里叶变换X(ej )之间的变换。X(ej )是连续频率的函数，称为频谱函数，且 X(ej )是周期的连续频率函数，其周期为2 。可见，非周期离散时间函数对应于频域中是一个连续的周期频率函数。对于有限长的离散时间序列，(16)可以表示为jninNX(ej )x(n)en nix( ni),x( n2),.,x( nN)jn2jnN上式可以方便地利用 MATLAB实现。三、实验内容(1) 已知x(t)是如下图所示的周期矩形脉冲信号x(t)I-c/2c/2

9、T1) 计算该信号的傅里叶级数;八2A .a0 A， ak sin(k Tk得傅里叶级数为：x(t)2 AA sin(k ) T kTcos(当k 1 T2) 利用MATLAB绘出由前N次谐波合成的信号波形，观察随着N的变化合成信号波形的变化规律；用MATLAB编写代码如下：N=i nput('N=');A=i nput('A=');c=i nput('c=');T=i nput('T=');t=:;x=A*c*o nes(size(t)/T;for n=1:Nx=x+(2*A/( n*pi)*si n(n *pi*c/T)*co

10、s(2*pi* n*t/T);endplot(t,x);xlabel('Time(sec)' )ti tle('N=' num2str( n)在命令窗口输入以下语句:>> subplot(223)>> DFTthird_2_1N=20A=1c=T=1>> subplot(224)>> DFTthird_2_1N=40A=1c=T=1Time(sec)N=201A"0.50k-1.5-0.5-2-10 12Time(sec)N=10-2-1012Time(sec)1.510.50-0.5-2-1012Tim

11、e(sec)>> subplot(221) >> DFTthird_2_1 N=5A=1c=T=1>> subplot(222) >> DFTthird_2_1 N=10A=1c=T=1图形如下：由以上四个图可知，随着 N的增大，合成信号的波形越来越接近原脉冲信号。3）利用MATLAB绘出周期矩形脉冲信号的频谱， 观察参数T和 变化时对频谱波形的影响。 可计算出傅里叶系数为：代k 0TkAsin(k Pk画出该信号频谱 MATLAB代码如下N=i nput('N=');c=in put('c=');A=i nput

12、('A=');T=i nput('T=');n1=-N:-1;c1=(A./(n 1*pi).*si n(n 1*pi*c/T);c0=c*A/T;n2=1:N;c2=(A./(n 2*pi).*si n(n 2*pi*c/T);cn=c1 c0 c2;n=-N:N;subplot(211);stem (n, abs(c n),'filled' );xlabel( 'omega/omega_0' ) title( 'Mag nitude of ck');subplot(212);stem (n,an gle(c n

13、),'filled' );xlabel( 'omega/omega_0' ); title( 卩 hase of ck' )命令窗口：>> DFTthird_2_2N=20 c=A=1T=4图形:Magn itude of ck43210-20-15-10-50510/ 01520N=20, A=1,改变T和时的波形变化:c=0.5;T=8;c/T=0.0625Mk k o ec=0.5;T=4;c/T=0.125/00.40.30.20.1c=1;T=4;c/T=0.25-1001020/ 03 2 10 0 0kcfo eaulrDyaM

14、-200c=0.5;T=2;c/T=0.25-20 -10 0 10 20/ 0综合分析以上四个图像可得：频谱的波形与占空比 一有关，对于T或 取不同的值时，当一TT不变，频谱波形不变，当变大，频宽（第一个过零点频率）减小，过零点频率增多，谱 线变得稀疏，各频率分量振幅增大。 观察实验结果，思考如下问题：1 1.什么是吉伯斯现象？产生吉伯斯现象的原因是什么？答：吉伯斯现象：将具有不连续点的周期函数（如矩形脉冲）进行傅里叶级数展开后选 取有限项进行合成。在不连续点附近，部分和有起伏，其峰值最大值是不连续点处高度的倍。不连续点处级数收敛于左右极限的平均值，t愈接近不连续点时，为将误差减小至低于某一

15、给定值，N必须取得很大。随着 N的增大，部分和的起伏就向不连续点处压缩，但是对有限的N值，起伏的峰值大小不变。产生吉伯斯现象的原因：当一个信号通过某一系统时，如果这个信号是不连续时间 函数，则因为一般的物理系统对信号的高频分量都有衰减作用，所以产生了吉伯斯 现象。1 2.以周期矩形脉冲为例，说明周期信号的频谱有什么特点？由结果可知：1. 周期性矩形脉冲信号的频谱是离散的，其频谱仅存在于0,2 0,3 °,等离散值处，谱线间隔距离为固定值0,各次谐波的频率均为基频0的整数倍。2. 随着谐波次数增高，谱线长度逐渐趋于收敛。所以离散型、收敛性以及谐波性是周期信号的共同特点。1 3.周期矩形

16、脉冲信号的有效频带宽度与信号的时域宽度之间有什么关系？答：根据B 2 /T（ B为有效频宽，T为脉冲的时域宽度），有效频带宽度与信号的时域宽度成反比。1 4.随着矩形脉冲信号参数/T的变化，其频谱结构（如频谱包络形状、过零点、谱线间隔等）如何变化？答：随着一的变换，频谱包络形状不变。一增大时，过零点增多，谱线间隔逐渐变大。TT（2）已知x（t）是如下图所示的矩形脉冲信号。x（t）屮-c/c/21)求该信号的傅里叶变换；>> syms t c A>> x=A*heaviside(t+c/2)-heaviside(t-c/2);ans =>> X=fourier

17、(x);>> collect(X)(2*A*si n(c*w)/2)/w所以傅里叶变换为：2AX( )一sin(y)2)利用MATLAB绘出矩形脉冲信号的频谱，观察矩形脉冲宽度 变化时对频谱波形的影响; 代码：A=in put( 'A=');c=i nput( 'c=');syms t wX=i nt(A*(heaviside(t+c/2)-heaviside(t-c/2)*exp(-j*w*t),t,-1,1) ezplot(abs(X),-6*pi,6*pi);grid on;xlabel( 'omega' ); ylabel(&

18、#39;Mag ni tude');title( '|X(omega)|')图形如下:|X( )|令A=1，改变的大小:-10 0 10c=0.252 50 10.40.30.20.1-10 0 10c=0.5e c - ncoy M-10 0 10e cs -CBU Me cs -CBU Mc=21 *. r -10 0 10增大时，频谱振幅增大，有效频宽减小。3) 让矩形脉冲的面积始终等于1，改变矩形脉冲宽度，观察矩形脉冲信号时域波形和频谱随矩形脉冲宽度的变化趋势。syms t wc=i nput('c=');x=(1/c).*(heaviside(

19、t+c/2)-heavi side(t-c/2);X=in t(x*exp(-j*w*t),t,-c/2,c/2); subplot(211);ezplot(x,-c,c);xlabel( 't');title( 'x(t)');subplot(212);ezplot(abs(X),-6*pi,6*pi); grid on ;xlabel( 'omega' ); ylabel( 'Mag ni tude');title('|X(w)|')axis tight时域波形及频谱在不同脉冲宽度下的图像:c=c=x(t)c=

20、1:c=2分析以上图形得：随着矩形脉冲宽度值的增加，时域波形幅值减小，信号频带宽度减小。观察实验结果，思考如下问题：2 1比较矩形脉冲信号和周期矩形脉冲信号的频谱，两者之间有何异同？ 答：异：矩形脉冲信号的频谱是连续的，周期矩形脉冲信号的频谱是离散的。同：有效频带宽度与脉冲宽度成反比当脉冲宽度相同时，二者频谱的包络相同。2 2.根据矩形脉冲宽度 变化时频谱的变化规律，说明信号的有效频带宽度与其时域宽度 之间有什么关系？当脉冲宽度0，脉冲的面积始终等于 1，其频谱有何特点？答：信号的有效频带宽度与其时域宽度成反比。由3）的四个图像可知，当减小时，信号频带宽度增大而最大幅值始终为1,当脉冲宽度0时

21、，信号的频谱无限趋近于一条幅值为1的直线。(3) 已知x(n)是如图所示的周期方波序列。x(n)厶* 止 * I 亠-L I 、丘-N-N1 0 N1nn利用MATLAB绘制出周期方波序列的频谱波形, 变化趋势。N=i nput('N=');N1=i nput( 'N1=');n=1:N;x=o nes(1,N1) zeros(1,N-1-2*N1)on es(1,N1+1);X=fft(x);stem( n,x,'filled');subplot(211)stem( n,x,'filled');改变参数N和N1的大小，观察频谱波

22、形的xlabel('n');title('x( n)N='n um2str(N)'N1='num2str(N1);subplot(212);stem( n,X,'filled');xlabel('k');title('X(k)N='n um2str(N)'N1='num2str(N1)N=9, 2=2时图像如下:x(n) N=9 N1=2dF!1.A1345nX(k) N=9 N1=26789n1s*4w1567890.50641234k由上面六组图可得：N1 一定的时候，随着 N的

23、增大，频谱的谱线逐渐变密，但包络线保持不变。N 一定的时候，随着N1的增大，频谱振幅增大，而且在以横轴为轴的两侧振动更为剧烈。543210-110203040k观察实验结果，思考如下问题：31.以周期方波序列为例，说明周期序列与连续周期信号的频谱有何异同。答：异：连续周期信号频谱随着时间变量t的绝对值的增大，振幅逐渐收敛于零;而周期序列频谱在一个周期内，由中间向两侧振幅逐渐增大。同：二者的频谱都是离散的，具有谐波性。32随着周期方波序列占空比的变化，其频谱如何随之变化？答：随着占空比的增大，谱线逐渐变得稀疏。(4 )已知一矩形脉冲序列。x(n)1,1 n| N10,| n| N1利用MATLAB绘制周期方波序列的频谱波