www.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;