cal_sinr.m

有限反馈的仿真程序

function[SINR]=cal_sinr(user,H,Quantized_H,P)
[K,M]=size(H);
for i=1:M
    userset_H(i,:)=Quantized_H(user(i),:);
end
T=userset_H'*inv( userset_H* userset_H');
for i=1:M
    temsum=0;
    for j=1:M
        if j~=i
            temsum=temsum+(abs(H(user(i),:)*T(:,j)/norm(T(:,j))))^2;
        end
    end
    SINR(i)=(P/M)*(abs(H(user(i),:)*T(:,i)/norm(T(:,i))))^2/(1+(P/M)*temsum);
end