蔣理++錢(qián)家好++陳曉東++王利軍

摘要：開(kāi)發(fā)藥代動(dòng)力學(xué)計(jì)算軟件可以利用現(xiàn)在非常普及的個(gè)人計(jì)算機(jī)解決二次微分方程計(jì)算和曲線(xiàn)繪圖的問(wèn)題。本項(xiàng)目開(kāi)發(fā)的軟件可以提高相關(guān)檢測(cè)部門(mén)在相關(guān)計(jì)算上的速度與準(zhǔn)確度。

國(guó)內(nèi)外比較多見(jiàn)的相關(guān)軟件主要有NONLIN，NONMEN，3P87/3P97，NDST，ABE等。均為非開(kāi)源軟件，使用的計(jì)算方法不一。其中3P97對(duì)樣本數(shù)據(jù)的擬合，計(jì)算房室模型參數(shù)計(jì)算部分可以做到免費(fèi)。本次研究的藥代動(dòng)力學(xué)計(jì)算軟件可以與3P97配合使用，可以滿(mǎn)足大部分相關(guān)科研院所，學(xué)校師生的需求。

關(guān)鍵詞：藥代動(dòng)力學(xué) 計(jì)算 微分方程 曲線(xiàn)圖

中圖分類(lèi)號(hào)：TP313 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1007-9416（2015）09-0000-00

1 項(xiàng)目背景

在畜牧業(yè)、漁業(yè)養(yǎng)殖領(lǐng)域，用藥物防治動(dòng)物疾病目前是一種最直接、最有效和最經(jīng)濟(jì)的措施，但藥物殘留對(duì)人類(lèi)健康和生態(tài)環(huán)境存在著嚴(yán)重的威脅。如何控制藥物在動(dòng)物體內(nèi)中的殘留，以保證產(chǎn)品在上市時(shí)能符合相關(guān)標(biāo)準(zhǔn)，是相關(guān)質(zhì)量檢測(cè)部門(mén)必須要解決的問(wèn)題。動(dòng)物在給藥后，藥物在動(dòng)物體會(huì)按一定的規(guī)律消除。掌握藥物在生物體內(nèi)的消除時(shí)間，是解決動(dòng)物產(chǎn)品質(zhì)量安全的重要科學(xué)方法，也是科學(xué)發(fā)展觀(guān)在畜牧業(yè)、漁業(yè)養(yǎng)殖上的具體體現(xiàn)。只有嚴(yán)格遵守休藥期，遵循藥物代謝動(dòng)力學(xué)，才能做到科學(xué)用藥、安全用藥。

描述某種藥物的代謝動(dòng)力學(xué)特征，最常用的方法是將機(jī)體看成由房室組成，藥物在這些房室間的運(yùn)轉(zhuǎn)及從房室中消除的速率符合一級(jí)過(guò)程或線(xiàn)性動(dòng)力學(xué)。只要機(jī)體內(nèi)某些部位吸收、運(yùn)轉(zhuǎn)或消除藥物的速率常數(shù)相似，即可歸入一個(gè)房室，而不受解剖學(xué)位置和學(xué)功能的限制。

本項(xiàng)目開(kāi)發(fā)的軟件可以提高相關(guān)檢測(cè)部門(mén)在相關(guān)計(jì)算上的速度與準(zhǔn)確度。國(guó)內(nèi)外比較多見(jiàn)的相關(guān)軟件主要有NONLIN，NONMEN，3P87/3P97，NDST，ABE等。均為非開(kāi)源軟件，使用的計(jì)算方法不一。其中3P97對(duì)樣本數(shù)據(jù)的擬合，計(jì)算房室模型參數(shù)計(jì)算部分可以做到免費(fèi)。本次研究的藥代動(dòng)力學(xué)計(jì)算軟件與3P97配合使用，可以滿(mǎn)足大部分相關(guān)科研院所，學(xué)校師生的需求。目前市面上實(shí)現(xiàn)相關(guān)功能的軟件也沒(méi)有實(shí)現(xiàn)免費(fèi)，相關(guān)技術(shù)的實(shí)現(xiàn)也沒(méi)有做到開(kāi)源，所以本軟件的開(kāi)發(fā)具有實(shí)際的市場(chǎng)價(jià)值。

2 主要算法分析

一般藥代動(dòng)力學(xué)的計(jì)算都涉及到二房室模型，休藥期二房室計(jì)算公式為：二次微分方程 ，手動(dòng)解指數(shù)方程無(wú)疑會(huì)給科研人員帶來(lái)比較大的麻煩，同時(shí)殘留量的衰減曲線(xiàn)圖的繪制也很麻煩，手工繪圖也難免帶來(lái)誤差。開(kāi)發(fā)藥代動(dòng)力學(xué)計(jì)算軟件可以利用現(xiàn)在非常普及的個(gè)人計(jì)算機(jī)解決休藥期計(jì)算的問(wèn)題。

解二次微分方程的算法主要思想是通過(guò)迭代的方法，選擇適當(dāng)?shù)牟介L(zhǎng)和精確度ε，將微分方程的計(jì)算值與漸出現(xiàn)檢出限進(jìn)行比較，使得出的理論休藥期x盡量趨近真值，同時(shí)也要兼顧到少或避免發(fā)生得不到值的情況或由于循環(huán)時(shí)間過(guò)長(zhǎng)而產(chǎn)生死循環(huán)的假象。同時(shí)在每次計(jì)算的時(shí)候描點(diǎn)作圖，將整個(gè)計(jì)算結(jié)果直接顯示在曲線(xiàn)圖上，使結(jié)果反映更加直觀(guān)。迭代求解流程如圖1。

