蒙特卡罗方法求冰淇淋.doc

蒙特卡罗方法求冰淇淋锥的体积 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>=R&z<=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;al