汪儉華，馮錫煒，張永攀，關(guān) 滄，畢福偉

（1.遼寧石油化工大學(xué) 計(jì)算機(jī)與通信工程學(xué)院，遼寧 撫順 113001；2.撫順石化公司石油三廠(chǎng) 遼寧 撫順 113001）

隨著現(xiàn)代科技的發(fā)展，蒙特卡羅模擬算法因?yàn)槟軌驅(qū)⒉淮_定性隨機(jī)事件以概率的形式表示出來(lái)，以這種事件出現(xiàn)的頻度估計(jì)這一隨機(jī)事件的概率來(lái)求得近似解，較好地解決了各種不確定性問(wèn)題，從而在物理科學(xué)、機(jī)械制造、金融計(jì)算、信息科學(xué)以及計(jì)算機(jī)科學(xué)等多個(gè)領(lǐng)域得到廣泛的重視和應(yīng)用[1]。蒙特卡羅模擬算法的主要思想為：當(dāng)所求解問(wèn)題是某種隨機(jī)事件出現(xiàn)的概率，或者是某個(gè)隨機(jī)事件的期望值時(shí)，通過(guò)某種“實(shí)驗(yàn)”的方法，以這種事件出現(xiàn)的頻率估計(jì)這一隨機(jī)事件的概率，或者得到這個(gè)隨機(jī)事件的某些數(shù)字特征，并將其作為問(wèn)題的解。

蒙特卡羅模擬分析結(jié)果的精度取決于采樣次數(shù)和計(jì)算次數(shù)，隨著采樣次數(shù)和計(jì)算次數(shù)的增加，得到的結(jié)果也就越準(zhǔn)確，但是相應(yīng)的計(jì)算速度就會(huì)降低，尤其是在大規(guī)模采樣下的計(jì)算量很大時(shí)，得到最終結(jié)果的計(jì)算速度會(huì)更慢。為了獲得精確、高效的結(jié)果，必須設(shè)計(jì)一個(gè)具有高性能計(jì)算能力的平臺(tái)來(lái)支撐。

本文以建立在Platform Symphony之上的云計(jì)算平臺(tái)為載體，通過(guò)Java編程對(duì)平臺(tái)功能進(jìn)行擴(kuò)展，并將Excel和COM組件與云計(jì)算平臺(tái)進(jìn)行集成，實(shí)現(xiàn)了基于云計(jì)算的蒙特卡羅模擬分析。

1 相關(guān)技術(shù)介紹

1.1 云計(jì)算

云計(jì)算是一種新興的網(wǎng)絡(luò)計(jì)算技術(shù)，它提供了通過(guò)網(wǎng)絡(luò)訪(fǎng)問(wèn)、可動(dòng)態(tài)配置的計(jì)算資源池。云計(jì)算涵蓋的范圍很廣泛，內(nèi)容也很豐富，通常將其分為以下3層[2-3]：

IaaS（Infrastructure as a Service，基礎(chǔ)設(shè)施即服務(wù)）。 將服務(wù)器、網(wǎng)絡(luò)設(shè)備、存儲(chǔ)設(shè)備等硬件資源集成為資源池供用戶(hù)按需使用。

PaaS（Platform as a Service，平臺(tái)即服務(wù)）。為用戶(hù)提供一個(gè)完整的平臺(tái)，它包含應(yīng)用開(kāi)發(fā)、接口開(kāi)發(fā)、數(shù)據(jù)庫(kù)開(kāi)發(fā)等，保證信息平臺(tái)在實(shí)現(xiàn)技術(shù)方面的擴(kuò)展性和兼容性。

SaaS（Software as a Service，軟件即服務(wù)）。為用戶(hù)提供基于網(wǎng)絡(luò)的應(yīng)用軟件服務(wù)。

三層結(jié)構(gòu)在真實(shí)的云計(jì)算環(huán)境中是可以相對(duì)獨(dú)立存在的，可以逐層搭建云系統(tǒng)，也可以根據(jù)需要只搭建其中的一層或幾層。

1.2 Platform Sym phony

