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 235 records · Page 13

Evaluation of spacecraft navigation using the Tracking and Data Relay Satellite System (TDRSS)

The navigational capabilities of the TDRSS for the operational support of user spacecraft is evaluated. A batch-weighted least-squares algorithm is used to fit tracking measurements and to generate orbit solutions for the TDRS-E spacecraft and the following user spacecraft: (1) the SMM, (2) the Landsat-5, (3) the Earth Radiation Budget Satellite (ERBS), and (4) the Solar Mesosphere Explorer (SME). TDRS-E orbit accuracy was evaluated by using the consistency of consecutive orbit solutions as a measure; user orbit accuracy was estimated by comparing orbits determined from TDRSS tracking with orbits obtained from ground tracking for the same time periods. It was found that the predicted accuracy of the TDRS-E orbit propagation was improved by as much as seven times by refining the spherical harmonic expansion model of the earth's potential field. The results obtained for SMM were fairly consistent, while those obtained for Landsat-5 were consistent to better than 80 meters. The orbits obtained using ERBS tracking were consistent to better than 60 meters and the SME orbits for two arcs with 14 and 20 passes of TDRSS tracking had a maximum position difference of 140 meters.

Samii, M. V.↗

Optimal model reduction and frequency-weighted extension

In this paper the quadratically optimal model reduction problem for single-input, single-output systems is considered. The reduced order model is determined by minimizing the integral of the magnitude-squared of the transfer function error. It is shown that the numerator coefficients of the optimal approximant satisfy a weighted least squares problem and, on this basis, a two-step iterative algorithm is developed combining a least squares solver with a gradient minimizer. The existence of globally optimal stable solutions to the optimization problem is established, and convergence of the algorithm to stationary values of the cost function is proved. The formulation is extended to handle the frequency-weighted optimal model reduction problem. Three examples demonstrate the optimization algorithm.

Spanos, J. T.↗

An improved solution to the 'Rankine-Hugoniot' problem

This paper presents an extension of the nonlinear least squares fitting technique of Vinas and Scudder (1986) (VS), which finds the physical and geometrical properties of nondissipational magnetohydrodynamic (MHD) shocks. The new method incorporates plasma temperature observations in the form of normal momentum flux and energy density flux conservation as well as plasma density, velocity, and magnetic field data. The new technique is capable of using known standard deviations in the individual measurement points to properly weight the fitting procedure. The new fitting code is validated through the analysis of synthetic shocks with known physical and geometrical properties. Finally, it is compared to the original VS method and the preaveraged velocity coplanarity technique.

Szabo, A.↗

A Local Coordinate Approach in the MLPG Method for Beam Problems

System matrices for Euler-Bernoulli beam problems for the meshless local Petrov-Galerkin (MLPG) method deteriorate as the number of nodes in the beam models are consistently increased. The reason for this behavior is explained. To overcome this difficulty and improve the accuracy of the solutions, a local coordinate approach for the evaluation of the generalized moving least squares shape functions and their derivatives is proposed. The proposed approach retains the accuracy of the MLPG methods.

Raju, Ivatury S.↗

Long-period perturbations in Starlette orbit and tide solution

The long-period perturbations present in the orbit of the Starlette geodetic satellite is analyzed using a one-year continuous orbit from October 1976 through October 1977 and a three-year orbit during 1983 through 1985. A solution for 66 ocean tidal parameters for 14 low-degree tidal constituents is obtained in a simultaneous least squares estimation procedure with other geodynamic and orbit parameters. The yearly values of the second-degree S(a) tides, which have been observed to produce significant seasonal variations on the Starlette node, are included in the solution. Good agreement is found between the solutions, other satellite solutions, the oceanographic tide solutions, and other investigations using the computation of the earth's tidal response. The results indicate that the Starlette tide solutions produces significantly improved orbit fits for both Lageos and Starlette.

Cheng, M. K.↗

Acoustic transmission in non-uniform ducts with mean flow. I - The method of weighted residuals. II - The finite element method

