51单片机

51音乐频谱

2020-01-27 11:22发布

站内问答 / 51单片机
6907 3 3
单独使用5110可以显示，单独使用FFT()和showbar（）函数用串口调试可以，为什么我把showbar（）里的数值显示在5110里就不行了，郁闷啊，请教大神
  1. #include"STC12C5A60S2.h"
  2. #include"5110.h"
  3. #include"ADDA.h"
  4. #include"math.h"
  5. #define  uchar unsigned char
  6. #define  uint  unsigned int
  7. struct compx                                                                           //定义复数结构体
  8. {
  9.         float real;//实部
  10.         float imag;//虚部
  11. };
  12. xdata struct compx s[ 64 ];  //数据缓冲区                         //FFT数据缓存放在XDATA空间
  13. struct compx EE(struct compx a1,struct compx b2)         //复数乘法
  14. {
  15.         struct compx b3;
  16.         b3.real=a1.real*b2.real-a1.imag*b2.imag; //实部1乘实部2-
  17.         b3.imag=a1.real*b2.imag+a1.imag*b2.real;
  18.         return(b3);
  19. }
  20. void FFT(struct compx xin[],int N)//64点                                  
  21. {
  22.         int f;
  23.         int m;
  24.         int nv2; //
  25.         int nm1;
  26.         int i;
  27.         int k;
  28.         int j=1;
  29.         int l;
  30.         struct compx v,w,t;
  31.         nv2=N/2;
  32.         f=N;
  33.         for(m=1;(f=f/2)!=1;m++){;}
  34.         nm1=N-1;
  35.         for(i=0;i<nm1;i++)                                           //倒序操作
  36.         {
  37.                 if(i<j)
  38.                 {
  39.                         t=xin[j];
  40.                         xin[j]=xin[i];
  41.                         xin[i]=t;
  42.                 }
  43.                 k=nv2;                                                       //k为倒序中相应位置的权值
  44.                 while(k<j)
  45.                 {
  46.                         j=j-k;
  47.                         k=k/2;
  48.                 }
  49.                 j=j+k;
  50.         }
  51.          
  52.         {
  53.                 int le,lei,ip;
  54.                   float pi;
  55.                   for(l=1;l<=m;l++)
  56.                    {
  57.                 le=pow(2,l);                                                 //乘方
  58.                     lei=le/2;
  59.                     pi=3.14159265;
  60.                     v.real=1.0;
  61.                 v.imag=0.0;  
  62.                    w.real=cos(pi/lei);                                           //旋转因子
  63.                     w.imag=-sin(pi/lei);
  64.                    
  65.                 for(j=1;j<=lei;j++)                                           //控制蝶形运算的级数
  66.                      {
  67.                         for(i=j-1;i<N;i=i+le)                                         //控制每级蝶形运算的次数
  68.                               {
  69.                                 ip=i+lei;
  70.                                        t=EE(xin[ ip ],v);
  71.                                        xin[ ip ].real=xin[ i ].real-t.real;   //蝶形计算
  72.                                        xin[ ip ].imag=xin[ i ].imag-t.imag;
  73.                                        xin[ i ].real=xin[ i ].real+t.real;
  74.                                        xin[ i ].imag=xin[ i ].imag+t.imag;
  75.                               }
  76.                               v=EE(v,w);   
  77.                       }     
  78.                    }
  79.           }
  80. }
  81. void showbar(void)
  82. {
  83.         uchar i,j,dat1,dat2;
  84.         uint  tmp1=33,tmp2=33;

  86.         for(i=0;i<32;i+=4)
  87.         {
  88.                 float t0=0,t2=0;

  90.                 t0=sqrt(pow((s[i  ].real+s[i+1].real),2)+pow((s[i  ].imag+s[i+1].imag),2))/2;
  91.                 t2=sqrt(pow((s[i+2].real+s[i+3].real),2)+pow((s[i+2].imag+s[i+3].imag),2))/2;

  93.                 tmp1=(uint)(t0/3);
  94.                 tmp2=(uint)(t2/3);

  96.                 for(j=2;j>0;j--)//行处理
  97.                 {                 
  98.                         if (tmp1>=8)
  99.                             {
  100.                                      dat1=7;
  101.                                     tmp1=tmp1-8;
  102.                                 }                         //满格
  103.                         else  if(tmp1==0)
  104.                             dat1=20;                                        //空格
  105.                         else  
  106.                             {
  107.                                 dat1=tmp1-1;
  108.                                 tmp1=0;
  109.                                 }                                        //1-7格
  110.                                                
  111.                         if (tmp2>=8){dat2=7;tmp2=tmp2-8;}                         //满格
  112.                         else  if(tmp2==0) dat2=20;                                        //空格
  113.                         else  {dat2=tmp2-1;tmp2=0;}                                        //1-7格

  115.                         /*UART_Put_Num(tmp1);
  116.                         UART_Send_Str("  ");
  117.                         UART_Put_Num(tmp2);
  118.                         UART_Send_Enter(); */


  121.                         LCD_write_pinpu(i/2,j-1,dat1); //dat1为0-8之间的数
  122.                         LCD_write_pinpu(i/2+1,j-1,dat2);
  123.                 }                        
  124.                                                                      //\\\\\  这里是显示的部分，重点是确定I(频率)，J，dat（幅度） 的值
  125.                                                                         // \\\\\\\\\\\\\\\\\\\\\         \\\\\\///////////////////
  126.                                                                          //\\\\\\\\\\\\         //////////////////////////////////////////////////
  127.         }

  129. }


  132. void delay_ms()
  133. {
  134.   int i,j;
  135.   for(i=0;i<6000;i++)
  136.     for(j=0;j<200;j++);
  137. }
  138. void main()
  139. {

  141.      int N=64,i;                                                             //变量初始化，64点FFT运算
  142.          float offset;
  143.            LCD_Init();
  144.   LCD_clear();
  145.          AD_init();    //A/D转换初始化
  146.      //UART_Init(); //串口初始化
  147.         // UART_Send_Str("串口调试AD值");
  148.         // UART_Send_Enter();
  149.          offset = AD_get(1);

  151.         while(1)
  152.         {         
  153.                      for(i=0;i<N;i++)                                      //采集音频信号
  154.                           {
  155.                                 ADC_CONTR=0xC8;                                         //40.96K采样率
  156.                                 while(!(ADC_CONTR&0x10));
  157.                                    s[i].real=(AD_get(1)-offset)/4;//((((int)ADC_DATA-128)/2))*4;
  158.                                    s[i].imag=0;
  159.                           }
  160.                     FFT(s,N);                                       //调用FFT函数进行变换
  161.                           showbar();                                                        //显示频谱                  
  162.         }
  163.   /*int i;
  164.   LCD_Init();
  165.   LCD_clear();
  166.   while(1)
  167.   {
  168.      for(i=0;i<8;i++)
  169.          {
  170.        LCD_write_pinpu(1,1,i);
  171.        delay_ms();
  172.          }
  173.   }        */
  174. }
复制代码
友情提示: 此问题已得到解决，问题已经关闭,关闭后问题禁止继续编辑，回答。
3条回答
成就与价值
2020-01-27 16:11
求楼主分享！

一周热门 更多>