PLCapcity.m

C1=sum(log(1+((abs(Hf(34:1000))).^2)./((abs(h(34:1000)')).^2)/1.58e5))*30e6/1000/log(2)
C2=sum(log(1+((abs(Hf(34:1000))).^2)./((abs(h(34:1000)')).^2)/1.58e7))*30e6/1000/log(2)
C3=sum(log(1+((abs(Hf(34:1000))).^2)/5.01e-10))*30e6/1000/log(2)
C4=sum(log(1+((abs(Hf(34:1000))).^2)/5.01e-8))*30e6/1000/log(2)
Ctest=966*log(1+1.9952e7)*30e6/1000/log(2)