想用FFT算法計(jì)算電壓的有效值。試驗(yàn)用數(shù)據(jù)由Ui=1024*[1+sin（2*i*PI/N）]（i=0，1，...N-1）得到理想的32點(diǎn)數(shù)值，存入數(shù)組data[32].Re中，而data[32].Im=0。對(duì)data[]逆序排列后，進(jìn)行FFT運(yùn)算，預(yù)想結(jié)果是：預(yù)想的運(yùn)算結(jié)果data[]的實(shí)部之和應(yīng)為AD_full的值，虛部之和為0???。但是以下程序運(yùn)行結(jié)果data[]的實(shí)部之和為0。查了多半天，沒(méi)有找到原因，頭都大了。請(qǐng)大家?guī)兔纯矗e(cuò)在哪里。謝謝！*/???
??/*test_data（float???u_test[]）和revers_order（struct???complx???A[]）經(jīng)驗(yàn)證沒(méi)有問(wèn)題，錯(cuò)誤很可能在FFT（）中。*/???
??#define???PI???3.141593???
??#define???AD_full???1024???
??#define???N???32???
??#include??????
??#include??????
??struct???complx???
??{???
??????double???Re;???
??????double???Im;???
??};???
????
??void???test_data（float???u_test[]）;/*產(chǎn)生32個(gè)理想的電壓數(shù)據(jù)點(diǎn)*/???
??void???FFT（struct???complx???A[]）;?????????
??/*對(duì)復(fù)數(shù)進(jìn)行FFT運(yùn)算，A[]的實(shí)部為電壓，虛部為0。預(yù)想的運(yùn)算結(jié)果A[]的實(shí)部之和應(yīng)為AD_full的值，虛部之和為0???*/???
??void???revers_order（struct???complx???A[]）;???/*???按照數(shù)組下標(biāo)的二進(jìn)制逆序重新排列數(shù)據(jù)*/???
????
??main（）???
??{???
????int???i;???
????struct???complx???data[N];???
????float???sum_data=0.0u_test[N]re_test[N]im_test[N]sum_re=0.0sum_im=0.0;???
????test_data（u_test）;???
????printf（“\ntest_data???are:\n“）;???
????for（i=0;i????{???
????????data[i].Re=u_test[i];???
????????data[i].Im=0.0;???
????????sum_data=sum_data+data[i].Re*data[i].Re;???
????????printf（“d[%2d].Re=%f?????d[%2d].Im=%1.1f\n???“idata[i].Reidata[i].Im）;???
????}???
????getch（）;???
????revers_order（data）;???
????printf（“\ntest???revers???data???are:\n“）;???
????for（i=0;i??????{???
????????printf（“d[%2d].Re=%f?????d[%2d].Im=%1.1f???“idata[i].Reidata[i].Im）;???
????????if（（i+1）%2==0）???printf（“\n“）;???
??????}???
????getch（）;???
????FFT（data）;???
????printf（“\ntest???FFT???data???are:\n“）;???
????for（i=0;i??????{???
????????sum_re=sum_re+data[i].Re;???
????????sum_im=sum_im+data[i].Im;???
????????printf（“d[%2d].Re=%f?????d[%2d].Im=%f???“idata[i].Reidata[i].Im）;???
????????if（（i+1）%2==0）???printf（“\n“）;???
??????}???
????????
????
??printf（“\n\nsum_re=%fsum_im=%fdata_sum=%f“sum_resum_imsqrt（sum_data/N））;???
????getch（）;???
??}???
????
??void???FFT（struct???complx???A[]）???
??{???
????struct???complx???UWT;???
????int???LELE1IJIP;???
????int???M=5;???
????double???temp;???
????i