DSP實現(xiàn)癲癇腦電信號處理

　　式(2)的含義是：信號f(t)在第2j尺度(第j層)的近似部分的小波系數(shù)，即低頻部分的小波系數(shù)Aj是通過第2j+1尺度(第j+1層)的逼近部分的小波系數(shù)Aj+1隔點插零后與重構濾波器h卷積以及第2j+1尺度(第j+1層)的細節(jié)部分的小波系數(shù)Dj+1隔點插零后與重構濾波器g卷積，然后求和得到的。不斷重復這一過程，直到第2°尺度，得到重構信號。

　　3 小波變換的DSP實現(xiàn)

　　3.1 腦電信號在CCS 2.2上的輸入與輸出

　　CCS 2.2(Code Composer Studio)是由TI公司推出的一種針對標準TMS320調試接口的集成開發(fā)環(huán)境(IDE)，利用CCS集成開發(fā)環(huán)境，用戶可以完成工程定義、程序編輯、編譯鏈接、調試和數(shù)據(jù)分析等工作環(huán)節(jié)[4]。我們把十進制的浮點數(shù)用兩個十六進制數(shù)進行表示，采用C語言實現(xiàn)。

　　再利用CCS中的File->Load Data將十六進制的數(shù)據(jù)導入到DSP的相應內存中去。

　　反過來，DSP處理之后的數(shù)據(jù)利用CCS的數(shù)據(jù)導出File->Save以文本文件形式保存，再用C語言進行數(shù)據(jù)逆轉化，把兩個十六進制數(shù)進行轉化成十進制的浮點數(shù)。

　　其中的result數(shù)組就是十進制的浮點型，origin數(shù)組就是十六進制的浮點型。

3.2 核心匯編程序介紹

　　以下是以16位定點乘法實現(xiàn)32位浮點乘法的部分匯編程序：

　　3.3 實驗結果與分析

　　圖2(a)為待處理腦電信號，經小波分解后(b)～(f)依次為各級逼近波形圖，(g)～(k)則依次為對應的細節(jié)波形圖。濾掉j=3的細節(jié)波形即圖2(i)，再進行小波重構后得到圖2(1)，從中發(fā)現(xiàn)，原圖2(a)和重構后的圖2(1)幾乎看不出明顯差別。