【Python信号处理】FFT变换后的坐标变换以及绘图

$$假设：f1=10;f2=20;f3=30$$
$$y=sin(2\ast pi\ast f1\ast t)+2\ast sin(2\ast pi\ast f2\ast t)+3\ast sin(2\ast pi\ast f3\ast t)$$
$$FrequencySample(fs)=128,Length(T)=5；注意此处的采样点数要满足奈奎斯特准则$$
$$SummationLength(N)=fs\ast T=128\ast 5$$
FFT后的频域转化公式为
$$F_i=\sum _{1<=i<=N}\frac{i}{N}\ast fs$$
FFT后的幅值转化公式为, Ai为傅里叶变换后的第i个值
$$A_i=\sum _{1<=i<=N}|A_i|\frac{2}{N}$$

双边傅里叶变换频率

from matplotlib.pyplot import plot, show, subplot
from numpy import pi, sin, linspace
from numpy.fft import fft, fftfreq


fs = 128
T = 5
f1 = 10
f2 = 20
f3 = 30
t = linspace(0, T, T*fs)
y = sin(2*pi*f1*t) + 2*sin(2*pi*f2*t) + 3*sin(2*pi*f3*t)
x = fftfreq((len(t))) * fs
subplot(121)
plot(t, y)
subplot(122)
plot(x, abs(fft(y)) / len(t) * 2)
show()

单边傅里叶变换频率

from matplotlib.pyplot import plot, show, subplot
from numpy import pi, sin, linspace
from numpy.fft import fft, fftfreq


fs = 128
T = 5
f1 = 10
f2 = 20
f3 = 30
t = linspace(0, T, T*fs)
y = sin(2*pi*f1*t) + 2*sin(2*pi*f2*t) + 3*sin(2*pi*f3*t)
x = fftfreq((len(t))) * fs
subplot(121)
plot(t, y)
subplot(122)
plot(x[:len(x)//2], abs(fft(y))[:len(x)//2] / len(t) * 2)
show()