matlab fftshift_matlab下实现FFT信号分析

奈奎斯特定理：

对于信号分析，我们需要掌握一个定理：

在进行模拟/数字信号的转换过程中，当采样频率 fs.max大于信号中最高频率 fmax 的 2 倍时(fs.max > 2fmax)，采样之后的数字信号完整地保留了原始信号中的信息，一般实际应用中保证采样频率为信号最高频率的2.56～4倍；采样定理又称奈奎斯特定理或香农采样定理。

利用matlab做频谱分析前我们需要了解分析过程中的一些基础知识，

matlab中的 fft 函数用法、fftshift 函数的用法

函数 1  fft ：

作用：快速傅里叶变换。

语法：

Y = fft(X)

Y = fft(X,n)

Y = fft(X,n,dim)

语法：

Y = fft(X) 用快速傅里叶变换 (FFT) 算法计算 X 的离散傅里叶变换 (DFT)。

• 如果 X 是向量，则 fft(X) 返回该向量的傅里叶变换；

• 如果 X 是矩阵，则 fft(X) 将 X 的各列视为向量，并返回每列的傅里叶变换；

• 如果 X 是一个多维数组，则 fft(X) 将沿大小不等于 1 的第一个数组维度的值视为向量，并返回每个向量的傅里叶变换；

Y = fft(X,n) 返回 n 点 DFT。如果未指定任何值，则 Y 的大小与 X 相同。

• 如果 X 是向量且 X 的长度小于 n，则为 X 补上尾零以达到长度 n；

• 如果 X 是向量且 X 的长度大于 n，则对 X 进行截断以达到长度 n；

• 如果 X 是矩阵，则每列的处理与在向量情况下相同；

• 如果 X 为多维数组，则大小不等于 1 的第一个数组维度的处理与在向量情况下相同；

函数 2 fftshift：

作用：将零频分量移到频谱中心。

语法：

Y = fftshift(X)

Y = fftshift(X,dim)

说明：

Y = fftshift(X) 通过将零频分量移动到数组中心，重新排列傅里叶变换 X。

• 如果 X 是向量，则 fftshift 会将 X 的左右两半部分进行交换；

• 如果 X 是矩阵，则 fftshift 会将 X 的第一象限与第三象限交换，将第二象限与第四象限交换；

• 如果 X 是多维数组，则 fftshift 会沿每个维度交换 X 的半空间；

应用举例

现有输入信号 2Hz，ADC 采样频率为 500Hz ，满足采样定理的要求，目标：对信号进行频谱分析，找出信号在频域上的噪声分量，并设计对应滤波器滤除噪声，提高信噪比。

新建脚本文件FFT.m，脚本内容如下图 1：

程序第 3 行定义了信号的采样频率为 500Hz，第 4 行定义了采样周期，第 6 行导入了脚本路径下的txt数据，第 7 行计算导入数据的长度，第 8 行是计算数据的频率范围，，第 9~11 行绘出数据的时域信号波形图；

程序第 14 ~20 行计算绘制 FFT频谱图。

图 1

程序运行结果如下图 2 所示，我们可以根据图 2 所展示的频谱信息设计相应的滤波电路或软件算法去除噪声。另外从图 2 中可以看出，FFT 变换得到的频谱是左右对称的，因此，我们只需要其中一边就能获得信号的所有信息。链接中演示了怎样将双边频谱图转换为单边频谱图，可参考学习 https://zhuanlan.zhihu.com/p/42893470

图 2

●贴片晶振的PCB layout需要注意哪些？

●CPU中的程序是怎么运行起来的？

●MOSFET栅极前 100 Ω 电阻有什么用？

●模拟电路设计的这12个细节，往往起着重要作用

●STM32：从菜鸟到牛人就是如此简单！

●什么是PID？讲个故事，秒懂！

●485通信自动收发电路 历史上最详细的解释