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 289 records · Page 16

M2, S2, K1 models of the global ocean tide

Ocean tidal signals appear in many geophysical measurements. Geophysicists need realistic tidal models to aid in interpretation of their data. Because of the closeness to resonance of dissipationless ocean tides, it is difficult for numerical models to correctly represent the actual open ocean tide. As an approximate solution to this problem, test functions derived by solving Laplace's Tidal Equations with ocean loading and self gravitation are used as a basis for least squares dynamic interpolation of coastal and island tidal data for the constituents M2, S2, and Kl. The resulting representations of the global tide are stable over at least a ?5% variation in the mean depth of the model basin, and they conserve mass. Maps of the geocentric tide, the induced free space potential, the induced vertical component of the solid earth tide, and the induced vertical component of the gravitational field for each contituent are presented.

Parke, M. E.↗

Polynomial compensation, inversion, and approximation of discrete time linear systems

The least-squares transformation of a discrete-time multivariable linear system into a desired one by convolving the first with a polynomial system yields optimal polynomial solutions to the problems of system compensation, inversion, and approximation. The polynomial coefficients are obtained from the solution to a so-called normal linear matrix equation, whose coefficients are shown to be the weighting patterns of certain linear systems. These, in turn, can be used in the recursive solution of the normal equation.

Baram, Yoram↗

A Combined HIPPARCOS and Multichannel Astrometric Photometer Study of the Proposed Planetary System of Rho Coronae Borealis

Hipparcos and Multichannel Astrometric Photometer (MAP) observations of rho Coronae Borealis independently display astrometric motion at the period of the proposed extrasolar planetary companion to the star. Individual least-squares fits to each astrometric data set yield independent estimates of the semimajor axis, inclination, and node angle that are in excellent agreement. A combined solution of the Hipparcos and MAP data yields an inclination of 0.5 deg, a node at 30.5 +/- 12.4, and a semimajor axis of 1.66 +/- 0.35 mas, indicating a companion mass of 0.14 +/- 0.05 solar masses over two orders of magnitude greater than the minimum mass for the companion as determined by radial velocity studies. This mass is approximately that of an M dwarf star, the companion cannot be a planetary object.

Gatewood, George↗

Considerations of open-loop, closed-loop, and adaptive multicyclic control systems

Four different types of self-tuning regulators were studied for multicyclic control of helicopter vibration. A numerical simulation of the helicopter is made, using a multivariable frequency-domain model, in terms of transfer function with six input control harmonics and six output harmonics. The model characteristics vary with flight speed. An off-line identification of model characteristics is made, using the least-squared-error method and using a succession of input and output measurements. The on-line identification of model characteristics is made using the Kalman filter solution. The optimal controls are calculated from the minimization of quadratic performance function based on response and multicyclic inputs. The performance of various regulators or controllers is judged from the stability, transient response, convergence time, and amplitude of the steady state.

Chopra, I.↗

Nearly-grazing optimal trajectories for aeroassisted orbital transfer

In the present treatment of optimal control problems arising in the study of coplanar aeroassisted orbital transfer, the hybrid combination of propulsive parameters in space and aerodynamic maneuvers employing lift modulation in the sensible atmosphere indicates that the optimal energy-viewpoint solution is the grazing trajectory; this trajectory is characterized by favorable values of the peak heating rate and the peak dynamic pressure. Numerical solutions are obtained by means of the sequential gradient restoration algorithm for optimal control problems. It is found that nearly-grazing trajectories yielding the least-square value of the path inclination have desirable characteristics from the standpoints of energy, heating rate, and dynamic pressure.

Miele, A.↗

Plate tectonics from VLBI and SLR global data

This study is based on data derived from fifteen years of observations of the SLR (side-looking radar) network and six years of the VLBI (very long baseline interferometry) network. In order to use all available information VLBI and SLR global data sets were combined in a least squares fashion to calculate station horizontal velocities. All significant data pertaining to a single site contribute to the station horizontal motion. The only constraint on the solution is that no vertical motion is allowed. This restriction does not greatly affect the precision of the overall solution given the fact that the expected vertical motion for most stations, even those experiencing post glacial uplift, is well under 1 cm/yr. Since the average baseline is under 4,000 km, only a small fraction of the station vertical velocity is translated into baseline rates so that the error introduced in the solution by restricting up-down station movement is minimal. As a reference, station velocities were then compared to the ones predicted by the NUVEL-1 geological model of DeMets et al. (1990). The focus of the study is on analyzing these discrepancies for global plate tectonics as well as regional tectonic settings. The method used also allows us not only to derive horizontal motion for individual stations but also to calculate Euler vectors for those plates that have enough stations located on the stable interior like North America, Pacific, Eurasia, and Australia.

Harrison, Christopher G. A.↗

Improved solution accuracy for Landsat-4 (TDRSS-user) orbit determination

