Code covered by the BSD License  

Highlights from
Oblique Spherical Triangle toolbox

Be the first to rate this file! 6 Downloads (last 30 days) File Size: 70.16 KB File ID: #23075
image thumbnail

Oblique Spherical Triangle toolbox

by Rody Oldenhuis

 

23 Feb 2009 (Updated 14 Sep 2009)

solves the general oblique spherical triangle

| Watch this File

File Information
Description

This toolbox can solve any of the 6 possible subproblems associated with an oblique spherical triangle, when only 3 of the 6 angles are known.

The toolbox basically is an implementation of the set of tools developed in [Wertz, 2001], which aimed to fully automize the procedure and do away with the need of user intervention.

EXAMPLE:

    [b1, c1, C1, b2, c2, C2] = aas(A, B, a)

gives both solutions for the sides b and c and the angle C, when given angles A, B and side a. This particular problem can be called the Angle-Angle-Side problem, hence the name.

Both degrees and radians are implemented, the difference is indicating by appending a 'd' to the function name. The above example expects and returns radians, whereas

    [b1, c1, C1, b2, c2, C2] = aasd(A, B, a)

expects and returns degrees.

[Wertz, 2001]
James R. Wertz, Mission Geometry: Orbit and Constellation Design and Management, 2001. Published by Microcosm and Kluwer Academic Publishers.

MATLAB release MATLAB 7.7 (R2008b)
Tags for This File  
Everyone's Tags
Tags I've Applied
Add New Tags Please login to tag files.
Comments and Ratings (2)
14 Jun 2009 Kamal Abdali

This comment is on the descriptive doc, not on the code. Specifically, it is about the auxiliary functions mal.m, mald.m, msl.m, and msld.m.

First, I think a correction is needed. In the equations for "cos c". the numerator has a plus sign which should be minus.

Second, a suggestion for a very minor improvement. The definitions of "sin c", "cos c", "sin C", and "cos C" all have the same denominator:

  1 - sin a sin b sin A sin B

This is equal to

 1 - (sin b sin A)^2

which contains only the given data b and A, not B which is a computed quantity (in ssa, the main user of the "middle side/angle laws"). So this expression has the advantage of avoiding any errors in the intermediate computation.

14 Sep 2009 Rody Oldenhuis

Whoah! I totally overlooked your comment! I actually found this comment by googling for one of my other submissions :s

Anyway, thanks for the comment. Actually, in MAL/MSL, that whole denominator is unneccesary, since the objective is to take the ATAN2. I posted an update which does just that. I also got rid of the whole MAL/MSL thing, and all dependencies I created back in a time I did not yet know DEAL() ^_^

Thanks for the feedback.

Please login to add a comment or rating.
Updates
14 Sep 2009

updated all files -- no dependencies should exist anymore.
Also got rid of H() and mal()/msl(); they were really just clutter.

Tag Activity for this File
Tag Applied By Date/Time
aerospace Cristina McIntire 23 Feb 2009 11:55:56
aeronautics Cristina McIntire 23 Feb 2009 11:55:56
spherical Cristina McIntire 23 Feb 2009 11:55:56
trigonometry Cristina McIntire 23 Feb 2009 11:55:56
spherical trig Rody Oldenhuis 14 Sep 2009 09:55:54
oblique spherical triangle Rody Oldenhuis 14 Sep 2009 09:55:54
trigonometry Rody Oldenhuis 14 Sep 2009 09:55:54
geometry Rody Oldenhuis 14 Sep 2009 09:55:54

Contact us at files@mathworks.com