feedback1.asv

有限反馈的仿真程序

clear all;
M=5;K=5;NFB=10;
H=channel(K,M);
W=codebook(M,2^NFB);
Quantized_H=rand_quantize(H,W);
FB_T=Quantized_H'*inv(Quantized_H*Quantized_H');
for series=1:100
    R(series,:)=once_sim(H,FB_T);
end
mean_R=mean(R);
x=1:50:1001;
plot(10*log10(x),mean_R,'b-');
hold on;
%%%%%%%%%%%%%%%%%%%%%%%%%%
ZF_T=H'*inv(H*H');
for series=1:100
    R(series,:)=once_sim(H,ZF_T);
end
mean_R=mean(R);
plot(10*log10(x),mean_R,'g-');
%%%%%%%%%%%%%%%%%%%%%%%%
C=M*log2(1+x);
plot(10*log10(x),C,'r-')