概率论与数理统计上机作业
概率论与数理统计上机实验 (6)
第四次上机内容安排财政收入决定因素的回归分析背景我国从1978年改革开放以来,国民经济一直保持着较高的增长速度,人民生活水平大幅度提高。
但近年来由于全球竞争加剧,中国对外开放程度的加深,国民经济的发展面临着很大的挑战。
由于我国改革以前是高度集中的计划经济,宏观调控政策运用很少,没有面临很多问题。
但改革以后,由于引入了市场经济,宏观调控对国民经济的发展起着至关重要的作用。
财政政策是宏观调控的重要手段,财政平衡对于国家之宏观经济有着决定性的影响力,所以政府都会有年度财政预算,以尽量实现财政平衡,避免财政赤字。
在财政政策的制定上,我国多年来基本上延续使用的是“量入为出”的方法。
在进行国家收入预算时,我们需要比较准确地测定次年度的财政收入,以此拟定预算支出,保障经济平稳发展。
因此,我们提出国家财政收入决定因素统计分析这一课题,通过统计的方法,研究多项因素对我国财政收入的影响,从而建立财政收入预测模型,以达到对现实经济活动进行指导的目的。
财政收入水平的高低是反映一国经济实力的重要标志。
在一定时期内,财政收入规模大小受许多因素的影响,如农业增加值、工业增加值、建筑业增加值、社会总人口、最终消费、是否受灾等因素。
我们认为,一个国家税收水平高低、国民生产总值规模大小、进出口额、社会从业人数的多少、其它收入的多少,是决定一个国家一定时期内财政收入规模的主要影响因素。
因此,在作业中,只取这五个变量作为解释变量,分析他们对财政收入的影响程度,从而达到预测未来财政收入的目的。
数据来源于中国历年统计年鉴2004年。
分析问题1、分析影响财政收入的有哪些因素,对财政收入影响最大的又是哪些因素。
2、分析各因素对财政收入的影响程度,说明各个影响因素影响重要度不同的原因。
3、分析比较财政收入影响因素模型的经济意义。
4、建立财政收入影响因素的模型。
5、一般而言,农业和建筑业的发展应该会使财政收入增加,请对你建立的模型进行分析、探讨,并修正你的模型。
概率论与数理统计作业
第一章随机事件与概率1. 将一枚均匀的硬币抛两次,事件C B A ,,分别表示“第一次出现正面”,“两次出现同一面”,“至少有一次出现正面”。
试写出样本空间及事件C B A ,,中的样本点。
解:{}反正正、正反、反正、反=Ω{}正正、正反=A ,{}正正=B ,{}正正、正反、反正=C2.设31)(=A P ,21)(=B P ,试就以下三种情况分别求)(A B P :(1)AB =∅,(2)B A ⊂,(3)81)(=AB P解:(1)5.0)()()()()(==-=-=B P AB P B P AB B P A B P(2)6/13/15.0)()()()()()(=-=-=-=-=A P B P AB P B P AB B P A B P (3)375.0125.05.0)()()()(=-=-=-=AB P B P AB B P A B P3.某人忘记了电话号码的最后一个数字,因而随机的拨号,求他拨号不超过三次而接通所需的电话的概率是多少?如果已知最后一个数字是奇数,那么此概率是多少?解: 记H 表拨号不超过三次而能接通。
Ai 表第i 次拨号能接通。
注意:第一次拨号不通,第二拨号就不再拨这个号码。
如果已知最后一个数字是奇数(记为事件B )问题变为在B 已发生的条件下,求H 再发生的概率。
4.进行一系列独立试验,每次试验成功的概率均为,试求以下事件的概率: (1)直到第r 次才成功;(2)在n 次中取得)1(n r r ≤≤次成功;解: (1)p p P r 1)1(--= (2)r n rr np p C P --=)1( 5. 设事件A ,B 的概率都大于零,说明以下四种叙述分别属于那一种:(a )必然对,(b )必然错,(c )可能对也可能错,并说明理由。
(1)若A ,B 互不相容,则它们相互独立。
(2)若A 与B 相互独立,则它们互不相容。
(3)()()0.6P A P B ==,则A 与B 互不相容。
概率论与数理统计上机实验报告
概率论与数理统计上机实验报告实验一【实验目的】熟练掌握 MATLAB 软件的关于概率分布作图的基本操作会进行常用的概率密度函数和分布函数的作图绘画出分布律图形【实验要求】掌握 MATLAB 的画图命令 plot掌握常见分布的概率密度图像和分布函数图像的画法【实验容】2 、设X : U (−1,1)(1 )求概率密度在 0 ,0.2 ,0.4 ,0.6 ,0.8,1 ,1.2 的函数值;(2 )产生 18 个随机数(3 行 6 列)(3 )又已知分布函数F ( x) = 0.45 ,求x(4 )画出X 的分布密度和分布函数图形。
【实验方案】熟练运用基本的MATLAB指令【设计程序和结果】1.计算函数值Fx=unifcdf(0, -1,1)Fx=unifcdf(0.2, -1,1)Fx=unifcdf(0.4, -1,1)Fx=unifcdf(0.6, -1,1)Fx=unifcdf(0.8, -1,1)Fx=unifcdf(1.0, -1,1)Fx=unifcdf(1.2, -1,1)结果Fx =0.5000Fx =0.6000Fx =0.7000Fx =0.8000Fx =0.9000Fx =1Fx =12.产生随机数程序:X=unifrnd(-1,1,3,6)结果:X =0.6294 0.8268 -0.4430 0.9298 0.9143 -0.7162 0.8116 0.2647 0.0938 -0.6848 -0.0292 -0.1565 -0.7460 -0.8049 0.9150 0.9412 0.6006 0.83153.求x程序:x=unifinv(0.45, -1,1)结果:x =-0.10004.画图程序:x=-1:0.1:1;px=unifpdf(x, -1,1);fx=unifcdf(x, -1,1);plot(x,px,'+b');hold on;plot(x,fx,'*r');legend('均匀分布函数','均匀分布密度');结果:【小结】运用基本的MATLAB指令可以方便的解决概率论中的相关问题,使数学问题得到简化。
概率论与数理统计上机题目及答案
实验项目一:数据整理中的统计计算一、实验要求:(1)掌握Excel中基本的数据处理方法;(2)学会使用Excel进行统计分组,能以此方式独立完成相关作业。
二、实验重点:了解数据整理的概念和内容。
掌握不同类型的统计图表。
三、实验难点:不同类型的统计图表四、实验要求:0、本实验课程要求学生已修《计算机应用基础》或类似课程。
此条为整门课程所要求,以后不再赘述。
1、已学习教材相关内容,理解数据整理中的统计计算问题;已阅读本次实验导引,了解Excel中相关的计算工具。
2、准备好一个统计分组问题及相应数据(可用本实验导引所提供问题与数据)。
3、以Excel文件形式提交实验报告(含:实验过程记录、疑难问题发现与解决记录(可选))。
此条为所有实验所要求,恕不赘述。
五、实验内容:1、在一批灯泡中随机抽取50只,测试其使用寿命,原始数据如下(单位:小时):700 716 728 719 685709 691 684 705 718706 715 712 722 691708 690 692 707 701708 729 694 681 695685 706 661 735 665668 710 693 697 674658 698 666 696 698706 692 691 747 699682 698 700 710 722进行等距分组,整理成频数分布表,并绘制频数分布图(直方图、折线图、曲线图)。
要求:(1)用MIN和MAX函数找出最小值和最大值,以50为组距,确定每组范围;(2)进行等距分组,整理成频数分布表,并绘制频数分布图(直方图、折线图、曲线图)。
3、温州市1978-2005年GDP(亿元)如下表要求:(1)作出趋势图(折线图或X-Y散点图);(2)用“添加趋势线”方法,找出一个最好的方程;(3)预测2006年、2007年温州市GDP。
4、书P140,6.4六、实验步骤与结果:1、实验项目二:数字特征的统计计算一、实验要求:学会使用Excel计算各种数字特征,能以此方式独立完成相关作业。
概率论与数理统计MATLAB上机实验报告
《概率论与数理统计》MATLAB上机实验实验报告一、实验目的1、熟悉matlab的操作。
了解用matlab解决概率相关问题的方法。
2、增强动手能力,通过完成实验内容增强自己动手能力。
二、实验内容1、列出常见分布的概率密度及分布函数的命令,并操作。
概率密度函数分布函数(累积分布函数) 正态分布normpdf(x,mu,sigma) cd f(‘Normal’,x, mu,sigma);均匀分布(连续)unifpdf(x,a,b) cdf(‘Uniform’,x,a,b);均匀分布(离散)unidpdf(x,n) cdf(‘Discrete Uniform’,x,n);指数分布exppdf(x,a) cdf(‘Exponential’,x,a);几何分布geopdf(x,p) cdf(‘Geometric’,x,p);二项分布binopdf(x,n,p) cdf(‘Binomial’,x,n,p);泊松分布poisspdf(x,n) cdf(‘Poisson’,x,n);2、掷硬币150次,其中正面出现的概率为0.5,这150次中正面出现的次数记为X(1) 试计算X=45的概率和X≤45 的概率;(2) 绘制分布函数图形和概率分布律图形。
答:(1)P(x=45)=pd =3.0945e-07P(x<=45)=cd =5.2943e-07(2)3、用Matlab软件生成服从二项分布的随机数,并验证泊松定理。
用matlab依次生成(n=300,p=0.5),(n=3000,p=0.05),(n=30000,p=0.005)的二项分布随机数,以及参数λ=150的泊松分布,并作出图线如下。
由此可以见得,随着n的增大,二项分布与泊松分布的概率密度函数几乎重合。
因此当n足够大时,可以认为泊松分布与二项分布一致。
4、 设22221),(y x e y x f +−=π是一个二维随机变量的联合概率密度函数,画出这一函数的联合概率密度图像。
《概率论与数理统计》课后习题及答案解析(上)
891010
1113
即有P(A)P(A1)P(A1A2)P(A1A2A3).
10101010
解法2沿用解法1的记号,知
P(A)1P(拨号3次都接不通)1P(A1A2A3)
1P(A|AA)P(A|A)P(A)17893.
312211
891010
五、某类产品毎百件成批,出厂验收时,规定从毎批中任意挑选5件为样品,若样品中发现有废品,则整批不予出厂.今有一批产品100件,其中有6件废品,问这批产品被拒绝出厂的概率有多大?
N(S)3
四、甲、乙两艘油轮驶向一个不能同时停泊两艘油轮的码头,它们都将在某日8时至
20时抵达码头,甲轮卸完油要1小时,乙轮要2小时,假设毎艘油轮在8时至20时的毎一时刻抵达码头的可能性相同.
1.求甲、乙两轮都不需要等候空出码头的概率;
2.设A表示甲、乙同一时刻抵达码头,问A是否为3C1C2C3C3C3
P(B)=P(Ai)P(B|Ai)393 983 9796
i0
C3C3C3C3C3C3C3C3
1212121212121212
121
2714277
211
(21)2
441
0.1458
220552205555445511553025
211
2.利用贝叶斯公式求得P(A
解:P(B|AUB)
P(AUB)
P(BAUBB)
P(A)P(B)P(AB)
P(AB)
P(A)P(B)P(AB)
P(A)P(AB)
0.60.5
0.10.125
P(A)P(B)P(AB)
0.60.70.50.8
2.某人有一笔资金,他投入基金的概率为0.58,购买股票的概率为0.28,两项投资都
概率统计上机作业二
《概率统计》上机作业(二)一、上机目的一、应用M a t l a b计算分布函数值;二、掌握M a t l a b计算随机变量的数字特征的计算方法.二、上机内容分布函数值的计算和随机变量的数字特征三、上机作业一、设一次试验中事件A发生的概率为p,那么在n次独立重复试验中,事件A恰好发生K次的概率为P_K。
试用MATLAB计算当n=100,p=0.6,k=20的概率值.y=binocdf(20,100,0.6)y =3.4204e-016二、设X~N(3, 22)(1)求(2)确定c,使得p=normcdf(5,3,22)-normcdf(2,3,22)p =0.0543p=normcdf(10,3,22)-normcdf(-4,3,22)p =0.2497p=1-(normcdf(2,3,22)-normcdf(-2,3,22))p =}3{},2{},104{},52{>><<-<<XPXPXPXP}{}{cXPcXP<=>0.9282p=1-normcdf(3,3,22)p =0.5000C=3三、已知某保险公司发现索赔要求中有25%是因被盗而提出的。
某年该公司收到10个索赔要求,试求其中包含不多于4个被盗索赔的概率.p=binocdf(4,10,0.25)p =0.9219四、假设一年中,某类保险者里面每个人死亡概率为0.05,现有1000人参加这类保险,试求在未来一年里,被保险者中有10人死亡的概率,并画泊松分布图.binopdf(10,1000,0.05)ans =2.2735e-012计算E(X)、E(2X+1)、E(X 2)-(E(X))2 .x=[-1 0 1 2];p=[0.4 0.2 0.1 0.3];EX=sum(x.*p)EX =0.3000>> x=[-1 0 1 2];p=[0.4 0.2 0.1 0.3];EX=sum(x.*p);y=x*2+1;EY=sum(y.*p)EY =1.6000x=[-1 0 1 2];p=[0.4 0.2 0.1 0.3];EX=sum(x.*p);y=x.^2;EY=sum(y.*p);t=EY-EX^2t =1.6100六、某公司年损失金额的概率分布列为:试计算该公司的期望值和标准差x=[500,1000,1500,2000];p=[0.82,0.15,0.02,0.01];EX=sum(x.*p)EX =610X=[500,1000,1500,2000];P=std(x,1)P =1.1180七、某保险公司1990年—1996年的保费收入如下表,年度 1990 1991 1992 1993 1994 1995 1996保费收入(万元) 104 162 188 264 320 400 442求:保费收入的平均值;样本方差和样本标准差,方差和标准差x=[104 162 188 264 320 400 442];mean(x)⎪⎪⎭⎫ ⎝⎛01.002.015.082.0200015001000500ans =268.5714x=[104 162 188 264 320 400 442];d=var(x,1)d =1.3564e+004x=[104 162 188 264 320 400 442];d1=var(x)d1 =1.5825e+004x=[104 162 188 264 320 400 442];s=std(x,1)s =116.4656x=[104 162 188 264 320 400 442];s=std(x)s =125.7973四、上机心得体会通过这次的上机,知道应用类问题到解决相对普通计算题的解答要相对难一些。
概率论与数理统计上机实验 (4)
第三次上机内容安排
广告策略、地理位置与销量之间的关系
作为对广告策略长期研究的一部分,某研究中心对BG公司的营销策略和广告策略进行了一项跟踪研究,以调查广告、地理位置与销量之间的关系。
本研究选择了60个规模相当的连锁超市组成了一个样本,其中20个位于广州,20个位于上海、20个位于北京。
对中选的每个超市的月度销量进行了一项标准化的数据统计,使得数据具有可对比性,整理后的数据资料如表所示,较高的得分表示较高的销量水平。
研究的第二部分考虑地理位置与某项广告策略之间是否存在明显的关系,这项广告策略包括电视、促销措施、报纸等。
这项研究也选择了60个规模相当的连锁超市组成了一个样本,其中20个位于广州,20个位于上海、20个位于北京。
这个研究记录的销售水平资料如表2:
实验内容
(1)资料1中的广州超市声称其月销量为5.8,且服从正态分布,据上述数据可否相信其表述?
(2)若上述广州超市的销量方差为5,在a=0.05的显著性水平下,能否认为其销量比较稳定?
(3)用描述统计方法概括说明两部分研究的资料,关于销量的初步观测结果是什么?
(4)对两个数据集进行方差分析,陈述每种情况下被检验的假设和研究结论?(5)说明上述两项资料中的单个均值处理分开的合理性。
(6)讨论这个研究的推广和你认为有用的其他分析。
概率论与数理统计上机题目及答案
867
988
849
958
934
1000
878
978
816
1001
918
1040
854
1098
900
869
949
890
1038
927
878
1050
924
905
954
890
1006
926
900
999
886
907
956
900
963
838
961
948
950
900
937
864
919
863
981
916
878
2、准备好一个或几个抽样推断问题及相应数据(可用本实验导引所提供问题与数据)。五、实验内容:
1、问题与数据
某厂用某机床加工某种零件,假设零件长度服从正态分布。现从一批该种中随机抽取10件,测得其长度如下(单位:cm):
6.1 5.7 6.5 6.0 6.3 5.8 6.3 6.1 5.9 6.4
要求:
(1)作出趋势图(折线图或X-Y散点图);
(2)用“添加趋势线”方法,找出一个最好的方程;
(3)预测2006年、2007年温州市GDP。
4、书P140,6.4
六、实验步骤与结果:
1、
实验项目二:数字特征的统计计算
一、实验要求:
学会使用Excel计算各种数字特征,能以此方式独立完成相关作业。
二、实验重点:掌握各种数字特征的函数求解方法。
891
870
986
913
850
911
886
950
926
(0198)《概率论与数理统计》网上作业题及答案
[0198]《概率论与数理统计》第二次作业[判断题]已知电子管的寿命X(小时)服从参数λ=1/1000的指数分布,则这电子管使用在1000小时以上的概率为1/e.参考答案:正确[单选题]当随机变量X服从( )分布时,DX=EX。
A:指数B:均匀C:泊松D:两点参考答案:C[单选题]设X的分布函数为F(x),则F(x)不满足的是( ).A:单调不减.B:x趋于正无穷时,F(x)以1为极限;C:x趋于负无穷时,F(x)以0为极限D:在每一点F(x)为左连续的。
参考答案:D[单选题]设随机变量为X与Y,已知DX=25,DY=36,相关系数ρ=0.4,则D(X-Y)=( ).A:85B:61C:11D:37参考答案:D[单选题]设随机变量X服从参数为2的泊松分布,用切比雪夫不等式估计P(|X-2|≥4)( ).A:≤1/8B:≥1/8C:≤7/8D:≥7/8参考答案:A[单选题]X为连续型随机变量,a为一常数,则P(X= a)的值( )A:必为零B:不一定为零.C:可能不为零。
参考答案:A[单选题]每次射击中靶的概率为0.7,现独立射击10次,用随机变量X表示命中的炮弹数,则X服从( )分布.A:二项.B:几何C:均匀.D:超几何参考答案:A[判断题]随机变量X的取值为不可列无穷多,则X必为连续型随机变量。
参考答案:错误[判断题]X服从二项分布B(n,p),Y服从二项分布B(m,p),且X与Y独立,则X+Y服从二项分布B(n+m,p)。
参考答案:正确[判断题]随机变量X、Y独立,则X与Y必不相关。
参考答案:正确[判断题]设X与Y独立,且X有概率密度函数为f(x), Y有概率密度函数为p(y),则(X,Y)的联合分布密度 f ( x , y )=f(x)p(y) .参考答案:正确第一次作业[判断题]袋内装有5个白球,3个黑球,从中任取两个球,则取到的两个球都是白球的概率为5/14。
参考答案:正确[判断题]设A、B、C表示三个事件, 用A、B、C的运算关系表示"A、B、C恰有一个不发生”为A+B+C. 参考答案:错误[判断题]从1,2,3,4,5,6这六个数中随机的、有放回的连续抽取4个,则"取到的4个数字完全不同”的概率为5/18.参考答案:正确[判断题]概率为零的事件一定是不可能事件.参考答案:错误[判断题]事件A、B相互独立,则有P(AB)=P(A)P(B) .参考答案:正确[判断题]设A、B为两个事件,P(A+B)=1,则A+B必为必然事件。
概率论与数理统计大作业二
班级:姓名:学号:分数:
《概率论与数理统计》大作业二(100分)
1、设二维随机变量(,)X Y 的联合概率密度为20112,(,)0,
y x y f x y ≤≤≤⎧=⎨⎩其他,求(1)(),(),()E X E Y E XY ;
(2)(),()D X D Y ;(3)判断,X Y 是否相关.(20分)
2、设随机变量(,)X Y 联合分布律为
讨论(1)Y X ,是否独立?
(2)Y X ,是否相关?(15分)3、设1210,,,X X X 是来自正态总体2
(0,0.3)X N 的一个简答随机样本。
求(1)(0)P X ≥;(2)求10
21( 1.44)i i P X =>∑.(已知20.1(10)16χ=)(10分)
班级:姓名:学号:
4、计算机在进行加法时,每个加数取整数(四舍五入),设所有取整误差是相互独立的,且它们都在[-0.5,0.5]上服从均匀分布。
(1)若将1500个数相加,问误差总和的绝对值超过15的概率是多少?(2)最少几个数相加在一起可使得误差总和的绝对值小于10的概率不超过90%?(15分)
5、设总体X 的概率密度函数1,0(),0,0x e x f x x θθ-⎧>⎪=⎨⎪≤⎩
12,,,n X X X 是取自总体X 的简单随机样本。
(1)求θ的矩估计量ˆθ
;(2)求θ的极大似然估计量ˆθ.(20分)6、设总体X 的概率分布为:其中θ是未知参数.
总X 有如下的样本值:3,1,3,1,3,1,2,3,求θ的矩估计值和极大似然估计值.(20
分)。
概率论上机实验报告
概率论上机实验报告概率论上机实验报告引言:概率论是数学中的一个重要分支,它研究的是随机现象的规律性。
概率论的应用十分广泛,涵盖了自然科学、社会科学、工程技术等各个领域。
为了更好地理解概率论的基本概念和方法,我们进行了一系列的上机实验,通过实际操作来探索概率事件的发生规律以及概率计算的方法。
实验一:硬币抛掷实验在这个实验中,我们使用了一枚标准的硬币,通过抛掷硬币的方式来研究硬币正反面出现的概率。
我们抛掷了100次硬币,并记录了每次抛掷的结果。
通过统计实验结果,我们可以得出硬币正反面出现的频率。
实验结果显示,硬币正面出现的次数为55次,反面出现的次数为45次。
根据频率的定义,我们可以计算出正面出现的概率为55%。
这个结果与我们的预期相符,说明硬币的正反面出现具有一定的随机性。
实验二:骰子掷掷实验在这个实验中,我们使用了一个六面骰子,通过投掷骰子的方式来研究各个面出现的概率。
我们投掷了100次骰子,并记录了每次投掷的结果。
通过统计实验结果,我们可以得出各个面出现的频率。
实验结果显示,骰子的六个面出现的次数分别为15次、18次、17次、16次、19次和15次。
根据频率的定义,我们可以计算出各个面出现的概率分别为15%、18%、17%、16%、19%和15%。
这个结果表明,在足够多次的投掷中,各个面出现的概率是相等的。
实验三:扑克牌抽取实验在这个实验中,我们使用了一副标准的扑克牌,通过抽取扑克牌的方式来研究各个牌面出现的概率。
我们随机抽取了100张扑克牌,并记录了每次抽取的结果。
通过统计实验结果,我们可以得出各个牌面出现的频率。
实验结果显示,各个牌面出现的次数相差不大,都在10次左右。
根据频率的定义,我们可以计算出各个牌面出现的概率都约为10%。
这个结果说明,在足够多次的抽取中,各个牌面出现的概率是相等的。
实验四:随机数生成实验在这个实验中,我们使用了计算机生成的随机数,通过生成随机数的方式来研究随机数的分布规律。
北京科技大学概率论与数理统计上机报告3
概率论与数理统计第三次上机报告专业:信息与计算科学班级:信计1502(35组)学生姓名:吕瑞杰陈炎睿何芝芝指导教师:***完成时间:2022年4月27日Matlab 概率论与数理统计上机练习(3)五、假设检验【例】(离散型分布检验)某工厂近五年发生了63起事故,按星期几可以分为[9 10 11 8 13 12],问该厂发生的事故数是有与星期几有关?-------------------------------------------------------------------------------------63605 1.666711.07050.8931--------------------------------------------------------------------------------------α=),对数学分析I【练习3.1】(基本计算,两个正态总体的假设检验,检验水平0.05(1)求课程中“专业(数学、信计)””的考试人数、平均分、最小值、最大值、极差、标准差、及格人数、及格率、优良人数(大于等于80)、优良率;写出标准差的计算公式。
(2)对“专业(数学、信计)”,检验方差、平均分是否相等。
(3)对“专业(数学、信计)”,检验及格率、优秀率是否相等。
(4)对“全体成绩”的分布进行检验,首先估计期望和方差,画出正态分布的密度函数曲线以及样本密度散点,对假设的正态分布进行检验。
Matlab程序实现:sy=[6060636340696560726762788290696072767893696895 7183607373607477718570896061776268607066847469 6160867369747174];se=[5081676577717662896565626260788166708053696661 4866696160608552686074606243616060647074657379 6043766663606068606060677464];alpha=0.05; %取显著水平为0.05sy1=length(sy);se1=length(se);%人数sy2max=max(sy);sy2min=min(sy);se2max=max(se);se2min=min(se);%最大值,最小值sy3=range(sy);se3=range(se);%极差sy4=mean(sy);se4=mean(se);%平均分sy5=sqrt(sum((sy-sy4).^2)/(sy1));se5=sqrt(sum((se-se4).^2)/(se1));%标准差%sy5=std(sy);se5=std(se);或sy6=length((find(sy>=60)));se6=length((find(se>=60)));%及格人数sy7=length((find(sy>=80)));se7=length((find(se>=80)));%优秀人数sy8=sy6/sy1;se8=se6/se1;%及格率sy9=sy7/sy1;se9=se7/se1;%优秀率fprintf('\t人数\t 平均分\t 最小值\t 最大值\t极差\t\t标准差\t\t及格人数及格率\t优秀人数优秀率\n');fprintf('--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------\n');fprintf('数学%4d\t%10.4f\t%4d\t\t%4d\t\t%4d\t\t%10.4f\t%4d\t %10.4f\t%4d\t%10.4f\n',sy1,sy4,sy2min,sy2max,sy3,sy5,sy6 ,sy8,sy7,sy9)fprintf('信计%4d\t%10.4f\t%4d\t\t%4d\t\t%4d\t\t%10.4f\t%4d\t %10.4f\t%4d\t%10.4f\n',se1,se4,se2min,se2max,se3,se5,sy6, sy8,se7,se9)fprintf('\n');%方法一fprintf('检验数学和信计的方差是否相等\n');[h1,p1,varci1,stats1]=vartest2(sy,se,alpha,'both');if(h1==0)disp('结果:方差相等');elsedisp('结果:方差不相等');endfprintf('\n');% %方法二% F=sy5^2/se5^2;%统计量F,满足F分布% alpha=0.05; %取显著水平为0.05% Fla1=finv(alpha/2,sy1-1,se1-1);Fla2=finv(1-alpha/2,sy1-1,se1-1);%求F的临界值% if (F>Fla1 && F<Fla2)% MM='数学分析1和数学分析2的方差无显著差异';% else% MM='数学分析1和数学分析2的方差有显著差异';% end% fprintf('检验数学分析1和数学分析2的方差是否相等\n');% fprintf('统计量F的值\t\t\t显著性水平\t\t临界值\t\t\t\t\t检验结果\n');% fprintf(' %.4f\t\t\t\t%.4f\t\t\t%.4f\t\t\t%15s\n',F,alpha,Fla1,MM);% fprintf('\n\n');%方法一fprintf('检验数学和信计的平均分是否相等\n');[h2,p2,muci2,stats2]=ttest2(sy,se,alpha,'both');if(h2==0)disp('结果:平均分相等');elsedisp('结果:平均分不相等');endfprintf('\n');%%方法二% %%%%方法三% sw=((sy1-1)*sy5^2+(se1-1)*se5^2)/(sy1+se1-2);% T=(sy4-se4)/sw/sqrt(1/sy1+1/se1);%统计量T,满T分布% Tla1=tinv(alpha/2,sy1+se1-2);Tla2=tinv(1-alpha/2,sy1+se1-2);%求出T的临界值% if (abs(T)<tinv(1-alpha/2,sy1+se1-2))% XX='数学分析1和数学分析2的平均分无显著差异';% else% XX='数学分析1和数学分析2的平均分有显著差异';% end% fprintf('检验数学分析1和数学分析2的平均分是否相等\n');% fprintf('统计量T的值\t\t\t显著性水平\t\t临界值\t\t\t\t\t检验结果\n ');% fprintf(' %.4f\t\t\t\t%.4f\t\t\t%.4f??%.4f\t\t%15s\n',T,alpha,Tla1,Tla2,XX);% fprintf('\n\n');p=(sy7+se7)/(sy1+se1);U=(sy9-se9)/sqrt((sy9+se9)*p*(1-p));Ua=norminv(1-alpha/2);if(abs(U)>Ua)disp('优秀率无显著差异');elsedisp('优秀率有显著差异');endp=(sy6+se6)/(sy1+se1);U=(sy8-se8)/sqrt((sy8+se8)*p*(1-p));Ua=norminv(1-alpha/2);if(abs(U)>Ua)disp('及格率无显著差异');elsedisp('及格率有显著差异');end[h3,p3,kstat3,critval3]=lillietest(sy,alpha);if(h3==1)disp('数学不是正态分布')elsedisp('数学是正态分布')end[h4,p4,kstat4,critval4]=lillietest(se,alpha);if(h4==1)disp('信计不是正态分布')elsedisp('信计是正态分布')end%%hist(sy)%直方图%[h5,p5,stats5]=chi2gof(sy)%可以检验分布%cdf=[sy,normcdf(sy,sy4,sy5)]%[h5,p5,ksstat,cv5]=kstest(sy,cdf)% a=0:1:100;% a=a';% CDF=[a,cdf(a,sy4,sy5)];% h = kstest(sy,CDF,0.05);S=[6565688174766869827774667372776062816668766074 8090696063686769626060606760776760606071726066 6186646060897374434061956970626663627874605076 6265847069837343717073716974606160707478489364 61797153606052636061656278606560858589696660]; [h,p,jbstat,critval]=jbtest(S,alpha);if(h==0)disp('服从正态分布');elsedisp('不服从正态分布');endsavg=mean(S);svar=var(S);x=20:130;y=normpdf(x,savg,sqrt(svar));d=5;a=20:d:130;pdf=hist(S,a)./length(S)./d;plot(x,y,'r');hold onscatter(a,pdf,'filled');hold off输出:>> lx3_1_lrj_41521335人数平均分最小值最大值极差标准差及格人数及格率优秀人数优秀率--------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------------数学5470.666740955510.0885530.9815 90.1667信计6065.51674389469.2313530.9815 50.0833检验数学和信计的方差是否相等结果:方差相等检验数学和信计的平均分是否相等结果:平均分不相等优秀率有显著差异及格率有显著差异数学不是正态分布Warning: P is less than the smallest tabulated value, returning 0.001.> In lillietest (line 206)In lx3_1_lrj_41521335 (line 99)信计不是正态分布服从正态分布六、方差分析【例1】(单因素方差分析)考虑温度对某化工产品得率的影响,选择五种不同温度进行试验,每一温度各做三次试验。
概率论与数理统计上机题目及答案
实验项目一:数据整理中的统计计算一、实验要求:(1)掌握Excel中基本的数据处理方法;(2)学会使用Excel进行统计分组,能以此方式独立完成相关作业。
二、实验重点:了解数据整理的概念和内容。
掌握不同类型的统计图表。
三、实验难点:不同类型的统计图表四、实验要求:0、本实验课程要求学生已修《计算机应用基础》或类似课程。
此条为整门课程所要求,以后不再赘述。
1、已学习教材相关内容,理解数据整理中的统计计算问题;已阅读本次实验导引,了解Excel中相关的计算工具。
2、准备好一个统计分组问题及相应数据(可用本实验导引所提供问题与数据)。
3、以Excel文件形式提交实验报告(含:实验过程记录、疑难问题发现与解决记录(可选))。
此条为所有实验所要求,恕不赘述。
五、实验内容:1、在一批灯泡中随机抽取50只,测试其使用寿命,原始数据如下(单位:小时):700 716 728 719 685709 691 684 705 718706 715 712 722 691708 690 692 707 701708 729 694 681 695685 706 661 735 665668 710 693 697 674658 698 666 696 698706 692 691 747 699682 698 700 710 722进行等距分组,整理成频数分布表,并绘制频数分布图(直方图、折线图、曲线图)。
要求:(1)用MIN和MAX函数找出最小值和最大值,以50为组距,确定每组范围;(2)进行等距分组,整理成频数分布表,并绘制频数分布图(直方图、折线图、曲线图)。
3、温州市1978-2005年GDP(亿元)如下表要求:(1)作出趋势图(折线图或X-Y散点图);(2)用“添加趋势线”方法,找出一个最好的方程;(3)预测2006年、2007年温州市GDP。
4、书P140,6.4六、实验步骤与结果:1、实验项目二:数字特征的统计计算一、实验要求:学会使用Excel计算各种数字特征,能以此方式独立完成相关作业。
《概率论与数理统计》作业
1、已知为三个事件,则中至少有两个发生这一事件可以表示为 ,,A B C ,,A B C .2、设事件互不相容,且,则,A B (),()P A p P B q ==()P AB = .3、设A ,B 为随机事件,且P (A )=0.7,P (A −B )=0.3,则P (AB )= .4、设A ,B ,C 为三事件,且P (A )=P (B )=1/4,P (C )=1/3且P (AB )=P (BC )=0, P (AC )=1/12,则A ,B ,C 至少有一个发生的概率为 .5、盒子里有4个红球鞋,5个白球,现从中任取两个,恰好有2个红球的概率为 . 二、计算题 1、从52张扑克牌中任意取出13张,问有5张黑桃,3张红心,3张方块,2张梅花的概率是多少?2、从一批由45件正品,5件次品组成的产品中任取三件,求其中恰的一件次品的概率.3、在电话号码簿中任取一电话号码,求后面4个数全不相同的概率(设后面4个数中的每一个数都等可能地取自0,1,2,……,9).4、一个袋内装有大小相同的7个球,其中4个是白球,3个是黑球,从中一次抽取3个,计算至少有两个是白球的概率.5、从5双不同的鞋子中任取4只,求这4只鞋子中至少有两只鞋子配成一双的概率.6、两人约定上午9∶00~10∶00在公园会面,求一人要等另一人半小时以上的概率.1、设事件,A B 满足:11(|)(|,()33P B A P B A P A ===,则 ()P B =.2、设()()()2212()P A B =U ,,53P A P A B P B A ===,则 .3、设随机事件满足,则 ,A B ()0.6,()0.9,(|)0.5P A P A B P B A ==U =()P B =.4、设事件两两独立,且,,A B C 1,()()()2ABC P A P B P C =∅==<,,则 ()9/16P A B C =U U ()P A =__________.5、设在一次试验中,事件A 发生的概率为p . 现进行次独立试验,则n A 至少发生一次的概率为 ,而事件A 至多发生一次的概率为 . 二、计算题1、某光学仪器厂制造的透镜,第一次落下时打破的概率是1/2,第一次落下未打破,第二次落下打破的概率为7/10,若前两次落下还未打破,第三次落下被打破的概率为9/10,试求透镜落下三次而未打破的概率.2、将两信息分别编码为A 和B 传递出来,接收站收到时,A 被误收作B 的概率为0.02,而B 被误收作A 的概率为0.01.信息A 与B 传递的频繁程度为2:1.问(1) 接收站收到信息A 的概率是多少?(2) 若接收站收到的信息是A ,试问原发信息是A 的概率是多少?3、加工某一零件需要经过四道工序,设第一、二、三、四道工序的次品率分别为0.02,0.03,0.05,0.03,假定各道工序是相互独立的,求加工出来的零件的次品率.4、掷一枚均匀硬币直到出现3次正面才停止. (1)问正好在第6次停止的概率;(2)问正好在第6次停止的情况下,第5次也是出现正面的概率.1、设随机变量X 的分布律为{},1,2,,aP X k k N N===L ,则 a =.2、设随机变量X 的分布函数为0,,0)(≥<−⎩⎨⎧=−x x e A x F x. (12)P 则= A ;ξ<≤= . 3、设随机变量X 的密度函数为()⎪⎩⎪⎨⎧≤≤⎟⎠⎞⎜⎝⎛−=其它0211122x x x f 则X 的分布函数 ()=x F .4、设,且,则 ~()X P λ(1)(2P X P X ===)(1)P X ≥=, 2(03)P X <<=.5、设随机变量X 的分布律为X 1− 1 02k p 2.01.03.04.0 则随机变量2X Y =的分布律为 . 6、设随机变量X 的分布函数为0,10.4,11()()0.8,131,3x x F x P X x x x <−⎧⎪−≤<⎪=≤=⎨≤<⎪⎪≥⎩则X 的分布律为 .二、计算题(请写在背面并标好题序,如果背面不够写请写在信笺上并装订好再交上来)1、设在15只同类型零件中有2只为次品,在其中取3次,每次任取1只,作不放回抽样,以X 表示取出的次品个数,求:(1)X 的分布律; (2)X 的分布函数;(3){}133,1,1,12222P X P X P X P X ⎧⎫⎧⎫⎧⎫≤<≤≤≤<⎨⎬⎨⎬⎨⎬⎩⎭⎩⎭⎩⎭<.2、有一繁忙的汽车站,每天有大量汽车通过,设每辆车在一天的某时段出事故的概率为0.0001,在某天的该时段内有1000辆汽车通过,问出事故的次数不小于2的概率是多少? 3、已知随机变量X 的密度函数为,求: ||(),x f x Ae x −=−∞<<+∞(1)A 值; (2){}0P X <<1; (3)随机变量X 的分布函数. ()F x 4、设()2~3,2X N (1)求{}{}{}{}25,410,2,P X P X P X P X <≤−<≤>>3; (2)确定使c {}{}P X c P X c >=≤. 5、设()~0,1X N (1)求的概率密度; (2)求的概率密度.X Y e =22Y X =+11、设二维离散型随机变量(,)X Y 具有概率分布律\36912151810.010.030.020.010.050.0620.020.020.010.050.030.0730.050.040.030.010.020.0340.030.090.060.150.090.02X Y则X 和Y 的边缘分布律分别为和 . 2、设二维随机变量(,)X Y 的联合分布函数为42(1)(1),0,0(,)0,x y e e x y F x y −−⎧−−>=⎨⎩其他> 则二维随机变量(,)X Y 的联合概率密度函数为 . 3、设二维随机变量(,)X Y 的联合分布函数为ππsin sin ,0,0(,)220,x y x y F x y ⎧≤≤≤≤⎪=⎨⎪⎩其他则πππ0,463P x y ⎧<≤<≤=⎨⎩⎭⎫⎬ .4、设随机变量,X Y 的概率密度分别为, 0,0(),0X xx f x e x −<⎧=⎨≥⎩0,0(),0Y y y f y e y −<⎧=⎨≥⎩且,X Y 相互独立,则二维随机变量(,)X Y 的联合密度函数为 . 5、设随机变量,X Y 相互独立,且,则()(~100,0.2,~50,0.2X b Y b )~Z X Y =+ . 二、计算题1、设随机变量(,)X Y 的概率密度为(6),02,24(,)0,k x y x y f x y −−<<<<⎧=⎨⎩其他(1)确定常数;(2)求k {}1,3P X Y ≤≤;(3)求{}1.5P X ≤;(4)求{}4P X Y +≤;(5)求边缘概率密度. 2、袋中有五个号码1从中任取三个,记其中最小的号码为,2,3,4,5.X ,最大的号码为Y .(1)求X 与Y 的联合概率分布;(2)X 与是否相互独立? Y 3、设某种型号的电子管的寿命(以小时计)近似地服从.随机地选取4只,求其中没有一只寿命小于180的概率.(160,400)N1、设随机变量,,X Y Z 相互独立,且,则 ()5,()11,()8E X E Y E Z ===(231)E X Y ++=; (4)E YZ X −=.2、设随机变量X 的分布律为123101X P p p p −且已知,则2()0.1,()0.9E X E X ==1p = ;2p = ;3p = .3、设随机变量(,)X Y 的概率密度为2,01,0(,)0,x y xf x y <<<<⎧=⎨⎩其他 则 ()E XY =.4、已知,,则= 4.1)(=X E 24)(2X E .0)(=X D .5、已知随机变量,且~(2)X P 22−=X Z ,则 _________, _________.()=Z E ()D Z =6、设连续型随机变量X 的密度函数为()1221−+−=x x e x f π()+∞<<∞−x 则_______.()=X D 7、设随机变量~(,)X U a b ,且,()3=X E ()34=X D ,则 a =, b =. 8、已知随机变量,则~(,),()12,()8X b n p E X D X ==p = ; n =.9、设且~(1,9),~(2,4)X N Y N ,X Y 相互独立,则(23)E X Y −= , (23)D X Y −=. 10、已知,则 ()2,()3,cov(,)1D X D Y X Y ===−cov(321,42)X Y X Y −++−=. 二、计算题1、设(,)X Y 的联合分布律为\1010.20.10.120.100.1300.30.1X Y −1求(),(),Y E X E Y E X ⎛⎞⎜⎟⎝⎠.2、设随机变量X 的概率密度为,0(,)2,120,1x x f x y x x ≤<⎧⎪=−≤≤⎨⎪⎩其他求.(),()E X D X 3、设随机变量,X Y 的概率密度分别为22e ,0()0,0x X x f x x −⎧>=⎨≤⎩ 44e ,0()0,0y Y y f y y −⎧>=⎨≤⎩求.2(),(23E X Y E X Y +−)4、设二维随机变量(,)X Y 在以为顶点的三角形区域上服从均匀分布,求(0,0),(0,1),(1,0)cov(,)X Y 和XY ρ.1、如果12,,,n X X X L 为来自总体X 的样本,X 的分布函数为,则()F x 12,,,n X X X L 的联合分布函数为 ;如果X 的概率密度为()f x ,则12,,,n X X X L 的联合概率密度为 .2、设1234,,,,5X X X X X 是来自总体的样本,则 (0,1)X N ∼521i i Z X ==∑∼.3、设2~(,)X N μσ,12,,,n X X X L 为来自总体X 的样本,则()E X = ,()D X = .4、设12~(),,,,n X P X X X λL 为来自总体X 的样本,则()E X =______,()D X =______.5、如果,且它们相互独立,则2(4),(5)X Y χχ∼∼2X Y +∼ .6、如果,则 2(10)X χ∼()E X =; ()D X =.7、 20.025(30)χ=, 20.05(61)χ=.8、设,且2(0,1),~(100)X N Y χ∼,X Y 相互独立,则统计量~t = .9、 0.01(20)t =, 0.25(50)t =.10、设,且相互独立,则统计量22(20),(30)U V χχ∼∼,U V 3~2UF V= .11、 0.05(9,12)F =,则 ()0.9512,9F =.12、设12,,,n X X X L 相互独立,2(,)i i X N iμσ∼,则 1ni i i a X η==∑∼.X 13、设12,,,n X X X L 是来自正态总体2~(,)X N μσ的样本,则X ∼ ,. ∼14、设12,,,n X X X L 相互独立,(,则 0i X N ∼,1)21~ni i T X ==∑.15、设两个随机变量X 与Y 相互独立,并且2(0,1),()X N Y n χ∼∼,则T =∼ . 二、计算题1、设总体,从总体X 中抽取一个容量为100的样本,求样本均值与总体均值之差的绝对值大于3的概率. (2~60,15X N )2、从一正态总体中抽取容量为10的样本,假定有2%的样本均值与总体均值之差的绝对值在4以上,求总体的标准差.1、设总体()~1,X b p ,12,,,n X X X L 是从总体X 中抽取的一个样本,则参数p 的矩估计量为=pˆ___________. 2、设总体X 的分布律为()()22123211X P θθθθ−−10<<θ是未知参数,12,,,n X X X L 是从中抽取的一个样本,则参数θ的矩估计量 =θˆ. 3、设,()~X E λ12,,,n x x x L 为X 的一组样本观察值,则的最大似然估计为 θ. 4、设2~(,)X N μσ,123,,,4X X X X 是来自总体X 的样本,设1230.10.30.24X X X α+++X 是(0)μμ≠的无偏估计,则 α=.5、设总体,根据来自2~(,0.04)X N μX 的容量为16的样本,测得样本均值为10.05x =,则总体均值μ的置信水平为的置信区间为 0.95.二、计算题1、设总体(1),0~()0,x x X f x θθ⎧+<=⎨⎩其他1<其中,1θ>−12,,,n X X X L 是X 的一个样本,求的矩估计量及极大似然估计量.θ2、设某种砖头的抗压强度(2,N )μσ,今随机抽取20块砖头,测得数据如下(kg.cm -2):64 69 49 92 55 97 41 84 88 99 84 66 100 98 72 74 87 84 48 81 (1)求μ的置信概率为0.95的置信区间. (2)求的置信概率为0.95的置信区间. 2σ三、证明题设X 1,X 2是从正态总体(2,N )μσ中抽取的样本 11221231211311ˆˆˆ;;3344222;X X X X X μμμ=+=+=+X 试证123ˆˆˆ,,μμμ都是μ的无偏估计量,并求出每一估计量的方差.1、已知某炼铁厂的铁水含碳量在正常情况下服从正态分布.现在测了5炉铁水,其含碳量(%)分别为()2~ 4.55,1.108X N 4.28 4.40 4.42 4.35 4.37问若标准差不改变,总体平均值有无显著性变化(=0.05)? α2、某种矿砂的5个样品中的含镍量(%)经测定为: 3.24 3.26 3.24 3.27 3.25设含镍量服从正态分布,问在=0.01下能否接收假设:这批矿砂的含镍量为3.25. α3、测量某种溶液中的水分,从它的10个测定值得出()()0.452%,0.037%x s ==.设测定值总体为正态,μ为总体均值,为总体标准差,试在水平下检验.σ0.05α=(1).()()01:0.5%;:0.5%H H μμ=<(2)()0:0.04%H σ=;()1:0.04%H σ<.。
概率论上机实验报告【范本模板】
《概率论与数理统计应用》实验报告班级:学号:姓名:实验目的:a.熟悉MATLAB的在概率计算方面的操作;b.掌握绘制常见分布的概率密度及分布函数图形等命令; c.会用MABLAB求解关于概率论与数理统计的实际应用题d.提高数据分析的能力实验题目与解答:1。
二项分布的泊松分布与正态分布的逼近设X ~B(n,p),其中np=21) 对n=101,…,105,讨论用泊松分布逼近二项分布的误差. 画处逼近的图形2) 对n=101,…,105, 计算 )505(≤<X P ,)9020(≤<X P 1)用二项分布计算 2)用泊松分布计算 3)用正态分布计算比较用泊松分布逼近与正态分布逼近二项分布的优劣。
问题分析:查询MATLAB 函数库可知泊松分布概率密度函数为(),poissdpf k lambda ,泊松分布概率函数为(),poisscpf k lambda .其中(k)(k,)!(k,)!k ifloor i poisspdf e k poisscdf e i λλλλλλ--===∑同时,二项分布概率密度函数为(),,binopdf x n p ,二项分布概率分布函数为(),,binocdf x n p 。
其中()()()()()()()()()0,1,,n 0,1,,n 0,,,,1n x x xn i i i n binopdf x n p p q I x x n binocdf x n p p p I i i --=⎛⎫= ⎪⎝⎭⎛⎫=- ⎪⎝⎭∑正态分布概率分布函数为(),,normcdf x μσ,其中()()222,,x normcdf x μσμσ--=利用,poissdpf binopdf 这两个函数,即可画出泊松分布和二项分布的概率密度曲线,设置变量Er 表示在每一点处,poissdpf binopdf 概率密度差值的绝对值,对Er 求平均值Aver ,并计算方差Var 。
概率论与数理统计上机实习题目
概率论与数理统计上机实习作业院系:**学院 班级:**班 姓名:** 学号:**一、 某人写了n 封信,又写了n 个信封,然后将这n 封信随机的装入这n 个信封中,用p n 表示至少有一封信装对的概率。
1.编制程序,用随机数模拟至少20000次,求当= 10时,p n 的值2.重复第一步,画出n=2,3,...,50时,p n 的散点图。
答:先在matlab 中编辑3个程序:rrank.m ,rrans.m ,thms.m (程序略)1. 输入 rrans(10,20000),按回车键。
输出结果如下: rrans(10,20000)ans =0.63385000000000(其中ans 为p n 的值)2. 输入thms(20000,2,1,50),结果如图所示:.二、 设X 1,X 2,…,Xn 相互独立且都服从区间[0,1]上的均匀分布,f(x)为区间[0,1]上的一个可积函数,由大数定律可知i ni X n 11=∑依概率收敛于⎰10)()(dx x f X Ef i .编制程序,用随机数模拟至少40000次,近似地求下列两个积分的值:dx e x ⎰102, dx x x ⎰10sin 答:运行结果如下:>> rand(200,200);>> fun=inline('exp(x.*x)','x');>> Isim=quad(fun,0,1)Isim =1.4627>> rand(200,200);>> fun=inline('sin(x)./(x+eps)','x');>> Isim=quad(fun,0,1)Isim =0.9461。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
概率论与数理统计上机
实习作业
制作人:
学号:
班级
1.假设车辆牌号的后两位数都是00~99中的某一个,任取m辆车(1060
≤≤),编制程序求出这m辆车中至少有两辆车后两位数相同m
的概率
P,画出m=10,11,…..,60时,m P的散点图。
m
MATLAB语言:
m=10:60;
p=zeros(size(m));
for i=1:length(m)
N=10000;
n=0;
for j=1:N
x=floor(100*rand(m(i),1));
v=zeros(100,1);
for k=1:m(i)
v(x(k)+1)=v(x(k)+1)+1;
end
a=sum(v>1);
if a>0
n=n+1;
end
end
p(i)=n/N;
end
scatter(m,p);
调试结果:
2.某人写了n封信,又写了n个信封,然后将这n封信随机地装入这n个信封中,用
P表示至少有一封信装对的概率。
n
(1)编制程序,用随机数模拟20000次,求当10
n=时,
P的值。
n (2)重复第一步,画出2,3,...,50
P的散点图。
n=时,
n
MATLAB语言:
t=20000;
q=2:50;
m=zeros(1,49);
for j=1:49
for i=1:t
k1=randperm(q(j));
k2=randperm(q(j));
k=k1-k2;
if sum(k==0)>=1
m(j)=m(j)+1;
end
end
end
h=m./t;
h(9)
scatter(q,m);
调试结果:
值为n=10时的Pn值。
3.设12,,...,n X X X 相互独立,且都服从区间[0,1]上的均匀分布,()f x 为
区间[0,1]上的一个可积函数,由大数定律可知1
1n
i i X n =∑依概率收敛于
1
()()i Ef X f x dx =⎰。
编制程序,用随机数模拟至少40000次,近似地求
下列两个积分的值。
(1)2
1
0x e dx ⎰ (2)1
0sin x
dx x
⎰ MATLAB 语言:
n=1000000; m=zeros(1,n); t=zeros(1,n); for i=1:n k=rand(1); m(i)=exp(k^2); t(i)=sin(k)/k; end p=sum(m)/n q=sum(t)/n
调试结果:
P 为2
1
0x e dx ⎰ 的值,q 为1
sin x
dx x
⎰的值。
4.假设男,女婴的出生率均为0.5,每个家庭只要有一个男婴出生就不再生下一胎,如果没有男婴,不论已经出生多少女婴,都可以继续要下一胎.假设某地有100万个家庭,按照上述生育政策,自行编制程序回答下一问题:
(1) 每个家庭大约有几个小孩? (2) 男女比例大约为多少?
(3)如果男婴与女婴的出生率为0.51与0.49,重新回答前两问。
MATLAB 语言:
max=1000000; m=zeros(max,1); n=zeros(max,1); p=zeros(max,1); q=zeros(max,1); while sum(m>0)<max for i=1:max k=rand(1); if m(i)==0 if k>=0.5 m(i)=m(i)+1; else n(i)=n(i)+1; end end
if p(i)==0
if k>=0.51
p(i)=p(i)+1;
else
q(i)=q(i)+1;
end
end
end
end
t=sum(m+n)/max;
l=sum(m)/sum(n);
w=sum(p+q)/max;
e=sum(p)/sum(q);
调试结果:
t为第一种情况每个家庭孩子的平均数l为第一种情况男女比列
w为第二种情况每个家庭孩子的平均数e为第二种情况男女比列。