%Values read from network analyzer (ohms)
Z175in = 25-j*272.8;
Z2in=4.2-j*102.6;
Z24in=2.8379+j*57.205;
Z26in=4.6289+j*84.121;
Z3in =186+j*363.6;

%rotate values as needed for TRL calibration adjustment (values below line
%on smith chart)
gam175in=(Z175in-50)/(Z175in+50);	% (Zin-Zo)/(Zin+Zo)
gam2in=(Z2in-50)/(Z2in+50);
gam175in=gam175in*exp(2*j*pi/2);
gam2in=gam2in*exp(2*j*pi/2);
Z175in=50*(1+gam175in)/(1-gam175in);
Z2in=50*(1+gam2in)/(1-gam2in);

% total length of line - Thru_length/2 (mils)
L=1040

%wavelengths at each frequency (mils)
lam175=4235.704;
lam2 =3705.246;
lam24=3086.292;
lam26=2848.206;
lam3 =2467.231;

%reflection coefficients at input at each frequency Zo=50 ohms
gam175in=(Z175in-50)/(Z175in+50);	% (Zin-Zo)/(Zin+Zo)
gam2in=(Z2in-50)/(Z2in+50);
gam24in=(Z24in-50)/(Z24in+50);
gam26in=(Z26in-50)/(Z26in+50);
gam3in=(Z3in-50)/(Z3in+50);

%reflection coefficients at the load at each frequency
gam175L=gam175in*exp(2*j*L*2*pi/lam175);
gam2L=gam2in*exp(2*j*L*2*pi/lam2);
gam24L=gam24in*exp(2*j*L*2*pi/lam24);
gam26L=gam26in*exp(2*j*L*2*pi/lam26);
gam3L=gam3in*exp(2*j*L*2*pi/lam3);

%impedance at load (ohms)
Z175L=50*(1+gam175L)/(1-gam175L);
Z2L=50*(1+gam2L)/(1-gam2L);
Z24L=50*(1+gam24L)/(1-gam24L);
Z26L=50*(1+gam26L)/(1-gam26L);
Z3L=50*(1+gam3L)/(1-gam3L);

%transmission coefficient at the load
T175L=1+gam175L;
T2L=1+gam2L;
T24L=1+gam24L;
T26L=1+gam26L;
T3L=1+gam3L;

r2d=180/pi;      % conversion for radians to degrees

magr = [abs(gam175L); abs(gam2L); abs(gam24L); abs(gam26L); abs(gam3L)];
angr = [angle(gam175L)*r2d; angle(gam2L)*r2d; angle(gam24L)*r2d; angle(gam26L)*r2d; angle(gam3L)*r2d];
magt = [abs(T175L); abs(T2L); abs(T24L); abs(T26L); abs(T3L)];
angt = [angle(T175L)*r2d; angle(T2L)*r2d; angle(T24L)*r2d; angle(T26L)*r2d; angle(T3L)*r2d];

%Mag[S11] Ang[S11(deg)] Mag[S21] Ang[S21(deg)] Mag[S12] Ang[S12(deg)] Mag[S22] ]
d = [magr angr magt angt zeros(size(magr)) zeros(size(magr)) ones(size(magr))]

% The formatting is not correct here.  Cut and paste the d table into a notebook editor.
%You need to add a column of frequencies at the front and a column of zeros (AngS22) at the end:

%frequencies are: 1.75, 2.0, 2.4, 2.6, 3.0