userselect.m

有限反馈的仿真程序

function[user]=userselect(H,Quantized_H,SINR,epsi)
[K,M]=size(H);
userset_num=zeros(1,M);
user=zeros(1,M);
userset_num(1)=K;
for i=1:K
    if SINR(i)==max(SINR)
        user(1)=i;
    end
end
for i=1:M-1
    num=0;
    for k=1:K
        flag=1;
        for j=1:i
            if abs(conj(Quantized_H(k,:))*Quantized_H(user(j),:).')>epsi
                flag=0;
            end
        end
        if flag==1;
            userset_num(i+1)=userset_num(i+1)+1;
            num=num+1;
            temp1(num)=k;
        end
    end
    if num==0
        num=1;
    end
    temp=zeros(1,num);
    for temp2=1:num
        temp(temp2)=SINR(temp1(temp2));
    end
    for temp2=1:num
        if temp(temp2)==max(temp)
            user(i+1)=temp1(temp2);
        end
    end
end