Search NASA⌕ Search

SEARCH · Search NASA

Results for “Matrix equations”

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

A diagonally inverted LU implicit multigrid scheme

A new Diagonally Inverted LU Implicit scheme is developed within the framework of the multigrid method for the three-dimensional unsteady Euler equations. The matrix systems that are to be inverted in the LU scheme are treated by local diagonalizing transformations that decouple them into systems of scalar equations. Unlike the Diagonalized ADI method, the time accuracy of the LU scheme is not reduced since the diagonalizing procedure does not destroy time conservation. Even more importantly, this diagonalization significantly reduces the computational effort required to solve the LU approximation and therefore transforms it into a more efficient method of numerically solving the three-dimensional Euler equations.

Yokota, J. W.↗

Modal methods in optimal control synthesis

Efficient algorithms for solving linear smoother-follower problems with quadratic criteria are presented. For time-invariant systems, the algorithm consists of one backward integration of a linear vector equation and one forward integration of another linear vector equation. Furthermore, the backward and forward Riccati matrices can be expressed in terms of the eigenvalues and eigenvectors of the Euler-Lagrange equations. Hence, the gains of the forward and backward Kalman-Bucy filters and of the optimal state-feedback regulator can be determined without integration of matrix Riccati equations. A computer program has been developed, based on this method of determining the gains, to synthesize the optimal time-invariant compensator in the presence of random disturbance inputs and random measurement errors. The program also computes the rms state and control variables of the optimal closed-loop system.

Bryson, A. E., Jr.↗

Linearized ray-trace analysis

A new, coordinate-free version of the exact ray-trace equations for optical systems consisting of conic reflecting, refracting and reference surfaces is presented. These equations are differentiated to obtain closed-form optical sensitivity dyadics. For computation, the sensitivities are evaluated in a single global coordinate frame and combined in linearized ray-trace matrix difference equations that propagate the rays and the sensitivities from element to element. One purpose of this analysis is to create optical models that can be directly integrated with models of the instrument structure and control systems for dynamic simulation.

Redding, David C.↗

Large-angle slewing maneuvers for flexible spacecraft

A new class of closed-form solutions for finite-time linear-quadratic optimal control problems is presented. The solutions involve Potter's solution for the differential matrix Riccati equation, which assumes the form of a steady-state plus transient term. Illustrative examples are presented which show that the new solutions are more computationally efficient than alternative solutions based on the state transition matrix. As an application of the closed-form solutions, the neighboring extremal path problem is presented for a spacecraft retargeting maneuver where a perturbed plant with off-nominal boundary conditions now follows a neighboring optimal trajectory. The perturbation feedback approach is further applied to three-dimensional slewing maneuvers of large flexible spacecraft. For this problem, the nominal solution is the optimal three-dimensional rigid body slew. The perturbation feedback then limits the deviations from this nominal solution due to the flexible body effects. The use of frequency shaping in both the nominal and perturbation feedback formulations reduces the excitation of high-frequency unmodeled modes. A modified Kalman filter is presented for estimating the plant states.

Chun, Hon M.↗

A computer package for the design and eigenproblem solution of damped linear multidegree of freedom systems

Systems described by the matrix differental equation are considered. An interactive design routine is presented for positive definite mass, damping, and stiffness matrices. Designing is accomplished by adjusting the mass, damping, and stiffness matrices to obtain a desired oscillation behavior. The algorithm also features interactively modifying the physical structure of the system, obtaining the matrix structure and a number of other system properties. In case of a general system, where the M, C, and K matrices lack any special properties, a routine for the eigenproblem solution of the system is developed. The latent roots are obtained by computing the characteristic polynomial of the system and solving for its roots. The above routines are prepared in FORTRAN IV and prove to be usable for the machines with low core memory.

Ahmadian, M.↗

Optimal filtering in the presence of unmodeled time correlated driving disturbances.

In many realistic data filtering problems, the cross correlation of the state estimation error and the state forcing function is unknown due to the poor knowledge of the time history of the forcing function. In this paper, the conservative and minimal approximation to the cross correlation terms is presented. It requires only the knowledge of the estimation error covariance and the forcing function covariance, with the choice of an associated free parameter left to the user. If the estimation error covariance and/or the forcing function covariance are bounded from above but not known exactly, the cross correlation approximation using those upper bounds remains conservative. This cross correlation approximation leads to a conservative approximation to the estimation error covariance matrix differential equation between measurement times. The free parameter is determined as the analytic solution to an associated optimal control problem. The procedure is expanded to include discrete linear measurement incorporation.

Fraser, D. C.↗

Perturbation guidance for minimum time flight paths of spacecraft.

