资源描述：
《北科数理统计与Matlab上机报告1》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、统计分析软件（matlab）实验报告1序号班级姓名学号日期时间地点1信计1502吕瑞杰415213352017.06.268:00-11:45实验楼102指导教师：李娜实验名称：一、matlab基本操作、概率计算实验任务：【练习1_01】生成(a,b)上的均匀分布的随机数；生成N(nu,sigma.^2)上的正态分布的随机数。【练习1_02】求个人中至少有两个人生日相同的概率【练习1_03】20个黑白棋子，随机抽出10个，求（1）10个一色（2）9个一色（3）8个一色（4）7个一色（5）6个一色（6）5个一色的概率。【练习1_04】将编号为
2、的本书任意地排列在书架上，求至少有一本书自左到右的排列序号与它的编号相同的概率。【练习1_05】求一个三位数之和为的概率。【练习1_06】（1）生成个二项分布的随机数，统计出每个可能取值的频数、频率；（2）计算其分布律，分别用数值和图形进行对比；（3）算出理论上最可能的取值点并于实际随机数进行比较。【练习1_07】画正态分布的概率密度函数曲线，产生个相应的随机数，画出直方图和带正态密度曲线的直方图。将随机数的频率曲线与概率密度函数曲线画在一起进行比对。【练习1_08】对不同的参数，画出分布的概率密度函数曲线，讨论的不同变化对曲线的影响。【练
3、习1_09】设，求，【练习1_10】求标准正态分布的临界值并画图。【练习1_11】（填充，二维均匀随机数）产生二维均匀分布和正态分布随机数，填充画图并将二维随机点画在一起。【练习1_12】求分布的双侧临界值并画图（变换不同的）。【练习1_13】选择两种离散型和两种连续性随机变量，（1）写出它们的分布，求出相应的数学期望、标准差、方差；（2）产生相应的随机数，求出样本的均值、标准差、方差。20实验目的：1.熟悉MATLAB在概率统计中的若干命令和使用格式。2.学会用Matlab填充画图的方法。3.熟悉Matlab的各种分布基本指令，并学会用M
4、atlab编程实现临界值的图形表示。求各种随机变量的各种数字特征。运行结果：【练习1_01】unifrnd(0,5,2,5)normrnd(0,1,2,5)ans=4.07360.63493.16181.39254.78754.52904.56690.48772.73444.8244ans=-1.34990.72540.7147-0.12411.40903.0349-0.0631-0.20501.48971.4172【练习1_02】a=[20:5:50];forn=1:7;p0(n)=prod(365:-1:365-a(n)+1)/365.
5、^a(n);endp1(n)=1-p0(n)p1=0.41140.56870.70630.81440.89120.94100.9704【练习1_03】2*nchoosek(10,10)/nchoosek(20,10)2*(nchoosek(10,9)*nchoosek(10,1))/nchoosek(20,10)202*(nchoosek(10,8)*nchoosek(10,2))/nchoosek(20,10)2*(nchoosek(10,7)*nchoosek(10,3))/nchoosek(20,10)2*(nchoosek(10,6
6、)*nchoosek(10,4))/nchoosek(20,10)(nchoosek(10,5)*nchoosek(10,5))/nchoosek(20,10)ans=1.0825e-050.00110.02190.15590.47740.3437【练习1_04】p=0;fori=1:100;p=p+(-1)^(i+1)*1/prod(1:i);endpp=0.6321【练习1_05】fork=1:27;n=0;forx=100:999;x1=fix(x/100);x2=rem(fix(x/10),10);x3=rem(x,10);20if
7、k==x1+x2+x3n=n+1;elsen=n;endendp=n/900endp=0.0011p=0.0033p=0.0067p=0.0111p=0.0167p=0.0233p=0.031120p=0.0400p=0.0500p=0.0600p=0.0678p=0.0733p=0.0767p=0.0778p=0.0767p=0.073320p=0.0678p=0.0600p=0.0500p=0.0400p=0.0311p=0.0233p=0.0167p=0.0111p=0.0067p=200.0033p=0.0011【练习1_06】p=
8、binornd(10,0.3,1,100);x=0:10;tabulate(p)a=ksdensity(p,x);plot(x,a,'ro')holdonplot(x,a)hol
