feedback2.m

有限反馈的仿真程序

clear all;
M=5;K=5;
Pdb=0:5:30;
P=10.^(Pdb/10);
for pwr=1:7
    NFB=((M-1)/3)*round(10*log10(P(pwr)));
    H=channel(K,M);
    W=codebook(M,2^NFB);
    Quantized_H=rand_quantize(H,W);
    T=Quantized_H'*inv(Quantized_H*Quantized_H');
    for series=1:1
        R(series)=once_sim(H,T,P(pwr));
    end
    mean_R(pwr)=mean(R);
end
x=1:20:1001;
plot(10*log10(x),mean_R)