bopp.m

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

clear;
clc;
close all;
powerpenalty=zeros(1,7);
%DD=zeros(1,6);
for kesi=1:1:7
    deta=2.85+0.3*0.2*kesi;
    save BER1 deta;

    powerpenalty(kesi)=klse822();
end
zz=min(powerpenalty);
for ii=1:7
    powerpenalty(ii)=powerpenalty(ii)-zz;
end
save CDBo_pi2 powerpenalty;
kesi=1.0:0.3:2.8;
p=polyfit(kesi,powerpenalty,7);
kesi1=1:0.01:2.8;
plot(kesi,powerpenalty,'r*',kesi1,polyval(p,kesi1),'b--');
grid;
title('Bo引起的功率损失');

xlabel('Bo (Rb)');
ylabel('Power Penalty');