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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1973

An application of a linear programing technique to nonlinear minimax problems

A differential correction technique for solving nonlinear minimax problems is presented. The basis of the technique is a linear programing algorithm which solves the linear minimax problem. By linearizing the original nonlinear equations about a nominal solution, both nonlinear approximation and estimation problems using the minimax norm may be solved iteratively. Some consideration is also given to improving convergence and to the treatment of problems with more than one measured quantity. A sample problem is treated with this technique and with the least-squares differential correction method to illustrate the properties of the minimax solution. The results indicate that for the sample approximation problem, the minimax technique provides better estimates than the least-squares method if a sufficient amount of data is used. For the sample estimation problem, the minimax estimates are better if the mathematical model is incomplete.

Schiess, J. R.↗
NASA NTRS1978

Estimate of procession and polar motion errors from planetary encounter station location solutions

Jet Propulsion Laboratory Deep Space Station (DSS) location solutions based on two JPL planetary ephemerides, DE 84 and DE 96, at eight planetary encounters were used to obtain weighted least squares estimates of precession and polar motion errors. The solution for precession error in right ascension yields a value of 0.3 X 10 to the minus 5 power plus or minus 0.8 X 10 to the minus 6 power deg/year. This maps to a right ascension error of 1.3 X 10 to the minus 5 power plus or minus 0.4 X 10 to the minus 5 power deg at the first Voyager 1979 Jupiter encounter if the current JPL DSS location set is used. Solutions for precession and polar motion using station locations based on DE 84 agree well with the solution using station locations referenced to DE 96. The precession solution removes the apparent drift in station longitude and spin axis distance estimates, while the encounter polar motion solutions consistently decrease the scatter in station spin axis distance estimates.

Pease, G. E.↗
NASA NTRS1995

The origin of spurious solutions in computational electromagnetics

The origin of spurious solutions in computational electromagnetics, which violate the divergence equations, is deeply rooted in a misconception about the first-order Maxwell's equations and in an incorrect derivation and use of the curl-curl equations. The divergence equations must be always included in the first-order Maxwell's equations to maintain the ellipticity of the system in the space domain and to guarantee the uniqueness of the solution and/or the accuracy of the numerical solutions. The div-curl method and the least-squares method provide rigorous derivation of the equivalent second-order Maxwell's equations and their boundary conditions. The node-based least-squares finite element method (LSFEM) is recommended for solving the first-order full Maxwell equations directly. Examples of the numerical solutions by LSFEM for time-harmonic problems are given to demonstrate that the LSFEM is free of spurious solutions.

Jiang, Bo-Nan↗
NASA NTRS1978

On least squares approximations to indefinite problems of the mixed type

A least squares method is presented for computing approximate solutions of indefinite partial differential equations of the mixed type such as those that arise in connection with transonic flutter analysis. The method retains the advantages of finite difference schemes namely simplicity and sparsity of the resulting matrix system. However, it offers some great advantages over finite difference schemes. First, the method is insensitive to the value of the forcing frequency, i.e., the resulting matrix system is always symmetric and positive definite. As a result, iterative methods may be successfully employed to solve the matrix system, thus taking full advantage of the sparsity. Furthermore, the method is insensitive to the type of the partial differential equation, i.e., the computational algorithm is the same in elliptic and hyperbolic regions. In this work the method is formulated and numerical results for model problems are presented. Some theoretical aspects of least squares approximations are also discussed.

Fix, G. J.↗
NASA NTRS1988

Accuracy study of the Upper Atmosphere Research Satellite (UARS) definitive attitude determination

The Upper Atmosphere Research Satellite (UARS) has two definitive attitude determination requirements: the definitive attitude of the Modular Attitude Control Subsystem (MACS) and the definitive attitude of the gimbaled Solar-Stellar Pointing Platform (SSPP). The onboard computer (OBC) will compute the MACS attitude using a Kalman filter and will transform this attitude solution through the SSPP gimbals to calculate the SSPP attitude. The attitude ground support system (AGSS) will compute the MACS attitude using a batch least-squares differential corrector algorithm and will also transform this solution through the gimbals to obtain the SSPP attitude. This paper reports the results of a prelaunch study to predict the accuracy of the OBC attitude solutions and the accuracy of the AGSS attitude solutions. The OBC and AGSS solution accuracies are then compared to establish the relative quality. The effects of star observability, sensor noise, and sensor misalignment uncertainties on attitude determination accuracy are analyzed for each case.

Snow, Frank↗
NASA NTRS1993

Analyzing Martian winds and tracer concentrations using Mars Observer data

