反演计算生物组织电特性参数.pdfVIP

反演计算生物组织电特性参数 张敏，黄卡玛，汤元九，’贺明寿 (四川大学屯子信息学院 610064) 摘要：本文采用了遗传算法求解生物组织的电特性参数。研究结果表明：利用GA可有 效地反演出生物组织的分层厚度信息和生物组织的某些德拜参数．研究中为模拟实际测 量结果．对正问题计算结果加入误差后再进行反演．从而找出了误差对反演结果的影响。 关键词：遗传算法德拜方程微波测量(FD)叶D 一、 引言 现有的生物组织介电常数和电导率的测量方法：卢，常采用开口同轴线测量生物组织 的反射系数，然后进行计算，得出其电特性参数。传统的计算方法中，有采用查反射系 数与介屯常数、电导率的关系图得到介电常数和电导率，或者根据天线理论，计算阻抗 变换，解出介电常数和电导率【l，2j。遗传算法是新近发展迅速扮全局最优搜索算法， 它根据自然界进化的原理，采用简单的编码技术和繁殖机制表达复杂的现象，从而解决 非常困难的问题，在许多领域得到广泛的应用．【3，4】 本文提出采用遗传算法搜索最优解的方法，反演计算生物组织的电特性参数．重点 讨论了目标函数的确定，适应度函数的选取对遗传结果的影响，而且针对实际情况，模 拟分析了测量误差对计算结果的影响，为实际应用遗传算法作出指导。 二、 生物组织屯特性反演计算模型 采用开口同轴线作测量，生物组织的模型分为三层：皮肤、脂肪和肌肉。为了得到 更多频率点的介电信息，采用脉冲作为测量的电磁波。由于生物组织为色散媒质，介电 常数与频率的关系用德拜方程表示。气表示频率趋于无穷对的相对介电常数．占。表示 准静态的相对介电常数，f表示驰豫时间。 占，=靠+二L=旦 (1) “ ’ 1+_，国f 在进行数值模拟时，必须采用(FD)叶D方法计算[5．6]。入射脉冲采用高斯脉冲 的时间导数波形．由于三层组织的介电常数和电导率不同。因此入射脉冲经过每两层交 界处产生反射。将入射脉冲与反射脉冲进行离散傅立叶变换，得到频谱图，发现两者的 频谱曲线相似，但主要频点的幅值不同。 1．0 O．8 O．6 04 02 OO O 1 O00E+0005．OO暑·01 1．00E·01 时问 蛹搴lHz) 田】入射脉冲的波形 圈2^射辣冲和反射脉冲的频谱 根据GA的计算步骤，本文的问题中，将德拜方程的未知参数用一系列二进制串表 达，这就是染色体串。然后将染色体的值代入正问题的计算．得到反射脉冲在某些重要 频点上的幅值，这些幅值是GA中的目标函数，将计算得到的目标函数与测量的这些频 点韵值相比较，根据适应度函数，计算出每一个二进制串的适应度值．根据适应度值选 择遗传算子计算产生新一代的染色体串。然后重复进行这几个步骤．直到满足停止条件。 GA优化的关键是选取恰当的适应度函数，因为它直接关系到搜索解空间的收敛速 度．本文采用如下函数作为适应度函数： fitness：P一。√R R：釜阶)一，r(f)】2 faI 式中c取常数0 1，Nr是所选择的频点数目，f’(i)是测量的各频点幅值，《i)是根据 染色体确定的参数所计算的相应频点幅值。 三、计算结果分析 假设肌肉的德拜参数己知，皮肤和脂肪的德拜参数中，已知占0，占。以