最近学校电赛院队招新,出的招新题就是低频示波器的。之前一直没有弄懂FFT,借着这次机会实现了一下。

  • FFT原理详解

  FFT,就是快速傅里叶变换,这个操作能够将时域信号转化成频域信号,然后对信号进行分析

  这样说可能有点抽象。讲细点就是指能够直观的看出来目标信号的频率是多少。x轴坐标本来是表示时间,FFT之后变成了表示频率,就是这个意思

  对于信号处理,FFT之后的结果,波峰一般会出现在我们希望测得信号的频率附近(十分相近)

  • 官方文件解释

stm32官方给了几个用于处理FFT的文件,如图所示:

其中有两个汇编文件两个头文件:汇编文件是定义了FFT的计算函数,我们直接调用即可

cr4_fft_1024_stm32.s是包含了计算1024个点的FFT的函数的汇编文件,另一个汇编文件同理

stm32_dsp.h里面有关于FFT处理函数的声明,我们包含了这个头文件之后直接调用函数即可

  • 算法解释
 1 //进行FFT运算等操作
2 void FFT_Wave(void)
3 {
4 u16 i;
5 float mid_value;
6 while(!ADC_flag)
7 {
8 LED1 = !LED1;
9 delay_ms(100);
10 }
11 ADC_flag = 0;
12
13 //获取最大值最小值
14 adc_value_max = adc_value_min = ADC_buff[1];
15 for(i = 0;i < NPT;i++)
16 {
17 //寻找最大值最小值
18 if(ADC_buff[i] >= adc_value_max)
19 {
20 adc_value_max = ADC_buff[i];
21 }
22 if(ADC_buff[i] <= adc_value_min)
23 {
24 adc_value_min = ADC_buff[i];
25 }
26 //先清空数组
27 fftin[i] = 0;
28 //移位,让后面16位为虚部
29 fftin[i] = ((s16)ADC_buff[i] << 16);
30 }
31 cr4_fft_1024_stm32(fftout,fftin,1024);//FFT
32 GetPowerMag();
33 //计算电压值
34 Vpp_true = (adc_value_max - adc_value_min) * 3.3 / 4096.0;//获得Vpp值
35 mid_value = (adc_value_max + adc_value_min) / 2;
36 for(int i = 0;i < NPT;i++)
37 {
38 if(ADC_buff[i] > mid_value)
39 {
40 rect_duty++;
41 }
42 }
43 rect_duty = rect_duty / 1024 * 100;
44 }

这是FFT的主体函数

第一步我们先要等待ADC采集完成,将数据存入数组当中准备进行处理

第二步是在采样值当中寻找最大值和最小值(遍历数组即可)

第三步是对数组进行移位处理(前面的是实部,后面的是虚部,由于我们采集到的电压都是实数,所以虚部都置0)

第四步是使用ST官方提供的函数进行FFT运算,得到运算之后的数组

第五步是根据频谱查找我们信号所对应的频率,也就是对频谱图当中所有的频率进行幅值的比较,找出幅值最大时所对应的频率,即为我们所需要测量的频率,其他的都可以看作噪声

在我们找到该频率之后,不能立刻输出,要与ADC的采样率相乘再除以1024,之后才能得到我们想要的信号频率

GetPowerMag函数定义如下:

 1 void GetPowerMag(void)
2 {
3 s16 lX,lY;
4 u32 i;
5 float maxmag;
6 for(i = 0;i < NPT / 2;i++)
7 {
8 lX = (fftout[i] << 16) >> 16;
9 lY = (fftout[i] >> 16);
10 float X = 1024 * ((float)lX) / 32768;
11 float Y = 1024 * ((float)lY) / 32768;
12 float mag = sqrt(X * X + Y * Y) / 1024;
13 FFT_Mag[i] = (u32)(mag * 65536);
14 }
15 FFT_Mag[0] >>= 1;//频谱图第一个是直流分量,无需乘2
16 for(int i = 0;i < NPT / 2;i++)
17 {
18 if((maxmag < FFT_Mag[i]) && (i != 0))
19 {
20 maxmag = FFT_Mag[i];
21 temp = i;
22 }
23 }
24 F_hz = temp * sampling_rate / 1024.0;
25 }

至此,我们就得到了我们所需信号的频率

鉴于本小白能力有限,如果有纰漏或改进之处,欢迎指正

特别提醒:ADC采样率应遵循奈奎斯特采样定理!采样率不是越高越好!

