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

实验3信号的频域分析（综合型实验）一、实验目的1）深入理解信号频谱的概念，掌握信号的频域分析方法。2）观察典型周期信号和非周期信号的频谱，掌握其频谱特性。二、实验原理与方法1.连续周期信号的频谱分析如果周期信号满足Dirichlet条件，就可展开为傅里叶级数的形式，即0(t)jktkkxce（1）0001(t)ejktkTcxdtT（2）其中0T表示基波周期，002/T为基波频率，0(...)T表示任一个基波周期内的积分。上面两式为周期信号复指数形式的傅里叶级数，系数kc成为x(t)的傅里叶系数。周期信号的傅里叶级数还可由三角函数的线性组合来表示，即00011(t)cossinkkkkxaaktbkt（3）其中000000000122(t),(t)cosk,(t)sinkkkTTTaxdtaxtdtbxtdtTTT（4）（3）式中同频率的正弦、余弦项合并可以得到三角函数形式的傅里叶级数，即001(t)cos(kt)kkkxAA（5）其中2200,,arctankkkkkkbAaAaba（6）任何满足Dirichlet条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数的叠加。周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号，但在实际应用中常采用有限项级数代替，所选级数项越多就越接近原信号。2.连续非周期信号的频谱分析对于非周期连续时间信号，信号的傅里叶变换和傅里叶逆变换定义为()(t)ejtXxdt（7）1(t)()e2jtxXd（8）以上两式把信号的时频特性联系起来，确立了非周期信号(t)x和频谱()X之间的关系。利用MATLAB可以方便地求出非周期连续时间信号的傅里叶变换，几种常见方法如下：1）符号运算法MATLAB的符号数学工具箱提供了直接求解傅里叶变换和反变换的函数，fourier函数和ifourier函数，基本调用格式为(x)Xfourierx(X)ifourier默认的时域变量为t，频域变量为。例：求2|t|(t)ex的傅里叶变换，代码及运行结果如下：symstx=exp(-2*abs(t));X=fourier(x)X=4/(w^2+4)所以傅里叶变换结果为24()4X也可利用int函数直接根据式（7）求傅里叶变换。2）数值积分法除了采用符号运算的方法外，还可以采用MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析。quad函数是一个用来计算数值积分的函数。利用quad函数可以计算非周期连续时间信号的频谱。Quad函数的一般调用格式为：y=quad(fun,a,b)y=quad(fun,a,b,TOL,TRACE,p1,p2,...)其中fun指定被积函数，可以用inline命令来创建，也可通过传递函数句柄的形式来指定，a、b表示定积分的上下限，TOL表示允许的相对或绝对积分误差，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值，如果TOL和TRACE为空矩阵，则使用缺省值，“p1，p2，...”表示被积函数除时间t之外所需的其他额外输入参数。3）数值近似法还可以利用MATLAB的数值计算的方法近似计算连续时间傅里叶变换。傅里叶变换()X可以由（9）式近似计算0()(t)elim(k)ejtjkkXxdtx（9）当x(t)为时限信号，且足够小，则（9）可以演变成()(k)ebjkkaXx（10）（10）中求和部分又可以表示成一个行向量和一个列向量的乘积(1)ee(k)e[x(),x((a1)),...,x(b)]...ejajabjkkajbxa（11）上式可以很方便地利用MATLAB实现。3.离散周期信号的分频域分析基波周期为N的周期序列(n)x可以用N个成谐波关系的复指数序列的加权和表示，即(2/N)n(n)jkkkNxce（12）这里k=N表示求和仅需包括一个周期内的N项，周期序列在一个周期内的求和与起点无关。将周期序列表示成式（12）的形式，称为离散傅里叶级数，而系数kc则称为傅里叶系数。离散傅里叶系数kc可由（13）式确定。(2/N)n1(n)ejkkkNcxN（13）傅里叶系数kc也称为(n)x的频谱系数，而且可以证明kc是以N为周期的离散频率序列。这说明了周期的离散时间函数对应于频域为周期的离散频率。这里，我们用周期N与傅里叶系数kc的乘积来表示周期离散时间信号的频谱，即(2/N)n(k)N(n)ejkkkNXcx（14）X(k)可以利用MATLAB提供的函数fft用来计算，调用格式为X=fft(x)该函数返回X(k)一个周期内的值，其中x表示x(n)一个周期内的样本值。4.离散非周期信号的频域分析非周期序列x(n)可以表示成一组复指数序列的连续和21(n)(e)2jjnxXed（15）其中(e)(n)ejjnnXx（16）式（16）成为x(n)的离散时间傅里叶变换，式（15）和（16）确立了非周期离散时间信号x(n)及其离散时间傅里叶变换(e)jX之间的变换。(e)jX是连续频率的函数，称为频谱函数，且(e)jX是周期的连续频率函数，其周期为2。可见，非周期离散时间函数对应于频域中是一个连续的周期频率函数。对于有限长的离散时间序列，（16）可以表示为12112(e)(n)e[x(n),x(n),...,x(n)]...NNjnjnnjjnNnnjneeXxe（17）上式可以方便地利用MATLAB实现。三、实验内容（1）已知x(t)是如下图所示的周期矩形脉冲信号。1）计算该信号的傅里叶级数；0aAT，2sin(k)kAakT，0kb得傅里叶级数为：122(t)sin(k)cos()kAktxATkTT2）利用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;forn=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)])在命令窗口输入以下语句：x(t)-c/2c/2A-TT......subplot(221)DFTthird_2_1N=5A=1c=0.5T=1subplot(222)DFTthird_2_1N=10A=1c=0.5T=1subplot(223)DFTthird_2_1N=20A=1c=0.5T=1subplot(224)DFTthird_2_1N=40A=1c=0.5T=1图形如下：由以上四个图可知，随着N的增大，合成信号的波形越来越接近原脉冲信号。3）利用MATLAB绘出周期矩形脉冲信号的频谱，观察参数T和变化时对频谱波形的影响。可计算出傅里叶系数为：A,0sin(k),k0kkTcAkT画出该信号频谱MATLAB代码如下：N=input('N=');c=input('c=');A=input('A=');T=input('T=');n1=-N:-1;c1=(A./(n1*pi)).*sin(n1*pi*c/T);c0=c*A/T;n2=1:N;c2=(A./(n2*pi)).*sin(n2*pi*c/T);cn=[c1c0c2];n=-N:N;subplot(211);stem(n,abs(cn),'filled');xlabel('\omega/\omega_0');title('Magnitudeofck');subplot(212);-2-1012-0.500.511.5Time(sec)N=5-2-1012-0.500.511.5Time(sec)N=10-2-1012-0.500.511.5Time(sec)N=20-2-1012-0.500.511.5Time(sec)N=40stem(n,angle(cn),'filled');xlabel('\omega/\omega_0');title('Phaseofck')命令窗口：DFTthird_2_2N=20c=0.5A=1T=4图形：N=20，A=1，改变T和时的波形变化：综合分析以上四个图像可得：频谱的波形与占空比T有关，对于T或取不同的值时，当T不变，频谱波形不变，当T变大，频宽（第一个过零点频率）减小，过零点频率增多，谱线变得稀疏，各频率分量振幅增大。-20-15-10-50510152000.050.10.150.2/0Magnitudeofck-20-15-10-50510152001234/0Phaseofck-20-100102000.020.040.060.08/0Maghitudeofckc=0.5;T=8;c/T=0.0625-20-100102000.050.10.150.2/0Maghitudeofckc=0.5;T=4;c/T=0.125-20-100102000.10.20.30.4/0Maghitudeofckc=1;T=4;c/T=0.25-20-100102000.10.20.30.4/0Maghitudeofckc=0.5;T=2;c/T=0.25观察实验结果，思考如下问题：1—1.什么是吉伯斯现象？产生吉伯斯现象的原因是什么？答：吉伯斯现象：将具有不连续点的周期函数（如矩形脉冲）进行傅里叶级数展开后选取有限项进行合成。在不连续点附近，部分和有起伏，其峰值最大值是不连续点处高度的1.09倍。不连续点处级数收敛于左右极限的平均值，t愈接近不连续点时，为将误差减小至低于某一给定值，N必须取得很大。随着N的增大，部分和的起伏就向不连续点处压缩，但是对有限的N值，起伏的峰值大小不变。产生吉伯斯现象的原因：当一个信号通过某一系统时，如果这个信号是不连续时间函数，则因为一般的物理系统对信号的高频分量都有衰减作用，所以产生了吉伯斯现象。1—2.以周期矩形脉冲为例，说明周期信号的频谱有什么特点？由结果可知：1.周期性矩形脉冲信号的频谱是离散的，其频谱仅存在于000,2,3,...等离散值处，谱线间隔距离为固定值0，各次谐波的频率均为基频0的整数倍。2.随着谐波次数增高，谱线长度逐渐趋于收敛。所以离散型、收敛性以及谐波性是周期信号的共同特点。1—3.周期矩形脉冲信号的有效频带宽度与信号的时域宽度之间有什么关系？答：根据2/BT（B为有效频宽，T为脉冲的时域宽度），有效频带宽度与信号的时域宽度成反比。1—4.随着矩形脉冲信号参数/T的变化，其频谱结构（如频谱包络形状、过零点、谱线间隔等）如何变化？答：随着T的变换，频谱包络形状不变。T增大时，过零点增多，谱线间隔逐渐变大。（2）已知x(t)是如下图所示的矩形脉冲信号。1）求该信号的傅里叶变换；symstcAx=A*[heaviside(t+c/2)-heaviside(t-c/2)];X=fourier(x);collect(X)ans=(2*A*sin((c*w)/2))/w所以傅里叶变换为：2()sin()2AXx(t)-c/2c/2A2）利用MATLAB绘出矩形脉冲信号的频谱，观察矩形脉冲宽度变化时对频谱波形的影响；代码:A=input('A