基于Java的药代动力学计算软件的实现.docVIP

基于Java的药代动力学计算软件的实现 　　摘要：开发药代动力学计算软件可以利用现在非常普及的个人计算机解决二次微分方程计算和曲线绘图的问题。本项目开发的软件可以提高相关检测部门在相关计算上的速度与准确度。 　　国内外比较多见的相关软件主要有NONLIN，NONMEN，3P87/3P97，NDST，ABE等。均为非开源软件，使用的计算方法不一。其中3P97对样本数据的拟合，计算房室模型参数计算部分可以做到免费。本次研究的药代动力学计算软件可以与3P97配合使用，可以满足大部分相关科研院所，学校师生的需求。 　　关键词：药代动力学 计算 微分方程 曲线图 　　中图分类号：TP313 文献标识码：A 文章编号：1007-9416（2015）09-0000-00 　　1 项目背景 　　在畜牧业、渔业养殖领域，用药物防治动物疾病目前是一种最直接、最有效和最经济的措施，但药物残留对人类健康和生态环境存在着严重的威胁。如何控制药物在动物体内中的残留，以保证产品在上市时能符合相关标准，是相关质量检测部门必须要解决的问题。动物在给药后，药物在动物体会按一定的规律消除。掌握药物在生物体内的消除时间，是解决动物产品质量安全的重要科学方法，也是科学发展观在畜牧业、渔业养殖上的具体体现。只有严格遵守休药期，遵循药物代谢动力学，才能做到科学用药、安全用药。 　　描述某种药物的代谢动力学特征，最常用的方法是将机体看成由房室组成，药物在这些房室间的运转及从房室中消除的速率符合一级过程或线性动力学。只要机体内某些部位吸收、运转或消除药物的速率常数相似，即可归入一个房室，而不受解剖学位置和学功能的限制。 　　本项目开发的软件可以提高相关检测部门在相关计算上的速度与准确度。国内外比较多见的相关软件主要有NONLIN，NONMEN，3P87/3P97，NDST，ABE等。均为非开源软件，使用的计算方法不一。其中3P97对样本数据的拟合，计算房室模型参数计算部分可以做到免费。本次研究的药代动力学计算软件与3P97配合使用，可以满足大部分相关科研院所，学校师生的需求。目前市面上实现相关功能的软件也没有实现免费，相关技术的实现也没有做到开源，所以本软件的开发具有实际的市场价值。 　　2 主要算法分析 　　一般药代动力学的计算都涉及到二房室模型，休药期二房室计算公式为：二次微分方程 ，手动解指数方程无疑会给科研人员带来比较大的麻烦，同时残留量的衰减曲线图的绘制也很麻烦，手工绘图也难免带来误差。开发药代动力学计算软件可以利用现在非常普及的个人计算机解决休药期计算的问题。 　　解二次微分方程的算法主要思想是通过迭代的方法，选择适当的步长和精确度ε，将微分方程的计算值与渐出现检出限进行比较，使得出的理论休药期x尽量趋近真值，同时也要兼顾到少或避免发生得不到值的情况或由于循环时间过长而产生死循环的假象。同时在每次计算的时候描点作图，将整个计算结果直接显示在曲线图上，使结果反映更加直观。迭代求解流程如图1。 　　在求解二次微分方程时，先将x设为0，然后进入循环，循环内计算 的值，将得到的值与最低残留限值作比较，误差小于设置误差值时，判定X为所求值，否则X+=b，b为步长，继续循环，在下一期计算前，判断X是否小于3000*b，如小于继续计算，否则判定X无解，根据大量的实际数值的计算，3000次循环足够得到X值。 　　图1 迭代求解二阶微分方程流程图 　　3 系统实现 　　药代动力学计算软件主要有四个模块：计算与绘图，数据管理，数据报表与打印，用户管理。 　　在计算与绘图模块，用户输入分布相的零时截距A， 消除相的零时截距B， 分布速率常数α，消除速率常数β，最低残留限量Y，允许误差大小ε，置信度σ，可以计算出理论置信时间范围，并画出图像。同时要求输入平均温度，以方便数据分析。对不能算出置信时间范围的数据会提示无法算出，输入数据有明显问题的，程序也会给予提醒。 　　数据管理模块主要完成数据从Excel表格中导入导出。经调研，科研机构现阶段数据模型的数量并不是特别巨大，而且研究人员也习惯于将数据存储在Excel表格内，所以软件不再另外建立数据库，直接对Excel表格中的数据进行处理。通过数据导入功能，可以直接将Excel表中的数据导入到计算软件的面板，点击某一行时会将该行的数据导入到参数输入栏中并计算绘制曲线。如果新算的数据要加入到Excel表中，只需要使用新增功能即可。对于行列不符合规范的表格，系统会提醒数据格式不符合规范，无法打开。 　　数据报表与打印主要支持数据计算后生成报告，帮助科研人员迅速生成文档，除了参数，计算结果和图片外，用户可以自己编辑备注信息加到报表中。 　　在用户管理模块，为了方便用户安装使用，也不需要加装其他数据库软件，