蒙特卡洛法

其實(shí)所謂的蒙特卡洛法就是我們的隨機(jī)化算法啊，最簡(jiǎn)單的應(yīng)用就是撒豆子求面積了。

那就舉兩個(gè)栗子，一個(gè)就是撒豆子，一個(gè)是非線性整數(shù)規(guī)劃問(wèn)題。

撒豆子

?

clc,clear x = unifrnd(0,12,[1,1e7]); y = unifrnd(0,9,[1,1e7]); ? cnt = sum(y<x.^2&x<=3) + sum(y<12-x&x>=3); s = 12*9*cnt/1e7;

?

?

非線性整數(shù)規(guī)劃

可以算一個(gè)概率，假設(shè)目標(biāo)函數(shù)落在高值區(qū)域的概率為? 、? ，則當(dāng)枚舉? 個(gè)點(diǎn)后，至少有一個(gè)點(diǎn)能落在高值區(qū)域的概率為：

?

% mente.m function [f,g] = mente(x) f = x(1)^2+x(2)^2+3*x(3)^2+4*x(4)^2+2*x(5)^2-8*x(1)-2*x(2)-3*x(3)-x(4)-2*x(5); g = [sum(x)-400x(1)+2*x(2)+2*x(3)+x(4)+6*x(5)-8002*x(1)+x(2)+6*x(3)-200x(3)+x(4)+5*x(5)-200];% test.m rand('state',sum(clock)); p0 = 0; tic for i = 1:10^6x = randi([0,99],1,5);[f,g] = mente(x);if all(g<=0)if p0<fx0 = x;p0 = f;endend end toc

?

轉(zhuǎn)載于:https://www.cnblogs.com/TreeDream/p/8297803.html

與50位技術(shù)專家面對(duì)面20年技術(shù)見(jiàn)證，附贈(zèng)技術(shù)全景圖

總結(jié)

以上是生活随笔為你收集整理的蒙特卡洛法的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。