數(shù)據(jù)采集及處理演示文稿

數(shù)據(jù)采集及處理演示文稿現(xiàn)在是1頁\一共有104頁\編輯于星期五數(shù)據(jù)采集及處理現(xiàn)在是2頁\一共有104頁\編輯于星期五5.1數(shù)字信號處理概述

對連續(xù)信號而言，隨著數(shù)字處理技術的發(fā)展，越來越迫切地要求連續(xù)信號的離散化。看似連續(xù)的信號是可以由其離散的樣本值來表征的。采樣——從連續(xù)時間信號中提取離散樣本的過程。抽樣即時間軸上離散化的過程。抽樣若按抽樣間隔來分，可分為均勻抽樣與非均勻抽樣。我們討論的是均勻抽樣。在沒有任何約束的條件下，離散時間樣本不能唯一地表示連續(xù)時間信號。因為有無限多個信號都可以產(chǎn)生一組給定的樣本值。一個連續(xù)時間信號必須在某一種條件下才能由其樣本來表示。現(xiàn)在是3頁\一共有104頁\編輯于星期五1、數(shù)字信號處理的主要研究內(nèi)容

數(shù)字信號處理主要研究用數(shù)字序列來表示測試信號，并用數(shù)學公式和運算來對這些數(shù)字序列進行處理。內(nèi)容包括數(shù)字波形分析、幅值分析、頻譜分析和數(shù)字濾波。0AtX(0)X(1)X(2)X(3)X(4)現(xiàn)在是4頁\一共有104頁\編輯于星期五2、測試信號數(shù)字化處理的基本步驟

物理信號對象傳感器電信號放大調(diào)制電信號A/D轉(zhuǎn)換數(shù)字信號計算機顯示D/A轉(zhuǎn)換電信號控制物理信號現(xiàn)在是5頁\一共有104頁\編輯于星期五3、數(shù)字信號處理的優(yōu)勢1)用數(shù)學計算和計算機顯示代替復雜的電路和機械結(jié)構(gòu)現(xiàn)在是6頁\一共有104頁\編輯于星期五2)計算機軟硬件技術發(fā)展的有力推動a)多種多樣的工業(yè)用計算機。

現(xiàn)在是7頁\一共有104頁\編輯于星期五b)靈活、方便的計算機虛擬儀器開發(fā)系統(tǒng)現(xiàn)在是8頁\一共有104頁\編輯于星期五5.2模數(shù)(A/D)、數(shù)模(D/A)與采樣定理采樣――利用采樣脈沖序列，從信號中抽取一系列離散值，使之成為采樣信號x(nTs)的過程．編碼――將經(jīng)過量化的值變?yōu)槎M制數(shù)字的過程。

量化――把采樣信號經(jīng)過舍入變?yōu)橹挥杏邢迋€有效數(shù)字的數(shù)，這一過程稱為量化．1、A/D轉(zhuǎn)換

現(xiàn)在是9頁\一共有104頁\編輯于星期五4位A/D:XXXXX(1)0101X(2)0011X(3)0000現(xiàn)在是10頁\一共有104頁\編輯于星期五2)A/D轉(zhuǎn)換器的技術指標

(3)模擬信號的輸入范圍;

如，5V，+/-5V，10V，+/-10V等。(1)分辨率;

用輸出二進制數(shù)碼的位數(shù)表示。位數(shù)越多，量化誤差越小，分辨力越高。常用有8位、10位、12位、16位等。(2)轉(zhuǎn)換速度;

指完成一次轉(zhuǎn)換所用的時間，如:1ms(1KHz)；10us(100kHz)現(xiàn)在是11頁\一共有104頁\編輯于星期五2、D/A轉(zhuǎn)換過程和原理

D/A轉(zhuǎn)換器是把數(shù)字信號轉(zhuǎn)換為電壓或電流信號的裝置。

D/A轉(zhuǎn)換器的技術指標

分辨率;轉(zhuǎn)換速度;模擬信號的輸出范圍;現(xiàn)在是12頁\一共有104頁\編輯于星期五A/D、D/A轉(zhuǎn)換過程中的量化誤差實驗：

現(xiàn)在是13頁\一共有104頁\編輯于星期五離散時間信號可認為是對模擬信號的采樣?！嘤斜匾懻摬蓸雍?，（1）信號的內(nèi)容是否有丟失？

（2）信號的頻譜有何變化？

（3）采樣信號能否恢復為原信號？

3采樣定理

