数字信号处理在生物医学工程中的应用

收稿日期:2005-09-21;修回时间:2006-03-28

数字信号处理在生物医学工程中的应用

娄 智

(安徽交通职业技术学院,合肥 230051)

摘 要:数字信号处理技术一诞生就显示了强大的生命力,展现了极为广阔的应用前景。主要讨论数字信号处理技术中小波分析、人工神经网络、维格纳分布在生物医学工程中的应用,并对数字信号处理技术在生物医学工程中的应用前景进行了展望。

关键词:数字信号处理;小波分析;人工神经网络;维格纳分布中图分类号:R 318

文献标识码:B

文章编号:1008-9632(2006)04-0038-03

1 引言

自20世纪60年代以来,随着计算机和信息学科

的飞速发展,大量的模拟信息被转化为数字信息来处理。于是就逐步产生了一门近代新兴学科)))数字信号处理(D i g ital Signal Processi n g ,简称DSP)技术。经过几十年的发展,数字信号处理技术现已形成了一门以快速傅里叶变换和数字滤波器为核心,以逻辑电路为基础,以大规模集成电路为手段,利用软硬件来实现各种模拟信号的数字处理,其中包括信号检测、信号变换、信号的调制和解调、信号的运算、信号的传输和信号的交换等各种功能作用的独立的学科体系。

而生物医学工程就是应用物理学和工程学的技术去解决生物系统中所存在的问题,特别是人类疾病的诊断、治疗和预防的科学。它包括工程学、医学和生命科学中的许多学科。本文主要讨论数字信号处理技术中小波分析、人工神经网络、维格纳分布在生物医学工程中的应用。

2 数字信号处理在生物医学工程中的应用211 小波分析在生物医学工程中的应用

[1~4]

近年来,小波的研究受到数学家,理论物理学家和工程学家们的关注,特别是在信号处理,图象处理,语音分析,模式识别,量子物理及众多非线性科学等应用领域,被认为是近年来在工具及方法上的重大突破。

所谓的小波变换是指把某一被称为基本小波(mot h er wave let)的函数W (t)作位移S 后,再在不同尺度A 下与待分析信号x (t)作内积:

WT x (A ,S )=1A

Q

x (t)W *(t-S

A )d t A >0

其中A 为伸缩因子,S 为平移因子。等效的频域表示:

WT x (A ,S )=A 2P

Q

X (X )7*(A X )e +j XS

d X

式中X (X ),7(X )分别是x (t),W (t)的傅里叶变换。

任何变换都必须能进行反变换才有实际意义,但反变换未必一定存在,对小波变换而言,所采用的小波必须满足允许条件C W =

Q

]

|W (X )|2

X

d X <]。由此推

论出7(X )=0或Q

W (t)dt =0,即小波变换必须具有带通性质。本来满足允许条件的W (t)便可作为基本小波,但考虑到频域上的局域要求,条件就更苛刻一些:即要求小波在频域上局域性能好,应要求W (t)的前n 阶矩为零,也就是Q

t n

W (t)d t =0,且n 越大越好。在频域上这相当于要求7(X )在X =0处有n 阶零点。

小波分析方法具有以下特点:(1)时频局部化特点,即可以同时提供时域和频域局部化信息。(2)多分辨率,即多尺度的特点,可以由粗到细逐步观察信号。(3)带通滤波的特点,可以根据中心频率的变化调节带宽,中心频率的高低与带宽成反向变化,可以观测出信号的低频缓变部分和高频突变部分。这种变焦特性决定了它对非平稳信号处理的特殊功能。在生物医学工程中的信号处理,信号压缩,医学图象处理中,小波变换均有用武之地。

适当地选择小波基,可以方便地检测出信号的奇异点,观测信号的瞬态变化以及时域分析中信号不见的信息;此外利用带通特性,将信号分解成不同频带低频分解波和高频分解波,并提取出信号中的非平稳信号。在生物医学工程中的人体电信号,如心电信号(ECG),脑电信号(EEG ),肌电信号(M UAP),视觉诱发电位信号(VEP)等均为非平稳的弱电信号,但对于这些信号的提取常因各信号的频谱相互交迭,以及信噪比较低加之工频及谐波干扰严重等而产生困难,而小波对非平稳信号的突出的处理能力,给人体电信号的提取带来了较以往各种滤波方法更为方便的手段。

38

此外,在CT成象方面,如何既增强边缘又平滑噪声一直是图象处理的难题,而小波变换,由于它可以同时在时域和频域内局部化,因而可以较好地处理图象局部细节,提高边缘分辨度,因此小波方法的出现给生物医学工程中信号图象分析提供了有力的手段。

