数学实验ly7概率2015.12剖析.ppt

* * * * 数据的统计分析 随机试验 参数估计 假设检验 随机试验 古典概率： 事件A发生的概率 在100人的团体中，如果不考虑年龄的差异，研究是否有 两个以上的人生日相同。假设每人的生日在一年365天中 的任意一天是等可能的，那么随机找 个人（365） 问这些人生日各不相同的概率是多少？ 至少有两个人生日相同的概率为多少？ P169 for n=1:100 p0(n)=prod(365:-1:365-n+1)/365^n; p1(n)=1-p0(n); end n=1:100; plot(n,p0,n,p1,--) xlabel(‘人数’),ylabel(‘概率’) legend(‘生日各不相同的概率’,‘至少两人相同的概率’) axis([0 100 -0.1 1.1]),grid on 重要的概率分布 ? 分布 ? ? 离散型随机变量 连续型随机变量 均匀分布 二项分布 泊松分布 几何分布 均匀分布 指数分布 正态分布 分布 t分布 ?F分布 字符 unid bino poiss geo unif exp norm chi2 t f 功能 概率密度 分布函数 逆概率分布 均值与方差 随机数生成 字符 pdf cdf inv stat rnd 表8-2 概率分布的命令字符 表8-3 运算功能的命令字符 正态分布 概率密度为 y=normpdf(x, mu, sigma) 的正态分布的密度函数 y=normpdf(x) 标准正态分布的密度函数 y=normcdf(x, mu, sigma) 的正态分布的分布函数 y=normcdf(x) 标准正态分布的分布函数 x=-6:0.01:6; y1=normpdf(x); z1=normcdf(x); y2=normpdf(x, 0, 2); z2=normcdf(x, 0, 2); subplot(1, 2, 1),plot(x, y1, x, y2); legend(N(0,1),N(0,2^2)) subplot(1, 2, 2),plot(x, z1, x, z2); legend(N(0,1),N(0,2^2)) 例：某厂生产一种设备，其平均寿命为10年，标准差为2年.如该设备的寿命服从正态分布，求寿命不低于9年的设备占整批设备的比例？ 解： 设随机变量 为设备寿命，由题意 clear p1=normcdf(9,10,2) p1 =0.3085 1-p1 ans = 0.6915 统计作图 在数据较小、较少的情况下输入 Matlab交互环境 境下输入 M文件的形式输入 数据 数据量较大，且不以计算机可读 形式存在 load *.M 读数据文件的命令 读入 load *.txt P180 函数名称 功能简介 Max(x) 求最大值 Min(x) 求最小值 Median(x) 求中值 Range(x) 求极差 Mean(x) 求算术平均值 Std(x) 求样本标准差 Var(x) 求样本方差 Cov(x) 求协方差矩阵 基本统计函数表 例8-19 某班（共有120名学生）的高等数学成绩如下： 74 63 78 76 89 56 70 97 89 94 76 88 65 83 72 41 39 72 73 68 14 76 45 70 90 46 54 61 75 76 49 57 78 66 64 74 78 87 86 73 47 67 21 66 79 67 68 65 56 84 66 73 68 72 76 65 70 94 53 65 77 78 53 74 59 50 98 67 89 78 63 92 54 87 84 80 63 64 85 66 69 69 60 54 75 33 30 62 74 65 84 73 55 85 75 76 81 71 83 72 56 84 76 75 67 65 35 94 59 47 45 67 75 36 78 82 94 70 84 75 根据以上数据作出该门课程成绩的频数表和直方图。 解：（1）数据输入： 方法1：在Matlab的交互环境下直接输入; 方法2：将以上数据以一列的形式存为A.txt文件,用 (2) 用hist命令作频数表和直方图：(区间个数为5，可省略) [N,X]=hist(A,5) 120名学生高数成绩的频数表； hist(A,5) 120名学生高数成绩的直方图； load A.txt 命令读入数据。 matleb命令： load A.txt↙ disp(‘高数成绩的频数表’),[N,X]=hist(A,5)%N为频数 频数是如何计算的？ a1=min(A);a2=max(