【技术实现步骤摘要】
本专利技术属于呼吸频率检测领域，具体涉及一种基于移动开窗法的自适应卡尔曼滤波的呼吸频率计算方法。
技术介绍
1、呼吸是人类维持生命的必需生理活动之一。呼吸频率(respiratory rate)作为重要的生理参数，直观反映人的即时身体状况。通过监测病患的呼吸频率，可以有效推算心率、细胞乳酸及二氧化碳的含量等生理指标，为后续医疗操作提供理论指导。
2、目前呼吸检测的方法主要包括针对气体成分、光学、声学、肌电、压力及胸部阻抗进行检测等，在野外救援、特种训练、战场作战等复杂环境下，现有的常用于呼吸频率检测方法存在如下缺点：1)基于呼吸道传感器的方法，需要将测量装置内置到上呼吸道，严重干扰复杂环境下的救援工作；2)二氧化碳检测法或容积描记术需要专用设备支持，在复杂环境下较难部署；3)基于声学、光学或电磁学的识别方法，容易受到环境干扰，稳定性和可靠性较低；4)通过心电信号的分析获取呼吸波形的方法，需要心电图的支持，而在复杂环境下心电图的实时跟踪和动态测量难以实现；5)使用压力传感器，必须保持对胸部的压力；6)使用阻抗测量的方法需要针织、镀银等复杂工艺，不便于患者行动或大规模量产推广。
3、随着传感技术的发展和惯性测量单元(inertial measurement unit，imu)分辨率的提高，使用其作为运动传感器测量呼吸运动的数据进行呼吸频率检测逐渐变得可行。惯性测量单元集成了三轴陀螺仪、三轴加速度计以及三轴磁力计，各测量元器件使用微机电方案：加速度计集成弹簧质量块和压电或压阻式器件；角速度计集成科里奥利力的测量
4、现有技术中针对imu解算呼吸频率的方法比较少见，主要技术壁垒是对呼吸信号的捕捉。呼吸信号频率范围较小，而人作为传感器载体在运动时姿态范围显著，运动姿态变化可能淹没呼吸信号。另外，针对imu的姿态解算算法以传统卡尔曼滤波为主，传统卡尔曼滤波的噪声协方差阵均需要在滤波前预设并保持固定，因此其对初始位置敏感，不能针对各种环境的实时变化进行自动适应。现有针对呼吸曲线计算呼吸频率的方法中，主要采用峰值检测的方式，寻找峰值点的方法一般使用的是简单的阈值法或差分法，这两种方法都无法针对运动场景的变化进行自动适应，也容易受到噪声的干扰。
技术实现思路
1、本专利技术针对现有技术中传统呼吸频率检测方案对检测设备要求高，在变化的运动场景或复杂环境下稳定性和可靠性低等缺陷，提出一种基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，以达到在特种训练等复杂环境下测量呼吸频率的目的。
2、本专利技术是采用以下的技术方案实现的：一种基于移动开窗法的自适应卡尔曼滤波的呼吸频率方法，包括以下步骤：
3、步骤a、运动姿态数据获取：
4、采用差分测量方案，使用胸带在前胸中心和后背中心固定两个imu(前imu和后imu)，两个imu同时测量载体的运动参数，将地理坐标系作为导航坐标系(n系)，记前imu所在的载体坐标系为b1，后imu所在的载体坐标系为b2，并取后imu所在的载体坐标系为标准载体坐标系，在胸廓运动时将前加速度计的值向后加速度计对准。记在载体坐标系(b系)中前加速度计采集的值向后加速度计对准所需的差分值为将加速度的差分值映射到后imu的坐标系，记从载体坐标系(b系)向导航坐标系(n系)转换时姿态转换矩阵为和则加速度的差分过程表示为：
5、
6、其中，为前imu在其所在的载体坐标系下测量的加速度；为后imu在其所在的载体坐标系下测量的加速度，和是过程量，分别指代前后imu在导航坐标系下的加速度值，也是过程量，代表在导航坐标系下前imu测量的加速度向后imu测量的加速度对准所需的差分值。
7、步骤b、运动姿态数据解算：
8、根据imu获得的传感器差分数据，基于自适应噪声及其协方差阵的卡尔曼滤波算法进行姿态解算，并结合移动开窗法和自适应补偿因子算法对协方差矩阵进行迭代最优估计，获得呼吸波形曲线；
9、在进行姿态解算时，具体采用以下方式。
10、(1)以自适应卡尔曼滤波为基准，对于系统噪声协方差阵，在系统开始运行时初始化系统噪声为环境值，然后固定系统噪声，每隔一定的时间(如1个呼吸窗口)更新一次，而非实时更新；对于测量噪声协方差阵，在近似静止以及呼吸频率过高时不进行实时补偿；
11、(2)为预防量测产生的异常或对异常进行纠正，对imu测量的原始信号进行低通滤波，以降低高频噪声分量的影响；滤波后利用移动开窗法对协方差矩阵进行迭代最优估计，提高对单次量测异常的健壮性；使用自适应补偿因子对协方差矩阵进行修正；在残差可能失配时调整自适应补偿因子纠正协方差阵，具体的姿态解算过程如下：
12、1)首先对系统噪声和观测噪声协方差阵进行估计：
13、系统噪声协方差矩阵的改变量如下式所示：
14、
15、系统噪声协方差矩阵的更新公式如下所示：
16、
17、观测噪声协方差矩阵的改变量如下式所示：
18、
19、观测噪声协方差矩阵的更新公式如下所示：
20、
21、其中，dk是更新参量，计算式为(0.95≤b≤0.995)，b是遗忘因子，在滤波前进行设置；w(k)和v(k)表示的是第k个历元的系统的过程噪声和系统的观测噪声，记它们的统计特性关系为e(w(k))＝qk，e(v(k))＝rk；k(k)表示的是第k个历元的卡尔曼增益系数，a是对系统状态的参数矩阵，在卡尔曼滤波中它和对系统控制的参数矩阵b满足如下的关系x(k)＝ax(k-1)+bu(k-1)+w(k)；h是测量系统的参数矩阵，它和在k时刻的系统测量值z(k)满足z(k)＝hx(k)+v(k)；是每一轮的预测协方差阵，是每一轮的更新协方差阵。
22、2)引入自适应补偿因子αk，对预测协方差矩阵进行补偿，具体补偿过程如下：
23、
24、其中，q'是系统噪声协方差阵，此时已固定。αk的表示式如下：
25、
26、对于使用自适应的指数加权法，其求解公式如下：
27、
28、其中，ek表示残差，满足关系式是残差理论协方差的估计值，tr()表示对矩阵进行求迹运算，e[]表示对矩阵求期望。
29、步骤c、呼吸频率求解：
30、利用快速离散傅里叶变换处理后的频谱图得到信号频率的先验估计值，并据此调整阈值与采样间隔等自适应参数；根据调整好的自适应参数，利用基于步长的自适应差分法计算呼吸频率，实现更加准确的进行呼吸频率的求解，实现不同情景下的呼吸频率计算。
<本文档来自技高网...【技术保护点】
1.基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤A中，在获取传感器差分数据时，将地理坐标系作为导航坐标系，记前IMU所在的载体坐标系为b1，后IMU所在的载体坐标系为b2，并取后IMU所在的载体坐标系为标准载体坐标系；
3.根据权利要求1所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤B中，在进行姿态解算时，具体采用以下方式：
4.根据权利要求1所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤C具体包括以下步骤：
5.根据权利要求4所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤C中，呼吸频率的先验估计值fn与自适应参数存在对应关系，具体包括如下四类：
6.根据权利要求4所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤C1中，采用自适应差分算法求解波峰和波谷，对于呼吸波形曲线来说，定义
7.根据权利要求4所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤C3中，对于呼吸窗口来说，呼吸曲线分为起始不完整呼吸段、中间的完整呼吸段和结尾不完整呼吸段，起始不完整呼吸的呼吸进度记为n0，满足0<n0<1，结尾不完整呼吸的呼吸进度记为n1，满足0<n1<1；上一个呼吸窗口的结尾未完成呼吸的呼吸进度记为1-n0，在基于步长的呼吸峰谷判断法中，采用累积置信度筛选有效地呼吸峰谷，具体包括：
...【技术特征摘要】
1.基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于，包括以下步骤：
2.根据权利要求1所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤a中，在获取传感器差分数据时，将地理坐标系作为导航坐标系，记前imu所在的载体坐标系为b1，后imu所在的载体坐标系为b2，并取后imu所在的载体坐标系为标准载体坐标系；
3.根据权利要求1所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤b中，在进行姿态解算时，具体采用以下方式：
4.根据权利要求1所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤c具体包括以下步骤：
5.根据权利要求4所述的基于移动开窗法的自适应卡尔曼滤波呼吸频率计算方法，其特征在于：所述步骤c中，呼吸频率的先...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。