复数自运算注意--小提示

         在编图像的FFT的递归实现时，没有使用STL的complex模板，直接对实部和虚部进行操作，发现频谱图像结果与OpenCV库实现的结果不一样，lena图像的FFT像是粗糙的布一样，下面左边是OpenCV生成的，右边是自己编的生成的。

 

      左边的可以看出左上--右下对角线方向有频带，而右边的看不出来，出错了。使用complex模板来计算则正常，那么应该是自己写的复数的运算出错了，仔细检查没看出来，单独写了几句程序测试了一下之后终于发现问题了。

      程序中有W=W*Wn的一个自乘运算，如果使用complex模板很简单，但是自己去写复数乘法的时候忽略了一点。

哪一点？

错误的程序：

		w_r = w_r*wn_r-w_i*wn_i;
		w_i = w_r*wn_i+w_i*wn_r;

正确的程序：

		T t_r = w_r;
		w_r = t_r*wn_r-w_i*wn_i;
		w_i = t_r*wn_i+w_i*wn_r;

      这下一目了然，做自乘运算的时候，计算虚部要用到原先的实部，而错误的程序中使用的是计算后的实部，所以。。。应该将原先的值暂存一下，计算时使用原先的值，下面的程序先计算的实部，所以只需要暂存实部就可以，如果先计算虚部那么很明显应该暂存虚部。

      复数的自加减由于实部虚部不混合所以不用担心，复数的除法的话也涉及到实部虚部的系数混合运算所以在自除时也要注意把原数据暂存的问题。