Platform Symphony是Platform Computing公司開(kāi)發(fā)的面向服務(wù)架構(gòu)的云計(jì)算平臺(tái)，用戶(hù)可以使用其提供的編程接口擴(kuò)展程序功能，從而得到個(gè)性化的服務(wù)[4]。Symphony基于主從架構(gòu)，由主服務(wù)器將任務(wù)發(fā)送給從服務(wù)器執(zhí)行計(jì)算，其架構(gòu)如圖1所示。

圖1 Symphony架構(gòu)圖Fig.1 Structure diagram of the symphony

SD（Session Director，服務(wù)注冊(cè)中心）。運(yùn)行在主服務(wù)器上，對(duì)通過(guò)編程接口開(kāi)發(fā)的應(yīng)用程序進(jìn)行注冊(cè)之后供客戶(hù)端使用。同時(shí)為各個(gè)子服務(wù)提供任務(wù)分配者模塊SSM，并與客戶(hù)端進(jìn)行連接。

SSM（Symphony Session Manager，任務(wù)分配者）。運(yùn)行在主服務(wù)器上，負(fù)責(zé)接受客戶(hù)端提交的任務(wù)，并根據(jù)相應(yīng)的調(diào)度配置策略（可以自行配置）將任務(wù)發(fā)送到合適的從服務(wù)器上執(zhí)行。

SIM（Service Instance Manager，任務(wù)執(zhí)行者）。可以運(yùn)行在主服務(wù)器與從服務(wù)器上，負(fù)責(zé)執(zhí)行具體的任務(wù)。它從SSM獲取需要執(zhí)行的任務(wù)，并將執(zhí)行結(jié)果返回給SSM，SSM接受后統(tǒng)一發(fā)送到客戶(hù)端。

資 源 管 理 模 塊 （Resource Orchestrator）。 負(fù) 責(zé) 管 理Symphony中的所有計(jì)算資源。

編程接口。供開(kāi)發(fā)人員開(kāi)發(fā)應(yīng)用程序來(lái)擴(kuò)展平臺(tái)的功能。

1.3 Excel

微軟出品的Excel是一款優(yōu)秀的數(shù)值計(jì)算軟件，功能強(qiáng)大、操作簡(jiǎn)單，可以用于數(shù)據(jù)處理、統(tǒng)計(jì)分析和輔助決策等操作。它還支持創(chuàng)建Excel宏命令以及VBA（Visual Basic For Application）編程來(lái)擴(kuò)展自身的功能，Excel經(jīng)過(guò)功能擴(kuò)展之后可以與云計(jì)算平臺(tái)進(jìn)行集成。集成的流程如下：1）創(chuàng)建Excel宏命令；2）使用 VBE（Visual Basic Editor）編程實(shí)現(xiàn)基本程序功能，然后通過(guò)API接口進(jìn)行功能擴(kuò)展；3）調(diào)用COM組件將Excel與云計(jì)算平臺(tái)集成。

1.4 COM組件

COM組件是微軟為了軟件的研發(fā)更適合人類(lèi)的行為邏輯而推出的軟件開(kāi)發(fā)技術(shù)[5]。在COM架構(gòu)下，開(kāi)發(fā)人員可以開(kāi)發(fā)出各種功能專(zhuān)一的組件，然后將其按照需要組合構(gòu)成復(fù)雜的應(yīng)用。它有以下優(yōu)點(diǎn)：1）可以隨時(shí)替換現(xiàn)有的組件，方便進(jìn)行功能擴(kuò)展；2）同一個(gè)組件可以在多個(gè)應(yīng)用中復(fù)用，利于模塊化開(kāi)發(fā)；3）可以將應(yīng)用部署到網(wǎng)絡(luò)環(huán)境。

2 基于云計(jì)算的蒙特卡羅模擬平臺(tái)

2.1 平臺(tái)的設(shè)計(jì)

基于云計(jì)算的蒙特卡羅模擬平臺(tái)運(yùn)用現(xiàn)代化的計(jì)算機(jī)技術(shù)以突破時(shí)間、空間的制約，通過(guò)精簡(jiǎn)、優(yōu)化、整合計(jì)算資源來(lái)進(jìn)行蒙特卡羅模擬分析，使所得的結(jié)果更加的精確、高效。根據(jù)蒙特卡羅模擬算法和云計(jì)算技術(shù)特點(diǎn)[6]，其結(jié)構(gòu)如圖2所示。

