


版權(quán)說(shuō)明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)
文檔簡(jiǎn)介
1、專業(yè):核工程與核技術(shù)實(shí)驗(yàn)一蒙特卡羅方法一、實(shí)驗(yàn)?zāi)康?、了解蒙特卡羅方法方法的基本思想;2、掌握蒙特卡羅方法計(jì)算面積、體積的方法;3、掌握由已知分布的隨機(jī)抽樣方法。二、實(shí)驗(yàn)原理Monte Carlo方法,又稱統(tǒng)計(jì)模擬方法或計(jì)算機(jī)隨機(jī)模擬方法,是一種基于“隨機(jī)數(shù)” 進(jìn)行數(shù)值模擬的方法,一種采用統(tǒng)計(jì)抽樣理論近似求解物理或數(shù)學(xué)問(wèn)題的方法。倘若待求量可以表述成某些特征量的期望值、某些事件出現(xiàn)的概率或兩者的函數(shù)形式, 那么可采用蒙特卡羅方法求解。在求解某些特征量的期望值或某些事件出現(xiàn)的概率時(shí),必 須構(gòu)建合符實(shí)際的數(shù)學(xué)模型。例如采用蒙特卡羅方法計(jì)算某函數(shù)所圍面積時(shí),構(gòu)建的數(shù)學(xué) 模型是構(gòu)造一已知面積的可均勻抽
2、樣區(qū)域,在該區(qū)域投點(diǎn),由伯努利定理大數(shù)定理可知, 進(jìn)入待求區(qū)域投點(diǎn)的頻率依概率 1收斂于該事件出現(xiàn)的概率(面積之比)。由已知分布的隨機(jī)抽樣方法指的是由已知分布的總體中抽取簡(jiǎn)單子樣。具體方法很多,詳見(jiàn)課本第三章。三、實(shí)驗(yàn)內(nèi)容1安裝所需計(jì)算工具(MATLAB等);以下內(nèi)容采用工具軟件中自帶偽隨機(jī)數(shù)發(fā)生器進(jìn)行計(jì)算。2、求解以下區(qū)域的面積、體積:2.1、 給定曲線y =2 - x2和曲線y3 = x2,曲線的交點(diǎn)為:Pi( - 1, 1 )、P2( 1, 1 )。曲線圍成平面有限區(qū)域,用蒙特卡羅方法計(jì)算區(qū)域面積;!22.2、 計(jì)算Z X y所圍體積z 1 Ji x2 y2其中 (x,y,z)| 1 x
3、 1, 1 y 1,0 z 2。3、對(duì)以下已知分布進(jìn)行隨機(jī)抽樣:三、實(shí)驗(yàn)報(bào)告編寫(xiě)1、給出各題的抽樣程序并解釋語(yǔ)句的含義;2、給出2.1和2.2抽樣結(jié)果誤差隨抽樣次數(shù)的關(guān)系圖,并解釋原因;表1實(shí)驗(yàn)記錄表3、給出3題的抽樣框圖、試驗(yàn)累積頻率與理論累積頻率關(guān)系圖,并給出抽樣次數(shù)(106) 與抽樣時(shí)間。2.1程序代碼編寫(xiě)如下:N=10A6;%總投點(diǎn)個(gè)數(shù)S=0; %記錄投點(diǎn)在所圍圖形中的個(gè)數(shù)ss=o;for i=1:Nx=2*rand-1;%產(chǎn)生的隨機(jī)變量x,yy=2*rand; ;%產(chǎn)生x和y的坐標(biāo)if(yv=2-xA2)&(yA3>=xA2)%判定是否落入所圍圖像中S=S+1; %進(jìn)入
4、則加1SS=SS+1A2;endendArea=4*S/N%計(jì)算面積Dev=SS/N-(S/NF2% 計(jì)算方差A(yù)=sqrt(Dev/N)%計(jì)算標(biāo)準(zhǔn)差toc實(shí)驗(yàn)數(shù)據(jù)如下:請(qǐng)輸入總投點(diǎn)個(gè)數(shù)150000-0.8-0.6-0.4-0.22.2實(shí)驗(yàn)代碼如下:clear;M=0;N= 5*10A4 ;for i=1:Nx=2*ra nd()-1;y=2*ra nd()-1;z=2*ra nd();t=xA2+yA2;s=zA2;if s>=tif t<=-s+2*zM=M+1;endendendtocMIANJI=M/N*8clear M N i x y;計(jì)算結(jié)果:N=50000時(shí)面積為3.135Q計(jì)算時(shí)間約0.282s,實(shí)驗(yàn)數(shù)據(jù)如下:2.3程序代碼編寫(xiě)如下:clear;clc;M = in put('輸入所需產(chǎn)生隨機(jī)變量的個(gè)數(shù):n');x = zeros(M,1);tic;for i=1:Mif(ran d()<=0.5)x(i) = max(ra nd(),ra nd()x(i) = max(x(i),ra nd();x(i) = max(x(i),ra nd();elsex(i) = min(ran d(),ra nd();x(i) = min( x(i),ra nd();endend plot(x,'g.')toe clear
溫馨提示
- 1. 本站所有資源如無(wú)特殊說(shuō)明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
- 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
- 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒(méi)有圖紙預(yù)覽就沒(méi)有圖紙。
- 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
- 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
- 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
- 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。
最新文檔
- 2025年電動(dòng)汽車(chē)技術(shù)考試試卷及答案
- 2025年城市公共設(shè)施管理知識(shí)測(cè)試卷及答案
- 物業(yè)客服個(gè)人工作總結(jié)
- 宴會(huì)廳音響設(shè)備調(diào)試標(biāo)準(zhǔn)
- 全膀胱切除術(shù)后造口護(hù)理
- 電力系統(tǒng)自動(dòng)化知識(shí)點(diǎn)習(xí)題集匯編
- 我身邊的英雄人物描寫(xiě)13篇范文
- 2025年江西汽車(chē)駕駛員客運(yùn)資格證考試題及答案
- 網(wǎng)絡(luò)平臺(tái)內(nèi)容服務(wù)提供合同
- 旅游行業(yè)旅游資源管理練習(xí)題
- 保潔學(xué)校管理制度
- 招聘渠道ROI評(píng)估模型-洞察及研究
- 2025春季學(xué)期國(guó)開(kāi)電大本科《人文英語(yǔ)4》一平臺(tái)機(jī)考真題及答案(第六套)
- 第七單元1認(rèn)識(shí)小數(shù)(課件)-三年級(jí)數(shù)學(xué)下冊(cè)(人教版)
- 2025年河北省中考麒麟卷生物(二)及答案
- 2024年民族出版社招聘事業(yè)編制專業(yè)技術(shù)人員真題
- 2025年食品安全管理員考試試題及答案
- 2025-2030骨科植入器材產(chǎn)業(yè)市場(chǎng)深度分析及發(fā)展趨勢(shì)與投資戰(zhàn)略研究報(bào)告
- T/SHPTA 071.1-2023高壓電纜附件用橡膠材料第1部分:絕緣橡膠材料
- 湖北省浠水縣聯(lián)考2025年七下數(shù)學(xué)期末質(zhì)量檢測(cè)試題含解析
- 生產(chǎn)基層管理培訓(xùn)課程
評(píng)論
0/150
提交評(píng)論