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

1、 标准实验 3 信号的频域分析（综合型实验）一、实验目的1）深入理解信号频谱的概念，掌握信号的频域分析方法。2）观察典型周期信号和非周期信号的频谱，掌握其频谱特性。二、实验原理与方法1.连续周期信号的频谱分析如果周期信号满足 dirichle条t 件，就可展开为傅里叶级数的形式，即1ckx(t)e dtjk t （2）0x(t)c ekjk t （1）0t0 tk0t其中 表示基波周期，2 /t为基波频率， (.)表示任一个基波周期内的积分。000t0c上面两式为周期信号复指数形式的傅里叶级数，系数 成为 x(t的) 傅里叶系数。周期信k号的傅里叶级数还可由三角函数的线性组合来表示，即x(t)

2、 aa cosk tb sink t（3）0k0k0k 1k 1122a0x( td)t,ax(t)cosk tdt,bx(t)sink tdt（4）其中t0 tt0 tt0 tk0k0000（3）式中同频率的正弦、余弦项合并可以得到三角函数形式的傅里叶级数，即x(t) aa cos(k t )（5）0k0kk 1ba a ,aa b ,arctank其中2k2k（6）a00kkk任何满足 dirichle条t 件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号，但在实际应用中常采用有限项级数代替，所选级数项越多就越接近原

3、信号。2.连续非周期信号的频谱分析对于非周期连续时间信号，信号的傅里叶变换和傅里叶逆变换定义为1x ( )x(t)e dtj tx(t)x ( )e dj t（7）（8）2以上两式把信号的时频特性联系起来，确立了非周期信号和频谱x(t) x ( )之间的关系。利用 matlab 可以方便地求出非周期连续时间信号的傅里叶变换，几种常见方法如下：1） 符号运算法文案 标准matlab 的符号数学工具箱提供了直接求解傅里叶变换和反变换的函数， fourier函数和 ifourier 函数，基本调用格式为x = fourier(x)x = ifourier(x)默认的时域变量为 t，频域变量为w 。=

4、 e例：求 x(t)的傅里叶变换，代码及运行结果如下：-2|t| syms t x=exp(-2*abs(t); x=fourier(x)x =4/(w2 + 4)4w=所以傅里叶变换结果为 x ( )w4 + 2也可利用 int 函数直接根据式（7）求傅里叶变换。2） 数值积分法除了采用符号运算的方法外，还可以采用 matlab 的 quad 函数，采用数值积分的方法来进行连续信号的频谱分析。quad 函数是一个用来计算数值积分的函数。利用 quad 函数可以计算非周期连续时间信号的频谱。quad 函数的一般调用格式为：y=quad(fun,a,b)y=quad(fun,a,b,tol,tr

5、ace,p1,p2,.)其中 fun 指定被积函数，可以用 inline 命令来创建，也可通过传递函数句柄的形式来指定，a、b 表示定积分的上下限，tol 表示允许的相对或绝对积分误差，trace表示以被积函数的点绘图形式来跟踪该函数的返回值，如果tol 和 trace 为空矩阵，则使用缺省值，“p1，p2，.”表示被积函数除时间 t 之外所需的其他额外输入参数。3） 数值近似法还可以利用 matlab 的数值计算的方法近似计算连续时间傅里叶变换。傅里叶变换 x (w ) 可以由（9）式近似计算+x (w ) = x(t)e w dt = limx(k d)e w d- j t- j kd（9

6、）d0k=-当 x(t)为时限信号，且d 足够小，则（9）可以演变成bx (w ) = d x(k )edw- jk d（10）k =a（10）中求和部分又可以表示成一个行向量和一个列向量的乘积文案 标准e- jadwe- j(a+1)dw bx(k )ed=x(a ), x(a 1) ),., x(b )+ d d d- jkwd（11）.k=ae- jbdw上式可以很方便地利用 matlab实现。3.离散周期信号的分频域分析基波周期为 n 的周期序列 x(n) 可以用 n 个成谐波关系的复指数序列的加权和表示，即x(n) =c ekjk (2p（12）/ n)nk =这里 k=表示求和仅需

