专利名称：脉搏波传导速度的混沌计算方法
技术领域：
本发明涉及一种医学上对脉搏波传导速度的混沌计算方法。
背景技术：
所谓PWV，指的是脉搏波在动脉系统的两个既定点间的传播速度，PffV是测量沿 着动脉段的压力波速度，测量的理论基础取决于Moens-Korteweg公式PWV2 =(杨氏模 具.h/2rp)，p为血液密度(通常在1.05左右)，h/2r为管壁厚度/动脉内直径。杨氏模 具代表弹性材料硬度的物理参数。在实践中，PWV是由距离/两个脉搏测量点的传输时间计算出的。对其深入的实 验研究中，PWV能够测量在同一动脉同一路线上的位点(例如，大动脉)。然而，在一个无创 的临床研究中，这两个测量位点是在人体体表的可触摸到的体表的动脉(通常在颈动脉和 股动脉之间测量大动脉的PWV，颈动脉和肱动脉之间测量上肢PWV，颈动脉和桡动脉测量手 臂的PWV，股动脉和颈骨动脉测量大腿的PWV)；严格地说，这些位点不总是在脉搏传输的同 一线路上。这种技术在测量上的局限性是这种方法主要问题。除了年龄，动脉粥样僵硬度和动脉僵硬度，PWV值也随着血压的升高而升高，取决 于包括杨氏模具的血管壁的自然弹性。这是PWV检测的第二个问题。在对血管性能的实际 效果进行估算之前已排除影响血压的变化。尽管有这些限制因素，因为它的原理简单，PffV 检测方法正在被广泛的应用，而且PWV不受反射波的影响。水锤公式特征性阻抗(Zc)是另一项测量动脉僵硬度的指标，它能通过测量脉压和流速 波及其相关的频率组分获得。当PWV和Zc表示为cm/s或dyne, s/cm3公式如下Zc = 1. 330*PP (mmHg) /max F (cm/s) = PWV* P (P = 1.05)。这个被称为水锤公式，它与 PWV 及 压力和流量有直接的关系。根据Moens-Korteweg方程，PffV与弹性系数的平方根成正比，由于动脉弹性的减 低，脉搏波在动脉系统的传播速度加快了。近年来，由于PWV与动脉僵硬度的这种正相关 性，再加上它的测量方法简便易行而又没有创伤，所以，被广泛应用来作为评估动脉僵硬度 的一个指标。Siebenhofer等的研究结果更显示，经过培训的不同研究人员之间所测得的PWV 值几乎相同，显示了极好的可重复性。Asmar教授1995年采用全自动测量的Complior检验 了 PWV测量的重复性，结果进一步说明PWV检测具有极好的重复性。目前，这种测量方法已 经成为PWV检测的“金标准”。PWV与各种疾病的关系1)PWV与动脉硬化动脉僵硬度与年龄呈正相关，这主要是由于血管壁中层退行 性变，中层胶原含量增加，导致弹力层随着年龄的增长而断裂，伴有中层纤维化和钙化。其 它相关因素包括性别、吸烟、高血压、糖尿病、高胆固醇血症等。而高血压、动脉退行性变、心 血管疾病又加剧了动脉硬化的进程，这种变化尤以中心弹力血管，如主动脉为显著。动脉硬化又导致了大动脉缓冲功能障碍，引起脉压(PP)以及收缩压(SBP)不成比例地增加，由此， 引发了左室后负荷的增加以及冠状动脉缺血的加重，冠状动脉灌注压降低，收缩期压力跃 升，增加收缩末期左室应力，导致左室肥厚。2)PWV与终末期肾病(ESRD)患者心血管疾病危险性的预测ESRD患者是心血管疾 病的高危人群。流行病学调查及临床研究表明，大血管损害是ESRD患者心血管疾病高发病 率及死亡率的主要辅助因子，大血管疾病在尿毒症患者中进展尤为迅速。PWV在ESRD患者 中的预测作用不容忽视。3)PWV与糖尿病动脉硬度的增加是糖尿病患者的另一个特征，这可能是由于循 环中糖基化终末产物增加的缘故。而且，在II型糖尿病患者中，连续动脉血压监测发现，BP 升高是一个普遍现象。动脉硬度增加与糖尿病、高血糖、高胰岛素血症以及糖耐量异常有 关。在年长一些的成年人中，研究显示，胰岛素抵抗、高胰岛素血症以及高血糖可能促进了 动脉硬化的发生与发展，这种作用独立于年龄而存在。高血糖可以激活胶原蛋白分解，导致 蛋白多糖沉积于动脉壁，从而引起动脉壁的破坏、硬化。II型糖尿病增加某些类型大血管 疾病的危险性。在西方人群中，CHD是II型糖尿病患者最常见的死亡原因，糖尿病人群死 于CHD的危险性约增加了 3倍，非致死性心梗和充血性心衰的危险性增长幅度相类似。肾 脏病变是II型糖尿病的一个常见的微血管并发症。在长期研究中，确诊为糖尿病20年以 上的患者肾脏病变的累积发生率约为27% -60%，出现蛋白尿后10年ESRD的发生率约为 10% -35%，糖尿病是大部分西方人ESRD的主要原因。伴有糖尿病的ESRD患者心血管疾病 死亡率更高，这可能是由于糖尿病加剧了患者动脉硬化的程度。PWV是一个极强的心血管病 死亡以及总死亡的预测因子。关于PWV对ESRD患者心血管疾病发病风险的预测因子。PWV 在心血管疾病防治方面的应用前景与展望。4)PffV与高血压PWV和动脉壁僵硬度相关，较高的PWV表明了较高的僵硬度。为 了降低高血压病人的死亡率，应该将血压降低到何种程度，这是一个争论中的问题。有研究 表明，高血压病人血压的深度降低与较低的心血管事件相关联；将血压控制在严格的范围 内，可以降低II型糖尿病人的大血管和微血管并发症的危险，这种降低与治疗用药的种类 没有关系。提示深度降压可能对高血压病人的动脉壁有好处。一项近期的研究观察到，187 个经降压治疗后血压达到150/87mmHg的高血压病人在6年的随访期间有主动脉加速僵硬 的现象，表明如果血压降得更低一点，或许可以防止主动脉僵硬化。理论上，血压的急剧下 降会延缓脉搏波沿动脉壁的传播，从而降低PWV。那么，让血压在较长期间内维持在下降后 的水平是否会更好地改善PWV呢？ Atsuhiro Ichihara等研究了深度降压和适度降压对无 糖尿病的高血压病人的动脉壁僵硬度的长期影响。该项研究用测量baPWV(肱踝脉搏波传 导速度)的方法来评估动脉壁僵硬度。评估和比较了在降压治疗后3个月、6个月、9个月 和12个月的PWV变化(APWV)和血压变化(ΔΒΡ)的比值。结果表明，适度降压没有能降 低高血压病人的PWV，而深度降压显著地降低了高血压病人的PWV。虽然PWV反映了功能性 和解剖学动脉壁弹性，本项研究的结果却提示，如果深度降压至接近正常水平，并维持至少 12个月，解剖学弹性是可以改善的。因此，PWV不仅可以预报高血压病人的心血管死亡率和 发病率，还可以作为高血压长期治疗血压控制的指标。在心血管疾病中，近来特别提出了动脉僵硬度增大是血管损伤的病理状态，与动 脉粥样硬化型心血管病紧密相连。而脉搏波传导速度(PWV)是动脉硬度的指示器和动脉损
4伤的标志量。测量臂踝PWV(baP-WV)作为一个替代标志量，可以用于动脉粥样硬化型心血 管病的诊断和管理。欧洲著名实验室Framingham的研究表明，高血压发病率随年龄增高而增加，年龄 超过60岁后，脉压增加和舒张压降低是大动脉弹性降低的表现，高血压实际上加速了大动 脉僵硬的发展过程。早期收缩压增加是外周血管阻力增加造成，而晚期收缩压增加则主要 是大动脉弹性降低所致。左室射血产生对血管的压力，形成沿血管壁传播的压力搏动，即 PWV。脉搏波在动脉系统的某一分叉处会产生波反射，主要取决于反射系数和动脉僵硬度。 正常时，PWV相对较低，反射波叠加于主动脉压力波形的舒张期，这对机体是有益的。而随 年龄增高，动脉壁结构发生改变，动脉弹性降低，PWV明显增加。血压升高时，使动脉管壁张 力增大，长期高血压状态会导致血管壁结构改变。因此，衰老和高血压共同作用的结果是动 脉弹性降低。PWV增加、发射位置距离心脏较近时，导致收缩压峰值增大，舒张压降低，脉压 水平升高。平均舒张压降低，导致冠状动脉灌注降低，收缩晚期压力跃升，导致左室肥厚。PffV是判断与心脑血管疾病有密切关系的外周动脉壁硬化程度的指标。运用“当 动脉硬化时由心脏输出的血液产生的波动(脉搏波)的传导速度会加快”这一原理，测量 2次心跳之间的波动(脉搏波)传导速度，判断血管的弹性程度。外周动脉硬化的进展有 先于冠状动脉硬化的倾向，因此PWV可以作为将来患冠状动脉粥样硬化疾病风险的评估指 标，是心脑血管疾病的最佳预测指标。PWV对所有死因尤其心血管疾病的死亡是非常有意义 的独立预测因素。鉴于PWV与动脉硬化及高血压的关系，PWV测量近年来被广泛使用。Blacher等对 约530名高血压患者所做的研究提示，根据Framingham方程式，PffV与AA (氨基酸)的存 在及程度具有显著相关性，而且，PWV > 13mPs作为心血管疾病发生危险性的强预测因子具 有很高的应用价值。在没有动脉粥样硬化的患者中，随着PWV的升高，心肌梗死(Ml)、冠心 病(CHD)、CHD死亡率、中风、心血管疾病(CVD)以及CVD死亡率都随之上升。所有心血管疾 病的危险因子与PWV的升高呈线性相关。而且，在任何一个年龄，大动脉PWV都是心血管疾 病死亡率最好的预测因子。60岁以上人群中，PWV > 13-15m/s发生心血管病死亡的相对危 险度(OR值)为7. 11。Asmar等对1470名原发性高血压患者的横断面研究也表明，PWV的 升高与CVD的存在具有显著相关性，这种相关性独立于抗高血压治疗而存在。Boutouyrie等对1045例临床无心血管病的高血压患者进行了研究，在基线测量 患者的PWV，并根据性别、年龄、血压、胆固醇、糖尿病和吸烟情况计算Framingham危险评 分，评估冠心病发生危险，平均随访期为5. 17年。结果显示，在随访期共发生了 53次冠脉 事件和97次心血管事件。经单变量分析，PWV传播速度增加者，冠脉事件和心血管事件发病 危险增加；PWV增加一个标准差(3. 15m/s)，发生冠脉事件的相对危险为11. 42 (P < 0. 01)， 发生心血管事件的相对危险为11. 41 (P < 0. 001)。多变量分析显示，对Framingham危险评 分因素进行校正后，PWV增加一个标准差，发生冠脉事件的相对危险为11. 34 (P = 0. 039)； 对众所周知的危险因素进行校正后，发生冠脉事件的相对危险为11. 39 (P = 0. 01)。研究者 认为，PWV测定有助于确定冠心病高危病人，使这部分病人从更加积极的治疗中获益。综上所述，脉搏波传导速度PWV是一个较为灵敏的显示血管弹性的指标，在没有 发展为高血压之前，在动脉粥样硬化发展之中，就可能影响血管弹性，此时应用PWV测定仪 进行早期监测及干预，对老年心血管疾病的预防和治疗起到积极的监控作用。
在1999年的法国，PffV最初被认为是心血管疾病还包括晚期肾病患者患者(ESRD) 的独立的危害因素的前瞻性研究。在这项研究中，PWV彡12m/s的比率相对于PWV彡9. 4m/ s的死亡率为5. 4。这项队列研究在同一所研究机构另一篇文献也验证了 PWV的准确性。据 有ESRD的患者比没有患此病的患者具有更高的死亡率和发病率；因此，有相当数量的人关 注PWV检测心血管疾病的可预知性。然而，最近的同龄研究证实PWV检测的预见性价值。另一个问题就是PWV的预见性是否是具有潜在的心血管疾病的患者或者是否高 的PWV值代表患有心血管疾病。若干的研究暗示这两种现象都有可能。Cruickshank等人 建议PWV可能是心血管的危害因素综合指数(如，高血压或者糖尿病的持续时间)，可以预 测先天心血管疾病。相反的，明显的心血管疾病患者具有高的PWV值。除了对大规模人口进行研究获得的证据，另一个问题是两种方法在临床环境的可 用性。医师相对于常规的血压测量不容易获得PWA和PWV值。这个测量需要特殊的设备和 训练有素的操作者同时需要新的医疗技术。他们通常需要病人仰卧着进行测量。尽管有这 些缺点，当通过大规模的人研究的有效的方法和在临床环境中有用的测量变得清晰，PWA和 PWV将会是心血管医师的必要的工具。多项研究表明正常的人随着年龄的增长其主动脉PWV值会逐步升高，表现为主动 脉的衰退和僵硬度。这样的升高会使高血压患者和心脏病患者的病情加重。研究还表现为 桡神经的后收缩压得扩增的增长，颈动脉和主脉压随着年龄的增长而增加。相反的，弗莱明 起初的研究数据，证实主动脉PWV的增长与年龄的关系。Adji等人报道说心血管病人AIX 指数的增长打折可能归功于日常多种心血管疾病药物的作用。他们还指出AIX指数在药物 治疗的情况下很容易变化；然而，PWV却不是这样。AIX也是依靠流动波外形，如果有相对的 后收缩流动的减少，它可能不会随着年龄增加而升高。这些观点暗示AIX和PWV检测在临 床环境中具有不同的作用。运用PWV检测，临床医师可以了解血管僵硬度，可以反映病人的 患病史。通过测量AIX值可以评估疗效者的药物治疗和正常心室射血。降低的AIX值可以 减速动脉僵硬度的进程。另一个仍待解决的问题是设备，这两个指数都影响着心血管系统。一个大动脉疾 病的恶性循环被建议，其中高中央大动脉造成动脉扩张并加速弹性纤维断裂(图10)。这 个理论基于工程原理和人类提供的供观察的研究成果；然而，他还不能通过实验和人类的 研究提供直接的证明，因为过去的20年在弹性动脉中有弹性破裂需要八亿搏动的延伸。此 外，他还不清楚这个理论可以应用于肌性动脉；举例来说，冠状动脉和周围动脉的靶器官。 最近的研究表明，脉压与冠状动脉的破裂有关，而且比以前所认为的更加深入；因此，将来 应该有关于动脉僵硬度的发病率的调查。现代科技已经带来先进的工具以评估临床医学中主动脉血流和中央动脉僵硬度。 尽管很多问题尚待解决，正如100年前袖带式血压计那样，新方法有对心血管研究和实践 创新的潜力。对PWV的前瞻性的研究，会提供新的策略以治疗心血管疾病。传统PWV计算方法如下脉搏波的起始点是测量PWV的关键，起始点为脉搏波启动的时间点，传播时间T = 起始点1的时刻-起始点2的时刻=ta-tb，起始点ta和tb随测量对象位置的变化而变化， 如果是测量下肢的PWV值，则起始点ta为心电图波形的起始时间点，而起始点tb为股动脉脉搏波的起始时间点，PWV = L/T，L是所测两脉搏波的距离。一般的测量方法如下峰值比例法是先求出脉搏波的最高点(峰值)，然后以峰值的10%作为起始点。导数法就是以脉搏波的导数最大值作为起始点。切线交叉法是以脉搏波基线与上升曲线的交点作为起始点。拐点法是求脉搏波二次导数来确定起始点。而上述这些检测方法的共同缺点是噪声处理不够，而起始点位置的信号小，噪声 污染严重。混沌技术具有对微弱信号敏感及能进行有效的统计处理的特性，Duffing系统方 程是由Duffing (1918)，Moon和Holmes (1979，1980)提出的，方程原型是非线性单摆的受迫 振动方程，而混沌基本函数采用Duffing Moon Holmes系统模型χ" (t)+kx' (t) -ax (t) +b3 (t) = rcos (ω t)(3.0)
发明内容
本发明提供的一种脉搏波传导速度的混沌计算方法，利用混沌技术具有对微弱信 号敏感及能进行有效的统计处理的特性，采用混沌技术，使PWV测量精度及重复一致性大 大提高。为了达到上述目的，本发明提供一种脉搏波传导速度的混沌计算方法，其包含以 下步骤步骤1、分别计算起始点位置；步骤1. 1、计算心脏处脉搏波的起始点ta位置；步骤1. 1. 1、进行混沌算法；Y' (t)=J(t)Y(t)(3.1)J(t)是式3. 0所构造的Duffing Moon Holmes系统线性化后的状态描述矩阵。对于一个已知动力系统，它的基本解通过矩阵表达可以写作Y(t) = Q(t)R(t)(3.2)Q是一个正交矩阵，R是一个对角元素都为正数的上三角矩阵，代入式(3. 1)，则 得Q' (t)R(t)+Q(t)R' (t) = J (t) Q (t) R (t)(3.3)由矩阵性质得R ‘ (t) = (QT (t) J (t) Q (t))(3. 4)步骤1. 1. 2、采用Lyapimov算法进行混沌判断，判断系统是否发生混沌；根据式(3. 2)和(3. 4)得到Lyapimov特性指数为若混沌系统的最大Lyapimov特性指数小于零，说明系统是非混沌的，最大 Lyapunov特性指数大于零，是系统处于混沌状态的充要条件，当系统的最大Lyapimov特性 指数大于零时，此时的时间点“t”就是心脏处脉搏波的起始点ta。步骤1. 2、计算远离心脏的肢端处脉搏波的起始点tb位置；
步骤1. 2. 1、实测肢端处的脉搏波波形；步骤1. 2. 2、对步骤1. 2. 1得到的波形进行功率谱分析，得到频谱图；步骤1.2.3、对频谱图中的峰值进行分析，得到峰值处的频率f(角频率ω = 2 π f)禾口幅值r ；步骤1. 2. 4、将步骤1. 2. 3得到的角频率和幅值代入混沌函数，计算得到肢端处脉 搏波的起始点tb = t ；x〃 (t)+kx' (t) -ax (t) +b3 (t) = T1Cos (ω jt) +r2cos (ω 2t) +r3cos (ω 3t)(3. 6)其中ri、r2、r3、Q1, ω2、ω3是根据PWV相应波形经功率谱分析求得，根据混沌函 数的特性就可快速算出起始点的正确位置；步骤2、计算脉搏波从心脏传导到肢端所用的时间T = ta-tb ；步骤3、测量心脏到肢端之间的距离，即两脉搏波之间的距离L ；步骤4、计算得到PWV;PffV = L/T。本发明采用混沌技术，使PWV测量精度及重复一致性大大提高。
图1是本发明的具体实施例中提供的实测的股动脉脉搏波的波形图
图2是对图1的波形进行功率谱分析得到的频谱波形图3是一组实施例中系统没有发生混沌时的实测波形图3A是对应式(3.6)的相平面图3B是脉搏波的波形图3C是对应式(3.6)的时序响应曲线；
图4是另一组实施例中系统没有发生混沌时的实测波形图4A是对应式(3.6)的相平面图4B是脉搏波的波形图4C是对应式(3.6)的时序响应曲线；
图5是实施例中系统发生混沌时的实测波形图5A是对应式(3.6)的相平面图5B是脉搏波的波形图5C是对应式(3.6)的时序响应曲线；
图6是本发明实的具体施例中提供的实测的脉搏波心电图波形。
具体实施例方式以下根据图1-图6，具体说明本发明的较佳实施方式一种脉搏波传导速度的混沌计算方法，其包含以下步骤步骤1、分别计算起始点位置；步骤1. 1、计算心脏处脉搏波的起始点ta位置；步骤1. 1. 1、进行混沌算法；
Y' (t)=J(t)Y(t)(3.1)J(t)是式3. 0所构造的Duffing Moon Holmes系统线性化后的状态描述矩阵。对于一个已知动力系统，它的基本解通过矩阵表达可以写作Y(t) = Q(t)R(t)(3.2)Q是一个正交矩阵，R是一个对角元素都为正数的上三角矩阵，代入式(3. 1)，则 得Q' (t)R(t)+Q(t)R' (t) = J(t)Q(t)R(t)(3.3)由矩阵性质得R ‘ (t) = (QT (t) J (t) Q (t))(3. 4)步骤1. 1. 2、采用Lyapimov算法进行混沌判断，判断系统是否发生混沌；根据式(3. 2)和(3. 4)得到Lyapimov特性指数为
(3.5)图6是实测的脉搏波心电图波形；如图3和图4所示，系统还未发生混沌；如图5所示，系统发生混沌，混沌发生点的位置如图5B所示，在图5A中，此时对应 的脉搏波的起始点为波值=0. 02(波值是脉搏波各点对应的电压值)，这时由式(3. 5)求出 的Lyapimov特性指数值为最大值，如图5B所示，相应时间点为18965，将其乘上采样时间间 隔 0. 003s,即为 56. 895 秒。即 ta = 56. 895 秒；步骤1. 2、计算远离心脏的肢端处脉搏波的起始点tb位置；步骤1. 2. 1、实测肢端处的脉搏波波形(如图1所示)；步骤1. 2. 2、对步骤1. 2. 1得到的波形进行功率谱分析，得到频谱图(如图2所 示)；步骤1.2. 3、对频谱图中的峰值进行分析，得到峰值处的频率f(角频率ω = 2 Jif)和幅值r (如图2所示)；Γι = 0. 0245 ；r2 = 0. 0175 ；r3 = 0. 009 ；= 9 ； f2 = 18 ；f3 = 26 ；ω 1 = 2 η T1 ； ω 2 = 2 η f2 ； ω 3 = 2 η f3 ；步骤1. 2. 4、将步骤1. 2. 3得到的角频率和幅值代入混沌函数x〃 (t)+kx' (t) -ax (t) +b3 (t) = T1Cos (ω jt) +r2cos (ω 2t) +r3cos (ω 3t)(3.6)计算得到肢端处脉搏波的起始点tb = 56. 782秒；步骤2、计算脉搏波从心脏传导到肢端所用的时间T = ta-tb = 56. 895-56. 782 = 0. 113 秒；步骤3、测量心脏到肢端之间的距离，即两脉搏波之间的距离L = 0. 9米；步骤4、计算得到PWV;PWV = L/T = 0. 9/0. 113 = 7. 96。
权利要求
一种脉搏波传导速度的混沌计算方法，其特征在于，包含以下步骤步骤1、分别计算起始点位置；步骤2、计算脉搏波从心脏传导到肢端所用的时间T＝ta tb；步骤3、测量心脏到肢端之间的距离，即两脉搏波之间的距离L；步骤4、计算得到PWV；PWV＝L/T。
2.如权利要求1所述的脉搏波传导速度的混沌计算方法，其特征在于，所述的步骤1包 含以下步骤步骤1. 1、计算心脏处脉搏波的起始点ta位置；步骤1. 2、计算远离心脏的肢端处脉搏波的起始点tb位置。
3.如权利要求2所述的脉搏波传导速度的混沌计算方法，其特征在于，所述的步骤1.1 包含以下步骤步骤1. 1. 1、进行混沌算法； Y' (t) = J(t)Y(t)(3. 1)J(t)是式3. 0所构造的Duffing Moon Holmes系统线性化后的状态描述矩阵。 对于一个已知动力系统，它的基本解通过矩阵表达可以写作 Y(t) = Q(t)R(t)(3.2)Q是一个正交矩阵，R是一个对角元素都为正数的上三角矩阵，代入式(3. 1)，则得 Q' (t)R(t)+Q(t)R' (t) = J(t)Q(t)R(t)(3.3)由矩阵性质得R' (t) = (QT(t)J(t)Q(t))(3.4)步骤1. 1. 2、采用Lyapimov算法进行混沌判断，判断系统是否发生混沌；根据式(3. 2)和(3. 4)得到Lyapunov特性指数为 若混沌系统的最大Lyapimov特性指数小于零，说明系统是非混沌的， 最大Lyapimov特性指数大于零，是系统处于混沌状态的充要条件，当系统的最大 Lyapimov特性指数大于零时，此时的时间点“t”就是心脏处脉搏波的起始点ta。
4.如权利要求2所述的脉搏波传导速度的混沌计算方法，其特征在于，所述的步骤1.2 包含以下步骤步骤1. 2. 1、实测肢端处的脉搏波波形；步骤1. 2. 2、对步骤1. 2. 1得到的波形进行功率谱分析，得到频谱图； 步骤1. 2. 3、对频谱图中的峰值进行分析，得到峰值处的频率f (角频率ω =2Jif)和幅值r ； 步骤1. 2. 4、将步骤1. 2. 3得到的角频率和幅值代入混沌函数，计算得到肢端处脉搏波 的起始点tb = t ；x〃 (t)+kx' (t) -ax (t) +b3 (t) = T1Cos (ω jt) +r2cos (ω 2t) +r3cos (ω 3t) (3. 6) 其中rpi^iv ωι、ω2、ω3是根据PWV相应波形经功率谱分析求得，根据混沌函数的 特性就可快速算出起始点的正确位置。
全文摘要
一种脉搏波传导速度的混沌计算方法，利用混沌技术具有对微弱信号敏感及能进行有效的统计处理的特性，采用混沌技术，使PWV测量精度及重复一致性大大提高。
文档编号G01R23/16GK101897574SQ20091005731
公开日2010年12月1日 申请日期2009年5月26日 优先权日2009年5月26日
发明者张琴舜 申请人:上海艾康菲医疗器械技术有限公司