标题: 五色LED音乐频谱源程序和原理图 原件清单 [打印本页]
作者: 旅人    时间: 2017-3-9 14:45
标题: 五色LED音乐频谱源程序和原理图 原件清单
12单片机频闪，五色LED音乐频谱源程序代码

原理图如下：

元件清单：

五色LED音乐频谱源程序：
  1. #include <stc12c5a60s2.h>
  2. #include<intrins.h>
  3. #define BIN(n) LongToBin(0x##n##)
  4. #define uchar  unsigned char
  5. #define uint  unsigned int
  6. #define SAMPLE_NUM 64
  7. #define NUM_2_LOG 6
  8. #define FFT_OUT_MIN 3
  9. #define LongToBin(n) (((n>>21)&0x80)|((n>>18)&0x40)|((n>>15)&0x20)|((n>>12)&0x10)|((n>>9)&0x08)|((n>>6)&0x04)|((n>>3)&0x02)|((n)&0x01))

  11. uchar code BRTable[SAMPLE_NUM] ={ 0, 32, 16, 48, 8, 40, 24, 56,4, 36, 20, 52, 12, 44, 28, 60, 2, 34, 18, 50, 10, 42, 26, 58, 6, 38, 22, 54, 14, 46, 30, 62, 1, 33, 17, 49, 9, 41, 25, 57,5, 37, 21, 53, 13, 45, 29, 61,3, 35, 19, 51, 11, 43, 27, 59,7, 39, 23, 55, 15, 47, 31, 63};
  12. char code sin_tabb[SAMPLE_NUM] = { 0 ,12 ,25 ,37 ,49 ,60 ,71 ,81 ,90 ,98 ,106 ,112 ,117 ,122 ,125 ,126 ,127 ,126 ,125 ,122 ,117 ,112 ,106 ,98 ,90 ,81 ,71 ,60 ,49 ,37 ,25 ,12 ,0 ,-12 ,-25 ,-37 ,-49 ,-60 ,-71 ,-81 ,-90 ,-98 ,-106 ,-112 ,-117 ,-122 ,-125 ,-126 ,-127 ,-126 ,-125 ,-122 ,-117 ,-112 ,-106 ,-98 ,-90 ,-81 ,-71 ,-60 ,-49 ,-37 ,-25 ,-12  };
  13.                                  
  14. char code cos_tabb[SAMPLE_NUM] = {127 ,126 ,125 ,122 ,117 ,112 ,106 ,98 ,90 ,81 ,71 ,60 ,49 ,37 ,25 ,12 ,0 ,-12 ,-25 ,-37 ,-49 ,-60 ,-71 ,-81 ,-90 ,-98 ,-106 ,-112 ,-117 ,-122 ,-125 ,-126 ,-127 ,-126 ,-125 ,-122 ,-117 ,-112 ,-106 ,-98 ,-90 ,-81 ,-71 ,-60 ,-49 ,-37 ,-25 ,-12 ,0 ,12 ,25 ,37 ,49 ,60 ,71 ,81 ,90 ,98 ,106 ,112 ,117 ,122 ,125 ,126 };
  15. uchar a[21];
  16. uchar keep,keepnum,anum,timernum,timernum2,lednum3,Ltime;//用于分离

  18. /*加入数组用于显示相应led灯数目*/
  19. uchar lednum[]={0x00,0x01,0x03,0x07,0x0f,0x1f,0x3f,0x7f,0xff};//0-7的显示数组  P2组控制
  20. int xdata FftReal[SAMPLE_NUM];
  21. int xdata FftImage[SAMPLE_NUM];
  22. sbit p30=P3^0;
  23. sbit p31=P3^1;
  24. sbit p32=P3^2;
  25. sbit p33=P3^3;
  26. sbit p34=P3^4;
  27. sbit p35=P3^5;//9-11的led控制
  28. sbit p36=P3^6;
  29. sbit p37=P3^7;

  31. void timerinit()//定时器 初始化函数
  32. {
  33.          TMOD=0x01;
  34.          TH0=(65536-6000)/256;
  35.          TL0=(65536-6000)%256;
  36.          EA=1;
  37.          ET0=1;
  38.          TR0=1;
  39. }
  40. void disp()
  41. {
  42.         timernum++;
  43.         if(timernum==6) timernum=1;
  44.         P2=0;//显示前先关闭
  45.         P3=P3&0x1f;
  46.         switch(timernum)
  47.         {
  48.                 case 1:anum=a[0];p34=0;p33=1;p32=1;p31=1;p30=1;break;
  49.                 case 2:anum=a[1];p34=1;p33=0;p32=1;p31=1;p30=1;break;
  50.                 case 3:anum=a[2];p34=1;p33=1;p32=0;p31=1;p30=1;break;
  51.                 case 4:anum=a[3];p34=1;p33=1;p32=1;p31=0;p30=1;break;
  52.                 case 5:anum=a[4];p34=1;p33=1;p32=1;p31=1;p30=0;break;
  53.         }
  54.         //anum=a[10];/*修改可以改变光柱高度 （anum值分开几个部分用定时器区分显示）(a[]内逐加) */
  55.         if(anum<=8){P2=lednum[anum];P3=P3&0x1f;}//屏蔽高三位
  56.         //else {P2=0xff;P3=lednum2[anum-9];}
  57.         if(anum==9){P2=0xff;p35=1;p36=0;p37=0;}
  58.         if(anum==10){P2=0xff;p35=1;p36=1;p37=0;}
  59.         if(anum==11){P2=0xff;p35=1;p36=1;p37=1;}

  61. }

  63. uchar STC_ADC()                   //！！根据数据手册写一个ad读取函数
  64.   {
  65.      uchar i;
  66.         ADC_RES   = 0;
  67.         ADC_RESL  = 0;
  68.   ADC_CONTR = BIN(10001000);
  69.   i=3;
  70.         while(i--);
  71.         while (1)                     
  72.      {
  73.          if (ADC_CONTR & BIN(10000))      
  74.          {
  75.             break;
  76.            }
  77.      }
  78.      ADC_CONTR = BIN(10000000);
  79.   return( ADC_RESL<<2) ;
  80. }
  81. short sqrt_16( unsigned long M)   
  82. {
  83.     unsigned int N, i;
  84.     unsigned long tmp, ttp;
  85.     if( M == 0 )            
  86.         return 0;
  87.    
  88.     N = 0;
  89.    
  90.     tmp = ( M >> 30 );        
  91.     M <<= 2;
  92.     if( tmp > 1 )            
  93.     {
  94.         N ++;               
  95.         tmp -= N;
  96.     }
  97.    
  98.     for( i=15; i>0; i-- )   
  99.     {
  100.         N <<= 1;           
  101.         
  102.         tmp <<= 2;
  103.         tmp += (M >> 30);  
  104.         
  105.         ttp = N;
  106.         ttp = (ttp<<1)+1;
  107.         
  108.         M <<= 2;
  109.         if( tmp >= ttp )   
  110.         {
  111.             tmp -= ttp;
  112.             N ++;
  113.         }      
  114.     }
  115.    
  116.     return N;
  117. }
  118. void FFT()
  119. {
  120. register    uchar i,bb,j,k,p,max;
  121. register short TR,TI,temp;
  122.     unsigned long ulReal;                             
  123.     unsigned long ulImage;
  124.    
  125.    
  126.                                                                  
  127. for(i=0; i<SAMPLE_NUM;i++)          //此处可以加入自动增益
  128. {
  129.                 FftReal[BRTable[i]] = STC_ADC()<<keep;//使显示保持在一定范围内
  130.         FftImage[i] = 0;
  131. }

  133.          keepnum=FftReal[2]/32;//提取等级数
  134.                    if((7<keepnum)&&(keepnum<=8)) {keep=1;}
  135.          else if((4<keepnum)&&(keepnum<=6)) {keep=2;}
  136.          else if((2<keepnum)&&(keepnum<=4)) {keep=3;}
  137.          else {keep=5;}       
  138.    
  139.     for( i=1; i<=NUM_2_LOG; i++)                          
  140.     {
  141.         bb=1;
  142.         bb <<= (i-1);                                      
  143.         for( j=0; j<=bb-1; j++)                           
  144.         {
  145.             p=1;
  146.             p <<= (NUM_2_LOG-i);            
  147.             p = p*j;
  148.             for( k=j; k<SAMPLE_NUM; k=k+2*bb)               
  149.             {
  150.                 TR = FftReal[k]; TI = FftImage[k]; temp = FftReal[k+bb];
  151.                 FftReal[k] = FftReal[k] + ((FftReal[k+bb]*cos_tabb[p])>>7) + ((FftImage[k+bb]*sin_tabb[p])>>7);
  152.                 FftImage[k] = FftImage[k] - ((FftReal[k+bb]*sin_tabb[p])>>7) + ((FftImage[k+bb]*cos_tabb[p])>>7);
  153.                 FftReal[k+bb] = TR - ((FftReal[k+bb]*cos_tabb[p])>>7) - ((FftImage[k+bb]*sin_tabb[p])>>7);
  154.                 FftImage[k+bb] = TI + ((temp*sin_tabb[p])>>7) - ((FftImage[k+bb]*cos_tabb[p])>>7);
  155.                
  156.                 FftReal[k]  >>= 1;            
  157.                 FftImage[k]  >>= 1;
  158.                 FftReal[k+bb]  >>= 1;                 
  159.                 FftImage[k+bb]  >>= 1;
  160.                                                                               
  161.             }  
  162.         }
  163.     }
  164.     max=0;
  165.     for( i=0; i<5; i++)//5
  166.     {  
  167.         ulReal = FftReal[i+1];
  168.         ulReal *= ulReal;
  169.         ulImage = FftImage[i+1];
  170.         ulImage *= ulImage;
  171.         
  172.         a[i] = sqrt_16( ulReal + ulImage );   //修改
  173.                         
  174.         if( a[i] < FFT_OUT_MIN )     
  175.             a[i] = 0;//修改
  176.         else
  177.           a[i] = a[i]-FFT_OUT_MIN;
  178.         if( a[i] >max)
  179.              max =a[i];
  180.                          //disp();                     
  181.     }
  182.     if(max>11) //11
  183.     {
  184.        max/=11;       
  185.         for( i=0; i<5; i++) //输出a的5个分离数值
  186.         {      
  187.              a[i]/=max;
  188.             
  189.         }  
  190.     }
  191. }            


  194. void main()
  195. {
  196.         P2M0=0xff;//  (11111111);//P2组设置为推挽输出
  197.         P2M1=0;

  199.         …………余下代码请下载附件…………
复制代码



可看得出用的是stc单片机
下载：
频谱.rar (181.27 KB, 下载次数: 122)
作者: wyj841224    时间: 2017-12-19 18:17
谢谢，学习了。
作者: queen123444    时间: 2019-4-2 20:33
这个程序能用吗？
作者: 1580418512    时间: 2019-5-2 09:50
学习学习
作者: 1915794157    时间: 2020-5-16 17:16
这个原理图是用什么软件画的



欢迎光临 (http://www.51hei.com/bbs/) Powered by Discuz! X3.1