function[sum_rate]=once_sim(H,T,P)
[K,M]=size(H);
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=sum(R);