qam16.m

星座图QAM4.大家参考参考啊！有更好的请邮箱联系

clear all
QAM16=[ ];
for w=1:4
    a=(-5+2*w)*2/3*(2*sqrt(2)-sqrt(5))
    for r=1:4
        b=(-5+2*r)*2/3*(2*sqrt(2)-sqrt(5))
        c=a+b*j
        QAM16=[QAM16,c]
    end
end

dB=20;
p=10^(dB/10);


NN=10000;
err=0;

for t=1:NN
    
    s=unidrnd(16);

    x=QAM16(s);
    h=1/2*randn(1,1)+sqrt(-1)*1/2*randn(1,1);
    w=1/2*randn(1,1)+sqrt(-1)*1/2*randn(1,1);
    
    y=p*h*x+w;
    
    dd=100;
    
    T=0;
    MM=100;
    for t=1:16
        z=QAM16(t);
        f=abs(y-p*h*z);
        if f<MM
            MM=f;
            T=t;
        else MM=MM;
            T=T;
        end
    end
    
    xx=QAM16(T);
   
    aa=abs(x-xx);
    if  aa == 0
        err=err;
    else err=err+1;
    end
    
end

pro_err=err/NN