C程序(new).docVIP

下载本文档

2
0
约3.25千字
约 5页
2016-12-31 发布于贵州
举报
版权申诉

C程序(new).doc

1、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。。
2、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
4、该文档为VIP文档，如果想要下载，成为VIP会员后，下载免费。
5、成为VIP后，下载本文档将扣除1次下载权益。下载后，不支持退款、换文档。如有疑问请联系我们。
6、成为VIP后，您将拥有八大权益，权益包括：VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
7、VIP文档为合作方或网友上传，每下载1次，网站将根据用户上传文档的质量评分、类型等，对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档

int nextpow2(int a) { int i; int mul;//用于计算2的幂次相乘的结果 if(a = 0) { i = 1; } else { i = 0; do { i+=1; mul = 1; //用于计算2的幂次 for(int j = 1; j=i; j++) mul*=2; } while (mul a); } return i; } void bitrevorder(int* input, const int num, double* output) { for(int k = 0; k num ;k++){ int a = *(input+k); //N为输入数据的比特的位数 unsigned int N = sizeof(a)*8; unsigned int reverseNum = 0; for(unsigned int i = 0; i N; i++){ if((a1) == 1){ int j = N-1-i; unsigned int mul = 1;//用于2的幂次的相乘的结果 while(j0){//while循环进行2的幂次的计算 mul*=2; j--; } reverseNum+=mul; } a=a1; } *(output+k) = reverseNum; } } void myFFT(int* xn, int length, double* outReal, double* outImge,int outLength) { //说明： //xn是输入数组的首项指针 //length是输入数组的长度 //outReal是输出数组的实部 //outImge是输出数组的虚部 //outLength是输出数组的长度 //xn1的长度为NN const int NN = length; int xn1[NN]; for (int i = 0; i NN; i++) { xn1[i] = *(xn+i); } int m = nextpow2(NN); //用于计算2的m次幂 const int N = 1 m; const int M = NN+N; if(NN N) { int xn2[M] ={0}; double yn2Real[M] = {0.0}; double yn2Imag[M] = {0.0}; for(int i = 0; i N; i++) xn2[i] = xn1[i]; bitrevorder(xn2,M,yn2Real); for(int L = 1; L =m; L++){ int D = 1L; double u_Real = 1.0; double u_Imag = 1.0; double WM_Real = cos(2*3.14/D); double WM_Imag = -sin(2*3.14/D); for(int j = 1; j=D/2; j++){ for(int k = j; j=N; j+=D){ int kp = k + D/2 -1; double t_Real = yn2Real[kp]*u_Real; double t_Imag = yn2Imag[kp]*u_Imag; yn2Real[kp] = yn2Real[k] - t_Real; yn2Imag[kp] = yn2Imag[k] - t_Imag; yn2Real[k] = yn2Real[k] + t_Real; yn2Imag[k] = yn2Imag[k] + t_Imag; } u_Real*=WM_Real; u_Imag*=WM_Imag; } } outReal = yn2Real; outImge = yn2Imag; outLength = M; } else { double yn1Real[N] = {0.0}; double yn1Imag[N] = {0.0}; bitrevorder(xn1,N,yn1Real); for(int L = 1; L =m; L++){ int D = 1L; double u_Real = 1.0; doub