This paper presents the results of a study to compare the orbit determination accuracy for a Tracking and Data Relay Satellite System (TDRSS) user spacecraft, Landsat-4, obtained using a Prototype Filter Smoother (PFS), with the accuracy of an established batch-least-squares system, the Goddard Trajectory Determination System (GTDS). The results of Landsat-4 orbit determination will provide useful experience for the Earth Observing System (EOS) series of satellites. The Landsat-4 ephemerides were estimated for the January 17-23, 1991, timeframe, during which intensive TDRSS tracking data for Landsat-4 were available. Independent assessments were made of the consistencies (overlap comparisons for the batch case and convariances for the sequential case) of solutions produced by the batch and sequential methods. The filtered and smoothed PFS orbit solutions were compared with the definitive GTDS orbit solutions for Landsat-4; the solution differences were generally less than 15 meters.

Oza, D. H.↗

Photometric properties of lunar terrains derived from Hapke's equation

Hapke's bidirectional reflectance equations provide the most rigorous available description of photometric behavior in terms of physically meaningful parameters. The primary objective of this study was to derive Hapke parameters for the lunar surface from both disk-integrated and disk-resolved photometric data. Hapke's equation was fit to the disk-integrated phase curves and disk-resolved data for dark, average, and bright terrain classes using an iterative, nonlinear least squares algorithm described by Helfenstein. Parameters were initially determined from the disk-integrated data, and the result was applied as a first guess to the iterative solution of parameters for individual terrain classes. Plots are presented of the disk-resolved data normalized to corresponding brightnesses predicted from the disk-integrated solution under the same illumination and viewing geometries. Systematic trends in disk-resolved parameters can be identified. Values for single scattering albedo (w) of the dominantly anorthositic average and bright terrains are significantly larger than the value for the basaltic dark terrains (mare). Values for surge brightness parameters, particle phase function, and average topographic slope angle are also discussed.

Helfenstein, Paul↗

Turbulent Plane Wakes Subjected to Successive Strains

Six direct numerical simulations of turbulent time-evolving strained plane wakes have been examined to investigate the response of a wake to successive irrotational plane strains of opposite sign. The orientation of the applied strain field has been selected so that the flow is the time-developing analogue of a spatially developing wake evolving in the presence of either a favourable or an adverse streamwise pressure gradient. The magnitude of the applied strain rate a is constant in time t until the total strain e(sup at) reaches about four. At this point, a new simulation is begun with the sign of the applied strain being reversed (the original simulation is continued as well). When the total strain is reduced back to its original value of one, yet another simulation is begun with the sign of the strain being reversed again back to its original sign. This process is done for both initially "favourable" and initially "adverse" strains, providing simulations for each of these strain types from three different initial conditions. The evolution of the wake mean velocity deficit and width is found to be very similar for all the adversely strained cases, with both measures rapidly achieving exponential growth at the rate associated with the cross-stream expansive strain e(sup at). In the "favourably" strained cases, the wake widths approach a constant and the velocity deficits ultimately decay rapidly as e(sup -2at). Although all three of these cases do exhibit the same asymptotic exponential behaviour, the time required to achieve this is longer for the cases that have been previously adversely strained (by at approx. equals 1). These simulations confirm the generality of the conclusions drawn in Rogers (2002) regarding the response of plane wakes to strain. The evolution of strained wakes is not consistent with the predictions of classical self-similar analysis; a more general equilibrium similarity solution is required to describe the results. At least for the cases considered here, the wake Reynolds number and the ratio of the turbulent kinetic energy to the square of the wake mean velocity deficit are determined nearly entirely by the total strain. For these measures the order in which the strains are applied does not matter and the changes brought about by the strain are nearly reversible. The wake mean velocity deficit and width, on the other hand, differ by about a factor of three when the total strain returns to one, depending on whether the wake was first "favourably" or "adversely" strained. The strain history is important for predicting the evolution of these quantities.

Rogers, Michael M.↗

Variability simulations with a steady, linearized primitive equations model

Solutions of the steady, primitive equations on a sphere, linearized about a zonally symmetric basic state are computed for the purpose of simulating monthly mean variability in the troposphere. The basic states are observed, winter monthly mean, zonal means of zontal and meridional velocities, temperatures and surface pressures computed from the 15 year NMC time series. A least squares fit to a series of Legendre polynomials is used to compute the basic states between 20 H and the equator, and the hemispheres are assumed symmetric. The model is spectral in the zonal direction, and centered differences are employed in the meridional and vertical directions. Since the model is steady and linear, the solution is obtained by inversion of a block, pente-diagonal matrix. The model simulates the climatology of the GFDL nine level, spectral general circulation model quite closely, particularly in middle latitudes above the boundary layer. This experiment is an extension of that simulation to examine variability of the steady, linear solution.

Kinter, J. L., III↗

Techniques in Doppler gravity inversion

The types of Doppler gravity data available for local as opposed to planetwide geophysical modeling are reviewed. Those gravity fields that are determined dynamically in orbit determination programs yield a smoothed representation of the local gravity field that may be used for quantitative modeling. An estimate of the difference between smoothed and true fields can be considered as a noise limitation in generating local gravity models. A nonlinear inversion for the geometry, depth, and density of the Mare Serenitatis mascon using an ellipsoidal model yielded a global least squares minimum in horizontal dimensions, depth, and thickness-density contrast product. It was subsequently found, by using a linear model, that there were an infinite number of solutions corresponding to various combinations of depth and lateral inhomogeneity. Linear modeling was performed by means of generalized inverse theory.

