46 传感器与微系统(Transducer and Microsystem Technologies) 2015年第34卷第6期 DOI:10.13873/J.1000-9787(2015)06--0046--04 全自动酶免分析仪定量分析系统研究 常海涛 一,祝连庆 ,王中宇 ,周哲海 ,郭阳宽 (1.北京航空航天大学仪器科学与光电工程学院,北京100191; 2.北京信息科技大学光电测试技术北京市重点实验室,北京100192) 摘要:为了提高酶免分析仪检测精度,建立了基于改进遗传算法的多组分定量分析校正模型。通过自 适应寻优空间技术和小生境技术,对遗传算法进行改进。采用全息凹面光栅作为分光元件,搭建扫描式紫 外可见分光光度实验系统,对食品工业中最常见的苋菜红和胭脂红二组分混合水溶液进行定量分析,并与 最小二乘回归方法进行了比较。实验结果表明:改进遗传算法预测精度与最小二乘法相近,对苋菜红和胭 脂红的预测均方根误差(RMSEP)分别为0.88,1.71 g/mL;收敛速度明显快于传统遗传算法,为多组分同 时测定提供了新的思路。 关键词:全息凹面光栅;多组分定量分析;改进遗传算法;自适应寻优空间;小生境技术 中图分类号:TH 744.1 文献标识码:A 文章编号:1000-9787(2015)06--0046--04 Research 0n quantitative analysis system for enzyme.1inked ● ' :l: lmmUn0aSSaV analVZer CHANG Hm.tao ,ZHU Lian.qing2。WANG Zhong.yu ,ZHOU Zhe—h ,GUO Yang.kuan (1.School of Instrumentation Science&Opto-Electronics Engineering,Beihang University, Beijing 100191,China;2.Beijing Key Laboratory for Optoelectronic Measurement Technology, Beijing Information Science&Technology University.Beijing 100192,China) Abstract:In order to improve detection precision of automatic enzyme—linked immunoassay analyzer,a correction model for muhieomponent quantitative analysis based on improved genetic algorithm is developed.Through adaptive search space technology and niche technology are proposed to improve performance of genetic algorithm. A UV-Vis spectrophotometer is built by using holographic concave grating as spectroscope.The comparison is made between the proposed method and least squares regression method by application to simultaneous determination of amaranth and carmine samples which arc commonly used in the food industry.Experimental results indicate that prediction precision of the improved genetic algorithm and LS are similar,for amaranth and carmine are 0.88,1.71 mL,respectively;and convergence rate is signiifcantly faster than ordinary GA.It also provides a new idea for muhicomponent simultaneous determination. Key words:holographic concave grating;multicomponent quantitative analysis;improved genetic algorithm; adaptive search space;niche technology 0引 言 间的定量校正模型 。 全自动酶免分析仪是临床检验必备的分析仪器,主要 近年来,遗传算法 (genetic algorithm,GA)作为一种全 用于传染病血清学标志物、肿瘤标志物及内分泌等各种免 局搜索方法,被广泛应用于组合优化、自适应控制 、图像 疫指标的检测 。传统的酶免分析仪受限于其单色器滤 处理 和光谱分析 等领域。特别是针对光谱分析中的 光片组数量,只能进行定性分析。为实现定量分析,需采集 波长选择问题,国内外学者的研究表明采用遗传算法进行 样本全光谱数据,并借助于化学计量学方法,如,偏最小二 波长选择,即可以剔除不相关变量、简化模型,又可提高模 乘(partial least square,PLS) 、主成分分析(principal com— 型的预测能力和稳健性 。但是遗传算法在使用过程中 ponent analysis,PCA)[3 J,建立光谱数据和待测样本属性之 仍存在一些问题,如,“早熟现象”和种群进化效率低等 。 收稿日期:2o14-1o-09 基金项目:北京市自然科学基金重点资助项目(K O1Ol0772O32);教育部“长江学者与创新团队”发展计划资助项目(IRT1212);国防技 术基础科研资助项目(J132012C001) 第6期 常海涛,等:全自动酶免分析仪定量分析系统研究 47 本文将多组分定量分析转化成组合最优化问题,提出 =O.2,为绝对位置差异的权重因子和辨识系数;A ---7/: 0.8为变化率的权重因子和辨识系数。 了一种基于改进遗传算法的定量分析方法。通过自适应寻 优空间和小生境技术加快种群的进化速度,提高算法的整 体性能。以全息凹面光栅作为分光元件,搭建扫描式紫外 可见分光光度计实验系统,并获得食品工业中最常见的苋 菜红和胭脂红二组分混合水溶液的吸收光谱数据。 1算法与原理 1.3遗传算法的改进 1)自适应寻优空间:首先随机产生初始种群,计算个 体适应度;然后以每代种群中最大适应度个体c 为中心生 成下一代解空间[c。+do,C。一 ],并在新的解空间中随机 产生新的种群;重复上述操作,经过固定遗传代数,锁定搜 索空间,引人选择、交叉、变异等遗传操作实现种群进化。’ 2)小生境技术:本文采用个体之间的广义海明距离作 为共享函数,对相似个体中适应度较小的个体施加一个较 1.1 回归模型的建立 可见光谱数据通常呈多变量、强相关性,并与样本组分 含量构成线性关系,符合朗伯比尔定律,即 A =C ・量,爪 ・b. (1) 强的惩罚函数,降低其适应度值,从而达到滤除相似个体, 增加种群的多样性,加快收敛速度的目的。 2实验系统搭建 其中, 为吸收系数矩阵,mL/( g・cm);6为光程, cm;C为样本浓度矩阵, mL;A为吸光度矩阵,n为样本 数,P为测量波长数,m为待测样本的组分数。由于建模用校 正集样本的吸光度和浓度矩阵已知,而且样本数通常大于组 分数,因此,可通过求超定方程最小范数最小二乘解的方法 得到各组分的吸收系数。对浓度矩阵c进行奇异值分解 C=U・S・VT. (2) 2.1分光光度系统 全息凹面光栅,具有像差校正、低杂散光和高信噪比等 优点。以全息凹面光栅为核心的分光光度系统的结构如 图1所示。 狭缝 狭缝 其中, 为n×n阶列正交矩阵, 为n×m阶对角阵;V 为m×m阶列正交矩阵。吸收系数矩阵,即超定方程解为 = ・(S /( ・A )). (3) 图1分光光度系统框图 Fig 1 Block diagram of spectrophotometer 对于每一个待测集样本,利用遗传算法在浓度解空间搜 索全局最优浓度解c ,带入到式(1),使得c .・K・b等 于或接近待测样本吸光度,C 即为待测样本的 州 。 从图1可以看出:系统主要由光源、单色器、探测器、控 制及处理电路和上位机组成。光源选用20 w碘钨灯,辐射 波长范围为320—800 rim,由于光源发热量大,采用外置结 1.2遗传算法 遗传算法将达尔文“适者生存”理论引入到种群进化过 程中,伴随着交叉、变异等遗传操作,实现种群中个体的信 息交换,并不断产生出新一代更适合环境的种群,最后收敛 于一个最适应环境的个体。本文采用的遗传算法操作参数 设定为:二进制编码,种群规模40,选择算子采用随机遍历 抽样,单点交叉(概率70%),变异概率1%,并采用精英主 义策略保留最优个体。需要特别指出的是,本文选取待测 样本预测吸光度曲线(即C ・K・6)和实测吸光度曲线 构,并配有散热片。单色器采用全息凹面光栅作为色散元 件,与传统的C—T结构光路相比,无需准直和聚焦成像系 统,简化系统结构的同时减少光能量的损失。样本杯为 1 cm见方石英材质的比色皿。探测器采用滨松R636—10光 电倍增管(光谱响应范围160—930 rim) 榨制和数据处理 部分,以TMS320F2812 DSP为核心处理器,实现光源供电、 光栅转台控制、探测器信号采集及处理和上位机数据传输 与通信。需要指出的是,全息凹面光栅固定于可微调的光 栅转接架上,确保光栅最凹处的切线与转台的中心轴线在 一的灰色综合关联度函数作为适应度函数。灰色综合关联度 既考虑了吸光度曲线间的绝对位置变化,又兼顾了曲线自 身的相对变化率,以综合集成的形式表征曲线间的相似程 度 ,公式如下 l+^ 个平面上。步进电机驱动蜗轮蜗杆机构来带动光栅转台 的连续转动,从而实现全谱段波长扫描。单色器实物如图2 所示。 1+A lAxm(k)I+A。I△ 。 (k)l+朋H ‘ (4) 其中,MAX 。maxl 。 ( )I+nm ̄IAx 。 ( )I; 图2单色器实物图 Fig 2 Monochrometer 为实测吸光度序列; 为预测吸光度序列; 。 (k)= (k)一 ( )为序列间的绝对位置差;Ax 。 (k)=( 。(k)一 。‰(k一1))一( ( )一 (k一1))为序列变化率差异。A = 传感器与微系统 第34卷 2.2实验材料与方法 待测样本选用食品工业中最常见的苋菜红、胭脂红两 种食用色素(国家标准物质研究中心提供)。采用精度为 10 手工加样器,吸取苋菜红和胭脂红纯组分标准溶液, 分别置于100 mL容量瓶中,用蒸馏水定容配制成两组分混 合溶液,其中各组分浓度在0 ̄200 ̄g/mL范围内均匀分布, 任意组合。选取蒸馏水为参比溶液,在波长381—630 llm (间隔1 itm)内,记录35组混合色素样本吸收光谱数据,其 中20组作为校正集,l5组作为预测集。经Savitzky—Golay 窗口移动多项式平滑后35组混合色素样本吸收光谱如图3 所示。从图3可以看出:光谱曲线平滑,滤噪效果明显,并 且苋菜红和胭脂红吸收峰接近,光谱重叠严重。 ≤ 图3平滑后苋菜红和胭脂红吸收光谱 Fig 3 Absorption spectra of Amaranth and Carmine wim smoothing 3实验结果与讨论 3.1 吸收系数计算 将校正集样本吸光度和浓度矩阵带人到式(2)和 式(3)中,通过解超定方程的方法计算苋菜红和胭脂红样本 在不同波长点的吸收系数。为进一步验证吸收系数的准确 性,将校正集样本浓度和两组分的吸收系数带入到公 式(1),反推出校正集样本的吸光度A。 (预测吸光度),计 算A 和A(实测吸光度)中各样本数据序列间的灰色综合 关联度。图4给出了校正集2O个样本预测吸光度A。 和 实测吸光度A的灰色综合关联度变化曲线,其中二者的灰 色综合关联度值均大于0.998,表明本文采用的解超定方程 法求解样本吸收系数的准确度很高。 岱 :H< 如 蜷 蝴 图4 Apred和A的灰色综合关联度曲线 Fig 4 Synthesized grey correlation(SGC)between Apredand A 3.2划分寻优空间 以胭脂红(100 mL)和苋菜红(80 mL)混合样本 为例,图5给出寻优空间随遗传代数变化示意,边界变化值 d。=5 mL。由于胭脂红和苋菜红吸收峰接近,吸收光谱 存在干扰,并且初始种群的解空间为[0,200],因此,生成的 第一代寻优空间,可能并不包含最优解。自适应寻优空问 技术能迅速定位全局最优解范围,本例中在第5代寻优空 间开始完全包含全局最优解,当遗传算法演化超过10代 时,寻优空间基本保持不变。最终胭脂红的寻优空间为 [95.19,105.19],苋菜红寻优空间为[74.35,84.35],均包 含全局最优解,且搜索空间缩小了20倍。 2O 寻优空间(1o代): +胭脂红 辍 15 器籍粒 19, 1 053.51J9 *苋莱 1o 鲴 5 40 6O 8O 1OO 120 样本浓度/p,g・mL 图5寻优空间随遗传代数变化曲线 Fig 5 Change curve of search space with evolutionary generations 3.3算法比较 分别采用遗传算法和改进遗传算法对预测集样本浓度 进行预测,以预测残差平方和(prediction residual error sum of square,PRESS)评价模型预测的准确性。图6(a),(b)分 别给出了浓度为50&80,60&60,80&100,lOO&8Op ̄g/mL的 苋菜红和胭脂红混合样本,两种方法每代最优个体的 PRESS值随遗传代数变化曲线。 从图6(a)可以看出:遗传算法预测结果存在不确定 性,由于早熟现象影响,a,b,d样本经过1000代进化仍收敛 于局部最优解;C样本经过69代遗传进化收敛于全局最优 解。改进遗传算法由于搜索空间进行了动态调整,每个样 本均能迅速收敛于全局最优解。由图6(b)可知,4个样本 均在30代以内收敛于全局最优解。 极 fa1遗传算法 g i 榭 第6期 常海涛,等:全自动酶免分析仪定量分析系统研究 49 系数和标准差基本相当。说明基于改进遗传算法的多组分 定量分析模型具有较高的预测精度。 表1 改进遗传算法与PLS预测性能比较 Tab 1 Comparison of prediction performance index by IGA and PLSmethod [4]Lee S,Kim K,Lee H,et a1.Improving the classiifcation accuracy f0r IR spectroscopic diagnosis of stomach and colon malinancy U—g uv叵 sing non-linear spectral feature extraction methods[J].Analyst, 2013,138:4076--4082. ∞ 舳 ∞ ∞ [5]Shao X G,Wang W,Hou Z Y,et a1.A new regression method based on independent component analysis[J].Talanta,2006, 69(3):676-680. [6]Clavaud M,Roggo Y,Daeniken R V,et a1.Chemometrics and in— line near infrared spectroscopic monitoring of a biopharmaceutical 4结论 Chinese hamster ovary cell culture:Prediction of multiple cuhiva— 本文提出了一种基于改进遗传算法的多组分定量分析 tion variables[J].Talanta,2013,111:28--38. 方法,将多组分定量分析转化成组合最优化问题,利用遗传 算法在浓度解空间搜索全局最优浓度解。采用自适应寻优 [7]Shen z H,Zhao Y K,wu w W.Niche pseudo—parallel genetic al— gorithms for path optimization of autonomous mobile robot[J]. Journal of Shanghai University:English Edition,2006,10(5): 449-453. 空间技术和小生境技术对传统遗传算法进行改进,一定程 度上加快了遗传算法收敛速度,克服“早熟现象”。搭建以 凹面全息光栅为核心的酶免分析仪分光光度系统,对食品 工业中最常见的苋菜红和胭脂红混合溶液定量分析。其中 苋菜红和胭脂红待测样本参考值与模型预测值间的相关系 数R分别达到了0.997 9,0.991 8,RMSEP分别为0.88, 1.71 g/mL。实验结果表明:本文所研究的方法具有较高 [8] 陈华,叶东,陈刚,等.遗传算法的数字图像相关搜索 法[J].光学精密工程,2007,15(10):1633—1637. [9] 曹晖,周延.多种群精英共享遗传算法在异常光谱识别 中的应用[J].光谱学与光谱分析,2011,31(7):1847--1851. [1O]Fei Q,Li M,Wang B,et a1.Analysis of cephalexin with NIR spec- tromety couplred to artiifcil neural netaworks with modified gene- 的预测精度,满足酶免分析仪定量分析的要求。 参考文献: [1]Dilorenzo M E,Timoney C F,Felder R A.Technological advance- ments in liquid handling robotics[J].Journal of Laboratory Auto— marion,2001,6(2):36-40. tic algorithm for wavelength selection[J].Chemometrics and In— telligent Laboratory Systems,2009,97:127--131. [11]Thakur M.A new genetic algorithm for global optimization of mul— timodal continuous functions[J].Journal of Computational Sci— ence,2014,5:298--311. [2] 乔晓艳,王艳景,李刚.偏最小二乘法荧光光谱预测啶虫脒 [12]刘思峰,党耀国,方志耕,等.灰色系统理论及其应用[M].北 京:科学出版社,2010:6--9. 农药残留[J].光学精密工程,2010,18(11):2369-2374. [3]Stanimimva I,Walczak B,Massart D L,et a1.A comparison be- tween two robust PCA algorithms[J].Chemometrics and Intelli— gent Laboratory Systems,2004,71:83-95. 、 作者简介: 常海涛(1982一),男,河北唐山人,博士研究生,主要从事光谱 分析、自动化仪表等方面的研究。 (上接第45页) 2011,48(8):516--519. [3]Hanse J G.Honeywell MEMS inertial technology&product sta- tus[C]∥ION NTM 2004,San Diego,California,USA,2004. [4] 何昆鹏.MEMS惯性器件参数辨识及系统误差补偿技术[D]. 哈尔滨:哈尔滨工程大学,2009. [5] Weinberg Harvey.Using the ADXL202 in pedometer and personal navigation applications[Z].Analog Devices Inc,2002. 向/m 图8系统生成的轨迹图 Fig 8 Trajectory generated by system [6] 宋丽君.基于MEMS器件的航向姿态测量系统的研究[D]. 西安:西北工业大学,2007. 本系统通过加速度计来计算行走距离,并通过磁强计和陀 螺仪来确定运动方向。系统硬件设备简单,成本低廉,能够 实时获得人的位置信息。最后通过匹配算法对位置信息进 行合理的修正,提高了测试轨迹的准确性。 参考文献: [1]汪苑,林锦国.几种常用室内定位技术的探讨【j].中国仪 [7]赵 琳.MEMS惯性器件的数据采集与误差分析技术研 究[D].哈尔滨:哈尔滨工程大学,2009. [8] 任晶晶,郭宇明,一种改进的进路搜索方法[J].仪器仪表用 户,2005(1):70. [9] 王 敏,赵晓雷.基于遍历搜索二叉树中最长路径的算法研 究[J].软件技术,2010(8):54--56. 作者简介: 张迅(1986一),男,山东潍坊人,硕士研究生,主要研究方向 为WLAN分析与优化。 器仪表,201I(2):54--57. [2]王淑华.MEMS传感器现状及应用[J].MEMS与传感器,