現(xiàn)在是14頁\一共有104頁\編輯于星期五顯然，采樣間隔T的選取是很重要的。直觀地看，T太大，則易丟失信息，使一些細節(jié)內(nèi)容無法反映出來，而T太小，顯然信息不易丟失，但卻使數(shù)據(jù)量明顯增大，處理費時。—采樣定理。∴應適當?shù)剡x擇現(xiàn)在是15頁\一共有104頁\編輯于星期五（1）正弦波采樣定理現(xiàn)在是16頁\一共有104頁\編輯于星期五X(0),X(1),X(2),……,X(n)

現(xiàn)在是17頁\一共有104頁\編輯于星期五每周期應該有多少采樣點？最少3點:三個方程求解三個變量現(xiàn)在是18頁\一共有104頁\編輯于星期五現(xiàn)在是19頁\一共有104頁\編輯于星期五

一個頻率有限信號如果頻譜只占據(jù) 的范圍，則信號可以用等間隔的抽樣值來唯一地表示。而抽樣間 隔不大于（其中），或者 說最低抽樣頻率為。

奈奎斯特頻率：（2）時域采樣定理現(xiàn)在是20頁\一共有104頁\編輯于星期五時域理想抽樣的傅立葉變換相乘相卷時域抽樣頻域周期重復現(xiàn)在是21頁\一共有104頁\編輯于星期五時域理想抽樣的傅立葉變換FTFT相乘

相卷積現(xiàn)在是22頁\一共有104頁\編輯于星期五周期矩形被沖激抽樣的頻譜先重復后抽樣現(xiàn)在是23頁\一共有104頁\編輯于星期五先抽樣時域抽樣頻域重復后重復時域重復頻域抽樣現(xiàn)在是24頁\一共有104頁\編輯于星期五

非理想抽樣信號的傅立葉變換乘卷現(xiàn)在是25頁\一共有104頁\編輯于星期五關于非理想抽樣現(xiàn)在是26頁\一共有104頁\編輯于星期五由抽樣信號恢復原連續(xù)信號取主頻帶：時域卷積定理：現(xiàn)在是27頁\一共有104頁\編輯于星期五卷積包絡相乘現(xiàn)在是28頁\一共有104頁\編輯于星期五（3）頻域抽樣定理

若信號為時限信號，它集中在的時間范圍內(nèi)，若在頻域中，以不大于的頻率間隔對的頻譜進行抽樣，則抽樣后的頻譜可以唯一地表示原信號?，F(xiàn)在是29頁\一共有104頁\編輯于星期五根據(jù)時域和頻域?qū)ΨQ性，可推出頻域抽樣定理偶函數(shù)變量置換現(xiàn)在是30頁\一共有104頁\編輯于星期五頻域抽樣后的時間函數(shù)相乘卷積現(xiàn)在是31頁\一共有104頁\編輯于星期五現(xiàn)在是32頁\一共有104頁\編輯于星期五抽樣定理小結(jié)時域?qū)Τ闃拥刃в陬l域?qū)χ貜?時域抽樣間隔不大于。頻域?qū)Τ闃拥刃в跁r域?qū)χ貜?頻域抽樣間隔不大于。滿足抽樣定理，則不會產(chǎn)生混疊。信號在時域以T離散化，→頻譜以1/T周期化，且幅度乘以1/T。時，若使采樣信號通過理想低通濾波器：即可恢復原信號。現(xiàn)在是33頁\一共有104頁\編輯于星期五

需注意，滿足采樣定理，只保證不發(fā)生頻率混疊，而不能保證此時的采樣信號能真實地反映原信號x(t)。工程實際中采樣頻率通常大于信號中最高頻率成分的3到5倍?，F(xiàn)在是34頁\一共有104頁\編輯于星期五

以上我們討論的都是滿足抽樣定理要求的情況。

如果對帶限信號抽樣時，抽樣頻率不夠高或抽樣間隔過大，就會出現(xiàn)頻譜的混迭，這一現(xiàn)象就稱為欠抽樣。欠抽樣使信號發(fā)生了頻譜的交叉。

但欠抽樣并不是百害而無一利的，在實際應用中，利用欠抽樣可使高頻變化的信息映射到低頻變化的信號上。這為高頻信號的測量帶來了便利。如頻閃儀和抽樣示波器等。4頻率混疊及防止

現(xiàn)在是35頁\一共有104頁\編輯于星期五不滿足抽樣定理時產(chǎn)生頻率混疊現(xiàn)象現(xiàn)在是36頁\一共有104頁\編輯于星期五頻混計算：

FsFsFsFs頻混正常Fs/2工程處理：混迭頻率=Fs-信號頻率現(xiàn)在是37頁\一共有104頁\編輯于星期五若各次諧波調(diào)制頻譜將相互交疊，觀察采樣信號的頻譜，可看到處就象一面鏡子，將信號頻譜超過的部分反射回來，造成混疊，故也將稱為折疊頻率。

