第十二讲随机变量的数字特征.pptVIP

中值定理与导数的应用 */22 Matlab实验（十二） 一、数据特征的Matlab实现 1.在MATLAB中趋势统计量的计算命令如下: mad (x) 平均绝对偏差 harmmean (x) 调和平均 iqr (x) 四分位极差 geomean (x) 几何平均 std (x) 标准差 trimmean (x,n) 切尾平均 var (x) 方差 median (x) 中位数 range (x) 极差 Mean(x) 算术平均 MATLAB命令 离中趋势 MATLAB命令 集中趋势 其中x是原始数据, n是上、下切除的百分数 , 例1. 已知以下100个数据： 459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680 926 653 164 487 734 608 428 1153 593 844 527 552 513 781 474 388 824 538 862 659 775 859 755 649 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638 699 634 555 570 84 416 606 1062 484 120 447 654 564 339 280 246 687 539 790 581 621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851 计算该组数据特征. 解：a=[459 362 624 542 509 584 433 748 815 505 612 452 434 982 640 742 565 706 593 680 926 653 164 487 734 608 428 1153 593 844 527 552 513 781 474 388 824 538 862 659 775 859 755 649 697 515 628 954 771 609 402 960 885 610 292 837 473 677 358 638 699 634 555 570 84 416 606 1062 484 120 447 654 564 339 280 246 687 539 790 581 621 724 531 512 577 496 468 499 544 645 764 558 378 765 666 763 217 715 310 851]; a1=a;b=a1(:); %将矩阵变成数列 T=[mean(b),median(b), trimmean(b,10), geomean(b), harmmean(b),range(b),var(b),std(b),iqr(b),mad(b)] 150.86 平均绝对偏差 499.06 调和平均 243.5 四分位极差 559.68 几何平均 196.629 标准差 600.64 切尾平均 38663.03 方差 599.5 中位数 1069 极差 600 算术平均 计算结果 离中趋势 计算结果 集中趋势 计算结果如下： 一般情况下中位数、切尾平均与四分位极差对异常值是稳健的。其他数据特征容易受异常值影响 2. 异常值的判别 在探索性数据分析时，有一种判别异常值的简单方法，首先计算数据的下、上截断点，数据中小于下截断点的数据为特小值，大于上截断点的数据为特大值，二者都是异常值. 1).数据的下、上截断点 计算下、上截断点的公式如下： 其中，R为四分位极差， 分别称为下四分位数与上四分位数 . 计算样本的P分位数： 显然 例2. 判别例1中的数据有无异常值. 对于0≤1/p1，样本容量为n的样本 其次序统计量(按照从小到大排序)记为： 则样本的P分位数 解：记原数据从小到大的次序统计量为x,命令为 x=sort(a(:)) 因为下(上)四分位数的np=25(或75)为整数，故有命令： Q1=(x(25)+x(26))/2 Q3=(x(75)+x(76))/2 R = iqr (x)