0 核数据处理预备知识（续） 核数据处理课件.ppt

1)设： 2)查正态概率积分表，找出给定显著水平α时的值kα 3) 用kα与实测计算的值k相比， 若kkα，认为差异不显著，数据可靠； 若kkα，认为差异显著，数据不可靠。 例2：分别测量10分钟得两个计数率1010cpm和1069cpm，问计数设备工作是否正常？ (取显著水平α=0.05，查表得：kα = 1.96) 解： Δ(k)为标准偏差，k个数据中，偏差Δ大于标准偏差Δ(k)的数据个数不得超过半个(即1/2k)。即在k个数据中，某一个数据与平均值的偏差Δ出现的概率小于1/2k时应舍弃。标准偏差Δ(k)由高斯分布来求，与k有关。 * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * 东华理工大学 张怀强 数据的预处理——网格变换 距离导数法（近点按距离加权平均） 计算机编程 步骤1：先统计出邻域中数据点的个数n，并将这些数据取出，放在另一新数组中。 步骤2：新数组做循环累加运算。 不足: 距离法是取离网格点最近的几个数据点值来计算网格点值, 不考虑方向, 所取的点很可能集中到某一侧, 其它方向会取不到点。 东华理工大学 张怀强 数据的预处理——网格变换 void distnet(float x[], float y[], float z[], float x1[], float y1[], float z1[], int nt, int nt1, float d) { float data[10], dist[10], dtemp; for (int i=0; int1; i++) { int n=0; for (int j=0; jnt; j++) { dtemp=srqt((x1[i]-x[j])*(x1[i]-x[j])+(y1[i]-y[j])*(y1[i]-y[j])); if (dtemp==0) { z1[i]=z[j]; n=-1; break; } else if (dtempd) { data[n]=z[j]; dist[n]=dtemp; n=n+1; } } if (n==0) z1[i]=0; else { float temp1=0, temp2=0; for (int k=0; kn; k++) { temp1=temp1+data[k]/dist[k]; temp2=temp2+1/dist[k]; } z1[i]=temp1/temp2; } } } 东华理工大学 张怀强 数据的预处理——网格变换 方位法(按方位取点加权平均) 按方位取点加权法: 以网格点(x,y)为中心，把区域划分成若干个象限，从每个象限内取点作加权平均, 这就克服了距离法偏向的缺点。 具体方法: 若求某个网格点(x, y) 的函数值时, 则以(x, y)为原点将平面分成四个基本象限； 再把每个象限等分成n0 份, 这样就把全平面分成4n0 等份。 东华理工大学 张怀强 数据的预处理——网格变换 方位法(按方位取点加权平均) 具体方法: 然后在每个等分角内寻找一个离(x, y)最近的数据点, 其值为Z i1, 它到(x, y)的距离为di1； 则网格(x, y)上的值为: 东华理工大学 张怀强 数据的预处理——网格变换 趋势面拟合法 趋势面分析是地学领域常用的数学方法, 它可以把长周期的趋势性变化和短周期的局部性变化分开。 根据趋势性变化可以总结规律, 根据局部性变化可以发现异常。 东华理工大学 张怀强 数据的预处理——网格变换 趋势面拟合法 从计算机制图角度看, 趋势面分