生物醫(yī)學信號處理ppt課件

1、l 5.1 生物醫(yī)學信號的特點生物醫(yī)學信號的特點l 5.2 生物醫(yī)學信號的檢測處理方法概述生物醫(yī)學信號的檢測處理方法概述l 5.3 信號及其描述信號及其描述l 5.4 信號處理的一般方法信號處理的一般方法l 5.5 應(yīng)用實例應(yīng)用實例l 生物醫(yī)學信號處理是生物醫(yī)學工程學的生物醫(yī)學信號處理是生物醫(yī)學工程學的一個重要研究領(lǐng)域，也是近年來迅速發(fā)一個重要研究領(lǐng)域，也是近年來迅速發(fā)展的數(shù)字信號處理技術(shù)的一個重要的應(yīng)展的數(shù)字信號處理技術(shù)的一個重要的應(yīng)用方面，正是由于數(shù)字信號處理技術(shù)和用方面，正是由于數(shù)字信號處理技術(shù)和生物醫(yī)學工程的緊密結(jié)合，才使得我們生物醫(yī)學工程的緊密結(jié)合，才使得我們在生物醫(yī)學信號特征的檢測

2、、提取及臨在生物醫(yī)學信號特征的檢測、提取及臨床應(yīng)用上有了新的手段，因而也幫助我床應(yīng)用上有了新的手段，因而也幫助我們加深了對人體自身的認識。們加深了對人體自身的認識。l人體中每時每刻都存在著大量的生命信人體中每時每刻都存在著大量的生命信息。由于我們的身體整個生命過程中都息。由于我們的身體整個生命過程中都在不斷地實現(xiàn)著物理的、化學的及生物在不斷地實現(xiàn)著物理的、化學的及生物的變化，因此所產(chǎn)生的信息是極其復(fù)雜的變化，因此所產(chǎn)生的信息是極其復(fù)雜的。的。l我們可以把生命信號概括分為二大類：我們可以把生命信號概括分為二大類：l化學信息化學信息 l物理信息物理信息 l化學信息是指組成人體的有機物在發(fā)生變化時化

3、學信息是指組成人體的有機物在發(fā)生變化時所給出的信息，它屬于生物化學所研究的范疇。所給出的信息，它屬于生物化學所研究的范疇。l物理信息是指人體各器官運動時所產(chǎn)生的信息。物理信息是指人體各器官運動時所產(chǎn)生的信息。物理信息所表現(xiàn)出來的信號又可分為電信號和物理信息所表現(xiàn)出來的信號又可分為電信號和非電信號兩大類。非電信號兩大類。l人體電信號，如體表心電人體電信號，如體表心電ECG信號、腦電信號、腦電EEG）、肌電）、肌電EMG）、眼電）、眼電EOG）、）、胃電胃電EGG等在臨床上取得了不同程度的應(yīng)等在臨床上取得了不同程度的應(yīng)用。人體磁場信號檢測近年來也引起了國內(nèi)外用。人體磁場信號檢測近年來也引起了國內(nèi)外

4、研究者和臨床的高度重視，我們把磁場信號也研究者和臨床的高度重視，我們把磁場信號也可歸為人體電信號?？蓺w為人體電信號。l人體非電信號，如體溫、血壓、心音、心輸出人體非電信號，如體溫、血壓、心音、心輸出量及肺潮氣量等，通過相應(yīng)的傳感器，即可轉(zhuǎn)量及肺潮氣量等，通過相應(yīng)的傳感器，即可轉(zhuǎn)變成電信號。變成電信號。l電信號是最便于檢測、提取和處理的信號。電信號是最便于檢測、提取和處理的信號。 l上述信號是由人體自發(fā)生產(chǎn)的，稱為上述信號是由人體自發(fā)生產(chǎn)的，稱為 “主動性信號。主動性信號。l另外，還有一種另外，還有一種“被動性信號，即人體被動性信號，即人體在外界施加某種刺激或某種物質(zhì)時所產(chǎn)在外界施加某種刺激或某

5、種物質(zhì)時所產(chǎn)生的信號。如誘發(fā)響應(yīng)信號，即是在刺生的信號。如誘發(fā)響應(yīng)信號，即是在刺激下所產(chǎn)生的電信號，在超聲波及激下所產(chǎn)生的電信號，在超聲波及X 射射線作用下所產(chǎn)生的人體各部位的超聲圖線作用下所產(chǎn)生的人體各部位的超聲圖象、象、X 射線圖象等也是一種被動信號。射線圖象等也是一種被動信號。這些信號是我們進行臨床診斷的重要工這些信號是我們進行臨床診斷的重要工具。具。 l我們這里所研究的生物醫(yī)學信號即是上我們這里所研究的生物醫(yī)學信號即是上述的包括主動的、被動的、電的和非電述的包括主動的、被動的、電的和非電的人體物理信息。的人體物理信息。l1信號弱信號弱 l2噪聲強噪聲強 l3頻率范圍一般較低頻率范圍一般

6、較低 l4隨機性強隨機性強 l1信號弱：直接從人體中檢測到的生理信號弱：直接從人體中檢測到的生理電信號其幅值一般比較小。如從母體腹電信號其幅值一般比較小。如從母體腹部取到的胎兒心電信號僅為部取到的胎兒心電信號僅為1050V，腦干聽覺誘發(fā)響應(yīng)信號小于腦干聽覺誘發(fā)響應(yīng)信號小于1V，自發(fā)，自發(fā)腦電信號約腦電信號約5150V，體表心電信號相，體表心電信號相對較大，最大可達對較大，最大可達5mV。l 因此，在處理各種生理信號之前要配置因此，在處理各種生理信號之前要配置各種高性能的放大器。各種高性能的放大器。l2噪聲強：噪聲是指其它信號對所研究噪聲強：噪聲是指其它信號對所研究對象信號的干擾。如電生理信號總

