蟻群算法在一種物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)上的調(diào)度研究與應(yīng)用

　　盧嘉軒，李 晉，周 延

　?。ㄉ虾４髮W(xué) 機(jī)電工程與自動(dòng)化學(xué)院 工程訓(xùn)練國(guó)家級(jí)教學(xué)實(shí)驗(yàn)示范中心，上海 200444）

　　摘要：針對(duì)一種物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)，提出了一套使用機(jī)器學(xué)習(xí)和蟻群算法的藥物調(diào)度方法。該方法使用閾值劃分的方式，對(duì)取藥記錄建立了嶺回歸、隨機(jī)森林回歸和神經(jīng)網(wǎng)絡(luò)的混合模型，以預(yù)測(cè)不同醫(yī)療箱中各類藥物的需求量。之后，根據(jù)硬件檢測(cè)的藥物實(shí)際數(shù)量，計(jì)算出各藥物的偏差值，在將調(diào)度問(wèn)題轉(zhuǎn)換為旅行商問(wèn)題以后，分別使用基本蟻群算法和最大最小蟻群算法對(duì)調(diào)度問(wèn)題進(jìn)行了求解和對(duì)比。實(shí)驗(yàn)表明，該方法能較好地預(yù)測(cè)藥物的需求量，規(guī)劃出一條合理的藥物調(diào)度路徑，為物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)的藥物調(diào)度提供了一種數(shù)據(jù)驅(qū)動(dòng)的解決方案。

　　關(guān)鍵詞：蟻群算法；機(jī)器學(xué)習(xí)；物聯(lián)網(wǎng)醫(yī)療箱；調(diào)度；旅行商問(wèn)題

　　0 引言

　　隨著智能醫(yī)療概念的興起，自動(dòng)售藥機(jī)和自助醫(yī)療箱在英美等國(guó)得到了廣泛推廣，在我國(guó)也逐漸普及。然而，傳統(tǒng)自動(dòng)售藥機(jī)在日常運(yùn)營(yíng)中需要進(jìn)行人為的藥物檢查和補(bǔ)給，在藥物的調(diào)度上也沒(méi)有一套較為合適的方案。除此之外，由于地理限制和宣傳力度欠佳，售藥機(jī)常常無(wú)人問(wèn)津，使用率較低。為此，一種基于物聯(lián)網(wǎng)的校園醫(yī)療箱系統(tǒng)應(yīng)運(yùn)而生，該系統(tǒng)由若干放置在學(xué)校不同位置的聯(lián)網(wǎng)的醫(yī)療箱所組成，可以在短時(shí)間內(nèi)為師生突發(fā)的創(chuàng)傷提供及時(shí)有效的救治藥物及必要的救援器具。該藥物箱系統(tǒng)已在某高校進(jìn)行了試運(yùn)營(yíng)，用戶可以通過(guò)微信小程序的終端查看附近的藥物箱位置和藥物余量，所有取藥記錄也將上傳至服務(wù)器，保存在數(shù)據(jù)庫(kù)中。

　　本文基于該新型物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)，提出了一種根據(jù)歷史取藥記錄和藥物余量進(jìn)行調(diào)度的通用方法。該方法分為藥物需求量預(yù)測(cè)和藥物調(diào)度路徑規(guī)劃兩大步驟。在需求量預(yù)測(cè)的實(shí)驗(yàn)中，針對(duì)歷史取藥記錄，建立了嶺回歸、隨機(jī)森林回歸和神經(jīng)網(wǎng)絡(luò)的混合模型，以預(yù)測(cè)各醫(yī)療箱中藥物的需求量。在調(diào)度路徑規(guī)劃的實(shí)驗(yàn)中，將預(yù)測(cè)的藥物需求量與實(shí)際數(shù)量進(jìn)行比較，將偏差量定義為藥物實(shí)際數(shù)量與預(yù)測(cè)需求量之差。該方法將藥物的調(diào)度問(wèn)題轉(zhuǎn)化為多個(gè)供應(yīng)商和多個(gè)需求者之間的特殊旅行商問(wèn)題(Traveling Salesman Problem,TSP)，并分別使用基本蟻群算法(Ant System,AS) ^[2] 和最大-最小蟻群算法(Max-Min Ant System,MMAS) ^[2-5] 進(jìn)行了求解和對(duì)比。該方法能為新型物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)的藥物分配與調(diào)度提供解決方案，有效降低日常運(yùn)營(yíng)成本。