7、包括一个周期内的 n 项，周期序列在一个周期内的求和与起点无关。将周期序列表示成式（12）的形式，称为离散傅里叶级数，而系数 c 则称为傅里叶k系数。离散傅里叶系数c 可由（13）式确定。k1c =kx(n)ep（13）- jk(2 / n)nnk=傅里叶系数c 也称为 x(n) 的频谱系数，而且可以证明c 是以 n 为周期的离散频率序列。kk这说明了周期的离散时间函数对应于频域为周期的离散频率。这里，我们用周期 n 与傅里叶系数c 的乘积来表示周期离散时间信号的频谱，即kx (k) = nc =x(n) e- jk (2p（14）/ n)nkk =x(k)可以利用 matlab提供的函数 f

8、ft 用来计算，调用格式为x=fft(x)该函数返回 x(k)一个周期内的值，其中 x 表示 x(n)一个周期内的样本值。4. 离散非周期信号的频域分析非周期序列 x(n)可以表示成一组复指数序列的连续和1+x(n) = x (e )e dw其中 x (e jw ) x(n) e- wn= （16）jwjwn（15）j2p2pn=-式（16）成为 x(n)的离散时间傅里叶变换，式（15）和（16）确立了非周期离散时间信号 x(n)及其离散时间傅里叶变换 x (e ) 之间的变换。 x (e ) 是连续频率w 的函数，称为jwjw频谱函数，且 x (e ) 是周期的连续频率函数，其周期为 2 。