在求解二次微分方程時(shí)，先將x設(shè)為0，然后進(jìn)入循環(huán)，循環(huán)內(nèi)計(jì)算 的值，將得到的值與最低殘留限值作比較，誤差小于設(shè)置誤差值時(shí)，判定X為所求值，否則X+=b，b為步長(zhǎng)，繼續(xù)循環(huán)，在下一期計(jì)算前，判斷X是否小于3000*b，如小于繼續(xù)計(jì)算，否則判定X無(wú)解，根據(jù)大量的實(shí)際數(shù)值的計(jì)算，3000次循環(huán)足夠得到X值。

圖1 迭代求解二階微分方程流程圖

3 系統(tǒng)實(shí)現(xiàn)

藥代動(dòng)力學(xué)計(jì)算軟件主要有四個(gè)模塊：計(jì)算與繪圖，數(shù)據(jù)管理，數(shù)據(jù)報(bào)表與打印，用戶(hù)管理。

在計(jì)算與繪圖模塊，用戶(hù)輸入分布相的零時(shí)截距A， 消除相的零時(shí)截距B， 分布速率常數(shù)α，消除速率常數(shù)β，最低殘留限量Y，允許誤差大小ε，置信度σ，可以計(jì)算出理論置信時(shí)間范圍，并畫(huà)出圖像。同時(shí)要求輸入平均溫度，以方便數(shù)據(jù)分析。對(duì)不能算出置信時(shí)間范圍的數(shù)據(jù)會(huì)提示無(wú)法算出，輸入數(shù)據(jù)有明顯問(wèn)題的，程序也會(huì)給予提醒。

數(shù)據(jù)管理模塊主要完成數(shù)據(jù)從Excel表格中導(dǎo)入導(dǎo)出。經(jīng)調(diào)研，科研機(jī)構(gòu)現(xiàn)階段數(shù)據(jù)模型的數(shù)量并不是特別巨大，而且研究人員也習(xí)慣于將數(shù)據(jù)存儲(chǔ)在Excel表格內(nèi)，所以軟件不再另外建立數(shù)據(jù)庫(kù)，直接對(duì)Excel表格中的數(shù)據(jù)進(jìn)行處理。通過(guò)數(shù)據(jù)導(dǎo)入功能，可以直接將Excel表中的數(shù)據(jù)導(dǎo)入到計(jì)算軟件的面板，點(diǎn)擊某一行時(shí)會(huì)將該行的數(shù)據(jù)導(dǎo)入到參數(shù)輸入欄中并計(jì)算繪制曲線(xiàn)。如果新算的數(shù)據(jù)要加入到Excel表中，只需要使用新增功能即可。對(duì)于行列不符合規(guī)范的表格，系統(tǒng)會(huì)提醒數(shù)據(jù)格式不符合規(guī)范，無(wú)法打開(kāi)。

數(shù)據(jù)報(bào)表與打印主要支持?jǐn)?shù)據(jù)計(jì)算后生成報(bào)告，幫助科研人員迅速生成文檔，除了參數(shù)，計(jì)算結(jié)果和圖片外，用戶(hù)可以自己編輯備注信息加到報(bào)表中。

在用戶(hù)管理模塊，為了方便用戶(hù)安裝使用，也不需要加裝其他數(shù)據(jù)庫(kù)軟件，用戶(hù)名和密碼經(jīng)加密后存儲(chǔ)在本地文件中。

在系統(tǒng)實(shí)現(xiàn)過(guò)程中，擴(kuò)展實(shí)現(xiàn)了藥代動(dòng)力學(xué)三房室模型與一房室模型的計(jì)算與繪制曲線(xiàn)圖的方法，三房室模型的實(shí)現(xiàn)方法與二房室類(lèi)似，增加了一個(gè)計(jì)算因子，一房室模型在時(shí)間軸上用對(duì)數(shù)表示，曲線(xiàn)圖為一條直線(xiàn)。在系統(tǒng)實(shí)現(xiàn)時(shí)重新繪制了相關(guān)面板，可以在計(jì)算模型選擇菜單欄中靈活切換這幾種模型。系統(tǒng)實(shí)現(xiàn)效果如圖2。

圖2 系統(tǒng)實(shí)現(xiàn)效果圖

參考文獻(xiàn)

[1]芮建中，盧建豐，凌樹(shù)森.NONMEN程序及其在臨床藥理學(xué)中的應(yīng)用.藥學(xué)進(jìn)展，1997.21.

[2]謝海棠，黃曉暉，孫瑞元.中國(guó)臨床藥理學(xué)與治療學(xué)，2001.4.

[3]Y. Daniel Liang Java.語(yǔ)言程序設(shè)計(jì).基礎(chǔ)篇（原書(shū)第8版），2011.6.

收稿日期：2015-08-05

作者簡(jiǎn)介：蔣理（1982—），男，江蘇南京人，南京信息職業(yè)技術(shù)學(xué)院教師，碩士研究生，研究方向：數(shù)據(jù)交換技術(shù)。錢(qián)家好（1995—），男，江蘇南京人，南京信息職業(yè)技術(shù)學(xué)院學(xué)生，專(zhuān)業(yè)：移動(dòng)互聯(lián)網(wǎng)應(yīng)用技術(shù)。陳曉東（1994—），男，江蘇宿遷人，南京信息職業(yè)技術(shù)學(xué)院學(xué)生，專(zhuān)業(yè)：移動(dòng)互聯(lián)網(wǎng)應(yīng)用技術(shù)。王利軍（1995—），女，河南商丘人，南京信息職業(yè)技術(shù)學(xué)院學(xué)生，專(zhuān)業(yè)：移動(dòng)互聯(lián)網(wǎng)應(yīng)用技術(shù)。