　　1 需求量預(yù)測(cè)

　　由于物聯(lián)網(wǎng)醫(yī)療箱具有藥物檢測(cè)、數(shù)據(jù)上傳的特點(diǎn)，因此所有的取藥記錄和藥物余量都會(huì)保存在云服務(wù)器中，這也為模型的建立提供了必要的數(shù)據(jù)支持。對(duì)于某類存放于醫(yī)療箱中的藥品，影響其需求量的因素可以大致概括為3類：所處時(shí)間，醫(yī)療箱地理位置和藥物自身特性。針對(duì)醫(yī)療箱系統(tǒng)而言，取藥記錄間接地反映了特定時(shí)間、特定地理位置上各類藥品的需求量。為此，可以通過(guò)歷史取藥記錄的數(shù)據(jù)，建立上述3類因素到藥品需求量的映射關(guān)系。

　　具體地，所處時(shí)間包括歷史取藥日期以及對(duì)應(yīng)的月份和星期；地理位置包括醫(yī)療箱所處的經(jīng)緯度數(shù)值和地理畫像，如地圖興趣點(diǎn)(Point of Interest,POI)密度、是否位于宿舍、食堂、運(yùn)動(dòng)場(chǎng)、體育館、實(shí)驗(yàn)室等特殊場(chǎng)所等；藥物自身特性，包括藥物所屬的類別以及由行業(yè)認(rèn)可的藥物評(píng)審數(shù)據(jù)庫(kù)提供的藥物評(píng)分等。在獲取并整合上述數(shù)據(jù)以后，即可以使用機(jī)器學(xué)習(xí)算法構(gòu)建藥品需求量與各因素之間的模型，即通過(guò)輸入上述特征來(lái)預(yù)測(cè)特定日期、特定地理位置下某藥品的日需求量或多日需求量。本文使用的模型包括嶺回歸、隨機(jī)森林回歸和神經(jīng)網(wǎng)絡(luò)。

　　嶺回歸(Ridge Regression) ^[9] 是線性回歸問(wèn)題中一種改良的最小二乘估計(jì)法，即在最小二乘估計(jì)法的損失函數(shù)中加上L2正則項(xiàng)，通過(guò)放棄最小二乘法的無(wú)偏性為代價(jià)以獲得更合適的回歸系數(shù)，防止過(guò)擬合。嶺回歸由于可供訓(xùn)練的參數(shù)較少，在數(shù)據(jù)量不太大、線性可分性強(qiáng)的問(wèn)題上表現(xiàn)較優(yōu)。

　　隨機(jī)森林回歸(Random Forest Regression) ^[9-11]是一種針對(duì)分類與回歸樹(shù)(CART,Classification AndRegression Tree)進(jìn)行集成的方法，通過(guò)取每棵CART回歸樹(shù)葉子結(jié)點(diǎn)的均值來(lái)得到預(yù)測(cè)值。相比于使用最小化基尼系數(shù)來(lái)選擇特征的CART樹(shù)，隨機(jī)森林往往擁有更強(qiáng)的泛化能力，在處理回歸問(wèn)題上也具有較好的通用性。

　　神經(jīng)網(wǎng)絡(luò)(Neural Networks,NNs) ^[9,10,12] 也稱人工神經(jīng)網(wǎng)絡(luò)，是一種模仿大腦神經(jīng)突觸聯(lián)接的結(jié)構(gòu)進(jìn)行信息處理的數(shù)學(xué)模型。該網(wǎng)絡(luò)使用大量基本神經(jīng)元進(jìn)行計(jì)算，并能通過(guò)反饋機(jī)制來(lái)優(yōu)化網(wǎng)絡(luò)參數(shù)，擁有學(xué)習(xí)的能力。深度神經(jīng)網(wǎng)絡(luò)在擁有大規(guī)模訓(xùn)練數(shù)據(jù)的任務(wù)上表現(xiàn)突出，但由于可供訓(xùn)練的參數(shù)數(shù)量龐大，在數(shù)據(jù)量過(guò)小的任務(wù)中可能遜于普通的機(jī)器學(xué)習(xí)算法。

