2011数学建模讲座.ppt

示例2 蒙特卡罗法 思考: 如何取c和d? 如何按均匀分布随机投入n个点? 如何统计落入Ω内的点的个数? 求由曲线y=1/x与直线y=x及x=2所围成图形的面积 %% 蒙特卡罗法演示 % 计算平面闭区域面积Ω %% 初始化 a=1; b=2; %设置定义域 %闭区域Ω的边界 f1=@(x) 1./x; f2=@(x) x; c=0; d=2; %设置覆盖矩形 x=linspace(a,b,100); plot(x,f1(x),x,f2(x)); %绘出闭区域Ω hold on %使得后面所绘图形在前面所绘图形的坐标系下 plot([a,b,b,a,a],[c,c,d,d,c],r); %绘出覆盖矩形S %% 随机投点 n=100; x=a+(b-a)*rand(1,n); y=c+(d-c)*rand(1,n); plot(x, y, .k); %绘出随机点 hold off %% 统计落入闭区域Ω内的点的个数 m=0; for i=1:n if y(i)=f1(x(i)) y(i)=f2(x(i)) m=m+1; end end S=(b-a)*(d-c)*m/n %% 计算近似面积 99创维杯全国大学生数学建模竞赛题目B题 钻井布局 勘探部门在某地区找矿。初步勘探时期已零散地在若干位置上钻井，取得了地质资料。进入系统勘探时期后，要在一