max_ratio.m

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

function [range]=Max_ratio(Ap,M1)
%Ap=39;              %for NRZ-OOK
%M1=1;               %for 2-OOK
Ndb=1;              %power penalty=1dB
range=zeros(0,50);
for ii=0:50
    ln_P_out=-3-0.1*ii;
    P_out=10.^(ln_P_out);
    range(ii+1)=4*M1*sqrt(Ndb)./sqrt( pi*Ap*log(1/P_out) );
end
%P_out=-3:0.1:-8;
%plot(P_out,range);