Kriging插值MATLAB代码

程序可以运行Kriging插值的算法过程，其中包含几种不同的半变异函数可供选择使用。

function [xfvalexitflagoutput] = fminsearch（funfcnxoptionsvarargin）
%FMINSEARCH Multidimensional unconstrained nonlinear minimization （Nelder-Mead）.
%   X = FMINSEARCH（FUNX0） starts at X0 and finds a local minimizer X of the
%   function FUN. FUN accepts input X and returns a scalar function value 
%   F evaluated at X. X0 can be a scalar vector or matrix. 
%
%   X = FMINSEARCH（FUNX0OPTIONS）  minimizes with the default optimization
%   parameters replaced by values in the structure OPTIONS created
%   with the OPTIMSET function.  See OPTIMSET for details.  FMINSEARCH uses
%   these options: Display TolX TolFun MaxFunEvals and MaxIter. 
%
%   X = FMINSEARCH（FUNX0OPTIONSP1P2...） provides for additional
%   arguments which are passed to the objective function F=feval（FUNXP1P2...）.
%   Pass an empty matrix for OPTIONS to use the default values.
%   （Use OPTIONS = [] as a place holder if no options are set.）
%
%   [XFVAL]= FMINSEARCH（...） returns the value of the objective function
%   described in FUN at X.
%
%   [XFVALEXITFLAG] = FMINSEARCH（...） returns a string EXITFLAG that 
%   describes the exit condition of FMINSEARCH.  
%   If EXITFLAG is:
%     1 then FMINSEARCH converged with a solution X.
%     0 then the maximum number of iterations was reached.
%   
%   [XFVALEXITFLAGOUTPUT] = FMINSEARCH（...） returns a structure
%   OUTPUT with the number of iterations taken in OUTPUT.iterations.
%
%   Examples
%     FUN can be specified using @:
%        X = fminsearch（@sin3）
%     finds a minimum of the SIN function near 3.
%     In this case SIN is a function that returns a scalar function value 
%     SIN evaluated at X.
%
%     FUN can also be an inline object:
%        X = fminsearch（inline（‘norm（x）‘）[1;2;3]）
%     returns a minimum near [0;0;0].
%
%   FMINSEARCH uses the Nelder-Mead simplex （direct search） method.
%
%   See also OPTIMSET FMINBND @ INLINE.

%   Reference: Jeffrey C. Lagarias James A. Reeds Margaret H. Wright
%   Paul E. Wright “Convergence Properties of the Nelder-Mead Simplex
%   Method in Low Dimensions“ SIAM Journal of Optimization 9（1）: 
%   p.112-147 1998.

%   Copyright 1984-2002 The MathWorks Inc. 
%   $Revision: 1.21 $  $Date: 2002/04/08 20:26:45 $
%
% Olivier Salvado Case Western Reserve University June04
%   Modified to work on Cost function smooth on a high scale but rough on a
%   small scale.
%   （1） the DiffMinChange option is taken into account to limit contraction
%   （2） two new options are added to initialize the starting trials （usual_deltazero_term_delta）.
%   （3） there is a possibility to display patches for the cases with 3 trials


defaultopt = struct（‘Display‘‘notify‘‘MaxIter‘‘200*numberOfVariables‘...
   ‘MaxFunEvals‘‘200*numberOfVariables‘‘TolX‘1e-4‘TolFun‘1e-4‘DiffMinChange‘1e-6...
   ‘usual_delta‘0.05...             % 5 percent deltas for non-zero terms
   ‘zero_term_delta‘0.00025

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件       49461  2014-04-09 15:09  kriging\ans.mat
     文件       10900  2014-04-21 16:05  kriging\fminsearch.asv
     文件       10917  2014-04-21 16:08  kriging\fminsearch.m
     文件        5066  2014-04-03 11:02  kriging\fminsearchbnd.m
     目录           0  2014-04-21 15:56  kriging\hehe\
     文件        5066  2014-04-03 11:02  kriging\hehe\fminsearchbnd.m
     文件        8434  2014-04-03 10:38  kriging\hehe\fminsearchbnd1.m
     文件        5066  2014-04-03 11:02  kriging\hehe\fminsearchbnd3.m
     文件       10636  2014-04-03 11:01  kriging\hehe\fminsearchbnd4.m
     文件          81  2014-04-03 11:10  kriging\hehe\readme.txt
     文件        6374  2012-10-15 11:11  kriging\kriging.m
     文件        1671  2014-04-19 13:43  kriging\run.asv
     文件        1717  2014-04-21 16:08  kriging\run.m
     文件       83290  2014-04-03 17:25  kriging\t11.mat
     文件        1671  2014-04-19 15:49  kriging\Untitled4.m
     文件       12718  2014-04-09 14:06  kriging\variogram.m
     文件       18297  2014-04-03 11:04  kriging\variogramfit.m