main.m

来自「有限反馈的仿真程序」· M 代码 · 共 33 行

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:50
        series
        for pwr=1:length(Pdb)
            H=channel(num_K,num_K);
            W=codebook(num_K,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
    num=0;
    for i=1:series
        if sum(isnan(R(i,:)))==0
            num=num+1;
            valid_R(num,:)=R(i,:);
        end
    end
    mean_R(num_K,:)=mean(valid_R);
end
hold on;
grid on;
plot(Pdb,mean_R(2,:),'r->');
plot(Pdb,mean_R(3,:),'g-o');
plot(Pdb,mean_R(4,:),'b-<');
legend('Num-User=2','Num-User=3','Num-User=4');
xlabel('SNR(dB)');
ylabel('Sum-rate(bps/Hz)');