峰值检测代码

峰值检测matlab，经测可以使用。喜欢就赶紧下载吧

function [tempajrainsnow]=adjmeteott（tempprecipaltmoylapsealtstatT...
    corrprecipaltstatPtcrValtcrIntVal）
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This code is part of the hydrological model GSM-SOCONT
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% MATLAB code written by Bettina Schaefli
% Ecole Polytechnique F閐閞ale de LausanneSwitzerland
% E-mail: bettina.schaefli@epfl.ch
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This function interpolates input for GSM-SOCONT
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Input:
%   For the names of input variables and constants see initgsmsocont.m
% Output:
%   tempaj rain snow: 
%       Matrices containing 1 time series per band
%       1 colonne per band format NxM where N the
%       number of time steps and M the number of bands
%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Note: the separation between snow- and rainfall is based on a fuzzy
% threshold not a single threshold as in the original paper Schaefli et
% al.2005
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%


noband=length（altmoy）;
N=length（temp）;

tempaj  =zeros（Nnoband）;
snow    =zeros（Nnoband）;
rain    =zeros（Nnoband）;

tcrinf=tcrVal-tcrIntVal/2;
tcrsup=tcrVal+tcrIntVal/2;


for j=1:noband
        tempaj（:j）     =temp+lapse*（altmoy（j）-altstatT）/100;
        precipaj（:j）   =precip+corrprecip/100*precip*（altmoy（j）-altstatP）/100;
        
        indneg=find（precipaj（:j）<0）;
        if ~isempty（indneg）
            ‘100% precipitation reduction for band‘
            j
            ‘see adjmeteott.m‘
        end
        precipaj（indnegj）=0;   
    
        %separation snow-rain
        kinf=find（tempaj（:j）<=tcrinf）;
        snow（kinfj）=precipaj（kinfj）;        
        kint=find（and（tcrinf        % the above is empty if tcrIntVal==0
        % in Matlab an empty index is not a problem nothing happens
        snow（kintj）=precipaj（kintj）.*（tcrsup-tempaj（kintj））/tcrIntVal;
        rain（:j）=precipaj（:j）-snow（:j）;
end

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% This code is part of the hydrological model GSM-SOCONT
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Citation:
%   Schaefli B. Hingray B. Niggli M. and Musy A. 2005. 
%       A conceptual glacio-hydrological model for high mountainous catchments.     
%       Hydrology and Earth System Sciences 9: 95 - 109.                                                                                             
%       http://www.hydrol-earth-syst-sci.net/9/95/2005/hess-9-95-2005.pdf
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% Copyright and Disclaimer
%%%%%%%%%%%%%%%%%%%%

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       4051  2013-09-11 10:02  GSMSOCONT_V2.5\GSMSOCONT_V2.5\adjmeteott.m

     文件       3759  2013-09-11 10:02  GSMSOCONT_V2.5\GSMSOCONT_V2.5\courbehypsoinv.m

     文件      34423  2003-06-03 18:26  GSMSOCONT_V2.5\GSMSOCONT_V2.5\exETP.txt

     文件       3030  2003-06-04 09:20  GSMSOCONT_V2.5\GSMSOCONT_V2.5\exhypsoglacier.txt

     文件       3030  2003-06-04 09:18  GSMSOCONT_V2.5\GSMSOCONT_V2.5\exhypsononglacier.txt

     文件      16358  2003-06-03 18:16  GSMSOCONT_V2.5\GSMSOCONT_V2.5\exPrec.txt

     文件      72324  2013-09-11 09:58  GSMSOCONT_V2.5\GSMSOCONT_V2.5\exqsim.txt

     文件      21530  2003-06-03 18:16  GSMSOCONT_V2.5\GSMSOCONT_V2.5\exTemp.txt

     文件       5855  2013-09-11 10:33  GSMSOCONT_V2.5\GSMSOCONT_V2.5\gsmsocont.m

     文件       3542  2013-09-11 09:42  GSMSOCONT_V2.5\GSMSOCONT_V2.5\gsmsocparprior.m

     文件       4551  2013-09-11 09:53  GSMSOCONT_V2.5\GSMSOCONT_V2.5\initgsmsocont.m

     文件       2484  2013-09-11 10:02  GSMSOCONT_V2.5\GSMSOCONT_V2.5\linres.m

     文件       1154  2013-09-11 09:47  GSMSOCONT_V2.5\GSMSOCONT_V2.5\nashbias.m

     文件       6128  2013-09-11 10:28  GSMSOCONT_V2.5\GSMSOCONT_V2.5\RunGsmSocont.m

     文件      12070  2013-09-11 10:39  GSMSOCONT_V2.5\GSMSOCONT_V2.5\socontplus.m

     文件       8689  2013-09-11 10:31  GSMSOCONT_V2.5\GSMSOCONT_V2.5\transpq.m

     目录          0  2013-09-11 04:57  GSMSOCONT_V2.5\GSMSOCONT_V2.5

     目录          0  2013-09-15 17:22  GSMSOCONT_V2.5

----------- ---------  ---------- -----  ----

               202978                    18