現(xiàn)在是38頁\一共有104頁\編輯于星期五實驗：

現(xiàn)在是39頁\一共有104頁\編輯于星期五頻混現(xiàn)象實驗：

現(xiàn)在是40頁\一共有104頁\編輯于星期五實際應用中，由于信號的時長總是有限的，時域有限→頻域無限。但高頻成分隨頻率的增高總是逐漸衰減的，∴采樣時，（1）近似地忽略的頻率分量；（2）選取采樣率；（3）為避免高于折疊頻率的雜散頻率分量進入采樣器造成頻譜混疊，常在采樣器前面加一個保護性的前置低通濾波器，以阻止高于的頻率分量進入采樣器。工程近似現(xiàn)在是41頁\一共有104頁\編輯于星期五今后在討論數(shù)字信號時，均假定已做過上述處理?，F(xiàn)在是42頁\一共有104頁\編輯于星期五A/D采樣前的抗混迭濾波：

物理信號對象傳感器電信號放大調(diào)制電信號A/D轉(zhuǎn)換數(shù)字信號展開低通濾波(0-Fs/2)放大現(xiàn)在是43頁\一共有104頁\編輯于星期五5.3信號的截斷、能量泄漏

為便于數(shù)學處理，對截斷信號做周期延拓，得到虛擬的無限長信號。

用計算機進行測試信號處理時，不可能對無限長的信號進行測量和運算，而是取其有限的時間片段進行分析，這個過程稱信號截斷?，F(xiàn)在是44頁\一共有104頁\編輯于星期五

周期延拓后的信號與真實信號是不同的，下面我們就從數(shù)學的角度來看這種處理帶來的誤差情況。

設有余弦信號x(t),用矩形窗函數(shù)w(t)與其相乘，得到截斷信號:y(t)=x(t)w(t)將截斷信號譜XT(ω)與原始信號譜X(ω)相比較可知，它已不是原來的兩條譜線，而是兩段振蕩的連續(xù)譜.原來集中在f0處的能量被分散到兩個較寬的頻帶中去了，這種現(xiàn)象稱之為頻譜能量泄漏?，F(xiàn)在是45頁\一共有104頁\編輯于星期五周期延拓信號與真實信號是不同的：能量泄漏誤差現(xiàn)在是46頁\一共有104頁\編輯于星期五能量泄漏實驗：

現(xiàn)在是47頁\一共有104頁\編輯于星期五克服方法之一：信號整周期截斷現(xiàn)在是48頁\一共有104頁\編輯于星期五

能量泄漏分主瓣泄漏和旁瓣泄漏，主瓣泄漏可以減小因柵欄效應帶來的譜峰幅值估計誤差，有其好的一面，而旁瓣泄漏則是完全有害的?，F(xiàn)在是49頁\一共有104頁\編輯于星期五泄露產(chǎn)生的問題:

1）降低了譜分析的頻率分辨力。2）引入了虛假的頻率分量。理想窗函數(shù)應具有的特點：1）主瓣寬度要小，即帶寬要窄。2）旁瓣高度與主瓣高度相比要小，且衰減要快?，F(xiàn)在是50頁\一共有104頁\編輯于星期五3常用的窗函數(shù)

1）矩形窗

現(xiàn)在是51頁\一共有104頁\編輯于星期五2）三角窗現(xiàn)在是52頁\一共有104頁\編輯于星期五3）漢寧窗現(xiàn)在是53頁\一共有104頁\編輯于星期五常用窗函數(shù)現(xiàn)在是54頁\一共有104頁\編輯于星期五5.4DFT與FFT

1、離散傅立葉變換

離散傅里葉變換（DiscreteFourierTransform）一詞是為適應計算機作傅里葉變換運算而引出的一個專用名詞。

x(t)截斷、周期延拓xT(t)周期信號xT(t)的傅里葉變換：現(xiàn)在是55頁\一共有104頁\編輯于星期五

對周期信號xT(t)采樣，得離散序列xT(n),將積分轉(zhuǎn)為集合：展開，得連續(xù)傅立葉變換計算公式：

現(xiàn)在是56頁\一共有104頁\編輯于星期五

采樣信號頻譜是一個連續(xù)頻譜，不可能計算出所有頻率點值，設頻率取樣間隔為：Δf=fs/N

頻率取樣點為｛0,Δf,2Δf,3Δf,....｝，有：

該公式就是離散傅立葉計算公式(DFT)現(xiàn)在是57頁\一共有104頁\編輯于星期五DFT的詳細解釋

