用DSP做AD采样然后做FFT变换

只看该作者 · 2015-5-26 14:39

出现了一个错误，给学长看了半天也没有解决

只看该作者 · 2015-5-26 21:41

你换个学长看看，采样后做FFT，你要先AD做出来啊，然后得到的是数字信号，你需要离散数字FFT。

2万 · 2015-5-26 21:43

ADC采集完成了是吧然后找个FFT的函数先做处理吧
结合理论在慢慢理解

只看该作者 · 2015-5-27 08:51

zhangmangui 发表于 2015-5-26 21:43
ADC采集完成了是吧然后找个FFT的函数先做处理吧
结合理论在慢慢理解 ...

AD是实验室以前学长们做的，拿过来直接用，但是在AD工程主函数里加入了我自己编的FFT，采不到数据，另一个同学也用的这个AD就采到了数据。昨天求助时候的问题已经解决了，也就是说，现在AD和FFT部分的程序都有，就是结合起来还是由问题

只看该作者 · 2015-5-27 08:51

玛尼玛尼哄发表于 2015-5-26 21:41
你换个学长看看，采样后做FFT，你要先AD做出来啊，然后得到的是数字信号，你需要离散数字FFT。 ...

AD已经做出来了，但是采不到数据

2万 · 2015-5-27 22:31

z641230919 发表于 2015-5-27 08:51
AD是实验室以前学长们做的，拿过来直接用，但是在AD工程主函数里加入了我自己编的FFT，采不到数据，另一 ...

看看你的ADC怎么触发采集的 FFT部分是不是进入了循环
ADC不再触发

只看该作者 · 2015-5-30 09:11

z641230919 发表于 2015-5-27 08:51
AD已经做出来了，但是采不到数据

都做出来了，还采集不到数据，莫非没有数据吗？

1万 · 2015-5-30 11:03

出来的ADC不是数据吗？你的意思是数据不能FFT，是不是FFT写的部队

只看该作者 · 2015-8-22 11:16

ADC的值怎么处理可以直接带入TI的fft lib，看资料有的是Q30 ,有的是Q31.可ccs中也没有Q31的格式。

1万 · 2015-8-23 12:12

那就一次性采集512组数据或者1024组数据的做下FFT变换的呗。

只看该作者 · 2015-8-25 12:56

我给共享个FFT算法的吧
/*¿ìËÙ¸µÀïÒ¶±ä»»*/
void fft(void)
{
int i=0,j=0,k=0,l=0;
    complex up,down,product;
    change();
    for(i=0;i<log(size_x)/log(2);i++)/*Ò»¼¶µûÐÎÔËËã*/
    {
    l=1<<i;
    for(j=0;j<size_x;j+=2*l)/*Ò»×éµûÐÎÔËËã*/
    {
    for(k=0;k<l;k++)/*Ò»¸öµûÐÎÔËËã*/
    {
mul(x[j+k+l],W[size_x*k/2/l],&product);
add(x[j+k],product,&up);
sub(x[j+k],product,&down);
x[j+k]=up;
x[j+k+l]=down;
      }
   }
   }
}

只看该作者 · 2015-8-25 12:57

这个是快速傅里叶逆变换的
/*¿ìËÙ¸µÀïÒ¶Äæ±ä»»*/
  void ifft()
  {
  int i=0,j=0,k=0,l=size_x;
  complex up,down;
  for(i=0;i< (int)( log(size_x)/log(2) );i++) /*Ò»¼¶µûÐÎÔËËã*/
  {
  l/=2;
  for(j=0;j<size_x;j+= 2*l ) /*Ò»×éµûÐÎÔËËã*/
  {
  for(k=0;k<l;k++) /*Ò»¸öµûÐÎÔËËã*/
  {
  add(x[j+k],x[j+k+l],&up);
  up.real/=2;up.img/=2;
  sub(x[j+k],x[j+k+l],&down);
  down.real/=2;down.img/=2;
  divi(down,W[size_x*k/2/l],&down);
  x[j+k]=up;
  x[j+k+l]=down;
  }
  }
  }
  change();
  }