detection/like1.m

   1 function [v0,v1]=like1(dw,t,stretch)
   2 % [v0,v1]=like1(dw,t,stretch)
   3 % dw > 0 if fire detected, <0 if not, 0 if cloud
   4 % t = time since fire arrival (days)
   5 Tmin=stretch(1);Tmax=stretch(2);Tneg=stretch(3);Tpos=stretch(4);
   6
   7
   8 tneg = (t-Tmin)./Tneg;
   9 tpos = (t-Tmax)./Tpos;
  10 v0y = -tneg.*tneg.*(t<Tmin)-tpos.*tpos.*(t>Tmax);
  11
  12 g0 = @(x) -0.25.*x.*x.*x+0.75*x-0.5;
  13 % v0n1 = (t<Tmin-2*Tneg);
  14 v0n2 = g0(-tneg-1).*(Tmin-2*Tneg<=t & t<=Tmin);
  15 v0n3 =  - (Tmin < t & t<Tmax);
  16 v0n4 = g0(tpos-1).*(Tmax <= t & t <= Tmax+2*Tpos);
  17 % v0n5 = (t>Tmax+2*Tpos);
  18 % v0n = v0n1 + v0n2 + v0n3 + v0n4 + v0n5;
  19 v0n =  v0n2 + v0n3 + v0n4 ;
  20 v0 = -dw.*(dw<0).*v0n + dw.*(dw>0).*v0y;
  21
  22 v1y = -2.*tneg.*(t<Tmin)-2.*tpos.*(t>Tmax);
  23
  24 g1 = @(x) -1.5*x.*x+1.5;
  25 v1n = - g1(-tneg-1).*(Tmin-2*Tneg<=t & t<=Tmin) ...
  26     + g1(tpos-1).*(Tmax <= t & t <= Tmax+2*Tpos);
  27 v1 = -dw.*(dw<0).*v1n + dw.*(dw>0).*v1y;
  28
  29 end
  30
  31
  32
  33