7、是伴對象信號的干擾。如電生理信號總是伴隨著由于肢體動作、精神緊張等帶來的隨著由于肢體動作、精神緊張等帶來的干擾，而且?；煊休^強的工頻干擾；誘干擾，而且?；煊休^強的工頻干擾；誘發(fā)腦電信號中總是伴隨著較強的自發(fā)腦發(fā)腦電信號中總是伴隨著較強的自發(fā)腦電；從母腹取到的胎兒心電信號常被較電；從母腹取到的胎兒心電信號常被較強的母親心電所淹沒。這給信號的檢測強的母親心電所淹沒。這給信號的檢測與處理帶來了困難。與處理帶來了困難。l 因此要求采用一系列的有效的去除噪聲因此要求采用一系列的有效的去除噪聲的算法。的算法。l3頻率范圍一般較低：經(jīng)頻譜分析可知，頻率范圍一般較低：經(jīng)頻譜分析可知，除聲音信號如心音頻譜成分較

8、高外，除聲音信號如心音頻譜成分較高外，其它電生理信號的頻譜一般較低。如心其它電生理信號的頻譜一般較低。如心電的頻譜為電的頻譜為0.0135Hz，腦電的頻譜分，腦電的頻譜分布在布在l30Hz之間。之間。l 因此在信號的獲取、放大、處理時要充因此在信號的獲取、放大、處理時要充分考慮對信號的頻率響應(yīng)特性。分考慮對信號的頻率響應(yīng)特性。l4隨機性強：生物醫(yī)學信號是隨機信號，一隨機性強：生物醫(yī)學信號是隨機信號，一般不能用確定的數(shù)學函數(shù)來描述，它的規(guī)律主般不能用確定的數(shù)學函數(shù)來描述，它的規(guī)律主要從大量統(tǒng)計結(jié)果中呈現(xiàn)出來，必須借助統(tǒng)計要從大量統(tǒng)計結(jié)果中呈現(xiàn)出來，必須借助統(tǒng)計處理技術(shù)來檢測、辨識隨機信號和估計它

9、的特處理技術(shù)來檢測、辨識隨機信號和估計它的特征。而且它往往是非平穩(wěn)的，即信號的統(tǒng)計特征。而且它往往是非平穩(wěn)的，即信號的統(tǒng)計特征如均值、方差等隨時間的變化而改變。征如均值、方差等隨時間的變化而改變。 這給生物醫(yī)學信號的處理帶來了困難。這給生物醫(yī)學信號的處理帶來了困難。l 因此在信號處理時往往進行相應(yīng)的理想化和因此在信號處理時往往進行相應(yīng)的理想化和簡化。當信號非平穩(wěn)性變化不太快時，可以把簡化。當信號非平穩(wěn)性變化不太快時，可以把它作為分段平穩(wěn)的準平穩(wěn)信號來處理；如果信它作為分段平穩(wěn)的準平穩(wěn)信號來處理；如果信號具有周期重復(fù)的節(jié)律性，只是周期和各周期號具有周期重復(fù)的節(jié)律性，只是周期和各周期的波形有一定程

10、度的隨機變異，則可以作為周的波形有一定程度的隨機變異，則可以作為周期平穩(wěn)的重復(fù)性信號來處理。更一般性的方法期平穩(wěn)的重復(fù)性信號來處理。更一般性的方法是采用自適應(yīng)處理技術(shù)，使處理的參數(shù)自動跟是采用自適應(yīng)處理技術(shù)，使處理的參數(shù)自動跟隨信號的非平穩(wěn)性而改變。隨信號的非平穩(wěn)性而改變。 1研究不同生物醫(yī)學信號檢測和提取的研究不同生物醫(yī)學信號檢測和提取的方法；方法； 2研究突出信號本身、抑制或除去噪聲研究突出信號本身、抑制或除去噪聲的各種算法；的各種算法； 3研究對不同信號的特征的提取算法；研究對不同信號的特征的提取算法； 4研究信號特征在臨床上的應(yīng)用。研究信號特征在臨床上的應(yīng)用。l5.2.1 生物醫(yī)學信號

11、檢測方法生物醫(yī)學信號檢測方法l5.2.2 生物醫(yī)學信號處理方法生物醫(yī)學信號處理方法l5.2.3 數(shù)字信號處理的特點數(shù)字信號處理的特點l生物醫(yī)學信號檢測是對生物體中包含的生物醫(yī)學信號檢測是對生物體中包含的生命現(xiàn)象、形狀、性質(zhì)和成分等信息進生命現(xiàn)象、形狀、性質(zhì)和成分等信息進行檢測和量化的技術(shù)。行檢測和量化的技術(shù)。l涉及到人機接口技術(shù)、低噪聲和抗干擾涉及到人機接口技術(shù)、低噪聲和抗干擾技術(shù)、信號拾取、分析與處理技術(shù)等工技術(shù)、信號拾取、分析與處理技術(shù)等工程領(lǐng)域，也依賴于生命科學如細胞生程領(lǐng)域，也依賴于生命科學如細胞生理、神經(jīng)生理等研究的進展。理、神經(jīng)生理等研究的進展。l信號檢測一般需要通過以下步驟：信號

12、檢測一般需要通過以下步驟：l生物醫(yī)學信號通過電極拾取或通過傳感生物醫(yī)學信號通過電極拾取或通過傳感器轉(zhuǎn)換成電信號，經(jīng)放大器及預(yù)處理器器轉(zhuǎn)換成電信號，經(jīng)放大器及預(yù)處理器進行信號放大和預(yù)處理，然后經(jīng)進行信號放大和預(yù)處理，然后經(jīng)A/D轉(zhuǎn)轉(zhuǎn)換器進行采樣，將模擬信號轉(zhuǎn)變?yōu)閿?shù)字換器進行采樣，將模擬信號轉(zhuǎn)變?yōu)閿?shù)字信號，輸入計算機，然后通過各種數(shù)字信號，輸入計算機，然后通過各種數(shù)字信號處理算法進行信號分析處理，得到信號處理算法進行信號分析處理，得到有意義的結(jié)果。有意義的結(jié)果。 心電電極、心音傳感器、導聯(lián)線心電電極、心音傳感器、導聯(lián)線心電、心音信號放大器心電、心音信號放大器數(shù)據(jù)采集卡數(shù)據(jù)采集卡A/D轉(zhuǎn)換卡）轉(zhuǎn)換卡

