A general precompiler for algebraic manipulation

A generalized precompiler for systems performing algebraic manipulation of Poisson series has been written. It accepts a trigonometric superset of FORTRAN IV similar to Jefferys' TRIGRUN language (Jefferys, 1972) and generates a valid FORTRAN IV program which drives an abstract formula manipulation machine. This machine is designed to be generally compatible with any manipulation system, and has been implemented with two such systems. The precompiler is written in standard FORTRAN IV and was designed to allow simple conversion for use on most computers.     

对数周期偶极子天线的电特性分析

正在发展中的云南天文台０．５ＧＨｚ—１．５ＧＨｚ快速射电频谱分析仪的接收天线部分是计划在现有的十米口径抛物面天线的焦点上安装两付互相垂直的对数周期偶极子天线作为馈源。本文介绍了设计一个ＦＯＲＴＲＡＮ语言程序对对数周期偶极子天线电特性进行分析的过程，它可以在整个工作频段范围内计算出任意频率点上一副对数周期偶极子天线两个主平面内的方向图、波瓣宽度、增益和输入阻抗，并把结果以数据文件的形式记录下来，结合绘图软件，就可给出直观的结果，从而为对数周期偶极子天线的设计方案提供了可靠的理论判据。本文最后针对为云南天文台抛物面所设计的馈源作具体的分析。     

A precise multipass method for satellite Doppler positioning

A precise and efficient algorithm is developed for determining the locations of radio beacons (e.g. of the ARGOS or COSPAS/SARSAT type) based on Doppler shift measurements in overflying satellites. The method distinguishes itself through: (1) the use of a very compact analytic orbital theory valid for all eccentricities, (2) autonomous orbit improvement based on Doppler data for one or more local reference beacons accessible at a single LUT, (3) simultaneous orbit improvement and calculation of beacon coordinates for an arbitrary number of satellites, satellite passes, and beacons, and (4) very efficient semi-analytic matrix inversion by partitioning into global, semi-global, and local parameters.The algorithm has been implemented in a FORTRAN program which can be run on a PC. Error statistics are presented from applications of the program to a large number of actual Doppler curves obtained with the ARGOS and COSPAS/SARSAT systems.     

CCD整体平差方法的FORTRAN语言归算程序

介绍了实现CCD整体平差的FORTRAN归算程序的详细流程,包括(a)数据收集、(b)预处理和(c)整体平差与结果输出共三个步骤,并给出了各步所调用的子程序名称及其功能说明。采用模拟数据对该程序进行了可靠性检验,并将该程序应用于4颗河外射电源的光学定位。结果表明, 由于利用整体平差方法参加归算的观测资料覆盖天区扩大,包含了较多的参考星,因而可以获得高于单底片归算的定位精度。     

Orbit classification in arbitrary 2D and 3D potentials

A method of classifying generic orbits in arbitrary 2D and 3D potentials is presented. It is based on the concept of spectral dynamics introduced by Binney &38; Spergel that uses the Fourier transform of the time series of each coordinate. The method is tested using a number of potentials previously studied in the literature and is shown to distinguish correctly between regular and irregular orbits, to identify the various families of regular orbits (boxes, loops, tubes, boxlets, etc.), and to recognize the second-rank resonances that bifurcate from them. The method returns the position of the potential centre and, for 2D potentials, the orientation of the principal axes as well, should this be unknown. A further advantage of the method is that it has been encoded in a FORTRAN program that does not require user intervention, except for 'fine tuning' of search parameters that define the numerical limits of the code. The automatic character makes the program suitable for classifying large numbers of orbits.     

A programming system for analytical series expansions on a computer

A programming system is described for the manipulation of Poisson series on a computer. The general structure is described together with the most important individual subroutines. The system exists in two versions, one with 3 angular and 3 polynomial variables and the other with twice as many variables. The programming has been done first in FORTRAN IV, but the most crucial subroutines have been rewritten in machine language for more efficiency. The storage mechanism is such that obsolete series can be erased from the main memory in order to make room for new series to be created.This paper presents the results of one phase of research carried out at the Jet Propulsion Laboratory, California Institute of Technology, under Contract No. NAS 7-100, sponsored by the National Aeronautics and Space Administration.     

