123.m

来自「一种新的光纤通信调制技术」· M 代码 · 共 11 行

Ap=39;              %for NRZ-OOK
M1=1;               %for 2-OOK
Ndb=1;              %power penalty=1dB
T=1e10;
ln_P_out1=-3:-0.1:-8;
P_out=10.^(ln_P_out1);
range_ratio1=( 400*4*M1*sqrt(Ndb)./sqrt( pi*Ap*log(1./P_out) ) ).^2;
ln_P_out2=-3:-1:-8;
P_out=10.^(ln_P_out2);
range_ratio2=( 400*4*M1*sqrt(Ndb)./sqrt( pi*Ap*log(1./P_out) ) ).^2;
plot(ln_P_out1,range_ratio1,'b-',ln_P_out2,range_ratio2,'ro');