13、）l生物醫(yī)學傳感器是獲取生物醫(yī)學信息并將其轉(zhuǎn)生物醫(yī)學傳感器是獲取生物醫(yī)學信息并將其轉(zhuǎn)換成易于測量和處理的信號的關(guān)鍵器件。生物換成易于測量和處理的信號的關(guān)鍵器件。生物醫(yī)學信號檢測技術(shù)的研究已涉及生物體各層次醫(yī)學信號檢測技術(shù)的研究已涉及生物體各層次的廣泛的生物信息。的廣泛的生物信息。l應(yīng)用電極可檢測心電、腦電、肌電、眼電和神應(yīng)用電極可檢測心電、腦電、肌電、眼電和神經(jīng)電等各種生物電信號；物理傳感器已用于血經(jīng)電等各種生物電信號；物理傳感器已用于血壓、血流、體溫，心音、脈搏、呼吸等各種生壓、血流、體溫，心音、脈搏、呼吸等各種生理量的測量；應(yīng)用化學傳感器可檢測血、尿等理量的測量；應(yīng)用化學傳感器可檢測血、尿

14、等體液中多種離子濃度；用于檢測酶、抗原、抗體液中多種離子濃度；用于檢測酶、抗原、抗體、神經(jīng)遞質(zhì)、激素、受體、體、神經(jīng)遞質(zhì)、激素、受體、DNA和和 RNA等等生物活性物質(zhì)的生物傳感器亦在研究及迅速發(fā)生物活性物質(zhì)的生物傳感器亦在研究及迅速發(fā)展之中；心磁、腦磁等生物磁信號的檢測方法展之中；心磁、腦磁等生物磁信號的檢測方法的研究正在受到重視。的研究正在受到重視。l生物醫(yī)學信號檢測技術(shù)已廣泛應(yīng)用于醫(yī)學研究、生物醫(yī)學信號檢測技術(shù)已廣泛應(yīng)用于醫(yī)學研究、臨床檢查、病人監(jiān)護、治療控制、以及人工器臨床檢查、病人監(jiān)護、治療控制、以及人工器官和運動醫(yī)學等領(lǐng)域，是一種基礎(chǔ)性技術(shù)。官和運動醫(yī)學等領(lǐng)域，是一種基礎(chǔ)性技術(shù)。l

15、生物醫(yī)學信號檢測的發(fā)展趨向是：發(fā)展微型化、生物醫(yī)學信號檢測的發(fā)展趨向是：發(fā)展微型化、多參數(shù)生物醫(yī)學傳感器，特別是加強化學傳感多參數(shù)生物醫(yī)學傳感器，特別是加強化學傳感器和生物傳感器的實用化研究；發(fā)展以生物電器和生物傳感器的實用化研究；發(fā)展以生物電和生物磁為代表的無創(chuàng)檢測技術(shù)；發(fā)展植入式、和生物磁為代表的無創(chuàng)檢測技術(shù)；發(fā)展植入式、動態(tài)監(jiān)測式技術(shù)和生物遙測技術(shù)；發(fā)展細胞和動態(tài)監(jiān)測式技術(shù)和生物遙測技術(shù)；發(fā)展細胞和分子水平的檢測技術(shù)。分子水平的檢測技術(shù)。l生物醫(yī)學信號處理是研究從被干擾和噪生物醫(yī)學信號處理是研究從被干擾和噪聲淹沒的信號中提取有用的生物醫(yī)學信聲淹沒的信號中提取有用的生物醫(yī)學信息的特征并作模

16、式分類的方法。息的特征并作模式分類的方法。l由于生物醫(yī)學信號具有隨機性強和噪聲背景強由于生物醫(yī)學信號具有隨機性強和噪聲背景強的特點，采用了諸多數(shù)字處理技術(shù)進行分析：的特點，采用了諸多數(shù)字處理技術(shù)進行分析：l如對信號時域分析的相干平均算法、相關(guān)技術(shù)；如對信號時域分析的相干平均算法、相關(guān)技術(shù)；l對信號頻域分析的快速傅立葉變換算法、各種對信號頻域分析的快速傅立葉變換算法、各種數(shù)字濾波算法；數(shù)字濾波算法；l對平穩(wěn)隨機信號分析的功率譜估計算法、參數(shù)對平穩(wěn)隨機信號分析的功率譜估計算法、參數(shù)模型方法；模型方法；l對非平穩(wěn)隨機信號分析的短時傅立葉變換、時對非平穩(wěn)隨機信號分析的短時傅立葉變換、時頻分布維格納分布

17、）、小波變換、時變參數(shù)頻分布維格納分布）、小波變換、時變參數(shù)模型、自適應(yīng)處理等算法；模型、自適應(yīng)處理等算法；l對信號的非線性處理方法如混沌與分形、人工對信號的非線性處理方法如混沌與分形、人工神經(jīng)網(wǎng)絡(luò)算法等。神經(jīng)網(wǎng)絡(luò)算法等。 l這些方法在生物醫(yī)學信號分析、醫(yī)學圖像技術(shù)和醫(yī)學這些方法在生物醫(yī)學信號分析、醫(yī)學圖像技術(shù)和醫(yī)學儀器中已得到了廣泛的應(yīng)用。例如：儀器中已得到了廣泛的應(yīng)用。例如：l采用相干平均技術(shù)已成功提取誘發(fā)腦電、希氏束電和采用相干平均技術(shù)已成功提取誘發(fā)腦電、希氏束電和心室晚電位等微弱信號；心室晚電位等微弱信號；l在心電和腦電體表標測中采用計算機進行多道信號同在心電和腦電體表標測中采用計算機

