11
void reflection_coefficient(struct FCOMPLEX zi, double *magnitude, double *phase)
13
struct FCOMPLEX zo, rho;
15
/* reflection coefficient rho = (Zi-Zo)/(Zi+Zo) */
18
rho=Cdiv(Csub(zi, zo), Cadd(zi,zo));
20
*phase=atan2(rho.i,rho.r);
24
fprintf(stderr,"Errno =%d in refco.c \n", errno);