matlab用蒙特卡罗法估算定积分.doc

matlab用蒙特卡罗法估算定积分.doc

《matlab用蒙特卡罗法估算定积分.doc》由会员分享，可在线阅读，更多相关《matlab用蒙特卡罗法估算定积分.doc（3页珍藏版）》请在冰豆网上搜索。

用蒙特卡罗法估算定积分

•求解示例：

程序是：

clc

n=10^6;

maxfx=1;

mimfx=0;

a=0;b=1;m=0;

d=1;

c=[];

fori=1:

n

x=a+rand*（b-a）;

y=rand*（maxfx-mimfx）+mimfx;

ify<=x^2

m=m+1;

end

end

fprintf（'积分值是：

%f\n',m/n*（maxfx-mimfx）*（b-a））

结果是：

积分值是：

0.333516

积分值是：

0.333116

积分值是：

0.334410

积分值是：

0.333132

积分值是：

0.333022

积分值是：

0.332899

积分值是：

0.333604

积分值是：

0.333535

积分值是：

0.333117

积分值是：

0.333312

积分值是：

0.334002

积分值是：

0.332966

积分值是：

0.332274

积分值是：

0.334320

积分值是：

0.333143

积分值是：

0.333144

积分值是：

0.333840

积分值是：

0.332857

积分值是：

0.333487

积分值是：

0.333332

积分值是：

0.333501

积分值是：

0.333358

积分值是：

0.333129

积分值是：

0.332874

积分值是：

0.333495

积分值是：

0.334282

积分值是：

0.333183

积分值是：

0.333594

积分值是：

0.333894

积分值是：

0.333766

积分值是：

0.333126

积分值是：

0.332993

积分值是：

0.333032

积分值是：

0.333090

积分值是：

0.332987

积分值是：

0.333853

积分值是：

0.333412

积分值是：

0.332419

积分值是：

0.333338

积分值是：

0.333622

积分值是：

0.334497

积分值是：

0.333565

积分值是：

0.333759

积分值是：

0.332973

积分值是：

0.332783

积分值是：

0.332990

积分值是：

0.333409

积分值是：

0.332452

积分值是：

0.334463

积分值是：

0.332927

积分值是：

0.332749

积分值是：

0.332269

积分值是：

0.333389

积分值是：

0.332691

积分值是：

0.333038

积分值是：

0.333918

积分值是：

0.333517

积分值是：

0.333182

积分值是：

0.333919

积分值是：

0.333463

积分值是：

0.332823

积分值是：

0.333679

积分值是：

0.332826

>>