【C言語】振幅スペクトルのグラフを出したい

Question

前提・実現したいこと

上図のような振幅スペクトルのグラフを指数関数(e^iみたいな)で作りたいです。
横軸が周波数、縦軸が振幅(の絶対値)です。
三角関数で式を表記した時はうまくいったのですが、指数関数で表記すると縦軸である振幅の絶対値(Fx[1])が0にしかならず、困っています。
式以外はすべて同じソースコードです。

うまくいったソースコード

#include <stdio.h>
#include <math.h>
#include <complex.h>

#define aM (10.0)    //時刻tにおける加速度a_max[m/s/s]
#define T  (0.0332)  //終端時刻[s]●

int main(void)
{
  double w = 0.0, dw = 1.0;  //角周波数
  
  do {
    double _Complex Fx[2];
    double f;

    f = w/(2.0*M_PI);  //周波数
    Fx[0] = aM*(cos(w*T) - I*sin(w*T) - 1.0)/(-I*w);  //加速度のフーリエ変換
    Fx[1] = cabs(Fx[0]);  //Fx[0]の絶対値(振幅の絶対値)
    
    printf("%g %g+%gi\n", f, creal(Fx[1]), cimag(Fx[1]));
       
    w += dw;
  } while (w <= 2000.0);
}

うまくいかなかったソースコード

#include <stdio.h>
#include <math.h>
#include <complex.h>

#define aM (10.0)    //時刻tにおける加速度a_max[m/s/s]
#define T  (0.0332)  //終端時刻[s]●

int main(void)
{
  double w = 0.0, dw = 1.0;  //角周波数
  
  do {
    double _Complex Fx[2];
    double f;

    f = w/(2.0*M_PI);
    Fx[0] = aM*(exp(-I*w*T) - 1.0)/(-I*w);
    Fx[1] = cabs(Fx[0]);
    
    printf("%g %g+%gi\n", f, creal(Fx[1]), cimag(Fx[1]));
       
    w += dw;
  } while (w <= 2000.0);
}

Comments

[jimbe]

式以外が同じであれば、式の各要素を表示して確認してみては如何でしょうか。

[jbpb0]

> Fx[0] = aM*(exp(-I*w*T) - 1.0)/(-I*w);

複素数で計算する場合は、「exp」ではなく「cexp」だったような

[im_a_student]

それで解決しました！ありがとうございます！！

Answer 1

Fx[0] = aM*(exp(-IwT) - 1.0)/(-I*w);

複素数で計算する場合は、「exp」ではなく「cexp」