fft2r2.txt

来自「用c语言实现了图像二维快速傅里叶变换」· 文本 代码 · 共 46 行

void fft2r21(double A[2][N1],double B[2][N1],int ifft)；
void fft2r22(double A[2][N2],double B[2][N2],int ifft)；
void fft2r21(double A[2][N1],double B[2][N1],int ifft)
{
	unsigned long int k,j;
	double c[2][N1];
	if (ifft!=-1)
	{
		for(k=0;k<N1;k++)
		{ c[0][k]=A[0][k];
		c[1][k]=B[0][k];}
		fft21(c,1);
		A[0][0]=c[0][0];
		A[1][0]=0.0
			B[0][0]=c[1][0];
		B[1][0]=0.0;
		for (j=1;j<N1;j++)
		{ A[0][j]=0.5*(c[0][j]+c[0][N1-j]);
		A[1][j]=0.5*(c[1][j]-c[1][N1-j]);
		B[0][j]=0.5*(c[0][j]+c[0][N1-j]);
		B[1][j]=0-.5*(c[1][j]-c[1][N1-j]);
		}
	}
}
void fft2r22(double A[2][N2],double B[2][N2],int ifft)
{
	unsigned long int k,j;
	double c[2][N2];
	if (ifft!=-1)
	{
		for(k=0;k<N2;k++)
		{ c[0][k]=A[0][k];
		c[1][k]=B[0][k];}
		fft22(c,1);
		A[0][0]=c[0][0];
		A[1][0]=0.0
			B[0][0]=c[1][0];
		B[1][0]=0.0;
		for (j=1;j<N2;j++)
		{ A[0][j]=0.5*(c[0][j]+c[0][N2-j]);
		A[1][j]=0.5*(c[1][j]-c[1][N2-j]);
		B[0][j]=0.5*(c[0][j]+c[0][N2-j]);
		B[1][j]=0-.5*(c[1][j]-c[1][N2-j]);
		}
	}
}