fft-music-display.docx

#include#define uchar unsigned char#define uint unsigned int #define channel 0x01 /设置AD通道为 P1.1/-sbit SDA_R=P12;sbit SDA_R_TOP=P13;sbit SDA_G=P14;sbit SDA_G_TOP=P15; sbit STCP=P16;sbit SHCP=P17;/-/-/放大128倍后的sin整数表（128）code char SIN_TAB128 = 0, 6, 12, 18, 24, 30, 36, 42, 48, 54, 59, 65, 70, 75, 80, 85, 89, 94, 98, 102, 105, 108, 112, 114, 117, 119, 121, 123, 124, 125, 126, 126, 126, 126, 126, 125, 124, 123, 121, 119, 117, 114, 112, 108, 105, 102, 98, 94, 89, 85, 80, 75, 70, 65, 59, 54, 48, 42, 36, 30, 24, 18, 12, 6, 0, -6, -12, -18, -24, -30, -36, -42, -48, -54, -59, -65, -70, -75, -80, -85, -89, -94, -98, -102, -105, -108, -112, -114, -117, -119, -121, -123, -124, -125, -126, -126, -126, -126, -126, -125, -124, -123, -121, -119, -117, -114, -112, -108, -105, -102, -98, -94, -89, -85, -80, -75, -70, -65, -59, -54, -48, -42, -36, -30, -24, -18, -12, -6 ;/放大128倍后的cos整数表（128）code char COS_TAB128 = 127, 126, 126, 125, 124, 123, 121, 119, 117, 114, 112, 108, 105, 102, 98, 94, 89, 85, 80, 75, 70, 65, 59, 54, 48, 42, 36, 30, 24, 18, 12, 6, 0, -6, -12, -18, -24, -30, -36, -42, -48, -54, -59, -65, -70, -75, -80, -85, -89, -94, -98, -102, -105, -108, -112, -114, -117, -119, -121, -123, -124, -125, -126, -126, -126, -126, -126, -125, -124, -123, -121, -119, -117, -114, -112, -108, -105, -102, -98, -94, -89, -85, -80, -75, -70, -65, -59, -54, -48, -42, -36, -30, -24, -18, -12, -6, 0, 6, 12, 18, 24, 30, 36, 42, 48, 54, 59, 65, 70, 75, 80, 85, 89, 94, 98, 102, 105, 108, 112, 114, 117, 119, 121, 123, 124, 125, 126, 126 ;/采样存储序列表code char LIST_TAB128 = 0, 64, 32, 96, 16, 80, 48, 112,8, 72, 40, 104, 24, 88, 56, 120,4, 68, 36, 100, 20, 84, 52, 116,12, 76, 44, 108, 28, 92, 60, 124,2, 66, 34, 98, 18, 82, 50, 114,10, 74, 42, 106, 26, 90, 58, 122,6, 70, 38, 102, 22, 86, 54, 118,14, 78, 46, 110, 30, 94, 62, 126,1, 65, 33, 97, 17, 81, 49, 113,9, 73, 41, 105, 25, 89, 57, 121,5, 69, 37, 101, 21, 85, 53, 117,13, 77, 45, 109, 29, 93, 61, 125,3, 67, 35, 99, 19, 83, 51, 115,11, 75, 43, 107, 27, 91, 59, 123,7, 71, 39, 103, 23, 87, 55, 119,15, 79, 47, 111, 31, 95, 63, 127;uchar COUNT=0,COUNT1=0,ADC_Count=0,LINE=15,G,T;uchar i,j,k,b,p; int Temp_Real,Temp_Imag,temp; / 中间临时变量 uint TEMP1;int xdata Fft_Real128; int xdata Fft_Image128; / fft的虚部 uchar xdata LED_TAB264;/记录 漂浮物 是否需要 停顿一下uchar xdata LED_TAB64;/记录红色柱状 uchar xdata LED_TAB164;/记录 漂浮点void Delay(uint a)while(a-);void FFT() /uchar X; 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; k7) + (Fft_Imagek+b*SIN_TABp)7); Fft_Imagek = Fft_Imagek - (Fft_Realk+b*SIN_TABp)7) + (Fft_Imagek+b*COS_TABp)7); Fft_Realk+b = Temp_Real - (Fft_Realk+b*COS_TABp)7) - (Fft_Imagek+b*SIN_TABp)7); Fft_Imagek+b = Temp_Imag + (temp*SIN_TABp)7) - (Fft_Imagek+b*COS_TABp)7); / 移位.防止溢出. 结果已经是本值的 1/64 Fft_Realk = 1; Fft_Imagek = 1; Fft_Realk+b = 1; Fft_Imagek+b = 1; /X=(Fft_Real1* Fft_Real1)+(Fft_Image1*Fft_Image1)7); Fft_Real0=Fft_Image0=0; /去掉直流分量/ Fft_Real63=Fft_Image63=0;for(j=0;j1);/求功率if(TEMP11)TEMP1-;else TEMP1=0;if(TEMP131)TEMP1=31; if(TEMP1(LED_TABj)LED_TABj=TEMP1;if(TEMP1(LED_TAB1j) LED_TAB1j=TEMP1;LED_TAB2j=18;/提顿速度=12 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速度为70周期一次； /0中断标志清零；1启动adc(ADC_START);001AD通道打开（这里为P1.1）;/- P2M0=1; P0M0=1; TMOD=0X12; TH0=0x30; /大约20K的采样率（要完整频段需40K以上。但音频中10k以下居多，故本人选择20K采样，美观些） TL0=0x30; TH1=0xEE; TL1=0XC0; ET0=1; /定时器0 打开 TR0=0; /关闭定时器 ET1=1; TR1=1; PT1=0; PT0=1; IPH=PADCH; IP=PADC; /中断优先级 EA=1;/总中断打开 void ADC_Finish() interrupt 5 ADC_CONTR &= !ADC_FLAG; Fft_RealLIST_TABADC_Count=(int)(ADC_RES)1)-256;/-512; /按LIST_TAB表里的顺序，进行存储 采样值,/ ADC_CONTR = ADC_POWER | ADC_SPEEDHH| ADC_START | channel;/ 为了采集负电压，采用 偏置采集。电压提高到1/2 vcc，所以要减去256 if(ADC_Count=127)ADC_Count+; else EADC=0;TR0=0; void LED_Display() interrupt 3 /中断一次 显示一行。 TH1=0xF3; TL1=0X00;for(G=0;G64;G+) /往点阵屏填充 一行的 数据if(LED_TABG=LINE+16)SDA_R_TOP=1;else SDA_R_TOP=0; if(LED_TABG0)LINE-;else LINE=15; / if(LED_TABCOUNT0)LED_TABCOUNT-; /柱状递减， COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; /柱状递减， COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; COUNT+; if(LED_TABCOUNT0)LED_TABCOUNT-; COUNT+; if(COUNT=64)COUNT=0; /漂浮物递减 if(LED_TAB2COUNT1=0)/判断是否需要停顿 if(LED_TAB1COUNT1LED_TABCOUNT1)LED_TAB1COUNT1-;/大于柱状则递减（保持漂浮物在柱状之上） else LED_TAB2COUNT1-; COUNT1+; if(LED_TAB2COUNT1=0) if(LED_TAB1COUNT1LED_TABCOUNT1)LED_TAB1COUNT1-; else LED_TAB2COUNT1-; COUNT1+; if(LED_TAB2COUNT1=0)/判断是否需要停顿 if(LED_TAB1COUNT1LED_TABCOUNT1)LED_TAB1COUNT1-;/大于柱状则递减（保持漂浮物在柱状之上） else LED_TAB2COUNT1-; COUNT1+; if(LED_TAB2COUNT1=0) if(LED_TAB1COUNT1LED_TABCOUNT1)LED_TAB1COUNT1-; else LED_TAB2COUNT1-; COUNT1+; if(C