18、進行多道信號同步處理并推求原始信號源的活動逆問題）；步處理并推求原始信號源的活動逆問題）；l在心電、腦電、心音、肺音等信號的自動識別分析中在心電、腦電、心音、肺音等信號的自動識別分析中應(yīng)用了多種信號處理方法進行特征提取與自動分類；應(yīng)用了多種信號處理方法進行特征提取與自動分類；l在生理信號數(shù)據(jù)壓縮和模式分類中引入了人工神經(jīng)網(wǎng)在生理信號數(shù)據(jù)壓縮和模式分類中引入了人工神經(jīng)網(wǎng)絡(luò)方法；絡(luò)方法；l在腦電、心電、神經(jīng)電活動、圖像分割處理、三維圖在腦電、心電、神經(jīng)電活動、圖像分割處理、三維圖像表面特征提取及建模等方面引入混沌與分形理論等，像表面特征提取及建模等方面引入混沌與分形理論等，已取得了許多重要的研究成

19、果并得到了廣泛的臨床應(yīng)已取得了許多重要的研究成果并得到了廣泛的臨床應(yīng)用。用。 l自自1960年以來，隨著計算機技術(shù)和現(xiàn)代年以來，隨著計算機技術(shù)和現(xiàn)代信息技術(shù)的飛速發(fā)展，產(chǎn)生了一門新的信息技術(shù)的飛速發(fā)展，產(chǎn)生了一門新的獨立學科體系：數(shù)字信號處理獨立學科體系：數(shù)字信號處理Digital Signal Processing, DSP）。）。l數(shù)字信號處理是利用計算機或?qū)Ｓ锰幚頂?shù)字信號處理是利用計算機或?qū)Ｓ锰幚硇酒詳?shù)值計算的方法對信號進行采芯片，以數(shù)值計算的方法對信號進行采集、分析、變換、識別等加工處理，從集、分析、變換、識別等加工處理，從而達到提取信息和便于應(yīng)用的目的。而達到提取信息和便于應(yīng)用的

20、目的。l數(shù)字信號處理技術(shù)主要是通過計算機算數(shù)字信號處理技術(shù)主要是通過計算機算法進行數(shù)值計算，與傳統(tǒng)的模擬信號處法進行數(shù)值計算，與傳統(tǒng)的模擬信號處理相比，具有如下特點：理相比，具有如下特點：l（1算法靈活算法靈活 l（2運算精確運算精確 l（3抗干擾性強抗干擾性強 l（4速度快速度快 l此外，數(shù)字系統(tǒng)還具有設(shè)備尺寸小，造此外，數(shù)字系統(tǒng)還具有設(shè)備尺寸小，造價低，便于大規(guī)模集成，便于實現(xiàn)多維價低，便于大規(guī)模集成，便于實現(xiàn)多維信號處理等突出優(yōu)點。信號處理等突出優(yōu)點。l在生物醫(yī)學信號處理領(lǐng)域，數(shù)字信號處在生物醫(yī)學信號處理領(lǐng)域，數(shù)字信號處理技術(shù)發(fā)揮著極其重要的作用。理技術(shù)發(fā)揮著極其重要的作用。l 5.3.

21、1 信號信號l 5.3.2 系統(tǒng)系統(tǒng)l 5.3.3 采樣采樣l信號信號(Signal)可以描述范圍極為廣泛的一類物可以描述范圍極為廣泛的一類物理現(xiàn)象。理現(xiàn)象。l在信號處理領(lǐng)域，信號被定義為一個隨時間變在信號處理領(lǐng)域，信號被定義為一個隨時間變化的物理量，例如心電監(jiān)護儀描記的病人的心化的物理量，例如心電監(jiān)護儀描記的病人的心電、呼吸等信號。電、呼吸等信號。l信號一般可以表示為一個數(shù)學函數(shù)式，以信號一般可以表示為一個數(shù)學函數(shù)式，以x(t)表示，自變量表示，自變量t為時間，為時間，x(t)表示信號隨時間表示信號隨時間t的的變化情況。如正弦波信號：變化情況。如正弦波信號：l一個實際信號除了用函數(shù)式表示外，

