[问答] fft傅里叶变换问题

[复制链接]

技术员

发表于 2017-9-14 15:20:07   171 查看 0 回复 显示全部楼层 倒序浏览
分享
现在固定频率8k采样,如果是50hz的信号相对应就是160个采样点,现在要进行傅里叶变换,有没有大神能帮我看看代码  pwcz函数是在干嘛,整个代码能不能帮我屡屡思路
  1. #include <stdio.h>   
  2. #include <stdlib.h>   
  3. #include <math.h>   
  4.    
  5. #define SAMP 160   
  6. #define frequency 47.5             //为方便仿真频率先固定   
  7. #define MN (1.0*SAMP*50/47.5)      
  8. float tmp[21];   
  9. float ratio[20];   
  10. struct compx   
  11. {      
  12.     float real;   
  13.     float imag;   
  14. };   
  15. struct compx sum[SAMP+1];   
  16. float base;                 /*-基波-*/   
  17. float component_O[20];      /* -各次谐波分量-*/   
  18. float component_O_rt[20];   
  19. float ftmp;   
  20. float angle[20];   
  21. int N = SAMP;   
  22. #define RN SAMP*2   
  23. float AD_RESULT[RN];   
  24. float CZ_RESULT[RN];   
  25.    
  26. /*float amend_mt[20]=                   //谐波修正系数  
  27.     {  
  28.         1.0,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,0.9752,  
  29.         0.9752*1.08,0.9752*1.015,0.9752*1.022,0.9752*1.025,0.9752*1.030,  
  30.         0.9752*1.032,0.9752*1.042,0.9752*1.055,0.9752*1.065,0.9752*1.083  
  31.     };  
  32. */   
  33.    
  34.    
  35.    
  36. struct compx EE(struct compx b1,struct compx b2)   
  37. {   
  38.     struct compx b3;      
  39.     b3.real=(b1.real*b2.real-b1.imag*b2.imag);   
  40.     b3.imag=(b1.real*b2.imag+b1.imag*b2.real);   
  41.     return(b3);   
  42. }   
  43.    
  44. float pwcz(float x)   
  45. {   
  46.     float y;   
  47.     float e1,e2,e3;   
  48.     unsigned int k;   
  49.    
  50.     if((unsigned int)(x*10)%10<=5)   
  51.     {      
  52.         k=(unsigned int)x;   
  53.         printf("k=%d\t",k);   
  54.     }   
  55.     else   
  56.     {   
  57.         k=(unsigned int)x+1;   
  58.         printf("k1=%d\t",k);   
  59.     }   
  60.     e1=(x-k)*(x-(k+1))/2.0;   
  61.     e2=(x-(k-1))*(x-(k+1))/-1.0;   
  62.     e3=(x-(k-1))*(x-k)/2.0;   
  63.     y=e1*AD_RESULT[k-1]+e2*AD_RESULT[k]+e3*AD_RESULT[k+1];   
  64.     printf("y=%f\n", y);   
  65.     return y;   
  66. }   
  67.    
  68. static void draw_data(float *buf,int len,int mode)   
  69. {   
  70.     FILE *stream;   
  71.     float u;   
  72.     float temp;   
  73.     stream = fopen( "data.txt", "w" );   
  74.     int n = 0;   
  75.     while(len> 0)   
  76.     {   
  77.         temp = *buf;   
  78.         if( temp < 0)   
  79.         {   
  80.             temp = temp/4 + 60;   
  81.         }   
  82.         else   
  83.         {   
  84.             temp = temp/4 + 60;   
  85.         }   
  86.    
  87.         fprintf(stream, "%d %f\t", n++, *buf);   
  88.         for (u = 0; u < temp-1; u++)   
  89.         {   
  90.             fprintf( stream, " ");   
  91.         }   
  92.         fprintf( stream, "*\n");   
  93.         buf += 1;   
  94.         len -= 1;   
  95.     };   
  96.     fclose( stream );   
  97. }   
  98.    
  99.    
  100.    
  101. void main()   
  102. {   
  103.     for(int n=1;n<=2*N;n++)      // 正弦波加谐波计算测试   
  104.     {   
  105.         AD_RESULT[n]=220*sin(2*3.1415926*(n-1)/MN)   
  106.                         +0*sin(2*2*3.1415926*(n-1)/MN)                          //2次   
  107. //                          +22*sin(3*2*3.1415926*(n-1)/MN)                     //3次   
  108. //                              +22*sin(4*2*3.1415926*(n-1)/MN)                 //4次   
  109. //                                  +22*sin(5*2*3.1415926*(n-1)/MN)             //5次   
  110. //                                      +22*sin(6*2*3.1415926*(n-1)/MN)         //6次   
  111. //                                      +22*sin(7*2*3.1415926*(n-1)/MN)             //5次   
  112. //                                      +22*sin(8*2*3.1415926*(n-1)/MN)             //5次   
  113. //                                      +22*sin(9*2*3.1415926*(n-1)/MN)             //5次   
  114. //                                      +6.6*sin(10*2*3.1415926*(n-1)/MN)         //10次   
  115. //                                      +6.6*sin(11*2*3.1415926*(n-1)/MN)             //5次   
  116. //                                      +6.6*sin(12*2*3.1415926*(n-1)/MN)             //5次   
  117. //                                      +6.6*sin(13*2*3.1415926*(n-1)/MN)             //5次   
  118. //                                      +6.6*sin(14*2*3.1415926*(n-1)/MN)             //5次   
  119. //                                      +6.6*sin(15*2*3.1415926*(n-1)/MN)             //5次   
  120. //                                      +6.6*sin(16*2*3.1415926*(n-1)/MN)        //16次   
  121.                                         +11*sin(17*2*3.1415926*(n-1)/MN)   
  122.                                         +11*sin(18*2*3.1415926*(n-1)/MN)   
  123.                                         +11*sin(21*2*3.1415926*(n-1)/MN);   
  124. //                                      +22*sin(19*2*3.1415926*(n-1)/MN);   //19次      
  125. //      printf("AD_RESULT[%d]=%f\n", n, AD_RESULT[n]);   
  126.    
  127.    
  128.     }   
  129.     for(n=1;n<=SAMP;n++)     // 正弦波加谐波计算测试   
  130.     {   
  131.         printf("AD_RESULT[%d]=%f\t", n, AD_RESULT[n]);   
  132.         CZ_RESULT[n]=pwcz(n*50.0/frequency);   
  133.     }   
  134.    
  135.     float tmp2;   
  136.     float temp;   
  137.    
  138.     for(int m=1; m<22;m++)   
  139.     {   
  140.         for(n=1; n<=SAMP; n++)   
  141.         {   
  142.             sum[m].real+=CZ_RESULT[n]*cos(2*3.1415926*m*n/SAMP);   
  143.         }   
  144.         for(n=1; n<=SAMP; n++)   
  145.         {   
  146.             sum[m].imag+=CZ_RESULT[n]*sin(2*3.1415926*m*n/SAMP);   
  147.         }   
  148.    
  149.         tmp2=2*sqrt((sum[m].real*sum[m].real+sum[m].imag*sum[m].imag))/SAMP;   
  150.         temp = tmp2;   
  151.         component_O[m]=temp;   
  152.         ftmp = component_O[m] * component_O[m];   
  153.         component_O_rt[m] = sqrt(ftmp);   
  154.         ratio[m] = component_O_rt[m]*100/component_O_rt[1];   
  155. //      printf("component_O_rt[%d]=%f  ratio[%d] = %f  angle[%d]=%f\n", n, component_O_rt[n], n, ratio[n], n, angle[n]);   
  156.         printf("component_O_rt[%d]=%f  ratio[%d] = %f  \n", m, component_O_rt[m], m, ratio[m]);   
  157.     }      
  158.     draw_data(AD_RESULT, 2*SAMP, 0);   
  159.     draw_data(CZ_RESULT, SAMP, 0);   
  160. }  