基礎(chǔ)設(shè)施：由計(jì)算節(jié)點(diǎn)、存儲(chǔ)設(shè)備、網(wǎng)絡(luò)設(shè)施以及基礎(chǔ)軟件組成的計(jì)算集群資源池。利用廉價(jià)的設(shè)備動(dòng)態(tài)擴(kuò)展資源池，以較低投資獲得高性能計(jì)算的能力。

軟件應(yīng)用：將蒙特卡羅模擬算法轉(zhuǎn)化為應(yīng)用服務(wù)，提供給客戶(hù)端使用。

2.2 平臺(tái)的實(shí)現(xiàn)

2.2.1 搭建計(jì)算節(jié)點(diǎn)集群

集群使用4個(gè)計(jì)算節(jié)點(diǎn)，每個(gè)節(jié)點(diǎn)配置相同，CPU是1.83 GHz，內(nèi)存2 G，千兆以太網(wǎng)卡，操作系統(tǒng)為Windows Server 2008，Java環(huán)境為 JDK 1.7，平臺(tái)軟件為 Symphony DE 5.0。

1）雙擊Windows平臺(tái)的安裝包SymphonyDE-5.0.0_win32.msi進(jìn)行安裝；

2）選取安裝了Symphony DE的主機(jī)Jeffw1為主服務(wù)器，對(duì)DE50confvem_resource.conf文件進(jìn)行以下配置，然后將經(jīng)過(guò)配置的文件復(fù)制到所有的從服務(wù)器相同目錄下；

3）關(guān)閉所有的服務(wù)器后進(jìn)行重啟，注意需要先啟動(dòng)主服務(wù)器。

2.2.2 擴(kuò)展平臺(tái)功能

在現(xiàn)代市場(chǎng)經(jīng)濟(jì)條件下，住房商品化、市場(chǎng)化引致住房金融的產(chǎn)生和繁榮。一方面，商品住房建設(shè)需要建設(shè)貸款，商品住房開(kāi)發(fā)需要開(kāi)發(fā)貸款，商品住房消費(fèi)需要租賃貸款和按揭貸款等，因此，住房金融的產(chǎn)生和繁榮反過(guò)來(lái)能夠促進(jìn)住房市場(chǎng)繁榮。但是，另一方面，住房金融也會(huì)因?yàn)橐宰》可唐窞槊浇榈耐顿Y和投機(jī)，容易衍生住房?jī)r(jià)格泡沫和金融泡沫，最終演化成住房危機(jī)和金融危機(jī)。因此，住房金融是本文研究的重點(diǎn)內(nèi)容。

Symphony DE支持C/C++、Java，C#等多種程序語(yǔ)言進(jìn)行功能擴(kuò)展。采用Java語(yǔ)言編寫(xiě)程序，在Eclipse中利用Ant工具編譯并生成ProjectAnalysisService程序包，然后將其部署到主服務(wù)器的Symphony DE中對(duì)平臺(tái)功能進(jìn)行擴(kuò)展。右擊Symphony運(yùn)行圖標(biāo)，選擇Add or Remove Application。在添加應(yīng)用向?qū)е刑砑优渲梦募rojectAnalysis.xm l，并選擇程序包ProjectAnalysisPackage.jar，單擊下一步，直到確認(rèn)即可。部署成功后，可以在管理界面看到如圖3所示內(nèi)容。

圖2 基于云計(jì)算的平臺(tái)結(jié)構(gòu)圖Fig.2 Structure diagram of the platform based on cloud

圖3 部署程序包到Symphony平臺(tái)圖Fig.3 Chart of deploy program package to the symphony platform

2.2.3 Excel與平臺(tái)集成

在Excel中調(diào)用Symphony DE的COM組件與平臺(tái)進(jìn)行集成，從Excel中啟動(dòng)VBE，在VBE中單擊“工具”→“引用”，找到DE50confwin32-vc7libCOMPlatform.Symphony.Soam.COM.dll文件，點(diǎn)擊確定。集成架構(gòu)由Excel客戶(hù)端、Symphony DE云計(jì)算控制器以及平臺(tái)上所運(yùn)行的應(yīng)用服務(wù)這3個(gè)部分組成，其工作流程如圖4所示。