22、還常常用一個實際信號除了用函數(shù)式表示外，還常常用曲線來表示。曲線來表示。)sin()(tAtxl信號分類：信號分類：l（1按信號取值的確定性與否按信號取值的確定性與否 ，分為：，分為：l確定性信號：確定性信號：x(t)可確切的表示成時間的可確切的表示成時間的函數(shù)函數(shù) l周期信號：周期信號：T為周期，為周期，n是任是任意整數(shù)意整數(shù) l非周期信號非周期信號l隨機信號：不能確定在某一給定時間的隨機信號：不能確定在某一給定時間的確切取值確切取值 l平穩(wěn)隨機信號平穩(wěn)隨機信號l非平穩(wěn)隨機信號非平穩(wěn)隨機信號)()(tnTxtxl（2按信號的時間取值特點，分為：按信號的時間取值特點，分為：l連續(xù)時間信號連續(xù)時

23、間信號 l離散時間信號離散時間信號 l如果如果t是定義在時間軸上的連續(xù)變化的量，稱是定義在時間軸上的連續(xù)變化的量，稱x(t)為連續(xù)時間信號連續(xù)信號），或模擬信為連續(xù)時間信號連續(xù)信號），或模擬信號。即連續(xù)信號是隨時間連續(xù)變化的，在一個號。即連續(xù)信號是隨時間連續(xù)變化的，在一個時間區(qū)間內(nèi)的任何瞬間都有確定的值。時間區(qū)間內(nèi)的任何瞬間都有確定的值。l如果如果t僅在時間軸上的離散點上取值，稱僅在時間軸上的離散點上取值，稱x(t)為為離散時間信號離散信號）。即離散信號只在離散時間信號離散信號）。即離散信號只在離散的時間點有確定的值。一般離散時間信號離散的時間點有確定的值。一般離散時間信號記為記為x(n)，n

24、取整數(shù)，這樣取整數(shù)，這樣x(n)表示為僅是整數(shù)表示為僅是整數(shù)n的函數(shù)，因此的函數(shù)，因此x(n)又稱為離散時間序列序又稱為離散時間序列序列）。如果列）。如果x(n)在幅度上也取離散值，即在時在幅度上也取離散值，即在時間和幅度上都取離散值的信號稱為數(shù)字信號。間和幅度上都取離散值的信號稱為數(shù)字信號。 l一般來說，離散信號的產(chǎn)生有兩種形式：一般來說，離散信號的產(chǎn)生有兩種形式：一種是信號源本身產(chǎn)生的就是離散信號；一種是信號源本身產(chǎn)生的就是離散信號；而更多的情況下，離散信號是通過對連而更多的情況下，離散信號是通過對連續(xù)信號的采樣得到的，例如在對病人監(jiān)續(xù)信號的采樣得到的，例如在對病人監(jiān)護時每隔半小時測一次體

25、溫，雖然病人護時每隔半小時測一次體溫，雖然病人的體溫是連續(xù)變化的，但采樣輸出的是的體溫是連續(xù)變化的，但采樣輸出的是離散信號，在一天內(nèi)得到離散信號，在一天內(nèi)得到48個采樣值，個采樣值，構(gòu)成了一個離散的體溫信號。構(gòu)成了一個離散的體溫信號。 0sinnl單位沖激函數(shù)與單位脈沖序列單位沖激函數(shù)與單位脈沖序列 ：l連續(xù)時間單位沖激函數(shù)連續(xù)時間單位沖激函數(shù)(t)定義為：定義為：l且且t0時，時，(t)0。 l離散時間單位脈沖序列離散時間單位脈沖序列(n)定義為：定義為： 1)(dtt0001)(nnnl系統(tǒng)系統(tǒng)System是指由若干相互作用和是指由若干相互作用和相互依賴的事物組合而成的具有某種特相互依賴的

26、事物組合而成的具有某種特定功能的整體。定功能的整體。l在信號處理領(lǐng)域，我們把系統(tǒng)定義為物在信號處理領(lǐng)域，我們把系統(tǒng)定義為物理器件的集合，它在受到輸入信號的激理器件的集合，它在受到輸入信號的激勵時，會產(chǎn)生輸出信號。輸入信號又稱勵時，會產(chǎn)生輸出信號。輸入信號又稱為激勵，輸出信號又稱為響應(yīng)。為激勵，輸出信號又稱為響應(yīng)。l對數(shù)字信號處理，系統(tǒng)可以抽象成一種對數(shù)字信號處理，系統(tǒng)可以抽象成一種變換，或一種運算，將輸入序列變換，或一種運算，將輸入序列x(n)變換變換成輸出序列成輸出序列y(n)。l對系統(tǒng)對系統(tǒng)T，輸入，輸入x(t)時輸出是時輸出是y(t)，我們稱，我們稱y(t)是是系統(tǒng)系統(tǒng)T對對x(t)的響

27、應(yīng)的響應(yīng)(Response)。l當輸入是單位沖激信號當輸入是單位沖激信號 時，系統(tǒng)的輸時，系統(tǒng)的輸出稱為系統(tǒng)的單位沖激響應(yīng)，用出稱為系統(tǒng)的單位沖激響應(yīng)，用h(t)表示。表示。h(t)反映了系統(tǒng)反映了系統(tǒng)T的固有的本質(zhì)，若系統(tǒng)的固有的本質(zhì)，若系統(tǒng)T是線性是線性時不變系統(tǒng)，只要知道了時不變系統(tǒng)，只要知道了h(t)，那么對于任意，那么對于任意的輸入的輸入x(t)，都可以通過公式求出其輸出：，都可以通過公式求出其輸出：l ll該公式稱為卷積積分。該公式稱為卷積積分。)()()()()()()(thtxhtxdthxty)()(ttxl對離散信號，當輸入是單位脈沖信號對離散信號，當輸入是單位脈沖信號(n

28、)時，時，系統(tǒng)的輸出稱為單位脈沖響應(yīng)系統(tǒng)的輸出稱為單位脈沖響應(yīng)h(n)。l例如，體表心電信號是心臟的電活動通過人體例如，體表心電信號是心臟的電活動通過人體傳到體表，通過電極拾取后得到的心電圖信號。傳到體表，通過電極拾取后得到的心電圖信號。心臟是心電圖的信號源，即心臟是心電圖的信號源，即x(t)，人體可以看，人體可以看作是一個容積導體，該導體可看作是系統(tǒng)作是一個容積導體，該導體可看作是系統(tǒng)T，x(t)經(jīng)過系統(tǒng)經(jīng)過系統(tǒng)T后的輸出后的輸出y(t)即是體表心電信號即是體表心電信號y(t)。然而由于人體這一容積導體對心電傳輸。然而由于人體這一容積導體對心電傳輸來說是非線性的，目前還無法得出該系統(tǒng)的單來說

29、是非線性的，目前還無法得出該系統(tǒng)的單位沖激響應(yīng)位沖激響應(yīng)h(t)。 人們也正在研究如何通過體人們也正在研究如何通過體表電位標測由反卷積來求解心電信號源的特征。表電位標測由反卷積來求解心電信號源的特征。l采樣是完成由連續(xù)時間信號轉(zhuǎn)換為離散時間信號的工采樣是完成由連續(xù)時間信號轉(zhuǎn)換為離散時間信號的工具，采樣一般由具，采樣一般由A/D轉(zhuǎn)換器完成。轉(zhuǎn)換器完成。A/D轉(zhuǎn)換器就如同一轉(zhuǎn)換器就如同一個電子開關(guān)，如果設(shè)定采樣頻率個電子開關(guān)，如果設(shè)定采樣頻率F，則開關(guān)每隔，則開關(guān)每隔T=1/F秒短暫閉合一次，將連續(xù)信號接通，得到一個離散點秒短暫閉合一次，將連續(xù)信號接通，得到一個離散點的采樣值。假設(shè)開關(guān)每次閉合的時

30、間為的采樣值。假設(shè)開關(guān)每次閉合的時間為秒秒 2倍的信號最高頻率倍的信號最高頻率Adequately sampledAliased due to undersamplingl 5.4.1 相干平均算法相干平均算法l 5.4.2 相關(guān)技術(shù)相關(guān)技術(shù)l 5.4.3 頻域分析技術(shù)頻域分析技術(shù)l 5.4.4 信號的濾波信號的濾波l 5.4.5 參數(shù)模型參數(shù)模型l在計算機普及應(yīng)用之前，信號處理裝置在計算機普及應(yīng)用之前，信號處理裝置或系統(tǒng)都是由模擬器件和電路組成，如或系統(tǒng)都是由模擬器件和電路組成，如RLC電路。這種系統(tǒng)的輸入輸出信號都電路。這種系統(tǒng)的輸入輸出信號都是模擬信號。隨著計算機和數(shù)字信號處是模擬信號。

31、隨著計算機和數(shù)字信號處理算法的發(fā)展，數(shù)字信號處理得到了飛理算法的發(fā)展，數(shù)字信號處理得到了飛速發(fā)展。我們在此簡要介紹一下數(shù)字信速發(fā)展。我們在此簡要介紹一下數(shù)字信號處理的基本方法。號處理的基本方法。l相干平均相干平均(Coherent Average)主要應(yīng)用于主要應(yīng)用于能多次重復(fù)出現(xiàn)的信號的提取。如果待能多次重復(fù)出現(xiàn)的信號的提取。如果待檢測的醫(yī)學信號與噪聲重疊在一起，信檢測的醫(yī)學信號與噪聲重疊在一起，信號如果可以重復(fù)出現(xiàn)，而噪聲是隨機信號如果可以重復(fù)出現(xiàn)，而噪聲是隨機信號，可用疊加法提高信噪比，從而提取號，可用疊加法提高信噪比，從而提取有用的信號。有用的信號。l疊加方法：按固定周期或固定觸發(fā)時刻

32、疊加方法：按固定周期或固定觸發(fā)時刻進行疊加。進行疊加。 l效果估計：效果估計：l l其中其中yi(t)為含有噪聲的待檢測信號，為含有噪聲的待檢測信號，其中其中s(t)為重復(fù)出現(xiàn)的有用信號，為重復(fù)出現(xiàn)的有用信號，ni(t)為為隨機噪聲。經(jīng)隨機噪聲。經(jīng)N次疊加后求平均，那么次疊加后求平均，那么 ：)()()(tntstyiiNiNiiNiiitnNtstntsNtyNty111)(1)()()(1)(1)(l若信號若信號s(t)的功率為的功率為P，噪聲，噪聲ni(t)的方的方差為差為2，那么對每一個，那么對每一個yi(t)，其信噪，其信噪比為比為P/ 2 。經(jīng)。經(jīng)N次平均后，噪聲的方次平均后，噪聲

33、的方差變?yōu)椴钭優(yōu)? /N，所以平均后信號的信噪，所以平均后信號的信噪比為比為N P/ 2 ，提高了，提高了N倍。倍。l例如心室晚電位為例如心室晚電位為V級，掩埋在噪聲級，掩埋在噪聲里，如按心動周期以里，如按心動周期以R峰點對齊，進峰點對齊，進行疊加、平均，則可檢出微弱的心室行疊加、平均，則可檢出微弱的心室晚電位信號。晚電位信號。l 信號的相關(guān)函數(shù)反映了兩個信號之間的信號的相關(guān)函數(shù)反映了兩個信號之間的相互關(guān)連的程度。相互關(guān)連的程度。l設(shè)有兩個信號設(shè)有兩個信號x(n)和和y(n)，定義它們的互相關(guān)，定義它們的互相關(guān)函數(shù)函數(shù)Across-correlation Funtionrxy為：為：ll它表示

