Floating point implementation: starting with an initial value of matrix P0 and vector U0, then proceeding to the algorithm described below with measurement y(??i) and basis r=(1 cos ??i sin ??i cos 2??i sin 2??i)T, where P is a matrix of dimension M???M and U and r are vectors of dimension M. M is the dimension of the basis function.