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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

State and model error estimation for distributed parameter systems

In-flight estimation of large structure model errors in order to detect inevitable deficiencies in large structure controller/estimator models is discussed. Such an estimation process is particularly applicable in the area of shape control system design required to maintain a prescribed static structural shape and, in addition, suppress dynamic disturbances due to the vehicle vibrational modes. The paper outlines a solution to the problem of static shape estimation where the vehicle shape must be reconstructed from a set of measurements discretely located throughout the structure. The estimation process is based on the principle of least-squares that inherently contains the definition and explicit computation of model error estimates that are optimal in some sense. Consequently, a solution is provided for the problem of estimation of static model errors (e.g., external loads). A generalized formulation applicable to distributed parameters systems is first worked out and then applied to a one-dimensional beam-like structural configuration.

Rodriguez, G.↗
NASA NTRS1978

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

An analysis of the shortening and lengthening of the phase of satellite-to-satellite (SST) data that passed within 40-700 km above the earth surface during its ATS-6 to GEOS-3 to ATS-6 path resulted in refractivity vs height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbations from the moon, sun, and a 15th-order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occulation period and a model ionosphere was estimated by a least-squares adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the two methods were caused by orbital errors, which propagated into the solution.

Liu, A. S.↗
NASA NTRS1978

On the determination and investigation of the terrestrial ionospheric refractive indices using GEOS-3/ATS-6 satellite-to-satellite tracking data

When the radio link between two satellites (GEOS-3/ATS-6) is intercepted by the earth's ionosphere and neutral atmosphere, a change in the Doppler frequency results. Travel through the atmosphere causes the Doppler phase to be advanced in the ionosphere's portion and retarded in the neutral portion of the atmosphere. Analysis of the shortening and lengthening of the phase of the Satellite-to-Satellite Tracking (SST) data that passed within 40-700 km above the earth's surface during its ATS-6 to GEOS-3 to ATS-6 path, caused by the atmosphere, results in refractivity versus height profiles. The SST Doppler data were used directly to adjust the GEOS-3 orbit. Perturbation from the Moon, Sun and a 15th order/degree earth gravity field were included in the orbit solution. This orbit was continued through the occultation period and a model ionosphere was estimated by a least-square adjustment of the Chapman ionosphere parameters from the SST data residuals. The refractivity profile obtained by this model ionosphere was compared to a refractivity profile obtained by a direct integral inversion of the SST data residuals. Systematic differences between the 2 methods were caused by orbital errors, which propagated into the solution. The SST data yielded refractive index profiles in a novel economical manner because no additional or special on-board equipment were required.

Liu, A. S.↗
NASA NTRS1976

Simplified solution for point contact deformation between two elastic solids

A linear-regression by the method of least squares is made on the geometric variables that occur in the equation for point contact deformation. The ellipticity and the complete eliptic integrals of the first and second kind are expressed as a function of the x, y-plane principal radii. The ellipticity was varied from 1 (circular contact) to 10 (a configuration approaching line contact). These simplified equations enable one to calculate easily the point-contact deformation to within 3 percent without resorting to charts or numerical methods.

Brewe, D. E.↗
NASA NTRS1977

Simplified solution for elliptical-contact deformation between two elastic solids

A linear regression by the method of least squares is made on the geometric variables that occur in the equation for elliptical-contact deformation. The ellipticity and the complete elliptic integrals of the first and second kind are expressed as a function of the x,y-plane principal radii. The ellipticity was varied from 1 (circular contact) to 10 (a configuration approaching line contact). The procedure for solving these variables without the use of charts or a high-speed computer would be quite tedious. These simplified equations enable one to calculate easily the elliptical-contact deformation to within 3 percent accuracy without resorting to charts or numerical methods.

Brewe, D. E.↗
NASA NTRS1990

Numerical experiments in geomagnetic modeling

Numerical tests were made, using least squares fitting of a spherical harmonic model, to a selection of Magsat data to determine the practical limits of this technique with modern computers. The resulting (M102189) model, whose coefficients were adjusted up to n = 50, was compared with M07AV6, a previous model which used least squares (on vector data) for coefficients up to n = 29, and Gauss-Legendre quadrature (on Z residuals) to adjust the coefficients up to n = 63. For the new least squares adjustment to n = 50 a condition number of 115 was obtained for the solution matrix, with a resulting precision of 11 significant figures. The M102189 model shows a lower and more Gaussian residual distribution than did M07AV6, though the Gaussian envelope fits to the residual distributions, even for the scalar field, gives "standard deviations' never lower than 6 nT, a factor of three higher than the estimated Magsat observational errors. Ionospheric currents are noted to have a significant effect on the coefficients of the internal potential functions.

