Eureka

Are you sure that Z(DUT)=30+j40 [ohm]? Can you tell by just observing the waveforms?

Let’s check with a python program.

# given waveforms
t=arange(0,1000)
v1=1.0*sin(t*(2.0*np.pi/444.0)+0.0)
v2=sqrt(5.0)/2.0*sin(t*(2.0*np.pi/444.0)+math.atan(1.0/2.0))
# Note that sqrt(5.0)/2.0=1.118 and atan(1.0/2.0)=26.56deg.

# start computing
z0=50.0+0.0j
cv1=1.0+0.0j
cv2=vratio*(cos(phase2advance)+sin(phase2advance)*1.0j)
cvr=2.0*cv1-cv2
cz=z0*(cv2/cvr)
gamma=(cz-z0)/(cz+z0)
swr=(1+abs(gamma))/(1-abs(gamma))

print "vratio=", vratio
print "phase2advance=", phase2advance*360.0/(2.0*np.pi)
print "Z(DUT)=", cz
print "gamma=",gamma
print "swr=", swr

And the results are:

vratio= 1.11803398875
phase2advance= 26.5650511771
Z(DUT)= (30+40j)
gamma= (-1.02119201983e-13+0.5j)
swr= 3.0

So you are correct!

Actually, the DUT was 30 [ohm] + 909.5 [nH], which is at 7026kHz 30 [ohm] + j40 [ohm].