51音乐频谱

晓毕8

单独使用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;
    
85. 
    
86.         for(i=0;i<32;i+=4)
    
87.         {
    
88.                 float t0=0,t2=0;
    
89. 
    
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;
    
92. 
    
93.                 tmp1=(uint)(t0/3);
    
94.                 tmp2=(uint)(t2/3);
    
95. 
    
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格
     
114. 
     
115.                         /*UART_Put_Num(tmp1);
     
116.                         UART_Send_Str("  ");
     
117.                         UART_Put_Num(tmp2);
     
118.                         UART_Send_Enter(); */
     
119. 
     
120. 
     
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.         }
     
128. 
     
129. }
     
130. 
     
131. 
     
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. {
     
140. 
     
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);
     
150. 
     
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. }

复制代码

semonpic

顶楼主，谢谢。

成就与价值

求楼主分享！