Search NASA⌕ Search

SEARCH · Search NASA

Results for “least squares solution”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Optimum data weighting and error calibration for estimation of gravitational parameters

A new technique was developed for the weighting of data from satellite tracking systems in order to obtain an optimum least-squares solution and an error calibration for the solution parameters. Data sets from optical, electronic, and laser systems on 17 satellites in GEM-T1 Goddard Earth Model-T1 (GEM-T1) were employed toward application of this technique for gravity field parameters. Also GEM-T2 (31 satellites) was recently computed as a direct application of the method and is summarized. The method employs subset solutions of the data associated with the complete solution to agree with their error estimates. With the adjusted weights the process provides for an automatic calibration of the error estimates for the solution parameters. The data weights derived are generally much smaller than corresponding weights obtained from nominal values of observation accuracy or residuals. Independent tests show significant improvement for solutions with optimal weighting. The technique is general and may be applied to orbit parameters, station coordinates, or other parameters than the gravity model.

Lerch, Francis J.↗
NASA NTRS1989

Eclipsing binary stars as tests of gravity theories - The apsidal motion of AS Camelopardalis

AS Camelopardalis is an 8th-magnitude eclipsing binary that consists of two main-sequence (B8 V and a B9.5 V) components in an eccentric orbit (e = 0.17) with an orbital period of 3.43 days. Like the eccentric eclipsing system DI Herculis, and a few other systems, AS Cam is an important test case for studying relativistic apsidal motion. In these systems, the theoretical general relativistic apsidal motion is comparable to that expected from classical effects arising from tidal and rotational deformation of the stellar components. Accurate determinations of the orbital and stellar properties of AS Cam have been made by Hilditch (1972) and Khalliulin and Kozyreva (1983) that permit the theoretical relativistic and classical contributions to the apsidal motion to be determined reasonably well. All the published timings of primary and secondary minima have been gathered and supplemented with eclipse timings from 1899 to 1920 obtained from the Harvard plate collection. Least-squares solutions of the eclipse timings extending over an 80 yr interval yield a smaller than expected apsidal motion, in agreement with that found by Khalliulin and Kozyreva from a smaller set of data. The observed apsidal motion for AS Cam is about one-third that expected from the combined relativistic and classical effects. Thus, AS Cam joins DI Her in having an observed apsidal motion significantly less than that predicted from theory.

Maloney, Frank P.↗
NASA NTRS1989

Automated generation of influence functions for planar crack problems

A numerical procedure for the generation of influence functions for Mode I planar problems is described. The resulting influence functions are in a form for convenient evaluation of stress-intensity factors for complex stress distributions. Crack surface displacements are obtained by a least-squares solution of the Williams eigenfunction expansion for displacements in a cracked body. Discrete values of the influence function, evaluated using the crack surface displacements, are curve fit using an assumed functional form. The assumed functional form includes appropriate limit-behavior terms for very deep and very shallow cracks. Continuous representation of the influence function provides a convenient means for evaluating stress-intensity factors for arbitrary stress distributions by numerical integration. The procedure is demonstrated for an edge-cracked strip and a radially cracked disk. Comparisons with available published results demonstrate the accuracy of the procedure.

Sire, Robert A.↗
NASA NTRS1991

Demonstration of precise estimation of polar motion parameters with the global positioning system: Initial results

Data from the Global Positioning System (GPS) were used to determine precise polar motion estimates. Conservatively calculated formal errors of the GPS least squares solution are approx. 10 cm. The GPS estimates agree with independently determined polar motion values from very long baseline interferometry (VLBI) at the 5 cm level. The data were obtained from a partial constellation of GPS satellites and from a sparse worldwide distribution of ground stations. The accuracy of the GPS estimates should continue to improve as more satellites and ground receivers become operational, and eventually a near real time GPS capability should be available. Because the GPS data are obtained and processed independently from the large radio antennas at the Deep Space Network (DSN), GPS estimation could provide very precise measurements of Earth orientation for calibration of deep space tracking data and could significantly relieve the ever growing burden on the DSN radio telescopes to provide Earth platform calibrations.

Lichten, S. M.↗
NASA NTRS1991

A new activity index for comets