The problem of transferring a rocket vehicle from a given circular orbit to a larger coplanar circular orbit in minimum time, using a constant low-thrust rocket engine, is considered. Parameters are chosen to correspond to a transfer from the earth's orbit in heliocentric space to the orbit of Mars. A path satisfying the first order necessary conditions of variational calculus is shown to be locally minimizing by application of a set of second order conditions. A physical explanation is offered to justify the retrothrust period occurring during the flight. A neighboring optimum feedback control law, based on estimated time-to-go, is applied to this problem. State variable and terminal constraint feedback gains are calculated while one of the second order conditions, involving the backward integration of a matrix Riccati equation, is being tested.

Wood, L. J.↗

The optimal control of merging aircraft - Implementation of the hybrid air traffic controller.

The control of merging aircraft is formulated as a finite-time, quadratic optimal control problem of a linear system with state and control constraints. The purpose of this paper is to demonstrate that the Hybrid Air Traffic Controller (HAC), which has been previously developed as a solution to this problem, may be easily implemented. Use is made of both the properties of the algebraic solution to the matrix Riccati equation and the structure of the linear model. This approach results in a real-time synthesis procedure for the HAC which does not rely on iterative numerical integration techniques.

Schatz, J. G.↗

Application of modern control theory to the design of optimum aircraft controllers

The synthesis procedure presented is based on the solution of the output regulator problem of linear optimal control theory for time-invariant systems. By this technique, solution of the matrix Riccati equation leads to a constant linear feedback control law for an output regulator which will maintain a plant in a particular equilibrium condition in the presence of impulse disturbances. Two simple algorithms are presented that can be used in an automatic synthesis procedure for the design of maneuverable output regulators requiring only selected state variables for feedback. The first algorithm is for the construction of optimal feedforward control laws that can be superimposed upon a Kalman output regulator and that will drive the output of a plant to a desired constant value on command. The second algorithm is for the construction of optimal Luenberger observers that can be used to obtain feedback control laws for the output regulator requiring measurement of only part of the state vector. This algorithm constructs observers which have minimum response time under the constraint that the magnitude of the gains in the observer filter be less than some arbitrary limit.

Power, L. J.↗

Evaluation of glide paths for landing a VTOL airplane using linear regulator theory.

A method of evaluating certain characteristics of approach paths for VTOL airplanes is presented which is based on the solution of the matrix Riccati equation to obtain an optimal state variable feedback controller. The longitudinal equations of motion of the airplane are linearized about a preselected path and the resulting system of equations is treated as a linear, time-varying regulator. The controller which minimizes a quadratic cost function is applied to the linearized system to determine the airplane's ability to return to the prescribed path given a specified initial error in altitude. The procedure is applied to the XC-142A, tilt-wing, V/STOL airplane, under decelerating approach conditions with a glide path consisting of two segments, the first having a smaller angle of descent than the second.

Reid, G. F.↗

Propagation of plane waves in a variable area duct carrying a compressible subsonic flow

A complete numerical solution is obtained for plane wave propagation in a variable-area duct which carries a mean compressible flow. The equation solved is a matrix formulation of the equations first derived by Tsien in 1952. Sample calculations are presented for four particular flow configurations: (1) a duct-nozzle-duct combination with the acoustic wave traveling with the flow, (2) a duct-nozzle-duct combination with the acoustic wave traveling against the flow, (3) an inlet-nozzle-duct configuration with the acoustic wave traveling against the flow, and (4) a duct-nozzle-exit configuration with the acoustic wave traveling with the flow. For each of these cases, the variation of energy flux with frequency is shown for particular geometry and flow condition.

Davis, S. S.↗

Ocean color spectra measured off the Oregon coast - Characteristic vectors

The ocean color spectrum is defined as the ratio of the spectrum of light upwelled from the sea to the spectrum of light incident on the sea surface. Ocean color spectra, observed from an airplane flown over waters off Oregon, are analyzed. The original spectra are resolved into fifty-five wavelength bands, each 5 nm wide. The shapes of these spectra are parameterized by, and shown to be accurately recoverable from, their first four principal components. These components are the scalar projections of each spectrum on the first four characteristic vectors of the sample covariance matrix. Regression equations are found with which phytoplankton pigment concentration and water transparency may be estimated as linear functions of the principal components. Pigment concentration estimates thus obtained are imprecise. The poor fit is at least partly due to the inappropriateness of the linear regression model and the neglect of other optically important substances typically present in sea water.

Mueller, J. L.↗

On reliable control system designs

A mathematical model for use in the design of reliable multivariable control systems is discussed with special emphasis on actuator failures and necessary actuator redundancy levels. The model consists of a linear time invariant discrete time dynamical system. Configuration changes in the system dynamics are governed by a Markov chain that includes transition probabilities from one configuration state to another. The performance index is a standard quadratic cost functional, over an infinite time interval. The actual system configuration can be deduced with a one step delay. The calculation of the optimal control law requires the solution of a set of highly coupled Riccati-like matrix difference equations. Results can be used for off-line studies relating the open loop dynamics, required performance, actuator mean time to failure, and functional or identical actuator redundancy, with and without feedback gain reconfiguration strategies.

Birdwell, J. D.↗