DSP實(shí)現(xiàn)癲癇腦電信號(hào)處理
式(2)的含義是:信號(hào)f(t)在第2j尺度(第j層)的近似部分的小波系數(shù),即低頻部分的小波系數(shù)Aj是通過第2j+1尺度(第j+1層)的逼近部分的小波系數(shù)Aj+1隔點(diǎn)插零后與重構(gòu)濾波器h卷積以及第2j+1尺度(第j+1層)的細(xì)節(jié)部分的小波系數(shù)Dj+1隔點(diǎn)插零后與重構(gòu)濾波器g卷積,然后求和得到的。不斷重復(fù)這一過程,直到第2°尺度,得到重構(gòu)信號(hào)。
3 小波變換的DSP實(shí)現(xiàn)
3.1 腦電信號(hào)在CCS 2.2上的輸入與輸出
CCS 2.2(Code Composer Studio)是由TI公司推出的一種針對(duì)標(biāo)準(zhǔn)TMS320調(diào)試接口的集成開發(fā)環(huán)境(IDE),利用CCS集成開發(fā)環(huán)境,用戶可以完成工程定義、程序編輯、編譯鏈接、調(diào)試和數(shù)據(jù)分析等工作環(huán)節(jié)[4]。我們把十進(jìn)制的浮點(diǎn)數(shù)用兩個(gè)十六進(jìn)制數(shù)進(jìn)行表示,采用C語(yǔ)言實(shí)現(xiàn)。
再利用CCS中的File->Load Data將十六進(jìn)制的數(shù)據(jù)導(dǎo)入到DSP的相應(yīng)內(nèi)存中去。
反過來,DSP處理之后的數(shù)據(jù)利用CCS的數(shù)據(jù)導(dǎo)出File->Save以文本文件形式保存,再用C語(yǔ)言進(jìn)行數(shù)據(jù)逆轉(zhuǎn)化,把兩個(gè)十六進(jìn)制數(shù)進(jìn)行轉(zhuǎn)化成十進(jìn)制的浮點(diǎn)數(shù)。
其中的result數(shù)組就是十進(jìn)制的浮點(diǎn)型,origin數(shù)組就是十六進(jìn)制的浮點(diǎn)型。
3.2 核心匯編程序介紹
以下是以16位定點(diǎn)乘法實(shí)現(xiàn)32位浮點(diǎn)乘法的部分匯編程序:
3.3 實(shí)驗(yàn)結(jié)果與分析
圖2(a)為待處理腦電信號(hào),經(jīng)小波分解后(b)~(f)依次為各級(jí)逼近波形圖,(g)~(k)則依次為對(duì)應(yīng)的細(xì)節(jié)波形圖。濾掉j=3的細(xì)節(jié)波形即圖2(i),再進(jìn)行小波重構(gòu)后得到圖2(1),從中發(fā)現(xiàn),原圖2(a)和重構(gòu)后的圖2(1)幾乎看不出明顯差別。
評(píng)論