复制代码


标签:fft变换
高级模式
您需要登录后才可以回帖 登录 | 注册

关闭

站长推荐 上一条 /7 下一条

快速回复 返回顶部 返回列表
-

推荐专区

技术干货集中营

专家问答

用户帮助┃咨询与建议┃版主议事

我的提问

工程师杂谈

工程师创意

工程师职场

论坛电子赛事

社区活动专版

发烧友活动

-

嵌入式论坛

ARM技术论坛

Android论坛

Linux论坛

单片机/MCU论坛

FPGA|CPLD|ASIC论坛

DSP论坛

嵌入式系统论坛

-

电源技术论坛

电源技术论坛

无线充电技术

-

硬件设计论坛

PCB设计论坛

电路设计论坛

电子元器件论坛

控制|传感

总线技术|接口技术

-

测试测量论坛

LabVIEW论坛

Matlab论坛

测试测量技术专区

仪器仪表技术专区

-

EDA设计论坛

multisim论坛

PADS技术论坛

Protel|AD|DXP论坛

Allegro论坛

proteus论坛|仿真论坛

EasyEDA-中国人自已的EDA工具

Orcad论坛

-

综合技术与应用

电机控制

智能电网

光电及显示

参考设计中心

汽车电子技术论坛

医疗电子论坛

-

开源硬件

DFRobot专区

树莓派论坛

智能硬件论坛

开发快智能硬件开发平台

Intel物联网开发者专区

Waveshare

乐美客SBC专区

Arduino论坛

BeagleBone论坛

机器人论坛

创客神器NanoPi

小钢炮CANNON

比派科技banana pi专区

-

无线通信论坛

无线通信技术专区

天线|RF射频|微波|雷达技术

-

IC设计论坛

芯片测试与失效分析

Mixed Signal/SOC[数模混合芯片设计]

Analog/RF IC设计

设计与制造封装测试

-

个人版区

阿东Verilog技术专版

直流马达驱动电路设计

LabVIEW英雄联盟

特权同学FPGA专区

-

厂商专区

灵动微电子 MM32

盈鹏飞嵌入式

TI论坛

TI Deyisupport社区

芯灵思嵌入式论坛

Tisan

米尔科技

庆科社区

WIZnet技术专区

Cypress技术论坛

飞凌嵌入式

Qualcomm技术论坛

英创嵌入式

机智云GoKit论坛

-

检测技术与质量

电磁兼容(EMC)设计与整改

安规知识论坛

检测与认证

-

消费电子论坛

手机技术论坛

平板电脑/mid论坛

音视/视频/机顶盒论坛

-

电子论坛综合区

聚丰众筹官方社区

新人报道区

聚丰供应链

-

论坛服务区

-

供求信息发布

供需广告

招聘┃求职发布区

电子展览展会专区