- 1、本文档共5页,可阅读全部内容。
- 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
蒙特卡罗方法求冰淇淋
蒙特卡罗方法求冰淇淋锥的体积
200820302021 赵越
问题叙述
蒙特卡罗方法是一种通过随机变量的统计试验求近似解的方法。在冰淇淋锥体和其外围六面体中随机产生N个坐标点,统计其中落在锥体中的点和总的点数的比例就可以近似地求出锥体的体积,在该实验中求出空间中x(-1,1),y(-1.1),z(0.2)区间中的冰淇淋锥体的体积大小。
问题分析
冰淇淋锥含于体积等于8的立方体中,
由于rand 产生0 到1之间的随机数,所以:
x=2*rand-1产生– 1到1之间的随机数;
y=2*rand-1产生– 1到1之间的随机数;
z=2*rand;产生0到2之间的随机数;
N个点均匀分布于六面体中,锥体中占有m个,则锥体与六面体体积之比近似为 m : N
多次反复试验,可以观察出其误差情况。
实验程序及注释
function data=icecream(L) %设置函数名
if nargin==0,L=10;end %设置循环次数最多为10
N=10000; %设置采样点的数量N为10000
for k=1:L %循环L次
P=rand(N,3); %随机取值函数N行3列
x=2*P(:,1)-1; %设置x在(-1,1)范围取值
y=2*P(:,2)-1; %设置y在(-1,1)范围取值
z=2*P(:,3); %设置z在(0,2)范围取值
R2=x.^2+y.^2;R=sqrt(R2); %设置锥体的函数表达
II=find(z=Rz=1+sqrt(1-R2)); %设置锥体的函数表达
m=length(II); %判断落在锥体内的点的数量
q(k)=8*m/N; %计算出锥体的体积大小
end
data=[q; q-pi]; %返回体积大小和计算误差
实验数据结果及分析
程序运行后,分别列出了取值N=10000的时候进行10次同样的蒙特卡罗方法计算出来的锥体体积大小,可以看出小数点后第2位的数值已经很不稳定,因此只能保证两位数字的有效数位,因此修改N的大小,进行相同的实验,当N分别取10^3,10^4,10^5,10^6时可以看出采样点的多少对数据精度的影响。
在N为10^3时,数据精度很差,为10^4时只能保证不到2位的有效数位,随着N的增加计算精度也随着增加,在N^6时,数据比较接近3.14。
实验结论
蒙特卡罗方法近似计算几何体体积,其随着采样点的数量提高,数据精度也逐渐提高,采样点的增加对计算量的影响成正比例关系增加,而不是成几何级数增加,因此使对图形体积的快速计算成为可能,但数据精度较低。
六、注记
蒙特卡罗方法可以用于快速计算函数曲线所围图形面积,即计算定积分,其优点是计算量小,采样点的提高对计算量提高不快,适合快速计算,计算面积的时候和计算体积的时候维度的增加对计算量的增加也是线性的,不像数值方法计算量提高很快,因此蒙特卡罗方法还适合多维度问题的计算。
对抛射曲线的改进实验
200820302021 赵越
一,问题叙述
在研究榴弹炮等类似的抛射曲线发射角度和射程的关系问题中,建立的物理模型为炮弹发射初速不变,炮弹发射后仅在重力的作用下作惯性飞行,飞行的轨迹为抛物线,但在实际问题中,用该模型计算出的射程会与实际的射程发生较大偏差,理论射程会大于实际射程,这主要是由于没有考虑到炮弹在空气中受到空气的助力,速度不断减小的原因。
问题分析
原抛物曲线公式为:
考虑到炮弹在空气中飞行受到空气助力,设空气阻力始终与飞行轨迹的切线方向相反,且阻力大小与速度呈正比例关系,侧有:,其中Z为阻力系数,再根据初值,当t=0时,V=Vo,可得到V=Vo*exp(-Zt),用V代替原抛物线公式中的Vo得到:
先假设Z=0.01
实验程序及注释
function fx=funx(x)
fx=exp(-0.01*x)-sqrt(2)*4.9*x/515; %定义求解抛物线零点的函数
v0=515;alpha
您可能关注的文档
- 考研无机化学 元素化学记忆内容汇总.doc
- 红杉树信息实时协同工作平台产品介绍.ppt
- 耐火材料电熔刚玉行业标准(060926再修改稿).doc
- 聆听往昔的音韵—鼓词调查结题报告.doc
- 聊城市地热回灌试验工作方案-改.doc
- 职业介绍机构资格认定流程图.doc.doc
- 职业健康监督管理.ppt
- 职业道德与职业精神.ppt
- 职称英语串讲点题班讲义一 1至180页.doc
- 职责康乐部经理岗位职责.doc
- 2024至2030年全球及中国乳品机械设备行业市场分析及投资建议报告.docx
- 2024至2030年中国白羽鸡市场调查与行业前景预测专题研究报告.docx
- 2024至2030年中国热镀锌产品制造行业市场预测与投资规划分析报告.docx
- 项目软件项目申请报告.docx
- 2024至2030年中国汽车副车架行业市场调查与前景预测研究报告.docx
- 2024至2030年中国传染病医院行业市场预测与投资规划分析报告.docx
- 2024至2030年全球与中国汽车竞赛轮胎市场现状及未来发展趋势.docx
- 伏特加酒项目申请报告.docx
- 2024至2030年中国上海市大健康产业发展前景预测与投资趋势分析报告.docx
- 2024至2030年中国辛醇行业深度调研及发展预测报告.docx
文档评论(0)