关于离散傅里叶变换DFT的导出

目录：

1.计算机进行数字信号处理的特点

2.离散复指数函数的特点

3.离散傅里叶级数(DFS)

4.离散傅里叶变换(DFT)

关于计算机处理数字信号的特点

1)计算机只能处理有限长度的信号

2)计算机只能处理离散的时域和频域信号

3)计算机进行频域分析时，把需要处理的时域信号当作一个周期内的信号

离散指数函数基的特点

1.$e^{j{w}_{0}n}$的角频域变换$2\pi$后将等于原信号

证：
$e^{j(w_0+2\pi )n}=e^{j{w}_{0}n}{e}^{j2\pi }=e^{j{w}_{0}n}$
这说明$w_0$的变换范围只能在$0-2\pi$，超出范围就变成原信号了

2.$e^{j{w}_{0}n}$在一定条件下是周期信号

证：
若$e^{j{w}_{0}n}$是周期信号，且周期是$N$,那么有

$e^{j{w}_0(n+N)}=e^{j{w}_{0}n}$=$e^{j{w}_{0}n}{e}^{j{w}_{0}N}$.

则$e^{j{w}_{0}N}=1$,那么$w_{0}N=2\pi m(m=0,1,2...)$

所以当且仅当$w_0=\frac{2\pi m}{N}$时，信号$e^{j{w}_{0}n}$为周期信号(即$w_0$是$\pi$的K倍，K为有理数)

且周期为$N$，基波频率为$\frac{2\pi m}{N}$

举例：
①$x(t)=e^{2\pi n/3}$

显然周期N=3，基波频率为$\frac{2\pi }{3}$，此时$m=1$；

②$x(t)=e^{3\pi n/4}$

此时周期N=8,因为$\frac{3\pi }{4}\ast 8=2\pi \ast 3$，基波频率为$\frac{3\pi }{4}$

以上性质对正余弦函数仍适应。

仿真分析：

n=0:1:69;
w_0=2*pi/3;

subplot(2,1,1);
plot(n,sin(w_0*n),'linewidth',1);
title('x=sin(w_0*n)');
subplot(2,1,2);
plot(n,sin((w_0+2*pi)*n),'linewidth',1);
title('x=sin((w_0+2pi)*n)');

2.离散复指数函数集

设基波频域为$w_0$，具有公共周期$N$(不是最小周期)的一组离散复指数函数集合为${\Phi }_k(n)$,$k$代表第k个函数。

显然${\Phi }_k(n)=e^{jk{w}_{0}n},k=0,\pm 1，\pm 2$

因为${\Phi }_{k+N}(n)=e^{j(k+N)w_{0}n}=e^{jk{w}_{0}n}{e}^{jN{w}_{0}n}$

根据周期复指数函数的特点，$e^{jN{w}_{0}n}=e^{j2\pi mn}=1$

所以${\Phi }_k(n)$的不同函数的个数为$N$,且${\Phi }_{k+N}={\Phi }_k$

离散傅里叶级数(DFS)

对于周期为$N$的离散信号$x(n)$，根据傅里叶级数进行展开有：

$x(n)=\sum _{k=0}^{N}X(k)e^{jk{w}_{0}n}$

$X(k)=\frac{1}{N}X(k)e^{-jk{w}_{0}n}$

因为只有$N$个不相同的信号，所以求和范围是0-N

这个变化对的特点是时域周期离散，频域也是周期离散

离散傅里叶变换(DFT)

DFT表达的是周期无限的离散信号，其频谱也是周期无限的，而计算机不能处理无限长的信号，因

此，常常把要处理的信号$x(n)$看成是某个周期信号$\overline{x(n)}$在一个周期$N$内的信号，

$N$为$x(n)$的长度。频谱也是$\overline{x(n)}$一个周期$N$的部分。

结果与matlab一致。