qub_feedback.m

有限反馈的仿真程序

clear all;
M=4;K=6400;
epsi=0.45;
Pdb=6:2:20;
P=10.^(Pdb/10);
for series=1:2
    for pwr=1:length(Pdb)
        pwr
        NFB=round((M-1)*log2(P(pwr))+8.68-log2(K));
        H=channel(K,M);
        W=codebook(M,2^NFB);
        Quantized_H=qub_quantize(H,W);
        SINR_sus=cal_sinr_sus(H,Quantized_H,P(pwr),epsi);
        user_sinr=userselect(H,Quantized_H,SINR_sus,epsi);
        SINR_sinr=cal_sinr(user_sinr,H,Quantized_H,P(pwr));
        R_sinr(series,pwr)=once_sim(SINR_sinr,M);
    end
end
mean_R_sinr=mean(R_sinr);
plot(Pdb,mean_R_sinr,'b-+');