%File: TRx_Calibrations.m
%Version: 1
%Date: 5 August 2018
%Author: David Ives
%email: di231@cam.ac.uk
%License: 



%UoC OSA_NSR = 0.9822 * TRx_NSR - 0.0081

OSA_NSR_UoC=[0.093428634;
0.08400981;
0.075592725;
0.068381717;
0.056298912;
0.050699071;
0.041720548;
0.031358841;
0.022778231;
0.014304741;
0.00642525;];

TRx_NSR_UoC=[0.102886881;
0.093549183;
0.085104005;
0.077739188;
0.067993907;
0.059274787;
0.050395296;
0.040144865;
0.031287438;
0.02274207;
0.014742801;];




%UCL OSA_NSR = 0.9457 * TRx_NSR - 0.0033

OSA_NSR_UCL=[0.020249046;
0.015990608;
0.01265145;
0.010081115;
0.007928848;
0.006324587;
0.005068278;
0.004064168;
0.003279379;
0.002670337;];

TRx_NSR_UCL=[0.025000454;
0.020382214;
0.01687765;
0.014011639;
0.012044611;
0.010220579;
0.008718781;
0.007849269;
0.006917147;
0.00649907;];



% UoB OSA_NSR = -3.709 * TRx_NSR^2 + 0.9837 * TRx_NSR - 0.0111

OSA_NSR_UoB=[3.14688E-05;
9.9513E-05;
0.000314688;
0.00099513;
0.001577174;
0.002499652;
0.003961682;
0.006278843;
0.009951295;
0.01577174;
0.024996523;
0.03961682;
2.66974E-05;
8.44247E-05;
0.000266974;
0.000844247;
0.001338042;
0.002120653;
0.003361009;
0.005326839;
0.008442472;
0.013380416;
0.02120653;
0.033610085;
2.46763E-05;
7.80332E-05;
0.000246763;
0.000780332;
0.001236743;
0.001960106;
0.003106559;
0.004923564;
0.007803323;
0.012367434;
0.019601062;
0.03106559;
2.93656E-05;
9.28623E-05;
0.000293656;
0.000928623;
0.001471769;
0.002332596;
0.003696916;
0.005859217;
0.009286234;
0.014717689;
0.023325965;
0.036969162;
2.52168E-05;
7.97425E-05;
0.000252168;
0.000797425;
0.001263833;
0.00200304;
0.003174604;
0.005031409;
0.007974246;
0.012638328;
0.0200304;
0.031746044;
2.6554E-05;
8.39711E-05;
0.00026554;
0.000839711;
0.001330852;
0.002109258;
0.003342948;
0.005298215;
0.008397106;
0.013308515;
0.021092575;
0.033429479;
2.53518E-05;
8.01696E-05;
0.000253518;
0.000801696;
0.001270602;
0.002013769;
0.003191608;
0.005058359;
0.008016958;
0.012706022;
0.020137688;
0.031916085;
2.84859E-05;
9.00803E-05;
0.000284859;
0.000900803;
0.001427676;
0.002262714;
0.00358616;
0.00568368;
0.009008026;
0.014276758;
0.022627137;
0.028485878;
0.035861596;];

TRx_NSR_UoB=[0.012271768;
0.012292735;
0.012484154;
0.013168517;
0.013739323;
0.014773978;
0.016305489;
0.019050205;
0.023446568;
0.030208269;
0.043896344;
0.068368617;
0.012407286;
0.012426797;
0.012570106;
0.013100158;
0.013473766;
0.013965476;
0.015369227;
0.017444958;
0.021012827;
0.027011181;
0.03784881;
0.058533134;
0.011898275;
0.01190692;
0.012094772;
0.012736974;
0.013332639;
0.013861674;
0.015252632;
0.017007519;
0.020364587;
0.026059606;
0.035438876;
0.053329952;
0.012243437;
0.012347516;
0.012477865;
0.013414363;
0.013894588;
0.014879415;
0.016459814;
0.018974807;
0.023095615;
0.030294042;
0.042675508;
0.064562738;
0.011790987;
0.011941072;
0.012031681;
0.012635072;
0.013067494;
0.013884761;
0.015172325;
0.017271229;
0.020813285;
0.026600858;
0.036361727;
0.055120092;
0.011949506;
0.01200743;
0.012299672;
0.012854881;
0.013418372;
0.01420458;
0.015536235;
0.017611658;
0.021128601;
0.027329805;
0.038125848;
0.058766728;
0.011705805;
0.011790987;
0.012125499;
0.012880752;
0.013675711;
0.014284203;
0.015625312;
0.017907584;
0.021606541;
0.027742569;
0.03882385;
0.05958288;
0.011081896;
0.011153904;
0.011360984;
0.012031681;
0.012539846;
0.013369785;
0.014845463;
0.017423642;
0.021503764;
0.028579519;
0.040999813;
0.05115108;
0.055475806;];

colors=[1 0 0 ; 0 0.5 0 ; 0 0 1; ];
set(groot,'defaultAxesColorOrder',colors);  


f1=figure(1);
f1.PaperSize=[10,7];
f1.PaperPosition=[0 0 10 7]; 
plot(TRx_NSR_UoC,OSA_NSR_UoC,'o',TRx_NSR_UoB,OSA_NSR_UoB,'o',TRx_NSR_UCL,OSA_NSR_UCL,'o')
ylim([-0.01 0.08])
xlim([0 0.08])
hold on
TRx1=(0.014:0.002:0.09);
TRx2=(0.006:0.002:0.026);
TRx3=(0.012:0.002:0.07);

OSA1=0.9822*TRx1 - 0.0081;
OSA2=0.9457*TRx2 - 0.0033;
OSA3=-3.709*TRx3.^2 +0.9837*TRx3 - 0.0111;

plot(TRx1,OSA1,'-',TRx3,OSA3,'-',TRx2,OSA2,'-')
hold off

xlabel('TRx constellation NSR')
ylabel('Optical equivalent AWGN NSR')
legend({'UoC: y = 0.9822 x - 0.0081', 'UoB: y = -3.709 x^2 + 0.9837 x - 0.0111', 'UCL: y = 0.9457 x - 0.0033'}, 'location', 'NorthWest','fontsize',7)

%saveas(gcf,'TRxCal','pdf')