During the courses of a day, the Mars Observer spacecraft will acquire globally distributed profiles of the martian atmosphere. It is highly desirable that this data be assembled into synoptic weather maps (complete specifications of the atmospheric pressure, temperature, and winds at a given time), which can in turn be used as starting points in the study of many meteorological phenomena. Unfortunately, the special nature of the Mars Observer data presents several challenges above and beyond the usual difficult problem of data initialization. Mars Observer atmospheric data will consist almost exclusively of asynoptic vertical profiles of temperatures (or radiances) and pressures, whereas winds are generally in balance with horizontal gradients of these quantities (which will not be observed). It will therefore be necessary to resort to dynamical models to analyze the wind fields. As a rule, data assimilation into atmospheric models can result in the generation of spurious gravity waves, so special steps must be taken to suppress these. In addition, the asynoptic nature of the data will require a four-dimensional (space and time) data assimilation scheme. The problem is to find a full set of meteorological fields (winds and temperatures) such that, when marched forward in time in the model, they achieve a best fit (in the weighted least-squares sense) to the data. The proposed solution is to develop a model especially for the Mars Observer data assimilation problem. Gravity waves are filtered from the model by eliminating all divergence terms from the prognostic divergence equation. This leaves a diagnostic gradient wind relation between the rotational wind and the temperature field. The divergent wind is diagnosed as the wind required to maintain the gradient wind balance in the presence of the diabatic heating. The primitive equations of atmospheric dynamics (with three principal dependent variables) are thus reduced to a simpler system with a single prognostic equation for temperature - the variable that will be best observed. (This balance system was apparently first derived by Charney as a first-order Rossby number expansion of the equations of motion). Experience with a full primitive equation model of the martian atmosphere indicates that a further simplification is possible: at least for short-term integrations, the model can be linearized about the zonally symmetric basic state.

Houben, Howard C.↗
NASA NTRS1981

Remarks on nonstandard least-squares problems

It is shown that after treating a least-squares problem in the standard fashion, the imposition of stochastic restrictions on certain parameters can be accomplished essentially by inverting a matrix of an order no larger than the total number of restricted parameters, and that generally there is an explicit term for correcting the unrestricted to the restricted solution. In addition, any known type of nonstandard least-squares adjustment problem can be reduced to an equivalent standard linear problem, that is, one in which there is exactly one observation per condition equation and in which the observations are not correlated.

Eichhorn, H.↗
NASA NTRS1996

A Systematic Solution Approach for Neutron Transport Problems in Diffuse Regimes

A systematic solution approach for the neutron transport equation, based on a least-squares finite-element discretization, is presented. This approach includes the theory for the existence and uniqueness of the analytical as well as of the discrete solution, bounds for the discretization error, and guidance for the development of an efficient multigrid solver for the resulting discrete problem. To guarantee the accuracy of the discrete solution for diffusive regimes, a scaling transformation is applied to the transport operator prior to the discretization. The key result is the proof of the V-ellipticity and continuity of the scaled least-squares bilinear form with constants that are independent of the total cross section and the absorption cross section. For a variety of least-squares finite-element discretizations this leads to error bounds that remain valid in diffusive regimes. Moreover, for problems in slab geometry a full multigrid solver is presented with V(1, 1)-cycle convergence rates approximately equal to 0.1, independent of the size of the total cross section and the absorption cross section.

Manteuffel, T. A.↗
NASA NTRS1976

Techniques for obtaining regional radiation budgets from satellite radiometer observations, phase 4 and phase 5

A scheme was developed which divides the earth-atmosphere system into 2060 elemental areas. The regions previously described are defined in terms of these elemental areas which are fixed in size and position as the satellite moves. One method, termed the instantaneous technique, yields values of the radiant emittance (We) and the radiant reflectance (Wr) which the regions have during the time interval of a single satellite pass. The number of observations matches the number of regions under study and a unique solution is obtained using matrix inversion. The other method (termed the best fit technique), yields time averages of We and Wr for large time intervals (e.g., months, seasons). The number of observations in this technique is much greater than the number of regions considered, and an approximate solution is obtained by the method of least squares.

Pina, J. F.↗
NASA NTRS1994

Preliminary solutions for the lunar gravity field from analysis of lunar orbiter tracking data

Knowledge of the gravitation field, in combination with surface topography, provides one of the principal means of inferring the internal structure of a planetary body. Previous analyses of the lunar gravitational field have been based on data from the Lunar Orbiters, the Apollo subsatellites, and the low altitude passes of the Apollo spacecraft. Recently, Konopliv et al. have reanalyzed all available Lunar Orbiter and Apollo subsatellite tracking data, producing a 60th degree and order solution. In preparation for the Clementine Mission to the Moon, we have also initiated a reanalysis of the Lunar Orbiter and Apollo subsatellite data. Our reanalysis takes advantage of advanced force and measurement modeling techniques as well as modern computational facilities. We applied the least squares collocation technique which stabilizes the behavior of the solution and high degree and order. The extension of the size of the field reduces the aliasing coming from the omitted portion of the gravitational field. This is especially important for the analysis of the tracking data from the Lunar Orbiters, as the periapse heights frequently ranged from 50 to 100 km.

Lemoine, F. G.↗
NASA NTRS1994

Theoretical study of the incompressible Navier-Stokes equations by the least-squares method

