实验四非周期信号频域分析1实验目的(1)掌握傅里叶变换的分析方法及其物理意义。(2)掌握各种典型的连续时间非周期信号的频谱特征以及傅里叶变换的主要性质。(3)学习掌握利用MATLAB语言编写计算CTFT的仿真程序，并能利用这些程序对一些典型信号进行频谱分析，验证CTFT的若干重要性质。2实验原理及方法2.1连续时间信号傅里叶变换——CTFT傅里叶变换在信号分析中具有非常重要的意义，它主要是用来进行信号的频谱分析的。傅里叶变换和其逆变换定义如下：4-14-2连续时间傅里叶变换主要用来描述连续时间非周期信号的频谱。任意非周期信号，如果满足狄里克利条件，它可以被看作是由无穷多个不同频率(这些频率都是非常的接近)的周期复指数信号ejt的线性组合构成的，每个频率所对应的周期复指数信号ejt称为频率分量，其相对幅度为对应频率的|X(j)|之值，其相位为对应频率的X(j)的相位。X(j)通常为关于的复函数，可以按照复数的极坐标表示方法表示为：X(j)=|X(j)|ejX(j)其中，|X(j)|称为x(t)的幅度谱，X(j)称为x(t)的相位谱。给定一个连续时间非周期信号x(t)，它的频谱是连续且非周期的。对于连续时间周期信号，也可以用傅里叶变换来表示其频谱，其特点是，连续时间周期信号的傅里叶变换是由冲激序列构成的，是离散的——这是连续时间周期信号的傅里叶变换的基本特征。2.2用MATLAB实现CTFT及其逆变换2.2.1用MATLAB实现CTFT的计算MATLAB进行傅里叶变换有两种方法，一种利用符号运算的方法计算，另一种是数值计算，本实验采用数值计算的方法。严格来说，用数值计算的方法计算连续时间信号的傅里叶变换需要有个限定条件，即信号是时限信号，也就是当时间|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)分别以图形的方式表现出来，并对图形加以适当的标注。例4-1：绘制非周期矩形脉冲信号时域波形及频谱函数。%Program4_1T=0.01;dw=0.1;t=-10:T:10;w=-6*pi:dw:6*pi;x=u(t+1/2)-u(t-1/2);X=x*exp(-j*t'*)*T;subplot(211),plot(t,x),title('非周期矩形信号')axis([-2,2,-0.2,1.5])subplot(212),plot(w,real(X))title('频谱函数'),axis([-20,20,-0.5,1.1])例4-2：绘制单边指数信号时域波形、信号的幅度频谱和相位频谱图。%Program4_2T=0.01;dw=0.1;%时间和频率变化的步长t=-10:T:10;w=-4*pi:dw:4*pi;x=exp(-t).*u(t);X=x*exp(-j*t'*)*T;%傅里叶变换X1=abs(X);%计算幅度谱phai=angle(X);%计算相位谱subplot(311),plot(t,x)title('单边指数信号信号')axis([-2,6,-0.2,1.5])subplot(312),plot(w,X1),title('幅度谱')axis([-15,15,-0.5,1.1])subplot(313),plot(w,phai),title('相位谱')axis([-15,15,-pi/2,pi/2])2.3.2用MATLAB实现傅里叶逆变换连续时间傅里叶逆变换用式4-2进行计算。式4-2重写如下：从定义式可看出，其计算方法与傅里叶变换是一样的，因此可以采用同样的矩阵运算的方法来计算，即x(t)=X(j)*exp(j*'*t)*d/(2*pi)例4-3：绘制非周期矩形脉冲信号和逆变换得到的时域信号波形图，并分析误差原因。MATLAB程序如下：