標(biāo)題: PIC單片機(jī)FFT音樂頻譜Proteus仿真程序 [打印本頁]
作者: feeling1791    時間: 2020-1-3 11:50
標(biāo)題: PIC單片機(jī)FFT音樂頻譜Proteus仿真程序
因為貼子附件大小問題，wav音樂文件沒有附上，請自行添加到目錄，打開FFT.pdsprj文件，然后雙擊Audio圖標(biāo)指定音樂wav位置，開始運(yùn)行即可，稍后再發(fā)下電腦里以前做的東東


單片機(jī)源程序如下:
  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);  

  64. #define LOG_NUMBER 5          

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

  81. unsigned char disp_time;
  82. unsigned char Led_scan_count;

  84. bit AD_readed_flag = 0;
  85. bit AD_getonce_flag = 0;

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

  91. bank2 unsigned char MAIN_TAB[16];
  92. bank1 unsigned char PEAK_TAB[16];
  93. bank3 unsigned char Peak_speed[16];

  95. bank1 int Fft_Real[32];
  96. bank2 int Fft_Image[32];

  98. unsigned int FFT_value;

  100. /********************************************************************/
  101. void delay_1ms(void)
  102. {                          
  103.     unsigned int i;
  104.     for(i = 0;i < 65;i++);
  105. }

  107. void delay_ms(unsigned int n)
  108. {
  109.     unsigned int i;
  110.     for(i = 0;i < n;i++)
  111.     {
  112.       delay_1ms();
  113.     }
  114. }


  117. /********************************************************************/
  118. void disp(void)            //2ms scan 16 segment
  119. {     
  120.      PORTD = 0xff;
  121.      PORTC = 0xff;       

  123.              
  124.         switch(Led_scan_count)
  125.         {
  126.                 case 1:
  127.                         
  128.                  PORTB = led_num[0] ;
  129.                  RD0 = 0;
  130.               break;

  132.                 case 2:
  133.        
  134.                  PORTB = led_num[1];
  135.                  RD1 = 0;  
  136.               break;

  138.                 case 3:

  140.                  PORTB = led_num[2] ;
  141.                  RD2=0;  
  142.               break;

  144.                 case 4:
  145.                
  146.                  PORTB = led_num[3];
  147.                  RD3 = 0;  
  148.               break;

  150.                 case 5:

  152.                  PORTB = led_num[4];
  153.                  RD4 = 0;
  154.               break;

  156.                 case 6:

  158.                  PORTB = led_num[5];
  159.                  RD5 = 0;  
  160.               break;

  162.                 case 7:
  163.                  PORTB = led_num[6];
  164.                  RD6 = 0;   
  165.               break;

  167.                 case 8:
  168.                  PORTB = led_num[7];
  169.                  RD7 = 0;   
  170.               break;

  172.                 case 9:
  173.                  PORTB = led_num[8];
  174.                  RC0 = 0;   
  175.               break;

  177.                 case 10:
  178.                  PORTB = led_num[9];
  179.                  RC1 = 0;   
  180.               break;
  181.                 case 11:
  182.                  PORTB = led_num[10];
  183.                  RC2 = 0;  
  184.               break;

  186.                 case 12:
  187.                  PORTB = led_num[11];
  188.                  RC3 = 0;  
  189.               break;
  190.                 case 13:
  191.                  PORTB = led_num[12];
  192.                  RC4 = 0;
  193.               break;

  195.                 case 14:
  196.                  PORTB = led_num[13];
  197.                  RC5 = 0;  
  198.               break;
  199.                 case 15:
  200.                  PORTB = led_num[14];
  201.                  RC6 = 0;   
  202.               break;

  204.                 case 16:
  205.                  PORTB = led_num[15];
  206.                  RC7 = 0;  
  207.               break;

  209.            }
  210. }


  213. /********************************************************************/
  214. void FFT(void)
  215. {
  216.   unsigned char i,b,j,k,p;
  217.   int Temp_Real,Temp_Imag,temp;

  219.     for( i = 1; i <= LOG_NUMBER; i++)                            /* for(1) */
  220.     {
  221.         b = 1;
  222.         b <<= (i-1);                                       
  223.         for( j = 0; j <= b-1; j++)                              /* for (2) */
  224.         {
  225.             p = 1;
  226.             p <<= (LOG_NUMBER-i);            
  227.             p = p*j;
  228.             for( k = j; k < 32; k = k + 2*b)                /* for (3)  */
  229.             {
  230.                 Temp_Real = Fft_Real[k];
  231.                 Temp_Imag = Fft_Image[k];
  232.                 temp = Fft_Real[k+b];

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

  237.                 Fft_Real[k+b] = Temp_Real - ((Fft_Real[k+b]*COS_TAB[p])>>7) - ((Fft_Image[k+b]*SIN_TAB[p])>>7);
  238.                 Fft_Image[k+b] = Temp_Imag + ((temp*SIN_TAB[p])>>7) - ((Fft_Image[k+b]*COS_TAB[p])>>7);     
  239.                                                                                                                      
  240.                 Fft_Real[k] >>= 1;            
  241.                   Fft_Image[k] >>= 1;
  242.                     Fft_Real[k+b] >>= 1;                 
  243.                       Fft_Image[k+b] >>= 1;
  244.                                                                               
  245.             }     
  246.         }
  247.     }
  248.    
  249.     for( j=0; j<16; j++)//16 frequency division  
  250.     {   

  252.             FFT_value=((((Fft_Real[j]* Fft_Real[j]))+((Fft_Image[j]*Fft_Image[j]))));
  253.             if(FFT_value < 4)
  254.                   FFT_value=0;
  255.             else if(FFT_value < 16)
  256.                   FFT_value=1;
  257.             else if(FFT_value < 36)
  258.                   FFT_value=2;
  259.             else if(FFT_value < 64)
  260.                   FFT_value=3;     
  261.             else if(FFT_value < 100)
  262.                   FFT_value=4;
  263.             else if(FFT_value < 144)
  264.                   FFT_value=5;
  265.             else if(FFT_value < 196)
  266.                   FFT_value=6;
  267.             else if(FFT_value < 256)
  268.                   FFT_value=7;
  269.             else FFT_value=7;           
  270.             /*
  271.             if(FFT_value < 9)
  272.                   FFT_value=0;
  273.             else if(FFT_value < 25)
  274.                   FFT_value=1;
  275.             else if(FFT_value < 49)
  276.                   FFT_value=2;
  277.             else if(FFT_value < 81)
  278.                   FFT_value=3;     
  279.             else if(FFT_value < 121)
  280.                   FFT_value=4;
  281.             else if(FFT_value < 169)
  282.                   FFT_value=5;
  283.             else if(FFT_value < 225)
  284.                   FFT_value=6;
  285.             else if(FFT_value < 289)
  286.                   FFT_value=7;
  287.             else FFT_value=7;
  288.             */
  289.                   
  290.             if(FFT_value > (MAIN_TAB[j]))
  291.             {
  292.                 MAIN_TAB[j] = FFT_value;
  293.             }   
  294.       
  295.         if(FFT_value > (PEAK_TAB[j]))
  296.         {  
  297.                  PEAK_TAB[j] = FFT_value;                
  298.              Peak_speed[j] = 12;    //peak keep time                                       
  299.         }                              
  300.     }
  301.               
  302. }                                                 

  304. void FFT_Process(void)
  305. {
  306.         unsigned char k;
  307.             if(AD_getonce_flag == 1)
  308.             {
  309.                      TMR2IE = 0;                        
  310.                      AD_getonce_flag = 0;
  311.                      
  312.              GODONE = 1;
  313.              while( GODONE == 1);
  314.              //AD_VALUE = (((unsigned int)ADRESH<<8)+ADRESL);   
  315.              Fft_Real[BRTable[AD_count]] = (((unsigned int)ADRESH<<8)+ADRESL);
  316.                         
  317.              AD_count++;
  318.              if(AD_count >= 32)
  319.              {
  320.                  AD_count= 0 ;                           
  321.                  AD_readed_flag = 1;                                
  322.              }
  323.              else
  324.              {
  325.                  TMR2 = 0x00;   
  326.                  TMR2IE = 1;                                            
  327.                  }
  328.              }            
  329.          if(AD_readed_flag == 1)
  330.          {
  331.              for(k = 0;k < 32;k++)
  332.              {
  333.                  Fft_Image[k] = 0;
  334.              }
  335.              FFT();
  336.              AD_readed_flag = 0;
  337.               
  338.              TMR2 = 0x00;   
  339.              TMR2IE = 1;                                    
  340.          }                  
  341. }



  345. /*******************************************************************************************************/
  346. void TMR0_Init(void)
  347. {   
  348.     OSCCON = 0B01110111;  // set timer0 as 8MS
  349.                           // "NC"----IRCF2----IRCF1----IRCF0----OSTS----HTS----LTS----SCS      
  350.                           // IRCF<2:0>   dis                          SCS = 1 --->inside

  352.         TMR0 = 256 - 250;          // 0.5us *250 *16 = 2ms
  353.                           // /RBPU  INTEDG  T0CS  T0SE  PSA  PS2  PS1  PS0                          
  354.     OPTION = 0B00000011;  // timer0,1:16   *16 = 2ms     
  355.                           
  356.     INTCON |= 0B10100000; //INT register:    GIE----PEIE----T0IE----INTE----RBIE----T0IF----INTF----RBIF                                                                                                                                              
  357. }

  359. /********************************************************************/
  360. void ADC_init(void)
  361. {       
  362.          TRISA0  =1;          //set it to input
  363.      //ANS0 = 1;
  364.    
  365.          ANSEL = 0B01110001;  //set AN0 to analog;

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


  373.          //ADCON1 = 0B10000000; //vdd & vss as vref   
  374.      VCFG0 = 0;  //VSS
  375.      VCFG1 = 0;  //VDD   
  376.      ADFM = 1;    // AD result right justified
  377.      //ADON = 1;     
  378.      //ADIE = 1;
  379.      //GODONE = 1;
  380.      
  381.      PEIE = 1;      
  382.      GIE = 1;               

  384. }

  386. /*
  387. void TMR1_Init(void)
  388. {   
  389.     //T1CON = 0X01;   //T1GINV  TMR1GE T1CKPS1 T1CKPS0 T1OSCEN /T1SYNC TMR1CS TMR1ON
  390.    
  391.     TMR1CS = 0;        //inside Fosc / 4
  392.     TMR1ON = 1;       //enable timer1
  393.     T1CKPS0 = 0;      //1:1
  394.     T1CKPS0 = 0;
  395.    
  396.    
  397.         TMR1IE = 1;
  398.         TMR1IF = 0;

  400.         TMR1H = 0XFF;
  401.         TMR1L = 0XE7;    //65535-25= 0XFFe7   0.025MS  4M/4
  402.        
  403.         PEIE = 1;
  404.         GIE = 1;
  405.        
  406. }
  407. */

  409. void PORT_Init(void)
  410. {      
  411.     TRISB = 0x00;         
  412.     TRISC = 0x00;
  413.     TRISD = 0x00;
  414.     TRISE = 0X00;
  415.       
  416. }
  417. void TMR2_Init(void)
  418. {
  419.     T2CON = 0x00;   //postscale 1:1
  420.     TMR2ON=1;
  421.    
  422.     PR2 = 0xff;       //30*0.25us +step = 25us
  423.    
  424.     TMR2 = 256-36;         
  425.     TMR2IE = 1;
  426.     TMR2IF = 0;
  427.    
  428.     //PEIE = 1;      
  429.     //GIE = 1;        

  431. }
  432.             
  433. void init(void)
  434. {
  435.     PORT_Init();
  436.     TMR0_Init();
  437.         ADC_init();
  438.     //TMR1_Init();
  439.     TMR2_Init();
  440. }

  442.    
  443. void interrupt ISR(void)
  444. {
  445.     unsigned char j;
  446.       
  447.     if(TMR2IF && TMR2IE)      //25us       
  448.     {     
  449.             TMR2IF = 0;
  450.             TMR2 = 256-36;
  451.             AD_getonce_flag = 1;
  452.                         
  453.     }
  454.   
  455.    
  456.     if(T0IF && T0IE)    //2ms
  457.     {
  458.         TMR0= 256 - 250;
  459.         T0IF=0;

  461.         Led_scan_count++;        //2ms
  462.             if(Led_scan_count > 16)
  463. ……………………

  465. …………限于本文篇幅 余下代碼請從51黑下載附件…………
復(fù)制代碼

所有資料51hei提供下載:
PIC_FFT V0.5.rar (166.49 KB, 下載次數(shù): 112)


作者: zwwhi    時間: 2020-1-3 17:01
真是高手，學(xué)習(xí)了
作者: 單片機(jī)真好玩呀    時間: 2020-6-13 00:49
樓主能做一個用89c51實現(xiàn)音樂頻譜的嗎
作者: hsl32    時間: 2020-7-9 21:37
單片機(jī)真好玩呀 發(fā)表于 2020-6-13 00:49
樓主能做一個用89c51實現(xiàn)音樂頻譜的嗎
  1. #include<STC12C5A60S2.H>
  2. typedef unsigned char uint8;
  3. typedef unsigned int  uint16;
  4. #define  channel  0x01                          //設(shè)置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整數(shù)表（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};

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

  32. //采樣存儲序列表
  33. 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,
  34.                            12,76,44,108,28,92,60,124,2,66,34,98,18,82,50,114,10,74,42,106,26,90,58,
  35.                                                    122,6,70,38,102,22,86,54,118,14,78,46,110,30,94,62,126,1,65,33,97,17,81,49,
  36.                                                    113,9,73,41,105,25,89,57,121,5,69,37,101,21,85,53,117,13,77,45,109,29,93,61,
  37.                                                    125,3,67,35,99,19,83,51,115,11,75,43,107,27,91,59,123,7,71,39,103,23,87,55,
  38.                                                    119,15,79,47,111,31,95,63,127};

  40. uint8 LED[32]={0x80,0x00,0xC0,0x00,0xE0,0x00,0xF0,0x00,
  41.                0xF8,0x00,0xFC,0x00,0xFE,0x00,0xFF,0x00,
  42.                0xff,0x80,0xff,0xC0,0xff,0xE0,0xff,0xF0,
  43.                            0xff,0xF8,0xff,0xFC,0xff,0xFE,0xff,0xFF,};
  44. uint8 COUNT=15,COUNT1=0,ADC_Count=0,LINE=15,G=0,T;
  45. uint8 i,j,k,b,p,anum,led=0;                 
  46. int Temp_Real,Temp_Imag,temp;                // 中間臨時變量  
  47. uint16 TEMP1,max=0;        
  48. int xdata Fft_Real[128];
  49. int xdata Fft_Image[128];               // fft的虛部
  50. uint8 xdata LED_TAB2[16];                                //記錄 漂浮物 是否需要 停頓一下
  51. uint8 xdata LED_TAB[16];                                //記錄紅色柱狀
  52. uint8 xdata LED_TAB1[16];                                //記錄 漂浮點(diǎn)
  53. /////////////shop60610364.taobao.com/////////
  54. void FFT()//基2fft
  55. {                 
  56.     for( i=1; i<=7; i++)                            /* for(1) */
  57.     {
  58.         b=1;
  59.         b <<=(i-1);                                 //碟式運(yùn)算，用于計算隔多少行計算例如 第一極 1和2行計算
  60.         for( j=0; j<=b-1; j++)                      /* for (2) */
  61.         {
  62.             p=1;
  63.             p <<= (7-i);            
  64.             p = p*j;
  65.             for( k=j; k<128; k=k+2*b)               /* for (3)  */
  66.             {
  67.                 Temp_Real=Fft_Real[k];
  68.                                 Temp_Imag=Fft_Image[k];
  69.                                 temp=Fft_Real[k+b];
  70.                 Fft_Real[k]=Fft_Real[k]+((Fft_Real[k+b]*COS_TAB[p])>>7)+((Fft_Image[k+b]*SIN_TAB[p])>>7);
  71.                 Fft_Image[k]=Fft_Image[k]-((Fft_Real[k+b]*SIN_TAB[p])>>7)+((Fft_Image[k+b]*COS_TAB[p])>>7);
  72.                 Fft_Real[k+b]=Temp_Real-((Fft_Real[k+b]*COS_TAB[p])>>7)-((Fft_Image[k+b]*SIN_TAB[p])>>7);
  73.                 Fft_Image[k+b]=Temp_Imag+((temp*SIN_TAB[p])>>7)-((Fft_Image[k+b]*COS_TAB[p])>>7);          // 移位.防止溢出. 結(jié)果已經(jīng)是本值的 1/64               
  74.                 Fft_Real[k] >>= 1;            
  75.                 Fft_Image[k] >>= 1;
  76.                 Fft_Real[k+b] >>= 1;                 
  77.                 Fft_Image[k+b] >>= 1;
  78.                                                                               
  79.             }     
  80.         }
  81.     }

  83.         for(j=0;j<16;j++)//16分頻                                                                                         
  84.         {                                                                                                                                                         
  85.                 TEMP1=((((Fft_Real[j+1]* Fft_Real[j+1]))+((Fft_Image[j+1]*Fft_Image[j+1])))>>1);//求各頻段幅值
  86.                 if(TEMP1<1)TEMP1=0;
  87.                 LED_TAB[j]=TEMP1;        
  88.                 if(LED_TAB[j]>max)max=LED_TAB[j];
  89.         }
  90.         if(max>16)//分級量化
  91.         {
  92.           max/=16;
  93.           for(j=0;j<16;j++)LED_TAB[j]/=max;
  94.         }      
  95. }


  98. void delay2us(void)   //2us延時
  99. {
  100.     unsigned char a,b;
  101.     for(b=1;b>0;b--)
  102.         for(a=3;a>0;a--);
  103. }



  107. void Init()
  108. {  
  109.         
  110. //-----------------------------------------------------------------------------------

  112.      P1ASF = 0x02;                 //0000,0010, 將 P1.1 置成模擬口
  113.      AUXR1 &=0xFB;                 //1111,1011, 令 ADRJ=0
  114.          EADC=1;                                           //AD中斷打開
  115.          ADC_CONTR = ADC_POWER | ADC_SPEEDHH | ADC_START | channel;//1110 1001   1打開 A/D （ADC_POWER）轉(zhuǎn)換電源；11速度為90周期一次；
  116.                                                                            //0中斷標(biāo)志清零；1啟動adc(ADC_START);001AD通道打開（這里為P1.1）;
  117. //-----------------------------------------------------------------------------------
  118.          
  119.          P2M0=1;
  120.      TMOD=0X12;                                          
  121.          TH0=0xb0;                                                                  
  122.          TL0=0xb0;
  123.          TH1=0xfe;                                                                                                   
  124.          TL1=0Xd4;
  125.          ET0=1;                                               //定時器0 打開
  126.          TR0=0;                                                      //關(guān)閉定時器
  127.          ET1=1;
  128.          TR1=1;
  129.          PT1=0;
  130.          PT0=1;
  131.      IPH=PADCH;
  132.          IP=PADC;                                           //中斷優(yōu)先級
  133.          EA=1;                                                   //總中斷打開
  134. }

  136. //----------------------------------------------------------------------------------------

  138. void ADC_Finish() interrupt 5
  139. {           ADC_CONTR &= !ADC_FLAG;
  140.             Fft_Real[LIST_TAB[ADC_Count]]=(int)((ADC_RES)<<1)+(ADC_RESL>>1)*2;//按LIST_TAB表里的順序，進(jìn)行存儲采樣值
  141.           if(ADC_Count<=127)ADC_Count++;
  142.           else {EADC=0;TR0=0;}                                                                                                                  
  143. }
  144. //---------------------------------------------------------------------------------------
  145. void LED_Display() interrupt 3                                         //中斷一次顯示一列
  146. {                  
  147.          TH1=0xfe;                                                                                                
  148.         TL1=0Xd4;                                                
  149.     G++;                                                         
  150.     if(G==33)G=1;
  151.         switch(G)                                                                     //往點(diǎn)陣屏填充一列的數(shù)據(jù)
  152.           {
  153.                 case 1:sw=1;anum=LED_TAB[0];P_3=0x00;break;
  154.                 case 2:anum=LED_TAB[0];P_3=0x01;break;
  155.                 case 3:anum=LED_TAB[1];P_3=0x02;break;
  156.                 case 4:anum=LED_TAB[1];P_3=0x03;break;
  157.                 case 5:anum=LED_TAB[2];P_3=0x04;break;
  158.                 case 6:anum=LED_TAB[2];P_3=0x05;break;
  159.                 case 7:anum=LED_TAB[3];P_3=0x06;break;
  160.                 case 8:anum=LED_TAB[3];P_3=0x07;break;
  161.                 case 9:anum=LED_TAB[4];P_3=0x08;break;
  162.                 case 10:anum=LED_TAB[4];P_3=0x09;break;
  163.                 case 11:anum=LED_TAB[5];P_3=0x0a;break;
  164.                 case 12:anum=LED_TAB[5];P_3=0x0b;break;
  165.                 case 13:anum=LED_TAB[6];P_3=0x0c;break;
  166.                 case 14:anum=LED_TAB[6];P_3=0x0d;break;
  167.                 case 15:anum=LED_TAB[7];P_3=0x0e;break;
  168.                 case 16:anum=LED_TAB[7];P_3=0x0f;break;
  169.                 case 17:sw=0;anum=LED_TAB[8];P_3=0x00;break;
  170.                 case 18:anum=LED_TAB[8];P_3=0x10;break;
  171.                 case 19:anum=LED_TAB[9];P_3=0x20;break;
  172.                 case 20:anum=LED_TAB[9];P_3=0x30;break;
  173.                 case 21:anum=LED_TAB[10];P_3=0x40;break;
  174.                 case 22:anum=LED_TAB[10];P_3=0x50;break;
  175.                 case 23:anum=LED_TAB[11];P_3=0x60;break;
  176.                 case 24:anum=LED_TAB[11];P_3=0x70;break;
  177.                 case 25:anum=LED_TAB[12];P_3=0x80;break;
  178.                 case 26:anum=LED_TAB[12];P_3=0x90;break;
  179.                 case 27:anum=LED_TAB[13];P_3=0xa0;break;
  180.                 case 28:anum=LED_TAB[13];P_3=0xb0;break;
  181.                 case 29:anum=LED_TAB[14];P_3=0xc0;break;
  182.                 case 30:anum=LED_TAB[14];P_3=0xd0;break;
  183.                 case 31:anum=LED_TAB[15];P_3=0xe0;break;
  184.                 case 32:anum=LED_TAB[15];P_3=0xf0;break;
  185.                     
  186.                 }

  188.            for(led=0;led<=anum;led++)                                 //遞增
  189.              { P0=LED[2*led+1];
  190.                    P2=LED[2*led];
  191.                    delay2us();
  192.                   
  193.                  }                                                                                 //遞減
  194.            for(;led>=0;led--)
  195.              { P0=LED[2*led+1];
  196.                    P2=LED[2*led];
  197.                    delay2us();
  198.                  }
  199. }
  200. //------------------------------------------------------------------------------------------------------------

  202. void Ad_Control() interrupt 1                                         //控制采樣率
  203. {
  204. ADC_CONTR = ADC_POWER | ADC_SPEEDHH| ADC_START | channel;         //開始AD采集
  205. }
  206.                                                                           
  207. //===============================================================================================================

  209. void main()
  210. {
  211.     uint8 i;
  212.          Init();
  213.         while(1)
  214.         {               
  215.                 ADC_Count=0;
  216.                 TR0=1;
  217.                 EADC=1;        //開啟定時器中斷0，開啟ADC
  218.                 while(ADC_Count<128);
  219.         for(i=0;i<128;i++)
  220.         {
  221.          Fft_Image[i]=0;
  222.         }
  223.                 FFT();//FFT運(yùn)算并轉(zhuǎn)換為各頻段幅值                                                   
  224.                  
  225.         }
  226. }
復(fù)制代碼
作者: 王tt    時間: 2020-7-15 11:10
hsl32 發(fā)表于 2020-7-9 21:37

采樣數(shù)據(jù)存儲列表內(nèi)數(shù)據(jù)是怎么計算的呢？
作者: gw168sh    時間: 2020-11-16 17:58
感謝樓主共享學(xué)習(xí)資料�。�！
作者: feeling1791    時間: 2022-12-20 15:14
王tt 發(fā)表于 2020-7-15 11:10
采樣數(shù)據(jù)存儲列表內(nèi)數(shù)據(jù)是怎么計算的呢？

把所有數(shù)據(jù)弄到EXECL，做個圖形就明白了
作者: 許玩童    時間: 2023-3-16 21:44
FFT自己編寫的C，不錯的內(nèi)容，學(xué)習(xí)一下
作者: chief1    時間: 2024-11-27 10:32
audio無法使用怎么辦呀，proteus8.17



歡迎光臨 (http://www.torrancerestoration.com/bbs/) Powered by Discuz! X3.1