An activity index (AI) is derived from observational data to measure the increase of activity in magnitudes for comets when brightest near perihelion as compared to their inactive reflective brightnesses at great solar distances. In all, AI determinations are made for 32 long-period (L-P) comets and for 14 short-period (S-P). The range of values of AI is the order of 3 to 10 magnitudes with a median about 6. An expected strong correlation with perihelion distance, q, was found to vary as approx. q(-2.3). Residuals from a least-square solution, delta AI, were used for comparing comets of different orbital classes, the standard deviation of a single value of delta AI was only plus or minus 1(m).1 for L-P comets and plus or minus 1(m)2 for S-P comets. Among the L-P comets, 19 of Period, P greater than 10(exp 4) yr yielded delta AI greater than or equal to -0(m).27 compared to +0(m).39 for 13 of 10(exp 4) greater than P greater than 10(exp 2) yr. This denies any effect of aging among the L-P comets. The 14 S-P comets yielded (delta AI) less 0(m).3 plus or minus 0(m).3 than the 32 L-P comets. The results suggest a common actively level, nature and probably origin for comets of all orbits, sizes and ages.

Whipple, Fred L.↗
NASA NTRS1991

The control network of Mars: April 1991

The modern geodetic control network of Mars was first established based on Mariner 9 images with 1-2 km/pixel resolutions and covered almost the entire Martian surface. The introduction of higher resolution (10-200 meter/pixel) Viking Orbiter images greatly improved the accuracy and density of points in the control network. Analysis of the Viking Lander radio tracking data led to more accurate measurements of Mars' rotation period, spin axis direction, and the lander coordinates relative to the inertial reference frame. The prime meridian on Mars was defined by the Geodesy/Cartography Group of the Mariner 9 Television Team as the crater Airy-0, located about 5 degrees south of the equator. The Viking 1 Lander site was identified on a high resolution Viking frame. The control point measurements form the basis of a least squares solution determined by analytical triangulation after the pixel measurements are corrected for geometric distortions and converted to millimeter coordinates in the camera focal plane. Photogrammetric strips encircling Mars at the equator and at 60 degree north south were used to strengthen the overall net and improve the accuracy of the coordinates of points. In addition, photogrammetric strips along 0, 90, 180, and 270 degrees longitude to the Viking 1 Lander site have all significantly strengthened the control network. Most recently, photogrammetric strips were added to the net along 30 degrees north latitude between 0 and 180 degrees, and along 30 degrees between 180 and 360 degrees. The Viking 1 Lander site and Airy-0 are linked through photogrammetric strips occurring along the 0 degree meridian from Airy-0 to 65 degrees north, from that point through the Viking 1 Lander site to the equator, and along the equator to 180 degrees longitude. The Viking 1 lander site is thus a well calibrated area with coordinates of points accurate to approximately 200 meters relative to the J2000 inertial coordinate system. This will be a useful calibration region for upcoming missions. The current status of the control network calculations is presented.

Davies, Merton E.↗
NASA NTRS1992

Estimation of Kalman filter gain from output residuals

This paper presents a procedure for estimating the Kalman filter gain from output residuals. The system state space model is assumed to be known, but the process and noise covariance are unknown. The proposed procedure consists of three basic steps. First, the output residuals are computed from the given model and a given set of input-output data. Second, a linear regression model for this part of the response is computed by a least squares solution. Third, the Kalman filter gain is then estimated from the coefficients of this model. Numerical results using experimental data are presented to illustrate the validity of the developed procedure.

Juang, Jer-Nan↗
NASA NTRS1993

What determines the stellar mass functions in globular clusters?

We analyze the dependence of stellar mass function slopes for a sample of 17 globular clusters on a variety of cluster parameters. The principal novelty of our approach is the use of appropriate multivariate statistical methods to disentangle the complex situation which is present in this problem: the slopes depend simultaneously on more than one variable, and many cluster parameters are mutually correlated. We find that the mass function slopes in the range M/M(solar) = 0.5-0.8 are largely determined by the position in the Galaxy and to a lesser extent by the cluster metallicity. Clusters closer to the Galactic center or plane have shallower mass function slopes. At a given distance to the Galactic center, clusters closer to the Galactic plane have shallower mass function slopes. At a given R(GC) and/or Z(GP), more metal-rich clusters have shallower mass function slopes. Thus, the monovariate correlations with the position or metallicity are both correct, but only partial, and in terms of slopes, biased descriptions of the situation. We present trivariate least-squares solutions where the mass function slopes can be predicted within the measurement accuracy. This relation can serve as a powerful observational constraint for theories of globular cluster formation and evolution, and it is one of the tightest correlations between globular cluster properties now known.

