phase_errorpp.m

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

clear;
clc;
close all;
powerpenalty=zeros(1,5);
for kesi=0:10:40
    detaPHI=0.76*pi*kesi/180;              

    save BER1 detaPHI;

    powerpenalty(0.1*kesi+1)=klse822();
end
zz=powerpenalty(1);
for ii=0:4
    powerpenalty(ii+1)=powerpenalty(ii+1)-zz;
end
save Phase_error powerpenalty;
kesi=0:10:40;
p=polyfit(kesi,powerpenalty,3);
kesi1=0:0.1:50;
plot(kesi,powerpenalty,'r*',kesi1,polyval(p,kesi1),'b-');
grid;
title('NRZ-DPSK phase-error 引起的功率损失');

xlabel('phase error (deg)');
ylabel('Power Penalty');