The Persian calendar for 3000 years

Using the analytical theory of the motion of the Earth around the Sun the times of the vernal (Spring) equinox has been calculated over the period from the Hijra (AD 622) to AD 3800. These data alone allow to decide whether a particular Persian (or Jalaali) calendar year is common or leap. Presented analysis shows that an algorithm implemented in the so called Khayam program is valid for the year 1799 to 2256 (1178 to 1634 Jalaali). A concise algorithm has been worked out that reconstructs the pattern of leap years over time span of about 3000 years. FORTRAN routines for conversion between the Jalaali, Gregorian and Julian calendars and the Julian Day Number are presented.     

Lamé functions of the first kind generated by computer

The definitions of the Lamé functions and the ellipsoidal coordinates inherent to them are introduced. Proceeding from Lamé's differential equation the four classes of Lamé functions of the first kind are generated by computer with the aid of formula manipulation techniques. For this purpose algebraic expressions for the coefficients in the Lamé polynomials are constructed by virtue of recurrence formulae and presented in tabular as well as machine readable form for further processing. The determination of a parameter appearing in the polynomial coefficients leads to an algebraic eigenvalue problem. It is shown how this eigenvalue problem is amenable to numerical treatment. The polynomial coefficients together with the eigenvalues are characteristic for the Lamé functions and define them completely.     

Automated,closed form integration of formulas in elliptic motion

In some perturbation theories it is possible to avoid expansion of the perturbations in powers of the eccentricity, obtaining results in closed form by using the true or eccentric anomaly instead of the mean anomaly. This paper describes an algorithm (which has been programmed for the 6600 computer using the formula manipulation system TRIGMAN) for automatically performing the integrals which arise in these theories.     

Literal solution for Hill's lunar problem

A computer program for the manipulation of power series in several variables is used to find the first order solution to Hill's lunar problem. The ratiom of the mean motion of the Sun to that of the Moon is kept as a formal parameter. The series inm are known to converge very poorly. It is shown how efficient algorithms among them the Lie transformation allow us to compute the series inm as far as they are needed. When the series are evaluated at Brown's numerical value form the results achieve or exceed his accuracy.     

Recursive calculation of hansen coefficients

Hansen coefficients are used in expansions of the elliptic motion. Three methods for calculating the coefficients are studied: Tisserand's method, the Von Zeipel-Andoyer (VZA) method with explicit representation of the polynomials required to compute the Hansen coefficients, and the VZA method with the values of the polynomials calculated recursively. The VZA method with explicit polynomials is by far the most rapid, but the tabulation of the polynomials only extends to 12th order in powers of the eccentricity, and unless one has access to the polynomials in machine-readable form their entry is laborious and error-prone. The recursive calculation of the VZA polynomials, needed to compute the Hansen coefficients, while slower, is faster than the calculation of the Hansen coefficients by Tisserand's method, up to 10th order in the eccentricity and is still relatively efficient for higher orders. The main advantages of the recursive calculation are the simplicity of the program and one's being able to extend the expansions to any order of the eccentricity with ease. Because FORTRAN does not implement recursive procedures, this paper used C for all of the calculations. The most important conclusion is recursion's genuine usefulness in scientific computing.     

A computer aided analysis of the Sitnikov Problem

We deal with the problem of a zero mass body oscillating perpendicular to a plane in which two heavy bodies of equal mass orbit each other on Keplerian ellipses. The zero mass body intersects the primaries plane at the systems barycenter. This problem is commonly known as theSitnikov Problem. In this work we are looking for a first integral related to the oscillatory motion of the zero mass body. This is done by first expressing the equation of motion by a second order polynomial differential equation using a Chebyshev approximation techniques. Next we search for an autonomous mapping of the canonical variables over one period of the primaries. For that we discretize the time dependent coefficient functions in a certain number of Dirac Delta Functions and we concatenate the elementary mappings related to the single Delta Function Pulses. Finally for the so obtained polynomial mapping we look for an integral also in polynomial form. The invariant curves in the two dimensional phase space of the canonical variables are investigated as function of the primaries eccentricity and their initial phase. In addition we present a detailed analysis of the linearized Sitnikov Problem which is valid for infinitesimally small oscillation amplitudes of the zero mass body. All computations are performed automatically by the FORTRAN program SALOME which has been designed for stability considerations in high energy particle accelerators.     

