Cap100basetx.m

a=7.26e-7;
b=4.56e-12;
f=[1:1024]/1024*2.5e8;
d4=385;
H4=-8.68*d4*(a*sqrt(f)+b*f);
C4=-42.5+15*log10(f/16e6);
txmask
rcv=logpspc+H4;
crsn=logpspc+C4;
SNR=H4-C4;
figure(1)
plot(f(1:200),rcv(1:200),f(1:200),crsn(1:200))
grid
xlabel('Frequency (Hz)')
ylabel('PSD (dBm/Hz)')
figure(2)
plot(f(1:200),SNR(1:200))
grid
xlabel('Frequency (Hz)')
ylabel('SNR (dB)')
capci=sum(log2(1+10.^(SNR(1:200)/10)))*5e7/200