stm32f103zet6 FFT测相位不准确

玉米斑

我用定时器2控制ADC采样，然后用dma转换，之后对得到的数据进行fft测相位，但是得到的结果非常不准确，并且大范围跳动着，这是为什么呢？
代码如下：

arm_cfft_f32(&arm_cfft_sR_f32_len1024, FFTin, 0, 1);
arm_cmplx_mag_f32(FFTin, FFTout, NPT);

//得到真实幅值
FFTout[0] /= 1024;
for ( i = 1; i < NPT; i++)//输出各次谐波幅值
{
FFTout /= 512;
}

//相位

for(i = 1;i<NPT/2;i++)

{   

    Phase = atan2(FFTin[i*2+1], FFTin[i*2+0]);

    Phase = Phase*(360/PI2);//求相位

}

玉米斑

顶一下

rebirth9

兄弟 我最近也在弄这个 但是我的采集就有问题 而且这个能够采集到80khz的交流信号嘛

rebirth9

可以给我发一下你的代码嘛

正点原子

继续努力

玉米斑

rebirth9 发表于 2022-4-25 23:58
兄弟 我最近也在弄这个 但是我的采集就有问题 而且这个能够采集到80khz的交流信号嘛

我这个采不到80khz，试试sinc插值看看

arm51avr

1、最好的办法是把原始数据通过串口传到PC，用MATLAB看下波形，写两行代码分析下。

1. close all
   
2. fp = fopen('D:\2021\matlab\D001.dat');
   
3. a01 = fread(fp,'int16');
   
4. 
   
5. b01 = a01(1:1024);
   
6. 
   
7. plot(b01);
   
8. 
   
9. m01 = b01-mean(b01);
   
10. win = hann(length(m01));
    
11. h01 = m01'.*win';
    
12. 
    
13. 
    
14. figure;
    
15. c = fft(h01);
    
16. c2 = abs(c(1:end));
    
17. plot(2*c2);

复制代码

2、如果数据采集没有问题的话，注意下信号源的的频率和幅度。

3、注意下arctan( )函数的返回范围，是否需要Phase unwrapping？