实验四连续信号的傅立叶变换.doc

实验 4 非周期信号的傅立叶变换分析一、实验目的(1) 熟悉连续非周期信号频谱特点及其分析方法; (2) 掌握用 MATLAB 实现傅立叶变换的两种方法; (3) 了解常用傅立叶变换性质的 MATLAB 实现方法; 二、实验原理 1 、傅里叶变换和其逆变换定义如下: ??????dtetxjX tj??)()( ?????????dejXtx tj)(2 1)( 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。按照教材中的说法,任意非周期信号, 如果满足狄里克利条件, 那么, 它可以被看作是由无穷多个不同频率( 这些频率都是非常的接近) 的周期复指数信号 e j? t 的线性组合构成的, 每个频率所对应的周期复指数信号 e j? t 称为频率分量( ponent ) ,其相对幅度为对应频率的|X(j ?)|之值,其相位为对应频率的 X(j ?) 的相位。 X(j ?) 通常为复函数,可以按照复数的极坐标表示方法表示为: X(j ?)=| X(j ?)|e j? X(j ?) 其中, | X(j ?)| 称为 x(t) 的幅度谱,而? X(j ?) 则称为 x(t) 的相位谱。给定一个连续时间非周期信号 x(t) ,它的频谱也是连续且非周期的。 2 、用 MATLAB 实现方法 MATLAB 进行傅里叶变换有两种方法,一种利用符号运算的方法计算,另一种是数值计算。 采用数值计算的方法来进行傅里叶变换的计算严格来说, 用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件, 即信号是时限信号( Time limited signal ) ,也就是当时间|t| 大于某个给定时间时其值衰减为零或接近于零, 这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。计算机只能处理有限大小和有限数量的数。采用数值计算算法的理论依据是: ( ) ( ) j t X j x t e dt ???????????????? k T jk TTe kT x ?)( lim 0 若信号为时限信号,当时间间隔 T 取得足够小时,上式可演变为: ????? NNk T jke kT xTjX ??)()(Teeetxtxtx tjtjtjN],,,[ )](, ),( ),([ 1221 ???????????上式用 MATLAB 表示为: X=x*exp(-j*t ’*w)*T 其中 X 为信号 x(t) 的傅里叶变换, w 为频率, T 为取样间隔。相应的 MATLAB 程序: T= ; dw= ; % 时间和频率变化的步长 t= -10: T :10; w= -4*pi: d w:4*pi; X(j ?) 可以按照下面的矩阵运算来进行: X=x*exp(-j*t ’*?)*T;% 傅里叶变换 X1=abs(X); % 计算幅度谱 phai=angle(X); % 计算相位谱为了使计算结果能够直观地表现出来,还需要用绘图函数将时间信号 x(t) ,信号的幅度谱|X(j ?)| 和相位谱? X(j ?) 分别以图形的方式表现出来,并对图形加以适当的标注。这里还需要注意,由于在 MATLAB 运算中,必须对连续信号( ) x t 进行取样,为了不丢失原信号( ) x t 的信息, 即反变