convex.m

来自「fading 在无线中的应用」· M 代码 · 共 16 行

function y=convex(g1,Pmax,A,N);
[l,n]=size(g1);
g2=zeros(1,n);
gm=gx(Pmax,A,N);
for i=1:n
   if g1(i)<gm(1)
      g2(i)=Pmax/N(2)-g1(i)*(Pmax/N(2)-gm(2))/gm(1);
   else
      g2(i)=gm(2)/(gm(1)-Pmax/N(1))*(g1(i) -Pmax/N(1));
   end
   g3=Pmax/N(2)-N(1)/N(2)*g1(i);
   g2(i)=max(g3,max(0,g2(i)));
 end;
y=g2;