2.2.4 實(shí)現(xiàn)過(guò)程中的關(guān)鍵技術(shù)

1）在Symphony DE中擴(kuò)展云計(jì)算平臺(tái)

部署的Project Analysis Service程序，通過(guò)API接口連接到Symphony DE，其代碼如下：

2）Excel中調(diào)用云計(jì)算環(huán)境

在Excel中執(zhí)行云計(jì)算必須連接到Symphony DE服務(wù)器，通過(guò)初始化SOA中間件以及獲取服務(wù)來(lái)完成云計(jì)算環(huán)境的初始化，其代碼如下：

3 實(shí)驗(yàn)與分析

實(shí)例描述：某信息系統(tǒng)集成項(xiàng)目計(jì)劃投資1 000萬(wàn)元，該項(xiàng)目預(yù)期年銷(xiāo)售收入為650萬(wàn)元，年銷(xiāo)售稅金為20萬(wàn)元，年經(jīng)營(yíng)成本為320萬(wàn)元，年折舊費(fèi)50萬(wàn)元，收入與成本的波動(dòng)變化幅度都為20%，分別在單機(jī)環(huán)境和云計(jì)算平臺(tái)上用蒙特卡羅模擬算法計(jì)算收益率大于等于20%的概率。

在云計(jì)算平臺(tái)上進(jìn)行蒙特卡羅模擬分析的應(yīng)用，按照如下步驟進(jìn)行：1）配置Symphony DE上的計(jì)算資源；2）初始化任務(wù)數(shù)據(jù)；3）在Excel中創(chuàng)建本地分析和云計(jì)算分析；4）調(diào)用自定義的分析策略；5）啟動(dòng)分析；6）在分析結(jié)束后統(tǒng)計(jì)結(jié)果。建立的實(shí)驗(yàn)?zāi)Ｐ停鐖D5所示。

實(shí)驗(yàn)1：?jiǎn)螜C(jī)計(jì)算節(jié)點(diǎn)為本地主機(jī)，云計(jì)算節(jié)點(diǎn)由1個(gè)主服務(wù)器，3個(gè)從服務(wù)器組成。將計(jì)算次數(shù)設(shè)置為10次，采樣次數(shù)從10 000到10 000 000變化，觀(guān)察單機(jī)計(jì)算與云計(jì)算下運(yùn)行的情況。

實(shí)驗(yàn)2：?jiǎn)螜C(jī)計(jì)算節(jié)點(diǎn)為本地主機(jī)，云計(jì)算節(jié)點(diǎn)由1個(gè)主服務(wù)器，3個(gè)從服務(wù)器組成。將采樣次數(shù)設(shè)置為10 000次，計(jì)算次數(shù)從10到400變化，觀(guān)察單機(jī)計(jì)算與云計(jì)算下運(yùn)行的情況。

表1 實(shí)驗(yàn)1結(jié)果Tab.1 Resu lt of experiment 1

實(shí)驗(yàn)3：云計(jì)算節(jié)點(diǎn)數(shù)由1個(gè)到4個(gè)變化。設(shè)采樣次數(shù)置為10 000次，計(jì)算次數(shù)置為100次，觀(guān)察云計(jì)算下運(yùn)行的情況。

圖4 工作流程圖Fig.4 Flow chart of the work

圖5 實(shí)驗(yàn)?zāi)Ｐ蛨DFig.5 Chart of experimental model

表2 實(shí)驗(yàn)2結(jié)果Tab.2 Result of experiment 2

圖6 實(shí)驗(yàn)1時(shí)間趨勢(shì)圖Fig.6 Time trend diagram of experiment 1

圖7 實(shí)驗(yàn)2時(shí)間趨勢(shì)圖Fig.7 Time trend diagram of experiment 2

圖8 實(shí)驗(yàn)3時(shí)間趨勢(shì)圖Fig.8 Time trend diagram of experiment 3

