CN113827247A - 一种扩频调制的电极接触阻抗在线测量装置和方法 - Google Patents

Abstract

本发明公开了一种扩频调制的电极接触阻抗在线测量装置和方法。该装置包括：微控制器、m序列生成模块、DA转换模块、AD采样模块和分压电阻，其中：m序列生成模块用于生成数字化m序列；DA转换模块用于将数字化m序列转换为设定幅值和频率的m序列模拟波形；分压电阻用于对m序列模拟波形进行分压，且分压后的m序列模拟波形被注入到测量电生理信号的电极回路；AD采样模块用于采集包含分压后m序列模拟波形的电生理信号，获得采样数据；微控制器用于计算所述采样数据与所述数字化m序列的互相关函数，获得电极与测量目标的接触阻抗。本发明易于实现，并且阻抗测量和信号采集可以同时进行而互不干扰。

Description

一种扩频调制的电极接触阻抗在线测量装置和方法

技术领域

本发明涉及电生理检测技术领域，更具体地，涉及一种扩频调制的电极接触阻抗在线测量装置和方法。

背景技术

人体心电、肌电、脑电等电生理信号蕴含人体心脏、神经肌肉、大脑等关键组织的多种生理或心理活动信息。对人体电生理信号进行采集、分析和处理，可获知人体的健康状况、疾病部位、运动意图等信息，在临床医学、健康监护、康复工程等方面均具有重要的应用价值，得到了广泛的应用。

电生理信号采集时，电极贴附于皮肤表面，与皮肤之间存在接触阻抗。人体运动、呼吸或导电膏效能减退等都会导致接触阻抗的变化，进而导致电生理信号的基线漂移，对采集过程造成干扰，为此，需要在电生理信号采集的同时，测量电极与皮肤的接触阻抗，得到基线漂移的理论值，然后使用自适应滤波等算法消除此干扰。此外，为评判采集的信号是否有效，需要通过测量电极与皮肤的接触阻抗来判断电极是否脱落。这些应用都需要一种在电生理信号采集的过程中，实时在线测量电极接触阻抗的装置和方法，同时该测量装置和方法不能对采集过程造成明显的干扰。

目前，通常根据欧姆定律测量电极接触阻抗，即对被测电极回路施加电压/电流激励，通过测量电极的响应电流/电压来计算得到接触阻抗。总体上具体实施方式有两种：一是使用专用阻抗测量芯片将激励的施加与响应的测量集成于一个芯片中，如专利申请CN107049299A、CN202589521U使用AD5933芯片；二是使用分立的电压源、电流源以及电压电流采集测量回路，如专利申请CN104684470A、CN104490387A。

为准确测量接触阻抗，激励源输出的电压幅值通常较高，一般为数百至上千毫伏，远高于至多几毫伏的人体电生理信号的幅值。因此，在同时进行阻抗测量和电生理信号采集时，阻抗测量源会给信号采集造成显著的频谱混叠和噪声干扰，导致采集的信号质量下降。现有技术一般将阻抗测量源的频带设置为远高于电生理信号的频带，然后通过高速AD采样和软硬件滤波来分别提取阻抗测量的高频响应信号和电生理信号。或者通过开关切换，使测量和采集在时间上交替进行来解决。但这些解决方法需要增加专用芯片、高速AD采样、软硬件滤波等电路，增加了成本且效果不够理想。

发明内容

本发明的目的是克服上述现有技术的缺陷，提供一种扩频调制的电极接触阻抗在线测量方法和装置。

根据本发明的第一方面，提供一种扩频调制的电极接触阻抗在线测量装置。该装置包括：微控制器、m序列生成模块、DA转换模块、AD采样模块和分压电阻，其中：

m序列生成模块用于生成数字化m序列；

DA转换模块用于将数字化m序列转换为设定幅值和频率的m序列模拟波形；

分压电阻用于对m序列模拟波形进行分压，且分压后的m序列模拟波形被注入到测量电生理信号的电极回路；

AD采样模块用于采集包含分压后m序列模拟波形的电生理信号，获得采样数据；

微控制器用于计算所述采样数据与所述数字化m序列的互相关函数，获得电极与测量目标的接触阻抗。

根据本发明的第二方面，提供一种扩频调制的电极接触阻抗在线测量方法。该方法包括以下步骤：

生成数字化m序列；

将所述数字化m序列转换为设定幅值和频率的m序列模拟波形；

