華科機(jī)械工程測試信息信號分析課件ch6-01數(shù)字信號分析

*1MEASUREMENTINFORMATIONSIGNALANALYSISINMECHANICALENGINEERING

機(jī)械工程測試?信息?信號分析

機(jī)械科學(xué)與工程學(xué)院機(jī)械電子信息工程系*2課件資料下載：郵箱地址：

密碼：注意下載時不要刪除原始文件*3

第六章數(shù)字信號分析(I)

DFT與FFT*4第六章主要內(nèi)容§6-5現(xiàn)代譜分析方法-最大熵譜估計(jì)§6-3FFT§6-4譜分析與譜估計(jì)§6-2離散傅立葉變換DFT§6-1模擬信號離散化時域采樣定理頻域采樣定理周期序列的離散傅立葉級數(shù)*5

時間函數(shù)頻率函數(shù)連續(xù)時間、連續(xù)頻率—傅里葉變換連續(xù)時間、離散頻率—傅里葉級數(shù)離散時間、連續(xù)頻率—序列的傅里葉變換離散時間、離散頻率—離散傅里葉變換Fourier變換的幾種可能形式*600t時域連續(xù)函數(shù)造成頻域是非周期的譜，而時域的非周期造成頻域是連續(xù)的譜密度函數(shù)。連續(xù)時間、連續(xù)頻率-FT域連續(xù)性周期性時域連續(xù)非周期頻域連續(xù)非周期*7連續(xù)時間、離散頻率-FS當(dāng)x(t)為連續(xù)時間周期信號時，可展開為傅立葉級數(shù)域連續(xù)性周期性時域連續(xù)周期頻域離散非周期時域連續(xù)函數(shù)造成頻域是非周期的譜，頻域的離散對應(yīng)時域是周期函數(shù)。時域周期為T0,頻域譜線間隔為20/T0*8離散時間、連續(xù)頻率--序列的FT對離散序列x(n)，其傅立葉變換為：若x(n)是信號x(t)的采樣序列，采樣間隔為T,則有：*9序列的FT域連續(xù)性周期性時域離散非周期頻域連續(xù)周期時域的離散化造成頻域的周期延拓，而時域的非周期對應(yīng)于頻域的連續(xù)*10上述三種情況至少在一個變換域有積分(連續(xù))，因而不適合進(jìn)行數(shù)字計(jì)算。域連續(xù)性周期性時域離散周期頻域離散周期時域的離散造成頻域的延拓（周期性）。因而頻域的離散也會造成時域的延拓（周期性）。要想在時域和頻域都是離散的，那么兩域必須是周期的。離散傅立葉變換*11對序列的傅立葉變換在頻域上加以離散化，令d=0，從而*12x(n)*13四種形式歸納類型時間函數(shù)頻率函數(shù)關(guān)系傅立葉變換連續(xù)非周期連續(xù)非周期

傅立葉級數(shù)連續(xù)周期(T0)離散(Ω0)非周期

序列傅立葉變換離散(Ts)非周期連續(xù)周期(Ωs)

離散傅立葉變換離散(Ts)周期(T0)離散(Ω0)周期(Ωs)

