R語言 幾何概型計算積分 蒙特卡羅方法

2021-10-22 14:00:07 字數 775 閱讀 1500

設計幾何概型實驗,估算積分的值

**詳細解釋在**注釋中

# 定義乙個函式,該函式為所求的積分函式

# 因此在計算上述定積分和計算sqrt(2)同理

# 在此用例為sqrt(2)

f = function(x)

sqrt(2)

# 生成x

x = seq(1.5, 3.5, length=100)

y = rep(0, length(x))

j = 1

# 計算每個x對應的y值

for (i in x)

# 根據函式劃出積分曲線

plot(x, y, type='o')

# 確定積分邊界

abline(v = 2)

abline(v = 3)

# 簡單設定隨機點 (x,y)|2<=x<=3,0<=y<=100

# 即隨機點的分布面積為100

s = 100

a = 0

cnt

# 有興趣的朋友可以適當的多迴圈幾次,這樣結果更為精確,但是執行時間太長,並沒有跑太久

while (cnt <= 100000)

# 否則為綠色

else

points(x = xx, y = yy, pch = 20, cex = 1, col = "green")

cnt = cnt + 1

}# 計算積分

print((a / 10000) * s)

pdf 徐鍾濟 蒙特卡羅方法 蒙特卡羅計算積分

通常情況下,我們不能解析地求解積分,必須借助其他方法,其中就包括蒙特卡羅積分。你可能還記得,函式的積分可以解釋為函式曲線下的面積。蒙特卡羅積分的工作原理是在a和b之間的不同隨機點計算乙個函式,將矩形的面積相加,取和的平均值。隨著點數的增加,所得結果接近於積分的實際解。蒙特卡羅積分用代數表示 與其他數...

python求定積分的函式 蒙特卡羅計算積分

作者 cory maklin 編譯 vk towards datas science 通常情況下,我們不能解析地求解積分,必須借助其他方法,其中就包括蒙特卡羅積分。你可能還記得,函式的積分可以解釋為函式曲線下的面積。蒙特卡羅積分的工作原理是在a和b之間的不同隨機點計算乙個函式,將矩形的面積相加,取和...

R語言由Monte Carlo方法計算積分

蒙特 卡羅方法 monte carlo method 也稱統計模擬方法,是二十世紀四十年代中期由於科學技術的發展和電子計算機的發明,而被提出的一種以概率統計理論為指導的一類非常重要的數值計算方法。是指使用隨機數 或更常見的偽隨機數 來解決很多計算問題的方法。f function o sum rexp...