9、可见，非周期离散时间函数pjw对应于频域中是一个连续的周期频率函数。对于有限长的离散时间序列，（16）可以表示为文案 标准e- jn w1- jn wnx (e ) =x(n)e=x(n ), x(n ),., x(n )- jwn2ne（17）jw12n.n=n1e- jn wn上式可以方便地利用 matlab实现。三、实验内容（1）已知 x(t)是如下图所示的周期矩形脉冲信号。x(t)a.-t-c/2c/2t1）计算该信号的傅里叶级数；t2akpta = a， a =sin(kp=) ，b 0ttk0k得傅里叶级数为：tx(t) = a +ttp2 kt2akp+p) cos(sin(k)

10、ttk =12）利用 matlab绘出由前 n 次谐波合成的信号波形，观察随着n 的变化合成信号波形的变化规律；用 matlab编写代码如下：n=input(n=);a=input(a=);c=input(c=);t=input(t=);t=-1.5:0.001:1.5;x=a*c*ones(size(t)/t;for n=1:nx=x+(2*a/(n*pi)*sin(n*pi*c/t)*cos(2*pi*n*t/t);endplot(t,x);xlabel(time(sec)title(n= num2str(n)在命令窗口输入以下语句：文案 标准 subplot(221) subplot(2

11、23) dftthird_2_1 dftthird_2_1n=5n=20a=1a=1c=0.5c=0.5t=1t=1 subplot(222) subplot(224) dftthird_2_1 dftthird_2_1n=10n=40a=1a=1c=0.5t=1c=0.5t=1图形如下：n=5n=101.51.5110.500.50-0.5-0.5-2-1012-2-1012time(sec)n=20time(sec)n=401.51.5110.500.50-0.5-0.5-2-1012-2-1012time(sec)time(sec)由以上四个图可知，随着 n 的增大，合成信号的波形越来越

12、接近原脉冲信号。t3） 利用 matlab绘出周期矩形脉冲信号的频谱，观察参数 和 变化时对频谱波形的影响。t可计算出傅里叶系数为：ta,k = 0tc = tka psin(kp),k 0kt画出该信号频谱 matlab代码如下：n=input(n=);c2=(a./(n2*pi).*sin(n2*pi*c/t);cn=c1 c0 c2;c=input(c=);a=input(a=);n=-n:n;t=input(t=);subplot(211);n1=-n:-1;stem(n,abs(cn),filled);xlabel(omega/omega_0);title(magnitude of

13、ck);subplot(212);c1=(a./(n1*pi).*sin(n1*pi*c/t);c0=c*a/t;n2=1:n;文案 标准stem(n,angle(cn),filled);xlabel(omega/omega_0);title(phase of ck)命令窗口： dftthird_2_2n=20c=0.5a=1t=4图形：magnitude of ck0.20.150.10.050-20-15-10-505101520w/w0phase of ck43210-20-15-10-505101520w/w0tn=20，a=1，改变 和 时的波形变化：tc=0.5;t=8;c/t=0

14、.0625c=0.5;t=4;c/t=0.1250.080.060.040.0200.20.150.1kkccffooedueduttiihhgg0.05aamm0-20-1001020-20-1001020w/w0w/w0c=1;t=4;c/t=0.25c=0.5;t=2;c/t=0.250.40.30.20.100.40.30.20.10kkccffooedueduttiihhggaamm-20-1001020-20-1001020w/w0w/w0t综合分析以上四个图像可得：频谱的波形与占空比 有关，对于 或 取不同的值时，当tttttt不变，频谱波形不变，当 变大，频宽（第一个过零点频率

15、）减小，过零点频率增多，谱t线变得稀疏，各频率分量振幅增大。文案 标准观察实验结果，思考如下问题：11. 什么是吉伯斯现象？产生吉伯斯现象的原因是什么？答：吉伯斯现象：将具有不连续点的周期函数（如矩形脉冲）进行傅里叶级数展开后选取有限项进行合成。在不连续点附近，部分和有起伏，其峰值最大值是不连续点处高度的 1.09 倍。不连续点处级数收敛于左右极限的平均值，t 愈接近不连续点时，为将误差减小至低于某一给定值，n 必须取得很大。随着 n 的增大，部分和的起伏就向不连续点处压缩，但是对有限的 n 值，起伏的峰值大小不变。产生吉伯斯现象的原因：当一个信号通过某一系统时，如果这个信号是不连续时间函数，

16、则因为一般的物理系统对信号的高频分量都有衰减作用，所以产生了吉伯斯现象。12. 以周期矩形脉冲为例，说明周期信号的频谱有什么特点？由结果可知：w w w w= ,2 ,3 ,.1. 周期性矩形脉冲信号的频谱是离散的，其频谱仅存在于等离000ww散值处，谱线间隔距离为固定值 ，各次谐波的频率均为基频 的整数倍。002. 随着谐波次数增高，谱线长度逐渐趋于收敛。所以离散型、收敛性以及谐波性是周期信号的共同特点。13. 周期矩形脉冲信号的有效频带宽度与信号的时域宽度之间有什么关系？答：= 2 / tb根据 bwp（为有效频宽， 为脉冲的时域宽度），有效频带宽度与信号的tw时域宽度成反比。14. 随着

17、矩形脉冲信号参数t /t的变化，其频谱结构（如频谱包络形状、过零点、谱线间隔等）如何变化？答：tt随着 的变换，频谱包络形状不变。 增大时，过零点增多，谱线间隔逐渐变大。tt（2）已知 x(t)是如下图所示的矩形脉冲信号。-c/2c/21）求该信号的傅里叶变换； syms t c a x=a*heaviside(t+c/2)-heaviside(t-c/2); x=fourier(x);ans = collect(x)(2*a*sin(c*w)/2)/w所以傅里叶变换为：wtsin( )22awwx ( ) =文案 标准t2）利用 matlab绘出矩形脉冲信号的频谱，观察矩形脉冲宽度 变化时对

18、频谱波形的影响；代码:a=input(a=);c=input(c=);syms t wx=int(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(magnitude);title(|x(omega)|)图形如下：|x(w)|10.8edut0.60.40.20ingam-15-10-5051015wt令 a=1，改变 的大小：c=0.25c=0.50.40.30.20.10.2edutedutiinng0.150.1gaa

19、mm-100w10-100w10c=1c=20.80.60.40.21.51eedutdutiinnggaamm0.5-10010-10010wwt增大时，频谱振幅增大，有效频宽减小。3）.让矩形脉冲的面积始终等于 1，改变矩形脉冲宽度，观察矩形脉冲信号时域波形和频谱随矩形脉冲宽度的变化趋势。syms t wtitle(x(t);c=input(c=);subplot(212);x=(1/c).*(heaviside(t+c/2)-heaviside(t-c/2);ezplot(abs(x),-6*pi,6*pi);grid on;x=int(x*exp(-j*w*t),t,-c/2,c/2)

20、;subplot(211);xlabel(omega);ylabel(magnitude);title(|x(w)|)axis tightezplot(x,-c,c);xlabel(t);文案 标准时域波形及频谱在不同脉冲宽度下的图像：c=0.25c=0.5x(t)x(t)4321021.510.50-0.25 -0.2 -0.15 -0.1 -0.050t0.050.10.150.20.25-0.5-0.4-0.3-0.2-0.10t0.10.20.30.40.5|x(w)|x(w)|0.80.60.40.20.80.60.4edutedutingamigam-15-10-5051015-1

21、5-10-50w51015wc=1：c=2x(t)x(t)10.40.200.50-1-0.8-0.6-0.4-0.20t0.20.40.60.81-2-1.5-1-0.50t0.511.52|x(w)|x(w)|0.80.60.40.20.80.60.40.2edueduttiingamngam-15-10-5051015-15-10-5051015ww分析以上图形得：随着矩形脉冲宽度值的增加，时域波形幅值减小，信号频带宽度减小。观察实验结果，思考如下问题：21.比较矩形脉冲信号和周期矩形脉冲信号的频谱，两者之间有何异同？答：异：矩形脉冲信号的频谱是连续的，周期矩形脉冲信号的频谱是离散的。同

22、：有效频带宽度与脉冲宽度成反比当脉冲宽度相同时，二者频谱的包络相同。t22.根据矩形脉冲宽度 变化时频谱的变化规律，说明信号的有效频带宽度与其时域宽度t之间有什么关系？当脉冲宽度答： 0，脉冲的面积始终等于 1，其频谱有何特点？信号的有效频带宽度与其时域宽度成反比。t由 3）的四个图像可知，当 减小时，信号频带宽度增大而最大幅值始终为1，当脉冲宽度t 0时，信号的频谱无限趋近于一条幅值为 1 的直线。文案 标准（3）已知 x(n)是如图所示的周期方波序列。x(n)1-n-n10n1nn利用 matlab绘制出周期方波序列的频谱波形，改变参数n 和 n 的大小，观察频谱波形的1变化趋势。n=in

23、put(n=);xlabel(n);n1=input(n1=);n=1:n;title(x(n) n= num2str(n) n1=num2str(n1);x=ones(1,n1) zeros(1,n-1-2*n1)ones(1,n1+1);subplot(212);stem(n,x,filled);xlabel(k);x=fft(x);stem(n,x,filled);subplot(211)title(x(k) n= num2str(n) n1=num2str(n1)stem(n,x,filled);n=9，n1=2 时图像如下：x(n) n=9 n1=210.50123456789nx(

24、k) n=9 n1=26420-212345k6789改变nn 和 的大小，频谱波形变化如下：15498763524312100-1-105k1005k10文案 标准x(k) n=40 n1=2x(k) n=40 n1=4x(k) n=20 n1=2x(k) n=20 n1=454108541083636242412120000-1-2-1-20510k15200510k152001020k304001020k3040由上面六组图可得：n 一定的时候，随着n 的增大，频谱的谱线逐渐变密，但包络线保持1不变。n 一定的时候，随着 n 的增大，频谱振幅增大，而且在以横轴为轴的两侧振动更为剧烈。1观察实验结果，思考如下问题：31.以周期方波序列为例，说明周期序列与连续周期信号的频谱有何异同。答：异：连续周期信号频谱随着时间变量 t 的绝对值的增大，振幅逐渐收敛于零；而周期序列频谱在一个周期内，由中间向两侧振幅逐渐增大。同：二者的频谱都是离散的，具有谐波性。32.随着周期方波序列占空比的变化，其频谱如何随之变化？答：随着占空比的增大，