利用DSP库做FFT运算

qingyu111

使用stm32利用DSP库做FFT运算调用cr4_fft_256_stm32(lBUFOUT, lBUFIN, NPT)计算各点幅值，

	
[mw_shl_code=c,true]for (i=0; i < nfill; i++)
  {
    lX= (lBUFOUT<<16)>>16; 
    lY= (lBUFOUT >> 16);       
    {
      float X=  64*((float)lX)/32768;
      float Y = 64*((float)lY)/32768;
      float Mag = sqrt(X*X+ Y*Y)/nfill;  // 先平方和,再开方
      lBUFMAG = (u32)(Mag*65536);  //
    }    
  }
}
[/mw_shl_code]
[mw_shl_code=c,true]64*((float)lX)/32768中64是什么？为什么要除以32768？能详细分析一下吗？[/mw_shl_code]
[mw_shl_code=c,true]

[mw_shl_code=c,true]另外lBUFMAG = (u32)(Mag*65536)中为什么要乘以65536？谢谢！[/mw_shl_code]


[/mw_shl_code]



[/mw_shl_code]

unciauncia

float X    = NPT * ((float)lX) /32768;   
   float Y    = NPT * ((float)lY) /32768
目的就是把数据浮点化，至于为什么是除以32768 。可以这么说，浮点化 就好像10进制里面的科学计数法。。32768=2的15次。除以32768也就是去除了浮点数后面的那个基数，只剩下前面的。比如1991 改写成1.991*10的三次幂，再除以10的三次方，只剩下1.991，便于余下的运算。
之后的语句就不用解释了吧至于最后一句要乘以65536是因为，我们定义的数据和我们需要求得的数据都是无符号32位的，之前已经把32位的数据拆开又分别浮点化了又开了个根号，所以再把它变回来 只需要乘以2的16次，也就是65536.
（我也是看完别人的讲解才懂的http://zhidao.baidu.com/link?url=wBCyVbXpHTMc6UvoEwjb_0SE4aAO7FQWd7CqqhWIOq96_eeyV9WgKL6ggk3tmx8TCCIJ6uozjabJIWFFZAkprv3isLK2lWEsPlxBNjiuyzK）

正点原子

帮顶。。。。