*14DFT重要性DFT是重要的變換分析有限長序列的有用工具。在信號處理的理論上有重要意義。在運(yùn)算方法上起核心作用，譜分析、卷積、相關(guān)都可以通DFT在計(jì)算機(jī)上實(shí)現(xiàn)。DFT是現(xiàn)代信號處理橋梁DFT要解決兩個問題：一是離散與量化，二是快速運(yùn)算。信號處理DFT(FFT)傅氏變換離散量化*156-1A/D與D/A轉(zhuǎn)換物理信號對象傳感器電信號放大調(diào)制電信號A/D轉(zhuǎn)換數(shù)字信號計(jì)算機(jī)顯示D/A轉(zhuǎn)換電信號控制物理信號顯示記錄信號預(yù)處理信號采集分析計(jì)算x(t)傳感器和測試系統(tǒng)的標(biāo)定和校準(zhǔn)？法定計(jì)量單位，標(biāo)準(zhǔn)信號源*16測試系統(tǒng)基本模型組成傳感器調(diào)理電路 放大器數(shù)據(jù)采集系統(tǒng) ADC系統(tǒng)靈敏度傳感器靈敏度放大器增益ADC靈敏度S1S2S3xxuvyyH1(j)H2(j)模擬部分?jǐn)?shù)字部分測試系統(tǒng)基本模型(測控系統(tǒng)前向通道)*17基本參數(shù)的確定任務(wù) 以系統(tǒng)分辯力、量程為依據(jù)，依次確定分辯力、量程和環(huán)境條件傳感器類型及其靈敏度S1量程、精度系統(tǒng)分辯率ADC的分辯率(Q=1/S3)ADC量程、傳感器輸出范圍放大器增益S2說明先兩頭，后中間多量程，低位A/DS1S2S3xxuvyyH1(j)H2(j)模擬部分?jǐn)?shù)字部分*18動態(tài)性能的確定任務(wù)根據(jù)被測信號的最高頻率fm和允許動態(tài)幅值誤差m%，確定各環(huán)節(jié)的動態(tài)參數(shù)。方法模數(shù)分開，各自預(yù)估模擬部分誤差表達(dá)理想頻率特性*19模擬部分動態(tài)性能的確定1傳感器與放大器均為一階動態(tài)誤差傳感器時間常數(shù)為1放大器時間常數(shù)為2=1/(2fb)合理選擇1和2，使=2fm時，*20模擬部分動態(tài)性能的確定2傳感器為二階，放大器為一階動態(tài)誤差傳感器固有轉(zhuǎn)角頻率為0傳感器阻尼比為(未指明取=0)放大器時間常數(shù)為2=1/(2fb)合理選擇0、和2，使=2fm時，*21數(shù)字部分動態(tài)性能的預(yù)估ADC的轉(zhuǎn)換時間Tc在Tc之內(nèi)，輸入信號的變化的誤差<ADC量化誤差的一半，即Q/2S/H的孔徑時間TAP(aperture

time)和孔徑抖動時間TAJ(Aperture

Uncertainty(jitter))對于TAP，提前啟動“H”對于TAJ校核說明還應(yīng)考慮采樣定理的要求！系統(tǒng)的動態(tài)性能往往受制于傳感器。*22數(shù)字部分動態(tài)性能的預(yù)估孔徑時間：在模擬量輸入通道中，A/D轉(zhuǎn)換器將模擬信號轉(zhuǎn)換成數(shù)字量需要一定的時間，完成一次A/D轉(zhuǎn)換器所需的時間稱為孔徑時間?？讖蕉秳樱嚎讖蕉秳?或稱孔徑誤差)是指采樣與采樣之間孔徑延遲時間的變化，起因是調(diào)制系統(tǒng)時鐘相位時的噪聲，通過對內(nèi)部ADC時鐘抖動和外部采樣時鐘抖動進(jìn)行和方根(root-sum-square)計(jì)算得到孔徑抖動。如果要求測量準(zhǔn)確，數(shù)據(jù)采樣系統(tǒng)必須要有極低的相位噪聲。隨著模擬輸入斜率(dV/dt)的增加，孔徑抖動也增大。一般來講，使用輸入頻率為MHz級的ADC時，時鐘抖動應(yīng)為亞皮秒級?？讖窖舆t：孔徑延遲是指在保持命令發(fā)出之后到ADC采樣保持放大器(SHA)完全打開采樣開關(guān)所需的時間，即ADC采樣發(fā)出命令到采樣實(shí)際開始的時間。有效孔徑延遲時間te包括孔徑延遲和SHA中模擬、數(shù)字傳輸延遲的影響，其值可正可負(fù)。*23誤差極限的預(yù)估任務(wù)按總誤差的限定值，分配各環(huán)節(jié)的允許誤差。(誤差極限的預(yù)估)過程預(yù)分配綜合調(diào)整再分配再綜合直至滿意模型誤差合成方根和法(幾何合成)被測量x的影響不屬測量系統(tǒng)，因而忽略。*246-1A/D轉(zhuǎn)換過程采樣―利用采樣脈沖序列，從信號中抽取一系列離散值，使之成為采樣信號x(nTs)的過程．量化―把采樣信號經(jīng)過舍入變?yōu)橹挥杏邢迋€有效數(shù)字的數(shù)，這一過程稱為量化．編碼―將經(jīng)過量化的值變?yōu)槎M(jìn)制數(shù)字的過程。

