2

我正在使用 Ooura FFT 计算 1024 个样本窗口中加速度计数据的 FFT。该代码工作正常,但由于某种原因,它产生了非常奇怪的输出,即幅度为 10^200 的连续频谱。

这是代码:

 OouraFFT *myFFT=[[OouraFFT alloc] initForSignalsOfLength:1024 NumWindows:10]; //       had to allocate it

 UIAcceleration *tempAccel = nil;


 double *input=(double *)malloc(1024 * sizeof(double));
 double *frequency=(double *)malloc(1024*sizeof(double));

 if (input)

 {

 //NSLog(@"%d",[array count]);
 for (int u=0; u<[array count]; u++)
 {
  tempAccel = (UIAcceleration *)[array objectAtIndex:u];
  input[u]=tempAccel.z;
  //NSLog(@"%g",input[u]);
 }

 }

 myFFT.inputData=input; // specifies input data to myFFT


 [myFFT calculateWelchPeriodogramWithNewSignalSegment]; // calculates FFT


 for (int i=0;i<myFFT.dataLength;i++) // loop to copy output of myFFT, length of spectrumData is half of input data, so copy twice
 {

  if (i<myFFT.numFrequencies)
  {
   frequency[i]=myFFT.spectrumData[i]; // 
  }
  else 

  {
   frequency[i]=myFFT.spectrumData[myFFT.dataLength-i]; // copy twice
  }

 }





 for (int i=0;i<[array count];i++)

 {
  TransformedAcceleration *NewAcceleration=[[TransformedAcceleration alloc]init];  
  tempAccel=(UIAcceleration*)[array objectAtIndex:i];

  NewAcceleration.timestamp=tempAccel.timestamp;
  NewAcceleration.x=tempAccel.x;
  NewAcceleration.y=tempAccel.z;
  NewAcceleration.z=frequency[i];
  [newcurrentarray addObject:NewAcceleration]; // this does not work

  //[self replaceAcceleration:NewAcceleration];
  //[NewAcceleration release];
  [NewAcceleration release];
 }

 TransformedAcceleration *a=nil;//[[TransformedAcceleration alloc]init]; // object containing fft of x,y,z accelerations


 for(int i=0; i<[newcurrentarray count]; i++)
 {
  a=(TransformedAcceleration *)[newcurrentarray objectAtIndex:i];
  //NSLog(@"%d,%@",i,[a printAcceleration]);
  fprintf(fp,[[a printAcceleration] UTF8String]);  //this is going wrong somewhow
 }

 fclose(fp);

 [array release];
 [myFFT release];
 //[array removeAllObjects];

 [newcurrentarray release];

 free(input);
 free(frequency);
4

2 回答 2

1

您需要调用输入和输出数组。这是我的代码中的一个错误,无法修复,因为该库已被 Apple 的新 Accelerate 框架淘汰。

于 2010-05-03T16:32:33.043 回答
0

嘿,你试过绘制加速度计的值吗?如果有任何突然变化,FFT 将超出范围。

您是否尝试过绘制 FFT?你能张贴看起来像什么的图片吗?如果它在大多数情况下看起来正常,而不是偶尔超出范围,我打赌你应该事先平滑你的加速度计数据。

另外,为什么要复制spectrumData 两次?该数组已经设置为为您提供可绘制的数据。

最后,您是否重叠了您提供给 OouraFFT 的数据段?它的设计目的是提供大约 50% 重叠的数据段。否则你可能会得到奇怪的边缘效应。

于 2010-03-21T19:52:43.877 回答