FFT C语言实现遇到的问题

for(int i = 0; i< 16;i++)
	{
		temp1[i].real = i+1;
		temp2[i].real = i+6;
	}

这是长度为16的两组数据，分别是1-16 和 6-21 。虚部不管，全部赋0 。
在后面补了15个0，使 “两个信号通过补0，变成长度等于两个序列长度之和减一”

//fft
	fft(fft_point, len, temp1, fftemp1);
	fft(fft_point, len, temp2, fftemp2);
	//print fft
	printf("\n数组一FFT结果fft1:\n\t");
	for (int i = 0;i < fft_point;i++) {
		printf("%lf\t",fftemp1[i].real);
	}
	printf("\n数组二FFT结果fft2:\n\t");
	for (int i = 0;i < fft_point;i++) {
		printf("%lf\t",fftemp2[i].real);
	}

这是FFT计算输出的32个结果

//ifft array and print
	ifft(fft_point,fftemp1,ifftemp1);
	ifft(fft_point,fftemp2,ifftemp2);
	printf("\n数组一IFFT 验算ifft1:\n\t");
	for (int i = 0;i < len;i++) {
		printf("%lf\t",model(ifftemp1[i]));
	}
	printf("\n数组二IFFT 验算ifft2:\n\t");
	for (int i = 0;i < len;i++) {
		printf("%lf\t",model(ifftemp2[i]));
	}

经过FFT，IFFT，和原先数据一样，可以验证 正确性
还是1-16和6-21

conv(fft_point,fftemp1,fftemp2,cov12);

void conv(int length, complex* xk, complex* yk, complex* result)
{
	int i;
	for (i = 0;i < length;i++) {
		result[i].real = xk[i].real*yk[i].real - xk[i].imaginary*yk[i].imaginary;
		result[i].imaginary = xk[i].real*yk[i].imaginary + yk[i].real*xk[i].imaginary;
	}
}

虚部全部为0，应该不需要考虑共轭的事情。
频域直接相乘，输出结果为这样的32个数据。

//ifft
	ifft(fft_point,cov12,result12);
	//print ifft
	printf("\n结果IFFT到时域ifft result:\n\t");
	for (int i = 0;i < fft_point;i++) {
		printf("%lf\t",result12[i].real);
	}
	//find max
	double val = 0.0;
	int maxpos = 0;
	for(int i = 0;i < fft_point; i++)
	{
		if(val < result12[i].real)
		{
			val = result12[i].real;
			maxpos = i;
		}
	}
	printf("\nMAXpos = %d:\n",maxpos);
	printf("\nMAXval = %lf:\n",val);

又到了看似成功其实失败的时候，将频域相乘的结果转换到时域
！最大值对应的下标居然是20！！！
不是15！
这个都不需要到matlab看结果就知道是错的了，想想就知道不需要时间延迟，就是最大值对上最大值，应该延时0的时候就是互相关的最大值啊。

for(int i = 0; i< 16;i++)
	{
		temp1[i].real = i+1;
		temp2[i].real = i+1;
	}

这下两个信号一模一样了

但是结果还不是15！
是21了。。更离谱了。。
唯一的一点安慰是数值没有算错

感觉可能在补零，或者相乘的时候出了点小问题。。。但是现在思路停滞，找不到原因。

“由于前面在快速计算相关结果的时候，就利用了快速傅里叶变换，所以可以在最后一步进行反傅里叶变换的时候，先进行补零，然后在进行”
这句话也没有完全理解，在频域数据的哪里补零呢。。。？左右一半 or 两个数据之间 or 左边 or 右边？