1、快速傅里叶变换的实现
什么是傅里叶定理?
法国科学家傅里叶提出,任何一条周期性曲线,无论多么跳跃或不规则,都能表示成一组光滑正弦曲线叠加之和。
什么是傅里叶变换?
傅里叶变换即是把一条周期性曲线拆解成一组光滑正弦曲线的过程。
傅里叶变换的目的是可将时域(即时间域)上的信号转变为频域(即频率域)上的信号,随着域的不同, 对同一个事物的认知角度也随之改变,因此在时域中某些不好处理的地方,在频域就可以较为简单地处理。这就可以大量减少处理信号存储量。
假设有一时间域函数:y = f(x),根据傅里叶的理论它可以被分解为一系列正弦函数的叠加,他们的振幅A,频率ω或初相位φ不同:
所以傅里叶变换可以把一个比较复杂的函数转换为多个简单函数的叠加,看问题的角度也从时间域转到了频率域,有些问题处理起来就会比较简单。
傅里叶变换相关函数
导入快速傅里叶变换所需模块
import numpy.fft as nf
通过采样数域采样周期求得傅里叶变换分解所得曲线的频率序列
freqs = nf.fftfreq(采样熟练,采样周期)
通过原函数值得序列经过快速傅里叶变换得到一个复数数组,复数的模代表的是振幅,复数的辐角代表初相位
nf.fft(原函数序列值) -> 目标函数序列值(复数数组)
通过一个复数数组(复数的模代表振幅,复数的辐角代表初相位)经过逆傅里叶变换得到合成的函数值数组
nf.ifft(目标函数值序列(复数)) -> 原函数值序列