%本函數(shù)為截斷函數(shù)，對任意的不在(-pi,pi)內的位相均可進行截斷
%dPw為待截斷位相，M為抽樣點數(shù)
% dPww為位于(-pi,pi)內的截斷位相

function [dPww] = wrap0(dPw)
[M,N] = size(dPw);

for i = 1:M
    for j = 1:N 
        while dPw(i,j) >= pi
            dPw(i,j) = dPw(i,j)-2*pi;
        end
        while dPw(i,j) < -pi
            dPw(i,j) = dPw(i,j)+2*pi;
        end  
        dPww(i,j) = dPw(i,j);
    end
end