TP2LC Calculates the number of upcrossings from the turning points. CALL: lc = tp2lc(TP,def,plotflag,sa); lc = a two column matrix with levels and number of upcrossings. [mx2] TP = the turning points. [nx2] def = 1, only upcrossings. 2, upcrossings and maxima (default). 3, upcrossings, minima, and maxima. 4, upcrossings and minima. plotflag = 0, no plotting 1, plot the number of upcrossings overplotted with Rice formula for the crossing intensity for a Gaussian process (default). sa = standard deviation of the process (Default estimates it from the number of upcrossings) See also lcplot
Calculates the number of upcrossings from a cycle count | |
Calculates min2Max and Max2min cycles from a sequence of turning points | |
Display message and abort function. |
% CHAPTER4 contains the commands used in Chapter 4 of the tutorial | |
Script to computer exercises 1 | |
Quick test of the routines in module 'cycles' | |
Extrapolates a sequence of turning points. |
001 function lc=tp2lc(tp,def,plotflag,sa) 002 %TP2LC Calculates the number of upcrossings from the turning points. 003 % 004 % CALL: lc = tp2lc(TP,def,plotflag,sa); 005 % 006 % lc = a two column matrix with levels and number of upcrossings. [mx2] 007 % TP = the turning points. [nx2] 008 % 009 % def = 1, only upcrossings. 010 % 2, upcrossings and maxima (default). 011 % 3, upcrossings, minima, and maxima. 012 % 4, upcrossings and minima. 013 % 014 %plotflag = 0, no plotting 015 % 1, plot the number of upcrossings overplotted 016 % with Rice formula for the crossing intensity 017 % for a Gaussian process (default). 018 % 019 % 020 % sa = standard deviation of the process 021 % (Default estimates it from the number of upcrossings) 022 % 023 % See also lcplot 024 025 % Tested on Matlab 5.3 026 % 027 % History: 028 % Created by PJ (Pär Johannesson) 09-Jan-2000 029 030 % Check input arguments 031 032 ni = nargin; 033 no = nargout; 034 error(nargchk(1,4,ni)); 035 036 if ni<2, def=[]; end 037 if ni<3, plotflag=[]; end 038 if ni<4, sa=[]; end 039 040 % Get min-max cycles 041 mM = tp2mm(tp); 042 % Get level crossings 043 %lc = cc2lc(mM,def,plotflag,sa); 044 lc = cc2lc(mM,def,0,sa); 045
Comments or corrections to the WAFO group