starcat2002 2019-09-16 17:25 采纳率: 0%
浏览 792

求Dalao告诉我FFT如何求幅值相位等等

我用
#define PI 3.1415926

#define X2D (PI / 180.0)
#define A 97.14
#define Q (12.37 * X2D)
#define DD (360.0 / 1024.0)

while (1024 > iCounter)
{
    in[iCounter][0] = (double)(A * sin((((double)iCounter) * DD * X2D) + Q));
    in[iCounter][1] = 0.0;
    //pData[iCounter] = (double)(A * sin((((double)iCounter) * DD * X2D) + Q));
    iCounter++;
}

创建一个波形,然后用FFTW
    fftw_plan plan;

plan = fftw_plan_dft_1d(1024, in, out, FFTW_FORWARD, FFTW_ESTIMATE);

fftw_execute(plan);
进行FFT

之后
for (iCounter = 0; iCounter < 200; iCounter++)
{
    double dT = sqrt(out[iCounter][0] * out[iCounter][0] + out[iCounter][1] * out[iCounter][1]);
    printf("%lf\n", dT);
    if (dT > AM)
    {
        AM = dT;
    }
}

得到的时49735这样的值。请问我是不是哪里错了?
另怎么从FFT结果中求得A,Q等值?
  • 写回答

1条回答 默认 最新

  • zqbnqsdsmd 2019-09-17 13:01
    关注
    评论

报告相同问题?

悬赏问题

  • ¥15 visual studio2022中文乱码无法解决
  • ¥15 关于华为5g模块mh5000-31接线问题
  • ¥15 keil L6007U报错
  • ¥15 webapi 发布到iis后无法访问
  • ¥15 初学者如何快速上手学习stm32?
  • ¥15 如何自动更换布娃娃图片上的衣服
  • ¥15 心理学eprime编程
  • ¥15 arduino esp8266开发
  • ¥15 stm32单片机通过485发送命令给驱动器控制电机转动,同样的代码f103可以控制电机转动,换到f407不能动了,但是用串口助手调试407显示发送的命令都是正确的,卡了好久了这是发送规则
  • ¥15 stm32f103c8t6最小系统板+2.8寸TFTLCD板子