Rayleigh scattering in planetary atmospheres

The vector equation of radiative transfer is solved both for conservative and non-conservative planetary atmospheres using the method of discrete ordinates. The atmosphere, bounded by a Lambert bottom, is considered plane-parallel and homogeneous. The scattering in the atmosphere obeys the Rayleigh or Rayleigh-Cabannes law. The compiled package of FORTRAN codes allows us to find the Stokes parameters for such an atmosphere at arbitrary optical depth.     

Ground traces of artificial earth satellites

This paper considers the ground trace of an artificial earth satellite. It determines the effects of the trace caused by perturbations due to atmospheric drag, the oblateness of the earth, and the moon and the sun as a third body.The necessary mathematical relations giving these perturbations which are available in literature are utilized (Betz, 1967; Brouwer and Clemence, 1961; Brouwer and Hori, 1961; Danby, 1962; Escobal, 1965; Kentet al., 1963; Kozai, 1962). Those relations unavailable elsewhere are derived.The computation was done by programming in FORTRAN language and utilizing an IBM 360/65.Captain, USAF     

Quenching of auroral hydrogen line emission by collisional ionization

With the aid of the classical impulse approximation a simple formula is derived for the cross section for electron loss from fast excited hydrogen atoms passing through a neutral gas. Detailed computations relevant to the quenching of hydrogen line emission in auroras are performed. It is found, incidentally, that the population distribution amongst the states of a level is unlikely to be brought into statistical equilibrium by the collisions suffered by the hydrogen atoms.     

A method of computing the gravitational field of an axially symmetric flat galaxy

A method has been developed for computing the gravitational force field of an axially symmetric flat galaxy from its surface mass density. The method is simple to program, fast, and accurate. An inversion formula is derived that allows computation of surface density from rotation curves by use of any method that converts density to force. The method is compared with a method of Clutton-Brock that utilizes Hankel transforms of Laguerre functions.     

Rayleigh-Cabannes scattering in planetary atmospheres II. Constant internal sources in nonconservative atmospheres

The vector equation of radiative transfer is solved for non-conservative homogeneous plane-parallel atmosphere using the method of discrete ordinates. The scattering processes in the atmosphere bounded by a Lambert bottom are described by the Rayleigh-Cabannes phase matrix. The primary radiation field is generated by constant internal sources. A package of FORTRAN subroutines is compiled to find the axial radiation field for such an atmosphere at arbitrary optical depth.     

Motion of a satellite in the equatorial plane of a spheroid

An exact, closed-form solution of the problem of the motion of a satellite in the equatorial plane of an oblate body is obtained. It is shown that the classic formula for the motion of the perihelion is a first order approximation to the exact formula.     

A note on an attempt at more efficient Poisson series evaluation

Eliminating many of the trigonometric function calls by a suitable series transformation has resulted in a substantial reduction of on-line computing time while very long Poisson series are being evaluated. A further reduction has been realized by applying a short SNOBOL processor to the FORTRAN coding of the transformed series which eliminates many of the multiplication operations during the course of series evaluation.     

Creation of the Hubble Space Telescope

The Hubble Space Telescope has been the most successful space astronomy project to date, producing images that put the public in awe and images and spectra that have produced many scientific discoveries. It is the natural culmination of a dream envisioned when rocket flight into space was first projected and a goal set for the US space program soon after NASA was created. The design and construction period lasted almost two decades and its operations have already lasted almost as long. The capabilities of the observatory have evolved and expanded with periodic upgrading of its instrumentation, thus realizing the advantages of its unique design. The success of this long-lived observatory is closely tied to the availability of the Space Shuttle and the end of the Shuttle program means that the end of the Hubble program will follow before long.