212人工神经网络(A NN)在生物医学工程中的应用[5,6]人工神经网络是指由大量简单元件(即神经元,可以用电子元件,光学元件等模拟)广泛相互连接而成的复杂网络系统。神经网络有很多具体模型,其共同的基本特征是:(1)以大规模并行处理为主;(2)采用分布式存储,具有很强的容错性和联想功能;(3)强调自适应过程和学习(训练)过程。

人工神经网络的最新发展使其成为信号处理的强有力工具,对于那些用其它信号处理技术无法解决的问题,人工神经网络的应用开辟了新的领域,许多A NN的算法和它们的应用已广泛的在自然科学的各个领域被报道,在这些网络模型中,多层感知器被认为是最有用的学习模型,广泛应用于脑电信号,心电信号的处理中。

20世纪80年代末,90年代初,神经网络的研究在国际上形成一股热潮,其原因是由于神经网络可将人脑的智能原理应用来解决工程技术及社会管理的许多复杂问题。生物医学工程工作者采用神经网络的方法来解释许多复杂的生理、病理现象,例如:心电、脑电、肌电、胃肠电等信号的识别,心电信号的压缩,医学图像的识别和处理等。人工神经网络是由大量的简单处理单元连接而成的自适应动力学系统,具有巨量并行性,分布式存贮,自适应学习的自组织等功能,可以用来解决生物医学信号分析处理中,用常规方法难以解决或无法解决的问题,神经网络在生物医学信号检测与处理中的应用主要集中在对脑电信号的分析,听觉诱发电位信号的提取;用于H olter系统的心电信号数据的压缩算法;医学图像的数据压缩算法等等。这些应用大多数都是基于神经网络的多层前馈网络反向传播算法(即BP算法)训练三层网络,该方法能解决许多信号处理中的难题,如语言合成与识别,视觉模式识别。从输出层开始,连接到第一隐层的连接权用如下算法校正:

$W ij=G D j X j

W ij(t+1)=W ij(t)+G D j X j

其中,W ij(t)对应时刻t输出层i到隐层节点j的连接权,X j隐层第j个节点的输出;G为学习率控制常数,D j 是误差。

由于神经网络可以把专家知识和先验知识结合进一个数学框架来完成提取特征和分类识别等功能,而不需要任何对数据和噪声的先验统计假设,也不需要把专家知识和经验归纳成严密清晰的条文,所以最适应于研究和分析生物医学信号。

213维格纳(W i g ner)分布在生物医学工程中的应用[8,9] 1932年W i g ner描述了一个函数W i g ner分数,为加强时域函数重建为时域图提供了最佳途径。W i g ner分布同时提供信号时域及频域特征,更适合分析随时间变化的能量和。W i g ner分析不需要假设信号是静止的,比FFT及AR分析有更高的分辨率。限制W igner 分布分析应用的不利特征为它只适用于单一成分的信号,如果信号中两种或者多种成分同时存在,函数中将产生伪峰,成为交叉项。对于确定性的连续时间复值函数,其维格纳分布的定义是:

W D x(t,X)=Q]-]x(t+S/2)x*(t-S/2)e-j XS d S 它可以看成是函数r x(t,S)=x(t+S/2)x*(t-S/2)对S的傅氏变换,其等效频域表示是:

W D x(t,X)=(1/2P)Q]-]x(X+N/2)x*(X-N/2)e j N t d N 经过应用数学界几十年的努力,维格纳分布的理论已逐步趋于成熟。进入20世纪80年代以来,许多学者采用维格纳分布对多种非平稳信号进行了分析。由于生物医学信号的非平稳性比较突出,因近年来国内外都有人希望采用维格纳分布来较好地表现它们的频率特性随时间的变化,特别对较微弱的电生理信号。维格纳分布在生物医学信号分析中的应用及发展主要包括生理节律、心电信号、血流速率波信号、体感诱发电位、超声多普勒信号、听神经电活动信号、声音信号、脑电信号、第一心音、心室晚电位、心率,血压和呼吸信号等方面。

目前国内有许多研究机构从事这方面工作。他们的研究对象大部分是高频心电图,心室晚电位和心音信号,已有的结论表明:维格纳分布对这几种信号的特异性和灵敏度较为显著。

维格纳分布在应用中的主要困难是由于它的加法性质而引入的交叉项干扰。有人在利用DP WD对高频心电图和心室晚电位进行分析时,采用了在时轴上作滑动平均的方法来抑制交叉项干扰,取得比较好的效果。另外,对信号进行迭加平均预处理和高通滤波将会取得更好的效果。

维格纳分布在多种非平稳信号分析中所表现的特异性和灵敏度已经得到证明。但到目前为止,各种信号的维格纳分布的分析结果还只停留在定性的特征上,进一步提取能作为诊断的定量特征指标将是今后的主要研究方向。

39