block = r(l : min(l + blocksize - 1, length(r)));
sum_M = sum(block .^ M);
- phi_est = angle(sum_M) / M; % assume phase of 0 symbol is 0.
+ %% if phase of 0 symbol is 0, use:
+ phi_est = angle(sum_M) / M;
+ %% if phase of 0 symbol is pi/M, use:
+ %% phi_est = angle(-sum_M) / M;
if l > 1
%% phase unwrapping