Search NASA⌕ Search

SEARCH · Search NASA

Results for “extended least squares”

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 109 records · Page 6

Kalman Filtering USNO's GPS Observations for Improved Time Transfer Predictions

The Global Positioning System (GPS) Master Control Station (MCS) performs the Coordinated Universal Time (UTC) time transfer mission by uploading and broadcasting predictions of the GPS-UTC offset in subframe 4 of the GS navigation message. These predictions are based on only two successive daily data points obtained from the US Naval Observatory (USNO). USNO produces these daily smoothed data points by performing a least-squares fit on roughly 38 hours worth of data from roughly 160 successive 13-minute tracks of GPS satellites. Though sufficient for helping to maintain a time transfer error specification of 28 ns (1 Sigma), the MCS's prediction algorithm does not make the best use of the available data from from USNO, and produces data that can degrade quickly over extended prediction spans. This paper investigates how, by applying Kalman filtering to the same available tracking data, the MCS could improve its estimate of GPS-UTC, and in particular, the GPS-UTC A(sub 1) term. By refining the A(sub 1) (frequency) estimate for GPS-UTC predictions, error in GPS time transfer could drop significantly. Additional, the risk of future spikes in GPS's time transfer error could similarly be minimized, by employing robust Kalman filtering for GPS-UTC predictions.

Hutsell, Steven T.↗

Aspects of the 40-50 day oscillation during the northern winter as inferred from outgoing longwave radiation

One-point correlation statistics and extended empirical orthogonal function analysis are used to show the existence of a principal mode of oscillation in the intraseasonal variability (ISV) of tropical convection. The data analyzed were derived from NOAA Outgoing Longwave Radiation data covering seven 6 mo periods from 1974-83. The data is filtered to remove fluctuations on time scales under 5 days and by subtracting a least-squares parabola to eliminate seasonal means. The previously-known 40-50 day oscillation is found to be the most prevalent component of the ISV of tropical convection. An E-W dipole forms and propagates over the equatorial Indian Ocean, then stalls and intensifies over Indonesia and the equatorial Indian Ocean, then stalls and intensifies over Indonesia and the equatorial central Pacific. Various anomalies generated in association with the dipole and dependent on its exact location when stalled are discussed.

Lau, K.-M.↗

Locally-Based Kernal PLS Smoothing to Non-Parametric Regression Curve Fitting

We present a novel smoothing approach to non-parametric regression curve fitting. This is based on kernel partial least squares (PLS) regression in reproducing kernel Hilbert space. It is our concern to apply the methodology for smoothing experimental data where some level of knowledge about the approximate shape, local inhomogeneities or points where the desired function changes its curvature is known a priori or can be derived based on the observed noisy data. We propose locally-based kernel PLS regression that extends the previous kernel PLS methodology by incorporating this knowledge. We compare our approach with existing smoothing splines, hybrid adaptive splines and wavelet shrinkage techniques on two generated data sets.

Rosipal, Roman↗

Automation of orbit determination functions for National Aeronautics and Space Administration (NASA)-supported satellite missions

The Flight Dynamics Facility (FDF) at Goddard Space Flight Center (GSFC) provides spacecraft trajectory determination for a wide variety of National Aeronautics and Space Administration (NASA)-supported satellite missions, using the Tracking Data Relay Satellite System (TDRSS) and Ground Spaceflight and Tracking Data Network (GSTDN). To take advantage of computerized decision making processes that can be used in spacecraft navigation, the Orbit Determination Automation System (ODAS) was designed, developed, and implemented as a prototype system to automate orbit determination (OD) and orbit quality assurance (QA) functions performed by orbit operations. Based on a machine-resident generic schedule and predetermined mission-dependent QA criteria, ODAS autonomously activates an interface with the existing trajectory determination system using a batch least-squares differential correction algorithm to perform the basic OD functions. The computational parameters determined during the OD are processed to make computerized decisions regarding QA, and a controlled recovery process isactivated when the criteria are not satisfied. The complete cycle is autonomous and continuous. ODAS was extensively tested for performance under conditions resembling actual operational conditions and found to be effective and reliable for extended autonomous OD. Details of the system structure and function are discussed, and test results are presented.

