DFT

傅里叶分析方法是信号分析的最基本方法傅里叶变换是傅里叶分析的核心经过它把信号从时间域变换到频率域进而研究信号的频谱结构和变化规律。

可是它的致命缺点是计算量太大时间复杂度过高当采样点数过高的时候计算缓慢由此出现了DFT的快速实现即下面的快速傅里叶变换FFT。

1.2.快速傅里叶变换FFT

计算量更小的离散傅里叶的一种实现方法。

快速傅氏变换FFT是离散傅氏变换的快速算法它是根据离散傅氏变换的奇、偶、虚、实等特性对离散傅立叶变换的算法进行改进获得的。

采样频率也称为采样速度或者采样率定义了每秒从连续信号中提取并组成离散信号的采样个数它用赫兹Hz来表示。

采样频率的倒数是采样周期或者叫做采样时间它是采样之间的时间间隔。

通俗的讲采样频率是指计算机每秒钟采集多少个信号样本。

采样定理

又称香农采样定理奈奎斯特采样定理是信息论特别是通信与信号处理学科中的一个重要基本结论。

采样定理指出若是信号是带限的而且采样频率高于信号带宽的两倍那么原来的连续信号能够从采样样本中彻底重建出来。

1.4.如何理解采样定理

在对连续信号进行离散化的过程当中不免会损失不少信息就拿一个简单地正弦波而言若是我1秒内就选择一个点很显然损失的信号太多了光着一个点我根本不知道这个正弦信号究竟是什么样子的天然也没有办法根据这一个采样点进行正弦波的还原很明显我采样的点越密集那越接近原来的正弦波原始的样子天然损失的信息越少越方便还原正弦波。

采样定理说明采样频率与信号频率之间的关系是连续信号离散化的基本依据。

它为采样率创建了一个足够的条件该采样率容许离散采样序列从有限带宽的连续时间信号中捕获全部信息。

test_fft_1.py

mpl.rcParams[axes.unicode_minus]

False

采样点选择1400个由于设置的信号频率份量最高为600赫兹根据采样定理知采样频率要大于信号频率2倍

因此这里设置采样频率为1400赫兹即一秒内有1400个采样点同样意思的

fft_y

normalization_y[range(int(N/2))]

#因为对称性只取一半区间单边频谱

plt.title(双边振幅谱(未求振幅绝对值),fontsize9,colorblack)

plt.subplot(233)

plt.title(双边振幅谱(未归一化),fontsize9,colorred)

plt.subplot(234)

plt.title(双边相位谱(未归一化),fontsize9,colorviolet)

plt.subplot(235)

plt.title(双边振幅谱(归一化),fontsize9,colorgreen)

plt.subplot(236)

plt.plot(half_x,normalization_half_y,blue)

plt.title(单边振幅谱(归一化),fontsize9,colorblue)

plt.show()运行

傅里叶定理指出任何频率为f0的周期信号都可以通过将频率为f0,2f0,3f0,4f0,5f0等的“正弦波”正弦波相加而精确地构建。

将周期时域信号分割为正弦波称为傅里叶分析。

“傅里叶级数”中的每个正弦曲线的特征在于频率振幅以及阶段f0被称为基频。

test_sawtooth.py

sawtooth()生成了频率为f200Hz的锯齿信号持续时间为2秒。

import