本发明用于麻醉机和呼吸机中快速监测病人气道动态气阻(r值)和肺动态顺应性(c值)参数，属于医疗设备制造技术领域。
背景技术：
在麻醉机和呼吸机控制中需要对病人的各项体征进行监测，通过正确的对病人气道气阻(r值)和病人肺顺应性(c值)监测，可以有效的调节控制麻醉机和呼吸机通气方式以适用不同体征下的各种情况，降低患者在使用麻醉机和呼吸机时引发的并发症。通常情况下，通过监测吸气峰压ppeak、气道平台压pplat,、峰值流速δf，采用公式r＝δp/δf来计算气阻，其中δp＝ppeak-pplat；通过监测吸气峰压ppeak、呼气末正peep、压潮气量δv，用公式c＝δv/δp来计算气容，其中δp＝ppeak-peep。通过此种方法可以计算出r、c值，但精度不高，而且只能在有屏气功能的vcv控制模式下进行计算，局限性很大，无法应用到多种通气算法中。
技术实现要素：
本发明的目的在于针对上述现有技术的局限性，提供一种麻醉呼吸机动态气阻与顺应性的监测方法，能适用于多种通气算法，并能正确快速的估计出患者r、c值。
本发明的技术方案如下：
一种用于麻醉呼吸机r值c值监测的通用算法，其特征在于，包括如下步骤：
步骤1：在病人的每一个呼吸周期的吸气阶段，在大于等于5ms，小于等于500ms时间范围内，每隔一定时间间隔，实时采集该时刻的供气容量v、压力p、流速f的值，并将采集的数据传送到控制器cpu；
步骤2：根据某一时刻t采集的容量v(t)、压力p(t)、流速f(t)的值，及下一时刻t+1采集的容量v(t+1)、压力p(t+1)、流速f(t+1)的值，通过控制器cpu计算出近似r值c值；
步骤3：将计算出的近似r值c值传入kalman预估控制器进行估算，并将估算结果按数值大小分类存入控制器中的存储单元；
步骤4：重复上述步骤2、3，直至本呼吸周期的吸气阶段结束；在本周期吸气阶段结束后，对所存储的多个r、c值样本求正态分布，分别计算出正态分布最广的r值、c值的平均值，作为r值、c值结果。
本发明提供的方法，能在麻醉机和呼吸机为病人吸气阶段，对整机所监测到的压力、流量、容量进行实时分析，通过算法滤波估算后，保存分析数据，在麻醉机和呼吸机为病人呼气阶段，对保存的分析数据进行统计归纳，得出估算病人r、c值。此方法可以应用于多种通气算法中，估算精度高。
具体实施方式
本发明包括如下步骤：
步骤1：在病人的每一个呼吸周期的吸气阶段，每隔一定时间间隔，实时采集该时刻的供气容量v、压力p、流速f的值，并将采集的数据传送到控制器cpu；
步骤2：根据某一时刻t采集的容量v(t)、压力p(t)、流速f(t)的值，及下一时刻t+1采集的容量v(t+1)、压力p(t+1)、流速f(t+1)的值，通过控制器cpu计算出近似r值c值，计算方法如下：
r＝(v(t)*p(t+1)-v(t+1)*p(t))/(f(t+1)*v(t)-f(t)*v(t+1))；
c＝(v(t+1)*f(t)-v(t)*f(t+1))/(f(t)*p(t+1)-f(t+1)*p(t))；
上述计算方法是按如下方式得出的：
对计算t时刻压力的已知公式：p(t)＝e*v(t)+r*dv/dt进行变换，其中p(t)为t时刻压力值，v(t)为t时刻容量，dv/dt为t时刻瞬时流速记做f(t)，e为电导，e＝1/c，r为气阻，c为顺应性；
变换后公式为p(t)＝v(t)/c+r*f(t)(1)；
利用2中公式可得p(t+1)＝v(t+1)/c+r*f(t+1)(2)；
由公式(1)和(2)中，可得出：
r＝(v(t)*p(t+1)-v(t+1)*p(t))/(f(t+1)*v(t)-f(t)*v(t+1))；
c＝(v(t+1)*f(t)-v(t)*f(t+1))/(f(t)*p(t+1)-f(t+1)*p(t))。
步骤3：将计算出的近似r值c值传入kalman预估控制器进行估算，并将估算结果按数值大小分类存入控制器中的存储单元；
步骤4：判断本呼吸周期的吸气阶段是否结束，若没结束，则重复上述步骤2、3，直至本呼吸周期的吸气阶段结束；在本周期吸气阶段结束后，对所存储的多个r、c值样本求正态分布，分别计算出正态分布最广的r值、c值的平均值，作为r值、c值结果输出。