Usually the theoretical analysis of the Navier-Stokes equations is conducted via the Galerkin method which leads to difficult saddle-point problems. This paper demonstrates that the least-squares method is a useful alternative tool for the theoretical study of partial differential equations since it leads to minimization problems which can often be treated by an elementary technique. The principal part of the Navier-Stokes equations in the first-order velocity-pressure-vorticity formulation consists of two div-curl systems, so the three-dimensional div-curl system is thoroughly studied at first. By introducing a dummy variable and by using the least-squares method, this paper shows that the div-curl system is properly determined and elliptic, and has a unique solution. The same technique then is employed to prove that the Stokes equations are properly determined and elliptic, and that four boundary conditions on a fixed boundary are required for three-dimensional problems. This paper also shows that under four combinations of non-standard boundary conditions the solution of the Stokes equations is unique. This paper emphasizes the application of the least-squares method and the div-curl method to derive a high-order version of differential equations and additional boundary conditions. In this paper, an elementary method (integration by parts) is used to prove Friedrichs' inequalities related to the div and curl operators which play an essential role in the analysis.

Jiang, Bo-Nan↗
NASA NTRS1973

Sparse matrix methods based on orthogonality and conjugacy

A matrix having a high percentage of zero elements is called spares. In the solution of systems of linear equations or linear least squares problems involving large sparse matrices, significant saving of computer cost can be achieved by taking advantage of the sparsity. The conjugate gradient algorithm and a set of related algorithms are described.

Lawson, C. L.↗
NASA NTRS1983

On obtaining the forward phase functions of Saturn ring features from radio occultation observations

It is noted that the near-forward scattering functions of particles in Saturn ring features are related to 3.6 cm radio occultation power spectra by a Fredholm integral equation of the first kind. The equation reduces to an algebraic system of equation whose solution by usual inversion techniques (that is, least mean squares) is ruled out by the near singularity of the forward transformation matrix. A combination of constrained linear inversion and a filtering algorithm based on eigenvector decomposition of the matrix reduces the instabilities; this yields derived phase functions valid over the range of zero to about 12 mrad. It is noted that these functions represent the collective forward diffraction lobe of particles greater than about 1 m in radius. Since multiple scattering of the signal is a significant effect, the measured phase functions must be adjusted to obtain the singly scattered component. This single-scattering correction is examined for two physical models, namely the monolayer and the classical discrete random slab; in addition, the fraction of opacity in submeter particles for each model for particular ring features is estimated.

Zebker, H. A.↗
NASA NTRS1987

Source parameters for 11 earthquakes in the Tien Shan, central Asia, determined by P and SH waveform inversion

The style and the distribution of faulting occurring today in the Tien Shan region were studied, by digitizing long-period World-Wide Standard Seismograph Network P and SH waveforms of 11 of the largest Tien Shan earthquakes between 1965 and 1982 and then using a least squares inversion routine to constrain their fault plane solutions and depths. The results of the examination indicate that north-south shortening is presently occurring in the Tien Shan, with the formation of basement uplifts flanked by moderately dipping thrust faults. The present-day tectonics of the Tien Shan seem to be analogous to those of the Rocky Mountains in Colorado, Wyoming, and Utah during the Laramide orogeny in Late Cretaceous and Early Tertiary time.

Nelson, Michael R.↗
NASA NTRS1976

Additional comments on 'Multistage least-squares parameter estimation'

The purpose of this correspondence is to call attention to other solutions to the problem addressed in the above paper, that of adding or deleting parameters or measurements from a given linear model. The solution algorithms proposed in the article are numerically ill-conditioned and are less efficient, in terms of both computational and storage requirements, than the methods described here. Several numerically stable, compact, flexible, and efficient algorithms for analyzing the parameter estimation problem are described and others are referenced.

Bierman, G. J.↗
NASA NTRS1987

Simplified solution for a class of fading memory filters

A fading memory filter is a least squares estimator (LSE) that applies an exponentially decaying weight to past measurements. When compared with a standard Kalman filter, its key advantages are asymptotic stability and reduced sensitivity to modeling errors. This paper derives a simple solution for a class of fading memory filters, resulting in a reduction in computational complexity. Steady state filter solutions are obtained for second- and third-order filters used in a global positioning system (GPS) receiver for high dynamic vehicles.

Statman, Joseph I.↗
NASA NTRS1975

Four methods of attitude determination for spin-stabilized spacecraft with applications and comparative results

The attitude of a spacecraft is determined by specifying independent parameters which relate the spacecraft axes to an inertial coordinate system. Sensors which measure angles between spin axis and other vectors directed to objects or fields external to the spacecraft are discussed. For the spin-stabilized spacecraft considered, the spin axis is constant over at least an orbit, but separate solutions based on sensor angle measurements are different due to propagation of errors. Sensor-angle solution methods are described which minimize the propagated errors by making use of least squares techniques over many sensor angle measurements and by solving explicitly (in closed form) for the spin axis coordinates. These methods are compared with star observation solutions to determine if satisfactory accuracy is obtained by each method.

Smith, G. A.↗