pmdpp.m

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

clear;
clc;
close all;
powerpenalty=zeros(1,4);
%DD=zeros(1,6);
for kesi=0:2:6
    deta=0.5*0.1*kesi;
    save BER1 deta;

    powerpenalty(0.5*kesi+1)=klse822();
end
zz=powerpenalty(1);
for ii=0:3
    powerpenalty(ii+1)=powerpenalty(ii+1)-zz;
end
save PMD_0 powerpenalty;
kesi=0:0.2:0.6;
p=polyfit(kesi,powerpenalty,3);
kesi1=0:0.001:0.6;
plot(kesi,powerpenalty,'r*',kesi1,polyval(p,kesi1),'b--');
grid;
title('NRZ-2DPSK PMD引起的功率损失');

xlabel('V^2*L*D[10^4(GHz)^2]');
ylabel('Power Penalty');