ZL=100+1j*80;Z0=75;RL=real(ZL);XL=imag(ZL);Y0=1/ZL

 

t1=(XL+((RL*((Z0-RL)^2+XL^2))/Z0)^.5)/(RL-Z0)

t2=(XL-((RL*((Z0-RL)^2+XL^2))/Z0)^.5)/(RL-Z0)

 

% when RL=Z0

D1_TL=1/(2*pi)*atan(t1);if atan(t1)<0 D1_TL=D1_TL+.5; end

D2_TL=1/(2*pi)*atan(t2);if atan(t2)<0 D2_TL=D2_TL+.5; end

 

B1=(RL^2*t1-(Z0-XL*t1)*(XL+Z0*t1))/(Z0*(RL^2+(XL+Z0*t1)^2))

B2=(RL^2*t2-(Z0-XL*t2)*(XL+Z0*t2))/(Z0*(RL^2+(XL+Z0*t2)^2))

 

Y0=1/Z0;

Dstub1=1/(2*pi)*atan(B1/Y0);if atan(Y0/B1)<0 Dstub1=Dstub1+.5; end

Dstub2=1/(2*pi)*atan(B2/Y0);if atan(Y0/B2)<0 Dstub2=Dstub2+.5; end