views:

161

answers:

6

I need to compute imaginary exponential in C.

As far as I know, there is no complex number library in C. It is possible to get e^x with exp(x) of math.h, but how can I compute the value of e^(-j), where j = sqrt(-1)?

+3  A: 

Note that exponent of complex number equals:

e^(jx) = cos(x)+j*sin(x)

Then:

e^(-j) = cos(-1)+j*sin(-1)
Vladimir
+5  A: 

Using the Euler's Formula you have that e^-i == cos(1) - i*sin(1)

Jack
e^(-k) = cos(k) - i*sin(k), you mean.
Deep-B
yes, but in his case k = 1. In any case the general form is e^(-ki) = cos(k)-i*sin(k), not e^(-k)
Jack
+1  A: 

Just use the cartesian form

if z = m*e^j*(arg);

re(z) = m * cos(arg);
im(z) = m * sin(arg);
Tom
+2  A: 

e^-j is just cos(1) - j*sin(1), so you can just generate the real and imaginary parts using real functions.

Paul R
+11  A: 

In C99, there is a complex type. Include complex.h; you may need to link with -lm on gcc. Note that Microsoft Visual C does not support complex; if you need to use this compiler, maybe you can sprinkle in some C++ and use the complex template.

I is defined as the imaginary unit, and cexp does exponentiation. Full code example:

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

int main() {
    complex x = cexp(-I);
    printf("%lf + %lfi\n", creal(x), cimag(x));
    return 0;
}

See man 7 complex for more information.

Thomas
+1 Didn't know about that.
Tom
I needed to do this manually.. Thanks for the future reference, I'll keep this in mind.
Erkan Haspulat
A: 

Is calling a c++ function a solution for you? The C++ STL has a nice complex-class and boost also has to offer some nice options. Write a function in C++ and declare it "extern C"

extern "C" void myexp(float*, float*);

#include <complex>

using std::complex;

void myexp (float *real, float *img )
{
  complex<float> param(*real, *img);
  complex<float> result = exp (param);
  *real = result.real();
  *img = result.imag();
}

Then you can call the function from whatever C-code you rely on ( Ansi-C, C99, ...).

#include <stdio.h>

void myexp(float*, float*);

int main(){
    float real = 0.0;
    float img = -1.0;
    myexp(&real, &img);
    printf ("e^-i = %f + i* %f\n", real, img);
    return 0;
}
Lucas
Thanks for the explanation but C++ is no good for me with this
Erkan Haspulat