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

求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 素材场景中光线烘焙后灯光失效
  • ¥15 请教一下各位,为什么我这个没有实现模拟点击
  • ¥15 执行 virtuoso 命令后,界面没有,cadence 启动不起来
  • ¥50 comfyui下连接animatediff节点生成视频质量非常差的原因
  • ¥20 有关区间dp的问题求解
  • ¥15 多电路系统共用电源的串扰问题
  • ¥15 slam rangenet++配置
  • ¥15 有没有研究水声通信方面的帮我改俩matlab代码
  • ¥15 ubuntu子系统密码忘记
  • ¥15 保护模式-系统加载-段寄存器