exercise 5.14

Z1=100;Z2=150;ZL=225

Z3=ZL

g1=(Z2-Z1)/(Z2+Z1)

g3=(Z3-Z2)/(Z3+Z2)

g2=-g1

% theta=beta*L=2*pi*D

% where D is fraction of wavelength D between [0 1]

 

dD=.01;D=[0:dD:1]

g_exact=(g1+g3*exp(-1j*2*2*pi*D))./(1+g1*g3*exp(-1j*2*2*pi*D))

g_approx=g1+g3*exp(-j*2*2*pi*D)

err=abs(abs(g_exact)-abs(g_approx))

plot(D,err),grid on;

 

 

 

 

 

 

 

 

 

​

​

​

​

nmax=find(err==max(err))

max(unique(err(nmax)))