【正文】 表424頻率檢測(cè)誤差表（）輸入信號(hào)（KHz）頻率檢測(cè)結(jié)果（KHz）檢測(cè)誤差（KHz）141140112601260026232625+233943395+147574755251315130174397440+1103661036511276412765+11406414065+11629916230+118396183951 比較表423和表424可知，在輸入信號(hào)幅值改變的情況下，頻率檢測(cè)的誤差都滿足預(yù)先的設(shè)計(jì)指標(biāo)，即誤差值在177。第五章 總 結(jié)非均勻采樣理論是數(shù)字信號(hào)處理領(lǐng)域中一個(gè)全新的發(fā)展方向，是世界各國(guó)科研人員都非常關(guān)注的一項(xiàng)課題。同時(shí)，非均勻采樣理論和技術(shù)是一個(gè)非常龐大而復(fù)雜的研究領(lǐng)域，涉及到的基礎(chǔ)學(xué)科有測(cè)試策略、非線性理論、概率論與隨機(jī)過(guò)程、空間變換理論、信號(hào)分析和處理以及多種工程技術(shù)學(xué)科。由于計(jì)算機(jī)技術(shù)的迅速發(fā)展以及人們對(duì)信息處理技術(shù)愈來(lái)愈高的要求，使得人們對(duì)非均勻采樣理論和技術(shù)的研究往縱深方向發(fā)展。本文中我們主要探討了非均勻采樣的基本理論及應(yīng)用，結(jié)合實(shí)際工程實(shí)踐，給出了實(shí)現(xiàn)非均勻采樣的近似方法——非均勻周期采樣，并研制出一套基于DSP的非均勻采樣系統(tǒng)。目前，我們的主要工作是利用非均勻采樣的抗頻率混疊性能，進(jìn)行正弦信號(hào)的頻譜檢測(cè)，并對(duì)一些傳統(tǒng)信號(hào)處理方法應(yīng)用于非均勻采樣進(jìn)行了初步的嘗試，取得了一定的成果。在利用非均勻采樣優(yōu)點(diǎn)的同時(shí)，我們也發(fā)現(xiàn)了其具有的一些缺陷，如頻譜檢測(cè)中的頻譜噪聲和信號(hào)恢復(fù)很困難等。因此，下一步的研究工作應(yīng)主要集中于以下幾個(gè)方面：1）如何消除非均勻采樣頻譜檢測(cè)中出現(xiàn)的頻譜噪聲。非均勻采樣頻譜檢測(cè)中含有一定程度的頻譜噪聲，這種現(xiàn)象極大地限制了非均勻采樣理論的發(fā)展和應(yīng)用范圍。鑒于非均勻采樣理論的巨大應(yīng)用前景，必須找到處理頻譜噪聲的方法。2）找到非均勻采樣頻譜計(jì)算的快速算法。正如FFT極大地推動(dòng)了傳統(tǒng)數(shù)字信號(hào)處理理論的應(yīng)用和發(fā)展一樣，非均勻采樣理論能夠得到廣泛使用和飛速發(fā)展的前提是其必須存在頻譜計(jì)算的快速算法。但是，非均勻采樣與均勻采樣相比，有著很大的不同，如不滿足均勻采樣中的時(shí)移關(guān)系等等，因此，均勻采樣理論中的頻譜計(jì)算快速算法不能應(yīng)用于非均勻采樣，必須重新尋找非均勻采樣頻譜計(jì)算的快速算法。3）如何由非均勻采樣的離散信號(hào)重構(gòu)模擬信號(hào)。信號(hào)重建是信號(hào)處理領(lǐng)域中一個(gè)重要方面，無(wú)論采樣什么處理方法，都需要在處理后進(jìn)行信號(hào)的恢復(fù)。目前，對(duì)于非均勻采樣得到的離散信號(hào)，還沒(méi)有一個(gè)可以普遍使用的信號(hào)恢復(fù)方法。4）非均勻采樣信號(hào)在時(shí)域及頻域下的能量對(duì)應(yīng)關(guān)系。5）對(duì)基于DSP的非均勻采樣系統(tǒng)作進(jìn)一步的改進(jìn)，使之能夠檢測(cè)更高頻率的信號(hào)，同時(shí)，改進(jìn)PC機(jī)應(yīng)用程序，使之具有更多的功能，如波形顯示、頻譜處理等功能。上述五個(gè)方面中，前四個(gè)方面是當(dāng)前人們都很感興趣的理論研究方向，它們之間的關(guān)系不是獨(dú)立的，而是相互聯(lián)系、相互促進(jìn)的，對(duì)這些問(wèn)題進(jìn)行廣泛、深入和細(xì)致的研究，無(wú)論在理論上還是在實(shí)際應(yīng)用上都具有重大而深遠(yuǎn)的意義。致 謝 在畢業(yè)論文即將完成之際，借此機(jī)會(huì)我想感謝一直以來(lái)關(guān)心，幫助我的老師，同學(xué)，朋友當(dāng)然還有我的親人，是你們的幫助和關(guān)系讓我順利完成了這項(xiàng)任務(wù)。 感謝我的指導(dǎo)導(dǎo)師***老師，本文從選題、論證、研究到最后完成，無(wú)不凝聚著指導(dǎo)導(dǎo)師的心血。*老師在理論上學(xué)識(shí)淵博，獨(dú)具匠心，他誨人不倦的精神，嚴(yán)謹(jǐn)治學(xué)的態(tài)度和嚴(yán)格的要求深深地影響著我，使我終生難忘。從他的身上，我學(xué)到了許多書(shū)本上難以學(xué)到的東西，對(duì)我今后的人生之路產(chǎn)生深遠(yuǎn)的影響。參考文獻(xiàn)[1] 趙洪亮等. TMS320C55xDSP應(yīng)用系統(tǒng)設(shè)計(jì). 北京：北京航空航天大學(xué)出版社,2008[2] 薛定宇. 控制系統(tǒng)仿真與計(jì)算機(jī)輔助設(shè)計(jì). 北京：機(jī)械工業(yè)出版社,2005[3] 王正林，劉明. 精通Matlab7. 北京：電子工業(yè)出版社，2006[4] 何薇. 計(jì)算機(jī)圖形圖象處理技術(shù)與應(yīng)用. 北京：清華大學(xué)出版社, 2007[5] 汪安民. 非均勻采樣信號(hào)的頻譜分析及信號(hào)處理系統(tǒng). 電子技術(shù)應(yīng)用,2006[6] 張勇. TMS320C5000系列DSP匯編語(yǔ)言程序設(shè)計(jì). 陜西西安：西安電子科技大學(xué)出版社,2004[7] 汪安民. 基于非均勻采樣的信號(hào)頻率檢測(cè)方法及其實(shí)現(xiàn). 湖北武漢：華中科技大學(xué)博士學(xué)位論文, 2004[8] 李玉生. 超高速并行采樣模擬/數(shù)字轉(zhuǎn)換的研究. 安徽合肥：中國(guó)科學(xué)技術(shù)大學(xué)博士學(xué)位論文, 2007[9] 侯利軍. 基于DSP的非均勻采樣系統(tǒng)的設(shè)計(jì)與實(shí)現(xiàn). 湖北武漢：華中科技大學(xué)碩士學(xué)位論文, 2005 [10] 石強(qiáng). 基于非均勻采樣的高速數(shù)據(jù)采集系統(tǒng)設(shè)計(jì). 四川成都：電子科技大學(xué)碩士學(xué)位論文，2007[11] 初仁辛，孫圣和. 非均勻采樣信號(hào)理論及其發(fā)展 黑龍江哈爾濱：電子科技導(dǎo)報(bào)，2000附錄基于MATLAB語(yǔ)言仿真程序 程序開(kāi)始：clear all。 %清所有的變量clf。 %清屏n=1224。 %設(shè)置數(shù)據(jù)點(diǎn)數(shù)f1=200。 %設(shè)置信號(hào)1頻率f2=700。 %設(shè)置信號(hào)2頻率f3=1100。 %設(shè)置信號(hào)3頻率dt=1/1024。 %設(shè)置采樣時(shí)間sample=1/dt。 %設(shè)置采樣頻率randn(‘state’,sum(100*clock))。 %產(chǎn)生隨機(jī)噪聲dis=2*rand(1,N)。dis=dismean(dis)。 %將隨機(jī)噪聲的均值設(shè)置成0dis=2*dis。 %擴(kuò)大隨機(jī)噪聲的幅度f(wàn)or k=1:N。 %產(chǎn)生待檢測(cè)的信號(hào)w(k)=1*sin(2*pi*f1*k*dt)+1*sin(2*pi*f2*k*dt)+1*sin(2* pi*f3*k*dt)+dis(k)。endfigure(1)。 %計(jì)算待檢測(cè)信號(hào)的均勻采樣頻譜fw=fft(w)。s=abs(fw).^2。f=(1:N1)*sample/N。plot(f,s(2:N))。xlabel(‘Frequency/Hz’)。 %橫坐標(biāo)說(shuō)明ylabel(‘Amplitude’)。 %縱坐標(biāo)說(shuō)明title(‘DFT of uniformly sampling signal’)。t=tnonunif(N,1,1)。 %調(diào)用非均勻采樣時(shí)間函數(shù)for k=1:N。 %產(chǎn)生非均勻采樣數(shù)據(jù) w2(k)=1*sin(2*pi*f1*t(k)*dt)+ 1*sin(2*pi*f2*t(k)*dt)+ 1*sin(2*pi*f3*t(k)*dt)+dis(k)。endfor k=1:N。 %計(jì)算非均勻采樣數(shù)據(jù)的頻譜 s2(k)=0。 for n=1:N。 %計(jì)算非均勻采樣數(shù)據(jù)的傅里葉變換 s2(k)=s2(k)+w2(n)*(cos(2*pi*t(n)*(k1)/N)j*sin(2*pi*t(n)*(k1)/N))。 end s2(k)=abs(s2(k)).^2。 end figure(2)。 %畫(huà)出非均勻采樣數(shù)據(jù)的頻譜f=(1:N1)*sample/N。plot(f,s2(2:N))。xlabel(‘Frequency/Hz’)。ylabel(‘Amplitude’)。title(‘DFT of additive random sampling signal’)。f=(1:N1)*sample/N。plot(f,s2(2:N))。xlabel(‘Frequency/Hz’)。ylabel(‘Amplitude’)。title(‘DFT of additive random sampling signal’)。.error(nargchk(3,7,nargin))。 %該函數(shù)定義一段隨機(jī)時(shí)間信號(hào) if nargin4, mode=‘a(chǎn)dd’。 elseif~strcmp(mode, ‘ per’)amp。~strcmp(mode, ‘a(chǎn)dd’)。 error(‘invalid mode’)。 end if nargin5 if ~strcmp (option, ‘unif’)amp。~strcmp(option, ‘norm’)。 error(‘invalid distribution option’)。 end else option=‘unif’。 endif nargin6v=version。if v(1)4 rand(‘state’,J)。 randn(‘state’,J)。 else rand(‘seed’,J)。 randn(‘seed’,J)。 end endt(1)=0。if strcmp(option, ‘unif’)。if strcmp(mode, ‘per’)。 t=［1:N］./F_s+sqrt(6)*SM./F_s*(rand(1,N))。 t=tt(1)。 elseif strcmp(mode, ‘a(chǎn)dd’)。 for ii=2:N。 t(ii)=t(ii1)+1./F_s+sqrt(12)*SM./F_s*()。 end end elseif strcmp(option, ‘norm’)。if strcmp(mode, ‘per’)。 t=［1:N］./F_s+sqrt(2)*SM./F_s*(rand(1,N))。 t=tt(1)。 elseif strcmp(mode, ‘a(chǎn)dd’)。 for ii=2:N。 t(ii)=t(ii1)+1./F_s+SM./F_s*(randn)。 end end endif nargin4t=t+t0。e