　　由于不同醫(yī)療箱中的藥品索取量和需求量可能不在同一個(gè)量級(jí)上，使用單一模型進(jìn)行預(yù)測(cè)并不合適。為此，本文采用閾值劃分的方式，針對(duì)數(shù)據(jù)量小于某個(gè)閾值的藥品，采用簡(jiǎn)單的嶺回歸模型進(jìn)行預(yù)測(cè)；而針對(duì)數(shù)據(jù)量大于該閾值的藥品，采用較為復(fù)雜的隨機(jī)森林回歸和神經(jīng)網(wǎng)絡(luò)共同進(jìn)行預(yù)測(cè)，并把這兩個(gè)模型輸出的均值作為最終的預(yù)測(cè)值。實(shí)驗(yàn)表明，這種方式能獲得比單一模型更高的準(zhǔn)確率。

　　2 調(diào)度路徑規(guī)劃

　　2.1 問(wèn)題描述與建模

　　物聯(lián)網(wǎng)醫(yī)療箱的藥品調(diào)度是醫(yī)療箱系統(tǒng)日常運(yùn)營(yíng)中的關(guān)鍵問(wèn)題，主要包括數(shù)量和路徑兩大難題，即需要向各醫(yī)療箱投放多少藥品以及如何規(guī)劃投放和調(diào)度的路徑。在傳統(tǒng)的自動(dòng)售藥機(jī)模式中，運(yùn)營(yíng)方需要事先預(yù)估各售藥機(jī)的藥品需求規(guī)模，并人為地對(duì)藥品的余量進(jìn)行檢查和補(bǔ)充，是一種基于經(jīng)驗(yàn)的方法。而在數(shù)據(jù)驅(qū)動(dòng)的物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)之下，由于所有的藥品需求量可以由已訓(xùn)練的網(wǎng)絡(luò)計(jì)算得出，因此可以將預(yù)測(cè)值與硬件檢測(cè)的實(shí)際藥品數(shù)量進(jìn)行比對(duì)，確定各類藥物的供需關(guān)系。

　　為此，本文假設(shè)醫(yī)療箱中的藥品數(shù)量需要滿足該藥品N日內(nèi)的需求，并定義醫(yī)療箱k中藥品的偏差量β _k,i 為當(dāng)前該藥品的實(shí)際數(shù)量與該藥品N日需求總量之差，即

其中，r _k,i 為醫(yī)療箱k中當(dāng)前藥品i的實(shí)際數(shù)量，可以由硬件檢測(cè)并上傳數(shù)據(jù)庫(kù)獲取得到；p _k,i 為預(yù)測(cè)的醫(yī)療箱k中藥品i的日平均需求量。若β _k,i 大于0，表明該醫(yī)療箱中的這類藥品處于供大于求的狀態(tài)，可以調(diào)出；若β _k,i 小于0，則表明該藥品處于緊缺狀態(tài)，希望調(diào)入。由此，醫(yī)療箱系統(tǒng)的藥物調(diào)度問(wèn)題即可轉(zhuǎn)化為一個(gè)特殊的旅行商問(wèn)題(TSP)，即派一輛小車從特定站點(diǎn)出發(fā)，訪問(wèn)各醫(yī)療箱和藥房、醫(yī)院、倉(cāng)庫(kù)等藥物供應(yīng)站點(diǎn)，尋找出總路程最短的Hamilton圈，并在這個(gè)過(guò)程中完成藥物的調(diào)度。

