Equinoxes and Solstices
EQNSOL computes the date of equinoxes and solstices for the input year. This is based upon the formulas given by Jean Meeus in his "Astronomical Algorithms". The date is computed in Julian Date. This function take in JD2CAL function of "Geodetic Toolbox" by Mike Craymer (ID:15285) to convert Julian Date into gregorian calendar date.
Syntax: Y=EQNSOL(YEAR,VERBOSE)
Inputs:
YEAR (default=current Year).
VERBOSE (default=1)
Outputs:
- Y=datenum formatted array.
EQ=EQNSOL(...) will stores dates in Matlab format
Example:
y=eqnsol(2007)
All dates and times are referred to GMT without DST
Spring Equinox 21-Mar-2007 00:08:28
Summer Solstice 21-Jun-2007 18:07:12
Autumn Equinox 23-Sep-2007 09:52:05
Winter Solstice 22-Dec-2007 06:08:51
Created by Giuseppe Cardillo
giuseppe.cardillo-edta@poste.it
To cite this file, this would be an appropriate format: Cardillo G. (2007) Equinoxes and Solstices: compute the date and time of equinoxes and solstices. http://www.mathworks.com/matlabcentral/fileexchange/17977
Cite As
Giuseppe Cardillo (2021). Equinoxes and Solstices (https://github.com/dnafinder/eqnsol), GitHub. Retrieved .
Comments and Ratings (1)
MATLAB Release Compatibility
Platform Compatibility
Windows macOS LinuxCategories
Tags
Acknowledgements
Inspired by: Geodetic Toolbox
Inspired: Aggregate timeseries data to solstices and equinoxes.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!Discover Live Editor
Create scripts with code, output, and formatted text in a single executable document.
Select a Web Site
Choose a web site to get translated content where available and see local events and offers. Based on your location, we recommend that you select: .
Select web siteYou can also select a web site from the following list:
Americas
- América Latina (Español)
- Canada (English)
- United States (English)
Europe
- Belgium (English)
- Denmark (English)
- Deutschland (Deutsch)
- España (Español)
- Finland (English)
- France (Français)
- Ireland (English)
- Italia (Italiano)
- Luxembourg (English)
- Netherlands (English)
- Norway (English)
- Österreich (Deutsch)
- Portugal (English)
- Sweden (English)
- Switzerland
- United Kingdom (English)
Great function! While the level of accuracy does not affect my needs, you may want to let others know if they are in need to to the minute or second, may want to validate the results. Comparing this functions results to the US Naval Observatory, the results differ by approximately 1-minute. I'm guessing this is an artifact of matlab and complier precision. Anything dealing with astronomy and floating point precision becomes significant.
http://aa.usno.navy.mil/data/docs/EarthSeasons.php