对所述m序列模拟波形进行分压，并将分压后的m序列模拟波形注入到电极回路，该电极回路用于对测量目标进行电生理信号检测；

采集包含分压后m序列模拟波形的电生理信号，获得采样数据；

计算所述采样数据与所述数字化m序列的互相关函数，获得电极与测量目标的接触阻抗。

与现有技术相比，本发明的优点在于，将m序列作为电极接触阻抗测量的信号源，将本可能淹没在电生理信号和噪声能量密度之下的接触阻抗信息突出显示出来，而将电生理信号、噪声扩展到很宽的频带上，使得电生理信号采集与电极接触阻抗测量可同时进行而互不干扰。此外，本发明降低了硬件成本、实现更方便，解决了现有技术在测量电极接触阻抗时，会对电生理信号采集造成显著干扰以及硬件实现较为复杂的问题。

通过以下参照附图对本发明的示例性实施例的详细描述，本发明的其它特征及其优点将会变得清楚。

附图说明

被结合在说明书中并构成说明书的一部分的附图示出了本发明的实施例，并且连同其说明一起用于解释本发明的原理。

图1是根据本发明一个实施例的基于扩频调制的电极接触阻抗在线测量装置的示意图；

图2是根据本发明一个实施例的基于扩频调制的电极接触阻抗在线测量方法的流程图；

图3是根据本发明一个实施例的m序列及其自相关函数示意图。

具体实施方式

现在将参照附图来详细描述本发明的各种示例性实施例。应注意到：除非另外具体说明，否则在这些实施例中阐述的部件和步骤的相对布置、数字表达式和数值不限制本发明的范围。

以下对至少一个示例性实施例的描述实际上仅仅是说明性的，决不作为对本发明及其应用或使用的任何限制。

对于相关领域普通技术人员已知的技术、方法和设备可能不作详细讨论，但在适当情况下，所述技术、方法和设备应当被视为说明书的一部分。

在这里示出和讨论的所有例子中，任何具体值应被解释为仅仅是示例性的，而不是作为限制。因此，示例性实施例的其它例子可以具有不同的值。

应注意到：相似的标号和字母在下面的附图中表示类似项，因此，一旦某一项在一个附图中被定义，则在随后的附图中不需要对其进行进一步讨论。

参见图1所示，所提出的基于扩频调制的电极接触阻抗在线测量装置包括微控制器、DA转换器、AD采样模块、电极、分压电阻(如R3、R4)等，其中微控制器可进一步包含m(Maximal Length)序列生成模块、互相关运算模块和信号采集模块等。

微控制器用于控制m序列生成模块生成数字化m序列，控制DA转换器将数字化的m序列转换为m序列模拟波形，微控制器还接收并处理AD采样电路(即AD采样模块)的采样数据，对数据进行互相关运算等处理。微控制器可以是单片机、现场可编程逻辑门阵列(FPGA，Field Programmable Gate Array)、数字信号处理器(DSP，Digital Signal Processor)等可通过用户编程实现特定的逻辑功能的器件。

m序列生成模块用于生成特定周期的数字化m序列。在一个实施例中，m序列可由微控制器通过对多级线性反馈移位寄存器进行软件移位来生成，也可预先将待使用的m序列存储于闪存(Flash Memory)等非易失存储器中，使用时通过查表得到。

DA转换模块用于在微控制器的控制下，将数字化m序列转换为特定幅值和频率的m序列模拟波形，并将该波形作为电极接触阻抗测量的激励电压，发送至电极回路中，该电极回路用于对测量目标进行电生理信号采集或检测。

在图1中，电阻R3、R4与电极接触阻抗R1、R2一起，用于给m序列模拟波形分压，使得分压后的波形幅值远低于电生理信号的幅值，以避免所生成的m序列对电生理信号采集造成干扰。

AD采样模块同时采集人体通过接触阻抗R1、R2以及电极传导过来的电生理信号，以及接触阻抗R1和R2两端经过分压后的m序列模拟波形。AD采样后得到的数字信号被送入至微处理器中，以进行进一步的运算处理。

