feedback1.m

有限反馈的仿真程序

clear all;
M=4;K=4;NFB=2;
Pdb=0:2:30;
P=10.^(Pdb/10);
for series=1:2
    series
    for pwr=1:length(Pdb)
        H=channel(K,M);
        W=codebook(M,2^NFB);
        Quantized_H=rand_quantize(H,W);
        FB_T=Quantized_H'*inv(Quantized_H*Quantized_H');
        R(series,pwr)=once_sim(H,FB_T,P(pwr));
    end
end
mean_R=mean(R);
plot(Pdb,mean_R,'b-');
hold on;
%%%%%%%%%%%%%%%%%%%%%%%%
for series=1:500
    for pwr=1:length(Pdb)
        H=channel(K,M);
        ZF_T=H'*inv(H*H');
        R(series,pwr)=once_sim(H,ZF_T,P(pwr));
    end
end
mean_R=mean(R);
plot(Pdb,mean_R,'g-o');
%%%%%%%%%%%%%%%%%%%%%%%%
for pwr=1:length(Pdb)
    C(pwr)=M*log2(1+P(pwr));
end
plot(Pdb,C,'r-+')