式(16)可改写为： E(Z。-io}=E{2o}-》a.E2a}-∑∑1aEZ} a4■1 =m〔i-芝…门a.m豆， (17) aho前1 a是■1 要使式(17)等于0必须有如下K个条件成立： ah。=1 (18) 乃 1as=0 atel 程质0 式(18)是无偏性的充分条件组。当”。=0时条件∑a。=1不能成立，也就是说，在 0k0●1 应用协同克立格法估计时，主变量，至少应有一个观测值Z:。 在上述条件下的估计方差是： 品.=E27.-Z} =0)-222e,） ha2p C(mn +宫良点点 (19) 式(19)中C(?。,V)为待估域V。中变量乙，(x)的直接协方差函数的平均值；C (7o,va4)为待估域7。与信息域a4之间区域化变量Z:。(x)与Za(x)的互协方差函数的 平均值；C(vaa,vs)表示信息域VaA与vgs'之间区域化变量Za4(x)与ZBs'(x)的互 协方差函数的平均值。 在式(18)表示的K个无偏条件约束下，要求式(19)所表示的估计方差0品：。达到极小时 来求式(15)中的诸权系数1·，可根据求条件极值的拉格朗日乘数法。 令F=-2(克a-1）-2,（分a） (k牛k。) 将F对诸2a,o,4(k卡k。)求偏导并令其为0， 99式 ( 16 )可改 写为 : 二、 z二。 。 一 忘 ` 。 } = 二、 z 二 , 。 } 一 誉 ; 。 , 。 E 、z 。 * 。 一 习 习 又 a o E { Z · , } a 几。 一 1 汤专 轰。 汀 玉一 1 月 告。 · m 告。 〔 , 一 习 , a , 。 〕 a 几 。 一 1 要使式 ( 1了) 等于 。 必须有如下 尤 个条件成立 : 习 久 a ` 。 = ; 习 , a 。 = 0 n 乙 一 习 二 , 。 习 久 a , ` 专 白o a 几 一 1 目 ` o ( 1 7 ) ( 1 8 ) 式 ( 18 ) 是无 偏性的 充分条件组 。 当 。 , 。 = 。 时条件 习 而 。 。 = 1 不能成立 , 也 就是说 , 在 J 告 O 一 1 应用协同克立格法 估计时 , 主 变量 寿 。 至少应有一个观测值 Z , 。 。 在上 述条件下 的估计方差是 : 。导 。 。 二 二 { 2 7 。 。 一 z 导 。 。 } K 月 为 = c 。 。 , 。 ( v 。 。 , v * 。 ) 一 2 习 习 几a 。 c , 。 , ( v 。 。 , 。 a * ) 生 J 口 几一 1 尤 f 十 习 习 九 ` ” ` I 习 习 几a , 几, , , C 。 * , ( 。 a 。 , ” , 。 , ) a 。一 1 刀 , I 一 1 ( 1 9 ) ` 一 1 k l . z 式 ( r 。 ) 中乞 , 。 。 。 (、 。 。 , v , 。 ) 为待估域 v ; 。 中变量 Z , 。 ( x ) 的直接协方差函数的平均值; c , 。 。 ( v , 。 , z, a 。 )为待估域 V . 。 与信息域 。 。 * 之’lia 区域化变量 Z , 。 ( 二 ) 与 Z a 。 ( x ) 的 互协 方差函数的 平 均值 ; C , , , ( Z, 。 。 , z, , , , ) 表示信息域 v 。 * 与 。 , , , 之间区域化变量 Z 。 。 (二 ) 与 Z 刀 ` , (二 ) 的互 协方差函数的平均值 。 在式 ( 18 ) 表示的 K 个无偏条件约束下 , 要求式 (1 9 ) 所表示的估计方差 『 令 , 。 达到极小时 来求式 ( 15 ) 中的诸权 系数 只 。 , , 可根据求条件极值的拉格 朗 日乘数法 。 令 二 一 a 导 、 。 一 2 ; 。 。 ( 习 “ a , 。 一 i ) 一 2 ; , ( n ` 习 , a * ) ` k今 k 。 , 口 玉。一 1 将 F 对诸几a , “ . 。 , 拌 。 ( k粉k 。 ) 求偏导并令其为