function [y] = funcion_heading(u1,u2)
%u1: coordenada x del avion
%u2: coordenada y del avion
%u3: coordenada z del avion
%c1,c2;constantes multiplicadoras
c1=0.01;
c2=0.99;
%r: vector con origen el avion y fin el origen de pista, sin coordenada z
r=[u2 abs(u1)];
%p: vector con origen el avion, fin el punto mas cercano de la linea que
%contiene el eje central de la pista y que es perpendicular a esta
p=[0 abs(u1)];
%hd: vector resultado de la suma ponderada de r y p
hd=c1*r+c2*p;
%calculo del angulo que forma el vector suma con el vector perpendicular
cosang=dot(hd,p)/(norm(hd)*norm(p));
ang=(acos(cosang))*180/pi;
%Salida
y(1)=-ang;
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
