shiyan_u0.m

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

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

    powerpenalty(0.5*kesi+1)=klse822();
end
zz=powerpenalty(1);
for ii=0:4
    powerpenalty(ii+1)=powerpenalty(ii+1)-zz;
end
save PMD powerpenalty;
kesi=0:0.2:0.8;
p=polyfit(kesi,powerpenalty,4);
kesi1=0:0.001:0.8;
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');