detection/like3.m

   1 function [v0,v1]=like3(dw,t,stretch)
   2
   3 % stretch:
   4 Tmin=stretch(1);Tmax=stretch(2);Tneg=stretch(3);Tpos=stretch(4);
   5
   6 max_like=log(0.9);
   7
   8 % get p0 = data likelhood of 1 as a function of t, p1=derivative
   9
  10 slope_neg=0.5;
  11 slope_pos=0.1;
  12
  13 p0=max_like*(t>=Tmin & t<=Tmax) + ...
  14    (max_like+(t-Tmin)*slope_neg).*(t<Tmin) + ...
  15    (max_like+(Tmax-t)*slope_pos).*(t>Tmax);
  16
  17 p1= (t<Tmin)*slope_neg - (t>Tmax)*slope_pos;
  18
  19 % everything else is derived from that
  20
  21 % p0 + p1 = 1
  22 % p0 = 1 - p1
  23 % logp0 = log(1 - p1) = log(1 -exp(logp0))
  24
  25 n0=log(1-exp(p0));
  26 n1=-p1.*exp(p0)./(1-exp(p0)).^2;
  27
  28 v0 = dw.*p0.*(dw>0) - dw.*n0.*(dw<0);
  29 v1 = dw.*p1.*(dw>0) - dw.*n1.*(dw<0);
  30
  31
  32 end
  33
  34