地统计插值方法课件.ppt
- 【下载声明】
1. 本站全部试题类文档,若标题没写含答案,则无答案;标题注明含答案的文档,主观题也可能无答案。请谨慎下单,一旦售出,不予退换。
2. 本站全部PPT文档均不含视频和音频,PPT中出现的音频或视频标识(或文字)仅表示流程,实际无音频或视频文件。请谨慎下单,一旦售出,不予退换。
3. 本页资料《地统计插值方法课件.ppt》由用户(晟晟文业)主动上传,其收益全归该用户。163文库仅提供信息存储空间,仅对该用户上传内容的表现方式做保护处理,对上传内容本身不做任何修改或编辑。 若此文所含内容侵犯了您的版权或隐私,请立即通知163文库(点击联系客服),我们立即给予删除!
4. 请根据预览情况,自愿下载本文。本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
5. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007及以上版本和PDF阅读器,压缩文件请下载最新的WinRAR软件解压。
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 统计 方法 课件
- 资源描述:
-
1、天行健,君子以天行健,君子以自强不息自强不息;地势坤,君子以地势坤,君子以厚德载物厚德载物。摘自:周易第七章第七章 地统计插值地统计插值第七章第七章 地统计插值地统计插值一、克立格法简介理解二、简单克立格与普通克立格掌握三、泛克立格法掌握四、协同克立格法掌握五、其他克立格法理解+了解1、克里格法概念及分类克立格法(kringing),又称克里金法、空间局部估计或空间局部插值,是建立在变异函数理论及结构分析基础上,在有限区域内对区域化变量的取值进行无偏最优估计的一种方法一、克立格法简介分类:分类:n 简单克立格法n 普通克立格法普通克立格法(Ordinary Kriging)点克立格法(Puct
2、ual Kriging)块段克立格法(Block Kriging)n 泛克立格法(Universal Kriging)n 协同克立格法(Co-Kriging)n 对数正态克立格法(Logistic Normal Kriging)n 指示克立格法(Indicator Kriging)n 析取克立格法(Disjuctive Kriging)一、克立格法简介2、克里格估计量(Kriging estimator)显然,估计的好坏取决于权重系数显然,估计的好坏取决于权重系数i一、克立格法简介待估点或待估块段待估点或待估块段的估计值的估计值待估点或待估块段影响范围内的有待估点或待估块段影响范围内的有效样本
3、值效样本值权重系数权重系数克立格估计量克立格估计量(1 1)无偏估计)无偏估计(2 2)最优估计)最优估计第七章第七章 地统计插值地统计插值一、克立格法简介理解二、简单克立格与普通克立格掌握三、泛克立格法掌握四、协同克立格法掌握五、其他克立格法理解+了解1、简单克里格法、简单克里格法二、简单克立格与普通克立格设Z(x)为区域化变量,满足二阶平稳假设或本征假设,其数学期望为m,为已知常数,协方差函数为C(h),变异函数为(h)点:点:块段:块段:待估块段为V,中心点在x,平均值为在待估块段V附近有n个样点xi(i=1,2,n),其观测值为Z(xi)(i=1,2,n),则:二、简单克立格与普通克立
4、格EZ(x)=m已知,令Y(x)=Z(x)-m,则 EY(x)=EZ(x)-m=EZ(x)-m=0 CY(x),Y(y)=EY(x)Y(y)待估块段新变量,观测值新变量为:Y(xi)=Z(xi)-m目标目标:找出一组权重系数 i(i=1n),使得Zv#(x)成为ZV(x)的线性、无偏、最优估计量二、简单克立格与普通克立格Y(V)的估计值Yv#是Y(xi)(i=1,2,n)的线性组合,则:则估计Z(V)的问题转化为估计Y(V)的问题在满足以下两个条件时,Yv#是Y(V)的线性、无偏、最优估计量。(1)(1)无偏性无偏性二、简单克立格与普通克立格(2)(2)最优性最优性对系数 i求偏导,并令其为0
5、,令n=3则偏导数通式为:则n个未知量,n个方程,解出系数 i,即,整理得简单克立格方程组,二、简单克立格与普通克立格得Y(V)的简单克立格估计量,记为Yk#得到简单克立格估计方差:二、简单克立格与普通克立格2、普通克里格法、普通克里格法二、简单克立格与普通克立格设Z(x)为区域化变量,满足二阶平稳假设或本征假设,其数学期望为m,为未知未知常数,协方差函数为C(h),变异函数为(h)点:点:块段:块段:估计中心点为x,体积为V的待估块段的平均值ZV(x),在待估块段V附近有n个样点xi(i=1,2,n),其观测值为Z(xi)(i=1,2,n),则:二、简单克立格与普通克立格目标目标:找出一组权
6、重系数 i(i=1n),使得ZV#(x)成为ZV(x)的线性、无偏、最优估计量(1)(1)无偏性无偏性即普通克立格是即普通克立格是条件条件无偏无偏二、简单克立格与普通克立格(2)(2)最优性最优性补充:条件极值补充:条件极值 拉格朗日乘数法拉格朗日乘数法若求函数z=f(x,y)在条件(x,y)=0下的可能极值点,可先建立拉格朗日函数:拉格朗日函数:F(x,y)=f(x,y)+(x,y)。为拉格朗日乘子拉格朗日乘子,对函数求x,y变量的一阶偏导偏导,并使之为0,与条件条件联立建立方程组。fx(x,y)+x(x,y)=0fy(x,y)+y(x,y)=0(x,y)=0对自变量多于两个的函数z=f(x
7、1,x2xn)在条件i(x1,x2xn)=0下的可能极值点求法同上:先建立拉格朗日函数,然后建立方程组。确定函数:建立方程组:普通克立格方程组普通克立格方程组n+1个未知量(i(i=1n),),n+1个方程则普通克立格估计量为Zk#:则普通克立格估计方差:普通克立格方程组和估计方差的变异函数表达:在二阶平稳条件二阶平稳条件下,可采用协方差或变异函数的方程组或计算式进行求解计算;在本证假设条件本证假设条件下,则只可采用变异函数的表达式求解计算当样品不是点支撑,而是以xi为中心的小块段vi时,普通克立格方程组和估计方差的表达式为:你会发现变化的只是xi到vi普通克立格方程组和估计方差的矩阵矩阵表达
8、:展开后?则普通克立格方程组矩阵矩阵表达为:K称为普通克立普通克立格矩阵格矩阵估计方差的矩阵矩阵表达:普通克立格方程组和估计方差矩阵的变异函数表达:注意:当样品不是点支撑,而是以xi为中心的小块段vi时,以上矩阵适用,只需将xi换为vi,但含义不同。假设条件适用情况同前。作业作业4 4 (1)写出普通克里格方程组和估计方差的推导过程,并分别用协方差函数和变异函数表达;(2)写出普通克里格方程组和估计方差的矩阵表达式推导过程,并分别用协方差函数和变异函数表达;(1)点普通克里格法的计算规则格网采样数据3、普通克里格法的计算、普通克里格法的计算二、简单克立格与普通克立格不规则格网采样数据(2)块段
9、普通克里格法的计算点采样数据块段采样数据(1)点点普通克里格法的计算普通克里格法的计算规则格网采样数据规则格网采样数据二、简单克立格与普通克立格在一研究区内,Z(x)为区域化变量,满足二阶平稳假设,协方差函数为C(h),变异函数为(h),且是一个二维各向同性各向同性的球状球状模型普通克立格估计量为:=22-0.51812.66+0.022 4.98+0.089 1.72+0.371 9.84+0.916=12.44(1)利用excel软件,用协方差函数表达式将前述的普通克里格插值计算过程做出,并计算出相应的估计方差;(2)利用excel软件,用变异函数表达式将前述的普通克里格插值计算过程做出,
10、并计算出相应的估计方差。备注:excel矩阵求逆函数minverse(),转置函数transpose(),相乘函数mmult()核心内容,必须掌握核心内容,必须掌握作业作业5 50.5180.0220.0890.3710.1470.1710.6050.077“屏蔽”作用(1)点点普通克里格法的计算普通克里格法的计算 不规则格网采样数据不规则格网采样数据二、简单克立格与普通克立格在一研究区内,Z(x)为区域化变量,满足二阶平稳假设,协方差函数为C(h),变异函数为(h),且是一个二维各向同性各向同性的指数指数模型。模型参数为C0=0,C=10,a=10,模型为x1x2x3x4x5x6x7Vx0(
11、2)块段块段普通克里格法的计算普通克里格法的计算二、简单克立格与普通克立格点采样数据采样点为点支撑,待估为块段v:可将v离散成若干点,求采样点和离散点之间的协方差函数或变异函数值之和,后取平均01020304V01020304v101020304V(2)块段块段普通克里格法的计算普通克里格法的计算块段采样数据业精于业精于勤勤,荒于嬉;,荒于嬉;行成于行成于思思,毁于随。,毁于随。摘自:韩愈的进学解第七章第七章 地统计插值地统计插值一、克立格法简介理解二、简单克立格与普通克立格掌握三、泛克立格法掌握四、协同克立格法掌握五、其他克立格法理解+了解第四节第四节 泛克里格法泛克里格法泛克里格法产生的原
展开阅读全文