　　設(shè)G=(C,L)是一個(gè)有向圖，其中 C ={ c ₁ , c ₂ ,?，c _m }為m個(gè)醫(yī)療箱或供應(yīng)站點(diǎn)的集合，L={l _ij |c _i ,c _j ∈C為集合C中元素兩兩連線構(gòu)成的邊集， d _i,j ( i,j =1,2,..., m )為醫(yī)療箱 i 和j 之間 l ij 的行車距離，可以通過(guò)調(diào)用地圖軟件API接口獲取其數(shù)值。在傳統(tǒng)的TSP問(wèn)題中，所有的站點(diǎn)是沒(méi)有區(qū)分性的，即可以任意選擇站點(diǎn)作為路徑的延續(xù)。而在本問(wèn)題中，由于需要進(jìn)行藥物的調(diào)度，因此對(duì)于候選站點(diǎn)k的選擇有如下的約束條件：

其中，a _i 為當(dāng)前貨車上已有藥物i的數(shù)量，β _k,i 為醫(yī)療箱k中藥品i的偏差量，Maxload為貨車的負(fù)載量上限。當(dāng)結(jié)點(diǎn)k的偏差量大于等于0，即該藥品供大于求時(shí)，需要將其移上貨車，要求移上貨車后藥品總數(shù)量不超過(guò)貨車的最大負(fù)載量；當(dāng)結(jié)點(diǎn)k的偏差量小于0，即該藥品供不應(yīng)求時(shí)，要求貨車上有足夠多的該類藥物以補(bǔ)給醫(yī)療箱。

　　如果一個(gè)醫(yī)療箱中所有藥品都滿足上述約束條件，且該站點(diǎn)還未訪問(wèn)，則稱該站點(diǎn)為該時(shí)刻的有效候選站點(diǎn)。

　　調(diào)度路徑搜索的目標(biāo)就是不斷選擇有效候選站點(diǎn)，直到所有站點(diǎn)都已被訪問(wèn)。

　　2.2 基本蟻群算法

　　蟻群算法(Ant System,AS)是一種模擬進(jìn)化的算法，由意大利學(xué)者多里科(Marco Dorigo)于1991年提出，通過(guò)模擬螞蟻在覓食過(guò)程中釋放信息素優(yōu)化路徑的行為，構(gòu)建了一套人工的螞蟻系統(tǒng)。蟻群算法在求解TSP問(wèn)題中得到了廣泛應(yīng)用，也因此作為本文醫(yī)療箱藥物調(diào)度問(wèn)題的基本求解算法。

　　在算法的開(kāi)始，將q只螞蟻隨機(jī)地置于m個(gè)醫(yī)療箱站點(diǎn)上。對(duì)于每只螞蟻，都擁有一張屬于自己的禁忌表tabu _k ，用來(lái)表示螞蟻k已經(jīng)訪問(wèn)過(guò)的醫(yī)療箱站點(diǎn)。在本問(wèn)題中，由于候選站點(diǎn)的約束條件限制，因此需要額外設(shè)置一張無(wú)效候選站點(diǎn)表invalid _k 作為不滿足公式(2)的候選站點(diǎn)的集合。

　　在t時(shí)刻，在醫(yī)療箱 i 處的螞蟻 k 需要根據(jù)該時(shí)刻的有效候選站點(diǎn)集J _k (i)，依據(jù)某一概率函數(shù)選擇下一個(gè)站點(diǎn)j。其中，有效候選站點(diǎn)集J _k (i)為去除了禁忌表中元素和無(wú)效候選站點(diǎn)表invalid _k 中所有元素的站點(diǎn)集合，即

螞蟻從站點(diǎn)轉(zhuǎn)移到站點(diǎn)的轉(zhuǎn)移概率定義為：

否則其中， α 為信息啟發(fā)式因子，表示軌跡的相對(duì)重要性；β 為期望啟發(fā)式因子，表示能見(jiàn)度的相對(duì)重要性； η _ij (t)是啟發(fā)函數(shù)，表示螞蟻 k 從站點(diǎn) i 轉(zhuǎn)移到站點(diǎn) j 的期望程度，這里取站點(diǎn) i 和站點(diǎn) j 實(shí)際行車距離的倒數(shù)，即

式中，站點(diǎn) i 到站點(diǎn) j 實(shí)際行車距離 d _i,j 越小，則 η _ij (t) 越大，因此 η _ij (t) 可以表示為螞蟻從站點(diǎn)i轉(zhuǎn)移到站點(diǎn)j的期望程度。

　　在算法的起始階段，所有邊上的信息素量是相等的，即 τ _ij (0)=C ( C 為常數(shù))。當(dāng)所有螞蟻都尋找到一條Hamilton回路后，各路徑的信息素量根據(jù)下式來(lái)更新：

其中， ρ 為信息素的蒸發(fā)系數(shù)， Δτ _ij (t) 為本次迭代中邊 l _ij上信息素的增量， Δ 為螞蟻 k 在本次迭代中留在邊l _ij 上的信息素。在Ant-Cycle模型中，定義為：

式中， Q 為正常數(shù)，表示信息素的強(qiáng)度； L _k 為螞蟻 k 在本次迭代中所經(jīng)路徑的總長(zhǎng)度。在每一輪迭代中，所有的螞蟻都去尋找一條Hamilton回路，并更新信息素，開(kāi)始下一輪迭代。當(dāng)?shù)喆芜_(dá)到最大進(jìn)化次數(shù)時(shí)，算法結(jié)束，當(dāng)前的最優(yōu)Hamilton回路即為算法尋找到的最優(yōu)調(diào)度路徑。