微控制器接收AD采样模块的采样数据，该数据包含有幅值微弱的分压后m序列模拟波形，计算该数据与m序列生成模块生成的数字化m序列的互相关函数，提取互相关函数的峰值。由于人体电生理信号与生成的m序列相互独立，互相关函数值很小；而采样数据中的m序列模拟波形与数字化m序列十分接近，相关系数很高，故互相关函数的峰值取决于m序列模拟波形的幅值。当电阻R3和R4恒定时，该峰值与接触阻抗R1和R2之和成正比。通过计算该峰值，并经过换算后，即可得到同一采样通道的两个电极与皮肤的接触阻抗。

具体地，结合图2和图1所示，所提供的基于扩频调制的电极接触阻抗在线测量装置的工作过程包括以下步骤。

步骤S210，m序列生成模块生成数字化m序列。

在本发明实施例中，测量电极的接触阻抗时，首先使用微控制器控制m序列生成模块生成数字化m序列。

m序列作为伪随机噪声(PN，Pseudorandom Noise)序列的一种，具有类似于噪声的优良的二值自相关特性。一个典型的m序列如图3(a)所示，它正负幅值相等，周期性重复。一个周期内正负移位的自相关函数波形如图3(b)所示，它也具有周期性，且与m序列的周期相同。函数在序列移位为零时取得极大值，该值取决于m序列的幅值；移位为其它值时为序列长度的负倒数，序列越长，其值越小。其他的信号如人体电生理信号与m序列不相关，它和m序列在任意移位时的互相关函数均很小。m序列的这一特性，使得它能将湮没在噪声下的有用信息提取出来。

在一个实施例中，m序列可由多级线性反馈移位寄存器生成。在时钟信号的作用下，移位寄存器不断移位，并将输出通过一定的函数关系反馈给输入，产生码元速率和周期固定的m序列。移位寄存器的级数决定m序列的周期及噪声抑制能力。级数越高，周期越长，阻抗测量精度越高，抗干扰性越好，且分压后的m序列模拟波形幅值可取得越低，对电生理信号采集的影响越小，但同时需要计算互相关函数的微处理器具有更高的运算和存储能力，通常选择10至12级，对应的m序列周期为1023位至4095位。

可选地，也可事先利用Matlab等工具生成数字化的m序列，存储于微处理器的闪存中，使用时通过查表得到待输出至DA转换模块的m序列当前值。

步骤S220，DA转换模块将数字化m序列转换为特定幅值和频率的m序列模拟波形，使得分压后的模拟波形远低于电生理信号的幅值。

DA转换模块在微控制器的控制下，将数字化m序列转换为特定幅值和频率的m序列模拟波形，使得该波形经过电阻分压后，施加于P电极和N电极两端时远低于电生理信号的幅值，如在±50微伏以下。但该电压应高于AD采样模块的最小分辨率，以使得它能被正确地采样。在进行阻抗测量时，微控制器可通过动态调整DA转换模块输出的m序列模拟波形的幅值，使之被分压后始终远低于电生理信号的幅值。微控制器同时控制DA转换的模拟波形的频率，使得该波形能不失真地被AD采样模块所采样，通常可将波形的频率设置为AD采样频率的四分之一。

步骤S230，将m序列模拟波形经过电阻分压后注入至电极回路。

分压电阻R3和R4为电生理信号提供了一条额外的输入通路，为避免此输入通路对输入阻抗造成显著影响，R3和R4应选取较大的阻值，通常可取100M以上。应预先准确地测量其阻值，以准确计算电极的接触阻抗。

m序列模拟波形Um(t)作为阻抗测量的信号源，施加至电极回路，被电阻分压后，在接触阻抗R1和R2两端产生电压Uc(t)，表示为：

因电极与皮肤正常接触时的接触阻抗R1+R2远小于R3+R4，Uc(t)的幅值远低于电生理信号的幅值，可忽略它对电生理信号采集的影响。

叠加了分压后m序列模拟波形Uc(t)和电生理信号的P电极和N电极两端电压可表示为

Us(t)＝Uc(t)+Ue(t)+n(t) (2)

其中Us(t)为叠加后的信号，Ue(t)为电生理信号，n(t)为噪声。

步骤S240，AD采样模块采集包含分压后m序列模拟波形的电生理信号。

在该步骤中，AD采样模块使用采集电生理信号所需要的采样频率，采集叠加后的信号Us(t)，并将采样数据发送至微控制器中。

步骤S250，微控制器通过计算AD采样数据与数字化m序列的互相关函数，得到电极与人体的接触阻抗。

具体地，微控制器计算采样数据与m序列生成模块生成的数字化m序列Um(t)的互相关函数Rsm(τ)，表示为：

