【正文】 *********/include include unsigned ioport port8003。 double p,q,s,vr,vi,poddr,poddi。 is=0。 is=2*is+(m2*j)。 } pr[0]=。 pi[1]=sin(p)。 i++) { p=pr[i1]*pr[1]。 pi[i]=spq。 vi=fi[it]。 fi[it+1]=vifi[it+1]。 l0=0。 it=(m1)*nv。 q=pi[m*j]*fi[it+j+nv/2]。 poddi=spq。 fi[it+j]=fi[it+j]+poddi。 fi[i]=fi[i]/(*n)。 if (fabs(fr[i])*fabs(fi[i])) { if ((fi[i]*fr[i])0) pi[i]=。 double x[128],pr[128],pi[128],fr[128],fi[128],mo[128]。 n=128 for (。 i++) { *px++ = port8003。 i++) { xm=*px。 px++。i=n1。 pz++。 int ep,arg,mt,s0,s1。 if(c==NULL) exit(1)。 if(inv==0) { sign=。xi[i]=xi[i]/n。i++) { if(ij) { tra(amp。xi[i],amp。k=n/2。 while(i!=1) { ep=ep+1。in。 b=1。 s1=0。while(in){arg=i+b。pi=xr[arg]*s[s1]+xi[arg]*c[s1]。xi[i]=xi[i]+pi。 b=b*2。詳細(xì)探索了超聲波流量計的方案設(shè)計，知道了檢測方案大致可分為傳播速度差法（包括：時差法、相位差法、頻差法），多普勒法，相關(guān)法等。文章中作者采用DSP進(jìn)行相關(guān)算法的實現(xiàn)，畫出了流程圖，給出了系統(tǒng)整體電路硬件設(shè)計方案，基于此，并進(jìn)行了DSP的編程設(shè)計。在論文撰寫過程中，他提供了許多寶貴的思路和建議，結(jié)合工作體會和經(jīng)歷，提出了很多有價值的觀點，為完成本次設(shè)計和論文給予了極大的幫助。該程序基本上滿足了設(shè)計之初的設(shè)想要求。設(shè)計所涉及到的東西，是前所未有的。再次感謝所有支持和幫助過我的領(lǐng)導(dǎo)、老師和同學(xué)們！參 考 文 獻(xiàn)，:中國計量出版社，1991，流量測量技術(shù)及儀表，機(jī)械工業(yè)出版社，(第二版).水利電力出版社，1990，信號檢測理論(第二版)，哈爾濱工業(yè)大學(xué)出版社，軟件無線電原理與應(yīng)用，電子工業(yè)出版社，:科學(xué)技術(shù)出版社，1999，1993，:西安電子科技大學(xué)出版社，200010．，北京:北京化工大學(xué)碩士學(xué)位論，南京::天津大學(xué)碩士學(xué)位論13. Yao Zhenjing, Meng Qinghao, Li Genwang, et al.. Noncrosstalk realtime ultrasonic range system with optimized chaotic pulse positionwidth modulation excitation. in Ultrasonics Symposium, 2008. IUS 2008. IEEE. 2008.14. Fortuna, L., M. Frasca, and A. Rizzo. Chaotic Pulse Position Modulation to Improve the Efficiency of Sonar Sensors. IEEE Transactions on Instrumentation and Measurement., 2003. 52(6): p. 18091815.15. 程曉暢, 蘇紹景,王躍科等. 偽隨機(jī)碼超聲擴(kuò)頻測距系統(tǒng)設(shè)計與算法[J]. 測試技術(shù)學(xué)報, 2007. 21(1): p. 7984.16. 童峰, 許水源, 許天增. 一種高精度超聲波測距處理方法[J]. 廈門大學(xué)學(xué)報, 1998. 37(4): p. 507513.17. 梅璐璐, 林京. 信號群延遲估計方法研究. 儀器儀表學(xué)報, 2009(001): p. 5054. 18. ?？? 陳兵, 唐斌. 快速多分量 LFM 信號的檢測與參數(shù)估計方法. 電子測量與儀器學(xué)報, 2008. 22(001): p. 2529. 19. Zhao Xiaochuan, Luo Qingsheng, Han Baoling, et al.. A novel ultrasonic ranging system based on the selfcorrelation of pseudorandom sequence. in Information and Automation, 2009. ICIA 39。通過這次設(shè)計，培養(yǎng)了自學(xué)能力，為以后的繼續(xù)學(xué)習(xí)打下基礎(chǔ)。通過這次設(shè)計，使我受益非淺。和我同組的同學(xué)負(fù)責(zé)軟件的編寫，他花了大量的時間查閱資料，請教老師和同學(xué)。但是由于實驗條件的限制，對DSP硬件電路的聯(lián)調(diào)沒有實現(xiàn)。然后，詳細(xì)時差法測量原理及主要特點，分析了時差法計算流量的公式表達(dá)。 free(s)。}s1=s1+s0。xi[arg]=xi[i]pi。pi=xi[arg]。kb。mt=ep。 c[i]=cos(ph*i)。 } ph=2*M_PI/n。 } ep=0。 } k=n/2。xr[j])。 for(i=0。in。 if(s==NULL) exit(1)。 float *s,*c。 }IFFT只要先將FFT變換得到的X(k)取共軛，就可以直接利用FFT子程序，F(xiàn)FT運(yùn)算和IFFT運(yùn)算可以共用一個子程序。 zm = (int)(mo[i]*)。 pz = (int*)0x4080。 pr[i]=x[i]。 for (i=0。 /* 在這里采集數(shù)據(jù)到0x4000 */ if 0 for( i = 0。 int *px = (int*)0x4000。 } else pi[i]=atan(fi[i]/fr[i])*。 i=n1。 i=n1。 fi[it+j+nv/2]=fi[it+j]poddi。 s=s*(fr[it+j+nv/2]+fi[it+j+nv/2])。 j=(nv/2)1。 nv=2*nv。 nv=2。 fi[it]=vi+fi[it+1]。 it=n2。 s=(pr[i1]+pi[i1])*(pr[1]+pi[1])。 for (i=2。 p=(*n)。 } fr[it]=pr[is]。 i=k1。 it=n1。double pr[],pi[],fr[],fi[]。由于,是實序列，所以求得為，若直接幅值相關(guān)法的計算復(fù)雜度為O (N2)，而改進(jìn)后的FFT算法程序的計算復(fù)雜度為，將大大地提高運(yùn)算的速度。雖然傅立葉變換比較復(fù)雜，但現(xiàn)在已經(jīng)有高速DSP和成熟的FFT算法，而兩個頻譜的相乘只不過是一次多相式乘法。DFT離散傅立葉變換是連續(xù)傅里葉變換在離散系統(tǒng)中的表現(xiàn)，完成離散時域信號到