www.gusucode.com > EasyKrig_V3.0工具箱matlab源码程序 > EasyKrig_V3.0/variogram/get3dvariopara_theo.m
function [model,nugt,sill,L,p,b]=get3dvariopara_theo() %% get parameter from the panel edit field (set by user) of semi-variogram/correlogram %% to compute model-based variogram (theory) % %% Last Modified December 22, 2001 %% %% Kriging Software Package version 3.0, May 1, 2004 %% Copyright (c) 1999, 2001, 2004, property of Dezhang Chu and Woods Hole Oceanographic %% Institution. All Rights Reserved. global hdl para %%% read in parameters model=get(hdl.vario.model,'Value'); value=get(hdl.vario.correlogram,'value'); if value == 1 model=-model; % correlogram model end nugt=str2num(get(hdl.vario.nugt_edit,'String')); sill=str2num(get(hdl.vario.sill_edit,'String')); L=str2num(get(hdl.vario.lscl_edit,'String')); p=str2num(get(hdl.vario.powr_edit,'String')); b=str2num(get(hdl.vario.hole_edit,'String')); if nugt > para.vario.max_nugt para.vario.max_nugt=1.2*nugt;end if sill > para.vario.max_sill para.vario.max_sill=1.2*sill;end if L > para.vario.max_lscl para.vario.max_lscl=1.2*L;end if p > para.vario.max_powr para.vario.max_powr=1.2*p;end if b > para.vario.max_hole para.vario.max_hole=1.2*b;end para.vario.nugt=nugt; para.vario.sill=sill; para.vario.lscl=L; para.vario.powr=p; para.vario.hole=b; para.vario.model=abs(model); para.vario.range=str2num(get(hdl.vario.range_edit,'String')); para.vario.res=str2num(get(hdl.vario.res,'String')); para.vario.vario=get(hdl.vario.variogram,'value'); para.vario.corr=get(hdl.vario.correlogram,'value');