main.asv

有限反馈的仿真程序

clear all;
M=4;K=4;
total_NFB=24;
Pdb=0:5:20;
P=10.^(Pdb/10);
for num_K=2:M;   
    for series=1:2
        series
        for pwr=1:length(Pdb)
            H=channel(K,M);
            W=codebook(M,2^round(total_NFB/num_K));
            Quantized_H=rand_quantize(H,W);
            T=Quantized_H'*inv(Quantized_H*Quantized_H');
            R(series,pwr)=once_sim(num_K,H,T,P(pwr));
        end
    end
    mean_R(num_K,:)=mean(R);
end
num
hold on;
plot(Pdb,mean_R(2,:),'r->');
plot(Pdb,mean_R(3,:),'g-o');
plot(Pdb,mean_R(4,:),'b-<');