Matlab绘制彩色星座图_改进

function Get_Color_Constellation（wavein）
    N = 4;
    Nbeside = 20;%左右两侧共2*Nbeside个点
    indexdata = zeros（length（wavein）N）;
    angledata = linspace（0pi/2N）;
    for m = 1:N
        data = real（exp（1j*angledata（m）））*real（wavein）+...
            imag（exp（1j*angledata（m）））*imag（wavein）;
        [~sortdata] = sort（data）; 
        indexdata（:m） = sortdata;
    end
    dis = zeros（length（indexdata）N）;
    for kk = 1:N
        for k = 1:length（indexdata）
            indexbeg = k-Nbeside;
            if indexbeg < 1
                indexbeg = 1;
            end
            indexend = indexbeg+Nbeside*2+1;
            if indexend > length（data）
                indexend = length（data）;
                indexbeg = length（data）-2*Nbeside;
            end  
            datatemp = wavein（indexdata（indexbeg:indexendkk））-wavein（indexdata（kkk））;
            [t1t2]=sort（abs（datatemp））;
            dis（indexdata（kkk）kk） = t1（2）;
        e