高等数学-read.ppt

高等数学 一、数据的录入、保存和调用； 二、 Matlab统计工具箱中的基本统计命令； 三、用Matlab软件求解概率统计问题 教学内容 一、数据的录入、保存和调用 例1上海市社会商品零售总额(单位：亿元）和全民所有制职工工资总额(单位：亿元）的数据如下： 175.0 155.0 137.4 95.9 77.5 68.7 67.9 61.7 51.8 41.4 零售总额 73.4 63.8 52.8 43.2 34.9 33.7 32.4 31.6 27.6 23.8 工资总额 87 86 85 84 83 82 81 80 79 78 年份 t=78:87; %输入年份的数据 x=[23.8 27.6 31.6 32.4 33.7 34.9 43.2 52.8 63.8 73.4];%输入工资 y=[41.4 51.8 61.7 67.9 68.7 77.5 95.9 137.4 155.0 175.0];%输入零售总额 save data t x y %将变量t,x,y的数据保存在文件data中 load data 进行统计分析时，用以下命令调用数据文件data中的数据 load data y(2) ans = 51.8000 方法2：用矩阵形式输入和保存数据 1、输入 data=[78 79 80 81 82 83 84 85 86 87; 23.8 27.6 31.6 32.4 33.7 34.9 43.2 52.8 63.8 73.4; 41.4 51.8 61.7 67.9 68.7 77.5 95.9 137.4 155.0 175.0] 2、将矩阵data的数据保存在文件data1中 Save data1 data 3、进行统计分析时，先用命令： load data1 t=data(1,:) %矩阵data的第一行的数据赋给变量t x=data(2,:) %矩阵data的第一行的数据赋给变量x y=data(3,:) %矩阵data的第一行的数据赋给变量y 二、计算基本统计量的命令： kurtosis skewness var std median mean 命令 峰度 偏度 方差 标准差 中位数 均值 功能 三、计算常见概率分布的命令 泊松分布 poiss 分布 chi2 二项分布 bino 正态分布 norm F分布 f 指数分布 exp t分布 t 均匀分布 unif 分布 命令 分布 命令 工具箱对每一种分步都提供了五类函数，其命令字符为 随机数生成 rnd 均值与方差 stat 逆概率分布 inv 分布函数 cdf 概率密度 pdf 功能 字符 当需要一种分布的某一类函数时，将以上的分布命令 字符和函数命令字符接起来，并输入自变量和参数即可。 例2画出正态分布N(0,1)和N(0,22)的概率密度函数图像。 输入命令： x=-6:0.01:6; y=normpdf(x); z=normpdf(x,0,2); plot(x,y,x,z) 1、密度函数：p=normpdf(x,mu,sigma) （当mu=0,sigma=1时可以却省） 2、概率分布：P=normcdf(x,mu,sigma ) 例3 计算标准正态分布的概率P{-1x1} 命令： P=normcdf(1)-normcdf(-1) 结果：P = 0.6827 例4 计算x~N(2,32)的概率P{-1x1} 命令： P=normcdf(1,2,3)-normcdf(-1,2,3) 结果：P = 0.2108 3、逆概率分布x=norm(P,mu,sigma) 即求出x，使得P(x1)=P,此命令可用来求分位数。 例5 取α=0.05，求 α=0.05时，P=0.975 norminv(0.975)=1.96 例6 计算均值100，标准差10的正态分布随机变量的 概率0.01的分位数。 命令：norminv(0.01,100,10) 结果： ans = 76.7365 4、均值与方差：[m,v]=normstat(mu,sigma) 例7 求正态分布N(3,52)的均值与方差。 命令：[m,v]=normstat(3,5) 结果：m=3,v=25 5、随机数生成：normrnd(mu,sigma,m,n) 产生m×n阶的正态分布随机数生成。 例8 ：命令：m=normrnd([1 2 3;4 5 6],0.1,2,3) 结果：m = 0.9567 2.0125 2.8854 3.8334 5.0288 6.1191 注意：以上仅仅以正态分布为例，其他分布