摘要首次提出基于多分辨分析去除食道动态pH监 测信号中常见噪声的算法。文章首先从实验入手对监测信号进行分析，然后结合分析结果依 据基于离散二进小波变换的多分辨分析给出合适的去噪算法，最后通过对临床数据的处理表 明本算法具有较好的去噪效果。

Multiresolution Analysis Based Denoising Algorithm for

Dynamic Esophageal pH Monitoring Signal

Wang LeiZheng ChongxunYuan JunfengSong Lei

(Biomedical Engineering Institute,Xi'an Jiaotong University, Xi'an710 049)

AbstractA denoising algorithm for dynamic esoph ageal pH monitoring signal is first introduced in this paper. The algorithm was developed on the basis of multiresolution analysis. The pH monitoring signal was investigated in experiments, then the algorithm was formulated by multiresoluti on analysis based on discrete dyadic wavelet transform combining with the result s of the investigation mentioned above. Application of this algorithm to clinica l data processing has proved its satisfactory denoising effect.

Key wordsDynamic esophageal pH monitoringMultiresol ution analysisDigital signal processingDenoising algorithm

1监测信号噪声的来源

食道动态pH监测可作为判断胃食道返流的“金标准”^［1］，监测系统借助 传感器和便携式记录器，对患者食道的pH进行长时间(如24 h)的动态监测，然后将记录器内 的数据传输至计算机进行分析。传感器由测量电极和参比电极组成，测量电极一般选用结构 简单、精度较高、成本低的锑电极。根据参比电极的不同可将传感器分为外参比式和内参比 式两种，前者的参比电极为体表贴片电极，在监测过程中贴于体表的合适位置(如胸骨柄处 体表)；后者则是将参比电极与测量电极集成在传感器导管上置于患者食道内，参比电极可 选用柔性Ag-AgCl电极管。在监测过程中会不可避免的引入如下几种噪声：

记录器电噪声。由于空间杂散电磁场能从传感器和电路部分耦合入信号通道，且在模拟电路 处理后仍有一定残留并形成噪声。该噪声可认为是平稳白噪声，同时由于pH传感器所获得的 电信号较强(数百毫伏)，其影响较小。

传感器偶然脱离食道内容物而引起的散点噪声。在监测过程中传感器导管是置于食道内的， 因此当出现电极头脱离食道内容物时会引起测量回路的空载，导致信号迭加上散点状的干扰 噪声，由于目前临床上将食道pH＜4判断为发生胃食道返流，故这种噪声的出现将直接导致 错误的诊断。

外参比式传感器引入的皮肤-粘膜低频漂移噪声^［1，2］。由于内参比式传感的成本 高、对保存和运输的要求高、尚无完善的消毒措施等原因，目前临床上广泛使用的仍为外参 比式传感器。然而，因为其参比电极是皮肤贴片电极，相当于在测量回路内串联了皮肤和食 道粘膜等人体组织，故皮肤电位和跨食道粘膜电位的波动会引入噪声，由于该噪声随人体和 环境而变化，较难定量描述，但实验表明有相当比例的食道动态监测信号出现至少是短时期 的监测失效(指噪声引起的测量误差折合0.5pH以上)。

图1所示为一例搀杂有上述噪声的食道动态pH监测信号，可以看出在许多测量过程中极易引 入的高斯白噪声，对食道动态pH监测信号的影响并不大。

图1具有上述噪声的食道动态pH监测信号

Fig 1Noised dynamic esophageal pH monitorying signal

2监测信号和噪声的实验分析

我们通过模拟实验来分析散点噪声：将传感器置于装有pH＝4缓冲液的试管内静置 1 min，然后在缓冲液中摇动传感器3 min，之后搅动缓冲液约4 min，接着将传感器不断提 出/放入试管4 min，最后恢复初始状态。实验表明散点噪声的分布密度与摇动、晃动的程度 有关，同时也能观察到低幅的记录器电噪声(见图2)。

图2测量回路空载导致的散点噪声

Fig 2Pulse-noise led by the idler of measure loop

在进行外参比食道动态pH监测时，由于外力作用、汗腺分泌、贴片电极与皮 肤的离子交换、粘膜所处生理状态的变化等因素，皮肤电位和跨食道粘膜电位在监测过程中 是不断波动的，我们对此做了大量实验，篇幅所限不能一一列举，这些实验的结果表明电位 波动的特点是低频漂移夹杂局部扰动，漂移幅值最大可达50 mV左右(约折合1pH)，正、负 向漂移均存在，但以正向为多，漂移周期一般在几十分钟。图3为动态监测内、外参比电极 之间电位差的波形片段(图中的阶跃变化由A/D量化误差造成)。

