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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Optical systolic solutions of linear algebraic equations

The philosophy and data encoding possible in systolic array optical processor (SAOP) were reviewed. The multitude of linear algebraic operations achievable on this architecture is examined. These operations include such linear algebraic algorithms as: matrix-decomposition, direct and indirect solutions, implicit and explicit methods for partial differential equations, eigenvalue and eigenvector calculations, and singular value decomposition. This architecture can be utilized to realize general techniques for solving matrix linear and nonlinear algebraic equations, least mean square error solutions, FIR filters, and nested-loop algorithms for control engineering applications. The data flow and pipelining of operations, design of parallel algorithms and flexible architectures, application of these architectures to computationally intensive physical problems, error source modeling of optical processors, and matching of the computational needs of practical engineering problems to the capabilities of optical processors are emphasized.

Neuman, C. P.↗

NASA NTRS1972

Free adjustment of a geometric global satellite network (solution MPS-7)

The basic purpose of this experiment was to compute reduced normal equations from the observational data of several different systems described below to combine them eventually with the normal equations of the Wild BC-4 observations taken in the DOD/DOC cooperative worldwide geodetic satellite program and provide station coordinates from a single least squares adjustment. The solution described is a partial one obtained without the use of the BC-4 data. The observational systems combined were the Baker-Nunn simultaneous camera observational systems combined were the Baker-Nunn simultaneous camera observations from the SAO worldwide network, the MOTS and PC-1000 optical observations in North America, miscellaneous camera observations in Europe which were included in the SAO69 solution, and, lastly, a group of optical observations where Baker-Nunn cameras observed simultaneously with MOTS and/or PC-1000 cameras in the previously mentioned group.

Mueller, I. I.↗

NASA NTRS1984

Definitive orbit determination for the HEAO-2 spacecraft

Precise ephemerides for the High Energy Astronomy Observatory-2 (HEAO-2) were computed to assist in the Charles Stark Draper Laboratory development of an onboard orbit determination technique. Weighted least-squares, batch orbit solutions were calculated using a high-precision earth gravity model and an approximate model for intermittent spacecraft thrusting. With these improvements, orbit solution consistencies at the 50- to 100-meter level were attained using the Goddard Trajectory Determination System and NASA S-band tracking data.

Smith, R. L.↗

NASA NTRS1988

Development of the NASA/FLAGRO computer program

A preliminary version of a new computer program has been developed for calculating crack propagation in cyclically loaded structures. The new program, titled NASA/FLAGRO, has numerous enhanced features over previous programs used for safe life analysis on space systems. The essential features of the new program consist of modular design, promptive input using a cathode ray tube terminal, improved crack growth rate analytical methods, many new crack case solutions, and incorporation of user-defined files for material properties and fatigue spectrums. Also, five options made available in the program are to compute: (1) safe life of a part with a preexisting crack; (2) critical crack size for a specified stress level; (3) stress-intensity factor numerical values for making comparisons or checks of solutions; (4) least-squares fit of growth rate equations to crack growth rate data; and (5) the conversion of a versus N data to Delta K versus da/dN.

Forman, Royce G.↗

NASA NTRS1989

Optimal moving grids for time-dependent partial differential equations

Various adaptive moving grid techniques for the numerical solution of time-dependent partial differential equations were proposed. The precise criterion for grid motion varies, but most techniques will attempt to give grids on which the solution of the partial differential equation can be well represented. Moving grids are investigated on which the solutions of the linear heat conduction and viscous Burgers' equation in one space dimension are optimally approximated. Precisely, the results of numerical calculations of optimal moving grids for piecewise linear finite element approximation of partial differential equation solutions in the least squares norm.

Wathen, A. J.↗

NASA NTRS2013

Wind Retrieval Algorithms for the IWRAP and HIWRAP Airborne Doppler Radars with Applications to Hurricanes

Algorithms for the retrieval of atmospheric winds in precipitating systems from downward-pointing, conically-scanning airborne Doppler radars are presented. The focus in the paper is on two radars: the Imaging Wind and Rain Airborne Profiler(IWRAP) and the High-altitude IWRAP (HIWRAP). The IWRAP is a dual-frequency (Cand Ku band), multi-beam (incidence angles of 30 50) system that flies on the NOAAWP-3D aircraft at altitudes of 2-4 km. The HIWRAP is a dual-frequency (Ku and Kaband), dual-beam (incidence angles of 30 and 40) system that flies on the NASA Global Hawk aircraft at altitudes of 18-20 km. Retrievals of the three Cartesian wind components over the entire radar sampling volume are described, which can be determined using either a traditional least squares or variational solution procedure. The random errors in the retrievals are evaluated using both an error propagation analysis and a numerical simulation of a hurricane. These analyses show that the vertical and along-track wind errors have strong across-track dependence with values of 0.25 m s-1 at nadir to 2.0 m s-1 and 1.0 m s-1 at the swath edges, respectively. The across-track wind errors also have across-track structure and are on average, 3.0 3.5 m s-1 or 10 of the hurricane wind speed. For typical rotated figure four flight patterns through hurricanes, the zonal and meridional wind speed errors are 2 3 m s-1.Examples of measured data retrievals from IWRAP during an eyewall replacement cycle in Hurricane Isabel (2003) and from HIWRAP during the development of Tropical Storm Matthew (2010) are shown.

