MATLAB克里格工具箱中文翻译版

MATLAB克里格工具箱(4.0版：2001年7月)翻译：阿童木看星星伊夫格拉顿等加拉弗勒克里格工具箱是分布式的自由和技术支持。规格克里格工具箱4.0版MATLAB6.1兼容。这是一个2.0版本的升级，已matlab下实现4.2编译，编译matlab下实现5.1和3.0版本。请注意，此升级只使用2-D矩阵，即使新的MATLAB版本支持更大的矩阵维数。优化功能要求的Matlab优化工具箱。然而，站在替代自我功能还提供了人，谁没有购买优化工具箱(参见“fitvario.m”志工具箱正常的研究是必要的。它是提供与克里格工具箱。说明此工具箱的发展是基于使用2个或3个标量的客观分析的必要性在物理海洋学的尺寸。这种类型的插值通常比更好的结果插值方法。此外，它的不可忽略的优势，使插值误差的估计。这个工具箱的功能几乎完全是从书Deutsch和Journel(1992)和Marcotte(1991)。变差函数的功能是墨西哥文件编制前,而协同克里格法的功能发表后，在Matlab格式，在1991年Marcotte的论文。所有的参数和例子可以发现，在英国，在该两本刊物。Journel和Huijbregts(1992)的书是最好的书semivariograms。一个完整的例子在物理海洋学的最优估计可以登曼和弗里兰(1985)发现的文件。同时，kridemo显示2-D目标的轮廓分析。登曼，吉隆坡和HJ斐然，1985年。相关秤，客观的测绘和统计检验Geostrophy超过大陆架。研究月RES,43：517-539。德语，C.V和AGJournel，1992年GSLIB:地统计的软件库和用户旨南。牛津牛津大学出版社，340页。Journel，AG和Huijbregts终审法院首席法官，1992年，矿业统计学。学术出版社，纽约，600页。Marcotte，D.1991。Cokrigeage与MATLAB。计算机及地质。17(9)：1265年至1280年。，建议或问题？许多功能还没有完全测试。请报告任何错误或问题伊夫格拉顿，INRS，淡，电话：(418)654-3764传真：(418)654-2600HYPERLINK"mailto:yves_gratton@inrs-eau.uquebec.ca"yves_gratton@inrs-eau.uquebec.ca克里格工具箱克里格工具箱内容变异函数的功能Confint置信区间。Fitvario优化的variogr。fitvari2variogr2优化（没有优化工具箱）。Fun估计变异函数从fitvario调用。gam2从vario2dr称为MEX文件。gam3从vario3dr称为MEX文件。gamv2从vario2di称为MEX文件。gamv2uv从var2diuv称为MEX文件。gamv3从vario3di称为MEX文件。Mrqmin最小二乘拟合的Levenberg-Marquardt方法。Outvario输出变差功能。Variogr模型的变异函数和相关图。variogr2variogr拟合程序不使用优化工具箱。vario2di变异函数的间隔不规则的2-D数据。vario2dr变异函数的定期间隔的2-D数据。vario3di变异函数的间隔不规则的3-D数据。vario3dr变异函数的定期间距为3-D数据。var2diuv变异函数的间隔不规则的2-D向量。克里格功能Barnes克里金数据的2-D空间滤波。CokriD尺寸的点或块协同克里格法。cokri2cokri引发的协同克里格法的功能。Davis点克里格使用戴维斯的方程组。FilrespBarnes的滤波器响应在波长域。TintoreBarnes过滤器在Tintore参数中的应用。（1991）相关功能Covsrt变换协方差矩阵的mrqmin。Deplie矢量矩阵转换。Gaussj高斯约旦消除线性方程的解。Kregrid矩阵（x2）的2-D网格坐标。kregrid3矩阵（x3）的3-D网格坐标。Kridemo克里格工具箱演示。Ksone从kstest称为MEX文件。Kstest柯尔莫哥洛夫-斯米尔诺夫正态性检验。mat3dp获取一个伪的3-D矩阵的值。mat4dp获取一个伪2,4-D矩阵的值。Means平均cokri2从所谓的功能。Mrqcof从mrqmin称为M文件。Transcokri2所谓的翻译功能。变异函数选项可用的变异函数是：4.5.6,7.6.半变异函数,交叉变异函数:统计:相关图:一般相对变异函数:两两和对变异函数:半变异函数对数:Semirodogram：Semimadogram：y{h}=52(-T-一vf)2M2M7OUj呻、1N闵氐qYgrW~{/-V(A)呻〕'■-?-1l-l;3=棗膚冷1M向YlW二TTryvV[ln(.Yf)-\n(ys)]--TH")r=i]MA)1机的曲)二硕叫-必|i=l2ScmivariogrammetHW:旳Q-I—斗乞工mn™■0变量的描述：C（H）=协方差距离h的函数，C（0）给出的方差。Y（H）=半方差=%[C（0）-C（H）]。h=分离载体。N（H）=样本对数Xi，yi=向量h隔开的样本对价值：xi是在开始的价值（或尾部）；yi是在间隔H尾部（或头部）的价值.zi,zi'etyi,yi'=和交叉变异函数（XI，YI）—样:yi和zi开始的价值和yi'和zi'是在间隔H.尾部的值有关详细信息，请咨询Deutsch和Journel（1992）。克里格选项可用克里格的选项是：1。简单的协同克里格法2。一个nonbias条件与普通协同克里格法（Isaaks和Srivastava）3。与Pnonbias条件的普通协同克里格法4。1阶漂移普遍协同克里格法5。2阶漂移普遍协同克里格法99。协同克里格法不执行，只计算样本方差SV协同克里格法是指用多个变量的克里格。当协同克里格法的被称为只在一个时间变量，结果将是那些简单克里格，普通克里格法，普遍克里格，点克里格或块克里格。Marcotte（1991年）的文件中可以找到更多的细节。ChiToolboxContents智工具箱内容使用这个工具箱是由正态性检验。所包含的函数的计算卡方概率函数和概率函数百分点。m文件可从MATLAB公共网站下载：HYPERLINK"http://www.mathworks.com/ftp/statv4/shtml"http://www.mathworks.com/ftp/statv4/shtml概率函数12Chiprob观察一个给定的〔2的概率值。Chitablex对于一个给定的概率值。Chiaux函数调用chitable。工具箱作者：彼得R.肖伍兹霍尔海洋研究所，伍兹霍尔，马02543（508）457-2000分机。2473HYPERLINK"mailto:pshaw@whoi.edu"pshaw@whoi.edubarnes巴恩斯目的克里金数据的2-D空间滤波。简介F=barnes（xi,yi,zi,c,g）说明巴恩斯（barnes）的过滤器是一个低通的2-D滤波器，其数学描述：工匕伉力耳力EF；（•）（】；◎）-H/J））其中爲（叮）二•D&]）=顽厂£儿（阳）£元心）n=1ff=l=勰P[—札心尸/4c]=exp[-Zrn（fJ）2/4cg].输入变量：xi：列网格坐标yi：行网格坐标ci：网格数据c,g：滤波器参数输出变量：F：过滤后的数据示例tintore功能提供了一个很好的例子，在物理海洋学的空间滤波器。参考文献马多克斯，类风湿性关节炎，1980年。一个分离宏观和中尺度的目的技术气象数据的特点。每月天气牧师，108：O八年至1121年。Tintore等，1991年。中尺度动力学和阿尔沃兰海的垂直运动。研究物理学。海洋学报，21：811-823。另请参阅tintore，filrespcokri/cokri2目的cokri：点或块D尺寸的协同克里格法。cokri2：从cokri称为协同克里格法的功能。简介cokri/cokri2目的cokri：点或块D尺寸的协同克里格法。cokri2：从cokri称为协同克里格法的功能。简介[x0s,s,sv,id,l]=cokri(x,x0,model,c,itype,avg,block,nd,ival,nk,rad,ntok,d)[x0s,s,id,l,k0]=cokri2(x,x0,id,model,c,sv,itype,avg,ng,d)说明协同克里格法是指用多个变量的克里格。当cokri只有一个被称为在一个时间变量，其结果将是那些简单克里格，普通克里格法，普遍克里格，点克里格或块克里格。更多细节，可以发现在Marcotte的论文中(1991年)。可用克里格的选项是：1。简单的协同克里格法2。一个nonbias条件与普通协同克里格法(Isaaks和Srivastava)3。与Pnonbias条件的普通协同克里格法4。1阶漂移普遍协同克里格法5。2阶漂移普遍协同克里格法99.cokriging不执行，只计算SV可用的变异函数模型：1。金块的影响2。指数模型3。高斯模型4。球状模型5。线性模型6。二次模型电源(H3)模型对数的SINC(H)贝塞尔[JO(H)]Bessel[Jo(h)]exp(-h)*cos(dh)exp(-h)*Jo(dh)exp(-h2)*cos(dh)exp(-h2)*Jo(dh)exp(-h2)*(1-dh2)1-3*min(h,1)2+2*min(h,1)3h2*log(max(h,eps))可以添加新的模型，很容易因为使用eval函数模型计算。输入变量：【译文】X：数据矩阵[XYZVARIVAR2...]X0：网格坐标西医子]模式：[模型，(H=R/A),旋转角度]。无旋转角度为各向同性分布。例如：型号=[1104月30日]表示的分布是各向同性的和它是代表一个范围10％的金块效果加上一个球状模型范围30。C：波幅模型克里格itype:选项块：向量(1x深)，估计块大小;点协同克里格法：块=的(1，D)ND：矢量(1个D)，块协同克里格法的离散网格;点协同克里格法：ND=的(1，D)ival:用于交叉验证的代码。0：没有交叉验证1：交叉验证是由删除一个变量时间在一个给定的的位置。2：交叉验证是通过消除所有的变量在执行给定的位置。NK：最近的邻居在X矩阵中使用的协同克里格法的数量。RAD：搜索半径。ntok：X0点将由ntok网格点的克里金。D：模型系数。这一系数已被添加到原来的Marcotte功能。警告：在cokri，模型是定义在H=R/A.在variogr，因变量为R，因此d=B*A(见variogr)。Inputvariables:输入变量：【原文】x:datamatrix[xyzvar1var2...]x0:gridcoordinates[xiyizi]model:[models,a(h=r/a),rotationangles].Norotationangleisrequiredforanisotropicdistribution.Example:model=[110;430]meansthatthedistributionisisotropicandthatitisrepresentedbyanuggeteffectofrange10plusasphericalmodelofrange30.c:amplitudesofthemodelsitype:krigingoptionblock:vector(1xD),givingthesizeoftheblocktoestimate;forpointcokriging:block=ones(1,D)nd:Vector(1xD),givingthediscretizationgridforblockcokriging;forpointcokriging:nd=ones(1,D)ival:Codeforcross-validation.0:nocross-validation1:cross-validationisperformedbyremovingonevariableatatimeatagivenlocation.2:cross-validationisperformedbyremovingallvariablesatagivenlocation.nk:numberofnearestneighborsinxmatrixtouseinthecokriging.rad:searchradius.ntok:pointsinx0willbekrigedbygroupsofntokgridpoints.d:modelcoefficients.ThiscoefficienthasbeenaddedtotheoriginalMarcotte'sfunction.Warning:Incokri,modelsaredefinedintermsofh=r/a.Invariogr,thedependentvariableisrandhenced=b*a(seevariogr).输出变量【译文】xOs:克里金数据矩阵XO地点。S：克里金在xO处的数据方差矩阵。SV：每个变量的方差。ID，看到Marcotte升(1991)输出变量【原文】Outputvariablesx0s:krigeddatamatrixatx0locations.s:krigeddatavariancematrixatx0locations.sv:varianceofeachvariable.id,lseeMarcotte(1991)参考Marcotte，D.1991。Cokrigeage用matlab。计算机及地质。17(9)：1265年至1280年。另请参阅cokri2,variogr，反，是指confint目的置信区间。简介[K2，K1]=confint(G，M，S2)说明置信区间的结构功能CONF{K28方差8K1}(1)结构的功能是衡量一个给定变量的函数的方差距离。在这种情况下的置信区间的估计是：(1)。K1=(N-1)*S2/C1K2=(N-1)*S2/C2(2)其中n=样本大小=M+1M=自由度$2=样本方差Cl和C2是由方程的解F(C1)=(1-G)/2(C2)=(1+G)/2(3)G=置信水平(95％，99％或类似)F=X分酉己解决，获得了功能chitable(池工具箱)。参考文献登曼，吉隆坡和HJ斐然(1985年)。相关秤，目标定位和统计检验Geostrophy大陆架。研究的三月分辨率，43：517-539。Kreyszig，E，1988年，高等工程数学，第六版，威利父子，纽约，p.1252另请参阅chitabledavis戴维斯目的点克里格使用戴维斯的方程小"其中A他1)r(用1)y(亦1)17(/?12).../(Aitti)722)....yt//2m)yU加2)...y\h耳册I1...111.ir=~Wi~WiandB-贰心p)1wm0A1Y(hik)是由距离h分隔的半方差样本对hik。非偏见条件要求的Wi的总和等于1。在这种情况下，一人多度自由必须引入与使用拉格朗日乘子久为了尽量减少估计误差。简介［Zp,Sp］=davis(data,xO,model,a,d,c,A)说明可用的模型：1。金块的影响2。指数模型3。咼斯模型4。球状模型5。线性模型6。二次模型7。功率模型(H"8。对数模型9。SINC(H)10。贝塞尔［JO(H)］11。EXP(-H)*COS(DH)12。EXP(-H)*JO(DH)13。EXP(-H2)*COS(DH)14。EXP(-H2)*JO(DH)15。EXP(-H2)x(1-DH2)输入变量：data:数据［XY变量］X0：网格坐标［西医］model:变异函数模型a：半变异范围d：模型系数(系数variogrb不同，与D相同，在cokri)c：模型振幅A：Y(HIK)矩阵，如果已经计算;若否，忽略的输入。输出变量：Zp：克里金数据矩阵X0职位Sp：在xO位置克里金数据的方差参考1986年戴维斯，JC统计和数据分析，在地质，第2版，威利父子，纽约，289页。deplie目的矢量矩阵转换。简介mat二deplie（y,nx,ny）说明改造一个向量y成为一个矩阵大小为nyxnx的mat。输入变量：y：行或列向量nx：在矩阵垫中的列数ny：在矩阵垫中的行数输出变量：mat:矩阵示例y=［（i：io）（1:10）（110）］；deplie(y,10,3)=123456789101234567891012345678910另请参阅kregridfilresp目的巴恩斯的滤波器响应在波长域。简介R=filresp(c,g)说明巴恩斯在滤波器响应的波长域为：R二傀(1+2?畀一盅)R°=exp(-7T4c//I；)其中^,C和G滤波器的参数，是波长。示例f2=filresp(200,0.6);Filterresponseinthe-wavelengOidomacm0asuodsauglzWavefength{distano?unit、参考文献巴恩斯，SL,1973年。中尺度目的地图使用加权时间序列分析意见。NOAA的技术。备忘录。广深港高速铁路NSSL-62，60页马多克斯，类风湿性关节炎，1980年。一个分离宏观和中尺度的目的技术气象数据的特点。周一WEA。牧师，108，1108:1121。另请参阅巴恩斯，tintorefitvario/fitvari2/fun目的fitvario：variogr优化。fitvari2：variogr2优化(没有优化工具箱)。fun：从fitvario估计变异函数简介fitvario(modeldata,a)fitvario(modeldata,a,b)fitvari2(modeLdata,a,b,C)f=fun(lam,data)说明最小二乘拟合半变异函数模型系数A,B和C输入变量：model:模型的类型(见variogr)data(:,1):X轴(距离)(GAM(:,1))data(:,2):Y轴(方差)(GAM(:,2))A,B,C：variogr系数A,B,C的初始值原文如下：data(:,1):x-axis(distanee)(gam(:,1))data(:,2):y-axis(varianee)(gam(:,2))a,betC:startingvaluesofeoeffieientsa,bandCofvariogr输出：图形输出显示半方差作为距离的函数曲线。实验结果显示为符号。纯行给出最适合的模型选择。A,B和C的最佳值也显示在图上。示例Exampler=(o：ioy；x=1-exp(-r);fitvario(2.[rx],2)fitvari2(2.[rx],2F031.5)BQL1SLIBAJEd迫□iyiflrce另请参阅gam2,gamv2,gamv2uvgam3,gamv3目的vario2dr,vario2di,var2diuv,vario3dr和vario3di称为MEX文件。简介[np,gam,hmrtmrhv,tv]=gam2(nlag,nx,ny,ndir,ixd,iyd,vr;tmin;tmax,nvarg,ivtail,ivhead,ivtype.nvar)[np,dis,gam,hm,tm,hv:tv]=gannv2(nd,x,y:vr,tmin,tmax,nlag,xlag,xltol,ndirpazm,atol,bancb