クリギング
クリギングは、ランダム又は確率過程共分散関数に基づいて空間回帰モデルおよび補間アルゴリズムです。
クリギング補間式は次のとおりです。
各点の重金属汚染の値は右重み値の各点について重金属汚染ポイント値として知られている式に挿入されます。
重量BLUP理論とクリギング重みを解きます：
ランダム誤差を含む線形システムとして表される確率変数の推定値は、選択は、推定値と真値の分散を最小にする線形BLUPシステムパラメータとして表すことができます。
前記一般クリギング重みと呼ばれる重み係数のための不明点、{ランダムサンプル、、。推定値と同一の真の値は、両方の最小分散の数学的期待値によって定義された既知の分散
上記重み係数に溶液のBLUP重量は、サンプル点と未知点間の共分散関数を含みます。
クリギングは統計グリッド法の多くの分野において非常に有用であり、それは汚染物質の濃度分布の問題と一致し、そしてクリギングの結果は、より自然な産生、効果的アウトライアを回避することができますまた、標準誤差を与えることができ、感謝クリギングアルゴリズムはまた、既知のポイントとの間の関係を検討し、考慮に推定位置と、既知の地点の地点間の相互関係をとります。したがって、客観的応答の汚染物質のより良好な分布が、それは推定値の比較的高い精度です。
コードは以下の通りです
サンプル％地形マッピング
CLC、クリア
データ=関数xlsread（ 'zz.xls'、 '付録1'、 'B4：D322' ）;％は、 ファイルの読み取り
％をデータ取得
のS =データ（：、1：2である。）;
Y =データ（：、3）;
％使用クリギング
シータ= [10 10];ロブ = 1E-1 1E-1]; UPB = [20 20];％ パラメータ
％のクリギング補間アルゴリズムツールボックス呼び出し
％フィッティング操作
[dmodel、パフォーマンス・] = dacefit（S、Y、regpoly0 @、@corrspherical、シータ、LOB、UPB）
M = 100;
％補間演算
X = gridsamp（[0 0; 30000 20000]、M） ;
[YX MSE] =予測（X-、dmodel）;
％補間得る
X1を= RESHAPE（X - （：、1）、M、M）; X2 = RESHAPE（X-（：、2）、M、M）;
YX = RESHAPE（YX、大きさ（X1））;
％プロット
図;
サーフ（X1、X2、YX）は
、ON HOLD
ホールドオフ
は、xlabel（ 'X / M'）
ylabelの（ 'Y / M'）
zlabel（ '高度'）
タイトル（ 'サンプリング地形'）
図;
contourf（X1、X2、YX）％の計画を行う
[C、H] =輪郭（ X1、X2、YX）を、
clabelが（C、H）
は、xlabel（ 'X / M'）
ylabelの（ 'Y / M'）
zlabel（ '高度'）
タイトル（ 'サンプリング地形'）
％の汚染物質濃度分布
CLC、クリア
'としての' B = { 、 'CDの'、 'のNi'、 '銅'、 '水銀'、 '鉛'、 '亜鉛'} 'のCr';
のNd =関数xlsread（ 'zz.xls'、 '付録2'、 'B4：I322' ）;ファイルを％読み取る
S =関数xlsread（ 'zz.xls'、 ' 付録1'、 'B4：C322を' ）;％は、 ファイル読み取る
データ読み取りに％のループ
= Iために1 :. 8。
Y = FIX（のNd（： 、I））;
％使用クリギング
シータ= [10 10]。LOB = 1E-1 1E- 1]; UPB = [20 20];％のパラメーター
％コールクリギングアルゴリズムツールボックス
％嵌合操作
[dmodel、パフォーマンス・] = dacefit（S、Y、regpoly0 @、@corrspherical、シータ、LOB、UPB）
M = 100;
％計算補間
X-gridsamp =（[0 0; 30000 20000]、M）;
[YX MSE] =予測（X-、dmodel）;
％補間う
X1 = RESHAPE（X - （：、1）、M、M）; X2 = RESHAPE（X-（：、2）、M、M）;
YX = RESHAPE（YX、サイズ（X1））;
％プロット
図、
メッシュ（X1、X2、YX）は
、ON HOLD
ホールドオフ
は、xlabel（ 'X / M'）
ylabelの（ 'Y / M'）
zlabel（ '濃度'）
タイトル（[Bを（ I）]）
図;
contourf（X1、X2、YX）
％予定ですか
[C、Hを] =輪郭（X1、X2、YX）
は、xlabel（ 'X / M'）
ylabelの（ 'Y / M'
zlabelを（ "濃度「）
タイトル（[B（I）]）
端
あなたは、ツールキットのファイルを呼び出す必要があります
そこに私のBaiduの雲のディスクでは、私はあなたを助けるために願っています
リンクします。https：//pan.baidu.com/s/1O-mqKowNBJ06llEldHu90A
抽出コード：g4wl