Cain, Joseph C.↗
NASA NTRS1990

Polarization studies of Zeeman affected spectral lines using the MSFC magnetograph

The Marshall Space Flight Center (MSFC) vector magnetograph records polarization images of absorption lines that are sensitive to magnetic fields. A method is presented for analyzing the Stokes spectral-line profiles of a photospheric Fe I absorption line (5250.2 A) which is influenced by the Zeeman effect. Using nonlinear least-square optimization, the observed Stokes profiles are compared with those generated from the theoretical solution of the polarized radiative transfer equations. The optimization process accounts for the spectral convolution of the source and the MSFC vector magnetograph. The resulting physical properties of the active region producing the polarized light are discussed.

Balasubramaniam, K. S.↗
NASA NTRS1993

Calibration of a synthetic aperture radiometer

Calibration algorithms for a synthetic aperture microwave radiometer are presented. The calibration is geared to Earth remote sensing applications, and is demonstrated on an airborne prototype of thinned array imager. Two approaches to the system calibration are presented; the first utilizes commonly available reference brightness temperature scenes, such as open water, and the second utilizes data collected on the antenna range. Both algorithms yield spatial response information which is cast in matrix form and inverted to obtain the image reconstruction formula. Experimental results are examined, and errors in some reconstructed images are linked to the present prototype antenna design. Algorithms for improving the synthesized antenna pattern side lobe performance are also presented. In one solution, the pattern efficiency is optimized by minimizing the pattern outside a defined beam. In another solution, the patterns are matched to a desired model pattern by the method of least squared errors. Both techniques offer an attractive alternative to aperture weighting.

Tanner, Alan B.↗
NASA NTRS1994

Robust inverse kinematics using damped least squares with dynamic weighting

This paper presents a general method for calculating the inverse kinematics with singularity and joint limit robustness for both redundant and non-redundant serial-link manipulators. Damped least squares inverse of the Jacobian is used with dynamic weighting matrices in approximating the solution. This reduces specific joint differential vectors. The algorithm gives an exact solution away from the singularities and joint limits, and an approximate solution at or near the singularities and/or joint limits. The procedure is here implemented for a six d.o.f. teleoperator and a well behaved slave manipulator resulted under teleoperational control.

Schinstock, D. E.↗
NASA NTRS1999

A Method to Solve Interior and Exterior Camera Calibration Parameters for Image Resection

An iterative method is presented to solve the internal and external camera calibration parameters, given model target points and their images from one or more camera locations. The direct linear transform formulation was used to obtain a guess for the iterative method, and herein lies one of the strengths of the present method. In all test cases, the method converged to the correct solution. In general, an overdetermined system of nonlinear equations is solved in the least-squares sense. The iterative method presented is based on Newton-Raphson for solving systems of nonlinear algebraic equations. The Jacobian is analytically derived and the pseudo-inverse of the Jacobian is obtained by singular value decomposition.

Samtaney, Ravi↗
NASA NTRS2017

Performance Optimizing Adaptive Control with Time-Varying Reference Model Modification

This paper presents a new adaptive control approach that involves a performance optimization objective. The control synthesis involves the design of a performance optimizing adaptive controller from a subset of control inputs. The resulting effect of the performance optimizing adaptive controller is to modify the initial reference model into a time-varying reference model which satisfies the performance optimization requirement obtained from an optimal control problem. The time-varying reference model modification is accomplished by the real-time solutions of the time-varying Riccati and Sylvester equations coupled with the least-squares parameter estimation of the sensitivities of the performance metric. The effectiveness of the proposed method is demonstrated by an application of maneuver load alleviation control for a flexible aircraft.

Adaptive Control↗
NASA NTRS2018

Time-Varying Modification of Reference Model for Adaptive Control with Performance Optimization

This paper presents a new adaptive control ap- proach that involves a performance optimization objective. The control synthesis involves the design of a performance optimizing adaptive controller from a subset of control inputs. The resulting effect of the performance optimizing adaptive controller is to modify the initial reference model into a time-varying reference model which satisfies the performance optimization requirement obtained from an optimal control problem. The time-varying reference model modification is accomplished by the real-time solutions of the time-varying Riccati and Sylvester equations coupled with the least-squares parameter estimation of the sen- sitivities of the performance metric. The effectiveness of the proposed method is demonstrated by an application of maneuver load alleviation control for a flexible aircraft.

Adaptive Control↗
NASA NTRSDate not supplied

Time-Varying Reference Model Modification of Adaptive Control with Multi-Objective Performance Optimization

