论文部分内容阅读
电阻抗肌动描记术(Electrical Impedance Myography, EIM)是哈佛医学院SewardB.Rutkove教授提出的一种基于肌肉电阻抗测量与建模分析的神经肌肉疾病诊断新方法,具有无辐射、非侵入、操作简单、低成本等优点,在神经肌肉疾病的早期发现与诊断、长期监测、治疗以及药物研制等方面具有潜在的应用前景。
然而,现有的EIM技术无法同时满足快速性与准确性的要求。首先,基于50kHz正弦激励的单频EIM分析方法包含信息较少,无法对肌肉宽频阻抗特征进行准确的建模分析;其次,基于正弦扫频激励的EIM分析方法虽然可以在不同频率提取更多的肌肉阻抗特征,并用于建立对应的等效电路模型,但通常测量时间较长,易受被测介质极化影响,且模型参数评估需要迭代优化,易受迭代初值影响,无法实现快速准确的肌肉特征辨识。因此,如何进一步改进EIM测量技术,优化EIM建模理论,建立快速准确的EIM辨识方法,对其在神经肌肉疾病诊断中的应用推广意义重大。针对上述问题,论文主要研究内容如下:
(1)单频EIM全响应(瞬态响应+稳态响应)瞬时测量及常用EIM整数阶、分数阶模型参数时域分析新方法研究。为实现基于单频激励的整数阶Fricke-Morse阻抗等效模型和基于含偏置单频激励的分数阶Cole阻抗等效模型快速准确辨识,提出(含偏置)单频全响应瞬时测量方法,以及上述模型参数数值解析/时域拟合方法。首先,推导(含偏置)正弦电流激励在上述模型下对应的电压全响应数学表达式。随后,提取该全响应信号在瞬态和稳态过程中包含的时域、频域特征,并通过联立全响应各阶段对应的特征方程求解得到模型参数的解析表达式。最后,采用之前解析方法评估的模型参数作为迭代初值,通过迭代优化使测量的与基于模型计算的电压全响应信号误差平方和最小,从而有效的减小初值选取不当对模型参数迭代优化速度和精度的影响;
(2)宽频(同时包含多个频率分量)EIM稳态响应整周期同步测量及常用EIM整数阶、分数阶模型参数频域分析新方法研究。为实现基于宽频激励的常用EIM阻抗等效模型参数快速准确评估,提出一种基于多频正弦(Multisine)激励的整周期同步测量方法及对应的频域建模方法。首先,分析基于Multisine激励的宽频电阻抗谱整周期同步测量原理,从根本上消除频谱泄露对宽频EIM测量的影响。随后,设计一种稀疏频谱分布的低波峰因数Multisine信号生成优化算法,并以此为基础合成适用于宽频EIM整周期同步测量的近二值Multisine激励信号,从而提高宽频EIM测量的信噪比。最后,通过分析EIM常用阻抗等效模型在Nyquist图上对应曲线的特征,推导一种基于最小二乘矩阵运算的模型特征参数频域快速拟合方法,使测量的与基于模型计算的阻抗谱数据在Nyquist图上对应曲线的误差平方和最小,从而避免使用传统最小二乘迭代运算易陷入局部最小值且耗时的缺点;
(3)EIM测量与建模新方法仿真和实验分析。为验证上述各方法的有效性,分别进行不同类型的仿真实验,并搭建基于新型宽频镜像恒流源的EIM高速测量系统,对标准2R-1C三元件电路和肌肉组织进行单频全响应测量、宽频稳态响应整周期测量,以及对应常用整数阶与分数阶阻抗等效模型参数评估。仿真结果表明在40dB高斯白噪声下,新提出的单频和宽频EIM分析方法与传统扫频方法相比,模型参数评估相对误差均<1%,但单频和宽频测量时间分别为0.295ms和1ms,远小于扫频所需的19ms。此外,通过对比标准生物电阻抗测定仪实验结果表明,对于标准2R-1C电路,采用新仪器单频测量分析对应的模型参数相对误差<5%,略大于采用新仪器宽频测量与标准仪器扫频测量的分析结果,后两者对应模型参数相对误差均<2%。对于肌肉组织,采用Cole模型(平均RMSE<2%)比Fricke-Morse模型(平均RMSE<12%)建模更准确,且对于后者,采用单频分析比宽频或扫频分析提取的模型参数RMSE平均值更低(约减小3%)。
以上研究成果有助于发展和完善EIM分析技术与理论,对简化阻抗测量电路、降低阻抗测量仪器成本、提高数据分析精度与速度等方面具有重要实用价值,并为新型EIM测量分析仪器的研制以及在后续神经肌肉疾病诊断中的应用提供新思路。
然而,现有的EIM技术无法同时满足快速性与准确性的要求。首先,基于50kHz正弦激励的单频EIM分析方法包含信息较少,无法对肌肉宽频阻抗特征进行准确的建模分析;其次,基于正弦扫频激励的EIM分析方法虽然可以在不同频率提取更多的肌肉阻抗特征,并用于建立对应的等效电路模型,但通常测量时间较长,易受被测介质极化影响,且模型参数评估需要迭代优化,易受迭代初值影响,无法实现快速准确的肌肉特征辨识。因此,如何进一步改进EIM测量技术,优化EIM建模理论,建立快速准确的EIM辨识方法,对其在神经肌肉疾病诊断中的应用推广意义重大。针对上述问题,论文主要研究内容如下:
(1)单频EIM全响应(瞬态响应+稳态响应)瞬时测量及常用EIM整数阶、分数阶模型参数时域分析新方法研究。为实现基于单频激励的整数阶Fricke-Morse阻抗等效模型和基于含偏置单频激励的分数阶Cole阻抗等效模型快速准确辨识,提出(含偏置)单频全响应瞬时测量方法,以及上述模型参数数值解析/时域拟合方法。首先,推导(含偏置)正弦电流激励在上述模型下对应的电压全响应数学表达式。随后,提取该全响应信号在瞬态和稳态过程中包含的时域、频域特征,并通过联立全响应各阶段对应的特征方程求解得到模型参数的解析表达式。最后,采用之前解析方法评估的模型参数作为迭代初值,通过迭代优化使测量的与基于模型计算的电压全响应信号误差平方和最小,从而有效的减小初值选取不当对模型参数迭代优化速度和精度的影响;
(2)宽频(同时包含多个频率分量)EIM稳态响应整周期同步测量及常用EIM整数阶、分数阶模型参数频域分析新方法研究。为实现基于宽频激励的常用EIM阻抗等效模型参数快速准确评估,提出一种基于多频正弦(Multisine)激励的整周期同步测量方法及对应的频域建模方法。首先,分析基于Multisine激励的宽频电阻抗谱整周期同步测量原理,从根本上消除频谱泄露对宽频EIM测量的影响。随后,设计一种稀疏频谱分布的低波峰因数Multisine信号生成优化算法,并以此为基础合成适用于宽频EIM整周期同步测量的近二值Multisine激励信号,从而提高宽频EIM测量的信噪比。最后,通过分析EIM常用阻抗等效模型在Nyquist图上对应曲线的特征,推导一种基于最小二乘矩阵运算的模型特征参数频域快速拟合方法,使测量的与基于模型计算的阻抗谱数据在Nyquist图上对应曲线的误差平方和最小,从而避免使用传统最小二乘迭代运算易陷入局部最小值且耗时的缺点;
(3)EIM测量与建模新方法仿真和实验分析。为验证上述各方法的有效性,分别进行不同类型的仿真实验,并搭建基于新型宽频镜像恒流源的EIM高速测量系统,对标准2R-1C三元件电路和肌肉组织进行单频全响应测量、宽频稳态响应整周期测量,以及对应常用整数阶与分数阶阻抗等效模型参数评估。仿真结果表明在40dB高斯白噪声下,新提出的单频和宽频EIM分析方法与传统扫频方法相比,模型参数评估相对误差均<1%,但单频和宽频测量时间分别为0.295ms和1ms,远小于扫频所需的19ms。此外,通过对比标准生物电阻抗测定仪实验结果表明,对于标准2R-1C电路,采用新仪器单频测量分析对应的模型参数相对误差<5%,略大于采用新仪器宽频测量与标准仪器扫频测量的分析结果,后两者对应模型参数相对误差均<2%。对于肌肉组织,采用Cole模型(平均RMSE<2%)比Fricke-Morse模型(平均RMSE<12%)建模更准确,且对于后者,采用单频分析比宽频或扫频分析提取的模型参数RMSE平均值更低(约减小3%)。
以上研究成果有助于发展和完善EIM分析技术与理论,对简化阻抗测量电路、降低阻抗测量仪器成本、提高数据分析精度与速度等方面具有重要实用价值,并为新型EIM测量分析仪器的研制以及在后续神经肌肉疾病诊断中的应用提供新思路。