【正文】 產(chǎn)生波形 for ( i=0。iSAMPLENUMBER。i++ ) { fWaveR[i]=INPUT[i]。//實部 fWaveI[i]=。 //虛部 w[i]=。 } FFT(fWaveR，fWaveI)。 //調(diào)用FFT子程序計算功率譜 for ( i=0。iSAMPLENUMBER。i++ ) { DATA[i]=w[i]。 //傳送數(shù)據(jù) } while ( 1 )。 //斷點}void FFT(float dataR[SAMPLENUMBER]，float dataI[SAMPLENUMBER]){ int x0，x1，x2，x3，x4，x5，x6，xx。 //定義xn int i，j，k，b，p，L。 //中間變量 float TR，TI，temp。 //中間變量 /********** 反演序列代碼 ************/為蝶形算法做準備 for ( i=0。iSAMPLENUMBER。i++ ) { x0=x1=x2=x3=x4=x5=x6=0。 x0=iamp。0x01。 x1=(i/2)amp。0x01。 x2=(i/4)amp。0x01。 x3=(i/8)amp。0x01。x4=(i/16)amp。0x01。 x5=(i/32)amp。0x01。 x6=(i/64)amp。0x01。 xx=x0*64+x1*32+x2*16+x3*8+x4*4+x5*2+x6。 dataI[xx]=dataR[i]。 } for ( i=0。iSAMPLENUMBER。i++ ) { dataR[i]=dataI[i]。 dataI[i]=0。 } /************** FFT算法代碼 *******************/ for ( L=1。L=7。L++ ) { /*第一層 */運算出第二層的運算次數(shù) b=1。 i=L1。 while ( i0 ) { b=b*2。 i。 } /* b= 2^(L1) */ for ( j=0。j=b1。j++ ) /* 第二層 */ { p=1。 i=7L。 while ( i0 ) /* p=pow(2，7L)*j。 */pow(x，y)計算x的y次冪 { p=p*2。 i。 } p=p*j。 for ( k=j。k128。k=k+2*b ) /* 第三層 */ { TR=dataR[k]。 TI=dataI[k]。 temp=dataR[k+b]。 dataR[k]=dataR[k]+dataR[k+b]*cos_tab[p]+dataI[k+b]*sin_tab[p]。 dataI[k]=dataI[k]dataR[k+b]*sin_tab[p]+dataI[k+b]*cos_tab[p]。dataR[k+b]=TRdataR[k+b]*cos_tab[p]dataI[k+b]*sin_tab[p]。 dataI[k+b]=TI+temp*sin_tab[p]dataI[k+b]*cos_tab[p]。 } /* 第三層結(jié)束 */ } /* 第二層結(jié)束 */ } /* 第一層結(jié)束*/ for ( i=0。iSAMPLENUMBER/2。i++ ) { w[i]=sqrt(dataR[i]*dataR[i]+dataI[i]*dataI[i])。 //開方求出w[i]} /* END FFT */void InitForFFT() //初始化工作變量{ int i。 for ( i=0。iSAMPLENUMBER。i++ ) { sin_tab[i]=sin(PI*2*i/SAMPLENUMBER)。 //公式 cos_tab[i]=cos(PI*2*i/SAMPLENUMBER)。 }}void MakeWave() //波形發(fā)生子程序{ int i。 for ( i=0。iSAMPLENUMBER。i++ ) { INPUT[i]=sin(PI*2*i/SAMPLENUMBER*3)*1024。 //采樣 }}2. 顯示程序include include define uint unsigned int define uchar unsigned char sbit RS=P3^0。sbit RW=P3^1。sbit E=P3^2。sbit ADDA=P1^0。sbit ADDB=P1^1。sbit ADDC=P1^2。sbit START=P1^4。sbit OE=P1^5。sbit EOC=P3^3。sbit CLK=P3^5。unsigned char addata。unsigned int volt100。unsigned char volt[3]。 //***************************延時子程序***************************void delay(uint i){ uint j。 for(。i0。i) for(j=0。j125。j++)。}//***************************寫數(shù)據(jù)子程序*************************void write_(uchar ){ RS=0。 RW=0。 P2=。 E=1。 delay(1)。 E=0。}//***************************寫數(shù)據(jù)子程序*************************void write_dat(uchar dat){ RS=1。 RW=0。 P2=dat。 E=1。 delay(1)。 E=0。}//**********************光標定位函數(shù)******************************void LcdPos (uchar Xpos，uchar Ypos){ uchar tmp。 Xposamp。=0x0f。 Yposamp。=0x01。 tmp=Xpos。 if(Ypos==1) tmp|=0xc0。 tmp|=0x80。 write_(tmp)。}//***********************LCD初始化子程序**************************LCDRset(void){ E=0。//**************使能端禁止 write_(0X38)。 delay(1)。 write_(0X01)。//**清屏 delay(1)。 write_(0X0C)。//**開顯示光標不顯示 delay(1)。 write_(0x06)。//**光標加一 delay(1)。}//*******************在指定行列顯示字符子程序**********************void writeChar(uchar Xpos，uchar Ypos，char c){ LcdPos(Xpos，Ypos)。 write_dat(c)。}//********************在指定行列寫入指定字符串*********************void writeString(uchar Xpos，uchar Ypos，char s[]){ uchar p=0。 while(1) { writeChar(Xpos，Ypos，s[p])。 p++。 if(s[p]==0) break。 if(++Xpos=15) { Xpos=0。 Ypos^=1。//異或運算 } }}void Delay1ms(unsigned int count){ unsigned int i，j。 for(i=0。icount。i++) for(j=0。j120。j++)。} void main(void){ P1=0。 P2=0。 LCDRset()。while(1) { START=1。 START=0。 while(!EOC)。 OE=1。 addata=P0。int a[100],i,k,n,max。 //計算電壓最大值 scanf(%d,amp。n)。 for(i=0。in。i++) { scanf(%d,amp。a[i])。 } max=a[0]。 for(i=1。i=n1。i++) { if(a[i]max) { max=a[i]。 k=i+1。 } if(a[0]=a[i]) { max=a[0]。 k=1。 } }addata=max。 OE=0。 volt100=(unsigned int)1(addata*100)/51。 volt[0]=(volt100/100)%10+39。039。 volt[1]=(volt100/10)%10+39。039。 volt[2]=(volt100)%10+39。039。 // delay(1)。 writeString(0，0，OUT=)。 writeChar(4，0，volt[0])。 writeChar(5，0，39。.39。)。 writeChar(6，0，volt[1])。 writeChar(7，0，volt[2])。 writeChar(8，0，39。V39。)。 // delay(1)。 volt100=(unsigned int)(addata*100)/51/5。 volt[0]=(volt100/100)%10+39。039。 volt[1]=(volt100/10)%10+39。039。 volt[2]=(volt100)%10+39。039。 writeString(0，1，IN= a=)。 writeChar(11，1，volt[0])。 writeChar(12，1，39。.39。)。 writeChar(13，1，volt[1])。 writeChar(14，1，volt[2])。 }