This paper presents an adaptive control approach that modifies a reference model in order to satisfy a multi-objective optimization of two performance signals: a linear signal and a quadratic signal. The time varying reference model modification is accomplished by the real-time solutions of the time-varying Riccati and Sylvester equations coupled with the least-squares parameter estimation of the sensitivities of the performance signals to be minimized. The effectiveness of the proposed method is to be demonstrated by a flight control application of maneuver load alleviation and drag minimization.

Adaptive Control↗
NASA NTRS1984

2 Pallas - 1982 and 1983 lightcurves and a new pole solution

Photoelectric lightcurves of the asteroid 2 Pallas obtained in March 1982 and May 1983 display amplitudes of 0.04 and 0.10 mag respectively. The latter lightcurve shows that Pallas was at V(1.0) mag of 4.51 + or -0.02 when it occulted 1 Vulpeculae on May 29, 1983. A least-squares best fit to an amplitude-aspect relation for all available lightcurve observations of Pallas between 1951 and 1983 yields two solutions for its pole position: lambda = 200, beta = 40 and lambda = 220, beta = 15, where the uncertainty regions correspond to an overall error estimate of + or -10 deg. Use of phase angle bisector coordinates give lower residuals than geocentric coordinates. The (220, 15) pole position is favored since it is in very good agreement with an independent pole solution obtained by photometric astrometry. This pole position implies that the latitude of the sub-earth point at the time of occultation was 22 deg.

Binzel, R. P.↗
NASA NTRS1989

A least-squares finite element method for incompressible Navier-Stokes problem

Most finite element schemes for solving the Navier-Stokes equations can be categorized into the Galerkin mixed method and the penalty method. The mixed method leads to a saddle-point problem. In order to guarantee the existence of a solution, the combination of velocity and pressure interpolations requires satisfaction of the Ladyzhenskaya Babuska Brezzi (LBB) consistency condition which precludes the use of equal order interpolations and many seemingly natural pairs of velocity and pressure elements. In a previous paper a least-squares finite element method based on the first order velocity-pressure-vorticity formulation for the Stokes problem was proposed. This method leads to a minimization problem. The choice of combination of elements is thus not subject to the LBB condition. The numerical experiments exhibit the optimal rate of convergence for all variables with equal order interpolations. A theoretical error analysis supports the numerical results. In this paper the least-squares finite element method is extended to solving the incompressible Navier-Stokes problem.

Jiang, Bonan↗
NASA NTRS1987

A method of calculating the total flow from a given sea surface topography

Using a simple dynamical model of a wind-driven ocean circulation of the Stommel type, and an analytical basis developed to objectively analyze the sea surface height residuals from an altimeter and, in the process, to determine the total flow instead of just the near surface geostrophic component associated with the given sea surface topography. The method is based on first deriving the solution to the forced problem for a given wind stress required to develop a hypothetical true or perfect data field and to establishing the basis for the objective analysis. The stream function and the surface height field for the forced problem are developed in terms of certain characteristic functions with the same expansion coefficients for both fields. These characteristic functions are simply the solutions for a homogeneous elliptic equation for the stream function and the solutions of an inhomogeneous balance equation for the height field. For the objective analysis, using a sample of randomly selected height values from the true data field, the height field characteristic functions are used to fit the given topography in a least squares sense. The resulting expansion coefficients then permit the synthesis of the total flow field via the stream function characteristic modes and the solution is perfectly well behaved even along the equator. The method of solution is easily adaptable to realistic ocean basis by straight forward numerical methods. The analytical basis of the theory and the results for an ideal rectangular basin on a beta plane are described.

Rao, Desiraju B.↗
NASA NTRS1997

Geopotential Error Analysis from Satellite Gradiometer and Global Positioning System Observables on Parallel Architecture

The recovery of a high resolution geopotential from satellite gradiometer observations motivates the examination of high performance computational techniques. The primary subject matter addresses specifically the use of satellite gradiometer and GPS observations to form and invert the normal matrix associated with a large degree and order geopotential solution. Memory resident and out-of-core parallel linear algebra techniques along with data parallel batch algorithms form the foundation of the least squares application structure. A secondary topic includes the adoption of object oriented programming techniques to enhance modularity and reusability of code. Applications implementing the parallel and object oriented methods successfully calculate the degree variance for a degree and order 110 geopotential solution on 32 processors of the Cray T3E. The memory resident gradiometer application exhibits an overall application performance of 5.4 Gflops, and the out-of-core linear solver exhibits an overall performance of 2.4 Gflops. The combination solution derived from a sun synchronous gradiometer orbit produce average geoid height variances of 17 millimeters.

Schutz, Bob E.↗