如何在C ++中使用复数“i”

时间:2013-07-29 13:52:34

标签: c++ overloading complex-numbers complextype

我现在编写一个简单的DFT算法,我想在复指数中使用复数i。我看到有人使用#include<complex>#include<cmath>,然后他们使用重载的符号I,例如exp(2*I)。但它似乎在我的visual studio编译器中不起作用。那么,任何人都可以给出一个使用复指数的简单例子吗?谢谢!

6 个答案:

答案 0 :(得分:7)

我最近也提到了这个问题,为未来的读者找到了一个简单的方法:

只需使用<complex>库,如下所示

#include <iostream>
#include <complex>
using namespace std ;

int main(int argc, char* argv[])
{
    const   complex<double> i(0.0,1.0);    
    cout << i << endl ;

    return(0) ;
}

答案 1 :(得分:5)

这是一个简短的完整示例:

#include <iostream>
#include <complex>
#include <cmath>

using namespace std;
typedef complex<double> dcomp;

main() {
  dcomp i;
  dcomp a;
  double pi;
  pi = 2 * asin(1);
  i = -1;
  i = sqrt(i);
  a = exp(2*pi*i);
  cout << "i is " << i << "and Euler was right: e(i pi) = " << a << endl;
} 

使用g ++进行测试

答案 2 :(得分:3)

您可以找到详细信息here

一个简单的方法是

#include <complex>

using std::complex;
const double pi = 3.1415;
void foo()
{
    complex<double> val(polar(1, pi/2.0); Create a complex from its olar representation
}

答案 3 :(得分:1)

C ++中的以下代码显示了用于实现虚数j的宏。众所周知,在编程中,术语i和j通常用作计数器变量。我改为使用大写字母J代表虚数以避免混淆。

/ * dcomplex.h

#ifndef DCOMPLEX_H_
#define DCOMPLEX_H_
#define J dcomplex(0.0,1.0)
typedef std::complex<double> dcomplex;
#endif /* DCOMPLEX_H_ */

使用此宏,可以在主代码中使用虚数J [与复杂库一起]。其使用示例如下所示:

....
....
#include <complex>
#include "dcomplex.h"

....
....
 tmp = tmp + t[n]*exp( (2.0*PI*(double)n*(double)l/(double)tab_size)*J );
....

...

其中tmp,t [n]是复数类型的变量,J是虚数。变量n,l和tab_size是整数类型。常数PI是众所周知的常量3.14 ...函数exp()被重载以处理复数。 [注:此代码示例是简单DFT]

的一部分

使用此宏,代码更具可读性。

答案 4 :(得分:1)

另一种方法是在C ++ 14之后使用std::literals::complex_literals::operator""i

#include <iostream>
#include <complex>

int main() {
    using namespace std::complex_literals;
    auto c = 1.0 + 3.0i;
    std::cout << "c = " << c << '\n';
}

输出:

c = (1,3)

答案 5 :(得分:0)

作为无理数的

pi 无法用 double 完全表示。 pi的不精确近似的 cos 可能产生接近但可能不完全是1的结果。同样 sin 的不精确近似的不精确近似> pi 喜欢导致数字的幅度非常小,可能不完全为0.为什么不将I定义为 std :: complex(0.0,1.0)和避免无端的不精确。