[问答]

fft傅里叶变换问题

问答对人有帮助，内容完整，我也想知道答案 0 现在固定频率8k采样，如果是50hz的信号相对应就是160个采样点，现在要进行傅里叶变换，有没有大神能帮我看看代码 pwcz函数是在干嘛，整个代码能不能帮我屡屡思路 #include #include #include #define SAMP 160 #define frequency 47.5 //为方便仿真频率先固定 #define MN (1.0SAMP50/47.5) float tmp[21]; float ratio[20]; struct compx { float real; float imag; }; struct compx sum[SAMP+1]; float base; /-基波-/ float component_O[20]; /* -各次谐波分量-/ float component_O_rt[20]; float ftmp; float angle[20]; int N = SAMP; #define RN SAMP2 float AD_RESULT[RN]; float CZ_RESULT[RN]; /float amend_mt[20]= //谐波修正系数 { 1.0,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752, 0.97521.08,0.97521.015,0.97521.022,0.97521.025,0.97521.030, 0.97521.032,0.97521.042,0.97521.055,0.97521.065,0.97521.083 }; / struct compx EE(struct compx b1,struct compx b2) { struct compx b3; b3.real=(b1.realb2.real-b1.imagb2.imag); b3.imag=(b1.realb2.imag+b1.imagb2.real); return(b3); } float pwcz(float x) { float y; float e1,e2,e3; unsigned int k; if((unsigned int)(x10)%10<=5) { k=(unsigned int)x; printf("k=%dt",k); } else { k=(unsigned int)x+1; printf("k1=%dt",k); } e1=(x-k)(x-(k+1))/2.0; e2=(x-(k-1))(x-(k+1))/-1.0; e3=(x-(k-1))(x-k)/2.0; y=e1AD_RESULT[k-1]+e2AD_RESULT[k]+e3AD_RESULT[k+1]; printf("y=%fn", y); return y; } static void draw_data(float buf,int len,int mode) { FILE stream; float u; float temp; stream = fopen( "data.txt", "w" ); int n = 0; while(len> 0) { temp = buf; if( temp < 0) { temp = temp/4 + 60; } else { temp = temp/4 + 60; } fprintf(stream, "%d %ft", n++, buf); for (u = 0; u < temp-1; u++) { fprintf( stream, " "); } fprintf( stream, "n"); buf += 1; len -= 1; }; fclose( stream ); } void main() { for(int n=1;n<=2N;n++) // 正弦波加谐波计算测试 { AD_RESULT[n]=220sin(23.1415926(n-1)/MN) +0sin(223.1415926(n-1)/MN) //2次 // +22sin(323.1415926(n-1)/MN) //3次 // +22sin(423.1415926(n-1)/MN) //4次 // +22sin(523.1415926(n-1)/MN) //5次 // +22sin(623.1415926(n-1)/MN) //6次 // +22sin(723.1415926(n-1)/MN) //5次 // +22sin(823.1415926(n-1)/MN) //5次 // +22sin(923.1415926(n-1)/MN) //5次 // +6.6sin(1023.1415926(n-1)/MN) //10次 // +6.6sin(1123.1415926(n-1)/MN) //5次 // +6.6sin(1223.1415926(n-1)/MN) //5次 // +6.6sin(1323.1415926(n-1)/MN) //5次 // +6.6sin(1423.1415926(n-1)/MN) //5次 // +6.6sin(1523.1415926(n-1)/MN) //5次 // +6.6sin(1623.1415926(n-1)/MN) //16次 +11sin(1723.1415926(n-1)/MN) +11sin(1823.1415926(n-1)/MN) +11sin(2123.1415926(n-1)/MN); // +22sin(1923.1415926(n-1)/MN); //19次 // printf("AD_RESULT[%d]=%fn", n, AD_RESULT[n]); } for(n=1;n<=SAMP;n++) // 正弦波加谐波计算测试 { printf("AD_RESULT[%d]=%ft", n, AD_RESULT[n]); CZ_RESULT[n]=pwcz(n50.0/frequency); } float tmp2; float temp; for(int m=1; m<22;m++) { for(n=1; n<=SAMP; n++) { sum[m].real+=CZ_RESULT[n]cos(23.1415926mn/SAMP); } for(n=1; n<=SAMP; n++) { sum[m].imag+=CZ_RESULT[n]sin(23.1415926mn/SAMP); } tmp2=2sqrt((sum[m].realsum[m].real+sum[m].imagsum[m].imag))/SAMP; temp = tmp2; component_O[m]=temp; ftmp = component_O[m] * component_O[m]; component_O_rt[m] = sqrt(ftmp); ratio[m] = component_O_rt[m]100/component_O_rt[1]; // printf("component_O_rt[%d]=%f ratio[%d] = %f angle[%d]=%fn", n, component_O_rt[n], n, ratio[n], n, angle[n]); printf("component_O_rt[%d]=%f ratio[%d] = %f n", m, component_O_rt[m], m, ratio[m]); } draw_data(AD_RESULT, 2SAMP, 0); draw_data(CZ_RESULT, SAMP, 0); } 复制代码 0
2017-9-14 15:20:07　　评论淘帖0 邀请回答您可以邀请以下用户，快速回答问题 × dianzijie5 该类别下有 16 个回答。邀请回答嵌入式-H 该类别下有 16 个回答。邀请回答 hgimtk 该类别下有 16 个回答。邀请回答 hucc 该类别下有 14 个回答。邀请回答舍我其谁84 该类别下有 13 个回答。邀请回答 HCPcry 该类别下有 12 个回答。邀请回答冰箱洗衣机该类别下有 12 个回答。邀请回答 iettke 该类别下有 11 个回答。邀请回答学生物的程序猿该类别下有 11 个回答。邀请回答小刘该类别下有 11 个回答。邀请回答 sjjs001 该类别下有 11 个回答。邀请回答 YOYOOO 该类别下有 10 个回答。邀请回答 bigbangboom 该类别下有 10 个回答。邀请回答 siyugege 该类别下有 9 个回答。邀请回答 drakannie 该类别下有 9 个回答。邀请回答 wuli北该类别下有 9 个回答。邀请回答广州洋钒该类别下有 9 个回答。邀请回答 zengweishuai 该类别下有 9 个回答。邀请回答 mintsy 该类别下有 8 个回答。邀请回答 fanB 该类别下有 8 个回答。邀请回答举报阿政相关推荐 • 非2的n次方求快速傅里叶变换可以使用DSP库函数中的FFT吗？ 3242 • 二傅里叶变换是什么？如何求傅里叶变换？ 2086 • 学习傅里叶变换意义和方法 2458 • 有关傅里叶变换的相关知识 2299 • 傅里叶变换可分为哪几类 1804 • 请问怎么用labview对一个连续信号连续变化的加速度值进行傅里叶变换？ 2323 • 急急急！！大神进来看看，帮我设计一个振动信号文件通过快速傅里叶变换分析其频域的后面板！！ 4149 • fft的相位问题 1935 • 傅里叶变换和拉普拉斯变换有什么区别 2654 • DSP操作教程 4-7 快速傅立叶变换（FFT）算法（CCS显示） 11252