信号与系统的实验报告汇总.docxVIP

信号与系统实验报告 ——信号的频域分析 班级 学号：1120111487 姓名：蒋志科 实验三 信号的频域分析 ——1120111487 信息工程（实验班） 蒋志科 一、实验目的 ①深入理解信号频谱的概念，掌握信号的频域分析方法。 ②观察典型周期信号和非周期信号的频谱，掌握其频谱特性 二、实验原理和方法 1、连续周期信号的频谱分析 如果周期信号满足狄里赫利条件，那就可以展开为傅里叶级数的形式，即： xt=k=-∞+∞ckejkω0t 式① ck=1T0T0xte-jkω0tdt 式② 式中，T0表示基波周期，ω0=2πT0为基波频率，T0（.）表示任一个基波周期内的积分。 式①和式②定义为周期信号复指数形式的傅里叶级数，系数ck称为 xt称为傅里叶系数。周期信号的傅里叶级数还可以由三角函数的线性组合来表示，即： xt=a0+k=1+∞akcoskω0t+k=1+∞bksinkω0t a0=1T0T0xtdt ak=2T0T0xtcoskω0tdt bk=2T0T0xtsinkω0tdt 同频率的正弦项和余弦项可以合并，得到三角函数形式的傅立叶级数，即： xt=A0+k=1+∞Akcos?（kω0t+θk) 其中A0=a0，Ak=ak2+bk2 θk=-arctanbkak 任何满足狄里赫利条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数叠加。一般说周期信号表示傅里叶级数需要无限多项才能逼近原信号，实际中选项数越多就越逼近原信号。 2、连续非周期信号的频域分析 对于非周期连续时间信号，信号的傅里叶变换和傅里叶逆变换定义为： Xω=-∞+∞x(t)e-jωtdt x(t)=12π-∞+∞Xωejωtdω 两式子把信号的时域特性和频域特性联系起来，确立了非周期信号x(t)和频谱Xω之间关系。采用MATLAB可以方便求取非周期连续时间信号的傅里叶变换。采用几种方法。 1）符合运算法 傅里叶变换和反变换函数分别用fourier和ifourier函数，其格式为X=fourier（x）和x=ifourier(X)默认的时域变量t，频域变量为ω。 2）数值积分法 除了符合运算方法外，也可以利用MATLAB的quad函数，采用数值积分的方法来进行连续信号的频谱分析。也可以利用这个函数计算非周期连续时间信号的频谱。其一般调用格式为：y=quad(fun,a,b) y=quad(fun,a,b,TOL,TRACE,p1,p2,…)其中fun指定被积函数，a、b表示定积分的上下限，TOL表示允许的相对或者绝对积分误差，TRACE表示以被积函数的点绘图形式来跟踪该函数的返回值。TOL,TRACE 为空则默认为缺省值。“p1,p2,…”表示被积函数除时间t之外的的其他额外输入参数。 3）数值近似法 数值近似法近似计算：Xω=-∞+∞x(t)e-jωtdt=lim?→0k=-∞+∞x(k?)e-jkk?? 当x(t)为时限信号且?足够小，就可以演变成Xω=?k=abx(k?)e-jkk?又可以表示成一个行向量和一个列向量的乘积。 3、离散周期时间信号的频域分析 基波周期为N的周期序列x(n)可以用N个成谐波关系的复指数序列的加权和表示，即： xn=k=Nckejk(2πN)n ck=1Nk=Nxne-jk(2πN)n 式k=N表示求和仅需包括一个周期内的N项。我们用周期N与傅里叶系数ck的乘积来表示周期离散时间信号的频谱，即： Xk=N?ckk=Nx(n)e-jk(2πN)n xk可以利用MATLAB提供函数fft用来计算，调用格式为X=fft(x)该函数返回Xk一个周期内的值，其中x表示xn一个周期内的样本值。 4、离散非周期时间信号的频域分析 非周期序列xn可以表示成一组复指数序列的连续和 x(n)=12π2πXejΩejΩndΩ XejΩ=n=-∞+∞x(n)e-jΩn 这称为x(n)的离散时间傅里叶变换，这两式确立了非周期离散时间信号x(n)及其离散时间傅里叶变换XejΩ之间的关系。 三、实验内容 A x(t) -2T -T -τ/2 0 τ/2 T 2T t … … （1）已知x(t)如图所示的周期矩形脉冲信号 ①计算该信号的傅里叶级数； ②利用MATLAB绘出由前N次谐波合成信号波形，观察随着N的变化合成信号波形的变化规律； ③利用MATLAB绘出周期矩形脉冲信号的频谱，观察参数T和τ变化时对频谱波形的影响。 观察实验结果，思考如下问题： Q1-