随机模拟—蒙特卡洛方法 Matlab代码实现

蒙特卡洛方法

蒙特卡洛方法（Monte Carlo method），也称统计模拟方法，是二十世纪四十年代中期由于科学技术的发展和电子计算机的发明，而被提出的一种以概率统计理论为指导的一类非常重要的数值计算方法。是指使用随机数（或更常见的伪随机数）来解决很多计算问题的方法。

基本思想

当所求解问题是某种随机事件出现的概率，或者是某个随机变量的期望值时，通过某种“实验”的方法，以这种事件出现的频率估计这一随机事件的概率，或者得到这个随机变量的某些数字特征，并将其作为问题的解。

应用及matlab代码

我们以面积计算的问题讲解蒙特卡洛算法的实现步骤如下：

1.画出图像

求要计算图像的基本形状。此处以y1=3x；y2=8-x为例，画出图像如下：

2.确定边界

通常为矩形，要求矩形必须完整包括所求图像。

如上图，因为图像交点为（2,6），因此，确定边界如红框所示：（0,0），（8,0），（8,6），（0,6）

3.随机统计

在（矩形）边界范围内随机产生点，并统计落在所求图像中的点。所用到的函数为unifrnd函数。本例子以10^7（7次方在保证运行速度的情况下，基本可以满足准确度）为例。

4.确定面积

用比值法求面积，即落在图像上的点：整个范围内的点=所求图像面积：矩形边界面积。

即所有图像I面积=矩形边界面积*落在图像上的点/整个范围内的点。

相关的知识

本例中使用到的unifrnd函数的相关知识。

R = unifrnd(A,B,M,N,...) or R = unifrnd(A,B,[M,N,...])

即产生m*n阶［a，b］均匀分布U（a，b）的随机数矩阵：unifrnd (a,b,m, n)

如y=(0,8,[1,10000000])，即产生1*10000000阶位于（0,8）之间的数。

5.代码实现

%画函数图像x=0:0.25:10;y1=3*x;y2=8-x;plot(x,y1,x,y2) axis([0 10 0 10]);legend('y1=3x','y2=8-x');title('蒙特卡洛算法');text(2,6,'交点');grid onnumPoint = 1e8;%产生随机点x=unifrnd(0,8,[1,numPoint]);y=unifrnd(0,6,[1,numPoint]);%统计所在所求图形中的点frequency=sum(y<3*x & x<=2)+sum(y<8-x & x>2);%计算面积area=6*8*frequency/numPoint

6.结果测验

由函数图像可知，所求区域面积为：24。

算法运行结果如下：

由结果可知，算法每次运行结果都不一样，这是有概率统计的特性所致。但是结果都稳定在24左右，这与理论值是一样的。

以上实例以规则图形为例，对于不规则图像，本算法同样适用，且更便利。

蒙特卡洛模拟这个名字听起来很厉害，简单来说其实就是进行多次随机抽样，结合题目和代码很容易理解
话不多说，show code!

clear,clc;
count=100000;
a=[20 15 10 5 5 3 3 3 3 1 1 1];%十个数，可以是其他不同的十个
c=1:12;
r=zeros(size(a));%每个位置出现的次数，初始化全部为零
for i=1:count
    b=c(randperm(length(a)));%b就是a中10个数的随机顺序排列
    for j=1:12
        if sum(a(b(1:j)))>=50
            break;
        end
    end
    r(b(j))=r(b(j))+1;%对应位置记录加1个
end
r=r/count;

下面是运行结果

V2G汽车的蒙特卡洛模拟matlab代码，可以从电网充电也可以向电网放电。

充电开始时刻、放电时刻、日行驶里程等参数可以自己设置，易于修改。

ID:6920644246447334

浪迹天涯

技术交流 资源共享 莫如博客陪您进步！

 

本文转载自：莫如博客