用蒙特卡罗方法计算国土面积问题.docVIP

本科学生综合性、设计性 实验报告 姓名＿＿胡延旭＿＿ ＿ 学号＿20091808020036＿＿＿ 姓名＿＿廉俊子＿＿ ＿ 学号＿20091808020037＿＿＿ 姓名＿＿盖枫娇＿＿ ＿ 学号＿20091808020028＿＿＿ 专业＿数学类＿＿＿ ＿ 班级＿09级7班＿＿＿ ＿ 实验课程名称＿＿用蒙特卡罗方法计算国土面积数学模型＿＿＿＿＿＿ 指导教师及职称＿＿崔利宏 教授＿＿＿＿＿＿ ＿ ＿ 开课学期 2010 至＿2011 学年＿ 第二＿学期 上课时间 2010 年 4 月 20 日 辽宁师范大学教务处编印 一、实验方案 实验名称：用蒙特卡罗方法计算国土面积。 实验时间：2011-5-10 小组合作：是● 否○ 小组成员：胡延旭 廉俊子 盖枫娇 实验目的： 熟练掌握蒙特卡罗方法的基本思想。 应用与练习用计算机模拟数学模型。 c)锻炼学生MATLAB编程的能力。 d) 锻炼学生利用数学知识建立模型解决实际问题的能力。 e) 培养学生动手能力与创新意识，增进对数学的兴趣。 实验设备与材料： 硬件设备：微型计算机。 软件设备：WindowsXP操作系统，数学软件MATLAB6.5。 《数学模型与数学建模》教材及相关辅导用书 Matlab21张志勇 Matlab6. X编程 Matlab22辅助小波分析与应用 Matlab23 mat lab 6.5辅助优化计算与设计 Matlab24 matlab6.x图像处理 《mat lab宝典》 实验方法步骤及注意事项 产生随机点于矩形区域内，产生的点数设为N。 如果点产生在曲线所围区域内，则计数器m加1。 根据公式 Aream*S/N （S为矩形面积） 应用程序： Clear all X=[7.0 10.5 13.0 17.5 34.0 40.5 44.5 48.0 56.0 61.0 68.5 76.5 80.5 91… 96.0 101.0 104.0 106.5 111.5 118.0 123.5 136.5 142.0 146.0 150.0 157.0 158.0]; Y1=[44 45 47 50 38 30 30 34 36 34 41 45 46 43 37 33 28 32 65 55 54 52 50 66 66 68]; Y2=[44 59 70 72 93 100 110 110 110 117 118 116 118 118 121 124 121 121… 121 122 116 83 81 82 86 85 68]; L=max(x)-min(x); H=max(y2)-min(y1); for k=1:10 n=10000; m=0; for i=1:n u(i)=unifrnd(min(x),max(x)); v(i)=unifrnd(min(y1),max(y2)); f1=interp(x,y1,u(i),’cubic’); f2=interp(x,y2,u(i),’cubic’); if(v(i)=f1v(i)=f2) m=m+1; end end S=L*H*m/n/18^2*1600 end 注意事项: N取值越大越接近真实值,但不能过大,这样会减慢计算机的运行速度甚至会造成死机.所以n的取值要适当。 注意X，Y1，Y2 的区域 多次试验 注意for循环语句的书写 设计算法时要严谨，以免算法设计不完整 实验数据处理方法： 应用MATLAB软件编程，通过for循环产生随机点，if语句限定曲线所围区域内。 通过多次实验，利用蒙特卡罗方法计算国土面积 参考文献：MATLAB实用教程 郑阿奇 主编 电子工业出版社 数学模型与数学建模 陈汝栋 于延荣 编著 国防工业出版社 数学建模竞赛---获奖论文精选与点评 韩中庚编 科学出版社 MATLAB数学实验 胡良剑 孙小君编著 北京：高等教育出版社 2006.6 6、指导老师对实验设计方案的意见： 指导老师签名： 年 月 日 二、实验报告 1、实验目的、设备与材料、理论依据、实验方法步骤见实验设计方案 2、实验现象、数据及结果 实验现象：取N=10000运行两次得到10个面积，取平均值作为面积。 数据： 44099 41913 42142 42636