function cr = reflexion_horizontal(e,d,angulo_grados) %Función que calcula el coeficiente de reflexión para polarización %horizontal ang=angulo_grados*(pi/180); f=1*10^9; c=3*10^8; numerador = sin(ang) - sqrt(e - cos(ang)*cos(ang)); denominador = sin(ang) + sqrt(e - cos(ang)*cos(ang)); rh = numerador/denominador; numerador = 2*pi*f*sqrt(e - sin(ang)); k=numerador/c; numerador=1-exp(-i*2*k*d); denominador=1-(rh^2)*exp(-i*2*k*d); cr=rh*(numerador/denominador); end