概率论上机实验报告.docVIP

实验一 随机变量的概率分布 一 实验目的 掌握计算随机变量分布律或概率密度值的Matlab命令； 掌握计算分布函数的Matlab命令； 学习常见分布的随机变量的模拟与应用。 二 实验内容 作业一：某公司准备招聘320人，其中正式工280人，临时工40人，报考人数1821人，满分400分，知考试平均成绩u=166分，而360分以上有31人，某甲考256分，而问他是否会被录取？可否聘为正式工？（考试成绩～N） 分析: 利用（x-u）/s～N(0,1)通过计算P（x360）=31/1821计算出方差,那个人的分数设为b，根据以上结果只需要比较得280分以上的人数与正式工录用人数的大小，可以知道他是否被录取，能否成为正式工。 解： 编写程序： clear,clc; N=1821, M=166; a=360,pa=31; b=256; x=0; d=a-x; t=norminv(1-pa/N,0,1) i=(a-M)/t x=linspace(66,266,200); y1=normpdf(x,M,i); y2=normcdf(x,M,i); subplot(1,2,1); plot(x,y1); subplot(1,2,2); plot(x,y2); pb=normcdf(b,M,i); paiming=N*(1-pb) if paiming=280 fprintf(he can be offered this job) else fprintf(he can not be offered this job) end 所得结果： N = 1821 a = 360 t = 2.1195 i = 91.5305 paiming = 296.3407 他不会被录取。 作业二：从一批火箭推力装置抽出十个实验，得燃烧时间如下： 50.7 54.9 54.3 44.8 42.2 69.8 53.4 66.1 48.1 34.5， 设燃烧时间～(),求μ和σ的置信水平为0.9的置信区间。 分析: MATLAB中有可以直接求出总体平均值和方差置信区间的函数normfit，可以利用该函数求出相应的置信区间。 解： 编写程序， x=[50.7,54.9,54.3,44.8,42.2,69.8,53.4,66.1,48.1,34.5]; alpha=0.1; [mu,sig,muci,sigci]=normfit(x,alpha); muci,sigci 实验结果： muci = 45.7629 57.9971 sigci = 7.6964 17.3609 三、上机心得体会 第一个题只要将其转化成标准的正态分布，然后将方差解出来。第二个题我将上课所用的理论使用，就是用χ2和t分布，但得到的置信区间非常大，后来发现MATLAB中存在一个专门求置信区间的函数，然后非常简单的求出来了置信区间。感觉MATLAB真是很好用。