请问各位朋友，在C++中如何实现复数数组？

[ruixiang6]

请问各位朋友，在C++中如何实现复数数组？我想有标准库中的complex来实现，为什么不行呢？我的是VC2005，程序如下：



#include "stdafx.h"
#include "math.h"
#include <complex>

void DFT(complex<float> a[],complex<float> b[],int N)
{
for (int k = 0;k<N;k++)
{
b[k]=(0,0);

for(int n = 0;n<N;n++)
{
complex<float> mm = (0,(2*3.14/N)*k*n);
b[k]=b[k]+(a[n]*(exp(mm));
}
}
}


int main()
{
complex<float> x[3]={(1,0),(1,0),(1,0)};
complex<float> y[3];

DFT(x,y,3);

for(int k=0;k<3;k++)
{
printf("y[%d] = %f+i%f\n",k,y[k]);
}

return 0;
}

[Y___Y]

b[k]=b[k]+(a[n]*(exp(mm)));少了一个括弧

[ruixiang6]

改了，还是不行。。

[ruixiang6]

我加了括弧为什么还是不行呀？？？

提示：错误 1 error C2065: “complex”: 未声明的标识符 g:\work\visual studio 2005\projects\dosdft2\dosdft2\dosdft2.cpp 5

[手插兜兜]

在标准库中，用 vector,就是数组。

[xyxfly]

error C2065: “complex”: 未声明的标识符 g:\work\visual studio 2005\projects\dosdft2\dosdft2\dosdft2.cpp 5



using namespace std;

[ruixiang6]

且不管结果如何，为什么你加了括号就能编译，而我的有错误呢？

[ruixiang6]

我加了，也不行。。。

[ruixiang6]

加上using namespace std;行了，可是为何结果不对呀，程序哪里有错吗？

[xyxfly]

说实话，我都不知道你在算什么 ：）

[ruixiang6]

简单一点，看下面的程序，为何显示全是0.000000？

#include "math.h"
#include <complex>
using namespace std;

int main()
{
complex<float> x[3]={(1,0),(2,1),(3,4)};

for(int k=0;k<3;k++)
{
printf("x[%d] = %f",k,x[k]);
}

return 0;
}

[xyxfly]

你把一个复数 %f输出？

[xyxfly]

http://www.programfan.com/club/showbbs.asp?id=203772&page=last

http://lix66.spaces.live.com/blog/

自己重载

[ruixiang6]

谢谢小虾米呀，我来看看。。

[ruixiang6]

C++标准库中有<complex>，我想直接用。但不知道怎么定义复数数组，怎么向它赋值，并输出？
求助！

[菜牛]

printf( "x[%d] = %f,%f ",k,x[k].real, x[k].imag);
这样试试

[ruixiang6]

printf( "x[%d] = %f,%f ",k,x[k].real, x[k].imag);是错的，偶加了括号变成
printf( "x[%d] = %f,%f ",k,x[k].real（）, x[k].imag（）);

但是按照下面的程序输出却是，0,0 1,0 4,0 为什么实部变也了原虚部，虚部全0？

#include "math.h"
#include <complex>
using namespace std;

int main()
{
complex<float> x[3]={(1,0),(2,1),(3,4)};

for(int k=0;k<3;k++)
{
printf("x[%d] = %f",k,x[k]);
}

return 0;
}

[fuyifan]

#include "math.h"
#include <complex>
#include <iostream>

using namespace std;

void DFT(complex<float> a[],complex<float> b[],int N)
{
for (int k = 0;k<N;k++)
{
b[k]=complex<float>(0,0);
//b[k]=(0,0);

for(int n = 0;n<N;n++)
{
complex<float> mm = complex<float>(0,(2*3.14/N)*k*n);
//complex<float>mm = (0,(2*3.14/N)*k*n);
b[k]=b[k]+(a[n]*(exp(mm)));
}
}
}


int main()
{
complex<float> x[3]={complex<float>(1,0),complex<float>(1,0),complex<float>(1,0)};
complex<float> y[3];

DFT(x,y,3);

for(int k=0;k<3;k++)
{
cout<<"y["<<k<<"]="<<y[k];
//printf("y[%d] = %f+i%f\n",k,y[k]);
}

return 0;
}

[ruixiang6]

楼上的朋友太厉害了，好了，谢谢！

也非常感谢小虾火，相互中，还有其他帮助我的朋友！

[ruixiang6]

我来贴一下重要代码。一共两个
第一个，是我没有使用标准库中complex来表示复数，是用C写的，是编译通过的。

下面是一个函数
void CdftcDlg::DFT(float a1[], float a2[], float b1[], float b2[], int N)
{
for (int k = 0;k<N;k++)
{
b1[k]=0;
b2[k]=0;

for (int n = 0;n<N;n++)
{
b1[k]=b1[k]+(a1[n]*cos((2*3.14/N)*k*n)-a2[n]*sin((2*3.14/N)*k*n));
b2[k]=b2[k]+(a1[n]*sin((2*3.14/N)*k*n)+a2[n]*cos((2*3.14/N)*k*n));
}
}

}



下面是单击按扭控件产生的事件，以期在文本框中取得数据，并在上面函数的作用下显示到文本框中：

void CdftcDlg::OnBnClickedButton1()
{

float x1[3],x2[3];

UpdateData();

x1[0]=m_edit1;
x1[1]=m_edit3;
x1[2]=m_edit5;
x2[0]=m_edit2;
x2[1]=m_edit4;
x2[2]=m_edit6;

float y1[3],y2[3];

CdftcDlg::DFT(x1,x2,y1,y2,3);

m_edit7 = y1[0];
m_edit8 = y2[0];
m_edit9 = y1[1];
m_edit10= y2[1];
m_edit11= y1[2];
m_edit12= y2[2];

UpdateData(FALSE);

}



下面是第二个，是我想使用标准库存中complex来表示复数，是想方便一点的，没想到通不过呀，哪知道有错误。

函数：
void CdftcDlg::DFT(complex<float> a[],complex<float> b[],int N)
{

for (int k = 0;k<N;k++)
{
b[k]=complex<float>(0,0);

for(int n = 0;n<N;n++)
{
complex<float> mm = complex<float>(0,(2*3.14/N)*k*n);
b[k]=b[k]+(a[n]*(exp(mm)));
}
}

}


单击事件：


void CdftcDlg::OnBnClickedButton1()
{

complex<float> x[3];

UpdateData();

x[0] = complex<float>(m_edit1,m_edit2);

x[1] = complex<float>(m_edit3,m_edit4);

x[2] = complex<float>(m_edit5,m_edit6);



complex<float> y[3];

CdftcDlg::DFT(x,y,3);

m_edit7 = y[0].real();
m_edit8 = y[0].imag();
m_edit9 = y[1].real();
m_edit10= y[1].imag();
m_edit11= y[2].real();
m_edit12= y[2].imag();

UpdateData(FALSE);

}

[jijun1018]

共同学习，顶一下吧