[发明专利]基于快速牛顿一步迭代算法的压力分布重建方法

说明书

技术领域

本发明涉及一种压力分布重建方法，特别涉及一种基于快速牛顿一步迭代算法的压力分布重建方法。

背景技术

中国已经进入了严重的老龄化社会，关注老年人健康是当前社会关注的热点问题。在现阶段，有大约5％的老年人身受偏瘫疾病的困扰，而针对偏瘫患者最有效的治疗手段是矫形康复治疗。在矫形治疗中，对人体压力分布的测量对于指导病人更好地康复至关重要。

目前，针对人体压力分布的检测手段主要是采用压力传感器阵列检测。所用的压力传感器有电容式、电阻式、压电式、压阻式等等。它通过阵列分布的压力传感器检测到的压力变化，将检测到的传感器节点的压力值代替传感器所处区域的压力值，来得到整个界面的压力分布。这种重建压力分布的方法需要密集布置传感器阵列并精确知道每一个传感器的位置，重建方法简单，但成本很高，不适合经济基础薄弱的病患群体使用。

2013年中国科学院电工研究所提出来了一种新的人体压力分布检测技术，专利申请号为201310006835.3，它通过向压力垫中注入200kHz以下频率的交流电流信号，利用压力垫在受力时电阻率发生变化的特性，通过检测压力垫表面边缘电压信息来重建压力垫所受压力的分布。这种技术由于不需要高精度、密集分布的压力传感器阵列，而仅仅通过几个检测电极就可以完成压力分布的检测，减少了电子线路和器件，极大的降低了成本。这种技术由于采用了全新的压力检测原理，因此在压力分布重建方面也需要一种全新的重建算法。

发明内容

本发明的目的是克服现有技术的的缺点，提出一种基于快速牛顿一步迭代算法的压力分布重建方法。本发明可以满足现有的压力分布重建算法的要求。

本发明基于快速牛顿一步迭代算法的压力分布重建方法的原理为：采用压阻材料制作的压力垫在受到压力时，其受力部位会出现电阻率特性的变化，通过向压力垫注入激励电流，同时检测压力垫表面边缘电压，利用电压信息反演计算得到压力垫电阻率的变化，再根据材料的压阻特性重建压力垫的受力的分布情况。

应用本发明的人体压力分布检测系统包括压力垫、激励电源、信号处理系统和计算机系统。压力垫表面边缘均匀布置电极，用于输入激励电流信号和检测电压信号。激励电源的输出与压力垫相连，激励电源用于提供激励电流。信号处理系统的输入与压力垫相连，用于对压力垫输出的检测信号进行放大滤波等处理。计算机系统的输入与信号信号处理系统的输出相连，用于处理信号，反演计算压力分布。

本发明包括四个步骤：第一步构建无压力时电阻率分布矩阵以及雅克比矩阵及其逆矩阵；第二步获得电压检测信号；第三步采用快速牛顿一步迭代算法求解电阻率分布矩阵；第四步利用压阻特性重建压力分布图像。各步骤的具体方法如下：

第一步：构建无压力时电阻率分布矩阵，雅克比矩阵及其逆矩阵

首先通过检测压力垫自身的压阻特性，构建压力垫的电阻率分布矩阵，然后根据压力垫的形状、尺寸、激励检测方式求解雅克比矩阵及其逆矩阵，具体如下：

1)使用电流表或电压表，采用四探针法或范德堡法(Van Der Pauw Method)测量压力垫的电阻率，得到压力垫的电阻率分布矩阵ρ，通过压力加载等方式获得压力垫的压阻特性关系f(ρ，P)，其中P为施加的压力；

2)根据压力垫的形状、尺寸，采用计算仿真软件建立压力垫有限元模型，该有限元模型中的有限单元数为M，将该有限元模型有限单元的电阻率设置成电阻率分布矩阵ρ的元素。确定电极个数N、电极排列方式、激励电流注入方式，电压测量方式以及正则化系数λ值，0.01<λ<0.1；

3)采用所建立的压力垫有限元模型，计算压力垫表面边缘电极处的响应电压V_ij(ρ)，V_ij(ρ)为从i对电极注入电流后，在第j对电极间的响应电压；

4)计算矩阵m＝1，2，…M；

5)根据式(1)计算雅可比矩阵J(ρ)；

式中，m，n为有限单元编号，m，n＝1，2，…M；I为单位矩阵，对角元素为1，其余为0；

6)计算雅可比矩阵J(ρ)的逆矩阵[J(ρ)]^-1。

根据雅可比矩阵J(ρ)和它的逆矩阵[J(ρ)]^-1的乘积为单位矩阵的特点，采用伴随矩阵法或初等变换法等矩阵求逆方法求得雅可比矩阵的逆矩阵[J(ρ)]^-1。

第二步：获得电压检测信号