[问答]

FFT如何在STM32中的实现与应用

问答对人有帮助，内容完整，我也想知道答案 0 FFT如何在STM32中的实现与应用 0
2021-10-13 08:58:26　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × heks 该类别下有 50 个回答。邀请回答 hgimtk 该类别下有 44 个回答。邀请回答新星之火12138 该类别下有 42 个回答。邀请回答 wang21cj 该类别下有 41 个回答。邀请回答 chm5 该类别下有 40 个回答。邀请回答 hjfjsdgfjdsf 该类别下有 37 个回答。邀请回答 werywer 该类别下有 35 个回答。邀请回答 h1654155957.9520 该类别下有 35 个回答。邀请回答 fdjslkjd 该类别下有 35 个回答。邀请回答 uvysdfydad 该类别下有 35 个回答。邀请回答 jenny042 该类别下有 34 个回答。邀请回答 dfzvzs 该类别下有 34 个回答。邀请回答维生素B2 该类别下有 34 个回答。邀请回答凤毛麟角该类别下有 34 个回答。邀请回答 ggfvxv 该类别下有 33 个回答。邀请回答刘洋该类别下有 33 个回答。邀请回答 meihuacg 该类别下有 33 个回答。邀请回答 hfgdzc 该类别下有 33 个回答。邀请回答 h1654155275.6678 该类别下有 32 个回答。邀请回答河南顺之航该类别下有 32 个回答。邀请回答举报马刺牛逼相关推荐 • 如何在stm32里面去实现FFT呢 1557 • STM32的DSP库可以实现FFT变换输出吗？ 337 • 如何使用STM32F103汇编库去实现FFT呢 1147 • labview中怎么实现FFT和逆FFT并播放波形 9558 • 怎么实现基于STM32-AD7606的FFT交流采样？ 1085 • 如何在CubeIDE中的STM32上通过USB实现“printf”发送消息？ 228 • 如何使用STM32提供的DSP库函数进行FFT 1738 • 如何使用STM32提供的DSP库进行FFT呢 1110 • 怎么使用cordic旋转方法在VHDL中实现FFT？ 1357 • 怎么实现基于LabVIEW的数字变频FFT设计？ 1814 1个回答