stm32F103RCT6使用FFT运算分析波形详解(非常新手)的更多相关文章

  1. 【Java】HashMap源码分析——常用方法详解

    上一篇介绍了HashMap的基本概念,这一篇着重介绍HasHMap中的一些常用方法:put()get()**resize()** 首先介绍resize()这个方法,在我看来这是HashMap中一个非常 ...

  2. valgrind和Kcachegrind性能分析工具详解

    一.valgrind介绍 valgrind是运行在Linux上的一套基于仿真技术的程序调试和分析工具,用于构建动态分析工具的装备性框架.它包括一个工具集,每个工具执行某种类型的调试.分析或类似的任务, ...

  3. MATLAB实例:PCA(主成成分分析)详解

    MATLAB实例:PCA(主成成分分析)详解 作者:凯鲁嘎吉 - 博客园 http://www.cnblogs.com/kailugaji/ 1. 主成成分分析 2. MATLAB解释 详细信息请看: ...

  4. 16S 基础知识、分析工具和分析流程详解

    工作中有个真理:如果你连自己所做的工作的来龙去脉都讲不清楚,那你是绝对不可能把这份工作做好的. 这适用于任何行业.如果你支支吾吾,讲不清楚,那么说难听点,你在混日子,没有静下心来工作. 检验标准:随时 ...

  5. JVM性能分析工具详解--MAT等

    获得堆转储文件 巧妇难为无米之炊,我们首先需要获得一个堆转储文件.为了方便,本文采用的是 Sun JDK 6.通常来说,只要你设置了如下所示的 JVM 参数: -XX:+HeapDumpOnOutOf ...

  6. nmon分析与详解

    1.命令安装 1.查看liunx版本版本x86_64_14i 目录:cd /nmon/logs/ 版本x86_64_14i [root@localhost u06]# cd / [root@local ...

  7. RTP协议分析和详解

    一.RTP协议分析 第1章.     RTP概述 1.1.  RTP是什么 RTP全名是Real-time Transport Protocol(实时传输协议).它是IETF提出的一个标准,对应的RF ...

  8. webpack4配置详解之新手上路初探

    前言 经常会有群友问起webpack.react.redux.甚至create-react-app配置等等方面的问题,有些是我也不懂的,慢慢从大家的相互交流中,也学到了不少. ​ 今天就尝试着一起来聊 ...

  9. [拉格朗日反演][FFT][NTT][多项式大全]详解

    1.多项式的两种表示法 1.系数表示法 我们最常用的多项式表示法就是系数表示法,一个次数界为\(n\)的多项式\(S(x)\)可以用一个向量\(s=(s_0,s_1,s_2,\cdots,s_n-1) ...

随机推荐

  1. SpringBoot静态资源配置访问上传文件

    使用SpringBoot项目开发上传文件的代码时,如果想访问已上传的文件,但处于测试阶段,而不想配置Nginx服务并启动这么繁琐,那么配置以下代码即可 @Override public void ad ...

  2. IIS MVC 发布错误 403.14-Forbidden Web 服务器被配置为不列出此目录的内容

     转:http://blog.csdn.net/csethcrm/article/details/37820135 IIS MVC 发布错误 403.14-Forbidden Web 服务器被配置为不 ...

  3. Tcp三次握手四次挥手个人学习

    最近想跳槽,学习了tcp中的三次握手与四次挥手,特意记录下,加深记忆 SYN 代表请求创建连接 FIN 表示请求关闭连接 ACK 代表确认接受,不管是三次握手还是四次分手,在回应的时候都会加上ACK= ...

  4. Java_lambda表达式之"stream流学习,Map学习,collect学习,Conllectors工具类学习"

    Lambda表达式学习 对List<Integer> userIdList = UserList.stream().map(User::getUserId).collect(Collect ...

  5. 遇到问题之“postman报Unsupported Media Type: Content type 'text/plain;charset=UTF-8' not supported”

    postman报Unsupported Media Type: Content type 'text/plain;charset=UTF-8' not supported postman之所以报Uns ...

  6. Spring官宣网传大漏洞,并提供解决方案

    Spring沦陷了!这样的标题这几天是不是看腻了?然而,仔细看看都是拿着之前的几个毫不相干的CVE来大吹特吹.所以,昨天发了一篇关于最近网传的Spring大漏洞的文章,聊了聊这些让人迷惑的营销文.以及 ...

  7. asp.net 可视化操作(一)——asp.net安装与使用

    目录 安装 创建网页 设计网页 运行 vs 2019安装asp.net 1.安装 打开vs,选择继续但无需代码 -->工具–>获取工具和功能 勾选如下选项后,点击关闭 点击更新等待安装完成 ...

  8. 模拟web服务器http请求应答

    我们在浏览器打开网页,其实是向远端服务器提出页面发送请求,远端服务器在接到请求后,就开始执行请求页面的程序文件,然后将执行结果通过html格式,发送到你的浏览器,再显示出来.以下用百度(www.bai ...

  9. [译] 沙箱中的间谍 - 可行的 JavaScript 高速缓存区攻击

    原文 The Spy in the Sandbox – Practical Cache Attacks in Javascript 相关论文可在 https://github.com/wyvernno ...

  10. 有关placeholder在ie9中的一点折腾

    有关placeholder在ie9中的一点折腾. placeholder属性定义: placeholder 属性规定可描述输入字段预期值的简短的提示信息(比如:一个样本值或者预期格式的短描述). 问题 ...