　　2.3 最大最小蟻群算法

　　基本蟻群算法提供了一種求解TSP問(wèn)題的方案，但是其存在收斂速度較慢、易陷于局部最優(yōu)的缺點(diǎn)。為此，Stutzle和Hoos提出了最大-最小螞蟻系統(tǒng)(Max-MinAnt System, MMAS)。相較于基本蟻群算法(AS)，最大最小蟻群算法做了如下改進(jìn)：

　　1)采用全局信息素更新，即在每一輪迭代結(jié)束后，僅更新本輪或全局最優(yōu)解路徑上的信息素，以加強(qiáng)最優(yōu)解的影響力，即將式(6)和(7)修改為：

其中，為本次迭代中最優(yōu)路徑或全局最優(yōu)路徑上信息素的增量，其值等于最優(yōu)路徑的總長(zhǎng)度 L _best 的倒數(shù)。

　　2) 限制每條邊上的信息素在固定范圍[τ _min ,τ _max ]內(nèi)，避免某條路徑上的信息素量遠(yuǎn)大于其他路徑，造成算法過(guò)早收斂于局部最優(yōu)解。

　　3)將初始時(shí)刻各條邊上的信息素量τ _ij (0)設(shè)為τ _max ，而不再是一任意的常數(shù) C ，使算法在初始階段能擁有較高的隨機(jī)性，以提高發(fā)現(xiàn)全局最優(yōu)解的概率。

　　3 實(shí)驗(yàn)與結(jié)果分析

　　3.1 數(shù)據(jù)描述與預(yù)處理

　　本文的實(shí)驗(yàn)數(shù)據(jù)來(lái)源于正在上海大學(xué)寶山校區(qū)試運(yùn)營(yíng)的醫(yī)療箱系統(tǒng)后臺(tái)，由于用戶需要使用微信小程序進(jìn)行取藥，因而所有的取藥記錄都會(huì)保存在云服務(wù)器的MySQL數(shù)據(jù)庫(kù)中。取藥記錄包括取藥時(shí)間、醫(yī)療箱編號(hào)、藥物編號(hào)、藥物名稱等字段。本實(shí)驗(yàn)使用Python訪問(wèn)數(shù)據(jù)庫(kù)，讀取取藥記錄，并進(jìn)行了數(shù)據(jù)的處理和整合，統(tǒng)計(jì)出了各醫(yī)療箱中各類藥物的平均日索取量。

　　針對(duì)取藥日期，將其轉(zhuǎn)化為項(xiàng)目開(kāi)始運(yùn)營(yíng)到該日期的偏差天數(shù)，并加入該日期所對(duì)應(yīng)的月份和星期作為模型的輸入字段；針對(duì)醫(yī)療箱位置，調(diào)用經(jīng)緯度轉(zhuǎn)換函數(shù)將其編碼為可以用單個(gè)字段反映地理坐標(biāo)的GeoHash格式 ^[13] ；針對(duì)醫(yī)療箱所在位置的特殊地理類型，如運(yùn)動(dòng)場(chǎng)、游泳館、實(shí)驗(yàn)室、宿舍等，分別進(jìn)行One-Hot編碼將其轉(zhuǎn)化為8個(gè)字段；針對(duì)藥物評(píng)分，通過(guò)API接口調(diào)用了全球藥物評(píng)審數(shù)據(jù)庫(kù)SERMO的相關(guān)評(píng)分?jǐn)?shù)據(jù)。

　　為了更好地刻畫醫(yī)療箱的地理畫像，本實(shí)驗(yàn)調(diào)用Google地圖接口分類統(tǒng)計(jì)了各醫(yī)療箱地理坐標(biāo)附近的興趣點(diǎn)(POI)數(shù)量，包括餐飲類數(shù)量、公司企業(yè)類數(shù)量、購(gòu)物商場(chǎng)類數(shù)量、交通設(shè)施類數(shù)量、道路地名類數(shù)量等。最后，將所有上述字段和對(duì)應(yīng)的日索取量進(jìn)行歸一化處理以消除量綱的影響。

　　3.2 需求量預(yù)測(cè)

　　本文假設(shè)歷史取藥記錄中的索取量間接地反映了特定醫(yī)療箱中該藥品的需求量，希望能夠構(gòu)建一個(gè)由上述特征到藥物日平均需求量的映射關(guān)系。由于經(jīng)預(yù)處理后的特征向量維數(shù)過(guò)高，本文使用主成分分析(Principalcomponents analysis,PCA)的方法對(duì)特征向量進(jìn)行降維，并以8:2的比例劃分訓(xùn)練集和測(cè)試集。