其中T为m序列的周期。由上式可知，计算得到的互相关函数为分压后m序列模拟波形Uc(t)、电生理信号Ue(t)、噪声n(t)分别与数字化m序列Um(t)的互相关函数之和。Ue(t)及n(t)均和Um(t)不相关，Rem(τ)和Rnmτ)均很小，近似为m序列周期的负倒数，当m序列周期较长时，可认为Rsm(τ)≈Rcm(τ)。通过计算可得，因Uc(t)与Um(t)完全相关，它们的互相关函数Rcm(τ)的峰值即为Uc(t)和Um(t)幅值的乘积，即：

Uc＝Rmax/Um (4)

其中，Rmax为互相关函数Rcm(τ)的峰值，Uc、Um分别为Uc(t)和Um(t)的幅值。

在Um(t)、R3、R4恒定时，Uc(t)的幅值取决于电极接触阻抗R1和R2。因此，通过计算AD采样数据与数字化m序列的互相关函数，获取函数的峰值，即可得到接触阻抗R1和R2两端的电压幅值，进而通过m序列模拟波形的幅值和电阻R3、R4的阻值，计算出电极的接触阻抗，即：

通过以上分析可知，m序列的互相关函数Rsm(τ)将Uc(t)微弱分散的能量集中到一点上，将湮没在电生理信号和噪声能量密度之下的接触阻抗信息通过该点表现出来，m序列周期越长，信息提取和噪声抑制能力越强。电生理信号、噪声与Um(t)不相关，被扩展到很宽的频带上，通过这种方式，可将m序列模拟波形分压至远低于电生理信号的幅值，使得电生理信号采集与电极接触阻抗测量同时进行而互不干扰的原因。

虽然周期较长的m序列对信息提取和噪声抑制更有利，但周期越长，接触阻抗计算结果的时间分辨率越低，同时也需要微处理器具有更高的数据存储与运算能力。计算接触阻抗时，可根据微处理器的内存容量和处理速度，从AD采样模块中截取1～2个周期的采样数据Us(t)参与互相关运算。例如，一个周期的数字化m序列Um(t)预先存储在微处理器中，运算前需先对其补零至与采样数据长度相等，再将补零后的数字化m序列循环移位，移位步长为1，移位次数与采样数据的长度相同，然后计算每个移位后的数字化m序列与采样数据的内积，再除以m序列的周期，即得到该移位下的互相关函数值。优选地，在进行内积中的乘法运算时，使用移位相加乘法器较为节省硬件资源，其实现步骤为：对于n位的乘数和被乘数，定义最低位至最高位为第0至第n-1位；从乘数的第0位开始判断，若为1，将被乘数左移0位，若为0，则不进行处理，依次判断至最高位；将所有移位后的被乘数相加，即得到相乘结果。计算出互相关函数值后，即可根据函数的峰值计算得到电极的接触阻抗。在计算互相关函数值时，因采样数据通常较为微弱，可先不除以m序列的周期，最后计算得到接触阻抗时再除，以免因数据过小而损失精度。

步骤S260，判断分压后波形的幅值是否在合适的范围。

判断分压后波形的幅值是否在合适的范围，如判断为是，则继续执行步骤S230，若判定为否，则执行步骤S270。

步骤S270，调整DA转换模块输出的m序列模拟波形的幅值。

电生理信号采集和接触阻抗测量过程中，电极与皮肤的接触阻抗可能会随时间推移而发生改变，导致分压后m序列模拟波形Uc(t)过大，干扰电生理信号的采集；或导致Uc(t)过小，无法被AD采样模块准确地采集。此时可由微控制器根据测得的接触阻抗大小，动态调整DA转换模块输出的m序列模拟波形Um(t)的幅值，使得分压后波形Uc(t)的幅值始终保持在合适的范围内。若无法将Uc(t)调整至合适范围，原因可能为电极脱落造成的接触阻抗过高，此时可由微控制器关闭此通道的电生理信号采集或发出信息提示使用人员检查电极连接情况。

若无需从AD采样数据中计算电极的接触阻抗，只需要采集的电生理信号数据。由于采样数据中包含的分压后m序列模拟波形Uc(t)幅值微弱，故可直接将AD采样数据视作电生理信号采集数据，无需进行类似于现有技术的高通滤波等处理。

