gusucode.com > m4k20功率分配图 mimo功率优化,OFDM典型仿真程序，贪婪算法仿真程序 > code/m4k20power/m4k20功率分配图xin/dengxin.m

    function y=dengxin(nn,K)




M=4;
a=2*(sqrt(M)-1)/(sqrt(M)*log2(sqrt(M)));
b=3/(M-1);
Pmax=20;
K=20;
sum2=0;
sum1=0;
for k=1:K
    sum1=sum1+1/nn(k);
    sum2=sum2+1/nn(k)*log(1/nn(k));
end
for k=1:K
  %  if log(nn(k))*sum1+sum2<-(b*Pmax/2)
  %      flag=1;
   %     p(k)=0;
    %    sum1=sum1-1/nn(k);
     %   sum2=sum2-1/nn(k)*log(1/nn(k));
     %else flag=0;
     %end
     %   if flag==0;
        p(k)=Pmax/(sum1.*nn(k));
        %  end
    end
sum(p);
y=Pmax/sum(p).*p;