北京期刊网站建设,企业互联网整合营销,网站设计制做报价,app源码下载%%monte_carlo_ff.m%被积函数(二重)function ffmonte_carlo_ff(x,y)ffx*y^2;%函数定义处end%%monte_carlo.m%蒙特卡洛计算二重积分function resultmonte_carlo(a,b,c,d,n,m)%先y后x积分#xff0c;a是x积分下限#xff0c;b是x积分上限#xff0c;c是y积分下限#xff0c;d…%%monte_carlo_ff.m%被积函数(二重)function ffmonte_carlo_ff(x,y)ffx*y^2;%函数定义处end%%monte_carlo.m%蒙特卡洛计算二重积分function resultmonte_carlo(a,b,c,d,n,m)%先y后x积分a是x积分下限b是x积分上限c是y积分下限d是y积分上限nm是蒙特卡洛参数sumxff0;for i1:nsumyff0;xffa(b-a)*rand();for j1:myffc(d-c)*rand();sumyffsumyffmonte_carlo_ff(xff,yff);endaversumyffsumyff/m;sumxffsumxff(b-a)*aversumyff;endresultsumxff/n;end%%show.mclear all;clc;format long;%二重积分的int函数法用于校对蒙特卡洛法syms x y;zx^2*y^2;x11;x22;y13;y24;disp(‘int法结果‘)ans_intint(int(z,y,y1,y2),x,x1,x2)%蒙特卡洛法计算结果disp(‘monte_carlo法结果‘)ans_monte_carlomonte_carlo(1,2,3,4,10000,20000)%二者误差disp(‘误差‘)errorabs(ans_monte_carlo-ans_int)原文http://www.cnblogs.com/blog-3123958139/p/6128149.html