傅里葉變換的四種基本形式1．連續(xù)時間與連續(xù)頻率的傅里葉變換

—

連續(xù)傅里葉變換

現(xiàn)在是58頁\一共有104頁\編輯于星期五連續(xù)非周期信號的傅立葉變換例：從傅立葉積分得到從周期信號取單脈沖得到現(xiàn)在是59頁\一共有104頁\編輯于星期五2．連續(xù)時間與離散頻率的傅里葉變換

—

傅里葉級數(shù)……現(xiàn)在是60頁\一共有104頁\編輯于星期五例：周期重復連續(xù)周期信號的傅立葉變換現(xiàn)在是61頁\一共有104頁\編輯于星期五3．離散時間與連續(xù)頻率的傅里葉變換

—

序列傅里葉變換現(xiàn)在是62頁\一共有104頁\編輯于星期五例：從非周期信號抽樣得到離散非周期序列相乘相卷頻域周期重復時域抽樣現(xiàn)在是63頁\一共有104頁\編輯于星期五4．離散時間與離散頻率的傅里葉變換

—

離散傅里葉級數(shù)…………現(xiàn)在是64頁\一共有104頁\編輯于星期五例：從連續(xù)周期信號的抽樣得到現(xiàn)在是65頁\一共有104頁\編輯于星期五第一個域離散函數(shù)?第二個域周期函數(shù)連續(xù)函數(shù)非周期函數(shù)?且易證：一個域中的周期函數(shù)的周期現(xiàn)在是66頁\一共有104頁\編輯于星期五5.4.3從離散傅立葉級數(shù)(DFS)到離散傅立葉變換(DFT)效仿連續(xù)周期信號有傅立葉級數(shù)，記作：離散周期序列也有傅立葉級數(shù)，記作：周期性以N為周期現(xiàn)在是67頁\一共有104頁\編輯于星期五離散周期序列的傅立葉級數(shù)(DFS)的正負運算對周期序列的基頻是是K次諧波分量，諧波系數(shù)是諧波成分中只有N個是獨立的，是周期的現(xiàn)在是68頁\一共有104頁\編輯于星期五現(xiàn)在是69頁\一共有104頁\編輯于星期五有限長序列是周期序列的一個周期有限長序列x(n)只有的N個值x(n)可看成是周期序列的主值序列，記作周期序列當叫做的主值周期，記作有限長序列的以N為周期的周期延拓現(xiàn)在是70頁\一共有104頁\編輯于星期五

的主值序列

也是周期性的，相當于有限長序列周期延拓當時，其主值序列 相當于一個有限長序列現(xiàn)在是71頁\一共有104頁\編輯于星期五

和都取主值周期，得到離散傅立葉變換(DFT)對現(xiàn)在是72頁\一共有104頁\編輯于星期五周期為N和周期為2N的不同當主值周期為0—N-1時，N點的DFT為當主值周期為0—2N-1時，2NDFT（接下頁）現(xiàn)在是73頁\一共有104頁\編輯于星期五現(xiàn)在是74頁\一共有104頁\編輯于星期五小結(jié)

是的主值序列是嚴格按傅立葉分析的概念得來的只是一種借用形式，一種算法用計算信號的頻譜時，采樣頻率必須大于兩倍的信號最高截止頻率對周期信號要取一個整周期現(xiàn)在是75頁\一共有104頁\編輯于星期五DFSDFT周期性以N為周期共軛性以N/2共軛現(xiàn)在是76頁\一共有104頁\編輯于星期五2、快速傅立葉變換

快速傅立葉變換(FFT)是離散傅立葉變換的一種有效的算法，通過選擇和重新排列中間結(jié)果，減小運算量。展開各點的DFT計算公式：XR(1)=x(0).cos(2pi*0*1/N)+x(1).cos(2pi*1*1/N)+x(2).cos(2pi*2*1/N)…..XR(2)=x(0).cos(2pi*0*2/N)+x(1).cos(2pi*1*2/N)+x(2).cos(2pi*2*2/N)…..現(xiàn)在是77頁\一共有104頁\編輯于星期五有大量重復的cos、sin計算，F(xiàn)FT的作用就是用技巧減少cos、sin項重復計算。

當采樣點數(shù)為1024點,DFT要求一百萬次以上計算量，而FFT則只要求一萬次。

現(xiàn)在是78頁\一共有104頁\編輯于星期五5.5柵欄效應

1、柵欄效應

為提高效率,通常采用FFT算法計算信號頻譜，設數(shù)據(jù)點數(shù)為N，采樣頻率為Fs。則計算得到的離散頻率點為:

