概率论matlab实验报告参考模板.docx
- 文档编号:466675
- 上传时间:2022-10-10
- 格式:DOCX
- 页数:31
- 大小:451.94KB
概率论matlab实验报告参考模板.docx
《概率论matlab实验报告参考模板.docx》由会员分享,可在线阅读,更多相关《概率论matlab实验报告参考模板.docx(31页珍藏版)》请在冰豆网上搜索。
概率论matlab实验报告参考模板
概率论与数理统计matlab上机
实验报告
班级:
学号:
姓名:
指导老师:
实验一常见分布的概率密度、分布函数生成
[实验目的]
1.会利用MATLAB软件计算离散型随机变量的概率,连续型随机变量概率密度值。
2.会利用MATLAB软件计算分布函数值,或计算形如事件{X≤x}的概率。
3.会求上α分位点以及分布函数的反函数值。
[实验要求]
1.掌握常见分布的分布律和概率密度的产生命令,如binopdf,normpdf
2.掌握常见分布的分布函数命令,如binocdf,normcdf
3.掌握常见分布的分布函数反函数命令,如binoinv,norminv
[实验内容]
常见分布的概率密度、分布函数生成,自设参数
1、X~B(20,0.4)
(1)P{恰好发生8次}=P{X=8}
(2)P{至多发生8次}=P{X<=8}
(1)binopdf(8,20,0.4)
ans=
0.1797
(2)binocdf(8,20,0.4)
ans=
0.5956
2、X~P
(2)
求P{X=4}
poisspdf(4,2)
ans=
0.0902
3、X~U[3,8]
(1)X=5的概率密度
(2)P{X<=6}
(1)unifpdf(5,3,8)
ans=
0.2000
(2)unifcdf(6,3,8)
ans=
0.6000
4、X~exp(3)
(1)X=0,1,2,3,4,5,6,7,8时的概率密度
(2)P{X<=8}
注意:
exp(3)与教材中参数不同,倒数关系
(1)exppdf(0:
8,3)
ans=
Columns1through3
0.33330.23880.1711
Columns4through6
0.12260.08790.0630
Columns7through9
0.04510.03230.0232
(2)expcdf(8,3)
ans=
0.9305
5、X~N(8,9)
(1)X=3,4,5,6,7,8,9时的概率密度值
(2)X=3,4,5,6,7,8,9时的分布函数值
(3)若P{X<=x}=0.625,求x
(4)求标准正态分布的上0.025分位数
(1)normpdf(3:
9,8,3)
ans=
Columns1through3
0.03320.05470.0807
Columns4through6
0.10650.12580.1330
Column7
0.1258
(2)normcdf(3:
9,8,3)
ans=
Columns1through3
0.04780.09120.1587
Columns4through6
0.25250.36940.5000
Column7
0.6306
(3)norminv(0.625,8,3)
ans=
8.9559
(4)norminv(0.975,0,1)
ans=
1.9600
6、X~t(3)
(1)X=-3,-2,-1,0,1,2,3时的概率密度值
(2)X=-3,-2,-1,0,1,2,3时的分布函数值
(3)若P{X<=x}=0.625,求x
(4)求t分布的上0.025分位数
(1)tpdf(-3:
3,3)
ans=
Columns1through3
0.02300.06750.2067
Columns4through6
0.36760.20670.0675
Column7
0.0230
(2)tcdf(-3:
3,3)
ans=
Columns1through3
0.02880.06970.1955
Columns4through6
0.50000.80450.9303
Column7
0.9712
(3)tinv(0.625,3)
ans=
0.3492
(4)tinv(0.975,3)
ans=
3.1824
7、X~卡方(4)
(1)X=0,1,2,3,4,5,6时的概率密度值
(2)X=0,1,2,3,4,5,6时的分布函数值
(3)若P{X<=x}=0.625,求x
(4)求卡方分布的上0.025分位数
(1)chi2pdf(0:
6,4)
ans=
Columns1through3
00.15160.1839
Columns4through6
0.16730.13530.1026
Column7
0.0747
(2)chi2cdf(0:
6,4)
ans=
Columns1through3
00.09020.2642
Columns4through6
0.44220.59400.7127
Column7
0.8009
(3)chi2inv(0.625,4)
ans=
4.2361
(4)chi2inv(0.975,4)
ans=
11.1433
8、X~F(4,9)
(1)X=0,1,2,3,4,5,6时的概率密度值
(2)X=0,1,2,3,4,5,6时的分布函数值
(3)若P{X<=x}=0.625,求x
(4)求F分布的上0.025分位数
(1)fpdf(0:
6,4,9)
ans=
Columns1through3
00.44790.1566
Columns4through6
0.05950.02550.0122
Column7
0.0063
(2)fcdf(0:
6,4,9)
ans=
Columns1through3
00.54420.8218
Columns4through6
0.92110.96090.9788
Column7
0.9877
(3)finv(0.625,4,9)
ans=
1.1994
(4)finv(0.975,4,9)
ans=
4.7181
实验二概率作图
[实验目的]
1.熟练掌握MATLAB软件的关于概率分布作图的基本操作
2.会进行常用的概率密度函数和分布函数的作图
3.会画出分布律图形
[实验要求]
1.掌握MATLAB画图命令plot
2.掌握常见分布的概率密度图像和分布函数图像的画法
[实验内容]
任选四种分布,自设参数(已画八种分布图像,可熟悉各分布特点)
1、X~B(20,0.4)
代码:
x=0:
20;
y=binopdf(x,20,0.4)
plot(x,y,'.')
结果:
2、X~exp(3)
概率密度图像
代码:
x=0:
0.01:
15;
y=exppdf(x,3)
plot(x,y)
结果:
分布函数
代码:
x=-1:
0.01:
15;
y=expcdf(x,3)
plot(x,y)
结果:
3、X~P(4)
概率密度图形
代码:
x=0:
10;
y=poisspdf(x,4)
plot(x,y,'.')
结果:
分布函数图形
代码:
x=0:
0.01:
10;
y=poisscdf(x,4)
plot(x,y)
结果:
4、X~U(3,8)
概率密度图形
代码:
x=0:
0.01:
10;
y=unifpdf(x,3,8)
plot(x,y,'.')
结果:
分布函数图形
代码:
x=0:
0.01:
10;
y=unifcdf(x,3,8)
plot(x,y)
结果:
5、X~N(4,9)
概率密度图形
代码:
x=-10:
0.01:
18;
y=normpdf(x,4,3);
plot(x,y)
结果:
分布函数图形
代码:
x=-10:
0.01:
18;
y=normcdf(x,4,3);
plot(x,y)
结果:
同一坐标系,均值是4,标准差分别为1,2,3的正态分布概率密度图形
代码:
x=-5:
0.01:
15;
y1=normpdf(x,4,1);
y2=normpdf(x,4,2);
y3=normpdf(x,4,3);
plot(x,y1,x,y2,x,y3)
结果:
6、X~t(3)
概率密度图形
代码:
x=-10:
0.01:
10;
y=tpdf(x,3);
plot(x,y)
结果:
分布函数图形
代码:
x=-10:
0.01:
10;
y=tcdf(x,3);
plot(x,y)
结果:
7、X~卡方(4)
概率密度图形
代码:
x=0:
0.01:
15;
y=chi2pdf(x,4);
plot(x,y)
结果:
分布函数图形
代码:
x=0:
0.01:
15;
y=chi2cdf(x,4);
plot(x,y)
结果:
8、X~F(4,9)
概率密度图形
代码:
x=0:
0.001:
10;
y=fpdf(x,4,9);
plot(x,y)
结果:
分布函数图形
代码:
x=0:
0.001:
10;
y=fcdf(x,4,9);
plot(x,y)
结果:
实验三数字特征
[实验目的]
1加深对数学期望,方差的理解
2理解数学期望,方差的意义,以及具体的应用
3加深对协方差,相关系数的理解
4了解协方差,相关系数的具体的应用
[实验要求]
1概率与频率的理论知识,MATLAB软件
2协方差,相关系数的理论知识,MATLAB命令cov,corrcoef
[实验内容]
P101-11
代码:
exp=[];
price=[-200100];
exp
(1)=expcdf(1,4)
exp
(2)=1-exp
(1)
Ey=exp*price'
结果:
exp=
0.2212
exp=
0.22120.7788
Ey=
33.6402
即平均获利为Ey=e^(-1/4)*300-200=33.6402
p101-13
代码:
Symsxy
fxy=(x+y)/3;
Ex=int(int(fxy*x,y,0,1),x,0,2)
Ey=int(int(fxy*y,y,0,1),x,0,2)
Exy=int(int(fxy*x*y,y,0,1),x,0,2)
E=int(int(fxy*(x^2+y^2),y,0,1),x,0,2)
结果:
Ex=
11/9
Ey=
5/9
Exy=
2/3
E=
13/6
>>
P102-22
代码:
Symsxy
fxy=1;
Ex=int(int(fxy*x,y,-x,x),x,0,1)
Ey=int(int(fxy*y,y,-x,x),x,0,1)
Ex2=int(int(fxy*x^2,y,-x,x),x,0,1)
Ey2=int(int(fxy*y^2,y,-x,x),x,0,1)
Dx=Ex2-Ex^2
Dy=Ey2-Ey^2
结果:
Ex=
2/3
Ey=
0
Ex2=
1/2
Ey2=
1/6
Dx=
1/18
Dy=
1/6
>>
P103-26
代码:
Symsxy
fxy=2-x-y;
Ex=int(int(fxy
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 概率论 matlab 实验 报告 参考 模板
![提示](https://static.bdocx.com/images/bang_tan.gif)