为什么F407做FFT之后，求相位，相位值一直在变

yimengxiaoyao

利用内部adc采集交流信号，然后做用dsp的库做fft之后求相位，相位值一直在变是什么问题

Samouyang

因为你采样的参考点每次不一样的，而我们通常说的相位是对0点的值，所以是会不一样的。

yimengxiaoyao

这个是求相位部分的代码

1. #define pi 3.1415926f
   
2. #define FFT_LENGTH 4096
   
3. 
   
4. float fft_inputbuf[FFT_LENGTH*2];        //FFT输入数组 计算采样电阻前的电压
   
5. float fft_outputbuf[FFT_LENGTH];        //FFT输出数组
   
6. 
   
7. 
   
8. float FFT_ADC(u8 ch,u8 mode)
   
9. {
   
10.                 float ADC_temp;
    
11.                 float max;
    
12.                 int max2;
    
13.                 float v_ac;
    
14.                 float v_dc;
    
15.                 float phase;//相位
    
16.                 double real,im;
    
17.                 float hudu;
    
18.                 u16 i;
    
19. 
    
20.                         while(!adc_flag)
    
21.                         {
    
22.                         }
    
23.                                 TIM_Cmd(TIM3,ENABLE);
    
24.                                 adc_flag=0;
    
25.                           max=0;
    
26.                                 for(i=0;i<FFT_LENGTH;i++)
    
27.                                 {
    
28.                                         ADC_temp=(float)ADC1data[ch][i]*(3.3f/4095.0f);
    
29.                                         fft_inputbuf[2*i]=ADC_temp;//高位放置实数  1010 0000
    
30.                                         fft_inputbuf[2*i+1]=0;//低位放置虚数       
    
31.                                 }
    
32.                                 arm_cfft_f32(&arm_cfft_sR_f32_len4096,fft_inputbuf,0,1);
    
33.                                 arm_cmplx_mag_f32(fft_inputbuf,fft_outputbuf,FFT_LENGTH);          //FFT计算（基4）               
    
34. 
    
35.                                 for(i=0;i<(FFT_LENGTH/2+1);i++)
    
36.                                 {  
    
37.                                         if(max < fft_outputbuf[i+2])       
    
38.                                         {
    
39.                                                  max = fft_outputbuf[i+2];        //找fftout的最大值
    
40.                                                  max2=i+2;
    
41.                                         }
    
42.                                 }
    
43.                                 if(mode==1)
    
44.                                 {                                                               
    
45.                                
    
46.                                         v_ac = (max/FFT_LENGTH)*2.308834f;
    
47.                                         return v_ac;
    
48.                                 }
    
49.                          else if(mode==2)
    
50.                                 {
    
51.                                         v_dc =  (fft_outputbuf[0]/FFT_LENGTH);
    
52.                                         return v_dc;
    
53.                                 }       
    
54.                         else
    
55.                                 {
    
56.                                         real=(double)fft_inputbuf[2*max2];
    
57.                                         im=(double)fft_inputbuf[2*max2+1];
    
58.                                         hudu=atan2(im,real);
    
59.                                         phase=hudu*180.0f/pi;
    
60.                                         if(phase<0)
    
61.                                                 phase+=360.0f;
    
62.                         //                phase= atan2((double)fft_inputbuf[2*max2+1],(double)fft_inputbuf[2*max2])*180.0/pi;
    
63.                                         return phase;
    
64.                                 }                                                       
    
65. }

复制代码

yimengxiaoyao

这里是在主函数里做的显示

1.                 while(1)
   
2.                 {
   
3.                         //ctp_test();
   
4.                         jiaodu=FFT_ADC(1,3);
   
5.                         sprintf((char*)dd,"jiaodu=%f",jiaodu);
   
6.                         LCD_ShowString(400,200,200,24,24,dd);       
   
7.                         delay_ms(1000);
   
8.                 }       

复制代码

edmund1234

相位变是正常的

正点原子

帮顶

yimengxiaoyao

edmund1234 发表于 2019-12-25 20:55
相位变是正常的

我想检测一个信号经过一个放大电路的相位差，是不是两个两个相位相减的值是不变的

zbtc啊

edmund1234 发表于 2019-12-25 20:55
相位变是正常的

你确定吗

edmund1234

zbtc啊 发表于 2020-1-4 19:13
你确定吗

这个世界不存两个完全一样的信号， 频率些许的差别， 就自然产生相位差， 这有什么需要怀疑的啊

科安达AA

楼主，最近再学习这个，能发份源码给我参考下吗，谢谢 391588997@qq.com

h418452224

楼主问题解决了么？我的也是一直在变

xiaocaijii

楼主问题解决了吗，我的也是一直在变，不断递减

云边有个你

yimengxiaoyao 发表于 2019-12-27 21:42
我想检测一个信号经过一个放大电路的相位差，是不是两个两个相位相减的值是不变的

你好，我最近也需要求两个相位差，请问你解决了吗，416341641，可以交流下吗