图3动态监测内、外参比电极之间电位差的波形片段

Fig 3Potential difference between inner-reference and outer-r eference electrode during dynamic monitoring

食道动态pH内、外参比同步监测的实验表明，正常食道pH在5～7之间波动， 波动的波幅较小且持续期比漂移噪声长。胃食道返流波一般由较陡的下降沿(返酸过程)和相 对较缓的上升沿(廓酸过程)组成，持续时间以1～10 min之间较多，极小值在pH=2.5左右(即 比胃内的pH稍高)，胃食道返流波表现的是胃内容物返流入食道及食道廓酸机制将胃内容物 从食道中廓清这一生理过程，从生理角度上分析，胃内容物返流入食道的时间较短，即食道 pH会很快由正常生理值(pH=6左右)下降至胃内容物的pH，而之后发生的廓酸过程则持续期较 长，并且由于个体差异导致持续期的变化很大。另外在监测信号中还有可能存在较少见的十 二指肠胃食道返流^［1］，表现为食道pH向碱性漂移(最大可达pH=8左右)，其持续期 与正常食道pH波动类似。

3多分辨分析理论及去噪算法

基于离散二进小波变换的多分辨分析具有在时、频域内均能将信号很好的分解及 重构的优点^［3～5］。设｛V_j｝是一给定的多分辨分析，s和w分别为相应的尺 度和小波函数，再设观测到的实际信号f具有有限的分辨率，即f∈V₀，则存在如下小波变 换分解与重构的迭代方程，其中A₀f(x)=f(x)，h(n)=＜s_-1(u),s(u-n)＞称为由尺 度 函数和小波函数构成的滤波器的脉冲响应；g(n)=＜w_-1(u),s(u-n)＞称为h(n)的镜像 滤波器；h(n)与h^*(n)互为共轭转置。

A_jf(x)和D_jf(x)可理解为在j尺度上的连续逼近和连续细节，f被分解为 频率不超过2^-j的成分和介于2^-j与2^-j+1之间的成分，这相当于通过带 通滤波器组将信号在不同频率通道上展开，且频率越高频带越宽。另外，若w具有对称支集 ［K，-K］，则尺度-空间平面上的锥｜x-x₀｜≤K×2^j叫作X₀的影响锥，锥内点 的小波变换值都受到f(x₀)的影响^［4］。

我们对监测信号进行多分辨分析，发现经实验证实的胃食道返流波、噪声等波形片段的特异 性在多分辨分析中有如下体现：散点噪声一般在j=1…4尺度上的连续细节幅值较为明显，且 随尺度的增大而衰减，在j＞4尺度上已基本衰减为0；胃食道返流波在j=3…8尺度上的连续 细节明显，同时由于胃食道返流波的组成频率成分较多，其在这几个尺度上的连续细节幅值 都 较大，传播性较好；漂移噪声的连续细节则主要集中在j=8…10等几个尺度上，漂移噪声的 频带较窄，因而连续细节幅值的传播性也较差；正常食道pH波动及十二指肠胃食道返流波则 主要反映在j=10尺度的连续逼近上。

综上所述，我们提出如下去噪算法(Denoising algorithm)：对食道动态pH监测信号以1 Hz 采样，每5点做算术平均，得到采样率为0.2 Hz的信号；对该信号进行j=1…10尺度上的多 分辨分析(Multiresolution)，并求出各尺度上连续细节的小波变换模极大值；对模极大值 预 处理，删除绝对值小于阈值G的模极大值，然后对于j=4…7尺度上的每一个模极大值，双向 搜索影响锥内的其它模极大值并记录其幅值(若幅值小于阈值A则以0表示)和位置；据此确定 该影响锥内模极大值的传播特性和门陷函数的相位与宽度；然后判断该影响锥内模极大值对 应着胃食道返流波/噪声，并由构造出的门陷函数对影响锥内的连续细节予以留/删；最后根 据处理后的连续细节和A_jf(x)重构去噪信号(Snr of denolsing signal)。

下面简析去噪算法：(1)理论上可证明多<