　　本實(shí)驗(yàn)中，設(shè)定劃分閾值為 w ，即訓(xùn)練集中數(shù)據(jù)條目少于的 w 藥品，使用嶺回歸進(jìn)行建模，否則分別使用隨機(jī)森林回歸和神經(jīng)網(wǎng)絡(luò)進(jìn)行建模與預(yù)測(cè)。其中，嶺回歸和隨機(jī)森林由Python的機(jī)器學(xué)習(xí)庫(kù)Scikit-learn實(shí)現(xiàn)，隨機(jī)森林的子模型數(shù)n_estimators設(shè)為30；神經(jīng)網(wǎng)絡(luò)使用Keras框架下的反向傳播算法(BP)進(jìn)行實(shí)現(xiàn)，包括輸入層和輸出層共6層，每層的神經(jīng)元個(gè)數(shù)分別為[31,100,50,30,10,1]，激活函數(shù)使用ReLU函數(shù)，梯度優(yōu)化使用Adam算法，本文使用均方根誤差(RMSE)來(lái)衡量模型的預(yù)測(cè)精度，即

式中， p gt,i 為真實(shí)的藥品平均日需求量， p _model,i 為模型預(yù)測(cè)出的藥品日需求量。本實(shí)驗(yàn)根據(jù)上述參數(shù)設(shè)置，調(diào)整劃分閾值w，針對(duì)數(shù)據(jù)量小于w的藥品統(tǒng)一使用嶺回歸進(jìn)行建模，針對(duì)其他藥品分別使用隨機(jī)森林回歸、神經(jīng)網(wǎng)絡(luò)、隨機(jī)森林與神經(jīng)網(wǎng)絡(luò)取均值的算法進(jìn)行了實(shí)驗(yàn)，計(jì)算出了所有藥品的整體均方根誤差，如表1所示。

　　3.3蟻群算法路徑調(diào)度

　　在得到各醫(yī)療箱中各類藥品的日需求量，即可根據(jù)硬件檢測(cè)的實(shí)際數(shù)量計(jì)算出各類藥品的偏差值。實(shí)驗(yàn)中假設(shè) N =7，n =5,即每個(gè)藥物箱中共有5件藥品，要求調(diào)度結(jié)束后每件藥品能夠滿足7日的需求量。由于當(dāng)前物聯(lián)網(wǎng)醫(yī)療箱僅在上海大學(xué)寶山校區(qū)試點(diǎn)運(yùn)營(yíng)，數(shù)量不多，因此本文模擬了50個(gè)醫(yī)療箱和藥品提供商的站點(diǎn)，以檢驗(yàn)蟻群算法的調(diào)度效果。本實(shí)驗(yàn)使用Python的圖形化GUI庫(kù)Tkinter分別編寫了基本蟻群算法和最大最小蟻群算法，程序界面如圖1所示。

　　圖 中 每 個(gè) 站 點(diǎn) 上 方 的 五 元 組( β _k,1, β _k,2 ,β _k,3 ,β _k,4 ,β _k,5 )表示當(dāng)前醫(yī)療箱 k 中5件藥品的偏差值，其中黑色填充的結(jié)點(diǎn)代表貨車的起始結(jié)點(diǎn)。以使用基本蟻群算法，設(shè)置螞蟻種群數(shù)q=50，貨車最大負(fù)載量 Maxload =800，信息啟發(fā)因子 α =1.0，期望啟發(fā)因子 β =1.0，信息素強(qiáng)度 Q =100為例，得到的最優(yōu)調(diào)度總距離為4501，調(diào)度路徑如圖2所示。

　　為了比較基本蟻群算法(AS)和最大最小蟻群算法(MMAS)在該問(wèn)題上的效果，本文分別就不同的貨車最大負(fù)載量 Maxload ，分別使用兩種算法進(jìn)行了檢驗(yàn)，尋找到的最優(yōu)路徑長(zhǎng)度如表2所示。

　　可見(jiàn)，最大最小蟻群算法(MMAS)往往能夠得到比基本蟻群算法(AS)更優(yōu)的解；此外，貨車最大負(fù)載量Maxload 也會(huì)在一定程度上影響到最優(yōu)調(diào)度路徑。

