Matlab第二章-实验数据的处理-2ppt课件

.1/47,2.2实验数据错误值的剔除,错误值（过失误差）的处理,测量数据在其采集、传输与记录过程中，有可能造成数据丢失或个别数据产生不切合实际的偏差，这种数据称为错误值，或称为奇异项或坏值。,人为因素（读错、记录错、操作错）不符合实验条件的环境突变（突然振动、电磁干扰等）,.2/47,2.2实验数据错误值的剔除,剔除错误值的目的：主要是恢复数据的客观真实性，而不是为了提高精度。,剔除错误值基本思想：规定一个置信水平，确定一个置信限度，凡是超过这个限度的误差，就认为它不属于随机误差的范围，是错误值，从而予以剔除。由于置信限度与统计量的分布及测量次数有关，剔除了一次坏值后，还需对留下的数据重新检验，直到满足置信限度。,.3/47,2.2实验数据错误值的剔除,剔除错误值的常用方法拉依达方法：非等置信概率肖维勒方法：等置信概率一阶差分法：预估比较法,.4/47,2.2实验数据错误值的剔除,1拉依达方法,如果某测量值与平均值之差大于标准偏差的三倍，则予以剔除。,这种方法用于正态分布的数据，当n时，其置信概率为99.7%。实际中，即使是在正态分布的条件下，由于测量次数不同，其置信概率就不同，可见拉依达方法是一种非等置信概率方法。,.5/47,2.2实验数据错误值的剔除,2肖维勒方法,对于遵从正态分布的情况，在n次测量结果中，如果某误差可能出现的次数小于半次时，就予以剔除。这实质上是规定了置信概率为1-1/2n，根据这一置信概率，可计算出肖维勒系数，也可从表中查出，当要求不很严格时，还可按下列近似公式计算。,.6/47,2.2实验数据错误值的剔除,肖维勒系数表,.7/47,2.2实验数据错误值的剔除,如果某测量值与平均值之差的绝对值大于标准偏差与肖维勒系数之积，则该测量值xi被剔除。,肖维勒方法是一种等置信概率方法。,.8/47,2.2实验数据错误值的剔除,3一阶差分法,一阶差分法是一种预估比较法，是用前两个测量值来外推即预估新的测量值，然后用预估值与实际测量值比较，并事先给定其允差限值，称作误差窗，以此来决定该测量值的取舍。,预估值：,比较判别：,.9/47,2.2实验数据错误值的剔除,3一阶差分法,一阶差分法特点：,适合于实时数据采集与处理过程精度除了与误差窗的大小有关外，还与前两点测量值的精确度有关被测物理量的变化规律不是单调递增或单调递减函数，这一方法将在函数的拐点处产生较大的误差，严重时将无法使用,.10/47,2.2实验数据错误值的剔除,例2.2.1利用肖维勒方法对下列数据的错误值进行剔除,1.50341.50621.50341.50241.49852.50001.50071.50671.49931.4969,上列数据保存于文件erro.dat,.11/47,2.2实验数据错误值的剔除,std.m,erroeliminate.m,.12/47,2.3实验数据的平滑滤波,在物理实验中，不论是人工观测的数据还是由数据采集系统获取的数据，都不可避免叠加上噪声信号。为了提高数据的质量，必须去除噪声。,.13/47,2.3实验数据的平滑滤波,错误值剔除：,.14/47,2.3实验数据的平滑滤波,数字滤波,时域滤波,频域滤波,平滑滤波,对一个按某种等间隔规律（如等时测量，等电压间隔测量等）测得的数据序列，对某一位置处的数据，利用其前后的数据经某种运算后，得到一个新的平滑数据，从而得到一个新的平滑后的数据序列。,平滑滤波是依据随机噪声的概率统计性质，通过适当运算，使噪声抑制到满意的程度。,.15/47,2.3实验数据的平滑滤波,2.3.1单纯移动平均,单纯移动平均法是在时域序列中，某一数据点前后对称地取出2n+1个数据（平滑宽度或平滑窗口），求其平均值作为结果序列中的数据，而取代原始数据。,.16/47,2.3实验数据的平滑滤波,例2.3.1分别利用3、9、15点单纯移动平均对下列数据进行平滑。,.17/47,2.3实验数据的平滑滤波,.18/47,2.3实验数据的平滑滤波,.19/47,2.3实验数据的平滑滤波,优点：方法简单，计算方便。缺点：方法误差：信号失真。数据丢失：前后各n个数据无法平滑。适用性：适用于变化缓慢的数据。,平滑窗口越大，平滑效果越好，但失真也越大。,.20/47,2.3实验数据的平滑滤波,2.3.2加权移动平均,加权的基本思想：平均区间内中心处数据的权值最大，愈偏离中心处的数据权值愈小。这样就减小了对真实信号本身的平滑作用。,权重系数可以采用最小二乘原理，即使平滑后的数据以最小均方差逼近原始数据。,.21/47,2.3实验数据的平滑滤波,设自变量x以步长h作等距观测的数据y,作变换,.22/47,2.3实验数据的平滑滤波,数据变为,构造平滑公式：,平滑多项式的系数由最小二乘原理确定,.23/47,2.3实验数据的平滑滤波,即,t取最靠近i的2n+1个整数值（平滑点的数目），且应保证m2n+1N，t=-n,-n+1,0n-1,n,.24/47,2.3实验数据的平滑滤波,线性平滑：m=1平滑公式为,系数由下式确定：,.25/47,2.3实验数据的平滑滤波,三点线性平滑n=1，t=-1,0,1,解得系数：,.26/47,2.3实验数据的平滑滤波,平滑公式为：,简记为：,.27/47,2.3实验数据的平滑滤波,三点线性平滑权重系数表,.28/47,2.3实验数据的平滑滤波,二次平滑：m=2，平滑公式为,系数由下式确定：,.29/47,2.3实验数据的平滑滤波,五点二次平滑n=2，t=-2,-1,0,1,2,.30/47,2.3实验数据的平滑滤波,五点二次平滑权重系数表,.31/47,2.3实验数据的平滑滤波,2.3.3matlab平滑函数,格式：Z=smooth(Y,span,method),说明：Z：平滑后的数据向量Y：被平滑的数据向量span：平滑点数，缺省为5点method：平滑方法，缺省为移动平滑moving-Movingaverage(default)单纯移动平均lowess-Lowess(linearfit)线性加权平滑loess-Loess(quadraticfit)二次加权平滑sgolay-Savitzky-Golayrlowess-RobustLowess(linearfit)rloess-RobustLoess(quadraticfit),.32/47,2.3实验数据的平滑滤波,.33/47,2.3实验数据的平滑滤波,.34/47,2.3实验数据的平滑滤波,.35/47,上机1,1利用一阶差分法对下列数据的错误值进行剔除,1.50341.50621.50341.50241.4