*256-1A/D轉(zhuǎn)換過程-例*266-1D/A轉(zhuǎn)換過程*27模擬信號采樣設(shè)有兩個不同頻率的余弦信號，x1(t)=cos2(10)t,x2(t)=cos2(50)t，現(xiàn)以fs=40Hz進(jìn)行采樣，得：什么原因？脈沖序列采樣理想脈沖采樣1、一個連續(xù)信號經(jīng)過理想采樣后，其采樣信號頻譜是一個周期性連續(xù)頻譜，頻譜周期延拓，周期為-mS。幅值被系數(shù)Cn加權(quán)。2、采樣頻率S大于信號頻帶寬度的2倍時，取樣信號的頻譜在每一周期都完整地保留原來模擬信號的信息1.時域采樣：當(dāng)p(t)為脈沖序列時，Cn=1/Ts，所以頻域卷積定理：*29矩形脈沖采樣1.時域采樣矩形脈沖采樣：實(shí)際采樣脈沖采樣脈沖為周期矩形脈沖時沖*30時域采樣定理2.時域采樣定理時域采樣定理：采樣頻率必須大于或等于信號x(t)中的最高頻率的兩倍。物理解釋：頻譜受限信號，如果頻譜只占據(jù)-m—m范圍，則信號可用等間隔采樣值唯一的表示，采樣間隔不大于1/2fm頻帶有限信號，其頻率大小在時域內(nèi)就是它的波形變化速率，波形的最高變化速度將受最高頻率分量的限制。時域采樣定理表明，一個信號在滿足一定條件下，可通過它的時域采樣點(diǎn)準(zhǔn)確地恢復(fù)原來信號的波形。*31頻混現(xiàn)象頻混現(xiàn)象*32頻混現(xiàn)象實(shí)例1.時域采樣采樣信號恢復(fù)實(shí)例：頻率混疊低通濾波低通濾波WR()1/TS=71/TS=14采樣信號恢復(fù)-圖示1.時域采樣從采樣信號頻譜XS()中復(fù)原X()，采用頻域矩形窗函數(shù)H()與XS()相乘，即：實(shí)現(xiàn)方法：將采樣信號xs(t)通過理想低通濾波器，傳輸函數(shù)為H()，在濾波器的輸出端可以得到頻譜X()的連續(xù)信號x(t)。圖示采樣信號恢復(fù)-公式推導(dǎo)1.時域采樣理想低通濾波器傳輸函數(shù)：在采樣信號xs(t)的每個采樣值上畫一個峰值為x(nTS)的Sa型函數(shù)波形，合成的波形就是x(t)。對不在取樣時刻任意點(diǎn)的數(shù)值是無限加權(quán)采樣值的總和，但內(nèi)插函數(shù)是衰減的，濾波不是理想濾波器，該點(diǎn)值只能是附近的一組有限值之和逼近。*35采樣信號恢復(fù)-每步示例由采樣樣值恢復(fù)原連續(xù)信號的過程惠特克波形重構(gòu)法或理想內(nèi)插法。近似逼近時域采樣信號恢復(fù)：每步示例Sa函數(shù)叫內(nèi)插函數(shù)的原因*36頻域采樣定理-圖示若x(t)的頻譜X()被間隔為0的脈沖序列在頻域中采樣，則在時域中等效于x(t)以一定的時間間隔周期延拓，時間周期為T0=2/0，周期信號的頻譜是離散的。信號的時域與頻域呈采樣（離散）與重復(fù)（周期）關(guān)系。

３.頻域采樣定理：圖示*37頻域采樣定理-公式推導(dǎo)

３.頻域采樣定理連續(xù)頻譜X()，對應(yīng)的時間函數(shù)x(t)，被間隔為1的脈沖序列采樣上式可看出：一個時間受限信號的長度為2tm，在頻域采樣間隔f11/2tm條件下，能夠從采樣點(diǎn)集合完全恢復(fù)原信號的頻譜。*38頻域采樣信號恢復(fù)-公式推導(dǎo)

