function DEC_X=judge(y)%最大似然檢測子函數，y1，y2為待檢測信號，nt為，返回值為與待檢測信號相隔最近的調制信號星座圖中對應的信號值
DEC_X=zeros(1,4);
for i=1:4
    for m=1:4
        d(i,m)=abs(y(i)-exp(j*2*pi/4*(m-0.5))).^2;%求調制信號星座圖中各調制信號與待檢測信號的距離平方
    end
end
%進行判決
for i=1:4
    [z,i1]=min([d(i,1) d(i,2) d(i,3) d(i,4)]);%z1為最小值，i1為最小值出現的位置
    DEC_X(i)=exp(j*2*pi/4*(i1-0.5));
end