資源描述：
《蒙特卡羅法(Monte Carlo method)》由會員上傳分享，免費在線閱讀，更多相關(guān)內(nèi)容在行業(yè)資料-天天文庫。
1、實驗十二計算機仿真實驗實驗目的：1.掌握全概率公式與貝葉斯公式；2.了解計算機仿真方法；3.了解蒙特卡羅法(MonteCarlomethod),具有初級編程能力.實驗原理：全概率公式:設A1,A2,…,An為兩兩互斥事件，B是A1+A2+…+An的子事件，則P(B)=P(A1)P(B
2、A1)+P(A2)P(B
3、A2)+…+P(An)P(B
4、An).貝葉斯公式:P(Ak
5、B)=P(Ak)P(B
6、Ak)/P(B).計算機仿真:就是在計算機上模擬各種實際系統(tǒng)的運行過程.計算機仿真通常用來產(chǎn)生規(guī)定分布的隨機變量.對于任
7、意隨機變量x，其分布函數(shù)為F(x)，設h=F(x)的分布函數(shù)為G(y)，則G(y)=P{h≤y}=P{F(x)≤y}=P{x≤F-1(y)}=y，這說明h服從[0，1]的均勻分布.一般的編程語言都提供了均勻分布隨機數(shù)發(fā)生器.應用隨機數(shù)模擬試驗的方法通常稱為蒙特卡羅法(MonteCarlomethod).蒙特卡羅法不僅適用于處理隨機性問題,如存貯、排隊、質(zhì)量檢驗、市場營銷、社會救急、生態(tài)競爭和傳染病等問題；也可處理定性問題,如計算多重積分、解積分方程及微分方程、解整數(shù)規(guī)劃(特別是非線性整數(shù)規(guī)劃)等.應用蒙特卡羅法
8、解規(guī)劃問題的基本思想是：先估計各個變量的大致取值范圍，每次試驗從中隨機取出一個樣本點，然后判斷它是否為可行點.若是則將其目標函數(shù)值與上一次的目標函數(shù)值相比較，記錄下較優(yōu)目標函數(shù)值與其樣本點；否則重新抽樣。直到試驗次數(shù)達到指定值或可行點數(shù)達到指定值為止.實驗內(nèi)容：1.設有兩個口袋，甲袋中盛有兩個白球，一個黑球，乙袋中盛有一個白球，兩個黑球.由甲袋任取一個球放入乙袋，再從乙袋中取出一個球.若從乙袋中取出的球是白球，那么從甲袋中取出放入乙袋的球是白球還是黑球?用計算機模擬上述過程1000次，問理論判斷是正確的有多少次
9、?2.用計算機模擬隨機變量x～j(x)=的取值200次.3.用計算機模擬隨機變量x～N(120,102)的取值800次，并畫出統(tǒng)計直方圖.4.應用蒙特卡羅法解非線性規(guī)劃問題:maxz=-2x2-y2+xy+8x+3ys.t.3x+y=10x≥0,y≥0.C語言簡介(僅介紹本實驗所用到的)1.標識符標識符是由程序員定義的單詞，如函數(shù)名、變量名等.標識符是由大小寫字母、數(shù)字和下劃線組成的，并以字母和下劃線開始.2.關(guān)鍵字void(無值型)char(字符型)int(整數(shù)型)long(長整數(shù)型)float(浮點型)do
10、uble(雙精度浮點型)if(如果)else(或者)for(循環(huán))while(循環(huán))break(滿足一定的條件終止循環(huán))return(返回函數(shù)值)3.函數(shù)形式類型函數(shù)名(參數(shù)){}4.庫函數(shù)簡介doublesin(doublex)doublecos(doublex)doubleexp(doublex)doublelog(doublex)doublesqrt(doublex)doublepow(doublex,doubley)doublefabs(doublex)intabs(intx)intprintf(con
11、stchar*format,...);屏幕格式化輸出函數(shù)FILE*fopen(constchar*filename,constchar*mode);文件打開intfprintf(FILE*fp,constchar*format,...);文件格式化輸出函數(shù)intfgetc(FILE*fp);從文件中讀出一字符intfclose(FILE*fp);文件關(guān)閉voidfarsetcolor(intcolor);設置輸出顏色voidfarline(intx1,inty1,intx2,inty2);畫直線voidfarr
12、ectangle(intleft,inttop,intright,intbottom);畫矩形unsignedfargetpixel(intx,inty);讀出點(x,y)的顏色voidfarputpixel(intx,inty,intpixelcolor);畫點intrandom(intNum),均勻產(chǎn)生0到Num-1中的一個隨機數(shù)5.示例計算9！#includevoidmain(void){inti;longn;n=1;for(i=1;i<=9;i++)n*=i;printf("9!=
13、%ld",n);getch();}6.TurboC(2.0)編輯命令F3錄入文件F9編譯Ctrl+F9運行Ctrl+KB定義塊首Ctrl+KK定義塊尾Ctrl+KC塊粘貼Ctrl+KV塊移動Ctrl+Y刪除當前行7.部分源程序程序LAB1_1.C求出方程sinx-x=1在(-2，2)內(nèi)的近似根#include#includefloatf(float