34、它表示x(n)不動，將不動，將y(n)在時間軸上左移或右在時間軸上左移或右移移m為正數(shù)時左移，為正數(shù)時左移，m為負數(shù)時右移為負數(shù)時右移m個個時間間隔后分別與時間間隔后分別與x(n)逐點對應(yīng)相乘后求和，逐點對應(yīng)相乘后求和，得到該得到該m點時刻的相關(guān)函數(shù)值點時刻的相關(guān)函數(shù)值rxy(m)。以。以m為為橫軸，橫軸，rxy(m)為縱軸可畫出相關(guān)函數(shù)曲線，該為縱軸可畫出相關(guān)函數(shù)曲線，該曲線反映了曲線反映了x(n)和和y(n)的相似程度。的相似程度。nxymmnynxmr取任意整數(shù)),()()(l 一個信號一個信號x(n)的自相關(guān)函數(shù)的自相關(guān)函數(shù)l（Autocorrelation Functionrxx定義

35、為：定義為：ll其中，其中，rxx(0)反映了信號反映了信號x(n)自身的能量。自身的能量。rxx(m)是偶函數(shù)，是偶函數(shù)，rxx(0)是其中的最大值。是其中的最大值。自相關(guān)函數(shù)曲線可反映信號自身的周期性自相關(guān)函數(shù)曲線可反映信號自身的周期性和噪聲水平。和噪聲水平。 nxxmmnxnxmr取任意整數(shù)),()()(l相關(guān)技術(shù)應(yīng)用范圍很廣，例如，我們可以利用相關(guān)技術(shù)應(yīng)用范圍很廣，例如，我們可以利用相關(guān)判斷在一個含有噪聲的記錄中有無我們所相關(guān)判斷在一個含有噪聲的記錄中有無我們所希望的信號。設(shè)記錄到的信號：希望的信號。設(shè)記錄到的信號：ll其中其中s(n)為信號，為信號，(n)為白噪聲白噪聲是指為白噪聲白

36、噪聲是指其頻譜為一非零常數(shù)的噪聲），現(xiàn)在我們不知其頻譜為一非零常數(shù)的噪聲），現(xiàn)在我們不知道當前記錄到的道當前記錄到的y(n)中是否存在中是否存在s(n)，但我們，但我們根據(jù)以前的工作已知道關(guān)于根據(jù)以前的工作已知道關(guān)于s(n)的先驗知識，的先驗知識，因此我們可以做因此我們可以做y(n)與與s(n)的互相關(guān)：的互相關(guān)：ll通常我們認為信號與白噪聲是不相關(guān)的，因此通常我們認為信號與白噪聲是不相關(guān)的，因此rs(m)等于零，于是等于零，于是rys(m) = rss(m) 。因此我。因此我們可以根據(jù)互相關(guān)函數(shù)們可以根據(jù)互相關(guān)函數(shù)rys(m)與自相關(guān)函數(shù)與自相關(guān)函數(shù)rss(m)是否相等來判斷在是否相等來判斷

37、在y(n)中是否含有信號中是否含有信號s(n)。 )()()(nnsny)()()(mrmrmrsssysl對于信號對于信號x(t)或或x(n)，我們可以在時域直，我們可以在時域直接對其進行分析，如濾波、求相關(guān)函數(shù)、接對其進行分析，如濾波、求相關(guān)函數(shù)、相干平均、特征提取等，然而，對信號相干平均、特征提取等，然而，對信號特征的深入研究，往往轉(zhuǎn)換到頻域進行特征的深入研究，往往轉(zhuǎn)換到頻域進行分析，有助于加深對信號特征的認識。分析，有助于加深對信號特征的認識。l頻域分析的一個典型應(yīng)用即是對信號進頻域分析的一個典型應(yīng)用即是對信號進行傅立葉變換，研究信號所包含的各種行傅立葉變換，研究信號所包含的各種頻率成

38、分。頻率成分。 l我們知道，對于一個周期信號，如正弦波信號：我們知道，對于一個周期信號，如正弦波信號：y=sin(t)，具有一個單一的頻譜值，具有一個單一的頻譜值。而對于任。而對于任意一個周期信號意一個周期信號f(t)都可用傅立葉級數(shù)表示為：都可用傅立葉級數(shù)表示為：l其中，其中，l即任何一個周期函數(shù)都可以展開成為頻率值為即任何一個周期函數(shù)都可以展開成為頻率值為基頻基頻和其和其m次倍頻次倍頻m的三角函數(shù)和的形式，的三角函數(shù)和的形式，系數(shù)系數(shù)am即為信號即為信號f(t)所包含的該頻率成分的頻譜。所包含的該頻率成分的頻譜。 mtjmmeatf)(tmjtmetjmsincosl進一步推廣，若取實際的

39、有限長離散采樣信號進一步推廣，若取實際的有限長離散采樣信號x(n)，可以將該有限長信號看作是周期信號的，可以將該有限長信號看作是周期信號的一個基本周期，同樣可以應(yīng)用傅立葉級數(shù)理論，一個基本周期，同樣可以應(yīng)用傅立葉級數(shù)理論，計算計算x(n)的頻譜，得到離散傅立葉變換公式：的頻譜，得到離散傅立葉變換公式：l應(yīng)用該公式計算離散傅立葉變換有一個快速算應(yīng)用該公式計算離散傅立葉變換有一個快速算法，這就是著名的快速傅立葉變換法，這就是著名的快速傅立葉變換FFT）。）。 10,)()()(102NkenxnxDFTkXNnnkNjl傅立葉變換只能對確定性信號進行分析，而隨機信號在時間上是無限的，在樣本上是無窮

40、多，其傅立葉變換不存在，因此，對隨機信號只能計算信號的功率譜。信號的功率譜可以由信號的相關(guān)函數(shù)計算得到：ll l 因此，只要我們能求出信號的相關(guān)函數(shù)rxx(m)，即可求出信號x(n)的功率譜。但是，真正的rxx(m)也很難求出，要靠由x(n)估計出來，這就是功率譜估計。功率譜估計在生物醫(yī)學信號處理中應(yīng)用極為廣泛，如在心電、心音、腦電等處理中取得了良好的效果。mmjxxjemreP)()(l在對醫(yī)學信號進行檢測分析時，由于生在對醫(yī)學信號進行檢測分析時，由于生物醫(yī)學信號總是不可避免的伴隨著不同物醫(yī)學信號總是不可避免的伴隨著不同頻率的噪聲干擾，為了有效的提取信號，頻率的噪聲干擾，為了有效的提取信號，

41、抑制噪聲，需要使用相應(yīng)的濾波器進行抑制噪聲，需要使用相應(yīng)的濾波器進行濾波。數(shù)字濾波器是數(shù)字信號處理中使濾波。數(shù)字濾波器是數(shù)字信號處理中使用的最廣泛的一種線性系統(tǒng)，是數(shù)字信用的最廣泛的一種線性系統(tǒng)，是數(shù)字信號處理的重要基礎(chǔ)。號處理的重要基礎(chǔ)。l數(shù)字濾波器作為一個線性系統(tǒng)，系統(tǒng)的數(shù)字濾波器作為一個線性系統(tǒng)，系統(tǒng)的輸入輸入x(t)包含信號包含信號s(t)和干擾和干擾n(t)，如果，如果s(t)和和n(t)在頻譜上不重疊，即可通過一個特在頻譜上不重疊，即可通過一個特定的濾波器系統(tǒng)濾除干擾定的濾波器系統(tǒng)濾除干擾n(t)，得到的輸，得到的輸出出y(t)近似地等于近似地等于s(t)。l濾波器有四種頻率響應(yīng)特

42、性：濾波器有四種頻率響應(yīng)特性：l低通、高通、帶通、帶阻低通、高通、帶通、帶阻l 低通濾波器有一個截止頻率低通濾波器有一個截止頻率0，只允許頻率，只允許頻率低于低于0的頻率成分通過，凡是輸入信號中頻率的頻率成分通過，凡是輸入信號中頻率成分高于成分高于0的均被濾除，因此在輸出信號的均被濾除，因此在輸出信號y(t)中只含有低于中只含有低于0的頻率成分。的頻率成分。l高通濾波器正好相反，只允許頻率高于高通濾波器正好相反，只允許頻率高于0的頻率成分通過。的頻率成分通過。l帶通濾波器有上下邊帶截止頻率帶通濾波器有上下邊帶截止頻率1和和2，只允許只允許12的頻率成分通過。的頻率成分通過。l帶阻濾波器阻止帶阻

43、濾波器阻止12的頻率成分通過。的頻率成分通過。l低通低通l高通高通l帶通帶通l參數(shù)模型技術(shù)是近年來得到迅速發(fā)展的參數(shù)模型技術(shù)是近年來得到迅速發(fā)展的新技術(shù)，受到人們的普遍關(guān)注。在對隨新技術(shù)，受到人們的普遍關(guān)注。在對隨機過程的研究中，由于我們不能得到一機過程的研究中，由于我們不能得到一個隨機過程的完整描述，只能得到它們個隨機過程的完整描述，只能得到它們有限次的有限長的觀察值，因此我們希有限次的有限長的觀察值，因此我們希望能用一個數(shù)學模型來模擬我們要研究望能用一個數(shù)學模型來模擬我們要研究的隨機過程，使該模型的輸出等于或近的隨機過程，使該模型的輸出等于或近似該過程。似該過程。l我們用一個白噪聲我們用一

44、個白噪聲u(n)作為輸入去激勵一個系作為輸入去激勵一個系統(tǒng)統(tǒng)h(n)，得到輸出，得到輸出x(n)，如果滿足：，如果滿足：ll該系統(tǒng)稱為自回歸模型該系統(tǒng)稱為自回歸模型AR模型或線性模型或線性預(yù)測模型，其物理意義是：模型現(xiàn)在的輸出預(yù)測模型，其物理意義是：模型現(xiàn)在的輸出x(n)是由現(xiàn)在的輸入是由現(xiàn)在的輸入u(n)和過去的和過去的p個輸出的線個輸出的線性加權(quán)得到。只要我們能求出系數(shù)，即可確定性加權(quán)得到。只要我們能求出系數(shù)，即可確定模型參數(shù)。模型參數(shù)。l通過該模型，可以完成很多有價值的通過該模型，可以完成很多有價值的研究工作，例如，可以估計信號研究工作，例如，可以估計信號x(n)的功率譜、的功率譜、進行

45、各種特征分析、判別分析等工作。進行各種特征分析、判別分析等工作。)()()(1nuknxanxpkkl應(yīng)用應(yīng)用AR模型估計信號的功率譜模型估計信號的功率譜l已知采樣信號已知采樣信號x(n)l我們用一個白噪聲我們用一個白噪聲u(n)作為輸入去激勵一作為輸入去激勵一個系統(tǒng)個系統(tǒng)h(n)，使其能夠得到輸出，使其能夠得到輸出x(n)，建立系，建立系統(tǒng)的統(tǒng)的AR模型：模型：l 若可以求出模型的系數(shù)若可以求出模型的系數(shù)ak和常數(shù)和常數(shù)2w，則可，則可用下式計算信號的功率譜：用下式計算信號的功率譜：)()()(1nuknxanxpkk2121)(pkkjkxxeaPl以上我們簡要介紹了生物醫(yī)學信號處理以上我

46、們簡要介紹了生物醫(yī)學信號處理的幾種基本方法。當然，信號處理的內(nèi)的幾種基本方法。當然，信號處理的內(nèi)容非常豐富，例如多采樣率信號處理、容非常豐富，例如多采樣率信號處理、非平穩(wěn)信號的時頻分布、同態(tài)濾波、自非平穩(wěn)信號的時頻分布、同態(tài)濾波、自適應(yīng)濾波、小波變換、人工神經(jīng)網(wǎng)絡(luò)、適應(yīng)濾波、小波變換、人工神經(jīng)網(wǎng)絡(luò)、混沌與分形等方法，在生物醫(yī)學信號處混沌與分形等方法，在生物醫(yī)學信號處理領(lǐng)域得到了廣泛的應(yīng)用。理領(lǐng)域得到了廣泛的應(yīng)用。l 5.5.1 心電信號的計算機分析心電信號的計算機分析l 5.5.2 心音信號的計算機分析心音信號的計算機分析 分析步驟分析步驟心電數(shù)據(jù)采集：心電數(shù)據(jù)采集：500Hz采樣頻率采樣頻率

47、心電信號預(yù)處理：濾除干擾基線漂移、心電信號預(yù)處理：濾除干擾基線漂移、50Hz、肌電、肌電）特征點檢測：特征點檢測：P、QRS、ST、T波波自動診斷：心律失常分析與波形分類自動診斷：心律失常分析與波形分類lQRS波形檢測算法：波形檢測算法：l經(jīng)典的經(jīng)典的QRS波檢測算法包括三部分波檢測算法包括三部分 ：l 線性濾波；線性濾波；l 非線性變換；非線性變換；l 決策規(guī)則。決策規(guī)則。l線性濾波一般采用中心頻率在線性濾波一般采用中心頻率在1025Hz之之間帶寬為間帶寬為510Hz的帶通濾波器，用于減除的帶通濾波器，用于減除ECG信號中的非信號中的非QRS波頻率成分，提高信噪比。波頻率成分，提高信噪比。l非線性變換的目的是將每個非線性變換的目的是將每個QRS波信號變波信號變換為單向正波峰。換為單向正波峰。l決策規(guī)則一般用峰值檢測器或自適應(yīng)閾值決策規(guī)則一般用峰值檢測器或自適應(yīng)閾值檢測器來檢測檢測器來檢測QRS波。波。 l由于心電信號是非平穩(wěn)隨機信號，利用由于心電信號是非平穩(wěn)隨機信號，利用經(jīng)典的經(jīng)典的QRS波檢測算法往往受到以下兩波檢測算法往往受到以下兩個因素的制約：個因素的制約：l QRS復(fù)合波的信號帶寬對于不同的復(fù)合波的信號帶寬對于不同的病人乃至同一病人的不同心搏均有所不病人乃至同一病人的不同心搏均有所不同。同。 各種噪聲與各種噪聲與QRS波的通帶相互交波的通帶相互交疊。疊。