STM32的ADC配合DMA采样以及STM32利用DSP库实现快速傅里叶变换,而当二者真正的结合到一起实现一个信号采样以及频域分析才可以发挥出很强大的功能。
' y. o2 P0 \% l' p: P' }) I+ z( Z7 H/ m9 m% f8 ^& l
本期我们就来演示,利用DMA+ADC+FFT实现信号采集与频域分析。
% e8 L b" X- M" [1 M( v4 D( }' b9 ^2 s; N5 `2 H/ m2 f5 B! F( }
CUBEMX配置
! O6 b S0 Q8 H; {# z2 ~! X `1 `/ {8 l j; w
0 d- }& h0 [1 F1 T; G
1 i \% P4 y5 Y. uADC配置这边除了之前的配置之外有一点区别,就是我们把循环采样给关掉,因为我们打算利用定时器来触发采样,这样子的话可以精准的控制采样率。
! \7 |3 ^- n7 A' E, ?7 W1 I# d0 u( P7 R n( ~
+ r9 T( e" h4 h' e% v1 p! `
- F/ \" E0 P6 V% D2 R这里我们选用TIM1的输入捕获来触发,当TIM1输出PWM波的时候,会选取上升沿来触发ADC采样,将采样的数据存储到DMA中。
6 K( j- Y" w) a: R- }
+ {$ [, Q6 D7 q' `1 h
- d9 g& F$ a6 F3 e0 s
4 w- e4 F! p" j+ Q" D' N. R
定时器中我们设置PWM的频率为500HZ,这代表着我们的采样率为500HZ,之后生成我们的代码。(这里非常奇怪,我明明觉得这么算好像是1000HZ但实际只有500HZ,我哪里算错了?)5 V# ^. _2 b. \ T
`2 j, y. h9 t代码编写
) E& ~( l- W' n+ B- #define FFT_LENGTH 1024//采样长度
) _- E4 X6 h2 ~+ _ - arm_cfft_radix4_instance_f32 scfft;//定义scfft结构体
O. B2 p* c6 H( [* ^/ [ - float FFT_InputBuf[FFT_LENGTH*2]; //FFT输入数组+ S$ P/ @, m: J7 K% V
- float FFT_OutputBuf[FFT_LENGTH]; //FFT输出数组5 V9 l. r+ _2 T5 N
- uint16_t AD_Value[1024] = {0};//存放ADC的值
复制代码
+ D: h9 {( r: l+ z2 C5 k上述代码需要定义在全局,否则可能会因为栈溢出(因为临时变量的定义是存放在栈中的)导致代码崩溃。0 E+ C, A2 O7 S" q% f( ^9 U
- HAL_TIM_PWM_Start(&htim1, TIM_CHANNEL_1); //用来触发adc采样 5 X; l3 r4 z& U# r
- HAL_ADC_Start_DMA(&hadc1, (uint32_t *)AD_Value, FFT_LENGTH);//开启ADC
) c( Y% E' J- Z - arm_cfft_radix4_init_f32(&scfft, FFT_LENGTH,0,1);
S( w. h7 {1 B! N - HAL_Delay(1000);//稍微等待一会ADC一轮转换结束
1 \1 l, X/ y$ o - for(int i=0; i < FFT_LENGTH; i++)2 i7 {9 g |: ?& F: z9 \2 Y
- {
# U' w# F% W3 Y; i g+ k+ r" O - FFT_InputBuf[2*i]=AD_Value[i]*3.3/4096; //实部/ r: {. j, A; J- I6 r0 o. \6 h
- FFT_InputBuf[2*i+1]=0; //虚部
6 G9 [$ y6 S3 n1 W+ I8 f0 ^/ ? - }4 D% g1 E9 K: B# S1 Y. |- i! d
- arm_cfft_radix4_f32(&scfft,FFT_InputBuf); 8 Y0 y6 L+ u' c2 a% N
- //arm_cmplx_mag_f32(FFT_InputBuf,FFT_OutputBuf,FFT_LENGTH); //取模得幅值& t' e' ~# j& O" X3 s3 `4 l1 X3 o
- for(int i = 0;i<FFT_LENGTH;i++)
" C2 K. `. Z2 d# ]/ p - {0 k: l; F# [$ h) d
- float32_t real = FFT_InputBuf[2 * i];
8 P/ E) d$ M6 L r6 i/ Y - float32_t imag = FFT_InputBuf[2 * i + 1];
+ }7 G, y. I( ?7 O0 j - float32_t magnitude = sqrtf(real * real + imag * imag);
; e4 Y: ~0 W) ` -
# U# t Z+ E( e" T. Z, V. D - // 打印每个频率分量的模值# X1 K$ e( P3 X6 U; p7 b3 h; [$ r
- printf("Fre: %f \r\n", magnitude);
0 y# a. p4 K* I' K7 r3 A% s - }
复制代码 r' U! t2 x. a7 f3 `5 A1 c
生成代码,我们由定时器触发采样,原始信号来自于一个250HZ的方波。
' u: f6 I9 I" u0 u: `' C* W. `! D! j" w- u3 v' d
k. a$ l( {3 n# a' T2 n+ y# k
3 d+ R. Z. e# V9 j0 i之后我们执行我们的FFT代码,首先我们先看一下采集的波形。6 R& D. C2 o& `) Q, L) m- P4 d5 g3 T
2 }. M) R0 h# n
) c! B+ j# h! I5 k
9 ^0 G6 @9 H$ r% k这里可以看到外面采集的方波。# ]* H7 _/ |5 b2 R! y
- W* `8 R0 X$ V: r: U3 w B5 A
注意事项8 [) Q' O. ~* r1 v5 |
(这里我的问题,根据奈奎斯特采样定理,采样频率应该至少是信号频率的两倍,否则会出现采样失真,因此为了能够顺利的采样,我将PWM的频率修改为100HZ)
( ?; d; O ]$ e1 Z7 ]* E/ p" Z0 o; L% r: A
* @- e2 c8 `% R$ X( T2 p4 L
+ l: x2 V n7 K+ ^这里也非常奇怪,我觉得PWM的频率应该是100HZ的,但是不知道为什么测出来是50HZ(难不成我又算错了?)% n, l+ p# r f7 x7 u
# Z+ N% Z+ X6 K3 c) j8 _
后来发现是我的分时系数算错啦~
2 V$ r! T: [1 u' c
% Y4 T, {& n; f变换结果6 O3 |' }3 {: m0 v/ `
然后来看看傅里叶变换之后的结果。0 Y* U5 }. @' r7 k' a: X' a4 d
, E- Y8 \0 A6 P# O }4 V6 Q
% s1 K6 _& t# N1 ~5 u* O- R, Y% m: y4 H t
(这个软件是我自己写的串口示波器)) G2 ^& D: `8 E' M; {4 T
& q' |: c+ N) i( R7 I每个FFT点对应的频率(f)与采样频率(Fs)和FFT的点数(N)之间的关系可以用以下的数学公式表示:
, s' b' d& {, n$ tf = k × (Fs / N)9 ?4 P+ a, I. R7 S+ f/ _9 o. a# ~0 o2 r
2 m0 @8 Y5 f0 v# c5 F其中,k 是FFT结果的索引(从0到N-1),Fs是采样频率,N是FFT的点数(与采样点数相同)。
; p. \* j: R/ i* {
" f# r* ]9 _$ \5 h [! }9 [; ?这里外设定的采样长度是1024,采样频率是500HZ,因此每个点之间的频率差是500/1024(这里按理来说1000更加方便计算但是不知道为什么1000就会卡死,1024就不会,这个也等我研究一下)。
! c2 b/ `) c7 e* z5 j/ ^
0 H! U7 s8 b: b$ h. B1 S. _8 l8 ]然后用Vofa+来算一下频率(自己的串口示波器不能显示下标,真是失败,我这两天就加上去这个功能)
9 x8 N5 ^! w' X. C+ W6 D1 ~8 D9 ]1 H0 X5 |
9 Q8 T8 t5 X2 I) E8 E( K, T6 x
" {! Q; z: [$ z* K. q* P5 r; p首先是开始的这个大信号,毫无疑问就是是直流信号(低频)所带来的幅度。9 i- J7 U# E2 b- Z1 O; I/ w# o
3 }1 c- k" j: o9 X! p; X
这里我给数据加一个直流滤波,就可以把这个很大的直流量给滤掉了。
! g1 W0 |1 L# G9 b- for(int i = 0;i<FFT_LENGTH;i++)
Y2 x+ h5 y0 L5 @. p) ]1 L - {
6 D/ h2 ]$ S5 O5 | - sum+=AD_Value[i];
2 L3 G+ v! S8 q. a8 w - }/ T' U/ A" v+ ]$ @" V. B; |
- sum/=1024;
$ ]% p9 }& B* T0 m - for(int i=0; i < FFT_LENGTH; i++)
* F' h, F/ {5 u - {
. Q& G& F. i: J+ O! } - FFT_InputBuf[2*i]=(AD_Value[i]-sum)*3.3/4096; //实部7 m! r/ A8 ^3 q9 |3 p7 Z: r( l: o9 S
- FFT_InputBuf[2*i+1]=0; //虚部
& h* p7 R/ k) x - }
复制代码
. I& E; ^; ?6 j5 k从图中也可以看到,FFT的结果不是一个单边频谱,而是一个双边频谱。- C! F2 _4 }& X$ i
4 S: z; }5 ?' q' f. c9 dFFT(快速傅里叶变换)的结果默认是双边频谱,这主要是因为FFT是一种处理复数信号的方法,而复数信号在频域中既有正频率分量也有负频率分量。即使你输入的是实数信号,FFT算法也会将其视为复数信号来处理,并返回包括正频率和负频率在内的双边频谱。% [# I$ g: {8 S; |) {+ p( Q
1 l+ o3 @8 L$ W9 O$ n+ k3 Y
1 f# i% |1 l" G) x' K' u$ r3 j
( G& l" l% b4 P0 g, W" {& q) v两个峰值之间相差100个点,也就是存在50HZ的倍数波。" i; R, e! j& O$ A
* s2 b5 w2 k4 `8 C8 ]2 q# H
/ N8 b) T. h9 R) y' g" Y5 B0 ~. }7 J5 z2 R8 a* K
这是因为方波是由它的基波频率以及高次谐波所组成。
% m- d* R1 x$ t r; w. w
# d2 @8 g. c; b
- ?) h4 {3 j7 I! } V! R* u' J$ z4 c3 _ t8 ^; f
& B9 ]9 p a' j
* P6 m A, y& a+ R7 ?* ] z0 Z' Y& n6 k- J; j, b; {% o
转载自:电路小白: [. o/ L. C2 B2 ?6 k
如有侵权请联系删除, w: ]0 u3 q: Y
2 a. H3 z6 x' E3 G$ _$ E% d1 r
c9 H3 w; |7 g$ T$ t D1 q
6 B: z$ A+ `( k/ n5 M5 `4 d |
学习学习,正好需要FFT频谱分析