微波EDA网,见证研发工程师的成长!
首页 > 研发问答 > 嵌入式设计讨论 > MCU和单片机设计讨论 > 关于音乐频谱FFT程序的求助

关于音乐频谱FFT程序的求助

时间:10-02 整理:3721RD 点击:

sbit led0 = P2^0;
sbit led1 = P2^1;
sbit led2 = P2^2;
sbit led3 = P2^3;
sbit led4 = P2^4;
sbit led5 = P2^5;
sbit led6 = P2^6;
sbit led7 = P2^7;
uint8 COUNT=15,COUNT1=0,ADC_Count=0,LINE=15,G=0,T;
uint8 i,j,k,b,p,anum,led=0;                 
int Temp_Real,Temp_Imag,temp;                // 中间临时变量  
uint16 TEMP1,max=0;
int xdata Fft_Real[128];
int xdata Fft_Image[128];               // fft的虚部
uint8 xdata LED_TAB[64];    //记录FFT计算出的各点幅值大小
void FFT()//基2fft
{                 
    for( i=1; i<=7; i++)                            /* for(1) */
    {
        b=1;
        b <<=(i-1);                                 //碟式运算,用于计算隔多少行计算例如 第一极 1和2行计算
        for( j=0; j<=b-1; j++)                      /* for (2) */
        {
            p=1;
            p <<= (7-i);            
            p = p*j;
            for( k=j; k<128; k=k+2*b)               /* for (3)  */
            {
                Temp_Real=Fft_Real[k];
    Temp_Imag=Fft_Image[k];
    temp=Fft_Real[k+b];
                Fft_Real[k]=Fft_Real[k]+((Fft_Real[k+b]*COS_TAB[p])>>7)+((Fft_Image[k+b]*SIN_TAB[p])>>7);
                Fft_Image[k]=Fft_Image[k]-((Fft_Real[k+b]*SIN_TAB[p])>>7)+((Fft_Image[k+b]*COS_TAB[p])>>7);
                Fft_Real[k+b]=Temp_Real-((Fft_Real[k+b]*COS_TAB[p])>>7)-((Fft_Image[k+b]*SIN_TAB[p])>>7);
                Fft_Image[k+b]=Temp_Imag+((temp*SIN_TAB[p])>>7)-((Fft_Image[k+b]*COS_TAB[p])>>7);          // 移位.防止溢出. 结果已经是本值的 1/64               
                Fft_Real[k] >>= 1;            
                Fft_Image[k] >>= 1;
                Fft_Real[k+b] >>= 1;                 
                Fft_Image[k+b] >>= 1;                                                               
            }     
        }
    }

for(j=0;j<64;j++)//64分频            
{                     
  TEMP1=((((Fft_Real[j+1]* Fft_Real[j+1]))+((Fft_Image[j+1]*Fft_Image[j+1])))>>1);//求各频段幅值
  if(TEMP1<1)TEMP1=0;
  LED_TAB[j]=TEMP1;
  if(LED_TAB[j]>max)max=LED_TAB[j];
}
if(max>64)//分级量化
{
   max/=64;
   for(j=0;j<64;j++)LED_TAB[j]/=max;
}      
}
void delay2us(void)   //2us延时
{
    unsigned char a,b;
    for(b=1;b>0;b--)
        for(a=3;a>0;a--);
}
void Init()
{  
     P1ASF = 0x02;                 //0000,0010, 将 P1.1 置成模拟口
     AUXR1 &=0xFB;                 //1111,1011, 令 ADRJ=0
  EADC=1;        //AD中断打开
  ADC_CONTR = ADC_POWER | ADC_SPEEDHH | ADC_START | channel;//1110 1001   1打开 A/D (ADC_POWER)转换电源;11速度为90周期一次;
                                                      //0中断标志清零;1启动adc(ADC_START);001AD通道打开(这里为P1.1);
  
  P2M0=1;
    TMOD=0X12;        
  TH0=0xb0;                                       
  TL0=0xb0;
  TH1=0xfe;               
  TL1=0Xd4;
  ET0=1;            //定时器0 打开
  TR0=0;            //关闭定时器
  ET1=1;
  TR1=1;
  PT1=0;
  PT0=1;
    IPH=PADCH;
  IP=PADC;        //中断优先级
  EA=1;         //总中断打开
}

没办法,还是FFT没学好,对于修改看不懂。

难于上青天啊啊啊啊我也卡了

Copyright © 2017-2020 微波EDA网 版权所有

网站地图

Top