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 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↗
NASA NTRS1993

Least-squares finite element solutions for three-dimensional backward-facing step flow

Comprehensive numerical solutions of the steady state incompressible viscous flow over a three-dimensional backward-facing step up to Re equals 800 are presented. The results are obtained by the least-squares finite element method (LSFEM) which is based on the velocity-pressure-vorticity formulation. The computed model is of the same size as that of Armaly's experiment. Three-dimensional phenomena are observed even at low Reynolds number. The calculated values of the primary reattachment length are in good agreement with experimental results.

Jiang, Bo-Nan↗
NASA NTRS1992

Least-squares finite element solution of 3D incompressible Navier-Stokes problems

Although significant progress has been made in the finite element solution of incompressible viscous flow problems. Development of more efficient methods is still needed before large-scale computation of 3D problems becomes feasible. This paper presents such a development. The most popular finite element method for the solution of incompressible Navier-Stokes equations is the classic Galerkin mixed method based on the velocity-pressure formulation. The mixed method requires the use of different elements to interpolate the velocity and the pressure in order to satisfy the Ladyzhenskaya-Babuska-Brezzi (LBB) condition for the existence of the solution. On the other hand, due to the lack of symmetry and positive definiteness of the linear equations arising from the mixed method, iterative methods for the solution of linear systems have been hard to come by. Therefore, direct Gaussian elimination has been considered the only viable method for solving the systems. But, for three-dimensional problems, the computer resources required by a direct method become prohibitively large. In order to overcome these difficulties, a least-squares finite element method (LSFEM) has been developed. This method is based on the first-order velocity-pressure-vorticity formulation. In this paper the LSFEM is extended for the solution of three-dimensional incompressible Navier-Stokes equations written in the following first-order quasi-linear velocity-pressure-vorticity formulation.

Jiang, Bo-Nan↗
NASA NTRS1989

On damping mechanisms in beams

A partial differential equation model of a cantilevered beam with a tip mass at its free end is used to study damping in a composite. Four separate damping mechanisms consisting of air damping, strain rate damping, spatial hysteresis and time hysteresis are considered experimentally. Dynamic tests were performed to produce time histories. The time history data is then used along with an approximate model to form a sequence of least squares problems. The solution of the least squares problem yields the estimated damping coefficients. The resulting experimentally determined analytical model is compared with the time histories via numerical simulation of the dynamic response. The procedure suggested here is compared with a standard modal damping ratio model commonly used in experimental modal analysis.

Banks, H. T.↗
NASA NTRS1991

On damping mechanisms in beams

A partial differential equation model of a cantilevered beam with a tip mass at its free end is used to study damping in a composite. Four separate damping mechanisms consisting of air damping, strain rate damping, spatial hysteresis and time hysteresis are considered experimentally. Dynamic tests were performed to produce time histories. The time history data is then used along with an approximate model to form a sequence of least squares problems. The solution of the least squares problem yields the estimated damping coefficients. The resulting experimentally determined analytical model is compared with the time histories via numerical simulation of the dynamic response. The procedure suggested here is compared with a standard modal damping ratio model commonly used in experimental modal analysis.

Banks, H. T.↗
NASA NTRS1997

Least-Squares Spectral Element Solutions to the CAA Workshop Benchmark Problems

This paper presents computed results for some of the CAA benchmark problems via the acoustic solver developed at Rocketdyne CFD Technology Center under the corporate agreement between Boeing North American, Inc. and NASA for the Aerospace Industry Technology Program. The calculations are considered as benchmark testing of the functionality, accuracy, and performance of the solver. Results of these computations demonstrate that the solver is capable of solving the propagation of aeroacoustic signals. Testing of sound generation and on more realistic problems is now pursued for the industrial applications of this solver. Numerical calculations were performed for the second problem of Category 1 of the current workshop problems for an acoustic pulse scattered from a rigid circular cylinder, and for two of the first CAA workshop problems, i. e., the first problem of Category 1 for the propagation of a linear wave and the first problem of Category 4 for an acoustic pulse reflected from a rigid wall in a uniform flow of Mach 0.5. The aim for including the last two problems in this workshop is to test the effectiveness of some boundary conditions set up in the solver. Numerical results of the last two benchmark problems have been compared with their corresponding exact solutions and the comparisons are excellent. This demonstrates the high fidelity of the solver in handling wave propagation problems. This feature lends the method quite attractive in developing a computational acoustic solver for calculating the aero/hydrodynamic noise in a violent flow environment.

Lin, Wen H.↗