实验六傅里叶变换及其反变换.doc

实验六傅里叶变换及其反变换 6.1实验目的 ?学会运用MATLAB求连续时间信号的傅里叶变换； .学会运用MATLAB求连续时间信号的傅里叶反变换； 3.学会运用MATLAB求连续吋间信号的频谱图。 6. 2实验原理及实例分析 1.连续时间信号傅里叶变换?…CTFT 傅里叶变换在信号分析屮具有非常重耍的意义，它主要是用來进行信号的频谱分析的。 傅里叶变换和其逆变换定义如下： 00 X(je)=⑴茁购d/ ]00 6.2x(t) = — \x(jco)严 dco 6.2 —00 连续吋间傅里叶变换主要用来描述连续吋间非周期信号的频谱。按照教材小的说法，任 意非周期信号，如果满足狄里克利条件，那么，它可以被看作是由无穷多个不同频率(这些 频率都是非常的接近)的周期复指数信号』°’的线性组合构成的，每个频率所对应的周期复 指数信号』称为频率分赧(frequency component), K相对幅度为对应频率的|X(j?)|之值， 其相位为对应频率的X(jo)的相位。 X(j(o)通常为关于的复函数，可以按照复数的极他标表示方法表示为： X(jo)=| X(j(o)|ejzx(iw) 其中,|X0(o)|称为x(t)的幅度谱,而ZX0CO)则称为x⑴的相位谱。 给定一个连续时间非周期信号x⑴，它的频谱也是连续且非周期的。对于连续时间周期 信号，也可以用傅里变换来表示其频谱，其特点是，连续时间周期信号的傅里叶变换时冇冲 激序列构成的，是离散的——这是连续吋间周期信号的傅里叶变换的基本特征。 2 .用MATLAB实现CTFT的计算 MATLAB进行傅里叶变换有两种方法，-?种利用符号运算的方法计算，另一种是数值计算。 1) MATLAB符号运算求解法 MATLAB符号数学工具箱提供了直接求解傅里叶变换与傅里叶反变换的函数 fourier( ifourier( )o常用的是：F=fourier(f) 默认返回值是关于?的函数。 f=fourier(F,t) 返回值是关于t的函数 例：利用MATLAB求单边指数信号f(t) = 的傅里叶变换，画岀f⑴及具幅度谱和相位 谱图。 %定义符号变量%f (t) =exp (- %定义符号变量 %f (t) =exp (-2*t) *u (t) %求傅里叶变换 %绘制f(t)的时域波形 %绘制幅度谱 %计算F(w)的虚部 f = exp(-2*t)*sym(，Heaviside(t)?); Fw 二 fourier(f); subplot(311); ezplot (f); axis([-l 2. 5 0 1. 1]); subplot(312); ezplot(abs(Fw)); im = imag(Fw); re = real(Fw); phase = atan(im/re); subplot(313); ezplot(phase); %计算F(w)的实部%计算相位谱%绘制相位谱1/abs(2+i w)w-i atanh(1/2 (1/(2+i w)-1/(2-i conj(w)))/(1/2/(2+i w)+1/2/(2-i conj(w)))) %计算F(w)的实部 %计算相位谱 %绘制相位谱 1/abs(2+i w) w -i atanh(1/2 (1/(2+i w)-1/(2-i conj(w)))/(1/2/(2+i w)+1/2/(2-i conj(w)))) w 2) MATLAB数值计算求解法 符号运算求解法的局限性在于，如果返回函数中冇诸如6 (t)等项，则川ezplot()函数无法 作图。对某些信号求变换时，其返回函数可能包含一些不能直接用符号表达的式子，因而也 不返回函数作图。故有必要给出连续信号傅里叶变换的数值计算法。 采用数值计算算法的理论依据是： 00X(jco)= 00 X(jco)= —co 00 =lim工兀伙门厂畑7Y 若信号为吋限信号，当吋间间隔T収得足够小吋，上式可演变为: N X(jco) = T^x(kT)e-j^T k=-N =［心)，曲2 )，…，或2郴)］° ［严 0叽，…,厂3 V 上式用MATLAB表示为： X=x*cxp(-j *t, *w)*T 其中X为信号x(t)的傅里叶变换，w为频率Q, T为时间步长。 相应的MATLAB程序： T = 0.01;dw = 0.1; %时间和频率变化的步长 t = -10:T:10; w =?4*pi:dw:4*pi; %X(j3)可以按照下而的矩阵运算来进行： X=x *cxp(-j*f*(D)*T %傅里叶变换 Xl=abs(X); %计算幅度谱 phai=angle(X); %计算相位谱 为了使计算结果能够肓观地表现出来，还需要用绘图函数将时间信号x(t),信号的幅度谱 |X(jco)|和相位