atmospheric winds↗

NASA NTRS1977

Ambiguity resolution for satellite Doppler positioning systems

A test for ambiguity resolution was derived which was the most powerful in the sense that it maximized the probability of a correct decision. When systematic error sources were properly included in the least squares reduction process to yield an optimal solution, the test reduced to choosing the solution which provided the smaller valuation of the least squares loss function. When systematic error sources were ignored in the least squares reduction, the most powerful test was a quadratic form comparison with the weighting matrix of the quadratic form obtained by computing the pseudo-inverse of a reduced rank square matrix. A formula is presented for computing the power of the most powerful test. A numerical example is included in which the power of the test is computed for a situation which may occur during an actual satellite aided search and rescue mission.

Argentiero, P. D.↗

NASA NTRS1993

Non-oscillatory and non-diffusive solution of convection problems by the iteratively reweighted least-squares finite element method

A comparative description is presented for the least-squares FEM (LSFEM) for 2D steady-state pure convection problems. In addition to exhibiting better control of the streamline derivative than the streamline upwinding Petrov-Galerkin method, numerical convergence rates are obtained which show the LSFEM to be virtually optimal. The LSFEM is used as a framework for an iteratively reweighted LSFEM yielding nonoscillatory and nondiffusive solutions for problems with contact discontinuities; this method is shown to convect contact discontinuities without error when using triangular and bilinear elements.

Jiang, Bo-Nan↗

NASA NTRS1990

Local and global accuracy estimates for boundary element analysis

An a posteriori error estimation strategy for the boundary element solution of two-dimensional heat conduction problem is presented. The error estimate formulation is based on the boundary element equations and employs a least squares fit of the solutions obtained from these equations. A local error estimate is evaluated for the unknown quantity in each element. The global error estimates are found from local (element) quantities. The range and accuracy of the proposed strategy is very good as demonstrated by problems containing smooth and non-smooth solutions.

Rencis, Joseph J.↗

NASA NTRS1976

On estimating gravity anomalies: A comparison of least squares collocation with least squares techniques

The least squares collocation algorithm for estimating gravity anomalies from geodetic data is shown to be an application of the well known regression equations which provide the mean and covariance of a random vector (gravity anomalies) given a realization of a correlated random vector (geodetic data). It is also shown that the collocation solution for gravity anomalies is equivalent to the conventional least-squares-Stokes' function solution when the conventional solution utilizes properly weighted zero a priori estimates. The mathematical and physical assumptions underlying the least squares collocation estimator are described, and its numerical properties are compared with the numerical properties of the conventional least squares estimator.

Argentiero, P.↗

NASA NTRS1977

On estimating gravity anomalies - A comparison of least squares collocation with conventional least squares techniques

Argentiero, P.↗

NASA NTRS1996

First-Order System Least-Squares for the Navier-Stokes Equations

This paper develops a least-squares approach to the solution of the incompressible Navier-Stokes equations in primitive variables. As with our earlier work on Stokes equations, we recast the Navier-Stokes equations as a first-order system by introducing a velocity flux variable and associated curl and trace equations. We show that the resulting system is well-posed, and that an associated least-squares principle yields optimal discretization error estimates in the H(sup 1) norm in each variable (including the velocity flux) and optimal multigrid convergence estimates for the resulting algebraic system.

Bochev, P.↗

NASA NTRS1973

Constrained least-squares analysis of petrologic problems with an application to lunar sample 12040.

Many petrologic problems, which may be expressed as a set of linear equations, have been solved by least-squares analysis. In many cases insufficient attention has been paid to the physical conditions of the model resulting in incorrect application of the method. This paper presents a systematic treatment of the application of least-squares analysis to petrologic problems including the direct utilization of physical constraints and weighting factors in the problem, and the assessment of uncertainties in the solution. As an example, least-squares analysis is used to examine, in detail, the mass balance equations for lunar rock 12040 and to determine the consistency of the available analytical data.

Reid, M. J.↗

NASA NTRS1984

Acoustooptic linear algebra processors - Architectures, algorithms, and applications

