feature de centrado de ciclo de histéresis

[jsilveyra]

agregar feature de centrado del ciclo de histéresis, después de convertirlo a M vs H en A/M y antes de calcular la curva anhisterética. No debería correrse con una curva input anhisterética.

Podría ofrecerse con un box con tic que dijera "Center hysteresis loop" y que por default esté marcada:

Este es el algoritmo de centrado (M y H left corresponden a la rama de la izq del ciclo de histéresis y M y H right a la de la derecha)

% Center hysteresis loop

% Center vertical axis
Mneg = min(Mleft);
Mpos = max(Mright);
Mshift = (Mpos+Mneg)/2;
Mleft = Mleft - Mshift;
Mright = Mright - Mshift;

% Center horizontal axis
Hcneg = interp1(Mleft, Hleft, 0, 'linear'); % linear interpolation
Hcpos = interp1(Mright, Hright, 0, 'linear'); % linear interpolation
Hshift = (Hcneg+Hcpos)/2;
Hleft = Hleft - Hshift;
Hright = Hright - Hshift;

Hshifted = [flip(Hleft) Hright];
Mshifted = [flip(Mleft) Mright];