多峰曲线拟合..docxVIP

PAGE1 / NUMPAGES13 目录 TOC \o 1-4 \h \z \u HYPERLINK \l _Toc376003365 1 背景介绍 PAGEREF _Toc376003365 \h 2 HYPERLINK \l _Toc376003366 2 数据处理 PAGEREF _Toc376003366 \h 4 HYPERLINK \l _Toc376003367 2.1 高斯拟合 PAGEREF _Toc376003367 \h 4 HYPERLINK \l _Toc376003368 2.1.1 共轭梯度法 PAGEREF _Toc376003368 \h 5 HYPERLINK \l _Toc376003369 2.1.2 0.618法 PAGEREF _Toc376003369 \h 6 HYPERLINK \l _Toc376003370 2.1.3 计算效果 PAGEREF _Toc376003370 \h 7 HYPERLINK \l _Toc376003371 2.2 分段多项式拟合 PAGEREF _Toc376003371 \h 8 HYPERLINK \l _Toc376003372 2.2.1 多项式拟合的不足 PAGEREF _Toc376003372 \h 8 HYPERLINK \l _Toc376003373 2.2.2 分段标准 PAGEREF _Toc376003373 \h 8 HYPERLINK \l _Toc376003374 2.2.3 峰位、峰面积计算 PAGEREF _Toc376003374 \h 11 HYPERLINK \l _Toc376003375 3 软件介绍 PAGEREF _Toc376003375 \h 11 HYPERLINK \l _Toc376003376 4 结语 PAGEREF _Toc376003376 \h 12  多峰曲线拟合 1 背景介绍 工程上常采用电泳的方法来分离不同的DNA片段。带电微粒在电场作用下，向着与其电性相反的电极移动，称为电泳(electrophoresis)。利用带电粒子在电场中移动速度不同而达到分离的技术称为电泳技术。 常用微流控芯片来做为DNA混合液的载体。微流控芯片技术是把生物、化学医学分析过程的样品制备、反应、分离、检测等基本操作单元集成到一块微米尺度的芯片上，自动完成分析全过程。将DNA混合液放在特制的微流控芯片中，可以实现对溶液的精确控制，对DNA片段的检测液变得方便简单。图1为一片微流控芯片实物图。 实验室采用电泳技术来分离DNA片段。首先对溶液中的DNA片段进行染色处理，然后将溶液注入微流控芯片。在微流控芯片的沟道两端加载高压进行电泳，用激光照射经染色处理的DNA片段使发出荧光，采用光电倍增管来检测荧光的存在进而检测电泳过程,流程示意图如图2。用单片机采集光电倍增管的电流信号传入计算机得到时间-电压曲线（如图3），根据曲线中峰的位置的个数来判断溶液中DNA片段的种类。 图1 微流控芯片 图2 DNA检测分析系统示意图 图3 电泳曲线 在图3中，每一个峰代表一个DNA片段，峰的位置和面积代表着DNA片段的种类信息。为了从曲线中获得DNA信息，需要对曲线进行相应的处理。可以先对曲线作拟合处理，再由拟合得到的数学表达式来求得曲线中峰的位置和面积。 本文讨论两种曲线拟合方式：1、高斯曲线拟合；2、多项式拟合。 2 数据处理 2.1 高斯拟合 由于曲线中有多个峰，因此考虑用高斯函数和洛伦兹函数的线性组合来拟合曲线，预计会有好的结果[1]。记拟合函数为： (1) 可以看到，由洛伦兹函数和高斯函数组合而成。其中各组合系数的意义如下： ：洛伦兹函数和高斯函数的组合比例，； ：峰位置信息，； ：峰的半高宽，； ：峰强度信息，。 要对电泳曲线进行拟合，就是要求参数，，，使函数能够和曲线匹配的最好。本文采用下式来刻画这种匹配效果： (2) 其中为测得数据点。这样，拟合就转化为求的极小值问题： (3) 这是一个求多元非线性函数的极值问题，目前没有解析法直接求得，可以采用迭代方法逐步获取该函数的极小值。本文采用共轭梯度法来求解。 2.1.1 共轭梯度法 共轭梯度法是最优化中最常用的方法之一。[2]它具有算法简便、存储需求小等优点，十分适合于大规模优化问题。共轭梯度法是介于最速下降法与牛顿法之间的一个方法，它仅需利用一阶导数信息，但克服了最速下降法收敛慢的缺点，又避免了牛顿法需要存储和计算Hesse矩阵并求