- 浏览: 497109 次
- 性别:
- 来自: 深圳
文章分类
- 全部博客 (185)
- job (15)
- linux/windows/unix/bash/shell (31)
- JAVA/J2EE/spring/hibernate/struts (30)
- VC/C++ (48)
- mysql/postgresql (6)
- php/jsp/asp/pear (1)
- FMS/flex/openlaszlo/red5/openmeetings (34)
- apache/tomcat/ftp/svn (6)
- xen/vm/Hadoop/cloudcompute (6)
- visual studio/eclipse/zendstudi/ant (8)
- others (1)
- windows异常处理 __try __except (1)
- (1)
- matlab (4)
- android (0)
最新评论
-
hongzhounlfd:
很透彻,很详细
依赖注入和控制反转 -
jefferyqjy:
谢谢~言简意赅~很明了!
依赖注入和控制反转 -
elderbrother:
太好了,谢谢
依赖注入和控制反转 -
east_zyd_zhao:
终于搞明白了
依赖注入和控制反转 -
Dremeng:
完美,一看就懂理解透彻
依赖注入和控制反转
matlab
<div id="wrapper">
<div id="bodycontainer" class="document">
<div id="resources">
<div class="noindex">
<div id="flakes">
<div class="fnavbox">
<ul></ul>
</div>
</div>
</div>
</div>
<div id="primarycontentcontainer">
<div class="noindex">
<h1>利用LabVIEW进行心率变异分析</h1>
<div class="docdetails">
<div class="ratings">2 评级 | <strong>3.00</strong> out of 5</div>
<!-- Read in other language section --><strong>阅读语言 </strong><select><option selected value="http://zone.ni.com/devzone/cda/epd/p/id/6275">Chinese (Simplified)</option>
<option value="http://zone.ni.com/devzone/cda/epd/p/id/5832">English</option></select> | <a><img class="bottomalign" src="http://baitai.iteye.com/images/global/us/icons/print.gif" border="0" alt="" width="18" height="13"> 打印</a> </div>
</div>
<div id="pagearea">
<div class="noindex">
<h3>概览</h3>
</div>
<p>心率变异(HRV)分析是一项在研究和临床中广泛使用的重要应用。LabVIEW包含多个信号处理VI,使用户可以非常方便地创建自定义的HRV分析应用。本文介绍了利用LabVIEW进行HRV分析的应用与实践。</p>
<div class="noindex">
<h3>目录</h3>
<ol>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc0">引言</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc1">什么是HRV?</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc2">采集RR间隔</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc3">HRV分析方法</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc4">应用</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc5">获取HRV分析起步工具</a> </li>
<li><a href="http://baitai.iteye.com/admin/blogs/#toc6">参考文献</a></li>
</ol>
<h3>下载</h3>
</div>
<a name="toc0"></a>
<h2><strong>引言</strong></h2>
<p align="left"><strong>心电图(ECG)</strong>是一种由心电图仪产生的图形,提供了有关个人心脏健康的信息。ECG常在临床应用中帮助医生诊断心跳过速等心脏疾病。除了直接分析ECG信号,研究人员和医生们还可以从ECG信号中提取其他非直接的测量信息,最常见的测量之一便是<strong>心率变异(HRV)</strong>。本文描述了如何利用NI LabVIEW图形化开发环境提取并分析HRV测量结果,并同时介绍了LabVIEW HRV分析起步工具—— 一款基于LabVIEW的面向HRV分析的程序。本文所展示的HRV分析结果来自于该HRV分析起步工具。</p>
<a name="toc1"></a>
<h2><strong>什么是HRV?</strong></h2>
<p align="left">一个正常的单周期ECG信号由若干个波组成,如图1所示。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57496.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图1</strong>. ECG信号的R波峰与RR间隔</p>
<p align="left">具有最高峰值的波被称为R波。RR间隔是指两个相邻R波之间的时间间隔。波峰较小的波分别为P波、T波和U波。</p>
<p align="left">心率变异(HRV)测量可以分析出这些RR间隔(它们显示了连续心跳之间的变化)是如何随时间而变化的。例如,图2展示了一个典型的成人的HRV。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57536.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图2.</strong> 一个成人的HRV</p>
<p align="left">图2中的HRV分析显示了该RR间隔范围从0.6秒变化至1秒,其平均值为0.8秒。</p>
<p align="left">HRV分析提供了<strong>自主神经系统(ANS)</strong>的一个定量标记,因为HRV的调节机制源自交感和副交感神经系统[5]。迄今为止,许多研究人员和工程师们都致力于发掘HRV分析中对临床应用有价值的信息。在过去的十年中,已经发表了超过2000篇关于HRV的文章[1]。这些文章考察了HRV与血压、心肌梗塞、神经系统、心律失常、糖尿病、呼吸、肾衰竭、性别、年龄、疲劳程度、药物、吸烟、酒精等因素的关系[2]。目前,多个医学学科都已经在对HRV进行研究,相信HRV分析在未来的患者图形与数据分析中的应用将越来越多。</p>
<a name="toc2"></a>
<h2><strong>采集RR间隔</strong></h2>
<p align="left">为了进行HRV分析,首先必须获取RR间隔。图3展示了<a id="OLE_LINK2" name="OLE_LINK2"></a> <a id="OLE_LINK1" name="OLE_LINK1">RR间隔的采集过程</a>。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image22756.jpg" alt="" width="513" height="88"></p>
<p align="center"><strong>图3.</strong> RR间隔的采集过程</p>
<a name="h30"></a>
<h3><strong>采集原始ECG信号</strong></h3>
<p align="left">ECG信号可通过各种心电图仪采集获得。这些仪器的带宽并不要求很高。可以利用NI多通道DAQ设备,如NI PCI-6221,从ECG记录设备的输出端采集原始ECG信号。典型的采样率设置为125 Hz或250 Hz。所采集的ECG信号可以以NI TDMS文件类型存储以便进行离线分析。一些在线数据库(如<a href="http://www.physionet.org/physiobank/database/#ecg">MIT-BIH数据库database</a>)也提供许多典型的ECG信号。该HRV分析起步工具包含一个导入<a href="http://www.physionet.org/physiobank/database/#ecg">MIT-BIH</a> ECG<a href="http://www.physionet.org/physiobank/database/#ecg">数据</a>的向导。</p>
<a name="h31"></a>
<h3><strong>从ECG信号中提取RR间隔</strong></h3>
<p align="left">从原始的ECG信号中可以提取到RR间隔。该提取过程通常包括预处理步骤和波峰检测步骤。如果原始ECG信号具有噪声和显著的基线漂移趋势,就有必要对其进行预处理。然后,您可以通过阈值设置或利用基于小波的波峰检测方法,检测R波峰,从而计算RR间隔。</p>
<p align="left">请参阅<a href="http://zone.ni.com/devzone/cda/tut/p/id/10041">NI开发者园地</a>,获取关于ECG信号处理和特性提取的更多信息。</p>
<a name="toc3"></a>
<h2><strong>HRV分析方法</strong></h2>
<p align="left">HRV分析有不同的方法[1][2]。方法之一便是时域分析,即利用暂态的RR间隔信号提取一些特别的度量。另一个方法是频谱分析,该方法以某个特定的速率对RR间隔进行插值处理,并将该间隔信号变换到频域。针对这两种方法,目前有一些标准可供使用[1]。还有一些其他方法,如非线性分析和时频分析,也有相关讨论和研究。</p>
<p align="left">图4展示了HRV分析的整个过程。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image22758.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图4.</strong> HRV分析的整个过程。</p>
<p align="left">LabVIEW具有多个可用于HRV分析的工具包或函数。如图5所示的HRV分析起步工具,展现了利用LabVIEW进行HRV分析的简便性。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57500.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图5.</strong> HRV分析起步工具的前面板</p>
<a name="h32"></a>
<h3><strong>预处理</strong></h3>
<p align="left">RR间隔信号可能存在一些基线抖动,这会导致不精确的暂态分析和频率分析。您可以选择是否在分析前对RR间隔信号进行去漂移处理,如图6所示。对RR间隔信号进行去漂移处理的方法,与<a href="http://zone.ni.com/devzone/cda/tut/p/id/10041">使用LabVIEW进行ECG信号处理</a>文档中所提及的方法相似。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57502.jpg" alt="" width="364" height="203"></p>
<p align="center"><strong>图6.</strong> 预处理配置对话框</p>
<p align="left">图7展示了去漂移处理的过程。图中的红线表示该基线漂移。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57504.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图7.</strong> RR间隔信号的预处理</p>
<a name="h33"></a>
<h3><strong>分析RR间隔信号</strong></h3>
<p align="left">下面描述了用于分析RR间隔信号的不同方法。HRV分析起步工具支持所有这些方法。</p>
<a name="h40"></a>
<h4><strong>时域分析</strong></h4>
<p align="left">对于时间序列分析,通常使用时域量度。许多量度可以从最初的RR间隔信号中提取,以展现ANS中的变化。表1列出了利用HRV分析起步工具提取的量度。</p>
<p align="center"><strong>表1.</strong> HRV的时域量度[1][5]</p>
<table border="1" cellspacing="0" cellpadding="0" align="center"><tbody>
<tr>
<td width="139" valign="top">
<p align="center"><strong>变量</strong></p>
</td>
<td width="54" valign="top">
<p align="center"><strong>单位</strong></p>
</td>
<td width="397" valign="top">
<p align="center"><strong>描述</strong></p>
</td>
</tr>
<tr>
<td colspan="3" width="590" valign="top">
<p align="center">统计学量度</p>
</td>
</tr>
<tr>
<td width="139" valign="top">R<a id="OLE_LINK10" name="OLE_LINK10"></a> <a id="OLE_LINK9" name="OLE_LINK9">R均值与标准差</a>
</td>
<td width="54" valign="top">s</td>
<td width="397" valign="top">所有RR间隔的均值与标准差</td>
</tr>
<tr>
<td width="139" valign="top">HR均值与标准差</td>
<td width="54" valign="top">1/min</td>
<td width="397" valign="top">所有心率的均值与标准差</td>
</tr>
<tr>
<td width="139" valign="top">RMSSD</td>
<td width="54" valign="top">ms</td>
<td width="397" valign="top">各个相邻RR间隔的差值的均方根</td>
</tr>
<tr>
<td width="139" valign="top">NN50计数</td>
<td width="54" valign="top"> </td>
<td width="397" valign="top">在所有测量值中,差值超过50ms的相邻RR间隔对的数目</td>
</tr>
<tr>
<td width="139" valign="top">pNN50</td>
<td width="54" valign="top">%</td>
<td width="397" valign="top">NN50计数除以所有RR间隔的总数目</td>
</tr>
<tr>
<td colspan="3" width="590" valign="top">
<p align="center">几何学量度</p>
</td>
</tr>
<tr>
<td width="139" valign="top">HRV三角指数</td>
<td width="54" valign="top"> </td>
<td width="397" valign="top">所有RR间隔的总数除以所有RR间隔的直方图的高度</td>
</tr>
<tr>
<td width="139" valign="top">TINN</td>
<td width="54" valign="top">ms</td>
<td width="397" valign="top">对于间距为7.8125 ms(1/128 s)的离散标度所测得的所有RR间隔,对其直方图中最高峰的最小均方误差三角插值的基线宽度。</td>
</tr>
</tbody></table>
<p align="left">LabVIEW中包含概率与统计VI,您可以利用这些VI方便地获取统计学度量,如图8所示。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57538.jpg" alt="" width="430" height="341"></p>
<p align="center"><strong>图8.</strong> LabVIEW中的概率与统计VI</p>
<p align="left">图9展示了一个成人的RR间隔信号样本。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57512.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图9.</strong> 利用时域分析方法获得的RR间隔信号采样</p>
<p align="left">图10展示了利用HRV分析起步工具中的时域分析方法得到的这些信号的统计学结果和几何学结果。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57514.jpg" alt="" width="364" height="211"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57516.jpg" alt="" width="364" height="265"></p>
<p align="center"><strong>图10.</strong> 利用时域分析方法的HRV分析</p>
<a name="h41"></a>
<h4><strong>频域分析</strong></h4>
<p align="left">HRV分析的频域分析方法从RR间隔信号中提取频域参数,如峰值频率和带内功率。ANS的两个组成部分——交感神经系统和副交感神经系统——会决定升高或降低心率,并影响RR间隔频谱中的不同频带。因而,可以利用频域分析监测ANS的状态。</p>
<p align="left">表2展示了HRV的常见频域度量,其中包括峰值频率以及超低频(VLF)、低频(LF)和高频(HF)等频带的带内功率。</p>
<p> </p>
<p align="center"><strong>表2.</strong> HRV的频域量度[1][5]</p>
<table border="1" cellspacing="0" cellpadding="0" align="center"><tbody>
<tr>
<td width="139" valign="top">
<p align="center"><strong>变量</strong></p>
</td>
<td width="54" valign="top">
<p align="center"><strong>单位</strong></p>
</td>
<td width="397" valign="top">
<p align="center"><strong>描述</strong></p>
</td>
</tr>
<tr>
<td width="139" valign="top">峰值频率</td>
<td width="54" valign="top">Hz</td>
<td width="397" valign="top">针对VLF、LF和HF频带的功率谱密度(PSD)估计的峰值频率</td>
</tr>
<tr>
<td width="139" valign="top">VLF </td>
<td width="54" valign="top">ms<sup>2</sup>
</td>
<td width="397" valign="top">0–0.04 Hz<a id="OLE_LINK20" name="OLE_LINK20"></a> <a id="OLE_LINK19" name="OLE_LINK19">频段的功率</a>
</td>
</tr>
<tr>
<td width="139" valign="top">LF </td>
<td width="54" valign="top">ms<sup>2</sup>
</td>
<td width="397" valign="top">0.04–0.15 Hz频段的功率</td>
</tr>
<tr>
<td width="139" valign="top">HF</td>
<td width="54" valign="top">ms<sup>2</sup>
</td>
<td width="397" valign="top">0.15–0.4 Hz频段的功率</td>
</tr>
<tr>
<td width="139" valign="top">LF归一化</td>
<td width="54" valign="top">无量纲</td>
<td width="397" valign="top">
<a id="OLE_LINK30" name="OLE_LINK30"> </a><a id="OLE_LINK29" name="OLE_LINK29">LF功率(采用归一化单位):LF/(总功率–VLF)*100</a>
</td>
</tr>
<tr>
<td width="139" valign="top">HF 归一化</td>
<td width="54" valign="top">无量纲</td>
<td width="397" valign="top">HF功率(采用归一化单位):HF/(总功率–VLF)*100。</td>
</tr>
<tr>
<td width="139" valign="top">LF/HF比</td>
<td width="54" valign="top"> </td>
<td width="397" valign="top">LF/HF(均采用ms<sup>2</sup>作为单位)</td>
</tr>
</tbody></table>
<p align="left">频域方法通常涉及以下三个步骤:</p>
<ol>
<li>对RR间隔信号进行重采样。 </li>
<li>估算该RR间隔信号的功率谱密度(PSD)。 </li>
<li>根据该PSD计算频域参数。 </li>
</ol>
<p align="left">从ECG信号中抽取的RR间隔信号通常不是以一个规则间隔采样所得。因此,该RR间隔信号并不是延时间轴等距离间隔。可以利用位于<strong>数学>>内插与外推</strong>函数选板中的一维插值VI,将最初的RR信号转换为均匀采样的时间序列。</p>
<p align="left">对于PSD估计,LabVIEW包含可以执行非参数估计方法或参数估计方法的VI。非参数方法包括基于FFT的PSD,如周期图法和<a id="OLE_LINK34" name="OLE_LINK34"></a> <a id="OLE_LINK33" name="OLE_LINK33">沃尔什法</a>。参数方法包括基于模型的PSD,如自回归(<a id="OLE_LINK36" name="OLE_LINK36"></a> <a id="OLE_LINK35" name="OLE_LINK35">AR)方法</a>和自回归滑动平均(ARMA)方法。<a id="OLE_LINK38" name="OLE_LINK38"></a> <a id="OLE_LINK37" name="OLE_LINK37">HRV分析起步工具</a>包含沃尔什法和AR频谱方法。LabVIEW高级信号处理工具包中则包含有其他方法。</p>
<p align="left">图11展示了利用HRV分析起步工具中的频域分析方法得到的一个RR间隔信号分析结果。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57518.jpg" alt="" width="436" height="257"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57520.jpg" alt="" width="436" height="257"></p>
<p align="center"><strong>图11.</strong> 利用频域分析方法的HRV分析</p>
<a name="h42"></a>
<h4><strong>联合时频域分析</strong></h4>
<p align="left">频域分析方法适用于随时间变化不大的RR间隔信号分析。然而,在长期HRV分析(如24小时RR间隔与药物药效学研究)等应用中,RR间隔信号会随时间显著变化。在这类情况下,您可以使用联合时-频分析(JTFA)方法分析稳态和时-频行为。</p>
<p align="left">联合时频域方法(JTFA)没有标准的度量参数。因此可以利用短时傅立叶变换(STFT)、Gabor展开或连续小波VI得到时频谱图并进行定性分析。文献[8]中提到了一个利用JTFA方法的药物药效学应用案例。</p>
<p align="left">高级信号处理工具包支持JTFA方法。图12显示了利用该工具包所含的LabVIEW代码,该代码可以计算得到RR间隔信号的STFT谱图。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57522.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图12.</strong> 获取RR间隔信号的STFT谱图</p>
<p align="left">图13展示了利用HRV分析起步工具中的JTFA方法获得的联合时频分析结果。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57524.jpg" alt="" width="278" height="182"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57526.jpg" alt="" width="278" height="182"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57528.jpg" alt="" width="278" height="183"></p>
<p align="center"><strong>图13.</strong> 利用JTFA方法的HRV分析</p>
<a name="h43"></a>
<h4><strong>非线性分析</strong></h4>
<p align="left">RR间隔信号是非线性的,因为它们来源于血流动力学、电生理学和体液参数、以及自主神经调节与中枢神经调节之间的复杂交互关系[1]。关于非线性分析的研究基本处于空白,所以没有任何相关标准。然而,这是一个充满希望的研究领域,我们需要在该领域开发出一种工具,从而帮助我们更好地理解人体系统的复杂性。</p>
<p align="left">HRV分析起步工具支持<a id="OLE_LINK40" name="OLE_LINK40"></a> <a id="OLE_LINK39" name="OLE_LINK39">庞加莱图</a>和<a id="OLE_LINK50" name="OLE_LINK50"></a> <a id="OLE_LINK49" name="OLE_LINK49">去漂移波动分析</a>(DFA)方法。下面将分别介绍这些方法。</p>
<p align="left"><strong>庞加莱图</strong></p>
<p align="left">庞加莱图是一种XY图,其X轴和Y轴表示两个连续的RR间隔。该图所用的模型是一个椭圆,该椭圆的对称线与X轴成45◦,如图14所示。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57530.jpg" alt="" width="390" height="247"></p>
<p align="center"><strong>图14.</strong> RR间隔的庞加莱图</p>
<p align="left">图14显示了下列数值:</p>
<ul>
<li>
<strong>Sd1</strong>—各点到对称线距离的标准差。Sd1刻画出了短期变化,这主要是由呼吸性窦性心律不齐(RSA)引起的。 </li>
<li>
<strong>Sd2</strong>—各点在沿着对称直线方向上的标准偏差。Sd2刻画了长期变化[2]。 </li>
</ul>
<p align="left">LabVIEW包含概率与统计VI和拟合VI,以及一个前面板XY图控件,可以利用该控件创建并显示一幅庞加莱图。图15显示了<strong>拟合</strong>VI。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57532.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图15.</strong> LabVIEW所包含的拟合VI库</p>
<p align="left"><strong>去漂移波动分析</strong></p>
<p align="left">去漂移波动分析(DFA)方法定量表示了RR间隔信号的分形尺度属性。该方法通过剔除不稳定的人为趋势,避免了错误的相关性。因此,相比其他方法,DFA方法能更好展示HRV长期相关性的变化情况。</p>
<p align="left">DFA方法利用线性拟合对DFA图进行分析。拟合直线的斜率α表示了HRV的相关性。如果α=0.5,那么该信号为白噪声,没有任何相关性。在一些情况下,DFA图不是严格地呈线性,而是由相交于某点的两条直线组成。该情形表明存在两个因子:一个短期内的小尺度指数和一个长期时间内的大尺度指数。</p>
<p align="left">图16显示了利用HRV分析起步工具中的DFA方法获得的一个RR间隔信号分析结果。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57534.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图16.</strong> RR间隔信号的DFA图</p>
<a name="toc4"></a>
<h2><strong>应用</strong></h2>
<p align="left">HRV分析结果可以应用于临床和研究。对于统计与频域方法,HRV分析结果[1]提出了一些临床应用HRV分析以帮助诊断各种心脏疾病的标准。利用HRV分析起步工具,研究人员能够发现HRV与疾病以及药物药效学之间的有趣关联。例如,您可以发现一名吸烟者与一名不吸烟者的ANS之间的一些相关性。在频谱测量(如LF、HF或LF/HF)中,可以发现一些显著的差别,利用它们检测出吸烟对ANS的影响。如欲了解更多专题,请查阅文献[2]。</p>
<a name="toc5"></a>
<h2><strong>获取HRV分析起步工具</strong></h2>
<p align="left">本文中提到的HRV分析起步工具包含于<a href="http://zone.ni.com/devzone/cda/tut/p/id/9037" target="_blank">NI Biomedical Startup Kit 2.0</a>。请<span style="text-decoration: underline;">下载</span>该生物医学起步工具以使用HRV分析起步工具。</p>
<p align="left">请致信<a href="mailto:dsp.nish@ni.com">dsp.nish@ni.com</a>,以获取关于本文的更多信息。</p>
<a name="toc6"></a>
<h2><strong>参考文献</strong></h2>
<p>[1] 欧洲心脏病学会与北美起搏与电生理学会联合任务组,心率变异——测量、生理解释与临床使用标准,循环,1996(93),5,页:1043—1065。</p>
<p>[2] Rajendra Acharya U、Paul Joseph K、Kannathal N、Lim Choo Min和Jasjit Suri S,心率变异,<a id="OLE_LINK66" name="OLE_LINK66"></a> <a id="OLE_LINK65" name="OLE_LINK65">心脏信号处理进展,</a>Springer-Verlag Berlin Heidelberg,2007,页码:121-165。</p>
<p>[3] U. Rajendra Acharya、N. Kannathal、Lim Choo Min和Jasjit S. Suri,利用心率信号进行心脏健康的线性、非线性和小波分析,心脏信号处理进展,Springer-Verlag Berlin Heidelberg,2007,页码:355-375。</p>
<p>[4] MICHAEL J. LEWIS,科学学士、科学硕士、博士、英国皇家特许物理学家和英国物理研究所学术成员,心率变异分析:一种评估心脏自主神经功能的工具,CIN:计算机、信息学和护理,2005(23),6,页码:335-341。</p>
<p>[5] Juha-Pekka Niskanen、Mika P. Tarvainen、Perttu O. Ranta-aho和Pasi A. Karjalainen,面向高级HRV分析的软件,生物医学中的计算机方法与程序,2004,76,页码:73-81。</p>
<p>[6] Peng CK、Havlin S和Stanley HE等,非稳定心跳时间序列中的标度指数与交叉现象的量化表示,混沌学,1995,5(1),页码:82-87。</p>
<p>[7] D. Singh、K. Vinod与S. C. Saxena,用于心率变异频谱分析的RR间隔时间序列的采样频率,医学工程与技术期刊,2004(28),6,页码:263–272。</p>
<p>[8] H G Adelmann,采用时频联合分析方法的临床药理学中的心率变异分析,心脏病学计算机技术,1999,26,页码:643-645。</p>
<p> </p>
<div class="noindex">
<div class="docdetails">
<div class="ratings">2 评级 | <strong>3.00</strong> out of 5</div>
<!-- Read in other language section --><strong>阅读语言 </strong><select><option selected value="http://zone.ni.com/devzone/cda/epd/p/id/6275">Chinese (Simplified)</option>
<option value="http://zone.ni.com/devzone/cda/epd/p/id/5832">English</option></select> | <a><img class="bottomalign" src="http://baitai.iteye.com/images/global/us/icons/print.gif" border="0" alt="" width="18" height="13"> 打印</a> </div>
</div>
<div class="noindex">
<div class="legalnotice">
<p><strong>法律条款</strong><br>本范例程序(简称“程序”)由National Instruments公司(简称“NI”)应用工程师开发。 尽管该程序的技术支持由NI提供,但是该程序并非完全通过测试和验证。NI不以任何方式保证其质量,也不保证相关产品或驱动程序的新版本出现时继续为其提供技术支持。 本程序仅以其“现状”向用户提供,程序没有任何担保。程序使用受ni.com网站上《使用条款》的约束。 (<a href="http://ni.com/legal/termsofuse/unitedstates/us/">http://ni.com/legal/termsofuse/unitedstates/us/</a>)</p>
</div>
</div>
</div>
</div>
<div class="noindex">
<div class="clearfloat"> </div>
</div>
</div>
<a id="nibug_2010" href="http://www.ni.com/niweek/?metc=mtxrhy"></a>
</div>
<div class="clearfloat"> </div>
<div id="footer" class="noindex">
<a rel="nofollow" href="http://www.ni.com/profile/zhs/">我的个人信息</a> | <a rel="nofollow" href="http://www.ni.com/rss/">RSS</a> | <a rel="nofollow" href="http://www.ni.com/legal/privacy/unitedstates/us/">隐私权说明</a> | <a rel="nofollow" href="http://www.ni.com/legal/">法律信息</a> | <a rel="nofollow" href="http://sine.ni.com/apps/utf8/nicc.call_me">联系NI</a> <span class="copyright">© 2011 National Instruments Corporation. All rights reserved.</span> |
<script type="text/javascript"></script><a href="mailto:?subject=Link%20to%20National%20Instruments%20Web%20Site&body=http://zone.ni.com/devzone/cda/epd/p/id/6275"><img src="http://baitai.iteye.com/images/global/us/email.gif" border="0" alt="转发该页"> 转发该页</a> </div>
<p>
<script src="/javascript/analysis/us/tag-footer.js" type="text/javascript"></script></p>
<div id="bodycontainer" class="document">
<div id="resources">
<div class="noindex">
<div id="flakes">
<div class="fnavbox">
<ul></ul>
</div>
</div>
</div>
</div>
<div id="primarycontentcontainer">
<div class="noindex">
<h1>利用LabVIEW进行心率变异分析</h1>
<div class="docdetails">
<div class="ratings">2 评级 | <strong>3.00</strong> out of 5</div>
<!-- Read in other language section --><strong>阅读语言 </strong><select><option selected value="http://zone.ni.com/devzone/cda/epd/p/id/6275">Chinese (Simplified)</option>
<option value="http://zone.ni.com/devzone/cda/epd/p/id/5832">English</option></select> | <a><img class="bottomalign" src="http://baitai.iteye.com/images/global/us/icons/print.gif" border="0" alt="" width="18" height="13"> 打印</a> </div>
</div>
<div id="pagearea">
<div class="noindex">
<h3>概览</h3>
</div>
<p>心率变异(HRV)分析是一项在研究和临床中广泛使用的重要应用。LabVIEW包含多个信号处理VI,使用户可以非常方便地创建自定义的HRV分析应用。本文介绍了利用LabVIEW进行HRV分析的应用与实践。</p>
<div class="noindex">
<h3>目录</h3>
<ol>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc0">引言</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc1">什么是HRV?</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc2">采集RR间隔</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc3">HRV分析方法</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc4">应用</a> </li>
<li>
<a href="http://baitai.iteye.com/admin/blogs/#toc5">获取HRV分析起步工具</a> </li>
<li><a href="http://baitai.iteye.com/admin/blogs/#toc6">参考文献</a></li>
</ol>
<h3>下载</h3>
</div>
<a name="toc0"></a>
<h2><strong>引言</strong></h2>
<p align="left"><strong>心电图(ECG)</strong>是一种由心电图仪产生的图形,提供了有关个人心脏健康的信息。ECG常在临床应用中帮助医生诊断心跳过速等心脏疾病。除了直接分析ECG信号,研究人员和医生们还可以从ECG信号中提取其他非直接的测量信息,最常见的测量之一便是<strong>心率变异(HRV)</strong>。本文描述了如何利用NI LabVIEW图形化开发环境提取并分析HRV测量结果,并同时介绍了LabVIEW HRV分析起步工具—— 一款基于LabVIEW的面向HRV分析的程序。本文所展示的HRV分析结果来自于该HRV分析起步工具。</p>
<a name="toc1"></a>
<h2><strong>什么是HRV?</strong></h2>
<p align="left">一个正常的单周期ECG信号由若干个波组成,如图1所示。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57496.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图1</strong>. ECG信号的R波峰与RR间隔</p>
<p align="left">具有最高峰值的波被称为R波。RR间隔是指两个相邻R波之间的时间间隔。波峰较小的波分别为P波、T波和U波。</p>
<p align="left">心率变异(HRV)测量可以分析出这些RR间隔(它们显示了连续心跳之间的变化)是如何随时间而变化的。例如,图2展示了一个典型的成人的HRV。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57536.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图2.</strong> 一个成人的HRV</p>
<p align="left">图2中的HRV分析显示了该RR间隔范围从0.6秒变化至1秒,其平均值为0.8秒。</p>
<p align="left">HRV分析提供了<strong>自主神经系统(ANS)</strong>的一个定量标记,因为HRV的调节机制源自交感和副交感神经系统[5]。迄今为止,许多研究人员和工程师们都致力于发掘HRV分析中对临床应用有价值的信息。在过去的十年中,已经发表了超过2000篇关于HRV的文章[1]。这些文章考察了HRV与血压、心肌梗塞、神经系统、心律失常、糖尿病、呼吸、肾衰竭、性别、年龄、疲劳程度、药物、吸烟、酒精等因素的关系[2]。目前,多个医学学科都已经在对HRV进行研究,相信HRV分析在未来的患者图形与数据分析中的应用将越来越多。</p>
<a name="toc2"></a>
<h2><strong>采集RR间隔</strong></h2>
<p align="left">为了进行HRV分析,首先必须获取RR间隔。图3展示了<a id="OLE_LINK2" name="OLE_LINK2"></a> <a id="OLE_LINK1" name="OLE_LINK1">RR间隔的采集过程</a>。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image22756.jpg" alt="" width="513" height="88"></p>
<p align="center"><strong>图3.</strong> RR间隔的采集过程</p>
<a name="h30"></a>
<h3><strong>采集原始ECG信号</strong></h3>
<p align="left">ECG信号可通过各种心电图仪采集获得。这些仪器的带宽并不要求很高。可以利用NI多通道DAQ设备,如NI PCI-6221,从ECG记录设备的输出端采集原始ECG信号。典型的采样率设置为125 Hz或250 Hz。所采集的ECG信号可以以NI TDMS文件类型存储以便进行离线分析。一些在线数据库(如<a href="http://www.physionet.org/physiobank/database/#ecg">MIT-BIH数据库database</a>)也提供许多典型的ECG信号。该HRV分析起步工具包含一个导入<a href="http://www.physionet.org/physiobank/database/#ecg">MIT-BIH</a> ECG<a href="http://www.physionet.org/physiobank/database/#ecg">数据</a>的向导。</p>
<a name="h31"></a>
<h3><strong>从ECG信号中提取RR间隔</strong></h3>
<p align="left">从原始的ECG信号中可以提取到RR间隔。该提取过程通常包括预处理步骤和波峰检测步骤。如果原始ECG信号具有噪声和显著的基线漂移趋势,就有必要对其进行预处理。然后,您可以通过阈值设置或利用基于小波的波峰检测方法,检测R波峰,从而计算RR间隔。</p>
<p align="left">请参阅<a href="http://zone.ni.com/devzone/cda/tut/p/id/10041">NI开发者园地</a>,获取关于ECG信号处理和特性提取的更多信息。</p>
<a name="toc3"></a>
<h2><strong>HRV分析方法</strong></h2>
<p align="left">HRV分析有不同的方法[1][2]。方法之一便是时域分析,即利用暂态的RR间隔信号提取一些特别的度量。另一个方法是频谱分析,该方法以某个特定的速率对RR间隔进行插值处理,并将该间隔信号变换到频域。针对这两种方法,目前有一些标准可供使用[1]。还有一些其他方法,如非线性分析和时频分析,也有相关讨论和研究。</p>
<p align="left">图4展示了HRV分析的整个过程。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image22758.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图4.</strong> HRV分析的整个过程。</p>
<p align="left">LabVIEW具有多个可用于HRV分析的工具包或函数。如图5所示的HRV分析起步工具,展现了利用LabVIEW进行HRV分析的简便性。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57500.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图5.</strong> HRV分析起步工具的前面板</p>
<a name="h32"></a>
<h3><strong>预处理</strong></h3>
<p align="left">RR间隔信号可能存在一些基线抖动,这会导致不精确的暂态分析和频率分析。您可以选择是否在分析前对RR间隔信号进行去漂移处理,如图6所示。对RR间隔信号进行去漂移处理的方法,与<a href="http://zone.ni.com/devzone/cda/tut/p/id/10041">使用LabVIEW进行ECG信号处理</a>文档中所提及的方法相似。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57502.jpg" alt="" width="364" height="203"></p>
<p align="center"><strong>图6.</strong> 预处理配置对话框</p>
<p align="left">图7展示了去漂移处理的过程。图中的红线表示该基线漂移。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57504.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图7.</strong> RR间隔信号的预处理</p>
<a name="h33"></a>
<h3><strong>分析RR间隔信号</strong></h3>
<p align="left">下面描述了用于分析RR间隔信号的不同方法。HRV分析起步工具支持所有这些方法。</p>
<a name="h40"></a>
<h4><strong>时域分析</strong></h4>
<p align="left">对于时间序列分析,通常使用时域量度。许多量度可以从最初的RR间隔信号中提取,以展现ANS中的变化。表1列出了利用HRV分析起步工具提取的量度。</p>
<p align="center"><strong>表1.</strong> HRV的时域量度[1][5]</p>
<table border="1" cellspacing="0" cellpadding="0" align="center"><tbody>
<tr>
<td width="139" valign="top">
<p align="center"><strong>变量</strong></p>
</td>
<td width="54" valign="top">
<p align="center"><strong>单位</strong></p>
</td>
<td width="397" valign="top">
<p align="center"><strong>描述</strong></p>
</td>
</tr>
<tr>
<td colspan="3" width="590" valign="top">
<p align="center">统计学量度</p>
</td>
</tr>
<tr>
<td width="139" valign="top">R<a id="OLE_LINK10" name="OLE_LINK10"></a> <a id="OLE_LINK9" name="OLE_LINK9">R均值与标准差</a>
</td>
<td width="54" valign="top">s</td>
<td width="397" valign="top">所有RR间隔的均值与标准差</td>
</tr>
<tr>
<td width="139" valign="top">HR均值与标准差</td>
<td width="54" valign="top">1/min</td>
<td width="397" valign="top">所有心率的均值与标准差</td>
</tr>
<tr>
<td width="139" valign="top">RMSSD</td>
<td width="54" valign="top">ms</td>
<td width="397" valign="top">各个相邻RR间隔的差值的均方根</td>
</tr>
<tr>
<td width="139" valign="top">NN50计数</td>
<td width="54" valign="top"> </td>
<td width="397" valign="top">在所有测量值中,差值超过50ms的相邻RR间隔对的数目</td>
</tr>
<tr>
<td width="139" valign="top">pNN50</td>
<td width="54" valign="top">%</td>
<td width="397" valign="top">NN50计数除以所有RR间隔的总数目</td>
</tr>
<tr>
<td colspan="3" width="590" valign="top">
<p align="center">几何学量度</p>
</td>
</tr>
<tr>
<td width="139" valign="top">HRV三角指数</td>
<td width="54" valign="top"> </td>
<td width="397" valign="top">所有RR间隔的总数除以所有RR间隔的直方图的高度</td>
</tr>
<tr>
<td width="139" valign="top">TINN</td>
<td width="54" valign="top">ms</td>
<td width="397" valign="top">对于间距为7.8125 ms(1/128 s)的离散标度所测得的所有RR间隔,对其直方图中最高峰的最小均方误差三角插值的基线宽度。</td>
</tr>
</tbody></table>
<p align="left">LabVIEW中包含概率与统计VI,您可以利用这些VI方便地获取统计学度量,如图8所示。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57538.jpg" alt="" width="430" height="341"></p>
<p align="center"><strong>图8.</strong> LabVIEW中的概率与统计VI</p>
<p align="left">图9展示了一个成人的RR间隔信号样本。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57512.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图9.</strong> 利用时域分析方法获得的RR间隔信号采样</p>
<p align="left">图10展示了利用HRV分析起步工具中的时域分析方法得到的这些信号的统计学结果和几何学结果。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57514.jpg" alt="" width="364" height="211"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57516.jpg" alt="" width="364" height="265"></p>
<p align="center"><strong>图10.</strong> 利用时域分析方法的HRV分析</p>
<a name="h41"></a>
<h4><strong>频域分析</strong></h4>
<p align="left">HRV分析的频域分析方法从RR间隔信号中提取频域参数,如峰值频率和带内功率。ANS的两个组成部分——交感神经系统和副交感神经系统——会决定升高或降低心率,并影响RR间隔频谱中的不同频带。因而,可以利用频域分析监测ANS的状态。</p>
<p align="left">表2展示了HRV的常见频域度量,其中包括峰值频率以及超低频(VLF)、低频(LF)和高频(HF)等频带的带内功率。</p>
<p> </p>
<p align="center"><strong>表2.</strong> HRV的频域量度[1][5]</p>
<table border="1" cellspacing="0" cellpadding="0" align="center"><tbody>
<tr>
<td width="139" valign="top">
<p align="center"><strong>变量</strong></p>
</td>
<td width="54" valign="top">
<p align="center"><strong>单位</strong></p>
</td>
<td width="397" valign="top">
<p align="center"><strong>描述</strong></p>
</td>
</tr>
<tr>
<td width="139" valign="top">峰值频率</td>
<td width="54" valign="top">Hz</td>
<td width="397" valign="top">针对VLF、LF和HF频带的功率谱密度(PSD)估计的峰值频率</td>
</tr>
<tr>
<td width="139" valign="top">VLF </td>
<td width="54" valign="top">ms<sup>2</sup>
</td>
<td width="397" valign="top">0–0.04 Hz<a id="OLE_LINK20" name="OLE_LINK20"></a> <a id="OLE_LINK19" name="OLE_LINK19">频段的功率</a>
</td>
</tr>
<tr>
<td width="139" valign="top">LF </td>
<td width="54" valign="top">ms<sup>2</sup>
</td>
<td width="397" valign="top">0.04–0.15 Hz频段的功率</td>
</tr>
<tr>
<td width="139" valign="top">HF</td>
<td width="54" valign="top">ms<sup>2</sup>
</td>
<td width="397" valign="top">0.15–0.4 Hz频段的功率</td>
</tr>
<tr>
<td width="139" valign="top">LF归一化</td>
<td width="54" valign="top">无量纲</td>
<td width="397" valign="top">
<a id="OLE_LINK30" name="OLE_LINK30"> </a><a id="OLE_LINK29" name="OLE_LINK29">LF功率(采用归一化单位):LF/(总功率–VLF)*100</a>
</td>
</tr>
<tr>
<td width="139" valign="top">HF 归一化</td>
<td width="54" valign="top">无量纲</td>
<td width="397" valign="top">HF功率(采用归一化单位):HF/(总功率–VLF)*100。</td>
</tr>
<tr>
<td width="139" valign="top">LF/HF比</td>
<td width="54" valign="top"> </td>
<td width="397" valign="top">LF/HF(均采用ms<sup>2</sup>作为单位)</td>
</tr>
</tbody></table>
<p align="left">频域方法通常涉及以下三个步骤:</p>
<ol>
<li>对RR间隔信号进行重采样。 </li>
<li>估算该RR间隔信号的功率谱密度(PSD)。 </li>
<li>根据该PSD计算频域参数。 </li>
</ol>
<p align="left">从ECG信号中抽取的RR间隔信号通常不是以一个规则间隔采样所得。因此,该RR间隔信号并不是延时间轴等距离间隔。可以利用位于<strong>数学>>内插与外推</strong>函数选板中的一维插值VI,将最初的RR信号转换为均匀采样的时间序列。</p>
<p align="left">对于PSD估计,LabVIEW包含可以执行非参数估计方法或参数估计方法的VI。非参数方法包括基于FFT的PSD,如周期图法和<a id="OLE_LINK34" name="OLE_LINK34"></a> <a id="OLE_LINK33" name="OLE_LINK33">沃尔什法</a>。参数方法包括基于模型的PSD,如自回归(<a id="OLE_LINK36" name="OLE_LINK36"></a> <a id="OLE_LINK35" name="OLE_LINK35">AR)方法</a>和自回归滑动平均(ARMA)方法。<a id="OLE_LINK38" name="OLE_LINK38"></a> <a id="OLE_LINK37" name="OLE_LINK37">HRV分析起步工具</a>包含沃尔什法和AR频谱方法。LabVIEW高级信号处理工具包中则包含有其他方法。</p>
<p align="left">图11展示了利用HRV分析起步工具中的频域分析方法得到的一个RR间隔信号分析结果。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57518.jpg" alt="" width="436" height="257"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57520.jpg" alt="" width="436" height="257"></p>
<p align="center"><strong>图11.</strong> 利用频域分析方法的HRV分析</p>
<a name="h42"></a>
<h4><strong>联合时频域分析</strong></h4>
<p align="left">频域分析方法适用于随时间变化不大的RR间隔信号分析。然而,在长期HRV分析(如24小时RR间隔与药物药效学研究)等应用中,RR间隔信号会随时间显著变化。在这类情况下,您可以使用联合时-频分析(JTFA)方法分析稳态和时-频行为。</p>
<p align="left">联合时频域方法(JTFA)没有标准的度量参数。因此可以利用短时傅立叶变换(STFT)、Gabor展开或连续小波VI得到时频谱图并进行定性分析。文献[8]中提到了一个利用JTFA方法的药物药效学应用案例。</p>
<p align="left">高级信号处理工具包支持JTFA方法。图12显示了利用该工具包所含的LabVIEW代码,该代码可以计算得到RR间隔信号的STFT谱图。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57522.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图12.</strong> 获取RR间隔信号的STFT谱图</p>
<p align="left">图13展示了利用HRV分析起步工具中的JTFA方法获得的联合时频分析结果。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57524.jpg" alt="" width="278" height="182"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57526.jpg" alt="" width="278" height="182"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57528.jpg" alt="" width="278" height="183"></p>
<p align="center"><strong>图13.</strong> 利用JTFA方法的HRV分析</p>
<a name="h43"></a>
<h4><strong>非线性分析</strong></h4>
<p align="left">RR间隔信号是非线性的,因为它们来源于血流动力学、电生理学和体液参数、以及自主神经调节与中枢神经调节之间的复杂交互关系[1]。关于非线性分析的研究基本处于空白,所以没有任何相关标准。然而,这是一个充满希望的研究领域,我们需要在该领域开发出一种工具,从而帮助我们更好地理解人体系统的复杂性。</p>
<p align="left">HRV分析起步工具支持<a id="OLE_LINK40" name="OLE_LINK40"></a> <a id="OLE_LINK39" name="OLE_LINK39">庞加莱图</a>和<a id="OLE_LINK50" name="OLE_LINK50"></a> <a id="OLE_LINK49" name="OLE_LINK49">去漂移波动分析</a>(DFA)方法。下面将分别介绍这些方法。</p>
<p align="left"><strong>庞加莱图</strong></p>
<p align="left">庞加莱图是一种XY图,其X轴和Y轴表示两个连续的RR间隔。该图所用的模型是一个椭圆,该椭圆的对称线与X轴成45◦,如图14所示。</p>
<p align="center"><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57530.jpg" alt="" width="390" height="247"></p>
<p align="center"><strong>图14.</strong> RR间隔的庞加莱图</p>
<p align="left">图14显示了下列数值:</p>
<ul>
<li>
<strong>Sd1</strong>—各点到对称线距离的标准差。Sd1刻画出了短期变化,这主要是由呼吸性窦性心律不齐(RSA)引起的。 </li>
<li>
<strong>Sd2</strong>—各点在沿着对称直线方向上的标准偏差。Sd2刻画了长期变化[2]。 </li>
</ul>
<p align="left">LabVIEW包含概率与统计VI和拟合VI,以及一个前面板XY图控件,可以利用该控件创建并显示一幅庞加莱图。图15显示了<strong>拟合</strong>VI。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57532.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图15.</strong> LabVIEW所包含的拟合VI库</p>
<p align="left"><strong>去漂移波动分析</strong></p>
<p align="left">去漂移波动分析(DFA)方法定量表示了RR间隔信号的分形尺度属性。该方法通过剔除不稳定的人为趋势,避免了错误的相关性。因此,相比其他方法,DFA方法能更好展示HRV长期相关性的变化情况。</p>
<p align="left">DFA方法利用线性拟合对DFA图进行分析。拟合直线的斜率α表示了HRV的相关性。如果α=0.5,那么该信号为白噪声,没有任何相关性。在一些情况下,DFA图不是严格地呈线性,而是由相交于某点的两条直线组成。该情形表明存在两个因子:一个短期内的小尺度指数和一个长期时间内的大尺度指数。</p>
<p align="left">图16显示了利用HRV分析起步工具中的DFA方法获得的一个RR间隔信号分析结果。</p>
<p align="center"><a><img src="http://baitai.iteye.com/cms/images/devzone/epd/image57534.jpg" border="0" alt="" width="525"></a><br><a>[+] 放大图片</a></p>
<p align="center"><strong>图16.</strong> RR间隔信号的DFA图</p>
<a name="toc4"></a>
<h2><strong>应用</strong></h2>
<p align="left">HRV分析结果可以应用于临床和研究。对于统计与频域方法,HRV分析结果[1]提出了一些临床应用HRV分析以帮助诊断各种心脏疾病的标准。利用HRV分析起步工具,研究人员能够发现HRV与疾病以及药物药效学之间的有趣关联。例如,您可以发现一名吸烟者与一名不吸烟者的ANS之间的一些相关性。在频谱测量(如LF、HF或LF/HF)中,可以发现一些显著的差别,利用它们检测出吸烟对ANS的影响。如欲了解更多专题,请查阅文献[2]。</p>
<a name="toc5"></a>
<h2><strong>获取HRV分析起步工具</strong></h2>
<p align="left">本文中提到的HRV分析起步工具包含于<a href="http://zone.ni.com/devzone/cda/tut/p/id/9037" target="_blank">NI Biomedical Startup Kit 2.0</a>。请<span style="text-decoration: underline;">下载</span>该生物医学起步工具以使用HRV分析起步工具。</p>
<p align="left">请致信<a href="mailto:dsp.nish@ni.com">dsp.nish@ni.com</a>,以获取关于本文的更多信息。</p>
<a name="toc6"></a>
<h2><strong>参考文献</strong></h2>
<p>[1] 欧洲心脏病学会与北美起搏与电生理学会联合任务组,心率变异——测量、生理解释与临床使用标准,循环,1996(93),5,页:1043—1065。</p>
<p>[2] Rajendra Acharya U、Paul Joseph K、Kannathal N、Lim Choo Min和Jasjit Suri S,心率变异,<a id="OLE_LINK66" name="OLE_LINK66"></a> <a id="OLE_LINK65" name="OLE_LINK65">心脏信号处理进展,</a>Springer-Verlag Berlin Heidelberg,2007,页码:121-165。</p>
<p>[3] U. Rajendra Acharya、N. Kannathal、Lim Choo Min和Jasjit S. Suri,利用心率信号进行心脏健康的线性、非线性和小波分析,心脏信号处理进展,Springer-Verlag Berlin Heidelberg,2007,页码:355-375。</p>
<p>[4] MICHAEL J. LEWIS,科学学士、科学硕士、博士、英国皇家特许物理学家和英国物理研究所学术成员,心率变异分析:一种评估心脏自主神经功能的工具,CIN:计算机、信息学和护理,2005(23),6,页码:335-341。</p>
<p>[5] Juha-Pekka Niskanen、Mika P. Tarvainen、Perttu O. Ranta-aho和Pasi A. Karjalainen,面向高级HRV分析的软件,生物医学中的计算机方法与程序,2004,76,页码:73-81。</p>
<p>[6] Peng CK、Havlin S和Stanley HE等,非稳定心跳时间序列中的标度指数与交叉现象的量化表示,混沌学,1995,5(1),页码:82-87。</p>
<p>[7] D. Singh、K. Vinod与S. C. Saxena,用于心率变异频谱分析的RR间隔时间序列的采样频率,医学工程与技术期刊,2004(28),6,页码:263–272。</p>
<p>[8] H G Adelmann,采用时频联合分析方法的临床药理学中的心率变异分析,心脏病学计算机技术,1999,26,页码:643-645。</p>
<p> </p>
<div class="noindex">
<div class="docdetails">
<div class="ratings">2 评级 | <strong>3.00</strong> out of 5</div>
<!-- Read in other language section --><strong>阅读语言 </strong><select><option selected value="http://zone.ni.com/devzone/cda/epd/p/id/6275">Chinese (Simplified)</option>
<option value="http://zone.ni.com/devzone/cda/epd/p/id/5832">English</option></select> | <a><img class="bottomalign" src="http://baitai.iteye.com/images/global/us/icons/print.gif" border="0" alt="" width="18" height="13"> 打印</a> </div>
</div>
<div class="noindex">
<div class="legalnotice">
<p><strong>法律条款</strong><br>本范例程序(简称“程序”)由National Instruments公司(简称“NI”)应用工程师开发。 尽管该程序的技术支持由NI提供,但是该程序并非完全通过测试和验证。NI不以任何方式保证其质量,也不保证相关产品或驱动程序的新版本出现时继续为其提供技术支持。 本程序仅以其“现状”向用户提供,程序没有任何担保。程序使用受ni.com网站上《使用条款》的约束。 (<a href="http://ni.com/legal/termsofuse/unitedstates/us/">http://ni.com/legal/termsofuse/unitedstates/us/</a>)</p>
</div>
</div>
</div>
</div>
<div class="noindex">
<div class="clearfloat"> </div>
</div>
</div>
<a id="nibug_2010" href="http://www.ni.com/niweek/?metc=mtxrhy"></a>
</div>
<div class="clearfloat"> </div>
<div id="footer" class="noindex">
<a rel="nofollow" href="http://www.ni.com/profile/zhs/">我的个人信息</a> | <a rel="nofollow" href="http://www.ni.com/rss/">RSS</a> | <a rel="nofollow" href="http://www.ni.com/legal/privacy/unitedstates/us/">隐私权说明</a> | <a rel="nofollow" href="http://www.ni.com/legal/">法律信息</a> | <a rel="nofollow" href="http://sine.ni.com/apps/utf8/nicc.call_me">联系NI</a> <span class="copyright">© 2011 National Instruments Corporation. All rights reserved.</span> |
<script type="text/javascript"></script><a href="mailto:?subject=Link%20to%20National%20Instruments%20Web%20Site&body=http://zone.ni.com/devzone/cda/epd/p/id/6275"><img src="http://baitai.iteye.com/images/global/us/email.gif" border="0" alt="转发该页"> 转发该页</a> </div>
<p>
<script src="/javascript/analysis/us/tag-footer.js" type="text/javascript"></script></p>
相关推荐
MATLAB源代码MATLAB源码大集合220MB上千个源码文件: 2013全国大学生数学建模B题matlab代码.rar 802.11a OFDM MATLAB仿真代码!(推荐).rar MATLAB DCT水印源程序代码.rar MATLAB GUI 实例.rar matlab GUI实例.rar ...
谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换器matlab仿真文件谐振变换...
高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码高斯伪谱法matlab代码...
Matlab地铁牵引供电系统建模与仿真 Matlab地铁牵引供电系统建模与仿真 Matlab地铁牵引供电系统建模与仿真 Matlab地铁牵引供电系统建模与仿真 Matlab地铁牵引供电系统建模与仿真 Matlab地铁牵引供电系统建模与仿真 ...
在开发过程中,有时我们需要将不同工具的优势结合起来,例如在C#应用中利用Matlab的强大计算能力及图形绘制功能。本文将详细介绍如何将Matlab的图形窗口嵌入到C#窗体中,以便在C#应用中实时展示Matlab生成的图形。 ...
matlab在关联编译器的时候,会出现找不到编译器的问题,原因是matlab要找windows sdk的mt.exe,找不到报错的.因为vs2012在安装时,windows sdk的mt.exe位置变了,而matlab还是按照老的路径找,肯定找不到(原来在windows ...
【Matlab精简版】是针对那些只需要使用Matlab基础功能的用户设计的,它相比标准版具有更小的体积,降低了对系统资源的占用,是轻量级的计算工具。Matlab是一款强大的数学计算软件,广泛应用于工程计算、数据分析、...
标题“matlab xgboost安装调试”涉及到的是在MATLAB环境中安装和配置XGBoost库的过程,以及可能遇到的问题和调试方法。XGBoost是一种高效、灵活且可扩展的梯度提升框架,常用于机器学习任务,特别是分类和回归问题。...
基于Matlab的卫星轨道仿真(源代码) 基于Matlab的卫星轨道仿真(源代码) 基于Matlab的卫星轨道仿真(源代码) 基于Matlab的卫星轨道仿真(源代码) 基于Matlab的卫星轨道仿真(源代码) 基于Matlab的卫星轨道...
代码 马尔科夫链蒙特卡洛模拟的matlab源代码代码 马尔科夫链蒙特卡洛模拟的matlab源代码代码 马尔科夫链蒙特卡洛模拟的matlab源代码代码 马尔科夫链蒙特卡洛模拟的matlab源代码代码 马尔科夫链蒙特卡洛模拟的matlab...
基于MATLAB神经网络的股票价格预测基于MATLAB神经网络的股票价格预测基于MATLAB神经网络的股票价格预测基于MATLAB神经网络的股票价格预测基于MATLAB神经网络的股票价格预测基于MATLAB神经网络的股票价格预测基于...
在MATLAB环境中,可以利用特定的驱动程序与CAN接口设备进行通信,例如周立功的USBCAN,这是一种通过USB接口连接到计算机的CAN收发器。MATLAB Guide程序则是用于构建图形用户界面(GUI)的工具,使得用户能够更加直观地...
在IT领域,C++与Matlab的交互是一个常见的需求,特别是在数据分析、科学计算以及图形可视化等方面。本主题将深入探讨如何使用C++调用Matlab的函数进行绘图操作。以下是一个详细的步骤和相关知识点: 首先,为了在...
此时,通过Matlab调用Cplex,可以利用Cplex的强大性能,结合Matlab的易用性,实现更高效的问题求解。 首先,要实现Matlab调用Cplex,你需要确保已经安装了Cplex和Matlab,并且在系统路径中包含了Cplex的库文件。...
This package is used for calling r language from matlab Most of the files already exist on the Internet But I wrote five to complement this package 1 In R enter: >install packages "rscproxy...
《cplex12.9与MATLAB 2021a的集成与应用》 在IT领域,优化问题的解决是关键任务之一,IBM的Cplex是一款强大的优化求解器,广泛应用于数学规划、线性规划、整数规划等领域。Cplex12.9版本与MATLAB 2021a的结合,为...
MATLAB是一款强大的数学计算软件,广泛应用于工程计算、科学建模和数据分析等领域。2020a版本的官方PDF文档提供了全面的指南,帮助用户掌握MATLAB的各种功能和使用技巧。以下是对这些知识点的详细说明: 1. **...
《Matlab教程(从新手到骨灰级玩家).pdf》是一个全面且深入的Matlab学习资源,适合不同水平的Matlab使用者。Matlab是一种强大的计算环境,尤其在矩阵运算、图像处理、信号处理和机器学习等领域有着广泛的应用。本教程...
【标题】"matlab生成菲涅尔透镜"涉及到的核心技术是使用MATLAB进行图形绘制,特别是光学领域的菲涅尔透镜模拟。MATLAB是一种强大的数值计算和数据分析软件,同时也支持高级图形绘制,使得用户可以方便地创建各种复杂...
matlab在关联编译器的时候 会出现找不到编译器的问题 原因是matlab要找windows sdk的mt exe 本文件修改方法也是根据wiills的那篇“matlab找不到编译器的问题”(matlab无法找到vs2012)的修改方法进行的 请参考:...