Phillips, R. J.↗

Digital adaptive flight controller development

A design study of adaptive control logic suitable for implementation in modern airborne digital flight computers was conducted. Two designs are described for an example aircraft. Each of these designs uses a weighted least squares procedure to identify parameters defining the dynamics of the aircraft. The two designs differ in the way in which control law parameters are determined. One uses the solution of an optimal linear regulator problem to determine these parameters while the other uses a procedure called single stage optimization. Extensive simulation results and analysis leading to the designs are presented.

Kaufman, H.↗

A decentralized square root information filter/smoother

A number of developments has recently led to a considerable interest in the decentralization of linear least squares estimators. The developments are partly related to the impending emergence of VLSI technology, the realization of parallel processing, and the need for algorithmic ways to speed the solution of dynamically decoupled, high dimensional estimation problems. A new method is presented for combining Square Root Information Filters (SRIF) estimates obtained from independent data sets. The new method involves an orthogonal transformation, and an information matrix filter 'homework' problem discussed by Schweppe (1973) is generalized. The employed SRIF orthogonal transformation methodology has been described by Bierman (1977).

Bierman, G. J.↗

Polar decomposition for attitude determination from vector observations

This work treats the problem of weighted least squares fitting of a 3D Euclidean-coordinate transformation matrix to a set of unit vectors measured in the reference and transformed coordinates. A closed-form analytic solution to the problem is re-derived. The fact that the solution is the closest orthogonal matrix to some matrix defined on the measured vectors and their weights is clearly demonstrated. Several known algorithms for computing the analytic closed form solution are considered. An algorithm is discussed which is based on the polar decomposition of matrices into the closest unitary matrix to the decomposed matrix and a Hermitian matrix. A somewhat longer improved algorithm is suggested too. A comparison of several algorithms is carried out using simulated data as well as real data from the Upper Atmosphere Research Satellite. The comparison is based on accuracy and time consumption. It is concluded that the algorithms based on polar decomposition yield a simple although somewhat less accurate solution. The precision of the latter algorithms increase with the number of the measured vectors and with the accuracy of their measurement.

Bar-Itzhack, Itzhack Y.↗

A solution scheme for the Euler equations based on a multi-dimensional wave model

A scheme for the solution of scalar advection on an unstructured mesh has been developed, tested, and extended to the Euler equations. The scheme preserves a linear function exactly, and yields nearly monotone results. The flux function associated with the Euler scheme is based on a discrete 'wave model' for the system of equations. The wave model decomposes the solution gradient at a location into shear waves, entropy waves and acoustic waves and calculates the speeds, strengths and directions associated with the waves. The approach differs from typical flux-difference splitting schemes in that the waves are not assumed to propagate normal to the faces of the control volumes; directions of propagation of the waves are instead computed from solution-gradient information. Results are shown for three test cases, and two different wave models. The results are compared to those from other approaches, including MUSCL and Galerkin least squares schemes.

Powell, Kenneth G.↗

Development of MLPG and LBIE Methods for Nonlinear Problems of Fracture

The truly meshless local Petrov-Galerkin (MLPG) method holds a great promise in solving boundary value problems, using a local symmetric weak form as a natural approach. In the present paper, in the context of MLPG and the meshless interpolation of a moving least squares (MLS) type, a method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The sub-domains for the shape functions in the MLS approximation are defined only from the primary nodes, and the secondary nodes use the same sub-domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and/or moved without an additional mesh. The essential boundary conditions can be imposed exactly, and non-convex boundaries can be treated without special techniques. Several numerical examples are presented to illustrate the performance of the present method.

Atluri, Satya N.↗

Tracking and Data Relay Satellite (TDRS) tracking and orbit determination

Analysis results were derived using the Goddard Trajectory Determination System to determine the orbit of TDRS-East. A batch weighted least-squares algorithm is used to fit the Bilateration Ranging Transponder system tracking measurements. Six state parameters (position and velocity components) and the solar radiation pressure coefficient are estimated. The solutions are obtained by fitting the data over 34 hr arcs that overlap by 10 hr. Ephemerides are generated over the 34 hr arcs and then compared over the overlapping regions. The position differences indicate the consistency of the solutions and give a measure of their accuracy. The effect of the modeling of perturbative forces (solar radiation pressure, nonspherical geopotential field) on these results are studied. Error analysis is performed using the ORAN program to estimate the effect of force-model errors and measurement-related errors on the TDRS-East ephemeris.

Teles, J.↗

Improved Nystrom Integrators

New formulation for Nystrom integrators offers increased precision and speed of computation. Formulation takes advantage of capabilities of modern electronic computers. Independent and dependent variables and various partial and total derivatives expressed via compact tensor notation. In this approach, typical set of constraint equations derived in about 1 day. Solution of constraint equations speeded by SEARCH computer program. Is general-purpose optimization program, which maximizes or minimizes cost function defined by user. Used successfully, to make least-squares fit of about 4,000 equations in 12 unknowns.

Lear, William M.↗