傅里叶变换 - 1

快速傅里叶变换需要用到单位根的性质，在此先介绍下复数单位根及其性质

1. 复数单位根及其性质
$$w^n=1$$

   符合条件的w值就是 复数单位根。将单位根在图中表示，横轴表示实部，纵轴表示虚部，python代码如下

# 绘制复数单位根图像
import numpy as np
import matplotlib.pyplot as plot


plot.axes(projection='polar')
plot.title('root of unit')
n = 8  #单位根的个数
rads = np.arange(0, 2*np.pi, 2*np.pi/n)
for rad in rads:
    plot.polar(rad, 1, 'o')
plot.show()

图片为：

复数单位跟的性质
$$w_n^k=cos\left(\frac{2\pi k}{n}\right)-i\ast sin\left(\frac{2\pi k}{n}\right),k=0,1,2,...,n-1$$
    $w_n^k$ 中，n表示单位根的个数，k表示第k个单位根。i对应虚数单位
1）当n>0满足 $$w_{dn}^{dk}=w_n^k$$ 因为
$$w_{dn}^{dk}=cos(\frac{2\pi dk}{dn})-i\ast sin(\frac{2\pi dk}{dn})=cos(\frac{2\pi k}{n})-i\ast sin(\frac{2\pi k}{n})=w_n^k,k\ge 0,n>0,d>0$$
2) 如果n是大于0的偶数，那么就有$$w_n^{\frac{n}{2}}=-1$$因为
$$w_n^{\frac{n}{2}}=cos\left(\frac{2\pi \frac{n}{2}}{n}\right)-i\ast sin\left(\frac{2\pi \frac{n}{2}}{n}\right)=cos\left(\pi \right)-i\ast sin\left(\pi \right)=-1-0=-1$$
3) $${\left(w_n^k\right)}^{{}^2}=w_n^{2k}$$当n是偶数时，满足
$${\left(w_n^{k+\frac{n}{2}}\right)}^2=w_n^{2k+n}=w_n^{2k}={\left(w_n^k\right)}^2$$ 推导:
   a :
$$\begin{gathered} {\left(w_n^k\right)}^2={\left[cos\left(\frac{2\pi k}{n}\right)-i\ast sin\left(\frac{2\pi k}{n}\right)\right]}^2 \\ =co{s}^2\left(\frac{2\pi k}{n}\right)-si{n}^2\left(\frac{2\pi k}{n}\right)-2\ast i\ast sin\left(\frac{2\pi k}{n}\right)\ast cos\left(\frac{2\pi k}{n}\right) \\ =cos\left(\frac{4\pi k}{n}\right)-i\ast sin\left(\frac{4\pi k}{n}\right)=w_n^{2k} \end{gathered}$$
   b：
$$\begin{gathered} {\left(w_n^{k+\frac{n}{2}}\right)}^2=w_n^{2k+n} \\ =cos\left(\frac{2\pi (2k+n)}{n}\right)-i\ast sin\left(\frac{2\pi (2k+n)}{n}\right) \\ =cos\left(\frac{4\pi k}{n}\right)-i\ast sin\left(\frac{4\pi k}{n}\right)=w_n^{2k}={\left(w_n^k\right)}^2 \end{gathered}$$
4)
$$\sum _{j=0}^{n-1}{\left(w_n^k\right)}^j=0$$根据等比数列求和公式
$$\sum _{j=0}^{n-1}{p}^j=\frac{p^n-1}{p-1}$$
性质4）公式可展开为
$$\sum _{j=0}^{n-1}{\left(w_n^k\right)}^j=\frac{{\left(w_n^k\right)}^n-1}{w_n^k-1}=\frac{w_n^{kn}-1}{w_n^k-1}=\frac{1-1}{w_n^k-1}=0$$
5)$$w_n^{k+j}=w_n^k\ast w_n^j$$推导：
$$\begin{gathered} w_n^k\ast w_n^j=\left[cos\left(\frac{2\pi k}{n}\right)-i\ast sin\left(\frac{2\pi k}{n}\right)\right]\ast \left[cos\left(\frac{2\pi j}{n}\right)-i\ast sin\left(\frac{2\pi j}{n}\right)\right] \\ =cos\left(\frac{2\pi k}{n}\right)\ast cos\left(\frac{2\pi j}{n}\right)-sin\left(\frac{2\pi k}{n}\right)\ast sin\left(\frac{2\pi j}{n}\right) \\ -i\ast \left[sin\left(\frac{2\pi k}{n}\right)\ast cos\left(\frac{2\pi j}{n}\right)+cos\left(\frac{2\pi k}{n}\right)\ast sin\left(\frac{2\pi j}{n}\right)\right] \end{gathered}$$ 根据三角公式积化和差，有
            $cos\left(\frac{2\pi k}{n}\right)\ast cos\left(\frac{2\pi j}{n}\right)=\frac{1}{2}\left[cos\left(\frac{2\pi k+2\pi j}{n}\right)+cos\left(\frac{2\pi k-2\pi j}{n}\right)\right]$
            $sin\left(\frac{2\pi k}{n}\right)\ast sin\left(\frac{2\pi j}{n}\right)=-\frac{1}{2}\left[cos\left(\frac{2\pi k+2\pi j}{n}\right)-cos\left(\frac{2\pi k-2\pi j}{n}\right)\right]$
          $sin\left(\frac{2\pi k}{n}\right)\ast cos\left(\frac{2\pi j}{n}\right)=\frac{1}{2}\left[sin\left(\frac{2\pi k+2\pi j}{n}\right)+sin\left(\frac{2\pi k-2\pi j}{n}\right)\right]$
          $cos\left(\frac{2\pi k}{n}\right)\ast sin\left(\frac{2\pi j}{n}\right)=\frac{1}{2}\left[sin\left(\frac{2\pi k+2\pi j}{n}\right)-sin\left(\frac{2\pi k-2\pi j}{n}\right)\right]$
推导公式可化简为：

$$\begin{gathered} =\frac{1}{2}\left[cos\left(\frac{2\pi k+2\pi j}{n}\right)+cos\left(\frac{2\pi k-2\pi j}{n}\right)\right]+\frac{1}{2}\left[cos\left(\frac{2\pi k+2\pi j}{n}\right)-cos\left(\frac{2\pi k-2\pi j}{n}\right)\right] \\ -i\ast \left[\frac{1}{2}\left[sin\left(\frac{2\pi k+2\pi j}{n}\right)+sin\left(\frac{2\pi k-2\pi j}{n}\right)\right]+\frac{1}{2}\left[sin\left(\frac{2\pi k+2\pi j}{n}\right)-sin\left(\frac{2\pi k-2\pi j}{n}\right)\right]\right] \\ =cos\left(\frac{2\pi k+2\pi j}{n}\right)-i\ast sin\left(\frac{2\pi k+2\pi j}{n}\right)=w_n^{k+j} \end{gathered}$$

本节主要介绍复数单位根及其性质，接下来介绍傅里叶变换及用python实现傅里叶变换
第一次用csdn的公式编辑器，有的公式编辑出来不美观，不知道怎么解决。特别是公式换行怎么让等号对齐。