MATLAB期末复习课件.pptVIP

记 Y = X1 + X2 + X3 + X4 + X5 + X6 随机变量 Y 服从二项分布 Y ~ B( n, p ) k =0,1,2,…,6 二项分布概率函数: binopdf(x,n,p) x 是n重贝努里试验中事件A出现的次数. %计算Galton试验板分布律 n=6;x=0:n; Y=binopdf(x,n,0.5) bar(x,Y) ans = 0.0156 0.0938 0.2344 0.3125 0.2344 0.0938 0.0156 2.二项分布 function Y=mlab72(n) x=0:n; Y=binopdf(x,n,0.5); bar(x,Y),hold on segm2=n/4; t=0:0.1:n; y=exp(-0.5*((t-n/2).^2)/segm2)/sqrt(2*pi*segm2); plot(t,y,r) 二项分布与正态分布近似关系实验 n=6 n=8 例3.已知在2000个产品中有30个次品,现抽查100个，求次品数X的分布律，计算抽查不超过2个次品的概率 (1)不放回抽样: X服从超几何分布 (2)放回抽样: X服从二项分布 次品率按30/2000 =0.015计算 命令文件　cytest.m x=0:13; %随机变量X的可能取值 P1=hygepdf(x,2000,30,100), %超几何分布的分布律 P2=binopdf(x,100,0.015), %二项分布的分布律 subplot(2,1,1);bar(x,P1,‘r’) %绘超几何分布条形图 subplot(2,1,2);bar(x,P2, ‘b’) %绘二项分布条形图 Pg=sum(P1(1:3)) %计算几何分布前三项 Pb=sum(P2(1:3)) %计算二项分布前三项 9/16 (3).产品抽查的概率计算 P1 = 0.2122 0.3402 0.2609 0.1274 0.0445 0.0119 0.0025 0.0004 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 P2 = 0.2206 0.3360 0.2532 0.1260 0.0465 0.0136 0.0033 0.0007 0.0001 0.0000 0.0000 0.0000 0.0000 0.0000 Pg = 0.8132 Pb = 0.8098 10/16 (4).电子器件寿命的参数估计 1100 380 350 190 620 520 450 410 340 400 320 415 305 300 290 280 310 140 例4．电子管使用寿命 X (单位:小时)服从指数分布. 样本数据如下,试估计 11/16 data=[140 310 280 290 300 305 415 320 400 ... 340 410 450 520 620 190 350 380 1100]; syms x theta f=exp(-x/theta)/theta; theta1=mle(Exponential, data) f=subs(f,theta,theta1) ezplot(f,0,1200) theta1 = 395.5556 分布类型为指数分布选择 ‘Exponential’ 极大似然估计函数命令格式: P=mle(‘dist’,data): Dist——指定分布的名称， Data——样本数据.P——返回极大似然估计值. 12/16 极大似然估计(Maximum Likelihood Estimation) theta1=mle(‘Exponential’, data) Normal 正态分布 Geometric Exponential Binomial dist 几何分布 指数分布 二项分布 分布类型 Weibull Uniform Poisson dist 威布尔分布 均匀分布 泊松分布 分布类型 统计直方图绘制方法hist()的使用格式 N = hist(Y,m) 将Y中的最小数和最大数所成区间分为m个小区间, 返回Y中元素在各小区间中的元素个数. 当“m”省略 时,默认m=10.当“N=”省略时,在图形窗口绘直方图. (5)、定积分