两路正弦波信号通过ADC采集信号进行FFT计算，怎么利用dsp库的FFT程序计算啊

yuanzimi

想测量两个正弦波的相位差，目前正弦信号已用ADC+DMAl两路通道采集到了正弦信号，接下来怎么利用怎么利用dsp库的FFT程序计算两个正弦波的幅值，实部，虚部，相位差啊，有哪位大神有没有FFT计算相关的例子程序啊，非常感谢！

WZTENG

做过的东西：
g_DSO为结构体，下面为它的部分成员,_usFFTPoints为点数(64/256/1024) 
uint32_t FFTIN [1024];
uint32_t FFTOUT[1024];
uint16_t FFTMAG[512];
uint32_t FFTPHASE[512];

cr4_fft_1024_stm32(g_DSO.FFTOUT,g_DSO.FFTIN,_usFFTPoints);


/******************************************************************************************
* 每个单元数据高字(高16位)中存储采样数据的实部，低字(低16位)存储采样数据的虚部(总是为0)
* 进行FFT变换，并计算各次谐波相位
* The phase value is in the range [-pi: pi]. 
*******************************************************************************************/
void ower_Phase_Radians(uint16_t _usFFTPoints)
{
    int16_t lX,lY;
    uint16_t i;

    for (i=0; i <(_usFFTPoints >> 1); i++)
    {
        lX= (g_DSO.FFTOUT<<16)>>16; /* 虚部 取低16bit，sine_cosine --> cos */
        lY= (g_DSO.FFTOUT >> 16);   /* 实部 取高16bit，sine_cosine --> sin */    
        {
            float X = (float)lX;
            float Y = (float)lY;
            float phase = atan2(X,Y);
             if (X>0)
             {
//                if (X>=0);
//                  //phase+=PI2;
//                else
//                  phase+=PI;  
             }
             else
             {             
//                if (Y>=0)
//                  phase+=PI2;                  
//                else
// { 
//   phase = -phase;
                  phase+=PI2; 
// }                   
             }                   
            g_DSO.FFTPHASE = phase* 180.0/PI;
        }    
    }
}

/******************************************************************************************
* 每个单元数据高字(高16位)中存储采样数据的实部，低字(低16位)存储采样数据的虚部(总是为0)
* 进行FFT变换，并计算各次谐波幅值
*******************************************************************************************/
void owerMag(uint16_t _usFFTPoints)
{
int16_t lX,lY;
uint16_t i;

if(_usFFTPoints == 1024)
{
_usFFTPoints = 512; 
}
else if(_usFFTPoints == 256)
{
_usFFTPoints = 256;
}
else
{
_usFFTPoints = 64;
}

/* 计算幅值 */
for (i=0; i < _usFFTPoints; i++)
{
lX= (g_DSO.FFTOUT<<16)>>16;               /* sine_cosine --> cos */
lY= (g_DSO.FFTOUT>> 16);                  /* sine_cosine --> sin */    
{
//     float X=  1024*((float)lX)/32768;
//     float Y = 1024*((float)lY)/32768;
//     float Mag = sqrt(X*X+ Y*Y)/1024;  
float X= (float)lX/32;
float Y = (float)lY/32;
float Mag = sqrt(X*X+ Y*Y)/1024;
g_DSO.FFTMAG= (uint32_t)(Mag*65536);
}    
}
}

WZTENG

我也是从这个文档学到一些东西，自己往后找相关的内容，介绍的比较详细。(附件)

正点原子

帮顶。。。。

yuanzimi

回复【2楼】正点原子:
---------------------------------
谢谢原子哥！回复【2楼】正点原子:
---------------------------------
回复【2楼】正点原子:
---------------------------------

yuanzimi

回复【3楼】WZTENG:
---------------------------------
请问你有这个FFT计算有一整套的例程吗，可以发给我看看吗，不慎感激

yuanzimi

回复【2楼】正点原子:
---------------------------------
谢谢！

luanshifusheng

回复【7楼】WZTENG:
---------------------------
同样的问题，有没有完整点的程序参考，谢谢

WZTENG

回复【8楼】luanshifusheng:
---------------------------------
这么好的资料还想怎么详细？
cr4_fft_1024_stm32(g_DSO.FFTOUT,g_DSO.FFTIN,_usFFTPoints);这个函数把参数传进去就行了。然后用上面的两个函数就可以得到幅值，相位，相位差自己把它们相减就行了，，，幅值是用勾股定理，相位用atan2();再不会就没有办法了。文档后面有专门的FFT函数介绍，慢慢悟。

yuanzimi

---
回复【9楼】WZTENG:
---------------------------------
非常感谢！