Djorgovski, S.↗
NASA NTRS1993

Estimation of Kalman filter gain from output residuals

This paper presents a procedure to estimate the Kalman filter gain from input-output measurement data with a given system model. The system model can be a finite element model or an experimental model from any identification method. The procedure consists of three basic steps. First, the stochastic portion related to the residuals of the response is computed. Second, the coefficients of a linear difference model for the stochastic portion are estimated by a least-squares solution that minimizes the filter residual. Third, the Kalman filter gain is computed from these model coefficients. Experimental results are presented to illustrate the usefulness of the developed procedure.

Juang, Jer-Nan↗
NASA NTRS1954

Application of several methods for determining transfer functions and frequency response of aircraft from flight data

In the process of analyzing the longitudinal frequency-response characteristics of aircraft, information on some of the methods of analysis has been obtained by the Langley Aeronautical Laboratory of the National Advisory Committee for Aeronautics. In the investigation of these methods, the practical applications and limitations were stressed. In general, the methods considered may be classed as: (1) analysis of sinusoidal response, (2) analysis of transient response as to harmonic content through determination of the Fourier integral by manual or machine methods, and (3) analysis of the transient through the use of least-squares solutions of the coefficients of an assumed equation for either the transient time response or frequency response (sometimes referred to as curve-fitting methods). (author)

Eggleston, John M↗
NASA NTRS1994

A new determination of radii and limb parameters for Pluto and Charon from mutual event lightcurves

Over the past several years Pluto-Charon mutual events have yielded progressively more accurate estimates of Charon's orbital elements and the radii of Pluto and Charon (e.g., Buie, Tholen, and Horne, 1992). Analysis of the 1988 stellar occultation by Pluto indicates a radius for Pluto that is about 4%, or 50 km, larger than the mutual event radius of 1151 km. One possible explanation for the discrepancy is that the mutual event modeling treats Pluto and Charon as uniformly bright disks. If they are limb-darkened, the mutual event fits could underestimate their radii. In this paper we use an independent mutual event data set (Young and Binzel, 1992) to fit for Pluto and Charon's radii in a manner independent of either object's limb profile or albedo distribution. Our least-squares solution indicates that Pluto's radius is 1164 +/- 22.9 km and Charon's radius is 621 +/- 20.6 km.

Young, Eliot F.↗
NASA NTRS1996

A Frequency-Domain Substructure System Identification Algorithm

A new frequency-domain system identification algorithm is presented for system identification of substructures, such as payloads to be flown aboard the Space Shuttle. In the vibration test, all interface degrees of freedom where the substructure is connected to the carrier structure are either subjected to active excitation or are supported by a test stand with the reaction forces measured. The measured frequency-response data is used to obtain a linear, viscous-damped model with all interface-degree of freedom entries included. This model can then be used to validate analytical substructure models. This procedure makes it possible to obtain not only the fixed-interface modal data associated with a Craig-Bampton substructure model, but also the data associated with constraint modes. With this proposed algorithm, multiple-boundary-condition tests are not required, and test-stand dynamics is accounted for without requiring a separate modal test or finite element modeling of the test stand. Numerical simulations are used in examining the algorithm's ability to estimate valid reduced-order structural models. The algorithm's performance when frequency-response data covering narrow and broad frequency bandwidths is used as input is explored. Its performance when noise is added to the frequency-response data and the use of different least squares solution techniques are also examined. The identified reduced-order models are also compared for accuracy with other test-analysis models and a formulation for a Craig-Bampton test-analysis model is also presented.

Blades, Eric L.↗
NASA NTRS2000

Retrieving the Balanced Winds on the Globe as a Generalized Inverse Problem

A generalized inverse technique is applied to retrieve two types of balanced winds that characterize the large-scale dynamics of the atmosphere: rotational winds based upon the linear balance equation, and divergent winds based upon the vorticity budget equation. Both balance equations are singular at or near the equator. The balance equations are transformed in spherical harmonic function space to an under-determined system, for which the scale-weighed least-squares solution consists of a sum of principal and singular components. The principal components represent the response to the source function for the regular eigenmodes, while the singular components are determined by the projection of an independent measurement on the singular eigenmodes. The method was tested with the NCEP/NCAR reanalysis data in which a quasi-balance condition exists. A realistic balanced wind field is retrievable when the singular components are computed based upon the reanalyzed wind data.

