DSP实现癫痫脑电信号处理

　　式(2)的含义是：信号f(t)在第2j尺度(第j层)的近似部分的小波系数，即低频部分的小波系数Aj是通过第2j+1尺度(第j+1层)的逼近部分的小波系数Aj+1隔点插零后与重构滤波器h卷积以及第2j+1尺度(第j+1层)的细节部分的小波系数Dj+1隔点插零后与重构滤波器g卷积，然后求和得到的。不断重复这一过程，直到第2°尺度，得到重构信号。

　　3 小波变换的DSP实现

　　3.1 脑电信号在CCS 2.2上的输入与输出

　　CCS 2.2(Code Composer Studio)是由TI公司推出的一种针对标准TMS320调试接口的集成开发环境(IDE)，利用CCS集成开发环境，用户可以完成工程定义、程序编辑、编译链接、调试和数据分析等工作环节[4]。我们把十进制的浮点数用两个十六进制数进行表示，采用C语言实现。

　　再利用CCS中的File->Load Data将十六进制的数据导入到DSP的相应内存中去。

　　反过来，DSP处理之后的数据利用CCS的数据导出File->Save以文本文件形式保存，再用C语言进行数据逆转化，把两个十六进制数进行转化成十进制的浮点数。

　　其中的result数组就是十进制的浮点型，origin数组就是十六进制的浮点型。

3.2 核心汇编程序介绍

　　以下是以16位定点乘法实现32位浮点乘法的部分汇编程序：

　　3.3 实验结果与分析

　　图2(a)为待处理脑电信号，经小波分解后(b)～(f)依次为各级逼近波形图，(g)～(k)则依次为对应的细节波形图。滤掉j=3的细节波形即图2(i)，再进行小波重构后得到图2(1)，从中发现，原图2(a)和重构后的图2(1)几乎看不出明显差别。