Xs(Fi)，F(xiàn)i=i*Fs

/N,i=0，1，2，.....，N/2X(f)f0ΔfΔ

如果信號中的頻率分量與頻率取樣點不重合，則只能按四舍五入的原則，取相鄰的頻率取樣點譜線值代替。

現(xiàn)在是79頁\一共有104頁\編輯于星期五柵欄效應誤差實驗：

現(xiàn)在是80頁\一共有104頁\編輯于星期五2能量泄漏與柵欄效應的關系

頻譜的離散取樣造成了柵欄效應，譜峰越尖銳，產(chǎn)生誤差的可能性就越大。

例如，余弦信號的頻譜為線譜。當信號頻率與頻譜離散取樣點不等時，柵欄效應的誤差為無窮大。現(xiàn)在是81頁\一共有104頁\編輯于星期五

實際應用中，由于信號截斷的原因，產(chǎn)生了能量泄漏，即使信號頻率與頻譜離散取樣點不相等，也能得到該頻率分量的一個近似值。從這個意義上說，能量泄漏誤差不完全是有害的。如果沒有信號截斷產(chǎn)生的能量泄漏，頻譜離散取樣造成的柵欄效應誤差將是不能接受的。

現(xiàn)在是82頁\一共有104頁\編輯于星期五

能量泄漏分主瓣泄漏和旁瓣泄漏，主瓣泄漏可以減小因柵欄效應帶來的譜峰幅值估計誤差，有其好的一面，而旁瓣泄漏則是完全有害的。小結(jié)：現(xiàn)在是83頁\一共有104頁\編輯于星期五2、頻率分辨力與頻率分析范圍

頻率分辨力：隨T(N)的增大，分辨力提高。故理論上與fs無關。頻率分析范圍：小于奈奎斯特頻率?，F(xiàn)在是84頁\一共有104頁\編輯于星期五5.6常用的數(shù)字信號處理算法

數(shù)字信號處理是測試技術中最常用和最需要掌握的部分，無論開發(fā)簡單或復雜的測控系統(tǒng)或儀器，都會用到數(shù)字信號處理知識?，F(xiàn)在是85頁\一共有104頁\編輯于星期五時域波形分析是最常用的信號分析手段，用示波器、萬用表等普通儀器顯示信號波形就可以特征參數(shù)。

5.6常用的數(shù)字信號處理算法

1、時域波形參數(shù)計算

0At峰值/雙峰值均值有效值方差周期現(xiàn)在是86頁\一共有104頁\編輯于星期五波形分析的應用超門限報警

信號類型識別

信號基本參數(shù)識別

Pp-p現(xiàn)在是87頁\一共有104頁\編輯于星期五1)峰值P，雙峰值Pp-pP1=data(0)P2=data(0)ForK=0ToNIfP1<data(k)ThenP1=data(k)EndIf

IfP2>data(k)ThenP2=data(k)EndIfNextP=P1P2P=P1-P2PPp-p現(xiàn)在是88頁\一共有104頁\編輯于星期五2)均值U=0ForK=0ToNU=U+data(k)NextU=U/N現(xiàn)在是89頁\一共有104頁\編輯于星期五3)均方值E2=0ForK=0ToNE2=E2+data(k)*data(k)NextRMS=sqr(E2/N)現(xiàn)在是90頁\一共有104頁\編輯于星期五4)方差V2=0ForK=0ToNV2=V2+(data(k)-U)*(data(k)-U)NextV=V2/N大方差

小方差

U=0ForK=0ToNU=U+data(k)NextU=U/N現(xiàn)在是91頁\一共有104頁\編輯于星期五信號幅值報警系統(tǒng)設計實驗：

現(xiàn)在是92頁\一共有104頁\編輯于星期五5)周期Tn=0AT=0.8*PForK=2ToNIfdata(k-1)<ATAnddata(k-2)<ATAnddata(k+1)>ATAnddata(k+2)>ATThenti(n)=Kn=n+1EndIfNextT=(ti(2)-ti(1))*dtAtT

現(xiàn)在是93頁\一共有104頁\編輯于星期五數(shù)字轉(zhuǎn)速表設計實驗：

現(xiàn)在是94頁\一共有104頁\編輯于星期五2、數(shù)字相關函數(shù)計算

變量之間的相依關系稱為相關。信號之間的相似關系稱為相關函數(shù)。0At現(xiàn)在是95頁\一共有104頁\編輯于星期五x(t)y(t)時延器

乘法器

y(t-τ)X(t)y(t-τ)積分器

Rxy(τ)Fori=0ToNr(i)=0Forj=0