cal_sinr_sus.m

有限反馈的仿真程序

function[SINR_sus]=cal_sinr_sus(H,Quantized_H,P,epsi)
[K,M]=size(H);
rou=P/M;
cosfai=sqrt((1+epsi)*(1-(M-1)*epsi)/(1-(M-2)*epsi));
fai=acos(cosfai);
for i=1:K
    sita(i)=acos(abs(conj(H(i,:))*Quantized_H(i,:).'/norm(H(i,:))));
end
for i=1:K
    SINR_sus(i)=rou*(norm(H(i,:)))^2*cos(sita(i)+fai)^2/(1+rou*(norm(H(i,:)))^2*sin(sita(i))^2);
end