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

.. .. .. 实验 3 信号的频域分析 （综合型实验 ） 一、实验目的 1）深入理解信号频谱的概念 ，掌握信号的频域分析方法 。 2）观察典型周期信号和非周期信号的频谱 ，掌握其频谱特性 。 二、实验原理与方法 1.连续周期信号的频谱分析 如果周期信号满足 Dirichlet 条件 ，就可展开为傅里叶级数的形式 ，即 x(t) ck ejk 0 t （ 1 ） ck 1 x(t)e jk 0t dt （ 2） k T0 T0 其中 T0 表示基波周期 ， 0 2 / T0 为基波频率 ， (...) 表示任一个基波周期内的积 T0 分。 上面两式为周期信号复指数形式的傅里叶级数 ，系数 ck 成为 x(t)的傅里叶系数 。周期信 号的傅里叶级数还可由三角函数的线性组合来表示 ，即 x(t) a0 ak cosk 0t bk sin k 0t （ 3 ） k 1 k 1 其中 a0 1 2 2 x(t)sink 0tdt （ 4） x(t)dt, ak x(t)cosk 0tdt, bk T0 T T0 T T0 T 0 0 0 （ 3 ）式中同频率的正弦 、余弦项合并可以得到三角函数形式的傅里叶级数 ，即 x(t) A0 Ak cos(k 0 t k ) （5） k 1 其中 A0 a0 , Akak2 bk2 , k arctanbk （6） ak 任何满足 Dirichlet 条件的周期信号都可以表示成一组谐波关系的复指数函数或三角函数 . 专业 .专注 . .. .. .. 的叠加 。 周期信号表示为傅里叶级数时需要无限多项才能完全逼近原信号 ，但在实际应用 中常采用有限项级数代替 ，所选级数项越多就越接近原信号 。 2.连续非周期信号的频谱分析 对于非周期连续时间信号 ，信号的傅里叶变换和傅里叶逆变换定义为 X ( ) x(t)e j t dt （ 7） x(t) 1 X ( )e j t d （ 8） 2 以上两式把信号的时频特性联系起来 ，确立了非周期信号 x(t) 和频谱 X ( ) 之间的关 系。 利用 MATLAB 可以方便地求出非周期连续时间信号的傅里叶变换 ，几种常见方法如下 ： 1 ） 符号运算法 MATLAB 的符号数学工具箱提供了直接求解傅里叶变换和反变换的函数 ， fourier 函数和 ifourier 函数 ，基本调用格式为 X fourier (x) x ifourier (X) 默认的时域变量为 t，频域变量为。 例：求 x(t) e 2|t| 的傅里叶变换 ，代码及运行结果如下 ： syms t x=exp(-2*abs(t)); X=fourier(x) X = 4/(w^2 + 4) . 专业 .专注 . .. .. .. 所以傅里叶变换结果为 X ( )  4 2 4 也可利用 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）式近似计算 X ( ) x(t)e j t dt lim x(k )e j k （ 9） k 当 x(t)为时限信号 ，且 足够小 ，则（ 9 ）可以演变成 . 专业 .专注 . .. .. .. b X ( ) x(k )e jk （ 10） a 10）中求和部分又可以表示成一个行向量和一个列向量的乘积 e b e x(k )e jk [x( a ), x((a 1) ),..., x(b )] k a ... e  ja j ( a 1) （ 11） jb 上式可以很方便地利用 MATLAB 实现 。 3.离散周期信号的分频域