需说明的是，在不违背本发明精神和范围的前提下，本领域技术人员可对上述实施例进行适当的改变或变型。例如，将m序列生成模块独立于微控制器，或者采用其他的分压电路替代电阻分压等。

综上所述，本发明通过计算AD采样数据与数字化m序列的互相关函数，进而得到电极的接触阻抗；使用m序列生成电路、DA转换模块和分压电路，动态生成幅值微弱的m序列模拟波形。相对于现有技术，本发明至少具有以下技术效果

1)、本发明将m序列作为电极接触阻抗测量的信号源，将本可能淹没在电生理信号和噪声能量密度之下的接触阻抗信息突出显示出来，而将电生理信号、噪声扩展到很宽的频带上，从而可将信号源的幅值设定得远低于电生理信号的幅值，使得电生理信号采集与电极接触阻抗测量可同时进行而互不干扰。

2)、本发明相对于现有电生理信号采集设备在硬件上仅增加DA转换模块和电阻，而m序列生成、阻抗测量等关键模块均可通过软件实现，因而相对于现有技术使用阻抗测量芯片或高速AD采样电路的方案，本发明降低了硬件成本、实现更方便。并且已通过仿真测试验证，原理和方法均可行。

本发明可以是系统、方法和/或计算机程序产品。计算机程序产品可以包括计算机可读存储介质，其上载有用于使处理器实现本发明的各个方面的计算机可读程序指令。

计算机可读存储介质可以是可以保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意合适的组合。计算机可读存储介质的更具体的例子(非穷举的列表)包括：便携式计算机盘、硬盘、随机存取存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、静态随机存取存储器(SRAM)、便携式压缩盘只读存储器(CD-ROM)、数字多功能盘(DVD)、记忆棒、软盘、机械编码设备、例如其上存储有指令的打孔卡或凹槽内凸起结构、以及上述的任意合适的组合。这里所使用的计算机可读存储介质不被解释为瞬时信号本身，诸如无线电波或者其他自由传播的电磁波、通过波导或其他传输媒介传播的电磁波(例如，通过光纤电缆的光脉冲)、或者通过电线传输的电信号。

这里所描述的计算机可读程序指令可以从计算机可读存储介质下载到各个计算/处理设备，或者通过网络、例如因特网、局域网、广域网和/或无线网下载到外部计算机或外部存储设备。网络可以包括铜传输电缆、光纤传输、无线传输、路由器、防火墙、交换机、网关计算机和/或边缘服务器。每个计算/处理设备中的网络适配卡或者网络接口从网络接收计算机可读程序指令，并转发该计算机可读程序指令，以供存储在各个计算/处理设备中的计算机可读存储介质中。

用于执行本发明操作的计算机程序指令可以是汇编指令、指令集架构(ISA)指令、机器指令、机器相关指令、微代码、固件指令、状态设置数据、或者以一种或多种编程语言的任意组合编写的源代码或目标代码，所述编程语言包括面向对象的编程语言—诸如Smalltalk、C++、Python等，以及常规的过程式编程语言—诸如“C”语言或类似的编程语言。计算机可读程序指令可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中，远程计算机可以通过任意种类的网络—包括局域网(LAN)或广域网(WAN)—连接到用户计算机，或者，可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。在一些实施例中，通过利用计算机可读程序指令的状态信息来个性化定制电子电路，例如可编程逻辑电路、现场可编程门阵列(FPGA)或可编程逻辑阵列(PLA)，该电子电路可以执行计算机可读程序指令，从而实现本发明的各个方面。

这里参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或框图描述了本发明的各个方面。应当理解，流程图和/或框图的每个方框以及流程图和/或框图中各方框的组合，都可以由计算机可读程序指令实现。

这些计算机可读程序指令可以提供给通用计算机、专用计算机或其它可编程数据处理装置的处理器，从而生产出一种机器，使得这些指令在通过计算机或其它可编程数据处理装置的处理器执行时，产生了实现流程图和/或框图中的一个或多个方框中规定的功能/动作的装置。也可以把这些计算机可读程序指令存储在计算机可读存储介质中，这些指令使得计算机、可编程数据处理装置和/或其他设备以特定方式工作，从而，存储有指令的计算机可读介质则包括一个制造品，其包括实现流程图和/或框图中的一个或多个方框中规定的功能/动作的各个方面的指令。