Mardirossian, H.↗

Automation of orbit determination functions for National Aeronautics and Space Administration (NASA)-supported satellite missions

The Flight Dynamics Facility (FDF) at Goddard Space Flight Center (GSFC) provides spacecraft trajectory determination for a wide variety of National Aeronautics and Space Administration (NASA)-supported satellite missions, using the Tracking Data Relay Satellite System (TDRSS) and Ground Spaceflight and Tracking Data Network (GSTDN). To take advantage of computerized decision making processes that can be used in spacecraft navigation, the Orbit Determination Automation System (ODAS) was designed, developed, and implemented as a prototype system to automate orbit determination (OD) and orbit quality assurance (QA) functions performed by orbit operations. Based on a machine-resident generic schedule and predetermined mission-dependent QA criteria, ODAS autonomously activates an interface with the existing trajectory determination system using a batch least-squares differential correction algorithm to perform the basic OD functions. The computational parameters determined during the OD are processed to make computerized decisions regarding QA, and a controlled recovery process is activated when the criteria are not satisfied. The complete cycle is autonomous and continuous. ODAS was extensively tested for performance under conditions resembling actual operational conditions and found to be effective and reliable for extended autonomous OD. Details of the system structure and function are discussed, and test results are presented.

Mardirossian, H.↗

Aperture codes for sensors viewing extended objects from space

The paper describes aperture codes viewing extended objects form space which find application in imaging an extended object which may have relatively low contrast, and whose lateral limits extend beyond the FOV of the sensor. Those elements of an extended object lying near the FOV limits are only partially coded, i.e., flux from those elements cannot cast a shadow of the entire aperture code onto the detector, as can elements near the center of the FOV. This has consequences for the algorithms used to reconstruct the image. The object field is divided into a number of elements which is smaller than the number of detector pixels, and a least squares fit to the data is performed. The methods used for choosing the matrices representing the aperture codes are discussed, and computer simulations of the effects of noise are described.

Curtis, C. C.↗

Maximizing TDRS Command Load Lifetime

The GNC software onboard ISS utilizes TORS command loads, and a simplistic model of TORS orbital motion to generate onboard TORS state vectors. Each TORS command load contains five "invariant" orbital elements which serve as inputs to the onboard propagation algorithm. These elements include semi-major axis, inclination, time of last ascending node crossing, right ascension of ascending node, and mean motion. Running parallel to the onboard software is the TORS Command Builder Tool application, located in the JSC Mission Control Center. The TORS Command Builder Tool is responsible for building the TORS command loads using a ground TORS state vector, mirroring the onboard propagation algorithm, and assessing the fidelity of current TORS command loads onboard ISS. The tool works by extracting a ground state vector at a given time from a current TORS ephemeris, and then calculating the corresponding "onboard" TORS state vector at the same time using the current onboard TORS command load. The tool then performs a comparison between these two vectors and displays the relative differences in the command builder tool GUI. If the RSS position difference between these two vectors exceeds the tolerable lim its, a new command load is built using the ground state vector and uplinked to ISS. A command load's lifetime is therefore defined as the time from when a command load is built to the time the RSS position difference exceeds the tolerable limit. From the outset of TORS command load operations (STS-98), command load lifetime was limited to approximately one week due to the simplicity of both the onboard propagation algorithm, and the algorithm used by the command builder tool to generate the invariant orbital elements. It was soon desired to extend command load lifetime in order to minimize potential risk due to frequent ISS commanding. Initial studies indicated that command load lifetime was most sensitive to changes in mean motion. Finding a suitable value for mean motion was therefore the key to achieving this goal. This goal was eventually realized through development of an Excel spreadsheet tool called EMMIE (Excel Mean Motion Interactive Estimation). EMMIE utilizes ground ephemeris nodal data to perform a least-squares fit to inferred mean anomaly as a function of time, thus generating an initial estimate for mean motion. This mean motion in turn drives a plot of estimated downtrack position difference versus time. The user can then manually iterate the mean motion, and determine an optimal value that will maximize command load lifetime. Once this optimal value is determined, the mean motion initially calculated by the command builder tool is overwritten with the new optimal value, and the command load is built for uplink to ISS. EMMIE also provides the capability for command load lifetime to be tracked through multiple TORS ephemeris updates. Using EMMIE, TORS command load lifetimes of approximately 30 days have been achieved.

