实验11-数据统计和分析.pdfVIP

清华大学《数学实验》实验报告  实验十一 数据的统计与分析 姓名：芦琛璘 班级：化 33 学号：2013011934 实验目的： 1、 掌握数据的统计描述和参数估计、假设检验的基本概念与原理，以及用 MATLAB 实现的 方法。  2、 练习用这些方法解决实际问题  实验内容： 【问题1】炮弹射击的目标为一个半径为 100m的圆形区域，弹着点以圆心为中 心呈二维正态分布，其中 x方向的标准差为 80m，y 方向的标准差为 50m 相关系 数 r=0.4，求炮弹命中圆形区域的概率。 【模型建立】 设目标中心坐标 x=0;y=0,圆的半径为100，则圆形区域表示为： Ω：x y R 二维正态分布的概率密度函数为： 1 1 f x,y exp 2 2 √1 21 代入数据得： 1 1 f x,y exp 2∗0.4 2∗80∗50 10.4 210.4 80 80∗50 50 √ 所以概率为 P ∬ exp 2∗0.4 (积分区域Ω） √ . ∗ ∗∗ . 由于该积分难以求解，所以采用蒙特卡罗方法计算 P 【模型求解】 由于该积分难以求解，所以采用蒙特卡罗方法计算 P，用matlab 编写如下程序： daba.m: sx=80; sy=50; r=0.4; n x=unifrnd(-100,100,1,n); y=unifrnd(-100,100,1,n); m=0;z=0; for i=1:n; if x(i)^2+y(i)^2=100^2 u=exp(-(x(i)^2/sx^2-2*r*x(i)*y(i)/(sx*sy)+y(i)^2/sy^2)/(2*(1- r^2))); z=z+u; m=m+1; 化学工程系 化 33 芦琛璘  清华大学《数学实验》实验报告  end end p=200*200*z/(2*pi*sx*sy*sqrt(1-r^2))/n 【结果如下】 运行多次结果如下： 次数 P 1. 0.6981 2. 0.6981 3. 0.6977 4.