二维离散傅里叶变换最详手算

前置知识

1. 1 复数

​ 定义(符号表示)：

​ a + bi （这里a，b是实数，i是虚数单位，a叫做复数的实部，b叫做复数的虚部。可以认为实数a等价于复数a+0i）

​ 几何(复平面)：

​ 水平的实轴与垂直的虚轴建立起来的复数的几何表示

​ 其他理解：

 对于方程 x^2 = 1*x*x = -1   
 我们可以把*x看作是一种变换，经过两次变换，1变成了-1 
 如果我们把这种变换看成是“逆时针旋转90%”，就能够实现1变为-1
 也就是说我们可以把 *i 看作是逆时针旋转90%

​ 如图

​ 复平面上一个点乘以（-1+7i）的几何意义：

​ 模长增加\srqt(50)，方向逆时针旋转arctan(7/-1) = 98.17°

​ 两个复数相乘的结果就是：让它们的模长相乘得到最终的模长，让它们的幅角相加得到最终的幅角

参考文章：https://zhuanlan.zhihu.com/p/48392273

1.2 e的美妙

​ e = 2.718281828459…

​ 指数增长模型

​ 假设你有1元钱存在银行里，银行的利率为100%，若银行一年付一次利息，一年后你可以拿到1元本金+1元利息共两元

​ 假设银行利率不变但利息改为半年一付，那么在半年时候把0.5的利息再次存入银行(复利)，一年后能拿到1元本金+1+0.25共2.25元

​ 假设利率还是不是但利息改为四月一付，拿到利息后进行复利，那么一年后能拿到2.37元

​ 计算公式 // r为增长率，x为增长周期，Q为初始值的Q倍
$$Q=(1+r{)}^x$$
​ 复利模型

​ 假设利息连续不断的付，进行连续复利，那么利润的这个天花板就是e！如下图

​ 高等数学的重要极限
$$e=\underset{n\to \infty }{\lim }(1+\frac{1}{n}{)}^n$$
​ 参考文章：https://zhuanlan.zhihu.com/p/48391055

1. 2 欧拉公式

​ 定义：

$$e^{ix}=\cos x+i\sin x$$

​ 欧拉恒等式x=Π：
$$e^{i\pi }=-1,e^{i\pi }+1=0$$
​

​ 推导：

​ 采用复数cos(x)+i*sin(x)可以描述单位圆周上的点的位置或运动轨迹

​

​ 自然底数e的定义式
$$e=\underset{n\to \infty }{\lim }(1+\frac{1}{n}{)}^n=\underset{n\to \infty }{\lim }(1+\frac{100\%}{n}{)}^n$$
​ 当增长率为虚数时候
$$e^i=\underset{n\to \infty }{\lim }(1+\frac{100\%\ast i}{n}{)}^n=1\ast \underset{n\to \infty }{\lim }(1+\frac{100\%\ast i}{n}{)}^n$$
​ 我们可以把上式拆分，看作是1进行了n次的变换，每一次的变换的模长增量都是1，角度增量都是1/n。

​

​ 每一次的旋转如图

​

即最终转动的弧度为1。由此可得

​
$$e^i=1\ast \underset{n\to \infty }{\lim }(1+\frac{100\%\ast i}{n}{)}^n=\cos 1+i\ast \sin 1$$
​ 更一般地对于e^(ix)

​ 即最终旋转的弧度为x。由此可得

​
$$e^{ix}=1\ast \underset{n\to \infty }{\lim }(1+\frac{100\%\ast ix}{n}{)}^n=\cos x+i\ast \sin x$$
​ 证毕。

​ 其他等式
$$\begin{gathered} \sin x=\frac{e^{ix}-e^{-ix}}{2i} \\ \cos x=\frac{e^{ix}+e^{-ix}}{2} \end{gathered}$$
参考文章：https://zhuanlan.zhihu.com/p/48392958

一 定义

二 例题

2.1 定义展开法

​ 题目
$$f(x,y)=\cos (\frac{\pi }{4}x),x=[0,511],计算其二维离散傅里叶变换F(u,v)$$
​ 解：
$$\begin{array}{rl}& 根据欧拉公式有f=\frac{e^{i\frac{\pi }{4}x}+e^{-i\frac{\pi }{4}x}}{2},设M=N=512\\ & F(u,v)=\sum _{x=0}^{M-1}\sum _{y=0}^{N-1}f(x,y)e^{-j2\pi (\frac{u}{m}x+\frac{v}{N}y)}\\ & =\frac{1}{2}\sum _{x=0}^{M-1}\sum _{y=0}^{N-1}(e^{j\frac{\pi }{4}x}+e^{-j\frac{\pi }{4}x})\ast e^{-j2\pi (\frac{u}{m}x+\frac{v}{N}y)}\\ & =\frac{1}{2}\sum _{x=0}^{M-1}\sum _{y=0}^{N-1}(e^{j(\frac{\pi }{4}x-\frac{2\pi u}{M}x-\frac{2\pi v}{N}y)}+e^{-j(\frac{\pi }{4}x+\frac{2\pi u}{M}x+\frac{2\pi v}{N}y)})\\ & =\frac{1}{2}\sum _{x=0}^{M-1}\sum _{y=0}^{N-1}(e^{j\frac{2\pi }{512}[(64-u)x-vy]}+e^{-j\frac{2\pi }{512}[(64+u)x+vy]})\\ & =\frac{1}{2}\sum _{x=0}^{M-1}\sum _{y=0}^{N-1}(\cos \frac{2\pi }{512}[(64-u)x-vy]+j\sin \frac{\pi v}{512}[(64+u)x+vy]+\cos \frac{2\pi }{512}[(64+u)x+vy]+-(j)\sin \underset{512}{}\\ & \end{array}$$