资源描述：
《数学建模实验报告数据的统计分析》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、数据的统计分析一、实验目的及意义本实验旨在通过对一些常见分布的概率计算和概率密度函数、分布函数曲线的直观认识、对数据分布的形态猜测、对某些概率分布的密度函数的参数估计（以正态为例）以及进行简单的正态假设检验，来揭示生活中的随机数据的一些统计规律．二、实验内容1.常见的分布的概率计算、密度函数、分布函数及其图形；2.参数估计；3.正态假设检验。三、实验步骤1.开启软件平台——MATLAB，开启MATLAB编辑窗口；2.根据求解步骤编写M文件3.保存文件并运行；4.观察运行结果(数值或图形)；5.根据观察到的结果和体会写出实验
2、报告。四、实验要求与任务根据实验内容和步骤，完成以下实验，要求写出实验报告。1．某人向空中抛硬币100次，落下为正面的概率为0.5，设这100次中出现正面向上的次数为X，试分别计算X=45和X≤45的概率，并画出分布函数的图形．(用到的matlab函数：binopdf,binocdf)2．设，用MATLAB编程计算：（1）当时，求；（2）若；（3）分别绘制时的概率密度函数图形．(用到的matlab函数：norminv,normpdf,normcdf)3．随机产生1000个服从参数为的指数分布的样本数据，画出直方图，并求参数
3、的估计值和置信水平为99%的置信区间．(用到的matlab函数：hist，exprnd,expfit)4．已知数据如下表X234578111415161819Y106.42108.2109.58110109.93110.49110.59110.6110.9110.76111111.2试建立Y与X之间的函数关系，并检验残差r是否服从均值为0的正态分布．(用到的matlab函数：polyfit,polyval,normplot或ttest或lillietest)五.程序代码及运行结果（经调试后正确的源程序）1．某人向空中抛硬币
4、100次，落下为正面的概率为0.5，设这100次中出现正面向上的次数为X，试分别计算X=45和X≤45的概率，并画出分布函数的图形．(用到的matlab函数：binopdf,binocdf)程序代码：(prog1.m)x=0:100;y=binopdf(x,100,0.5);p1=binopdf(45,100,0.5);p2=binocdf(45,100,0.5);disp(['P(X=45)=',num2str(p1)])disp(['P(X≤45)=',num2str(p2)])plot(x,y,'b-','LineW
5、idth',2);title('X～b(100,0.5)');holdonplot(45,p1,'go','MarkerEdgeColor','k','LineWidth',2,'MarkerFaceColor','g','MarkerSize',8)str1='P(X=45)=';str2=num2str(p1);str=strcat(str1,str2);text(10,0.05,str);str1='P(X≤45)=';str2=num2str(p2);str=strcat(str1,str2);text(10,0.
6、04,str);运行结果：P(X=45)=0.048474P(X≤45)=0.18412．设，用MATLAB编程计算：（1）当时，求；（2）若；（3）分别绘制时的概率密度函数图形．(用到的matlab函数：norminv,normpdf,normcdf)程序代码：(prog2.m)fprintf('(1)X～N(2,0.25)')p1=normcdf(2.9,2,0.5)-normcdf(1.8,2,0.5);p2=1-normcdf(-3,2,0.5);p3=1-normcdf(3.5,2,0.5)+normcd
7、f(0.5,2,0.5);disp(['P(1.8＜X＜2.9)=',num2str(p1)])disp(['P(X＞-3)=',num2str(p2)])disp(['P(
8、X-2
9、＞1.5)=',num2str(p3)])fprintf('(2)X～N(2,0.25)')x=norminv(normcdf(1.8,2,0.5)+0.25,2,0.5);disp(['P(1.8＜X＜x)=2.5,x=',num2str(x)])fprintf('(3)如图')x=0:0.05:4;y1=normpdf(x,2,0
10、.2);y2=normpdf(x,2,0.5);y3=normpdf(x,2,0.9);holdonplot(x,y1,'b-',x,y2,'r-',x,y3,'g-','LineWidth',2);legend('σ=0.2','σ=0.5','σ=0.9');运行结果：(1)X～N(2,0.25)P