頻域采樣信號恢復(fù)時域窗函數(shù)：以F0=1/2或T0=2代入上式，則頻譜內(nèi)插公式為：上式說明恢復(fù)過程是內(nèi)插過程，需要無限多的采樣點(diǎn)從-+求和，實(shí)際上因所取的樣點(diǎn)有限、精度有限，濾波器也不具有理想頻率特性，信號恢復(fù)是近似恢復(fù)。*39頻域采樣信號恢復(fù)-圖示4.頻域采樣定理P160在頻域中對X()進(jìn)行采樣，等效于x(t)在時域中重復(fù)，采樣間隔不大于1/2tm，用矩形脈沖作選通信號，就可以恢復(fù)原信號x(t)柵欄效應(yīng)：頻域采樣后，只能獲得采樣點(diǎn)的頻率成分，其余的頻率成分一概舍去，就如透過柵欄觀景。*40復(fù)指數(shù)序列復(fù)指數(shù)序列：N為周期的序列x(t),x(n)都是周期性信號，但對連續(xù)時間信號來說，取值可在-+區(qū)間，而且任意選擇0都具有周期性，其周期T=2/0。但對離散時間信號，因?yàn)?，表示在?shù)字頻率上相差2整數(shù)倍的所有離散時間復(fù)指數(shù)序列都是一樣的。即離散域的頻率的有效取值是在02或-的任一間隔為2區(qū)間范圍*5.周期序列的時域分析DFS離散時間周期信號從連續(xù)時間周期信號均勻采樣例：正弦序列例：復(fù)指數(shù)序列TS為時域采樣時間間隔，0=0TS為離散域頻率連續(xù)時間信號FS表達(dá)式離散時間周期信號的時域分析t=nTS離散化*42離散時間周期信號的時域分析5.周期序列的時域分析DFS離散序列FS表達(dá)式：復(fù)指數(shù)形式x(n)是以N為周期的周期序列，表示一個周期的采樣點(diǎn)數(shù)，則：0為離散域的基本頻率，單位是弧度，k0是k次諧波的數(shù)字頻率；TS為采樣時間間隔；N為一個周期采樣點(diǎn)數(shù)。0為周期信號的基波頻率，T0為周期信號的周期。說明：一個周期序列可分解成一系列諧波(k0)關(guān)系的復(fù)指數(shù)序列之和。當(dāng)周期信號從連續(xù)域變到離散域后，頻率也從-+的無限范圍，映射到(0=0TS,0N=2)從0-2的有限范圍內(nèi)，連續(xù)域傅里葉級數(shù)可表示為具有無限多個諧波分量，離散域只含有有限個諧波分量，總共諧波數(shù)為k=2/0=N。*43離散時間周期信號的時域分析5.周期序列的時域分析DFS離散序列FS表達(dá)式：復(fù)指數(shù)形式上式為離散傅里葉級數(shù)DFS，是一個有限項(xiàng)級數(shù)，對離散時間周期信號，若用代表正交函數(shù)集的復(fù)指數(shù)序列來表示，該集合是有限的。N是它的最高頻率分量項(xiàng)，與CFS重要區(qū)別。X(k0)計(jì)算：*44離散時間周期信號的時域分析n為時間離散變量，k為頻率離散變量，2/N為離散域基頻(譜線間隔)。X(k0)是復(fù)指數(shù)序列各諧波分量的復(fù)振幅，反映了各諧波分量的幅度和相位，表示離散時間周期信號的頻譜。將k以k+N代入上式DFS變換對離散時間信號的頻譜。是一個以N(2)為周期的周期性離散頻譜，個譜線間的間隔為0=2/N，且存在著諧波的關(guān)系。*45周期信號DFS-舉例1已知一正弦離散時間信號x(n)=cos0n，分別求出當(dāng)

(1)0=20.5及(2)0=/3時的傅里葉級數(shù)表示式和相應(yīng)頻譜解：(1)正弦序列只有在滿足0/2=m/N=有理數(shù)的條件下才是周期序列，非周期序列，不能展開為DFS。*46周期信號DFS-舉例1已知一正弦離散時間信號x(n)=cos0n，分別求出當(dāng)

(1)0=20.5及(2)0=/3時的傅里葉級數(shù)表示式和相應(yīng)頻譜解：比較*47周期信號DFS-舉例2已知一周期序列x(n)，周期N=6，如圖，求該序列的頻譜及時域表示式*48周期信號DFS-舉例2已知一周期序列x(n)，周期N=6，如圖，求該序列的頻譜及時域表示式*49離散時間周期信號的頻域分析離散時間周期信號與連續(xù)時間周期信號頻譜間的內(nèi)在聯(lián)系舉例1：x(t)=6cost，現(xiàn)以采樣間隔T=0.25秒采樣，求取樣后周期序列的頻譜并與原始信號x(t)的頻譜進(jìn)行比較。解：求得*50離散時間周期信號的頻域分析在一個周期內(nèi)|X(k0)|=|X(k0)|，說明離散周期信號在-<<范圍內(nèi)的離散頻譜，可以準(zhǔn)確地等于連續(xù)時間周期信號的離散頻譜*51離散時間周期信號的頻域分析舉例2：(a)按fs1=16，Ts1=1/16，則每一基本周期取樣點(diǎn)數(shù)N=T0/Ts1=16，則：比較圖(a)和(c)，滿足采樣定理，不會出現(xiàn)頻率混疊*52離散時間周期信號的頻域分析舉例2：(b)按fs2=8，Ts2=1/8，N=T0/Ts2=8，則幅度頻譜為：比較圖(b)和(c)，不滿足采樣定理，出現(xiàn)頻率混疊，f2=5的正弦分量混疊到f1=3的余弦分量上，f1=3混疊到f2=5上。*53頻率混疊與能量泄漏混疊時域欠采樣時，出現(xiàn)頻率混疊無法恢復(fù)原信號頻譜，因而不能從時域采樣點(diǎn)準(zhǔn)確地重建原連續(xù)信號。同理，頻域欠采樣時，出現(xiàn)波形混疊無法恢復(fù)原頻譜對應(yīng)的信號，也不能從頻域采樣值重建原連續(xù)頻譜。改進(jìn)方法：提高采樣速率，增加采樣點(diǎn)數(shù)，減少混疊對頻譜分析的影響*54頻率混疊與能量泄漏舉例混疊：N=8，取樣間隔T=T0/N=0.5N=16，取樣間隔T=T0/N=0.25頻譜圖試分別求出如下圖(a)所示矩形周