detection/blending.m

   1 function T=blending(Tf,Ta,T0,T1)
   2 % T=blending(Tf,Ta,T0,T1)
   3 % fire arrival times Tf at T0 and Ta at T1 into T
   4
   5 if all(Tf>=T0),
   6     T=Ta;
   7 else
   8     p=1.1;
   9     f=max(Tf-T0,0).^p;
  10     a=max(T1-Ta,0).^p;
  11     T=(f.*Ta+a.*Tf)./(a+f);
  12 end
  13 if any(isnan(T(:))),warning('T0 too large'),end
  14 end