　　4 結(jié)論

　　本文針對(duì)一種新型物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)，提出了一套使用機(jī)器學(xué)習(xí)和蟻群算法的藥物調(diào)度方法，使用嶺回歸、隨機(jī)森林回歸和神經(jīng)網(wǎng)絡(luò)的混合算法建立了藥品需求量的預(yù)測(cè)模型，并使用基本蟻群算法和最大最小蟻群算法對(duì)藥品的調(diào)度問(wèn)題進(jìn)行了求解。該方法為物聯(lián)網(wǎng)醫(yī)療箱系統(tǒng)提供了一種有效的解決方案，簡(jiǎn)化了日常運(yùn)營(yíng)維護(hù)的過(guò)程。此外，在其他的調(diào)度問(wèn)題中，該方法也提供了一種可參考的解決思路。

　　致謝

　　最后感謝上海大學(xué)機(jī)電工程與自動(dòng)化學(xué)院“挑戰(zhàn)杯”大學(xué)生課外學(xué)術(shù)科技作品競(jìng)賽項(xiàng)目對(duì)本文的研究所提供的支持。以及上海大學(xué)“羅姆杯”大學(xué)生機(jī)電工程創(chuàng)新設(shè)計(jì)大賽對(duì)本文研究的支持。

　　參考文獻(xiàn)：

　　[1] 趙安琪.自動(dòng)售藥機(jī)破冰之路[J].中國(guó)藥店,2010(10):32-32.

　　[2] 陳雷毅,潘榮,許強(qiáng),等.高校校園自動(dòng)售藥機(jī)的再設(shè)計(jì)[J].大眾文藝,2016(11).

　　[3] Zhao G, Luo W, Sun R, et al. A Modified Max-Min Ant System for Vehicle RoutingProblems[C]// International Conference on Wireless Communications. IEEE, 2008.

　　[4]段海濱.蟻群算法原理及其應(yīng)用[M]. 北京 科學(xué)出版社,2005.

　　[5] Stützle T,Hoos H H.MAX–MIN ant system[J]. Future generation computer systems, 2000,16(8): 889-914.

　　[6] 陳昌敏,謝維成,范頌頌.自適應(yīng)和最大最小蟻群算法的物流車輛路徑優(yōu)化比較[J].西華大學(xué)學(xué)報(bào)（自然科學(xué)版）,2011,30(3).

　　[7] 李勇,段正澄.動(dòng)態(tài)蟻群算法求解TSP問(wèn)題[J].計(jì)算機(jī)工程與應(yīng)用,2003,39(17):103-106.

　　[8] 柴寶杰,劉大為.基于粒子群優(yōu)化的蟻群算法在TSP中的應(yīng)用[J].計(jì)算機(jī)仿真,2009,26(8):89-91.

　　[9] Peter Harrington. Machine learning in action [M].2013.

　　[10] Liaw A, Wiener M. Classification and regression by randomForest[J].R news,2002,2(3):18-22.

　　[11] Bose N K, Liang P.Neural Network Fundamentals with Graphs,Algorithms,and Applications(McGraw-Hill Series in Electrical Computer Engineering)[J].1996.

　　[12] Balki? Z, ?o?tari? D, Horvat G. GeoHash and UUID identifier for multi-agent systems[C]//KES International Symposium on Agent and Multi-Agent Systems: Technologies andApplications. Springer, Berlin, Heidelberg, 2012: 290-298.

　　[13] Cassel M, Lima F.Evaluating one-hot encoding finite state machines for SEU reliabilityin SRAM-based FPGAs[C]//On-Line Testing Symposium, 2006. IOLTS 2006. 12th IEEEInternational. IEEE, 2006: 6 pp.

　　作者簡(jiǎn)介:

　　通信作者：李晉，男，上海大學(xué)機(jī)電工程學(xué)院實(shí)驗(yàn)師，碩士，主要從事微控制器技術(shù)，人工智能算法應(yīng)用優(yōu)化。

　　第一作者：盧嘉軒，男，上海大學(xué)計(jì)算機(jī)學(xué)院本科生，主要從事人工智能算法研究。

　　第三作者：周延，男，上海大學(xué)計(jì)算機(jī)學(xué)院本科生，主要從事人工智能算法研究。

　　本文來(lái)源于科技期刊《電子產(chǎn)品世界》2019年第7期第80頁(yè)，歡迎您寫論文時(shí)引用，并注明出處