Brown, Aaron J.↗

Signal Prediction With Input Identification

A novel coding technique is presented for signal prediction with applications including speech coding, system identification, and estimation of input excitation. The approach is based on the blind equalization method for speech signal processing in conjunction with the geometric subspace projection theory to formulate the basic prediction equation. The speech-coding problem is often divided into two parts, a linear prediction model and excitation input. The parameter coefficients of the linear predictor and the input excitation are solved simultaneously and recursively by a conventional recursive least-squares algorithm. The excitation input is computed by coding all possible outcomes into a binary codebook. The coefficients of the linear predictor and excitation, and the index of the codebook can then be used to represent the signal. In addition, a variable-frame concept is proposed to block the same excitation signal in sequence in order to reduce the storage size and increase the transmission rate. The results of this work can be easily extended to the problem of disturbance identification. The basic principles are outlined in this report and differences from other existing methods are discussed. Simulations are included to demonstrate the proposed method.

Juang, Jer-Nan↗

STS-8 bet results

The final Best Estimate Trajectory (BET) products, i.e., the reconstructed trajectory, the Extended BET, AEROBET and MMLE input files, generated for the eighth NASA Space Shuttle flight are documented. The reconstructed trajectory (inertial BET) for this Challenger flight, the first night landing is discussed. State (position, velocity, and attitude) plus three accelerometer scale factors were determined from fitting the Guam S-band data, seven C-band passes, and pseudo Doppler and altimeter during rollout on Runway 22. The anchor epoch utilized for the batch weighted-least-squares determination was Sept. 5, 1983 7h1m50s.0 (25310 GMT seconds). The spacecraft altitude at epoch is approx. 617 kft. IMU2 data were selected for the reconstruction.

Findlay, J. T.↗

Magnetic field structure of interplanetary magnetic clouds at 1 AU

Interplanetary magnetic clouds emerge as a feature of the solar wind at 1 AU, exhibiting enhanced field strength and lower plasma temperature and density than the surrounding plasma. A least-squares program has been developed which fits magnetic field data within a cloud, while estimating such cloud properties as its size, maximum field strength, and axis inclination. The results obtained from a study of 12 clouds observed at 1 AU point to a probable cloud axis direction within 15 deg of the ecliptic plane and about 100 deg from the sun's direction, when projected into the ecliptic plane. A wide variety of orientations is observed; some extend to 80 deg from the ecliptic.

Lepping, R. P.↗

East-west inclination of large-scale photospheric magnetic fields

Sixteen years of Wilcox Solar Observatory (WSO) magnetogram data have been studied to determine the solar cycle variation and latitude dependence of the east-west inclination of photospheric magnetic field lines. East-west inclination is here defined as the angle between a field line and its local radial vector, as projected onto the plane of the latitude and line of sight. Inclination is determined by a least-squares fit of observed magnetic fields to a simple projection model, and is found to depend on polarity and to change with the solar cycle. Leading and following polarities are tipped towards each by about 9 deg and have an overall net tilt in the direction of rotation (to the west) of 0.6 deg. New cycles are seen to begin at high latitudes and to grow through the lower latitudes over approximately 5 years, providing evidence for an extended cycle length of 16-18 years.

