matlab第七讲.ppt

matlab第七讲.ppt

《matlab第七讲.ppt》由会员分享，可在线阅读，更多相关《matlab第七讲.ppt（45页珍藏版）》请在冰豆网上搜索。

MATLAB的统计分析的统计分析第七讲王文健MATLAB数据处理与应用数据处理与应用2011-2012学年选修课学年选修课西南交通大学摩擦学研究所主要内容点资料图绘制常用统计函数标准差与协方差相关系数随机数概率分布2TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析绘制资料点分布图统计分析中最常用的方法是绘制资料点分布图，通过资料点分布图可以显示数据的分布情况及范围绘制点资料图函数：

bar（）、hist（）、rose（）、stairs（）、stem（）等举例：

x=1:

10;y=rand（size（x）;bar（x,y）;3TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析4TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析统计分析函数max（x）找出x的最大值max（x,y）找出x和y的最大值y,i=max（x）找出x的最大值y及位置icumsum（x）计算数组x的累加值std（x）数组x的元素的标准差cumprod（x）计算数组x的累加连乘值mean（x）计算x的平均值median（x）计算x的中位数sum（x）计算数值x的和5TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析统计分析函数prod（x）计算数组x的连乘值y=sort（x）排序函数举例1：

rains=126.8148.5173.0148.4194.7208.9;328.8300.7268.3210.5278.4321.5;avg_rain=mean（rains）avg_rain=mean（avg_rain）max_rain=max（rains）max_rain,x=max（rains）6TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析统计分析函数min_rain=min（rains）s_sort=sort（rains）举例2：

x=12345;prod（x）cumsum（x）cumprod（x）7TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析标准差和协方差标准差的概念计算函数标准差函数stdrains=126.8148.5173.0148.4194.7208.9328.8300.7268.3210.5278.4321.5s_dev=std（rains）协方差的概念函数cov（x）计算向量x的方差；对于矩阵，cov（x）对应于协方差矩阵cov（rains）sqrt（cov（rains）8TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析标准差和协方差当x和y是等长度的观测值时，可用cov（x,y）计算其协方差，运算等同于cov（x（:

）y（:

）举例：

x=126.8148.5173.0148.4194.7208.9328.8300.7268.3210.5278.4321.5;y=120.8138.5158.0148.4194.7208.9328.8300.7268.3210.5378.4221.5;cov（x,y）9TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析综合分析函数datastats（）函数ds=datastats（x）对数据x做统计运算，x应该是一个实数列向量，如果以复数代入，其虚数部分将被忽略，ds的内容如下：

ds.num样本容量ds.max最大值ds.min最小值ds.mean平均值ds.median中值ds.range最大值-最小值ds.std标准差10TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析综合分析函数datastats（）函数xds,yds=datastats（x,y）对数据x和y同时做统计运算，结果赋给变量xds和yds举例：

x1=15141110974342;x2=13171511131713111111;da1=datastats（x1）da1,da2=datastats（x1,x2）11TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析相关系数函数r=corrcoef（x,y）计算数组x和y的相关系数矩阵举例：

x1=15141110974342;x2=13171511131713111111;r=corrcoef（x1,x2）12TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析相关系数函数r,p=corrcoef（x,y）同时输出一个矩阵P，其元素值就是不相关加上检验的p值，意义为在真实相关为0的前提下，通过随机抽取获得观察值，出现相关系数r值的可能性举例：

x1=15141110974342;x2=13171511131713111111;r,p=corrcoef（x1,x2）13TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析相关系数函数r,p,rlo,rup=corrcoef（x,y）同时返回矩阵rlo和rup，它们都是与R等大小，并包含有对每一系数的95%置信区间的上下边界举例：

x1=15141110974342;x2=13171511131713111111;r,p,rlo,rup=corrcoef（x1,x2）意义：

对于相关系数0.5232的95%置信区间的上下边界在-0.15870.8672之间14TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析随机数分析工程问题时，如果要仿真不可预测且不规律的现象，可利用随机数方式产生近似值随机数根据统计特征分为：

均匀随机数和常态随机数均匀随机数是指其值均匀分布于给定区间常态随机数则是呈现高斯分布均匀随机数MATLAB中rand（）函数可生成区间01上的均匀随机数rand（n）或rand（m,n）15TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析均匀随机数rand（1,600）;subplot（2,1,1）;hist（ans,20）subplot（2,1,2）;plot（ans）16TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析均匀随机数如果要产生不介于01区间的随机数，采用将随机数值从01区间转换到其他区间如果区间为a,b，a为下限值，b为上限值，公式为：

x=（b-a）*r+a举例：

生成一个107的随机数矩阵，其值再130之间的整数fori=1:

10;x（i,1:

7）=round（29-0）*rand（1,7）+1）;end;fori=1:

10;sort（x（i,1:

7）;end17TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大学摩擦学研究所MATLAB统计分析统计分析常态随机数用MATALB的randn（）函数可生成符合常态分布的随机数，符合高斯分布函数randn（n）和randn（n,m）分别可产生nn和nm随机数矩阵，其平均值为0和方差为1举例：

x=-2.9:

0.2:

2.9;y=randn（1,5000）;subplot（2,1,1）;hist（y,x）title（Histogramofnormalrandomdata）y1=rand（1,5000）;subplot（2,1,2）;hist（y1）title（Histogramofuniformrandomdata）18TribologyResearchInstituteTribologyResearchInstituteSOUTHWESTJIAOTONGUNIVERSITYSOUTHWESTJIAOTONGUNIVERSITY西南交通大