最近在看傅里叶转换(FT)和快速傅里叶转换(FFT),感觉很难,包括号称“小学生都能看懂的快速傅里叶转换”,也是一样晦涩难懂,当然,这不排除现在的小学生智力都很强,而我的智商已经在走下坡路了。言归正传,我也来实现一下傅里叶转换(FT),下一节再实现快速傅里叶转换(FFT)。
傅里叶转换是非常有用的一个算法,在声学、信号处理等领域有广泛的应用,例如核磁共振得到的数据是很难看出特征的,用傅里叶转换之后,就很容易看出谱图的特征。两个多项式的乘法也称之为“卷积”,“卷积”操作是AI中非常常见且重要的操作,傅里叶转换也可以用于两个多项式的乘法。傅里叶转换的时间复杂度是O(N2),快速傅里叶转换的时间复杂度为O(N * logN)。这两者的差距就相当于冒泡排序和快速排序的差别,所以我们很有必要学习一下这两种转换。在这里,我们就用两个多项式的乘法为例介绍傅里叶转换。
傅里叶转换 (FT)
有三个数a, b, c。其中a = 111, b = 11111,c = a * b,我们很容易用计算器求出c = 1233321。但是如果a和b都是非常大的数,例如都有100位,我们就没办法用计算器求解了。如果用C语言的编程,则会超出范围。
int a = 1111111111111111111111111111111111111111111; //假设有100位
int b = 1111111111111111111111111111111111111111111; //假设有100位
int c = a * b;
对于这样的问题,我们该怎么办呢?这就要用到傅里叶转换了:用数组arr1存放a中的每一位数据;用数组arr2存放b中的每一位数据。然后让arr1中的每一个元素与arr2中的每一个元素相乘,并累加。代码如下:
#include<stdio.h>
#include<stdlib.h>
int * arr;
//傅里叶转换,时间复杂度为O(N * M)
void ft(int * arr1, int * arr2, int n, int m)
{
int a = 0;
int b = 0;
for(int i = 0; i < n; i++)<