答案对人有帮助，有参考价值 0 0、前言音乐频谱显示说白了就是“儿童版”频谱仪。笔者平时比较喜欢听音乐，闲暇之余听音乐的时候如果有个频谱显示的小玩具在旁边跳来跳去的也挺有意思的，所以笔者去万能的某宝上搜索了一下，发现便宜的都很小，大一点的都很贵，而且都需要音频接头输入，很麻烦，所以笔者就自制了这个小玩具。效果图如下图1所示。效果视频:https://www.bilibili.com/video/BV1dZ4y1x77d 图1：实物展示其功能有一下3个： 1、实时显示音乐频谱 2、将小玩具倒立，显示当前时间和温度。 3、在显示时间和温度时，一段时间后屏幕自动熄灭。触动小玩具即可再次点亮 1、FFT（Fast Fourier transform）快速傅里叶变换是一种利用计算机进行高速快捷的离散傅里叶变换的一种方法，而什么是傅里叶变换？傅里叶变换的物理意义是什么呢？当年，著名的数学家傅里叶表示：任何连续测量的时序或信号，都可以表示为不同频率的正弦波信号的无限叠加。而通过傅里叶变换就可以将叠加的波形以正弦波为基本单位将其拆分开。图2：傅里叶变换而快速傅里叶变换为什么可以快速呢？首先先看一下离散傅里叶变换的公式：通过公式不难看出如果我们让计算机去计算时，我们需要计算N^2次运算。但是通过观察可以发现，在计算中有想当一般的数据都是冗余的，并且不会显示在频谱中。所以我们可以将这些数据保存下来计算的时候直接调用或在此数据基础上做简单的运算。通过这样的方法我们就可以将DFT的N^2次运算减少到Nlog10（N）次运算了。例如如果做1000采样点的傅里叶变换使用DFT的运算次数就是1000000次，而使用FFT就可以减少到3000次运算。 2、STM32实现FFT 因为笔者做的这个玩具使用的芯片是STM32单片机F10系列。所以在做FFT的时候没有官方库的支持（以前是有，但是F4系列上线后就取消了，真实万恶的资本主义dogs）虽然网上有很多教程是直接移植或者找之前的DSP库，但是笔者觉得着实有些麻烦，F10系列没有DSP运算核心，这个DSP库就犹如花瓶。这里主要借鉴了这位大神的文章：https://blog.csdn.net/yga_airspace/article/details/86688278 在单片机上做FFT时需要注意堆栈的大小，以免发生硬件错误。代码如下： double * FFT(double pr,int n) { static double amp_data[FFT_max2];//定义输出幅值 static double fr_data[FFT_max2];//定义输出实部 static double fi_data[FFT_max2];//定义输出虚部 static double pr_data[FFT_max];//缓存输入实部 static double pi_data[FFT_max];//缓存输入虚部 unsigned int k;//n=2^k; int i,j; int it,m,is,nv; double p,q,s,vr,vi,poddr,poddi; //===============初始化================================// k = log10(n)/log10(2);//n=2^k; for(i=0;i { pr_data = (pr+i);//读入缓存 pi_data = 0;//读入缓存 } //========将pr[0]和pi[0]循环赋值给fr_data[]和fi_data[]==========// for (it=0; it<=n-1; it++) { m=it; is=0; for(i=0; i<=k-1; i++) { j=m/2; is=2is+(m-2j); m=j; } fr_data[it]=pr_data[is]; fi_data[it]=pi_data[is]; } //==================================================// pr_data[0]=1.0; pi_data[0]=0.0; p=6.283185306/(1.0n); pr_data[1]=cos(p); //将w=e^-j2pi_data/n用欧拉公式表示 pi_data[1]=-sin(p); //================计算pr_data和pi_data========================// for (i=2; i<=n-1; i++) { p=pr_data[i-1]pr_data[1]; q=pi_data[i-1]pi_data[1]; s=(pr_data[i-1]+pi_data[i-1])(pr_data[1]+pi_data[1]); pr_data=p-q; pi_data=s-p-q; } //===================计算fr_data和fi_data=====================// for (it=0; it<=n-2; it=it+2) { vr=fr_data[it]; vi=fi_data[it]; fr_data[it]=vr+fr_data[it+1]; fi_data[it]=vi+fi_data[it+1]; fr_data[it+1]=vr-fr_data[it+1]; fi_data[it+1]=vi-fi_data[it+1]; } //=================蝴蝶操作=========================// m=n/2; nv=2; for (i=k-2; i>=0; i--) { m=m/2; nv=2nv; for (it=0; it<=(m-1)nv; it=it+nv) { for (j=0; j<=(nv/2)-1; j++) { p=pr_data[mj]fr_data[it+j+nv/2]; q=pi_data[mj]fi_data[it+j+nv/2]; s=pr_data[mj]+pi_data[mj]; s=s(fr_data[it+j+nv/2]+fi_data[it+j+nv/2]); poddr=p-q; poddi=s-p-q; fr_data[it+j+nv/2]=fr_data[it+j]-poddr; fi_data[it+j+nv/2]=fi_data[it+j]-poddi; fr_data[it+j]=fr_data[it+j]+poddr; fi_data[it+j]=fi_data[it+j]+poddi; } } } //=================计算幅值输出==========================// for (i=0; i<=n-1; i++) { amp_data=sqrt(fr_datafr_data+fi_datafi_data); //幅值计算 } //====================返回数据=========================// return amp_data;//返回幅值 } 直接通过KEIL仿真可以观察出FFT的输出波形与MATLAB对比图如下图3，图4 图3：Keil仿真波形图4：MATLAB波形方波的FFT测试如下图5所示：图5：方波测试 3、显示屏选用与点阵驱动笔者使用了6个88点阵组成的2416点阵作为显示屏。每个88点阵使用MAX7219驱动。MAX7219之间级联保证占用的I/O口最少。其电路原理图如下图6所示图6：MAX7219驱动电路因为MAX7219这款IC输入电压范围在4-5.5V故一般需要5V供电，信号端也需要5-3.3V的转换。虽然大部分STM32的I/O口都是可以容忍5V输入与输出的。但笔者为了保险起见加了3路的3.3V-5V电路。其电路原理图如下。 4、STM32与音频输入电路笔者将整个小玩具分成了两个块电路板，其中一块全是点阵和MAX7219一块是STM32、锂电池充放电电路、MPU6050外围电路和一路的音频采集电路，引出剩下部分引脚。音频采集电路主要使用咪头，经MAX4466做前置放大，放大后信号直接接入STM32内置ADC引脚。电路图如下图7所示图7：MAX4466驱动电路 5、点阵屏时间显示时间的计算采用STM32内置的RTC时钟即可，首先需要显示数字的字模，根据电路和字体大小（53），计算出数字0-9字模如下： unsigned char disp1[11][3]={ {0x1F,0x11,0x1F},//0 {0x00,0x1F,0x08},//1 {0x1D,0x15,0x17},//2 {0x1F,0x15,0x15},//3 {0x1F,0x04,0x1C},//4 {0x17,0x15,0x1D},//5 {0x17,0x15,0x1F},//6 {0x18,0x17,0x10},//7 {0x1F,0x15,0x1F},//8 {0x1F,0x15,0x1D},//9 {0x00,0x0A,0x00} //: }; 根据电路结构和扫描顺序，写出底层显示驱动程序，程序如下： //参数从左到右依次为：字体数据，字体长度，字体宽度，起始位置的x坐标（自右向左），起始位置y坐标（自下而上） void MAX7219_N_display(unsigned char cData,unsigned char Length,unsigned char Width,unsigned char x,unsigned char y) { static u16 iSbufWidth[24];//纵向显示缓存 u16 iSbuf;//用户数据缓存 u16 iSbuf_Middle;//缓存中间变量 unsigned int i; for(i=0;i { iSbuf = 0;//缓存清零 iSbuf \|= (cData+i)< iSbufWidth[i+x] \|= iSbuf;//向显示缓存写入有效位的“1” iSbuf = ~*(cData+i);//按位取反 iSbuf <<= y;//左移到显示位置 iSbuf = ~iSbuf;//相当于左移补1 iSbuf_Middle = 0xffff;//置位 iSbuf_Middle <<=(y+Width); iSbuf \|= iSbuf_Middle;//补全高位使得只有有效位存在"0"的可能，非有效位全为1 iSbufWidth[i+x] &= iSbuf;//向缓存写入有效位的“0” } for(i=0;i<8;i++) { Write_Max7219_N(i+1,iSbufWidth,iSbufWidth[i+8],iSbufWidth[i+16],iSbufWidth[i+16]>>8,iSbufWidth[i+8]>>8,iSbufWidth>>8);//显示 } } 6、总结这个装置整体来说比较简单，因为有姿态传感器和一个点阵屏，所以开发的自由度也比较高。本文附资源，资源包含两个电路板的PCB源文件和开发的源代码。如有好玩的想法或者改进欢迎留言讨论。

2021-10-13 16:25:17 评论举报钱杨静