一、介绍
本人在一段程序中想实现快速傅里叶变换(FFT)功能,但在CSDN上找了一圈,发现此类文章有点少。
不过,在查找了不少书本和网站的相关资料后,还是将功能实现了。现在,特意分享所写的部分代码,供各位参考参考。
二、代码
在放出代码之前,我得先说明一下。
你想要的结果决定你的傅里叶变换应该怎么写。每个人程序的傅里叶代码都是有所不同的,不过大致框架可以一样。这也是为什么不同文章所给出的代码都不同。
关于我这一段代码,主要的功能就是看出频率,以及对应的幅值,从而以此实现区分两类动作。话我就不多说了,到时我写另外一篇详细一点的,从易到难争取都覆盖到。
%载入数据
drinking_acc_x = drinking(:,1);
%FFT变换 X_d
S = drinking_acc_x;
Fs = 400; %采样频率Hz
N