劉 勇,高 娜,李志剛,朱 靈,李 飛,張 龍
(中國科學院安徽光學精密機械研究所光電子技術(shù)研究室,安徽合肥 230031)
聚合酶鏈式反應(polymerase chain reaction,PCR)是20世紀80年代發(fā)展起來的體外快速擴增特定基因或DNA序列的技術(shù)。它具有特異、敏感、快速、簡便、產(chǎn)率高等優(yōu)點,在分子生物學和基因工程實驗室、疾病檢測、臨床應用、商品檢疫、新醫(yī)藥開發(fā)和工農(nóng)業(yè)制品等方面顯示了巨大潛力和廣泛用途[1]。微流控PCR是在微流控芯片上加工的微反應室中進行PCR,利用芯片體積小、比表面積大和集成度高等特點,實現(xiàn)芯片上的DNA快速擴增,具有耗樣少、擴增效率高、速度快和可集成等優(yōu)勢,在國內(nèi)外已得到廣泛研究和商業(yè)化應用[2]。
PCR由高溫變性、低溫退火和適溫延伸組成一個周期,循環(huán)進行,使目的基因片段得以迅速擴增。其中,溫度控制的精度非常重要,直接影響PCR實驗成功與否,較大的溫度偏差會導致PCR實驗失?。?]。近年來研究表明[4],提高PCR熱循環(huán)中3個溫區(qū)之間的升降溫速率,不僅能夠縮短整個PCR擴增過程所需時間,而且特異性更好。因此,PCR反應中的溫度控制不僅涉及升溫、降溫和恒溫過程,還要求溫度控制精度高,超調(diào)小,調(diào)節(jié)時間短。在溫度控制中,微流控芯片周圍環(huán)境溫度也會對其造成干擾,因此,這種溫度控制具有時變、非線性等不利因素。如果采用傳統(tǒng)的PID算法,控制參數(shù)選擇較困難,系統(tǒng)存在抗干擾能力差、超調(diào)量大等缺點[5,6]。
本文采用系統(tǒng)辨識的方法建立了被控對象的數(shù)學模型,利用前饋控制對周圍環(huán)境溫度進行補償,并根據(jù)所建模型的動態(tài)特性提出了超前校正的控制方法,在實際系統(tǒng)調(diào)試中,又根據(jù)現(xiàn)場的條件采用了分段控制,成功實現(xiàn)了快速升降溫和精確控溫。
溫度控制結(jié)構(gòu)如圖1,溫度控制系統(tǒng)框圖如圖2所示,主要由以下幾部分組成:溫度信號檢測與輸入、控制器、執(zhí)行機構(gòu)和微流控芯片。采用Pt100熱敏電阻器測溫[7],利用恒壓驅(qū)動方式將Pt100隨溫度的阻值變化轉(zhuǎn)換成電壓信號,放大后經(jīng)過16位A/D轉(zhuǎn)換器MX7705將數(shù)字信號反饋到輸入端,和設定溫度值比較??刂破鞲鶕?jù)溫度誤差e調(diào)節(jié)超前校正和積分環(huán)節(jié)的控制參數(shù),輸出可變占空比的PWM信號至執(zhí)行機構(gòu)。溫度控制的執(zhí)行器采用了半導體制冷片,半導體制冷片既能實現(xiàn)加熱,又能實現(xiàn)制冷,通過輸入電流的控制可實現(xiàn)高精度的溫度控制[8]。輸出信號經(jīng)過驅(qū)動電路實現(xiàn)半導體制冷片加熱或制冷,從而使微流控芯片PCR反應室溫度達到所設定的溫度值。
圖1 溫度控制結(jié)構(gòu)示意圖Fig 1 Schematic diagram of temperature control structure
圖2 溫度控制系統(tǒng)框圖Fig 2 Block diagram of temperature control system
對于微流控PCR溫度控制系統(tǒng),采用系統(tǒng)辨識的方法建立了被控對象的數(shù)學模型,其傳遞函數(shù)為G(s)=7.5/s(0.7s+1)4。由于系統(tǒng)受周圍環(huán)境溫度的影響,會損失部分熱量,使溫度降低,利用前饋控制對其進行了溫度補償
其中,Tsetpoint為設定溫度,Tambient為周圍環(huán)境溫度,Kc為補償因子。通過前饋控制輸出一個附加量,補償了損失的熱量,提高系統(tǒng)對周圍環(huán)境的抗干擾能力。
對所建模型進行開環(huán)頻率分析,發(fā)現(xiàn)其相位裕量為負值,要使最小相位系統(tǒng)穩(wěn)定,相位裕量必須為正值,故設計了超前校正控制器。超前校正的傳遞函數(shù)[9]為
式中 α稱為分度系數(shù),T為時間常數(shù)。超前校正具有正的相角特性,能提供的最大正相角為
根據(jù)穩(wěn)態(tài)誤差確定系統(tǒng)增益K,繪制未校正系統(tǒng)的對數(shù)頻率特性,確定出未校正系統(tǒng)的相位裕量。根據(jù)關系式(4)確定超前校正提供的最大正相角
實際使用過程中發(fā)現(xiàn)系統(tǒng)穩(wěn)態(tài)時溫度存在誤差,增加積分環(huán)節(jié)可以改善系統(tǒng)的穩(wěn)態(tài)性能,減小穩(wěn)態(tài)誤差。但是當選用較大的積分因子時,系統(tǒng)易出現(xiàn)超調(diào)和振蕩;而選用較小的積分因子時,系統(tǒng)基本無超調(diào),但是調(diào)整時間較長,導致升降溫速率較慢。因此,采用積分—超前校正算法很難同時實現(xiàn)快速升降溫和精確控溫。
為了實現(xiàn)快速升降溫和精確的控溫,提出了分段控制,設置一個誤差閾值β,根據(jù)溫度誤差e大小采用不同的控制參數(shù),其規(guī)則如下:
1)當|e|>β時,為使系統(tǒng)具有較好的動態(tài)特性,提高系統(tǒng)響應速度,縮短調(diào)整時間,選取較大的增益;同時為避免系統(tǒng)響應出現(xiàn)較大超調(diào),積分因子為0。
2)當|e|≤β時,為使系統(tǒng)具有較好的穩(wěn)定性,避免引起振蕩,減小增益;同時為了減小穩(wěn)態(tài)誤差,增加積分環(huán)節(jié),選擇合適的積分因子。
分段控制使參數(shù)調(diào)節(jié)更靈活,不僅發(fā)揮了積分環(huán)節(jié)和增益的優(yōu)勢,同時還避免了引起系統(tǒng)不穩(wěn)定。
分別采用積分—超前校正算法和分段控制算法對系統(tǒng)進行測試,圖3的曲線是系統(tǒng)從72℃升溫到95℃的數(shù)據(jù),圖4的曲線是系統(tǒng)從95℃降溫到60℃的數(shù)據(jù)。由圖中可以看出:采用積分—超前校正算法時穩(wěn)態(tài)誤差較小,但是系統(tǒng)調(diào)整時間較長,升降溫速率較慢;采用分段控制算法時,系統(tǒng)調(diào)整時間變短,平均升溫速率4.46℃/s,平均降溫速率4.36℃/s,同時又保持了較小的穩(wěn)態(tài)誤差,穩(wěn)態(tài)誤差不大于0.5℃。
將系統(tǒng)用于微流控PCR實驗。圖5是系統(tǒng)實際的溫度控制過程中部分周期的溫度循環(huán)曲線,其中,變性階段95℃持續(xù)10 s,退火階段53℃持續(xù)15s,延伸階段72℃持續(xù)20s,一個溫控循環(huán)用時60s,升降溫速率較快,超調(diào)很小,完全可以滿足微流控PCR實驗的需求。
圖3 升溫曲線Fig 3 Heating curve
圖4 降溫曲線Fig 4 Cooling curve
圖5 PCR溫度循環(huán)曲線Fig 5 PCR temperature cycling curves
根據(jù)微流控PCR快速升降溫、精確控溫的要求,設計了一種分段式超前校正控制器,使系統(tǒng)不僅具有較好的動態(tài)特性,還具有較好的穩(wěn)態(tài)特性。實驗結(jié)果表明:分段控制算法能夠靈活的調(diào)整參數(shù),具有較好的控制品質(zhì),不僅縮短了系統(tǒng)調(diào)整時間,取得較好的升降溫速率,還保持了較小的穩(wěn)態(tài)誤差。目前,集成該算法的溫控系統(tǒng)已滿足微流控PCR實驗的需求,計劃用于致病菌PCR實驗。
[1] 陳文元,張衛(wèi)平.集成微流控聚合物PCR芯片[M].上海:上海交通大學出版社,2009:29-34.
[2] 姚李英,陳濤,王升啟,等.高聚物基PCR微流控芯片技術(shù)[J].中國生物工程雜志,2004,24(3):90 -93.
[3] 馬建崗.基因工程學原理[M].2版.西安:西安交通大學出版社,2007:82-96.
[4] Wittwer C T,Garling D J.Rapid cycle DNA amplication:Time and temperature optimization[J].BioTechniques,1991,10(1):76-83.
[5] 畢雪芹,倪 原,王麗娟.PCR溫度模糊PID控制器的設計與仿真[J].計算機工程與設計,2008,29(10):2604 -2606.
[6] 邱憲波,袁景淇.帶前饋補償?shù)幕驍U增反應高精度溫控系統(tǒng)[J].儀器儀表學報,2005,26(8):831 -832.
[7] 松井邦彥.傳感器應用技巧141例[M].梁瑞林,譯.北京:科學出版社,2006:77 -92.
[8] 陳振林,孫中泉.半導體制冷器原理與應用[J].微電子技術(shù),1999,27(5):63 -65.
[9] 龐國仲,自動控制原理[M].合肥:中國科學技術(shù)大學出版社,2007:322-326.