The problem of acoustic transmission through nonuniform ducts containing a high-speed subsonic flow is studied by means of the method of weighted residuals in the form of a modified Galerkin method and a Galerkin formulation of the finite element method. The method of weighted residuals is shown to employ the basis functions generated from eigenvalue calculations for the case of no flow, and is verified by comparison with exact eigenvalue calculations in the uniform duct case and numerical solutions of the one-dimensional form of the equations in the nonuniform duct case. The finite element scheme based on both the Galerkin method and the residual least squares method and employing eight-noded isoparametric elements is presented and used to investigate multimodal propagation by the coupling of the solution in the duct nonuniform section to modal expansions in uniform sections. Comparison of the results of the two methods reveals them to be in substantial agreement, and predicts the importance of multimodal interactions at high Mach numbers.

Eversman, W.↗

A least squares finite element scheme for transonic flow around harmonically oscillating airfoils

The present investigation shows that a finite element scheme with a weighted least squares variational principle is applicable to the problem of transonic flow around a harmonically oscillating airfoil. For the flat plate case, numerical results compare favorably with the exact solution. The obtained numerical results for the transonic problem, for which an exact solution is not known, have the characteristics of known experimental results. It is demonstrated that the performance of the employed numerical method is independent of equation type (elliptic or hyperbolic) and frequency. The weighted least squares principle allows the appropriate modeling of singularities, which such a modeling of singularities is not possible with normal least squares.

Cox, C. L.↗

Multistep modeling of protein structure: application towards refinement of tyr-tRNA synthetase

The scope of multistep modeling (MSM) is expanding by adding a least-squares minimization step in the procedure to fit backbone reconstruction consistent with a set of C-alpha coordinates. The analytical solution of Phi and Psi angles, that fits a C-alpha x-ray coordinate is used for tyr-tRNA synthetase. Phi and Psi angles for the region where the above mentioned method fails, are obtained by minimizing the difference in C-alpha distances between the computed model and the crystal structure in a least-squares sense. We present a stepwise application of this part of MSM to the determination of the complete backbone geometry of the 321 N terminal residues of tyrosine tRNA synthetase to a root mean square deviation of 0.47 angstroms from the crystallographic C-alpha coordinates.

NASA Discipline Exobiology↗

Orbit determination accuracies using satellite-to-satellite tracking

The uncertainty in relay satellite sate is a significant error source which cannot be ignored in the reduction of satellite-to-satellite tracking data. Based on simulations and real data reductions, it is numerically impractical to use simultaneous unconstrained solutions to determine both relay and user satellite epoch states. A Bayesian or least squares estimation technique with an a priori procedure is presented which permits the adjustment of relay satellite epoch state in the reduction of satellite-to-satellite tracking data without the numerical difficulties introduced by an ill-conditioned normal matrix.

Vonbun, F. O.↗

Efficient Implementation of Minimal Polynomial and Reduced Rank Extrapolation Methods

The minimal polynomial extrapolation (MPE) and reduced rank extrapolation (RRE) are two effective techniques that have been used in accelerating the convergence of vector sequences, such as those that are obtained from iterative solution of linear and nonlinear systems of equation. Their definitions involve some linear least squares problems, and this causes difficulties in their numerical implementation. Timewise efficient and numerically stable implementations for MPE and RRE are developed. A computer program written in FORTRAN 77 is also appended and applied to some model problems.

Sidi, Avram↗

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.↗

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.↗

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.↗

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.↗

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.↗

Stochastic average model methods

We consider the solution of finite-sum minimization problems, such as those appearing in nonlinear least-squares or general empirical risk minimization problems. We are motivated by problems in which the summand functions are computationally expensive and evaluating all summands on every iteration of an optimization method may be undesirable. Here we present the idea of stochastic average model (SAM) methods, inspired by stochastic average gradient methods. SAM methods sample component functions on each iteration of a trust-region method according to a discrete probability distribution on component functions; the distribution is designed to minimize an upper bound on the variance of the resulting stochastic model. We present promising numerical results concerning an implemented variant extending the derivative-free model-based trust-region solver POUNDERS, which we name SAM-POUNDERS.

97 MATHEMATICS AND COMPUTING↗