[生物学]matlab在科学计算中的应用8.ppt

  1. 1、本文档共113页,可阅读全部内容。
  2. 2、有哪些信誉好的足球投注网站(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
[生物学]matlab在科学计算中的应用8

第8章 概率论与数理统计问题的求解 概率分布与伪随机数生成 统计量分析 数理统计分析方法及计算机实现 统计假设检验 方差分析及计算机求解 8.1概率分布与伪随机数生成 8.1.1 概率密度函数与分布函数概述 通用函数计算概率密度函数值 函数 pdf 格式 P=pdf(‘name’,K,A) P=pdf(‘name’,K,A,B) P=pdf(‘name’,K,A,B,C) 说明 返回在X=K处、参数为A、B、C的概率密度值,对于不同的分布,参数个数是不同;name为分布函数名。 例如二项分布:设一次试验,事件Y发生的概率为p,那么,在n次独立重复试验中,事件Y恰好发生K次的概率P_K为:P_K=P{X=K}=pdf(bino,K,n,p) 例: 计算正态分布N(0,1)的随机变量X在点0.6578的密度函数值。 解: pdf(norm,0.6578,0,1) ans = 0.3213 例:自由度为8的卡方分布,在点2.18处的密度函数值。 解: pdf(chi2,2.18,8) ans = 0.0363 随机变量的累积概率值(分布函数值) 通用函数cdf用来计算随机变量的概率之和(累积概率值) 函数 cdf 格式 cdf(‘name’,K,A) cdf(‘name’,K,A,B) cdf(‘name’,K,A,B,C) 说明 返回以name为分布、随机变量X≤K的概率之和的累积概率值,name为分布函数名. 例: 求标准正态分布随机变量X落在区间(-∞,0.4)内的概率。 解: cdf(norm,0.4,0,1) ans = 0.6554 例:求自由度为16的卡方分布随机变量落在[0,6.91]内的概率。 解: cdf(chi2,6.91,16) ans = 0.0250 随机变量的逆累积分布函数 MATLAB中的逆累积分布函数是已知,求x。 命令 icdf 计算逆累积分布函数 格式 icdf(‘name’,F,A) icdf(‘name’,F,A,B) icdf(‘name’,F,A,B,C) 说明 返回分布为name,参数为a1,a2,a3,累积概率值为F的临界值,这里name与前面相同。 如果F= cdf(‘name’,X,A,B,C) , 则 X = icdf(‘name’,F,A,B,C) 例:在标准正态分布表中,若已知F=0.6554,求X 解: icdf(norm,0.6554,0,1) ans = 0.3999 例:公共汽车门的高度是按成年男子与车门顶碰头的机会不超过1%设计的。设男子身高X(单位:cm)服从正态分布N(175,6),求车门的最低高度。 解:设h为车门高度,X为身高。 求满足条件 F{Xh}=0.99,即 F{Xh}=0.01故 h=icdf(norm,0.99, 175, 6) h = 188.9581 8.1.2 常见分布的概率密度函数与分布函数 8.1.2.1 Poisson分布 例:绘制 l =1,2,5,10 时 Poisson 分布的概率密度函数与概率分布函数曲线。 x=[0:15]; y1=[]; y2=[]; lam1=[1,2,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) 8.1.2.2 正态分布 例: x=[-5:.02:5]; y1=[]; y2=[]; mu1=[-1,0,0,0,1]; sig1=[1,0.1,1,10,1]; sig1=sqrt(sig1); for i=1:length(mu1) y1=[y1,normpdf(x,mu1(i),sig1(i))]; y2=[y2,normcdf(x,mu1(i),sig1(i))]; end plot(x,y1), figure; plot(x,y2) 8.1.2.3 分布 例:绘制 为(1,1),(1,0.5),(2,1),(1,2),(3,1)时 x=[-0.5:.02:5]‘; %x=[-eps:

文档评论(0)

skvdnd51 + 关注
实名认证
内容提供者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档