CN101816556A量化心肺循环交互作用的频谱分析方法

CN101816556A量化心肺循环交互作用的频谱分析方法 (10)申请公布号101816556A (43)申请公布日xx.09.01101816556A*101816556A* (21)申请号xx10118641.6 (22)申请日xx.02.27A61B5/026(xx.01)A61B5/021(xx.01)G06F17/00(xx.01) (71)申请人李枝新地址中国台湾台北县 (72)发明人李枝新 (74)专利代理机构北京市浩天知识产权代理事务所11276代理人刘云贵 (54)发明名称量化心肺循环交互作用的频谱分析方法 (57)摘要本发明是一种用于量化心肺循环交互作用的频谱分析方法，利用与心脏相关的监测讯号进行心脏功能的评估，并透过撷取一固定频带内的频域讯号来进行，能避开频率较低且会影响心脏相关监测讯号的偶发干扰，因此无需将受测者进行全身麻醉或使之处于镇定状态下即能进行本发明，进而作为受测者心脏功能的评估指针。 (51)Int.Cl. (19)中华人民共和国国家知识产权局 (12)发明专利申请权利要求书2页说明书8页附图9页*101816556A*101816556A101816556A101816559A1/2页21.一种量化心肺循环交互作用的频谱分析方法，将一时间区间内的一动脉压讯号进行频谱分析，其步骤包含(a)将该动脉压讯号转换为一脉压差讯号；(b)将该脉压差讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带以得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值作为心脏功能的判断指针。 2.如权利要求1所述的频谱分析方法，其中该时间区间内至少包含10次的呼吸循环。 3.如权利要求1所述的频谱分析方法，其中于(a)步骤中，该脉压差讯号依下式所取得PPnorm(PP-PPmean)/PPmean其中，PPnorm为该脉压差讯号，PP为该动脉压讯号，PPmean为该动脉压讯号于该时间区间内的平均值。 4.如权利要求1所述的频谱分析方法，其中于(b)步骤中，利用快速傅立叶进行该脉压差讯号的频谱转换。 5.如权利要求1所述的频谱分析方法，其中于(b)步骤中，该功率频谱讯号经过一归一化步骤。 6.如权利要求1所述的频谱分析方法，其中于(c)步骤中，当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 7.一种量化心肺循环交互作用的频谱分析方法，将一时间区间内的一动脉压讯号进行频谱分析，其步骤包含(a)将该动脉压讯号转换为一射血量讯号；(b)将该射血量讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带以得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值作为心脏功能的判断指针。 8.如权利要求7所述的频谱分析方法，其中该时间区间内至少包含10次的呼吸循环。 9.如权利要求7所述的频谱分析方法，其中于(b)步骤中，利用快速傅立叶进行该射血量讯号的频谱转换。 10.如权利要求7所述的频谱分析方法，其中于(b)步骤中，该射血量讯号经过一归一化步骤。 11.如权利要求7所述的频谱分析方法，其中于(c)步骤中，当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 12.一种量化心肺循环交互作用的频谱分析方法，将一时间区间内的一血流量讯号进行频谱分析，其步骤包含(a)将该血流量讯号转换为一血量差讯号；(b)将该血量差讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带以得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频权利要求书101816556A101816559A2/2页3谱密度能量累加值作为心脏功能的判断指针。 13.如权利要求12所述的频谱分析方法，其中该时间区间内至少包含10次的呼吸循环。 14.如权利要求12所述的频谱分析方法，其中于(a)步骤中，该血量差讯号依下式所取得BFnorm(BF-BFmean)/BFmean其中，BFnorm为该血量差讯号，BF为该血流量讯号，BFmean为该血流量讯号于该时间区间内的平均值。 15.如权利要求12所述的频谱分析方法，其中于(b)步骤中，利用快速傅立叶进行该血量差讯号的频谱转换。 16.如权利要求12所述的频谱分析方法，其中于(b)步骤中，该功率频谱讯号经过一归一化步骤。 17.如权利要求12所述的频谱分析方法，其中于(c)步骤中，当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 18.一种量化心肺循环交互作用的频谱分析方法，将一时间区间内的一血流流速讯号进行频谱分析，其步骤包含(a)将该血流流速讯号转换为一流速差讯号；(b)将该流速差讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带而得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值作为心脏功能的判断指针。 19.如权利要求18所述的频谱分析方法，其中该时间区间内至少包含10次的呼吸循环。 20.如权利要求18所述的频谱分析方法，其中于(a)步骤中，该流速差讯号依下式所取得BFnorm(BF-BFmean)/BFmean其中，BFnorm为该流速差讯号，BF为该血流流速讯号，BFmean为该血流流速讯号于该时间区间内的平均值。 21.如权利要求18所述的频谱分析方法，其中于(b)步骤中，是利用快速傅立叶进行该流速差讯号的频谱转换。 22.如权利要求18所述的频谱分析方法，其中于(b)步骤中，该功率频谱讯号经过一归一化步骤。 23.如权利要求18所述的频谱分析方法，其中于(c)步骤中，当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 权利要求书101816556A101816559A1/8页4量化心肺循环交互作用的频谱分析方法技术领域0001本发明关于一种用于量化心肺循环交互作用的频谱分析方法，更特别的是关于一种将心脏相关的监测数据借由频谱分析方法来量化心肺循环交互作用。 背景技术0002血液动力学参数的相关监测是评估心脏功能的主要依据，也是受测者的生命征象之一，例如发生低血容量休克的始动环节主要就是血容量(BloodVolume，BV)减少、心输出量(Cardiac Output，CO)急剧减少或外周血管容量的扩大所造成的。 0003心脏的前负荷(preload)，亦即血容量，指心肌收缩之前所遇到的负荷，实际上是心室舒张末期容量或心室舒张末期室壁张力的反应。 前负荷是调节心搏出量(Stroke Volume，SV)的重要因素，而心搏出量亦与心输出量有关。 每分钟心脏的血液输出量等于每次心搏出量与每分钟心跳次数(Heart Rate，HR)的乘积，即COSV*HR。 0004前负荷增加，即心肌初长度增加，心肌收缩力增强，每搏量增加。 但如前负荷过度增加，超过心肌纤维最适初长度，则反而使心肌收缩力减弱，每搏量减少。 若使心脏扩张并增加心脏的血容量及增加足够的心肌纤维伸长量，根据史大林法则(Starlings Law)，下一次心脏的收缩力将增大，进而增加每分钟左心室射入主动脉的血量，即心输出量。 因此，血容量(即前负荷)是评估心脏功能的一项重要指针。 0005传统上用来评估血容量的心脏容积指针皆会有一些限制。 一种评估血容量的方法是利用压力来作为判断指标，主要将一中心静脉导管通过右心进入到肺动脉中，利用导管末端的压力计及连通管原理来量测左心房开口处的压力。 然而，实际上需要的是左心室的压力状况，以左心房的压力来预测左心室的压力，临床上并不准确，尤其是有合并心肺急慢性疾病的患者，其准确度会更差。 此外，在此种测量方式下，若要直接测量左心室的压力和容积就必需进入导管室进行心导管手术，此部份的费用不仅昂贵，在心导管手术过程中，操作者和受测者更将接受可观的放射线暴露量，且部分受测者因病况紧急更没有时间进行如此繁杂的手术程序。 0006传统上，也有利用顺流导管来间接测定左心房的压力以推估左心室的血容量，实际操作时，需将顺流导管由外围静脉漂入肺动脉的分枝，再以气球将肺动脉阻塞，此时所测得的压力即为脉动脉楔压(Pulmonary ArteryWedgePressure，PAWP)，可间接用来评估左心房的压力，但是在有肺动脉疾病或是肺病的患者可能会有较大的误差。 由于此操作方式，需将顺流导管送入脉动脉分支，导入的过程中需穿过两个心脏瓣膜，在没有放射线透视设备协助时，有一定的困难，而病房或是加护病房有放射线透视设备者非常少，即便有该设备，操作者和受测者也会接触大量的放射线暴露，因此应用在急重症患者上，执行顺流导管常遇到困难，且于临床上的研究显示，使用顺流导管治疗病患，不但对死亡率没有改善，更会增加病患的多种合并症，因此新的治疗指引，多已不再建议使用顺流导管。 此外，压力大并不一定代表血容量大，对于功能较差的心脏来说，较硬的心肌也会造成压力的增大。 再者，此种测量方式极易受待测者自律呼吸的干扰或是正压呼吸器中压力设定的干扰，因而说明书101816556A101816559A2/8页5仅适用于受测者在平静的状态下其量测结果才是正确的。 0007另一种评估的方法为利用热稀释技术(Transpulmonary thermodilution)来测定心脏容积，其主要利用一中心静脉导管和带有热敏探头的一股动脉导管透过一心肺容积监测仪(Peripherally IntroducedContinuous CardiacOutputMonitoring，PiCCO)来监测心输出量。 此种方式为一种心脏容积的测定，主要是于右心注入冰水，并于股动脉处量测血液温度，当心脏的容积越大，一定量的冰水对血液温度的影响就越小，所量测得的温度变化就越小。 然而此方法需要不断地打入冰水，平均每测定一次就需要在受测者的右心脏内打入1015毫升的冰水，对于长时间的监测来说其步骤相当繁琐。 此外，心脏容积大亦不一定代表血容量即越大，心肌收缩功能好的心脏虽然容积小，但血容量仍可达到正常值；反之，心肌收缩功能不佳且容积小的心脏尽管容积大，仍无法达到正常的血容量。 0008再一种评估的方法为利用动脉导管量测脉压差(收缩压与舒张压之间的差值)的变异量来作为心脏血容量大小的判断指针。 此一方法是利用呼吸时胸内压力的变化所造成心脏射血量(stroke volume)的变异(即心肺循环交互作用)来推估心脏功能的变化，心脏血容量较低时，心脏功能较差，心肺循环交互作用对心脏射血量的影响较大，反应在动脉的脉压差上，脉压差的变异量也会较大，故可用来推估心脏的功能和血容量，传统上是于监测时间的7.5秒内取动脉压的最大值与最小值以及动脉压的平均数，利用最大值减去最小值再除以平均数，当脉压差越大即代表心脏的血容量越小；反之，当脉压差越小即代表心脏的血容量越大。 然而此种方法最大的缺点就在于仅适用于全身麻醉或处于镇定状态下的受测者，因为此方法极易受到极端值的干扰，如咳嗽或自律呼吸，一旦极端值受此影响，对心肺循环交互作用的测定精确度将大幅降低。 0009此外其它的评估方法如使用心脏超音波，主要用来计算心脏腔室容积的大小，需要医师或是合格的超音波技术员来使用超音波并进行不断地量测，此种方式十分耗时和耗费人力，也难以达到连续性的监测。 还有一种评估方法如核医扫描，然其所运用的是放射线同位素，故部分受测者并不适宜使用。 0010现行作为心脏血容量的判断指针中，多需要将受测者全身麻醉或使之处于镇定状态，然而对处于休克的受测者进行全身麻醉或给予镇定又将是另一项高风险的行为。 目前使用的各种技术更有准确度、安全性、放射线暴露和执行困难等缺点。 因此，亟需一种适用于各种受测者而不需进行全身麻醉或镇定且能连续、快速又正确的心脏功能评估方法。 发明内容0011本发明主要目的在于提供一种用于量化心肺循环交互作用的频谱分析方法，利用与心脏相关的监测讯号进行心脏功能的评估，并利用突发性影响心肺循环交互作用的低频率特性(如咳嗽或自律呼吸，其发生的频率很低)，透过撷取一固定频带内的频域讯号来进行心脏功能的分析，而无需将受测者进行全身麻醉或使之处于镇定状态下即可进行。 0012本发明的另一目的在于利用心肺循环交互作用的量化结果来作为心脏功能的评估指针。 吸气时肺部会胀大而压迫到心脏，此时左心室受胸腔内正压影响会增加射血量，但右心室会因为此正压而减少射入肺动脉的血量，进而影响下一次左心室的射血量。 因此，呼气与吐气会造成胸内压的变化，而胸内压变化会影响心脏的收缩而造成每次心收缩搏出量的变化量发生变异，此变异会随着心脏心肌机能的强弱而有所差异，因此将心肺循环交互说明书101816556A101816559A3/8页6作用量化后即可作为相关心脏功能的评估指针。 0013为达上述目的，本发明提供一种用于量化心肺循环交互作用的频谱分析方法，将一时间区间内的一动脉压讯号进行频谱分析，其步骤包含(a)将该动脉压讯号转换为一脉压差讯号；(b)将该脉压差讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带以得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值是作为心脏功能的判断指针。 其中，该时间区间内至少包含10次的呼吸循环。 该脉压差讯号是依下式所取得PPnorm(PP-PPmean)/PPmean；PPnorm为该脉压差讯号，PP为该动脉压讯号，PPmean为该动脉压讯号于该时间区间内的平均值。 而当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 0014为达上述目的，本发明提供一种用于量化心肺循环交互作用的频谱分析方法，是将一时间区间内的一动脉压讯号进行频谱分析，其步骤包含(a)将该动脉压讯号转换为一射血量讯号；(b)将该射血量讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带以得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值是作为心脏功能的判断指针。 而当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 0015为达上述目的，本发明另提出一种量化心肺循环交互作用的频谱分析方法，将一时间区间内的一血流量讯号进行频谱分析，其步骤包含(a)将该血流量讯号转换为一血量差讯号；(b)将该血量差讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带以得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值是作为心脏功能的评估指针。 其中，该时间区间内至少包含10次的呼吸循环。 该血量差讯号是依下式所取得BFnorm(BF-BFmean)/BFmean；其中，BFnorm为该血量差讯号，BF为该血流量讯号，BFmean为该血流量讯号于该时间区间内的平均值。 而当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 0016为达上述目的，本发明再提出一种量化心肺循环交互作用的频谱分析方法，将一时间区间内的一血流流速讯号进行频谱分析，其步骤包含(a)将该血流流速讯号转换为一流速差讯号；(b)将该流速差讯号进行频谱转换以转换成一功率频谱讯号；(c)从该功率频谱讯号中撷取0.11.5Hz的频带而得到一功率频谱密度的分布，并将此频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值，该功率频谱密度能量累加值是作为心脏功能的评估指针。 其中，该时间区间内至少包含10次的呼吸循环。 该流速差讯号依下式所取得BFnorm(BF-BFmean)/BFmean；其中，BFnorm为该流速差讯号，BF为该血流流速讯号，BFmean为该血流流速讯号于该时间区间内的平均值。 而当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，心脏的血容量已充足。 附图说明0017图1为本发明于一实施例中的流程图。 0018图2为本发明第一实施例中动脉压讯号对时间的关系图。 说明书101816556A101816559A4/8页70019图3A至3C为本发明第一实施例中不同情况下脉压差的功率频谱密度分布图。 0020图4为本发明第一实施例中动脉血压讯号与射血量讯号对时间的关系图。 0021图5为本发明第一实施例中射血量功率频谱密度与脉压差功率频谱密度的关系图。 0022图6为本发明第二实施例中血流量讯号对时间的关系图。 0023图7A至7C为本发明第二实施例中不同情况下脉压差的功率频谱密度分布图。 0024图8为本发明第二实施例中6头小猪在放血实验中血流量功率频谱密度能量值的平方根与失血量的关系图。 0025图9为本发明第三实施例中血流流速讯号对时间的关系图。 0026图10A至10C为本发明第三实施例中不同情况下流速差的功率频谱密度分布图。 0027图11为本发明第三实施例中9头小猪在放血实验后的血流流速功率频谱密度能量值的平方根与心输出变化比例的关系图。 具体实施方式0028为充分了解本发明的目的、特征及功效，现借由下述具体的实施例，并配合所附的图式，对本发明做一详细说明，说明如后0029本发明的量化心肺循环交互作用的频谱分析方法可用于各种与心脏机能相关的心输出讯号，如动脉压(Arterial BloodPressure，ABP)、血流量、血流流速.等讯号，凡任何监测此等讯号的监测仪器皆可附加本发明的频谱分析方法，而可成为精确且多功能化的监测仪器。 0030请参阅图1，是本发明于一实施例中的流程图。 首先，输入一心输出讯号；接着将该心输出讯号归一化为一心输出差值讯号；再进行频谱转换；将转换后的一功率频谱讯号进行归一化；最后，撷取0.11.5Hz的频带得到一功率频谱讯号的分布，并将此频带内功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值而作为心脏功能的评估指针。 0031频谱转换的方式有很多种，可为离散傅立叶转换(Discrete FourierTransform，DFT)、快速傅立叶转换(Fast FourierTransform，FFT)、离散余弦转换(Discrete CosineTransform，DCT)、离散哈特利转换(Discrete HartleyTransform，DHT)、或离散小波转换(Discrete WaveletTransform，DWT)，上述转换方式并非本发明的限定，凡可用于时域转换成频域的转换方式皆可用于本发明。 于本发明的实施例中，利用韦尔奇法(Welchs Method)作为示例，于Matlab计算程序中进行快速傅立叶转换(FFT)以获得讯号的功率频谱密度(Power SpectrumDensity，PSD)的分布曲线。 Matlab计算程序中的函式运算于此不再赘述，仅列出韦尔奇法于Matlab中的函式，如式 (1)所示0032Pxx，fpwelch(xn，nfft，fs，window，noverlap) (1)0033其中，xn是讯号序列；nfft是快速傅立叶转换(FFT)的长度；fs是采样频率区间；window是选用的窗函数，须小于或等于nfft；noverlap是估算功率频谱密度时每一段重迭的长度，必须小于nfft；而输出的Pxx是估算出的功率频谱密度，f是频率坐标，Pxx与f分别为功率频谱密度分布图的纵坐标与横坐标。 0034在利用Matlab计算程序进行快速傅立叶转换以获得讯号的功率频谱密度的分布说明书101816556A101816559A5/8页8曲线之前，须先对输入的讯号进行一处理步骤，使输入的讯号中各讯号点间的时间间隔一致，而利于进行快速傅立叶转换。 于本发明的所有实施例中可利用三次曲线作内插法使该脉压差讯号、血流量差讯号、血流流速差讯号中每一点的时间间隔一致。 然而此三次曲线内插法仅为一种示例，其它如邻近点插值法、线性插值法、云形插值法.等皆可适用本发明。 上述内插法应用于Matlab计算程序的函式中，为Matlab程序的一般操作，于此不再赘述。 0035在取得该功率频谱讯号后，撷取约0.11.5Hz的频带来分析，以作为心脏功能的评估指针。 而就一般成人来说，较佳是撷取0.150.75Hz的频带，此乃每分钟呼吸9个循环至每分钟呼吸45个循环下的频带范围。 0036在本发明的第一实施例中将量化心肺循环交互作用的频谱分析方法应用于动脉压的监测上，动脉压的取得可利用现有的心肺容量监测仪或其它任何的动脉压监测装置。 0037首先，可透过任一动脉压监测仪器取得一时间区间内的一动脉压讯号，将该动脉压讯号利用归一化的方式转换为一脉压差(Pulse PressureVariation，PPV)讯号，转换的方式依式 (1)进行0038PPnorm(PP-PPmean)/PPmean (1)0039其中，PPnorm为经归一化后的该脉压差讯号；PP为该动脉压讯号；PPmean为该动脉压讯号于该时间区间内的平均值。 而为取得足够的动脉压讯号以供分析，该时间区间可包含2次以上的呼吸循环，而于一较佳实施态样下为至少包含10次的呼吸循环，于实际量测时大约为1分钟以上。 0040请参阅图2，为本发明第一实施例中动脉压讯号对时间的关系图。 图2仅为动脉压的一种示例，图中实线表示动脉压讯号，虚线表示经内插法处理的该脉压差讯号。 由图中可知该动脉压讯号具有规则的高低起伏，此现象即是受到胸内压变化的影响，而造成心收缩及心输出强度的变异。 0041接着，将该脉压差讯号利用前述的快速傅立叶转换进行频谱转换以转换成一功率频谱讯号，并经过一归一化步骤使功率频谱密度的能量大小归一化。 于一实施态样中，可将功率频谱密度对取样的该时间区间作归一化。 最后将所撷取频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值。 0042请参阅图3A至3C，为本发明第一实施例中不同情况下脉压差的功率频谱密度分布图。 图3A为经麻醉的受测者由呼吸器给予正压呼吸的脉压差功率频谱密度分布图；图3B为受测者未接呼吸器且未经麻醉并具有自律呼吸的脉压差功率频谱密度分布图；图3C为具有较大血容量的受测者的脉压差功率频谱密度分布图。 0043由图3B中可知即使受测者未受全身麻醉或未受镇定剂作用的情况下亦能明确呈现出一峰值。 图3A与图3B皆代表受测者的胸内压对心输出造成较大的影响，而图3C所示未包含任何明显的峰值。 0044该功率频谱讯号是量化心肺循环交互作用的结果，透过该功率频谱讯号的分布曲线可评估出受测者心脏的血容量。 由于本发明频谱分析方法系针对特定频率范围，因此对于频率很低的偶发性咳嗽或突发性的动脉压变化将不会落在如图3A至3C中所示的0.150.75Hz的频带范围内。 当呼吸频率较高时，于功率频谱密度图上即代表峰值越往高频率方向移动。 说明书101816556A101816559A6/8页90045由于相同的胸内压会对血容量较小的心脏产生较大的影响，因此在胸内压起伏间，脉压差的变异量较大者即代表心血容量较小，亦即前负荷较小；相反的，当心脏的血容量充足时，心肺循环交互作用的影响会减低，脉压差变异量也会下降。 功率频谱密度的能量值会随着脉压差变异量的大小而累加，即脉压差变异量越大则功率频谱密度的能量的累加值越大，再借由频带的选取，可为0.11.5Hz，而对于成人来说，于一较佳实施态样下是选取0.150.75Hz，即可获得该动脉压讯号的功率频谱密度分布图。 当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，代表脉压差变异量的累加值很小，即胸内压对该受测者的心输出影响较小，代表心脏血容量已相对充足，不随呼吸压力变异，心肺循环交互作用的影响较低。 0046在本发明的第一实施例中，该动脉压可转换为一射血量讯号，使该射血量讯号亦可进行本发明的频谱分析方法，其转换的公式如式 (2)00470048式 (2)中，SV是射血量，SA是动脉压图形下的积分值，SVR是外围血管阻力(systemic vascularresistance)，C(p)是血管系统(主要主动脉)的顺应性(pliance)，而dP/dt是动脉压图形在收缩期对时间的微分，cal是针对个别受测者的校正值。 由这个公式可将该动脉压图形转化为随着每次心搏(beat-by-beat)变化的射血量图形。 0049请参阅图4，为本发明第一实施例中动脉压讯号与射血量讯号对时间的关系图。 图中实线表示动脉压讯号，虚线表示经转换后再经内插法处理得到的该射血量讯号。 接着再经由上述的频谱分析方法亦可得到图3的结果。 0050接着请参阅图5，为本发明第一实施例中射血量功率频谱密度与脉压差功率频谱密度的关系图。 图中是以9只小猪的动脉压与射血量的频谱分析结果来比较，其相关系数(r-square)为0.95，代表动脉压与射血量间具有明显的线性关系，亦显示其可互相替换。 0051在本发明的第二实施例中是将量化心肺循环交互作用的频谱分析方法应用于血流量(blood flowrate)的监测，血流量的取得可利用现有的红外线监测仪或其它任何的血流量监测装置。 红外线监测仪利用血红素对红外线的吸收来监测血量。 0052如前所述，首先取得一时间区间内的一血流量讯号将该血流量讯号利用归一化的方式转换为一血量差讯号，转换的方式依式 (3)进行0053BFnorm(BF-BFmean)/BFmean (3)0054其中，BFnorm为经归一化后的该血量差讯号；BF为该血流量讯号；BFmean为该血流量讯号于该时间区间内的平均值。 而为取得足够的血流量讯号以供分析，该时间区间可包含2次以上的呼吸循环，而于一较佳实施态样下为至少包含10次的呼吸循环，于实际量测时大约为1分钟以上。 0055如图6所示，为本发明第二实施例中血流量讯号对时间的关系图。 红外线监测仪会将血红素对红外线的吸收转为电压讯号，如图6所示的纵坐标。 图中实线表示该血流量讯号，虚线表示经内插法处理的该血流量差讯号。 由图中可知该血流量讯号具有规则的高低起伏，此现象即是受到胸内压变化的影响，而造成心收缩及心输出强度的变异。 说明书101816556A101816559A7/8页100056同样地，接着将该血流量差讯号利用前述的快速傅立叶转换进行频谱转换以转换成一功率频谱讯号，并经过一归一化步骤使功率频谱的能量大小归一化。 于一实施态样中，可将经由能量值累加的功率频谱能量密度对取样的该时间区间作归一化，再借由频带的选取，可为0.11.5Hz，而就一般成人来说，较佳是撷取0.150.75Hz的频带，即可获得该血流量讯号的功率频谱能量密度分布图。 最后将所撷取频带内的功率频谱密度的所有能量值累加以得到一功率频谱密度能量累加值。 0057接着，请参阅图7A至7C，为本发明第二实施例中不同情况下血流量差的功率频谱密度分布图。 图7A是经麻醉的受测者由呼吸器给予正压呼吸的血流量差功率频谱密度分布图；图7B是受测者未接呼吸器且未经麻醉并具有自律呼吸的血流量差功率频谱密度分布图；图7C是具有较大血容量的受测者的血流量差功率频谱密度分布图。 而如前所述，当该功率频谱密度能量累加值小于或等于4.6210-4(分钟-1)时，代表血量差变异量的累加值很小，即胸内压对该受测者的心输出影响较小，代表心脏血容量已相对充足，不随呼吸压力变异，心肺循环交互作用的影响较低。 于图7中，在二倍主频处皆有一较小的能量波形(peak)，其乃频谱分析下产生的谐波能量。 0058接着请参阅图8，是本发明第二实施例中6头小猪在放血实验中血流量功率频谱密度能量值的平方根与失血量的关系图。 其中，将功率频谱密度能量值取平方根将有助于彰显血流量与失血量的关系。 于此，6头小猪在呼吸器支持及全身静注麻醉下接受放血实验并撷取0.150.75Hz频带的功率频谱密度，由图中可得知在0.150.75Hz频带下，功率频谱密度的能量值的平方根和失血量为呈线性相关。 而失血量可用来推估心脏的血容量(前负荷)，失血量愈多，即心脏的前负荷愈小时，心脏射出血量也随之愈小，血流量功率频谱密度的能量值平方根愈大，亦可反应心脏功能在放血过程中的变化。 因此，由图8可知，血流量功率频谱密度的能量值平方根的变化亦可显示失血程度的信息。 0059在本发明的第三实施例中是将量化心肺循环交互作用的频谱分析方法应用于血流流速(blood flowvelocity)的监测，血流流速的取得可利用血液在两个电极间流过所产生的电位差，或利用超音波多普勒监测仪，或其它任何的血流流速监测装置。 前述方法可不接触血液直接测量血流速。 一种超音波多普勒监测仪是对股动脉进行血流流速的监测。 0060如前所述，首先取得一时间区间内的一血流流速讯号将该血流流速讯号利用归一化的方式转换为一流速差讯号，转换的方式依式 (4)进行0061BFVnorm(BFV-BFVmean)/BFVmean (4)0062其中，BFVnorm为经归一化后的该流速差讯号；BFV为该血流流速讯号；BFVmean为该血流流速讯号于该时间区间内的平均值。 而为取得足够的血流流速讯号以供分析，该时间区间可包含2次以上的呼吸循环，而于一较佳实施态样下为至少包含10次的呼吸循环，于实际量测时大约为1分钟以上。 0063如图9所示，为本发明第三实施例中血流流速讯号对时间的关系图。 图中实线表示该血流流速讯号，虚线表示经内插法处理的该流速差讯号。 由图中可知该血流流速讯号具有规则的高低起伏，此现象即是受到胸内压变化的影响，而造成心收缩及心输出强度的变异。 00