Shrauner, J. A↗

Analysis of altimetry over inland seas

Satellite-based altimetric data taken by GEOS-3 and SEASAT over the Black Sea and Caspian Sea are analyzed and a least squares collocation technique is used to predict the geoid undulation on a .25-degree by .25-degree grid and to transform these geoid undulations to free air gravity anomalies. This project entailed processing satellite altimeter data over inland seas for recovery of area mean gravity information. Gravity information in this area of the world is not readily available, so the possibility of obtaining it from the processing of altimeter observations is attractive. The principal objective was to complete and extend analyses done in a previous study, verify those results, and document the results and techniques. A secondary objective was to improve the algorithms and results, if possible. The approach used involved editing geoid height data to remove overland data; evaluating geoid height differences at crossover points; removing orbit errors from geoid heights using crossover differences; gridding geoid height data at .25-degree by .25-degree intervals; and estimating the gravity anomalies from gridded geoid heights using the collocation technique.

Au, A. Y.↗

Uncertainty Analysis of Inertial Model Attitude Sensor Calibration and Application with a Recommended New Calibration Method

Statistical tools, previously developed for nonlinear least-squares estimation of multivariate sensor calibration parameters and the associated calibration uncertainty analysis, have been applied to single- and multiple-axis inertial model attitude sensors used in wind tunnel testing to measure angle of attack and roll angle. The analysis provides confidence and prediction intervals of calibrated sensor measurement uncertainty as functions of applied input pitch and roll angles. A comparative performance study of various experimental designs for inertial sensor calibration is presented along with corroborating experimental data. The importance of replicated calibrations over extended time periods has been emphasized; replication provides independent estimates of calibration precision and bias uncertainties, statistical tests for calibration or modeling bias uncertainty, and statistical tests for sensor parameter drift over time. A set of recommendations for a new standardized model attitude sensor calibration method and usage procedures is included. The statistical information provided by these procedures is necessary for the uncertainty analysis of aerospace test results now required by users of industrial wind tunnel test facilities.

Tripp, John S.↗

Standard remote manipulator system docking target augmentation for automated docking

A docking target is provided for use in automated docking of a first vehicle on which the target is located. The target comprises a pair of laterally extending arm portions lying in substantially the same plane and a central post extending outwardly from the plane of the arm portions. At least three reflectors are located on the target. Two of the reflectors are located at the outboard ends of the arms portions and another reflector is located at the end of the central post. In an important embodiment, the reflectors comprise individual pieces of retroreflective tape. The reflectors, when viewed from the front of the target, are aligned along the longitudinal center line of the target, and can take a number of different shapes including circular or square.

Dabney, Richard W.↗

Chemical kinetic studies of atmospheric reactions using tunable diode laser spectroscopy

IR absorption using tunable diode laser spectroscopy provides a sensitive and quantitative detection method for laboratory kinetic studies of atmospheric trace gases. Improvements in multipass cell design, real time signal processing, and computer controlled data acquisition and analysis have extended the applicability of the technique. We have developed several optical systems using off-axis resonator mirror designs which maximize path length while minimizing both the sample volume and the interference fringes inherent in conventional 'White' cells. Computerized signal processing using rapid scan (300 kHz), sweep integration with 100 percent duty cycle allows substantial noise reduction while retaining the advantages of using direct absorption for absolute absorbance measurements and simultaneous detection of multiple species. Peak heights and areas are determined by curve fitting using nonlinear least square methods. We have applied these techniques to measurements of: (1) heterogeneous uptake chemistry of atmospheric trace gases (HCl, H2O2, and N2O5) on aqueous and sulfuric acid droplets; (2) vapor pressure measurements of nitric acid and water over prototypical stratospheric aerosol (nitric acid trihydrate) surfaces; and (3) discharge flow tube kinetic studies of the HO2 radical using isotopic labeling for product channel and mechanistic analysis. Results from each of these areas demonstrate the versatility of TDL absorption spectroscopy for atmospheric chemistry applications.

