%本函數(shù)為截?cái)嗪瘮?shù)，對(duì)任意的不在(-pi,pi)內(nèi)的位相均可進(jìn)行截?cái)?%dPw為待截?cái)辔幌?，M為抽樣點(diǎn)數(shù)
% dPww為位于(-pi,pi)內(nèi)的截?cái)辔幌?
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