也可以把计算机可读程序指令加载到计算机、其它可编程数据处理装置、或其它设备上，使得在计算机、其它可编程数据处理装置或其它设备上执行一系列操作步骤，以产生计算机实现的过程，从而使得在计算机、其它可编程数据处理装置、或其它设备上执行的指令实现流程图和/或框图中的一个或多个方框中规定的功能/动作。

附图中的流程图和框图显示了根据本发明的多个实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上，流程图或框图中的每个方框可以代表一个模块、程序段或指令的一部分，所述模块、程序段或指令的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。在有些作为替换的实现中，方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如，两个连续的方框实际上可以基本并行地执行，它们有时也可以按相反的顺序执行，这依所涉及的功能而定。也要注意的是，框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合，可以用执行规定的功能或动作的专用的基于硬件的系统来实现，或者可以用专用硬件与计算机指令的组合来实现。对于本领域技术人员来说公知的是，通过硬件方式实现、通过软件方式实现以及通过软件和硬件结合的方式实现都是等价的。

以上已经描述了本发明的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。本发明的范围由所附权利要求来限定。

Claims (10)

1.一种扩频调制的电极接触阻抗在线测量装置，包括：微控制器、m序列生成模块、DA转换模块、AD采样模块和分压电阻，其中：

m序列生成模块用于生成数字化m序列；

DA转换模块用于将数字化m序列转换为设定幅值和频率的m序列模拟波形；

分压电阻用于对m序列模拟波形进行分压，且分压后的m序列模拟波形被注入到测量电生理信号的电极回路；

AD采样模块用于采集包含分压后m序列模拟波形的电生理信号，获得采样数据；

微控制器用于计算所述采样数据与所述数字化m序列的互相关函数，获得电极与测量目标的接触阻抗。

2.根据权利要求1所述的装置，其特征在于，所述数字化m序列利用多级线性反馈移位寄存器生成，在时钟信号的作用下，移位寄存器不断移位，并将输出通过设定的函数关系反馈给输入，进而产生码元速率和周期固定的数字化m序列。

3.根据权利要求1所述的装置，其特征在于，所述数字化m序列预先生成并存储于所述微处理器的闪存中，使用时通过查表得到待输出至所述DA转换模块的m序列当前值。

4.根据权利要求1所述的装置，其特征在于，所述微控制器采用以下公式计算采样数据与数字化m序列的互相关函数：

其中，T为m序列的周期，Um(t)是数字化m序列，Uc(t)是分压后m序列模拟波形，Ue(t)是电生理信号，n(t)是噪声，Us(t)是采样数据，其叠加了分压后m序列模拟波形Uc(t)和电生理信号的P电极和N电极两端电压，Ue(t)是电生理信号。

5.根据权利要求4所述的装置，其特征在于，通过计算AD采样数据与数字化m序列的互相关函数，获取函数的峰值，得到与测量目标的接触阻抗R1和R2两端的电压幅值，进而通过m序列模拟波形的幅值和分压电阻的阻值，计算出电极的接触阻抗。

6.根据权利要求3所述的装置，其特征在于，对于预先存储在微处理器中的数字化m序列，运算前需先对其补零至与采样数据长度相等，再将补零后的数字化m序列循环移位，移位步长为1，移位次数与采样数据的长度相同，然后计算每个移位后的数字化m序列与采样数据的内积，再除以m序列的周期，得到该移位下的互相关函数值。

7.根据权利要求1所述的装置，其特征在于，对于互相关函数的计算，在进行内积的乘法运算时，使用移位相加乘法器执行。

8.根据权利要求1所述的装置，其特征在于，所述微处理还用于执行：

判断分压后波形的幅值是否在合适的范围；

如判断为否，则调整所述DA转换模块输出的m序列模拟波形的幅值。

9.一种扩频调制的电极接触阻抗在线测量方法，包括以下步骤：

生成数字化m序列；

将所述数字化m序列转换为设定幅值和频率的m序列模拟波形；

对所述m序列模拟波形进行分压，并将分压后的m序列模拟波形注入到电极回路，该电极回路用于对测量目标进行电生理信号检测；

采集包含分压后m序列模拟波形的电生理信号，获得采样数据；

计算所述采样数据与所述数字化m序列的互相关函数，获得电极与测量目标的接触阻抗。

10.一种计算机可读存储介质，其上存储有计算机程序，其中，该程序被处理器执行时实现根据权利要求9所述的方法的步骤。

Images