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

实验4 非周期信号的傅立叶变换分析 实验目的 熟悉连续非周期信号频谱特点及其分析方法； 掌握用MATLAB实现傅立叶变换的两种方法； 了解常用傅立叶变换性质的MATLAB实现方法； 实验原理 1、傅里叶变换和其逆变换定义如下： 4.1 4.2 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。按照教材中的说法，任意非周期信号，如果满足狄里克利条件，那么，它可以被看作是由无穷多个不同频率（这些频率都是非常的接近）的周期复指数信号ej(t的线性组合构成的，每个频率所对应的周期复指数信号ej(t称为频率分量（frequency component），其相对幅度为对应频率的|X(j()|之值，其相位为对应频率的X(j()的相位。 X(j()通常为复函数，可以按照复数的极坐标表示方法表示为： X(j()=| X(j()|ej( X(j() 其中，| X(j()|称为x(t)的幅度谱，而(X(j()则称为x(t)的相位谱。 给定一个连续时间非周期信号x(t)，它的频谱也是连续且非周期的。 2、用MATLAB实现方法 MATLAB进行傅里叶变换有两种方法，一种利用符号运算的方法计算，另一种是数值计算。 采用数值计算的方法来进行傅里叶变换的计算 严格来说，用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件，即信号是时限信号（Time limited signal），也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零，这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。计算机只能处理有限大小和有限数量的数。 采用数值计算算法的理论依据是： 若信号为时限信号，当时间间隔T取得足够小时，上式可演变为： 上式用MATLAB表示为： X=x*exp(-j*t’*w)*T 其中X为信号x(t)的傅里叶变换，w为频率，T为取样间隔。 相应的MATLAB程序： T = 0.01; dw = 0.1; %时间和频率变化的步长 t = -10:T:10; w = -4*pi:dw:4*pi; X(j()可以按照下面的矩阵运算来进行： X=x*exp(-j*t’*()*T; %傅里叶变换 X1=abs(X); %计算幅度谱 phai=angle(X); %计算相位谱 为了使计算结果能够直观地表现出来，还需要用绘图函数将时间信号x(t)，信号的幅度谱|X(j()|和相位谱( X(j()分别以图形的方式表现出来，并对图形加以适当的标注。 这里还需要注意，由于在MATLAB运算中，必须对连续信号进行取样，为了不丢失原信号的信息，即反变换后能不失真地恢复原来信号，取样间隔T的确定必须满足取样定理的要求，即取样间隔T必须小于奈奎斯特频率。 2.1.1 用MATLAB实现傅里叶逆变换 连续时间傅里叶逆变换可用式4.2进行计算。式4.2重写如下： 从定义式可看出，其计算方法与傅里叶变换是一样的，因此可以采用同样的矩阵运算的方法来计算，即 x(t)=X(j()*exp(j(’*t)*d( 具体的MATLAB函数如下： t = -5:0.01;5; % 指定信号的时间范围，此范围应根据信号的持续时间确定。 dw = 0.1; w = -4*pi:d(:4*pi; X = input(‘Type in the expression of X(jw)’); x = X* exp(jw’*t)*dw; 然后用绘图函数就可以绘制出逆变换得到的时域信号波形图。 符号运算的计算方法 MATLAB提供了两个函数fourier()和ifourier()， 它们可分别用来计算傅立叶正变换和 反变换，这是计算傅立叶变换的符号计算方法。 （1）傅里叶变换 在Matlab中，傅里变换变换由函数fourier()实现。fourier()有三种调用格式： ① F=fourier(f) 求时间函数f(t)的傅里叶变换fourier(f,v) 求时间函数f(t)的傅里叶变换fourier(f,u,v) 对自变量为u的函数f(u)求傅里叶变换Matlab中，傅里变换逆变换由函数ifourier()实现。与函数fourier()相类似，ifourier()也有三种调用格式： ① f=ifourier(F) 求函数F(j?)的傅里叶变换fourier(F,u) 求函数F(j?)的傅里叶变换fourier(F,v,u) 求函数F