Koskie-Gajic算法

function [pnextrcount]=koskie_gajic1
sir_tar=9;%目标信干比
k=5000;%bi/2ci
w=5e+6;%扩频带宽
eps=1e-15;%迭代误差
v=1e+4;%传输速率
q=5e-15;%背景噪声
n=1:1:50;
d=301:14:1000;%随机50个用户离基站的距离

for i=1:50
    dd=d（1:i）;
    h=0.0075*dd.^（-3.6）;
    power0=5e-15*ones（1i）;
    ponce=power0;
    
j=sum（ponce.*h）-ponce.*h+q;
r=w/v*ponce.*h./j;
pnext=sir_tar*（ponce./r）-k*（ponce./r）.^2;

while sum（abs（pnext-ponce）>=eps）>0
    ponce=pnext;
    j=sum（ponce.*h）-ponce.*h+q;
    r=w/v*ponce.*h./j;
    pnext=sir_tar*（ponce./r）-k*（ponce./r）.^2;
    [a]=find（pnext<0）