Worsnop, Douglas R.↗

Towards a Third-Order Accurate, Second-Derivative-Free, Shock-Capturing Finite-Volume Method for Hypersonic Flows on Tetrahedral Grids

In this paper, we report progress in the development of a third-order accurate, second-derivative-free, shock-capturing finite-volume solver for three-dimensional unstructured grids. The method is economical in the sense that the computation and storage of second derivatives are not required for third-order accuracy. It is based on point-valued numerical solutions stored at cells, gradients computed and stored at nodes, and an efficient projected-derivative formula that eliminates the need for second derivatives in a quadratic solution interpolation. The projected-derivative formula is also used to eliminate second derivatives from a high-order flux quadrature formula, so that it can be implemented conveniently in the form of a numerical flux at a face center plus a correction term. Similarly, a high-order source quadrature formula can also be implemented in the form of a cell-center point evaluation plus a similar correction term. These features make it relatively straightforward to extend an existing second-order finite-volume code to third-order. This paper reports progress of implementing the method in the NASA VULCAN-CFD code and discusses the implementation of a high-order accurate limiter for shock capturing.

Weighted Least-Squares↗

Radar altimeter waveform modeled parameter recovery

Satellite-borne radar altimeters include waveform sampling gates providing point samples of the transmitted radar pulse after its scattering from the ocean's surface. Averages of the waveform sampler data can be fitted by varying parameters in a model mean return waveform. The theoretical waveform model used is described as well as a general iterative nonlinear least squares procedures used to obtain estimates of parameters characterizing the modeled waveform for SEASAT-1 data. The six waveform parameters recovered by the fitting procedure are: (1) amplitude; (2) time origin, or track point; (3) ocean surface rms roughness; (4) noise baseline; (5) ocean surface skewness; and (6) altitude or off-nadir angle. Additional practical processing considerations are addressed and FORTRAN source listing for subroutines used in the waveform fitting are included. While the description is for the Seasat-1 altimeter waveform data analysis, the work can easily be generalized and extended to other radar altimeter systems.

Source record↗

Wind-tunnel investigation of six shielded total-pressure tubes at high angles of attack subsonic speeds

The effect of inclination of the air stream on the measured pressures of six shielded total-pressure tubes (all designed for end-mounting on a horizontal boom) has been determined for an angle-of-attack range of approximately -30° to 65°. The tests were conducted at a Mach number of 0.26 in the Langley stability tunnel and at Mach numbers of 0.50, 0.70, 0.90, and 0.95 in the Langley high-speed 7- by 10-foot tunnel. The results of the tests showed that curved venturi entries are less sensitive to inclination of the air stream than the conical entry of the standard Kiel design. Of the three curved entries tested, the tube having the bluntest entry proved the least sensitive to inclination. The critical angle (that is, the sngle.at which the total-pressure error equals 1 percent of the impact pressure) for this tube was ±63.2° at a Mach number of 0.26. The tests also showed that the critical angle of a shielded tube could be extended to higher positive angles of attack by means of a slant profile. The critical angle of a tube having a conical entry with a 10° slant profile, for example, was found to average 9° greater than that of a similar tube having a square profile. Although the critical angle of this tube at angles of yaw was less than that at angles of attack, the critical angle in yaw was equal to that of the tube with the square profile. Tests of the effect of varying the position of the total-pressure probe in the shield showed that, in general, no advantage was to be gained by changing from the position used in the standard Kiel design. Limited tests of the effect of varying the internal diameter of the shield showed that an increase in the diameter of the throat corresponding an increase in the internal area of 50 percent resulted in a relatively small increase (2.5°) in the critical angle. The effect of Mach number on the sensitivity of the tubes to inclination was shown to be appreciable. An increase in Mach number from 0.26 to 0.90, for example, resulted in a decrease in the critical angle of most of the tubes of about 4°.

Walter R Russell↗