once_sim.asv

有限反馈的仿真程序

function[sum_rate]=once_sim(H,T)
P=1;
[K,M]=size(H);
for pwr=1:21
    for i=1:M
        temsum=0;
        for j=1:K
            if j~=i
                temsum=temsum+(abs(H(i,:)*T(:,j)/norm(T(:,j))))^2;
            end
        end
        SINR(i)=(P/M)*(abs(H(i,:)*T(:,i)/norm(T(:,i))))^2/(1+(P/M)*temsum);
        R(i)=log2(1+SINR(i));
    end
    sum_rate(pwr)=sum(R);
    P=P+100;
end