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
>>