View Single Post
  #4  
Old 17-03-2014, 04:04 AM
pdalek's Avatar
pdalek (Patrick)
Registered User

pdalek is offline
 
Join Date: Mar 2013
Location: Melbourne, Australia
Posts: 112
Using the ancient notation
d : declination of star, -south +north
L : latitude of observer, -south +north
LHA : local hour angle of star, -east +west
Hc : altitude
Zn : azimuth

Hc = asin( sin(d) sin(L) + cos(d) cos(L) cos(LHA) )
Z = acos( (sin(d) - sin(L) sin(Hc)) / (cos(L) cos(Hc)) )
If sin(LHA) < 0 then Zn = Z else Zn = 360 - Z

Other ways to write formula but his gets signs right without having to think.

Images 1-3 show azimuth for three objects at differing declinations as observed from Melbourne

Max and min Zn (i.e. direction changes) when partial derivative of Zn wrt LHA is zero.
Messy to work out closed form solution. Easier numerically.

Image 4 shows maximum azimuth vs declination for Melbourne observer. Minimum = 360 - maximum.

Generating a plot or table of values takes four lines of code using Mathematica.
Attached Thumbnails
Click for full-size image (b oct.jpg)
4.9 KB15 views
Click for full-size image (can.jpg)
5.1 KB20 views
Click for full-size image (w cen.jpg)
5.1 KB9 views
Click for full-size image (max.jpg)
6.8 KB10 views
Reply With Quote