STM32官方DSP库进行FFT运算的问题

smileqq

使用STM32官方的DSP库进行FFT运算，然后求每一个频率点的幅值，并把幅值用printf在串口打印出来，问题是我把计算幅值的函数放到main.c文件中，结果是正确的，但把它放到其他文件中，在main函数中调用，结果就不正确了，哪位大神知道问什么？
这是main函数，作用就是计算FFT后每一点的幅值，其中GetPowerMag(FFT_Out_Array, FFT_Mag_Array);函数是在fft.c中定义的，但结果不对，如果拿到main.c中就正确了
9 {7 ^' G0 J. q" X$ y& y# F#include "config.h"
#include "led.h"
#include
, x* G7 y( v( x6 v' y5 C#include "fft.h"
& p& I+ V5 _, n1 |#define Fs 6400 //采样率
#define PI2 6.28318530717959
( r' o# W, s4 w1 ?long FFT_In_Array[NPT];
long FFT_Out_Array[NPT];
unsigned long FFT_Mag_Array[NPT/2];
void fft_Init(void);
+ L9 H8 I6 G2 k# Zint main(void)
{
, L, t( D5 B7 |. O) [3 {" C( `u8 i;
$ K" K4 T0 Q4 V# N5 k: {+ Q: HRCC_Configuration();
GPIO_Configuration();
8 L9 V; d; ~6 W1 D$ `NVIC_Configuration();
USART1_Configuration();

//ADC_Configuration();
! _5 U. |7 K0 z1 w6 |1 y9 H
0 D! T8 R, h2 x; T8 ^& }
//DMA_Configuration();
" K) `5 D! R; \5 [! W; \: z7 `+ X// initLED();
9 Q" j5 q+ N# J( D2 [1 {fft_Init();
1 `6 e" P9 z& J* {/ o0 N
" p# X# X! @' Y8 q, W* Z; e' f* v//FFT运算
+ {, |# r. M) Ucr4_fft_256_stm32(FFT_Out_Array, FFT_In_Array, NPT);
; s0 B8 t2 i  V/ H$ s& ?1 _
. @2 u& s, {& `  J) G//计算幅值
$ w+ Z) ^$ d; B+ |0 n0 v$ l& G+ ZGetPowerMag(FFT_Out_Array, FFT_Mag_Array);

//打印每一点的幅值
for(i=0;i<NPT/2;i++)
{      
& V# P5 x* Y' m1 W    printf("%d=%ld\n", i, FFT_Mag_Array);
. `$ I3 {, l- f4 o8 k}

while (1)
1 A' t1 y+ l; q% t7 a% d: J{                           
% l3 {8 Z# ?4 |- s5 r- {
& o0 Z6 \6 y4 \+ l) @0 R  T: c
}

return 0;
}
//用指定频率的波形，初始化FFT数组
void fft_Init()
4 `, \+ z6 e1 K7 X{
: g% n' N, ^6 m3 }unsigned short i;
3 b' w9 \: r& jfloat fx;
for(i=0; i<NPT; i++)
3 e- X1 h" W8 j- I{
fx = 1500 * sin(PI2 * i * 350.0 / Fs) +
2700 * sin(PI2 * i * 8400.0 / Fs) +
  s- c, {6 L$ |7 P* A; w4000 * sin(PI2 * i * 18725.0 / Fs);
2 r# v9 @/ P$ B; V- j/ _: wFFT_In_Array = ((signed short)fx) << 16;
}
* D$ I# C$ y' M: u. k( \2 \}
以下是fft.c---------------------------------------
#include "fft.h"
; z/ {# x& x' i. m1 l) ivoid GetPowerMag(long *pIn, unsigned long *pMag)
8 e- n  [: o4 j/ V/ ?3 I! s  \{
! B9 T; o! Z9 b- }8 u  b  Ksigned short lX,lY;
) I) Y0 a1 K: ~+ b7 O2 efloat X,Y,Mag;
& O% X$ `+ T+ qunsigned short i;
) r7 k& c3 |( o6 {& ffor(i=0; i {
lX = (pIn > 16;
" j. L9 r# C3 {, a# JlY = (pIn >> 16);
X = NPT * ((float)lX) / 32768;
  l0 p# B: \9 ^1 |Y = NPT * ((float)lY) / 32768;
# t( s# K" c" B" kMag = sqrt(X * X + Y * Y) / NPT;
if(i == 0)
# x6 H7 A& |. C/ V5 wpMag = (unsigned long)(Mag * 32768);
& ]  Y2 h2 V9 i) ]2 K! Helse
: D' ^0 r9 M' n5 B' @pMag = (unsigned long)(Mag * 65536);
}
# Z. j3 b; E  ~2 Y( z# g" O) N9 q, ^}

yxlk

从源代码的管理上来说，应该不存在你所描述的问题。
像你这样出现的问题，可能是变量或其它未同步引起的。确认一下两种传入参数与各全局变量是否相同吧！