Lu, Huei-Iin↗
NASA NTRS2004

Least-Squares Self-Calibration of Imaging Array Data

When arrays are used to collect multiple appropriately-dithered images of the same region of sky, the resulting data set can be calibrated using a least-squares minimization procedure that determines the optimal fit between the data and a model of that data. The model parameters include the desired sky intensities as well as instrument parameters such as pixel-to-pixel gains and offsets. The least-squares solution simultaneously provides the formal error estimates for the model parameters. With a suitable observing strategy, the need for separate calibration observations is reduced or eliminated. We show examples of this calibration technique applied to HST NICMOS observations of the Hubble Deep Fields and simulated SIRTF IRAC observations.

Arendt, R. G.↗
NASA NTRS2005

Regional Deformation Studies with GRACE and GPS

GRACE data indicate large seasonal variations in gravity that have been shown to be to be related to climate-driven fluxes of surface water. Seasonal redistribution of surface mass deforms the Earth, and our previous study using GRACE data demonstrate that annual radial deformations of +/-13 mm in the region of Amazon River Basin were observed by both GRACE and ten GPS sites in the region. For the GRACE determinations, we estimate in a least-squares solution for each Stokes coefficient parameters that represent the amplitudes of the annual variation. We then filter these parameters based on a statistical test that uses the scatter of the postfit residuals. We demonstrate by comparison to the GPS amplitudes that this method is more accurate, for this region, than Gaussian smoothing. Our model for the temporal behavior of the gravity coefficients includes a rate term, and although the time series are noisy, the glacial isostatic adjustment signal over Hudson s Bay can be observed. .

Davis, J. L.↗
NASA NTRS2017

Precision Asteroid Astrometry

Among the methods used to guide spacecraft to their destinations, Optical Navigation (OpNav) remains an effective option. OpNav makes use of star fields and small body ephemerides to precisely locate spacecraft. To facilitate accurate OpNav, the small body ephemerides must be constantly updated; asteroid orbits accumulate errors each year of a few milliarcseconds. Through extended exposures, taken with strategic offsets, a least-squares solution can be found that determines updated ephemeris data. This updated data can also be used by occultation astronomers to gain further information about the small bodies, including their size and shape. Using the 24-inch telescope at the Caltech Table Mountain Observatory (TMO), we capture two or more 180 second exposures of each target. These images, combined with a file for the predicted background star field and two reference files, are then processed through a series of scripts and programs. Starting with a prediction file and two to five exposures of the asteroid, the data is processed. This original data is about 32MB per observation. Once the data are reduced to only Right Ascension and Declination for each target, the data are ready for delivery. This consists of text only, and for each target takes about 80 bytes; this resulting data reduction is about five orders of magnitude. This method produces observed positions that are refined by about 12 milliarcseconds, a refinement that is accomplished almost nowhere else. The occultation observations that are facilitated by the ephemerides being refined also produce results that are not possible in any other way from ground-based observations.

Dial, Jason↗
NASA NTRSDate not supplied

Constrained Field Construction Using Bernstein Polynomial Derived Operators

This method recasts the mass-conservation equation as a Sylvester equation employing high-order derivative operators derived from modified Bernstein polynomial expansions. Given prescribed velocity fields, the algorithm yields a constrained, least-squared solution for the associated density field. To demonstrate the practical utility of this methodology, it is applied to a computationally-derived two-dimensional isolator dataset. This reconstruction method is envisioned to be used in conjunction with diagnostic techniques to aid in the quantification of isolator flow fields, since obtaining highly characterized datasets within this engine component is exceedingly difficult.

Isolator↗
NASA NTRSDate not supplied

Improvements in VLBI Data Analysis for Advanced in Geodesy

I will provide an overview of advances in VLBI data analysis implemented at NASA GSFC. That include the use of B-splines for modeling non-linear station position evolution, the use of all non-tidal mass loadings in data reduction, the use of slant a priori path delays from the output of numerical weather models, the use of B-splines for evaluation of the Earth orientation parameters as a continuous function, the use of elevation-dependent weighting, processing both 24 hr and 1 hr experiments in a single least square solutions and other innovations. I will show how this changes in data analysis affect space geodesy results, such as estimation of station positions, Earth orientation parameters, and source coordinates.

Leonid Petrov↗