实验四连续信号的傅立叶变换.docx
《实验四连续信号的傅立叶变换.docx》由会员分享,可在线阅读,更多相关《实验四连续信号的傅立叶变换.docx(12页珍藏版)》请在三一办公上搜索。
1、实验四连续信号的傅立叶变换实验4 非周期信号的傅立叶变换分析 一、 实验目的 (1) 熟悉连续非周期信号频谱特点及其分析方法; (2) 掌握用MATLAB实现傅立叶变换的两种方法; (3) 了解常用傅立叶变换性质的MATLAB实现方法; 二、 实验原理 1、傅里叶变换和其逆变换定义如下: X(jw)=-x(t)e-jwtdt 4.1 x(t)=12pX(jw)e-jwtdw 4.2 连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。按照教材中的说法,任意非周期信号,如果满足狄里克利条件,那么,它可以被看作是由无穷多个不同频率的周期复指数信号ew的线性组合构成的,每个频率所对应的周期复jt
2、指数信号ew称为频率分量,其相对幅度为对应频率的|X(jw)|之jt值,其相位为对应频率的X(jw)的相位。 X(jw)通常为复函数,可以按照复数的极坐标表示方法表示为: X(jw)=| X(jw)|ej X(jw)其中,| X(jw)|称为x(t)的幅度谱,而X(jw)则称为x(t)的相位谱。 给定一个连续时间非周期信号x(t),它的频谱也是连续且非周期的。 2、用MATLAB实现方法 MATLAB进行傅里叶变换有两种方法,一种利用符号运算的方法计算,另一种是数值计算。 2.1 采用数值计算的方法来进行傅里叶变换的计算 严格来说,用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件,
3、即信号是时限信号,也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零,这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。计算机只能处理有限大小和有限数量的数。 采用数值计算算法的理论依据是: X(jw)=-x(t)e-jwtdt=limT0k=-x(kT)e-jkwTT 若信号为时限信号,当时间间隔T取得足够小时,上式可演变为: NX(jw)=Tk=-Nx(kT)e-jkwT=x(t1),x(t2),L,x(t2N+1)e-jwt1,e-jwt2,L,e-jwt2N+1T 上式用MATLAB表示为: X=x*exp(-j*t*w)*T 其中X为信号x(t)的
4、傅里叶变换,w为频率,T为取样间隔。 相应的MATLAB程序: T = 0.01; dw = 0.1; %时间和频率变化的步长 t = -10:T:10; w = -4*pi:dw:4*pi; X(jw)可以按照下面的矩阵运算来进行: X=x*exp(-j*t*w)*T; %傅里叶变换 X1=abs(X); %计算幅度谱 phai=angle(X); %计算相位谱 为了使计算结果能够直观地表现出来,还需要用绘图函数将时间信号x(t),信号的幅度谱|X(jw)|和相位谱 X(jw)分别以图形的方式表现出来,并对图形加以适当的标注。 这里还需要注意,由于在MATLAB运算中,必须对连续信号x(t)
5、进行取样,为了不丢失原信号x(t)的信息,即反变换后能不失真地恢复原来信号x(t),取样间隔T的确定必须满足取样定理的要求,即取样间隔T必须小于奈奎斯特频率。 2.1.1 用MATLAB实现傅里叶逆变换 连续时间傅里叶逆变换可用式4.2进行计算。式4.2重写如下: x(t)=12pX(jw)e-jwtdw 从定义式可看出,其计算方法与傅里叶变换是一样的,因此可以采用同样的矩阵运算的方法来计算,即 x(t)=X(jw)*exp(jw*t)*dw 具体的MATLAB函数如下: t = -5:0.01;5; % 指定信号的时间范围,此范围应根据信号的持续时间确定。 dw = 0.1; w = -4*
6、pi:dw:4*pi; X = input(Type in the expression of X(jw); x = X* exp(jw*t)*dw; 然后用绘图函数就可以绘制出逆变换得到的时域信号波形图。 2.2 符号运算的计算方法 MATLAB提供了两个函数fourier和ifourier, 它们可分别用来计算傅立叶正变换和 反变换,这是计算傅立叶变换的符号计算方法。 傅里叶变换 在Matlab中,傅里变换变换由函数fourier实现。fourier有三种调用格式: F=fourier(f) 求时间函数f(t)的傅里叶变换,返回函数F的自变量默认为w,即F(jw)=Ff(t); F=fou
7、rier(f,v) 求时间函数f(t)的傅里叶变换,返回函数F的自变量为v,即F(jv)=Ff(t); F=fourier(f,u,v) 对自变量为u的函数f(u)求傅里叶变换,返回函数F的自变量为v,即F(jv)=Ff(u)。 傅里叶逆变换 在Matlab中,傅里变换逆变换由函数ifourier实现。与函数fourier相类似,ifourier也有三种调用格式: f=ifourier(F) 求函数F(jw)的傅里叶逆变换,返回函数f的自变量默认为x,即f(x)=FF f=ifourier(F,u) 求函数F(jw)的傅里叶逆变换,返回函数f的自变量为u,即f(u)=FF f=ifourier
8、(F,v,u) 求函数F(jv)的傅里叶逆变换,返回函数f的自变量为u,即f(u)=FF-1-1(jw); -1(jw)。 (jv) 这里要注意的是,在调用上述两个函数之前,先要用syms命令对所用到的变量(如t、u、v、w)等进行说明,也就是要将这些变量说明成符号变量。对于fourier中的函数f或ifourier中的F,也要用syms将f或F说明成为符号表达式。另外,在采用fourier及ifourier得到的返回函数,仍然是符号表达式。若需要对返回函数作图时,只能用ezplot绘图命令,而不能用plot命令。如果返回函数中含有d(w)等项,用ezplot也无法作图。 fourier函数的
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 实验 连续 信号 傅立叶 变换
![提示](https://www.31ppt.com/images/bang_tan.gif)
链接地址:https://www.31ppt.com/p-3436283.html