matlab三角函数的逆,matlab实现ctft及其逆变换的计算

matlab实现ctft及其逆变换的计算

用 MATLAB 实现 CTFT 及其逆变换的计算5.1 用 MATLAB 实现 CTFT 的计算MATLAB 进行傅里叶变换有两种方法，一种利用符号运算的方法计算，另一种是数值计算，本实验要求采用数值计算的方法来进行傅里叶变换的计算。严格来说，用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件，即信号是时限信号(Time limited signal) ，也就是当时间|t|大于某个给定时间时其值衰减为零或接近于零，这个条件与前面提到的为什么不能用无限多个谐波分量来合成周期信号的道理是一样的。计算机只能处理有限大小和有限数量的数。采用数值计算算法的理论依据是：dtetxjXj)((kTjkTex)(lim0若信号为时限信号，当时间间隔 T 取得足够小时，上式可演变为：NkTjkexj)()( Tettx NtjtjtjN],,[],,[ 12211221 上式用 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)分别以图形的方式表现出来，并对图形加以适当的标注。5.2 用 MATLAB 实现傅里叶逆变换 连续时间傅里叶逆变换可用式 9.7 进行计算。式 9.7 重写如下：dejXtxtj)(21)(从定义式可看出，其计算方法与傅里叶变换是一样的，因此可以采用同样的矩阵运算的方法来计算，即x(t)=X(j)*exp(j’*t)*d具体的 MATLAB 函数如下：t = -5:0.01;5; % 指定信号的时间范围，此范围应根据信号的持续时间确定。dw = 0.1; w = -4*pi:d:4*pi;X = (‘Type in the expression of X(jw)’);x = X* exp(jw’*t)*dw;然后用绘图函数就可以绘制出逆变换得到的时域信号波形图。-20 -15 -10 -5 0 5 10 15 2000.20.40.60.811.2三、实验内容和要求实验前，必须首先阅读本实验原理，读懂所给出的全部范例程序。实验开始时，先在计算机上运行这些范例程序，观察所得到的信号的波形图。并结合范例程序应该完成的工作，进一步分析程序中各个语句的作用，从而真正理解这些程序。实验前，一定要针对下面的实验项目做好相应的实验准备工作，包括事先编写好相应的实验程序等事项。1. 编写程序，绘制下面的信号的波形图：)5cos(1)3cos(1)cs() 000 ttttx 10)cos()2sint其中， 0 = 0.5π，要求将一个图形窗口分割成四个子图，分别绘制 cos( 0t)、cos(3 0t)、cos(5 0t) 和 x(t) 的波形图，给图形加 title，网格线和 x 坐标标签，并且程序能够接受从键盘输入的和式中的项数。