清华大学数学实验_实验11数据的统计和分析.docVIP

数学实验 实验11 数据的统计和分析 实验11 数据的统计和分析 实验目的： 1.???掌握数据的统计描述和参数估计、假设检验的基本概念与原理，以及用MATLAB实现的方法； 2.???练习用这些方法解决实际问题。 实验内容： 5，炮弹射击的目标为一半径100m的圆形区域，弹着点以圆心为中心呈二维正态分布。设在密度函数式中σx=80m和σy=50m,相关系数r=0.4，求炮弹命中圆形区域的概率。 解： 问题的建模： 设目标中心坐标为x=0，y=0，圆的半径R=100m，则圆形区域表示为 Ω:x2+y2≤R2 二维正态分布的概率密度为 px,y=12πσxσy1-r2exp?{-121-r2[x2σx2-2rxyσxσy+y2σy2]} 由于X,Y方向不相关，炮弹命中圆形区域的概率为 P=Ωpx,ydxdy=Ω12πσxσy1-r2exp?{-121-r2[x2σx2-2rxyσxσy+y2σy2]}dxdy-0.5* 化简得 P=Ω1fx,ydxdy≈1nk=1mfxk,yk fx,y=12πσxσy1-r2exp?{-12[x2σx2-2rxyσxσy+y2σy2]} 其中Ω1是整个圆（引入了相关系数，所以不能用第一象限代替），xk,yk是n个点中落在Ω1内的点的坐标， R=100,σx=80, σy=50，而随机点xi，yi(I=1,2,…,n)分别为[-1，1]和[-1，1]区间上的均匀分布随机数。 模型的求解： Matlab代码如下： r=0.4; sx=0.8;sy=0.5 n=100000;m=0;z=0; x=unifrnd(-1,1,1,n); y=unifrnd(-1,1,1,n); for i=1:n if x(i)^2+y(i)^2=1 u=exp(-0.5*(x(i)^2/sx^2-2*r*x(i)*y(i)/sx/sy+y(i)^2/sy^2)); z=z+u; m+m+1; end end P=4*z/2/pi/sx/sy/n %考虑到正方形边长变成了2，所以还应乘以4 运行四次，结果如下： P = 0.6984 P = 0.6983 P =0.7000 P =0.7002 因此可认为概率在0.7左右。 11，轧钢有两道工序：粗轧和精轧。粗轧钢坯时由于各种随机因素的影响，得到的钢材的长度呈正态分布，其均值可由轧机调整，而方差是设备精度决定的，不能改变；精轧时将粗轧得到的钢材轧成规定的长度(可以认为没有误差)。如果粗轧后的钢材长度大于规定长度，精轧时要把多出的部分轧掉，造成浪费；如果粗轧后的钢材长度已经小于规定长度，则整根报废，浪费更大。问题是已知钢材规定的长度l和粗轧后的钢材长度的均方差δ，求可以调整的粗轧时钢材长度的均值m，使总的浪费最小。试从以下两种目标函数中选择一个，在l=2m，δ=20cm条件下求均值m： (1)每粗轧一根钢材的浪费最小； (2)每得到一根规定长度钢材的浪费最小。 解： 问题的建模： 已知粗轧后钢材的长度符合正态分布U（m，δ）。 其中，m是待求量，δ=0.2m是一定量。 符合规格的长度是l=2m 令每一根总浪费量是W。由数学知识可得： W=l∞x-lpxdx+-∞lxpxdx=-∞∞xpxdx-l∞lpxdx=m-lP 假设粗轧N根钢材，成品材有PN根，则共浪费长度为mN-lPN，粗轧一根平均浪费（mN-lPN)/N=m-lP，成品材一根平均浪费S=(mN-lPN)/PN=m/P-l。 以每得到一根规定长度钢材的浪费最小为目标函数，令 dSdm=P-mp(l)P2=0 则有P-mpl=0 模型的求解： Matlab代码如下： function W=Lab_11(m); s=0.2,n=100000; x=normrnd(m,s,1,n); d1=0,d2=0; for i=1:n; if x(i)=2; d1=d1+x(i)-2; else d2=d2+x(i); end end d=d1+d2; W=d/n; [mm,fv]=fminsearch(@Lab_11,2); 运行结果是： 2. 3312，2.3316， 2.3297，2.3250，因此当m=2.33m时，每粗轧一根浪费最小。