- 1、本文档共18页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
第3章概率统计实分析及MatlAb求解
概率统计实例分析及MatlAb求解
第3章 概率统计实例分析及MatlAb求解 1
3.1 随机变量分布与数字特征实例及MATLAB求解 1
3.1.1 MATLAB实现 1
3.1.2 相关实例求解 2
3.2 数理统计实例分析及MATLAB求解 4
3.1.1 MATLAB实现 4
3.1.2 相关实例求解 4
3.3参数估计与假设检验实例分析及MATLAB求解 5
3.1.1 MATLAB实现 5
3.1.2 相关实例求解 5
3.4 方差分析实例求解 10
3.1.1 MATLAB实现 10
3.1.2 相关实例求解 10
3.5 判别分析应用实例及求解 14
3.1.1 MATLAB实现 14
3.1.2 相关实例求解 14
3.6 聚类分析应用实例及MATLAB求解 16
3.1.1 MATLAB实现 16
3.1.2 相关实例求解 16
3.1 随机变量分布与数字特征实例及MATLAB求解
3.1.1 MATLAB实现
用mvnpdf和mvncdf函数可以计算二维正态分布随机变量在指定位置处的概率和累积分布函数值。
利用MATLAB统计工具箱提供函数,可以比较方便地计算随机变量的分布律(概率密度函数)分布函数中矩阵元素求期望的函数mean和var,若要求整个矩阵所有元素的均方差,则要使用std2函数
伪随机数生成函数:
A=gamrnd(a,lambda,n,m) % 生成n*m的分布的伪随机矩阵
B=raylrnd(b,n,m) %生成rayleigh的伪随机数
3.1.2 相关实例求解
例2-1 计算服从二维正态分布的随机变量在指定范围内的累积分布函数值并绘图。
程序:
%二维正态分布的随机变量在指定范围内的累积分布函数图形
mu=[0 0];
sigma=[0.25 0.3;0.3 1];%协方差阵
x=-3:0.1:3;y=-3:0.2:3;
[x1,y1]=meshgrid(x,y);%将平面区域网格化取值
f=mvncdf([x1(:) y1(:)],mu,sigma);%计算累积分布函数值
F=reshape(f,numel(y),numel(x));%矩阵重塑
surf(x,y,F);
caxis([min(F(:))-0.5*range(F(:)),max(F(:))]);%range(x)表示最大值与最小值的差,即极差。
axis([-3 3 -3 3 0 0.5]);
xlabel(x);
ylabel(y);
zlabel(Probability Density);
图1 二维正太分布累积分布函数值图
例 设的概率密度为求。
解程序syms x
f1=x/1500^2;
f2=(3000-x)/1500^2;
Ex=int(x*f1,0,1500)+int(x*f2,1500,3000)
运行结果:
Ex =1500
Ex =1/3
=0.5,1,3,5,10 时 Poisson 分布的概率密度函数与概率分布函数曲线。
代码如下:
x=[0:15]; y1=[]; y2=[]; lam1=[0.5,1,3,5,10];
for i=1:length(lam1)
y1=[y1,poisspdf(x,lam1(i))]; y2=[y2,poisscdf(x,lam1(i))];
end
plot(x,y1), figure; plot(x,y2)
图2 泊松分布概率密度函数图
图3 泊松分布概率分布函数
3.2 数理统计实例分析及MATLAB求解
3.1.1 MATLAB实现
在MATLAB中各种随机数可以认为是独立同分布的,即简单随机样本。常用分布的随机数产生方法,可用分布英文名称缩写加上rnd,例如:
x=betarnd(a,b,m,n) 参数为a,b的beta分布;
x=binornd(N,p,m,n) 参数为N,p的二项分布;
3.1.2 相关实例求解
例2-4:设总体密度函数
试从该总体中抽取容量为1000的简单随机样本。
解 利用MATLAB编辑窗口保存以下程序,保存为ex11.m
n=1000;
x=zeros(1,n);
k=0;
while kn
a=rand*pi-pi/2;
b=rand/2;
if b(cos(a)/2)
k=k+1;
x(k)=a;
end
end
hist(x,-pi/2:0.2:pi/2)
保存完成之后,在命令窗口执行ex11,则x被赋值,且可以得到这个容量为1000的样本的直方图。
图7 直方图
3.3参数估计与假设检验实例分析及MATLAB求解
3.1.1 MATLAB实现
3.1.2 相
您可能关注的文档
最近下载
- 隔离放大器与仪用放大器.ppt
- 百家讲坛文字稿——《屈原》.doc
- 上虞市临江化工有限公司年产500吨3,4,5-三氟溴苯及2000吨2,4-二氯-5-氟苯乙酮改造项目环境影响报告书.doc VIP
- 基于协同过滤的房源推荐系统的设计与实现.docx VIP
- 临床护理技术操作并发症预防和处理规范.pdf VIP
- T/ZZB 1390-2019 涤纶低弹丝油剂.pdf VIP
- GBT50964-2014 小型水电站运行维护技术规范.docx
- 浙江省杭州市萧山区2022-2023学年四年级下学期数学期末检测试卷.docx VIP
- 太阳能路灯技术规范[共20页].pdf VIP
- FMC6000移动式破碎机操作及维修手册.pdf
文档评论(0)