用matlab验证傅里叶变换的基本性质

任务2

1. 线性叠加性

ax1t+bx2t↔aX1f+bX2f

构造信号2sin(2π50t)+4sin(2π100t) ↔2X1f+4X2f 满足线性叠加性质。

1. 时间尺度改变性

xat↔1aX(fa)

该式表明若信号xt 在时间轴上被压缩成原信号的1a ，则其频谱函数在频率轴上将展宽a 倍，而其幅值相应地减至原信号幅值地1a 。时间尺度改变特性表面，信号的持续时间与信号占有的频带宽成反比。

构造矩形脉冲信号f(t) 、f(2t) 、f(0.5t) ，将f(t) 在时域上压缩成原信号的12 （即信号f(2t) ），绘图发现f(2t) 的频带宽是f(t) 的2倍，而其幅值相应地减至原信号幅值地12 ，符合时间尺度改变特性。

1. 时移特性

xt±t0↔Xfe±j2πft0

由公式可得信号在时域上平移，则在频域上只是附加了一个线性相位，例如正弦信号 sin（t） 和 sin（t+45°）， 两个信号幅度、频率是相同的，区别仅仅是在初始相位上，画出两信号的波形图，我们知道后一信号相对于前一信号仅仅是发生了平移，这就比较好理解了。

用Matlab分别绘出 x1= 12e -2tu(t) ；x2= 12e -2t+t1u(t+t1) ; x3= 12e -2t-t1u(t-t1) 的幅度频谱和相位频谱，绘图结果如图所示，上图代表信号的时域波形，可以看出信号向左平移t1 个单位、向右平移 t1 个单位；中图是信号的幅值谱，可以看到平移后信号的幅值谱未变，三个信号的幅值谱重叠；下图是信号的相位谱，可以看出三个信号的相位谱是不同的。

1. 频移特性

xte±j2πf0t↔X(f∓f0)

结论:将信号xt 乘以因子ej2πf0t ,对应于将频谱函数沿频率轴右移f0 ;将信号xt 乘以因子e-j2πf0t ，对应于将频谱函数沿频率轴左移f0 。

构造信号xt=2sin(2π50t)，乘以因子ej2π10t ，绘图发现频域图右移了10Hz,如图所示：

1. 时域卷积特性

x1t*x2t↔X1(f)X2(f)

对于信号的时域卷积有: (1)离散时间信号的卷积和:它的定义为离散时间信号x1(n) 和x2(n) 的卷积和为:

xn=x1n*x2n=m=-∞∞x1m*x2n-m

设序列x1n 在区间n1~n2 ,非零，序列x2n 在m1~m2 非零，那么就有xn=x1n*x2n 的非零区间就为(n1+m1) ~(n2+m2) 并且区间长度为(n1+m1) ~(n2+m2)+1 ，则只需计算序列x(n)的非零区间就可以表示整个序列。

构造离散序列x1t =[5 2 8 8 1 1 2]; x2t =[8 2 3 1 5 8 9 3 5 5];分别做x1t 与x2t 的卷积时域图和幅值谱图，将所得到的频域图做快速傅里叶逆变换（FFT）得到时域图，发现该时域图与原时域相同，满足时域卷积特性。

(2)连续时间信号的卷积：构造信号f1t=ut-ut-3,f1t=e-tut,ft=f1t*f1t 。时域卷积对应频域相乘，满足时域卷积特性。将得到的频域图做快速傅里叶逆变换（FFT）可以得到原卷积信号的时域图。

1. 频域卷积特性

x1tx2t↔X1f*X2(f)

构造离散序列x1t =[5 2 8 8 1 1 3 2]; x2t =[8 2 3 1 5 8 9 5];分别作x1tx2t 的时域图和幅值谱图，再分别对x1t 和x2t 做快速傅里叶变换（FFT）得到幅值谱X1(f)和X2(f),画出X1f*X2(f) 幅值谱图，再将X1f*X2(f) 作快速傅里叶逆变换得到时域图，发现与x1tx2t 的时域图相同，符合频域卷积特性。