表1顯示了在單位計(jì)算次數(shù)下，隨著采樣次數(shù)的增加，在單機(jī)計(jì)算與云計(jì)算情況下的運(yùn)行數(shù)據(jù)。表2顯示了單位采樣次數(shù)下，隨著計(jì)算次數(shù)的增加，在單機(jī)計(jì)算與云計(jì)算下的運(yùn)行數(shù)據(jù)。隨著采樣次數(shù)與計(jì)算次數(shù)的增加，所獲的數(shù)據(jù)精度越高，但所需時(shí)間呈上升趨勢(shì)，云計(jì)算相對(duì)本機(jī)計(jì)算所需的時(shí)間大幅減少，得到的結(jié)果更加準(zhǔn)確，在大規(guī)模計(jì)算時(shí)的優(yōu)勢(shì)更為明顯，消耗的時(shí)間趨勢(shì)分別如圖6、圖7所示。圖8為不同數(shù)量的云計(jì)算節(jié)點(diǎn)進(jìn)行蒙特卡羅模擬分析所消耗的時(shí)間趨勢(shì)。隨著計(jì)算節(jié)點(diǎn)數(shù)的增加，所需的時(shí)間呈指數(shù)下降趨勢(shì)。

4 結(jié) 論

云計(jì)算作為IT技術(shù)的第四次革命，將給大規(guī)模計(jì)算帶來(lái)新的機(jī)遇。本文充分應(yīng)用現(xiàn)有的計(jì)算資源搭建云計(jì)算環(huán)境，通過(guò)Excel與云計(jì)算的集成，既解決了大規(guī)模計(jì)算效率的問(wèn)題，也將云計(jì)算技術(shù)由企業(yè)和組織推向了普通使用者。實(shí)驗(yàn)表明，用云計(jì)算技術(shù)來(lái)實(shí)現(xiàn)蒙特卡羅模擬分析，所得結(jié)果精確、高效，具有不可比擬的優(yōu)越性和實(shí)用價(jià)值。

[1]陳鄧安，叢中旖，王子明.在效能分析中應(yīng)用Excel數(shù)據(jù)分析模塊進(jìn)行蒙特卡羅模擬[J].科學(xué)技術(shù)與工程，2012，12（2）:397-401.CHEN Deng-an，CONG Zhong-yi，WANG Zi-ming.Research on monte carlo simulation method in effectiveness analysis by data analysis modules in Excel[J].Science Technology and Engineering，2012，12（2）:397-401.

[2]趙廣才，張雪萍.云計(jì)算技術(shù)分析及其展望[J].電子設(shè)計(jì)工程，2011，19（22）:4-10.ZHAO Guang-cai，ZHANG Xue-ping.The analysis of cloud computing technology and its prospect[J].Electronic Design Engineering，2011，19（22）:4-10.

[3]陳康，鄭緯民.云計(jì)算:系統(tǒng)實(shí)例與研究現(xiàn)狀[J].軟件學(xué)報(bào)，2009，20（5）:1337-1348.CHEN Kang，ZHENG Wei-min.Cloud computing:system instance and research[J].Journal of Software，2009，20 （5）:1337-1348.

[4]Platform Computing.Symphony使用文檔[EB/OL].（2011-09-20）/[2012-5-20].http://www.platform.com/workload-manage ment/SOA-grid-computing/.

[5]喬磊，馬火林，陳海紅，等.基于COM組件計(jì)算巖層真厚度的軟件設(shè)計(jì)[J].工程地球物理學(xué)報(bào)，2010，7（3）:368-371.QIAO Lei，MA Huo-lin，CHEN Hai-hon，et al.The software designed for true thickness of the layer based on COM components[J].Chinese Journal of Engineering Geophysics，2010，7（3）:368-371.

[6]鄭廣成.一種SOA云服務(wù)平臺(tái)架構(gòu)研究與應(yīng)用 [J].計(jì)算機(jī)應(yīng)用與軟件，2011，28（12）:228-231.ZHENG Guang-cheng.An soa cloud service platform architecture research and application[J].Computer Applications and Software，2011，28（12）:228-231.

[7]HPC Community.About Symphony DE. （2010-08-11）/[2012-05-20].http://www.hpccommunity.org/content/about-symphony-de-145/.