function kpc=align(f)
%ALIGN Real alignment
% ALIGN(F) return a constant matrix KPC which minimizes
% the least squares error of
% f*kpc - diag(exp(sqrt(-1)*ai))
% over b
function A=align(V)
% A=ALIGN(V) returns a constant matrix A which is the real alignment
% of the INVERSE of the complex input matrix V, i.e.
%
% A=arg min_{A,theta}||V*A-diag(ex
