找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 5727|回复: 7
打印 上一主题 下一主题
收起左侧

PIC单片机FFT音乐频谱Proteus仿真程序

  [复制链接]
跳转到指定楼层
楼主
因为贴子附件大小问题,wav音乐文件没有附上,请自行添加到目录,打开FFT.pdsprj文件,然后双击Audio图标指定音乐wav位置,开始运行即可,稍后再发下电脑里以前做的东东


单片机源程序如下:
  1. /********************************************************************************
  2. ;*  File:FFT v1.0                                                               *  
  3. ;*  Date:2015-10-22                                                             *
  4. ;*  Version:1.0                                                                 *
  5. ;*  Author:Ben                                                                  *   
  6. ;*  Wechat:feeling1791                                                          *
  7. ;********************************************************************************/
  8. /*
  9. // Configuration Mask Definitions
  10. #define CONFIG_ADDR        0x2007
  11. // Oscillator
  12. #define EXTCLK                0x3FFF        // External RC Clockout
  13. #define EXTIO                0x3FFE        // External RC No Clock
  14. #define INTCLK                0x3FFD        // Internal RC Clockout
  15. #define INTIO                0x3FFC        // Internal RC No Clock
  16. #define EC                0x3FFB        // EC
  17. #define HS                0x3FFA        // HS
  18. #define XT                0x3FF9        // XT
  19. #define LP                0x3FF8        // LP
  20. // Watchdog Timer
  21. #define WDTEN                0x3FFF        // On
  22. #define WDTDIS                0x3FF7        // Disabled / SWDTEN control
  23. // Power Up Timer
  24. #define PWRTDIS                0x3FFF        // Off
  25. #define PWRTEN                0x3FEF        // On
  26. // Master Clear Enable
  27. #define MCLREN                0x3FFF        // MCLR function is enabled
  28. #define MCLRDIS                0x3FDF        // MCLR functions as IO
  29. // Code Protect
  30. #define UNPROTECT        0x3FFF        // Code is not protected
  31. #define CP                0x3FBF        // Code is protected
  32. // Data EE Read Protect
  33. #define DUNPROTECT        0x3FFF        // Do not read protect EEPROM data
  34. #define CPD                0x3F7F        // Read protect EEPROM data
  35. // Brown Out Detect
  36. #define BORDIS                0x3CFF        // BOD and SBOREN disabled
  37. #define SWBOREN                0x3DFF        // SBOREN controls BOR function (Software control)
  38. #define BORXSLP                0x3EFF        // BOD enabled in run, disabled in sleep, SBOREN disabled
  39. #define BOREN                0x3FFF        // BOD Enabled, SBOREN Disabled
  40. // Internal External Switch Over Mode
  41. #define IESOEN                0x3FFF        // Enabled
  42. #define IESODIS                0x3BFF        // Disabled
  43. // Monitor Clock Fail-safe
  44. #define FCMEN                0x3FFF        // Enabled
  45. #define FCMDIS                0x37FF        // Disabled
  46. // Low Voltage Programming
  47. #define LVPDIS                0x2FFF        // Disabled
  48. #define LVPEN                0x3FFF        // Enabled
  49. // In-Circuit Debugger Mode
  50. #define DEBUGEN                0x1FFF        // Enable ICD2 debugging
  51. #define DEBUGDIS        0x3FFF        // Disable ICD2 debugging
  52. #define CONFIG_ADDR2        0x2008
  53. // Brown-out Reset Voltage
  54. #define BORV21                0x3EFF        // 2.1 Volts
  55. #define BORV40                0x3FFF        // 4.0 Volts
  56. // Flash Memory Write Protection
  57. #define WP0                0x3DFF        // Protect 0h-0FFh
  58. #define WP1                0x3BFF        // Protect lower half of flash
  59. #define WP2                0x39FF        // Protect all of flash
  60. */
  61. #include <pic.h> //pic16f887
  62. __CONFIG(INTIO&WDTDIS&PWRTDIS&MCLREN&BOREN&LVPDIS&DUNPROTECT&DEBUGDIS&UNPROTECT&IESOEN&FCMEN);  

  63. #define LOG_NUMBER 5          

  64. unsigned char AD_count;
  65. const unsigned char BRTable[32] = {0,16,8,24,4,20,12,28,2,18,10,26,6,22,14,30,
  66.                                   1,17,9,25,5,21,13,29,3,19,11,27,7,23,15,31};
  67. const  char SIN_TAB[32] =
  68.                          { 0 ,25 ,49 ,71 ,90 ,106 ,117 ,125 ,127 ,
  69.                            125 ,117 ,106 ,90 ,71 ,49 ,25 ,0 ,
  70.                            -25 ,-49 ,-71 ,-90 ,-106 ,-117 ,-125 ,-127 ,
  71.                            -125 ,-117 ,-106 ,-90 ,-71 ,-49 ,-25 };
  72.                        
  73. const  char COS_TAB[32] =
  74.                          {127 ,125 ,117 ,106 ,90 ,71 ,49 ,25 ,0 ,
  75.                           -25 ,-49 ,-71 ,-90 ,-106 ,-117 ,-125 ,-127 ,
  76.                          -125 ,-117 ,-106 ,-90 ,-71 ,-49 ,-25 ,0 ,
  77.                            25 ,49 ,71 ,90 ,106 ,117 ,125  };

  78. unsigned char disp_time;
  79. unsigned char Led_scan_count;

  80. bit AD_readed_flag = 0;
  81. bit AD_getonce_flag = 0;

  82. bank3 const unsigned char Disp_Number[8]={0x01,0x03,0x07,0x0f,0x1f,0x3f,0x7f,0xff};//0-7
  83. bank3 const unsigned char Disp_peak[8]={0x01,0x02,0x04,0x08,0x10,0x20,0x40,0x80};//0-7
  84. bank3 unsigned char led_num[16];

  85. bank2 unsigned char MAIN_TAB[16];
  86. bank1 unsigned char PEAK_TAB[16];
  87. bank3 unsigned char Peak_speed[16];

  88. bank1 int Fft_Real[32];
  89. bank2 int Fft_Image[32];

  90. unsigned int FFT_value;

  91. /********************************************************************/
  92. void delay_1ms(void)
  93. {                          
  94.     unsigned int i;
  95.     for(i = 0;i < 65;i++);
  96. }

  97. void delay_ms(unsigned int n)
  98. {
  99.     unsigned int i;
  100.     for(i = 0;i < n;i++)
  101.     {
  102.       delay_1ms();
  103.     }
  104. }


  105. /********************************************************************/
  106. void disp(void)            //2ms scan 16 segment
  107. {     
  108.      PORTD = 0xff;
  109.      PORTC = 0xff;       

  110.              
  111.         switch(Led_scan_count)
  112.         {
  113.                 case 1:
  114.                         
  115.                  PORTB = led_num[0] ;
  116.                  RD0 = 0;
  117.               break;

  118.                 case 2:
  119.        
  120.                  PORTB = led_num[1];
  121.                  RD1 = 0;  
  122.               break;

  123.                 case 3:

  124.                  PORTB = led_num[2] ;
  125.                  RD2=0;  
  126.               break;

  127.                 case 4:
  128.                
  129.                  PORTB = led_num[3];
  130.                  RD3 = 0;  
  131.               break;

  132.                 case 5:

  133.                  PORTB = led_num[4];
  134.                  RD4 = 0;
  135.               break;

  136.                 case 6:

  137.                  PORTB = led_num[5];
  138.                  RD5 = 0;  
  139.               break;

  140.                 case 7:
  141.                  PORTB = led_num[6];
  142.                  RD6 = 0;   
  143.               break;

  144.                 case 8:
  145.                  PORTB = led_num[7];
  146.                  RD7 = 0;   
  147.               break;

  148.                 case 9:
  149.                  PORTB = led_num[8];
  150.                  RC0 = 0;   
  151.               break;

  152.                 case 10:
  153.                  PORTB = led_num[9];
  154.                  RC1 = 0;   
  155.               break;
  156.                 case 11:
  157.                  PORTB = led_num[10];
  158.                  RC2 = 0;  
  159.               break;

  160.                 case 12:
  161.                  PORTB = led_num[11];
  162.                  RC3 = 0;  
  163.               break;
  164.                 case 13:
  165.                  PORTB = led_num[12];
  166.                  RC4 = 0;
  167.               break;

  168.                 case 14:
  169.                  PORTB = led_num[13];
  170.                  RC5 = 0;  
  171.               break;
  172.                 case 15:
  173.                  PORTB = led_num[14];
  174.                  RC6 = 0;   
  175.               break;

  176.                 case 16:
  177.                  PORTB = led_num[15];
  178.                  RC7 = 0;  
  179.               break;

  180.            }
  181. }


  182. /********************************************************************/
  183. void FFT(void)
  184. {
  185.   unsigned char i,b,j,k,p;
  186.   int Temp_Real,Temp_Imag,temp;

  187.     for( i = 1; i <= LOG_NUMBER; i++)                            /* for(1) */
  188.     {
  189.         b = 1;
  190.         b <<= (i-1);                                       
  191.         for( j = 0; j <= b-1; j++)                              /* for (2) */
  192.         {
  193.             p = 1;
  194.             p <<= (LOG_NUMBER-i);            
  195.             p = p*j;
  196.             for( k = j; k < 32; k = k + 2*b)                /* for (3)  */
  197.             {
  198.                 Temp_Real = Fft_Real[k];
  199.                 Temp_Imag = Fft_Image[k];
  200.                 temp = Fft_Real[k+b];

  201.                 Fft_Real[k] = Fft_Real[k] + ((Fft_Real[k+b]*COS_TAB[p])>>7) + ((Fft_Image[k+b]*SIN_TAB[p])>>7);
  202.                 Fft_Image[k] = Fft_Image[k] - ((Fft_Real[k+b]*SIN_TAB[p])>>7) + ((Fft_Image[k+b]*COS_TAB[p])>>7);

  203.                 Fft_Real[k+b] = Temp_Real - ((Fft_Real[k+b]*COS_TAB[p])>>7) - ((Fft_Image[k+b]*SIN_TAB[p])>>7);
  204.                 Fft_Image[k+b] = Temp_Imag + ((temp*SIN_TAB[p])>>7) - ((Fft_Image[k+b]*COS_TAB[p])>>7);     
  205.                                                                                                                      
  206.                 Fft_Real[k] >>= 1;            
  207.                   Fft_Image[k] >>= 1;
  208.                     Fft_Real[k+b] >>= 1;                 
  209.                       Fft_Image[k+b] >>= 1;
  210.                                                                               
  211.             }     
  212.         }
  213.     }
  214.    
  215.     for( j=0; j<16; j++)//16 frequency division  
  216.     {   

  217.             FFT_value=((((Fft_Real[j]* Fft_Real[j]))+((Fft_Image[j]*Fft_Image[j]))));
  218.             if(FFT_value < 4)
  219.                   FFT_value=0;
  220.             else if(FFT_value < 16)
  221.                   FFT_value=1;
  222.             else if(FFT_value < 36)
  223.                   FFT_value=2;
  224.             else if(FFT_value < 64)
  225.                   FFT_value=3;     
  226.             else if(FFT_value < 100)
  227.                   FFT_value=4;
  228.             else if(FFT_value < 144)
  229.                   FFT_value=5;
  230.             else if(FFT_value < 196)
  231.                   FFT_value=6;
  232.             else if(FFT_value < 256)
  233.                   FFT_value=7;
  234.             else FFT_value=7;           
  235.             /*
  236.             if(FFT_value < 9)
  237.                   FFT_value=0;
  238.             else if(FFT_value < 25)
  239.                   FFT_value=1;
  240.             else if(FFT_value < 49)
  241.                   FFT_value=2;
  242.             else if(FFT_value < 81)
  243.                   FFT_value=3;     
  244.             else if(FFT_value < 121)
  245.                   FFT_value=4;
  246.             else if(FFT_value < 169)
  247.                   FFT_value=5;
  248.             else if(FFT_value < 225)
  249.                   FFT_value=6;
  250.             else if(FFT_value < 289)
  251.                   FFT_value=7;
  252.             else FFT_value=7;
  253.             */
  254.                   
  255.             if(FFT_value > (MAIN_TAB[j]))
  256.             {
  257.                 MAIN_TAB[j] = FFT_value;
  258.             }   
  259.       
  260.         if(FFT_value > (PEAK_TAB[j]))
  261.         {  
  262.                  PEAK_TAB[j] = FFT_value;                
  263.              Peak_speed[j] = 12;    //peak keep time                                       
  264.         }                              
  265.     }
  266.               
  267. }                                                 

  268. void FFT_Process(void)
  269. {
  270.         unsigned char k;
  271.             if(AD_getonce_flag == 1)
  272.             {
  273.                      TMR2IE = 0;                        
  274.                      AD_getonce_flag = 0;
  275.                      
  276.              GODONE = 1;
  277.              while( GODONE == 1);
  278.              //AD_VALUE = (((unsigned int)ADRESH<<8)+ADRESL);   
  279.              Fft_Real[BRTable[AD_count]] = (((unsigned int)ADRESH<<8)+ADRESL);
  280.                         
  281.              AD_count++;
  282.              if(AD_count >= 32)
  283.              {
  284.                  AD_count= 0 ;                           
  285.                  AD_readed_flag = 1;                                
  286.              }
  287.              else
  288.              {
  289.                  TMR2 = 0x00;   
  290.                  TMR2IE = 1;                                            
  291.                  }
  292.              }            
  293.          if(AD_readed_flag == 1)
  294.          {
  295.              for(k = 0;k < 32;k++)
  296.              {
  297.                  Fft_Image[k] = 0;
  298.              }
  299.              FFT();
  300.              AD_readed_flag = 0;
  301.               
  302.              TMR2 = 0x00;   
  303.              TMR2IE = 1;                                    
  304.          }                  
  305. }



  306. /*******************************************************************************************************/
  307. void TMR0_Init(void)
  308. {   
  309.     OSCCON = 0B01110111;  // set timer0 as 8MS
  310.                           // "NC"----IRCF2----IRCF1----IRCF0----OSTS----HTS----LTS----SCS      
  311.                           // IRCF<2:0>   dis                          SCS = 1 --->inside

  312.         TMR0 = 256 - 250;          // 0.5us *250 *16 = 2ms
  313.                           // /RBPU  INTEDG  T0CS  T0SE  PSA  PS2  PS1  PS0                          
  314.     OPTION = 0B00000011;  // timer0,1:16   *16 = 2ms     
  315.                           
  316.     INTCON |= 0B10100000; //INT register:    GIE----PEIE----T0IE----INTE----RBIE----T0IF----INTF----RBIF                                                                                                                                              
  317. }

  318. /********************************************************************/
  319. void ADC_init(void)
  320. {       
  321.          TRISA0  =1;          //set it to input
  322.      //ANS0 = 1;
  323.    
  324.          ANSEL = 0B01110001;  //set AN0 to analog;

  325.          ADCON0 = 0B10000001; //fosc/2,AN0 on,adc enable
  326.                              // adcs1 adcs0 chs3 chs2 chs1 chs0 go/done adon
  327.      //ADCS1 = 1;   //INIT 500k/8 62KHz/32point
  328.      //ADCS0 = 0;


  329.          //ADCON1 = 0B10000000; //vdd & vss as vref   
  330.      VCFG0 = 0;  //VSS
  331.      VCFG1 = 0;  //VDD   
  332.      ADFM = 1;    // AD result right justified
  333.      //ADON = 1;     
  334.      //ADIE = 1;
  335.      //GODONE = 1;
  336.      
  337.      PEIE = 1;      
  338.      GIE = 1;               

  339. }

  340. /*
  341. void TMR1_Init(void)
  342. {   
  343.     //T1CON = 0X01;   //T1GINV  TMR1GE T1CKPS1 T1CKPS0 T1OSCEN /T1SYNC TMR1CS TMR1ON
  344.    
  345.     TMR1CS = 0;        //inside Fosc / 4
  346.     TMR1ON = 1;       //enable timer1
  347.     T1CKPS0 = 0;      //1:1
  348.     T1CKPS0 = 0;
  349.    
  350.    
  351.         TMR1IE = 1;
  352.         TMR1IF = 0;

  353.         TMR1H = 0XFF;
  354.         TMR1L = 0XE7;    //65535-25= 0XFFe7   0.025MS  4M/4
  355.        
  356.         PEIE = 1;
  357.         GIE = 1;
  358.        
  359. }
  360. */

  361. void PORT_Init(void)
  362. {      
  363.     TRISB = 0x00;         
  364.     TRISC = 0x00;
  365.     TRISD = 0x00;
  366.     TRISE = 0X00;
  367.       
  368. }
  369. void TMR2_Init(void)
  370. {
  371.     T2CON = 0x00;   //postscale 1:1
  372.     TMR2ON=1;
  373.    
  374.     PR2 = 0xff;       //30*0.25us +step = 25us
  375.    
  376.     TMR2 = 256-36;         
  377.     TMR2IE = 1;
  378.     TMR2IF = 0;
  379.    
  380.     //PEIE = 1;      
  381.     //GIE = 1;        

  382. }
  383.             
  384. void init(void)
  385. {
  386.     PORT_Init();
  387.     TMR0_Init();
  388.         ADC_init();
  389.     //TMR1_Init();
  390.     TMR2_Init();
  391. }

  392.    
  393. void interrupt ISR(void)
  394. {
  395.     unsigned char j;
  396.       
  397.     if(TMR2IF && TMR2IE)      //25us       
  398.     {     
  399.             TMR2IF = 0;
  400.             TMR2 = 256-36;
  401.             AD_getonce_flag = 1;
  402.                         
  403.     }
  404.   
  405.    
  406.     if(T0IF && T0IE)    //2ms
  407.     {
  408.         TMR0= 256 - 250;
  409.         T0IF=0;

  410.         Led_scan_count++;        //2ms
  411.             if(Led_scan_count > 16)
  412. ……………………

  413. …………限于本文篇幅 余下代码请从51黑下载附件…………
复制代码

所有资料51hei提供下载:
PIC_FFT V0.5.rar (166.49 KB, 下载次数: 109)


评分

参与人数 1黑币 +50 收起 理由
admin + 50 共享资料的黑币奖励!

查看全部评分

分享到:  QQ好友和群QQ好友和群 QQ空间QQ空间 腾讯微博腾讯微博 腾讯朋友腾讯朋友
收藏收藏4 分享淘帖 顶 踩
回复

使用道具 举报

沙发
ID:679307 发表于 2020-1-3 17:01 | 只看该作者
真是高手,学习了
回复

使用道具 举报

板凳
ID:745514 发表于 2020-6-13 00:49 来自手机 | 只看该作者
楼主能做一个用89c51实现音乐频谱的吗
回复

使用道具 举报

地板
ID:767423 发表于 2020-7-9 21:37 | 只看该作者
单片机真好玩呀 发表于 2020-6-13 00:49
楼主能做一个用89c51实现音乐频谱的吗
  1. #include<STC12C5A60S2.H>
  2. typedef unsigned char uint8;
  3. typedef unsigned int  uint16;
  4. #define  channel  0x01                          //设置AD通道为 P1.0
  5. #define ADC_POWER   (1<<7)
  6. #define ADC_SPEEDHH (0x03<<5)
  7. #define ADC_START   (1<<3)
  8. #define PADCH       (1<<5)
  9. #define ADC_FLAG    (1<<4)
  10. #define P_3 P3 //32列扫描控制
  11. sbit sw=P1^2;//74HC154组切换控制位
  12. //放大128倍后的sin整数表(128)        
  13. code char SIN_TAB[128]={0,6,12,18,24,30,36,42,48,54,59,65,70,75,80,85,89,94,98,102,
  14.                         105,108,112,114,117,119,121,123,124,125,126,126,126,126,126,
  15.                                                 125,124,123,121,119,117,114,112,108,105,102,98,94,89,85,80,75,
  16.                                                 70,65,59,54,48,42,36,30,24,18,12,6,0,-6,-12,-18,-24,-30,-36,-42,
  17.                                                 -48,-54,-59,-65,-70,-75,-80,-85,-89,-94,-98,-102,-105,-108,-112,
  18.                                                 -114,-117,-119,-121,-123,-124,-125,-126,-126,-126,-126,-126,-125,
  19.                                                 -124,-123,-121,-119,-117,-114,-112,-108,-105,-102,-98,-94,-89,-85,
  20.                                                 -80,-75,-70,-65,-59,-54,-48,-42,-36,-30,-24,-18,-12,-6};

  21. //放大128倍后的cos整数表(128)
  22. code char COS_TAB[128]={127,126,126,125,124,123,121,119,117,114,112,108,105,102,98,94,
  23.                         89,85,80,75,70,65,59,54,48,42,36,30,24,18,12,6,0,-6,-12,-18,-24,
  24.                                                 -30,-36,-42,-48,-54,-59,-65,-70,-75,-80,-85,-89,-94,-98,-102,
  25.                                                 -105,-108,-112,-114,-117,-119,-121,-123,-124,-125,-126,-126,-126,
  26.                                                 -126,-126,-125,-124,-123,-121,-119,-117,-114,-112,-108,-105,-102,
  27.                                                 -98,-94,-89,-85,-80,-75,-70,-65,-59,-54,-48,-42,-36,-30,-24,-18,
  28.                                                 -12,-6,0,6,12,18,24,30,36,42,48,54,59,65,70,75,80,85,89,94,98,102,
  29.                                                 105,108,112,114,117,119,121,123,124,125,126,126};

  30. //采样存储序列表
  31. code char LIST_TAB[128] = {0,64,32,96,16,80,48,112,8,72,40,104,24,88,56,120,4,68,36,100,20,84,52,116,
  32.                            12,76,44,108,28,92,60,124,2,66,34,98,18,82,50,114,10,74,42,106,26,90,58,
  33.                                                    122,6,70,38,102,22,86,54,118,14,78,46,110,30,94,62,126,1,65,33,97,17,81,49,
  34.                                                    113,9,73,41,105,25,89,57,121,5,69,37,101,21,85,53,117,13,77,45,109,29,93,61,
  35.                                                    125,3,67,35,99,19,83,51,115,11,75,43,107,27,91,59,123,7,71,39,103,23,87,55,
  36.                                                    119,15,79,47,111,31,95,63,127};

  37. uint8 LED[32]={0x80,0x00,0xC0,0x00,0xE0,0x00,0xF0,0x00,
  38.                0xF8,0x00,0xFC,0x00,0xFE,0x00,0xFF,0x00,
  39.                0xff,0x80,0xff,0xC0,0xff,0xE0,0xff,0xF0,
  40.                            0xff,0xF8,0xff,0xFC,0xff,0xFE,0xff,0xFF,};
  41. uint8 COUNT=15,COUNT1=0,ADC_Count=0,LINE=15,G=0,T;
  42. uint8 i,j,k,b,p,anum,led=0;                 
  43. int Temp_Real,Temp_Imag,temp;                // 中间临时变量  
  44. uint16 TEMP1,max=0;        
  45. int xdata Fft_Real[128];
  46. int xdata Fft_Image[128];               // fft的虚部
  47. uint8 xdata LED_TAB2[16];                                //记录 漂浮物 是否需要 停顿一下
  48. uint8 xdata LED_TAB[16];                                //记录红色柱状
  49. uint8 xdata LED_TAB1[16];                                //记录 漂浮点
  50. /////////////shop60610364.taobao.com/////////
  51. void FFT()//基2fft
  52. {                 
  53.     for( i=1; i<=7; i++)                            /* for(1) */
  54.     {
  55.         b=1;
  56.         b <<=(i-1);                                 //碟式运算,用于计算隔多少行计算例如 第一极 1和2行计算
  57.         for( j=0; j<=b-1; j++)                      /* for (2) */
  58.         {
  59.             p=1;
  60.             p <<= (7-i);            
  61.             p = p*j;
  62.             for( k=j; k<128; k=k+2*b)               /* for (3)  */
  63.             {
  64.                 Temp_Real=Fft_Real[k];
  65.                                 Temp_Imag=Fft_Image[k];
  66.                                 temp=Fft_Real[k+b];
  67.                 Fft_Real[k]=Fft_Real[k]+((Fft_Real[k+b]*COS_TAB[p])>>7)+((Fft_Image[k+b]*SIN_TAB[p])>>7);
  68.                 Fft_Image[k]=Fft_Image[k]-((Fft_Real[k+b]*SIN_TAB[p])>>7)+((Fft_Image[k+b]*COS_TAB[p])>>7);
  69.                 Fft_Real[k+b]=Temp_Real-((Fft_Real[k+b]*COS_TAB[p])>>7)-((Fft_Image[k+b]*SIN_TAB[p])>>7);
  70.                 Fft_Image[k+b]=Temp_Imag+((temp*SIN_TAB[p])>>7)-((Fft_Image[k+b]*COS_TAB[p])>>7);          // 移位.防止溢出. 结果已经是本值的 1/64               
  71.                 Fft_Real[k] >>= 1;            
  72.                 Fft_Image[k] >>= 1;
  73.                 Fft_Real[k+b] >>= 1;                 
  74.                 Fft_Image[k+b] >>= 1;
  75.                                                                               
  76.             }     
  77.         }
  78.     }

  79.         for(j=0;j<16;j++)//16分频                                                                                         
  80.         {                                                                                                                                                         
  81.                 TEMP1=((((Fft_Real[j+1]* Fft_Real[j+1]))+((Fft_Image[j+1]*Fft_Image[j+1])))>>1);//求各频段幅值
  82.                 if(TEMP1<1)TEMP1=0;
  83.                 LED_TAB[j]=TEMP1;        
  84.                 if(LED_TAB[j]>max)max=LED_TAB[j];
  85.         }
  86.         if(max>16)//分级量化
  87.         {
  88.           max/=16;
  89.           for(j=0;j<16;j++)LED_TAB[j]/=max;
  90.         }      
  91. }


  92. void delay2us(void)   //2us延时
  93. {
  94.     unsigned char a,b;
  95.     for(b=1;b>0;b--)
  96.         for(a=3;a>0;a--);
  97. }



  98. void Init()
  99. {  
  100.         
  101. //-----------------------------------------------------------------------------------

  102.      P1ASF = 0x02;                 //0000,0010, 将 P1.1 置成模拟口
  103.      AUXR1 &=0xFB;                 //1111,1011, 令 ADRJ=0
  104.          EADC=1;                                           //AD中断打开
  105.          ADC_CONTR = ADC_POWER | ADC_SPEEDHH | ADC_START | channel;//1110 1001   1打开 A/D (ADC_POWER)转换电源;11速度为90周期一次;
  106.                                                                            //0中断标志清零;1启动adc(ADC_START);001AD通道打开(这里为P1.1);
  107. //-----------------------------------------------------------------------------------
  108.          
  109.          P2M0=1;
  110.      TMOD=0X12;                                          
  111.          TH0=0xb0;                                                                  
  112.          TL0=0xb0;
  113.          TH1=0xfe;                                                                                                   
  114.          TL1=0Xd4;
  115.          ET0=1;                                               //定时器0 打开
  116.          TR0=0;                                                      //关闭定时器
  117.          ET1=1;
  118.          TR1=1;
  119.          PT1=0;
  120.          PT0=1;
  121.      IPH=PADCH;
  122.          IP=PADC;                                           //中断优先级
  123.          EA=1;                                                   //总中断打开
  124. }

  125. //----------------------------------------------------------------------------------------

  126. void ADC_Finish() interrupt 5
  127. {           ADC_CONTR &= !ADC_FLAG;
  128.             Fft_Real[LIST_TAB[ADC_Count]]=(int)((ADC_RES)<<1)+(ADC_RESL>>1)*2;//按LIST_TAB表里的顺序,进行存储采样值
  129.           if(ADC_Count<=127)ADC_Count++;
  130.           else {EADC=0;TR0=0;}                                                                                                                  
  131. }
  132. //---------------------------------------------------------------------------------------
  133. void LED_Display() interrupt 3                                         //中断一次显示一列
  134. {                  
  135.          TH1=0xfe;                                                                                                
  136.         TL1=0Xd4;                                                
  137.     G++;                                                         
  138.     if(G==33)G=1;
  139.         switch(G)                                                                     //往点阵屏填充一列的数据
  140.           {
  141.                 case 1:sw=1;anum=LED_TAB[0];P_3=0x00;break;
  142.                 case 2:anum=LED_TAB[0];P_3=0x01;break;
  143.                 case 3:anum=LED_TAB[1];P_3=0x02;break;
  144.                 case 4:anum=LED_TAB[1];P_3=0x03;break;
  145.                 case 5:anum=LED_TAB[2];P_3=0x04;break;
  146.                 case 6:anum=LED_TAB[2];P_3=0x05;break;
  147.                 case 7:anum=LED_TAB[3];P_3=0x06;break;
  148.                 case 8:anum=LED_TAB[3];P_3=0x07;break;
  149.                 case 9:anum=LED_TAB[4];P_3=0x08;break;
  150.                 case 10:anum=LED_TAB[4];P_3=0x09;break;
  151.                 case 11:anum=LED_TAB[5];P_3=0x0a;break;
  152.                 case 12:anum=LED_TAB[5];P_3=0x0b;break;
  153.                 case 13:anum=LED_TAB[6];P_3=0x0c;break;
  154.                 case 14:anum=LED_TAB[6];P_3=0x0d;break;
  155.                 case 15:anum=LED_TAB[7];P_3=0x0e;break;
  156.                 case 16:anum=LED_TAB[7];P_3=0x0f;break;
  157.                 case 17:sw=0;anum=LED_TAB[8];P_3=0x00;break;
  158.                 case 18:anum=LED_TAB[8];P_3=0x10;break;
  159.                 case 19:anum=LED_TAB[9];P_3=0x20;break;
  160.                 case 20:anum=LED_TAB[9];P_3=0x30;break;
  161.                 case 21:anum=LED_TAB[10];P_3=0x40;break;
  162.                 case 22:anum=LED_TAB[10];P_3=0x50;break;
  163.                 case 23:anum=LED_TAB[11];P_3=0x60;break;
  164.                 case 24:anum=LED_TAB[11];P_3=0x70;break;
  165.                 case 25:anum=LED_TAB[12];P_3=0x80;break;
  166.                 case 26:anum=LED_TAB[12];P_3=0x90;break;
  167.                 case 27:anum=LED_TAB[13];P_3=0xa0;break;
  168.                 case 28:anum=LED_TAB[13];P_3=0xb0;break;
  169.                 case 29:anum=LED_TAB[14];P_3=0xc0;break;
  170.                 case 30:anum=LED_TAB[14];P_3=0xd0;break;
  171.                 case 31:anum=LED_TAB[15];P_3=0xe0;break;
  172.                 case 32:anum=LED_TAB[15];P_3=0xf0;break;
  173.                     
  174.                 }

  175.            for(led=0;led<=anum;led++)                                 //递增
  176.              { P0=LED[2*led+1];
  177.                    P2=LED[2*led];
  178.                    delay2us();
  179.                   
  180.                  }                                                                                 //递减
  181.            for(;led>=0;led--)
  182.              { P0=LED[2*led+1];
  183.                    P2=LED[2*led];
  184.                    delay2us();
  185.                  }
  186. }
  187. //------------------------------------------------------------------------------------------------------------

  188. void Ad_Control() interrupt 1                                         //控制采样率
  189. {
  190. ADC_CONTR = ADC_POWER | ADC_SPEEDHH| ADC_START | channel;         //开始AD采集
  191. }
  192.                                                                           
  193. //===============================================================================================================

  194. void main()
  195. {
  196.     uint8 i;
  197.          Init();
  198.         while(1)
  199.         {               
  200.                 ADC_Count=0;
  201.                 TR0=1;
  202.                 EADC=1;        //开启定时器中断0,开启ADC
  203.                 while(ADC_Count<128);
  204.         for(i=0;i<128;i++)
  205.         {
  206.          Fft_Image[i]=0;
  207.         }
  208.                 FFT();//FFT运算并转换为各频段幅值                                                   
  209.                  
  210.         }
  211. }
复制代码
回复

使用道具 举报

5#
ID:801012 发表于 2020-7-15 11:10 | 只看该作者

采样数据存储列表内数据是怎么计算的呢?
回复

使用道具 举报

6#
ID:512651 发表于 2020-11-16 17:58 | 只看该作者
感谢楼主共享学习资料!!!
回复

使用道具 举报

7#
ID:651494 发表于 2022-12-20 15:14 | 只看该作者
王tt 发表于 2020-7-15 11:10
采样数据存储列表内数据是怎么计算的呢?

把所有数据弄到EXECL,做个图形就明白了
回复

使用道具 举报

8#
ID:1066888 发表于 2023-3-16 21:44 | 只看该作者
FFT自己编写的C,不错的内容,学习一下
回复

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

手机版|小黑屋|51黑电子论坛 |51黑电子论坛6群 QQ 管理员QQ:125739409;技术交流QQ群281945664

Powered by 单片机教程网

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