LC2DPLUS Upper bound for total damage from level crossings. Calculates the upper bound for the total damage/damage intensity using the crossing spectrum/dowcrossing intensity. CALL: D = lc2dplus(cross,beta,n); where D = the upper bound for the total damage/damage intensity, cross = a two column matrix with levels u and corresponding downcrossing intensities, beta = a vector with beta-values, n = (optional input argument) the number of slice levels (default = 500).
Calculates the cycle count which has the highest damage |
001 function D=lc2dplus(cross,beta,num_slice) 002 %LC2DPLUS Upper bound for total damage from level crossings. 003 % 004 % Calculates the upper bound for the total damage/damage intensity using 005 % the crossing spectrum/dowcrossing intensity. 006 % 007 % CALL: D = lc2dplus(cross,beta,n); 008 % 009 % where 010 % D = the upper bound for the total damage/damage intensity, 011 % cross = a two column matrix with levels u and corresponding 012 % downcrossing intensities, 013 % beta = a vector with beta-values, 014 % n = (optional input argument) the number of slice levels 015 % (default = 500). 016 017 % Copyright 1993, Mats Frendahl, Dept. of Math. Stat., University of Lund. 018 019 if nargin<3 020 num_slice=500; 021 end 022 023 [cc,delta]=down2cc(cross,num_slice); 024 025 amplitudes=cc(:,1)-cc(:,2)/2; 026 cc_length=length(cc); 027 028 delta=max(cross(:,2))/num_slice; 029 030 deltas=[delta*ones(1,cc_length-1) .75*delta]; 031 for i=1:length(beta) 032 D(i)=deltas*amplitudes.^beta(i); 033 end 034
Comments or corrections to the WAFO group