=2之2i1C4(4,a)-2,-2C(W,7a=0 0入ak k/1 Va&=1,2,…,nik=1,2,…,K -2(空 a。-1)=0 眼.-2（克a)=0 Vk卡k。 经整理得到协同克立格方程组为(20)式所示： ∑Bx1C(wi,a)-4=Ca(。va) Va=1,2,…k5k=1,2,…,K ho at。=1 (20) d南0●1 a=0 (20)式中的 (空：+K)个未知数(2 n个权系数元a及K个拉格朗日乘子4：) 可以由（公”，+K)个线性方程组求解，即用(20）式求出的，才能保证Z.是 Z口。的最优无偏线性估计量。 相应的协同克立格方差为， .a C (m va (21) 在普通克立格法中，变异函数y()与协方差C()之间存在着如(9)式所示的关系，因此可 以很方便的用y(h)代替C(h)表示普通克立格方程组，但在协同克立格法中，y'(h)与 C(h)之间的关系如(10)式所示，从而不能简单地用y4'.(h)来表示协同克立格方程组 (式(20))，只有当 C(h)=C(h) (22) 时，才能使 ys'a(h)=Cs'(0)-C'(h) (23) 成立，这时，协同克立格方程组既可用互协方差函数也可用互变异函数表示，即用(-y (h))代替方程组(20)式中的(C'()导出用互变异函数表示的协同克立格方程组。 要使协同克立格方程组（或(20）)具有唯一解的条件是： (1)互协方差矩阵〔Cs(vw',vα)〕严格正定，为此，必须采用正定的点协同区域 100“ 云 走 , 一 i n ` I 习 人刀 。 , C * , , ( 。 , , , , 。 a , ) 一 2 ; * 一 Z C 。 。 , ( v 。 。 , v a , ) = 口 。 , . 1 V a , = 1 , 2 , … , n , ; 吞= 1 , 2 , … , K n 为。 一 2 ( 习 “ a , 。 一 i 口石 。 一 1 ) 一 刀 舀 一 2 ( 名 “ a 。 ) = 。 V k今掩 。 汀 一 一。 即一、解叭 经整理得到协 同克立格方程组为 (2 0) 式所示 : 兀 打 ` I 习 习 久, 。 ` c , , ` ( v , 。 ` , t, a 。 ) 一 “ * 二 C ` 。 ` ( v , 。 U a . k , 一 l 口 , , 一 i 刀& o 习 久a 肠 “ 1 口 为。 一 1 V a , = 1 , 2 , … k * ; k = 1 , 2 , … , K ( 2 0 ) 久 一 名心 少、厂| | (20 , 式 中的 l(P 一 ` ) 个 ” 数 (恩 二 个权系数 “ 。 及 ` 个拉格 朗 日乘子 动 可以 由 (习 , 。 十 K ) 个线性方程 组 求 解 , 即 用 ( 20 ) 式 求 出 的 肠 , 才 能 保 证 Z认 。 是 曲. 1 Z 甲 , 。 的最优 无偏线性估计量 。 相应的协 同克立格方差为 , a 2 7 吞。 尤 月 舀 = c , 。 * 。 ( v 。 。 , v , 。 ) + ; , 。 一 习 习 “ 。 , c , 。 , ( v ` 。 , v a ; ) 孟七 1 口盖 . 1 ( 2 1 ) 在普通克立格法 中 , 变异函数 夕 (h) 与协方差 C ( h) 之 间存在着如 ( 9) 式所示的关系 , 因 此 可 以很方便的用 夕h( ) 代替 C ( h) 表示普通克立格方程组 , 但在协同 克立格法 中 , 夕, , . (h) 与 C , , , h( ) 之间的关系如 (l 0) 式所示 , 从而 不能简单地 用 夕; , * ( h) 来表示协同 克 立 格 方 程 组 (式 ( 2 0 ) ) , 只 有当 C * , 。 (h ) = C 。 。 , ( h ) ( 2 2 ) 时 , 才能使 夕 , , 。 (h ) = C 、 , , ( o ) 一 C 。 , 。 ( h ) ( 2 3 ) 成 立 , 这时 , 协同克立格方程 组既可 用互协方差函数也可用 互变异函 数表示 , 即用 ( 一 夕。 ` . ( h) ) 代替方程组 ( 2 0) 式 中的 (c , , , ( h ) 导 出用互变异函数表示的协同克立格 方程组 。 要使 协 同克立格方程组 (或 (2 0 ) 具有 唯一解 的条件是 : (1 ) 互协方差矩 阵 〔 C , , , (巧 。 , , 。 a * ) 〕严 格正定 , 为此 , 必须采用 正定的点 协 同 区 域 1 0 0