-
Notifications
You must be signed in to change notification settings - Fork 0
/
regDenom.m
20 lines (16 loc) · 836 Bytes
/
regDenom.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% regDenom.m %%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
% function par=regDenom(point,step,par,tune)
% construct regularize denomiator
function par=regDenom(point,step,par,tune)
par.denom=par.gamma-step.q'*step.r;
% regularize denominator
e2f=abs(point.fnew-point.f)+step.agood*(abs(point.g)'*abs(step.p));
dcor=tune.DeltaH*(e2f/step.agood^2+abs(step.q)'*abs(step.r));
if par.denom>=0, par.denom=par.denom+dcor;
else, par.denom=par.denom-dcor;
end;
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%