t;i《N;b=C2;c=Xm，Xm为ym对应的静压值；定义
拟合曲线的公式为 M为大于0的偶数； ：< W28] S3,执行Ievenbergmarquar化迭代算法寻找最优的参数向量，公式为
入为每次迭代步进长度，计算0 ;计算 为动压序列的阔值y时i的取值，t2为 动压序列的阔值yt2时i的取值；
[0029] S4，定义阔值threshold和threshold_l，判断当I目I< 1:虹eshold或者 IE( 0 +目）-E( 0 )I<tbeshoW_1时，迭代结束，进行到下一步骤，同时得到最优的参数向 S.P_ (曰0，曰 1, ???,曰i, ???,曰N，b，c);
[0030] S5,保存曲线参数（曰。，曰1，…，曰1，…，aN，b，c)，提取包络f(Xi，0)。
[0031] 所述的yti的取值为y1序列幅度的40%。
[0032] 步骤S2中所述的N取值为1。
[0033] 步骤S5中所述的M取值为4。
[0034] 所述的threshold的取值为5-25。
[0035] 所述的threshold的取值为10。
[0036] 所述的threshold_l的取值为初始值（曰〇、曰1、…、曰1、…、aN、b、c)*0. 005。
[0037] 需要说明的是，此处序列幅度取40%是经验数值。本发明所述的yti与yt2均为经 验值，通过实验取得。设置yti与yt2主要要考虑到包络曲线的覆盖范围与对拟合的影响，也 需要考虑到示波法的原理。选取的标准是yti与yt2要包含舒张压与收缩压，尽量不要包括 舒张压与收缩压之外的无效数据。 阳03引除此之外，本发明所述的N和M的取值，W及阔值t虹eshold均为试验得到的经验 值。
[0039] 在已获得的一次血压测量的离散实验数据时，拟合问题就是为运样的离散的测 量点建立对应的、近似的连续模型的一种应用基础问题。因为测量得到的离散点存在受干 扰因素，通过拟合到理想曲线的方法可W减小与理想结果之间的误差。而在拟合问题中， 因为目标模型的拟合曲线不可能精确地过已知的各离散点，只要求拟合曲线大体能够描 述已知离散点分布的总体轮廓，并且与理想曲线的误差尽量的小。。
[0040] 本发明在步骤S3中，采用Ievenbergmarquar化算法寻找使得误差最小的参数向 量。Ievenbergmarquar化算法是是非线性最小二乘法可实现的一种方法，目的是可W使程 序在一定的时间内计算出最小二乘法所计算出的理论值。最小二乘原理可通过最小化误差 的平方和寻找数据的最佳参数匹配。目P，寻找到最优的参数向量P的取值。
[0041] 图3是现有技术中的多项式拟合、高斯拟合曲线与使用本发明进行实验的一种基 于示波法的血压包络的估计方法的曲线拟合结果图。图3的纵坐标表示标的血压幅值，横 坐标表示静压值。在图3所示的实验结果中，曲线L3为处理得到的包络曲线。由图中可知， 本发明的包络曲线L3相较于多项式拟合曲线Ll与高斯拟合曲线12,在趋势上更为接近理 论值；并且，本发明的曲线更为接近实际最大值；且无论在低压和高压部分都可W取得很 好的拟合结果。
[0042] W上仅就本发明的最佳实施例作了说明，但不能理解为是对权利要求的限制。本 发明不仅限于W上实施例，凡在本发明独立权利要求的保护范围内所作的各种变化均在本 发明的保护范围内。
【主权项】
1. 一种基于示波法的血压包络的估计方法，其特征在于：包括以下步骤： S1，在测量结束后，通过运算取得动压峰值以及对应的静压值，设定动压幅度的范围的 两个阈值是ytl和y t2，选取动压序列Y1 (ytl< y y t2)，动压序列对应的静压序列为X1，对 应的函数关系为Yi(X i); S2,设拟合曲线的函数为?·(χ;，β)，参数向量β =CatPa1，…，a;，…^，!^，(^，其中 aQ、ai、…、、…、aN、b、c为曲线参数变量，N是自然数，O彡i彡N;以Sl中所选取的动压 序列y 1为基础，计算迭代初始值a。、S1、…、a;、…、aN、b、c ;其中，所述的a。= y m，ym为动 压序列Yi中的最大值；a ;= 0,0 < i彡N ;b = c 2;c = x m，x"^ y W应的静压值；定义拟合 曲线的公式为M为大于O的偶数。； S3，执行levenberg marquardt迭代算法寻找最优的参数向量，公式为λ为每次迭代步进长度，计算β ;计算tl为动压序列yi的阈值y u时i的取值，t2为 动压序列Y1的阈值y t2时i的取值； S4，定义阈值 threshold 和 threshold_l，判断当 I Θ I < threshold 或者 E ( β + Θ ) -E ( β ) I < thresholcLl时，迭代结束，进行到下一步骤，同时得到最优的参数向 里 β (&〇，&1，···，···，，b，C); S5,保存曲线参数（aQ，ai，…，，…，aN，b，c)，提取包络。2. 根据权利要求1所述的一种基于示波法的血压包络的估计方法，其特征在于：所述 的yu的取值为y 1序列幅度的40%。3. 根据权利要求1所述的一种基于示波法的血压包络的估计方法，其特征在于：步骤 S2中所述的N取值为1。4. 根据权利要求1所述的一种基于示波法的血压包络的估计方法，其特征在于：步骤 S5中所述的M取值为4。5. 根据权利要求1所述的一种基于示波法的血压包络的估计方法，其特征在于：所述 的threshold的取值为5-25。6. 根据权利要求5所述的一种基于示波法的血压包络的估计方法，其特征在于：所述 的threshold的取值为10。7. 根据权利要求1所述的一种基于示波法的血压包络的估计方法，其特征在于：所述 的 threshold_l 的取值为初始值 Cat^a1、…、a;、...、aN、b、c)*0.005〇
【专利摘要】本发明一种基于示波法的血压包络的估计方法，在测量结束后先选取一个范围内的动压序列，定义拟合函数并设定参数向量的初始值，采用levenberg？marquardt算法的最小二乘原理寻找使得误差最小的参数向量，最终确定包络曲线。本发明提供了一种误差小、精度高的基于示波法的血压包络的估算方法。
【IPC分类】A61B5/0225
【公开号】CN105342590
【申请号】CN201510732538
【发明人】胡英俊, 袁哲明
【申请人】浙江漫思网络科技有限公司
【公开日】2016年2月24日
【申请日】2015年11月2日