dengxin.m

来自「ofdm系统中功率自适应分配算法研究的matlab源代码（k＝20）」· M 代码 · 共 30 行

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;