Upload New File

Code/inverseGC.m

+%**************************************************************************
+% * Programmiersprache MATLAB
+% * Dateiname: inverseGC.m
+%
+% * Kurs- und Distanzberechnung auf einem Grokreis
+%  
+%
+% ******************************** Inputs *********************************
+% * lat1: Breitengrad des Startpunktes []
+% * lon1: Lngengrad des Startpunktes []
+% * lat2: Breitengrad des Endpunktes []
+% * lon2: Lngengrad des Endpunktes []
+% 
+%
+% ******************************** Outputs ********************************
+% * s_inv: Grokreisdistanz zwischen Start- und Endpunkt [m]
+% * a12:   Startkurs (rechtsweisender Kurs gegenber geograph. Nordpol, 
+%          Uhrzeigersinn pos.) []
+% 
+%
+% * Einrichtung/Autoren: DLR-LY/F.Linke/B.Lhrs
+% * Bearbeitungszustand: freigebenen
+% * Version: 1.0
+% * Datum: 10032013
+% * Letzte nderung: Kommentare eingefgt
+%**************************************************************************
+%**************************************************************************
+%****************************Copyright (C) 2011****************************
+%************DLR Deutsches Zentrum fuer Luft- und Raumfahrt e.V.***********
+%************************German Aerospace Center e.V. *********************
+%********************Institut fuer Lufttransportsysteme/*******************
+%******************Institute of Air Transportation Systems*****************
+%***************Tel. +49 531 295 3801, Fax: +49 531 295 2979***************
+%
+% www.dlr.de/ly
+% www.tu-harburg.de/ilt
+%
+%*  These coded instructions, statements, and computer programs contain****
+%*  unpublished proprietary information of the German Aerospace Center*****
+%*  e.V. and are protected by the copyright law. They may not be***********
+%*  disclosed to third parties or copied or duplicated in any form,********
+%*  in whole or in part, without the prior written consent of the German***
+%*  Aerospace Center e.V.**************************************************
+%**************************************************************************
+
+function [s_inv,a12]=inverseGC(lat1, lon1, lat2, lon2)
+% Konvertierung von Grad zu Rad
+lat1=lat1*pi/180;
+lon1=-lon1*pi/180;
+lat2=lat2*pi/180;
+lon2=-lon2*pi/180;
+
+% Winkel zwischen Start- und Endpunkt
+dist = 2*asin(sqrt((sin((lat1-lat2)/2))^2+...
+    cos(lat1)*cos(lat2)*(sin((lon1-lon2)/2))^2));
+
+% Umrechnung in Meter
+s_inv=dist*180/pi*60./0.0005399568;
+
+if (dist < 1E-16) % Wenn Start und Endpunkt (numerisch) identisch
+  brg = 0;        % dann Kurs = 0
+  
+elseif (abs(cos(lat1)) < 1E-16)   % Wenn Startpunkt (numerisch) am Pol 
+  brg = (3 - sin(lat1)) * pi / 2; % liegt, kann man nur Richtung 
+                                  % Sden/Norden
+  
+else  %sonst
+  brg = modcrs(atan2(sin(lon1-lon2)*cos(lat2),cos(lat1)*sin(lat2)...
+      -sin(lat1)*cos(lat2)*cos(lon1-lon2)));
+end
+
+a12=brg*180/pi;  % Umrechnung in 
+
+
+% crs in Weterbereich 0<crs<=2*pi konvertieren
+function course=modcrs(crs)
+twopi = 2 * pi;
+course = twopi - mod(twopi - crs,twopi);