Architectures, algorithms, and applications for systolic processors are described with attention to the realization of parallel algorithms on various optical systolic array processors. Systolic processors for matrices with special structure and matrices of general structure, and the realization of matrix-vector, matrix-matrix, and triple-matrix products and such architectures are described. Parallel algorithms for direct and indirect solutions to systems of linear algebraic equations and their implementation on optical systolic processors are detailed with attention to the pipelining and flow of data and operations. Parallel algorithms and their optical realization for LU and QR matrix decomposition are specifically detailed. These represent the fundamental operations necessary in the implementation of least squares, eigenvalue, and SVD solutions. Specific applications (e.g., the solution of partial differential equations, adaptive noise cancellation, and optimal control) are described to typify the use of matrix processors in modern advanced signal processing.

Casasent, D.↗

NASA NTRS2003

Lunar Rotation and the Lunar Interior

Variations in rotation and orientation of the Moon are sensitive to solid-body tidal dissipation, dissipation due to relative motion at the fluid-core/ solid-mantle boundary, and tidal Love number k2. There is weaker sensitivity to flattening of the core-mantle boundary (CMB) and fluid core moment of inertia. Accurate Lunar Laser Ranging (LLR) measurements of the distance from observatories on the Earth to four retroreflector arrays on the Moon are sensitive to lunar rotation and orientation variations and tidal displacements. Past solutions using the LLR data have given results for dissipation due to solid-body tides and fluid core plus Love number. Past detection of CMB flattening has been marginal but is improving, while direct detection of the core moment has not yet been achieved. Three decades of Lunar Laser Ranging (LLR) data are analyzed using a weighted least-squares approach. The lunar solution parameters include dissipation at the fluid-core/solid-mantle boundary, tidal dissipation, dissipation-related coefficients for rotation and orientation terms, potential Love number k2, a correction to the constant term in the tilt of the equator to the ecliptic which is meant to approximate the influence of core-mantle boundary flattening, and displacement Love numbers h2 and l2. Several solutions, with different combinations of solution parameters and constraints, are considered.

Williams, J. G.↗

NASA NTRS1990

Global configuration of a magnetic cloud

A magnetic cloud associated with a 2N flare on January 1, 1978 was observed by IMP-8, Helios A, Helios B, and Voyager 2. The variation of the magnetic field observed at each spacecraft is represented to good approximation by Lundquist's solution for a cylindrically symmetric force-free magnetic field with constant alpha. A least-squares fit of Lundquist's solution to the data from each spacecraft gives the local orientation of the axis of the magnetic cloud. The times of the estimated boundaries of the magnetic cloud at each spacecraft, together with the speeds of the boundaries and the spacecraft position, give the positions of the boundaries at a given time. From these results the magnetic cloud is determined to resemble a flux rope whose minor radius is approximately 0.15 AU at 1 AU, and whose radius of curvature at 1 AU is approximately 1/3 AU.

Burlaga, L. F.↗

NASA NTRS1993

Comparison of TOPEX/Poseidon orbit determination solutions obtained by the Goddard Space Flight Center Flight Dynamics Division and Precision Orbit Determination Teams

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using the Goddard Trajectory Determination System (GTDS) and a real-time extended Kalman filter estimation system to process Tracking Data and Relay Satellite (TDRS) System (TDRSS) measurements in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. GTDS is the operational orbit determination system used by the FDD, and the extended Kalman fliter was implemented in an analysis prototype system, the Real-Time Orbit Determination System/Enhanced (RTOD/E). The Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generates an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the Geodynamics (GEODYN) orbit determination system with laser ranging tracking data. The TOPEX/Poseidon trajectories were estimated for the October 22 - November 1, 1992, timeframe, for which the latest preliminary POD results were available. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch cases were assessed using overlap comparisons, while the sequential cases were assessed with covariances and the first measurement residuals. The batch least-squares and forward-filtered RTOD/E orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 10 meters (m) for the batch least squares and less than 18 m for the sequential estimation solutions. The differences among the POD, GTDS, and RTOD/E solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C.↗

NASA NTRS1979

Ambiguity resolution for satellite Doppler positioning systems

The implementation of satellite-based Doppler positioning systems frequently requires the recovery of transmitter position from a single pass of Doppler data. The least-squares approach to the problem yields conjugate solutions on either side of the satellite subtrack. It is important to develop a procedure for choosing the proper solution which is correct in a high percentage of cases. A test for ambiguity resolution which is the most powerful in the sense that it maximizes the probability of a correct decision is derived. When systematic error sources are properly included in the least-squares reduction process to yield an optimal solution the test reduces to choosing the solution which provides the smaller valuation of the least-squares loss function. When systematic error sources are ignored in the least-squares reduction, the most powerful test is a quadratic form comparison with the weighting matrix of the quadratic form obtained by computing the pseudoinverse of a reduced-rank square matrix. A formula for computing the power of the most powerful test is provided. Numerical examples are included in which the power of the test is computed for situations that are relevant to the design of a satellite-aided search and rescue system.

Argentiero, P.↗