概率论与数理统计上机作业
概率论与数理统计上机实验 (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)讨论这个研究的推广和你认为有用的其他分析。
- 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为第二种情况男女比列。