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