Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “LEAST SQUARE FIT”

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 721 records · Page 40

Empirical source noise prediction method with application to subsonic coaxial jet mixing noise

A general empirical method, developed for source noise predictions, uses tensor splines to represent the dependence of the acoustic field on frequency and direction and Taylor's series to represent the dependence on source state parameters. The method is applied to prediction of mixing noise from subsonic circular and coaxial jets. A noise data base of 1/3-octave-band sound pressure levels (SPL's) from 540 tests was gathered from three countries: United States, United Kingdom, and France. The SPL's depend on seven variables: frequency, polar direction angle, and five source state parameters: inner and outer nozzle pressure ratios, inner and outer stream total temperatures, and nozzle area ratio. A least-squares seven-dimensional curve fit defines a table of constants which is used for the prediction method. The resulting prediction has a mean error of 0 dB and a standard deviation of 1.2 dB. The prediction method is used to search for a coaxial jet which has the greatest coaxial noise benefit as compared with an equivalent single jet. It is found that benefits of about 6 dB are possible.

Zorumski, W. E.↗

Analysis of 12-700 micron reflectance spectra of three optical black samples

Normal incidence, specular reflectance spectra of three samples of optical black coatings have been measured in 10 filter passbands at effective wavelengths between 12 and 700 microns. Strong absorption, due to amorphous silicate material in the coatings, is found in the 15-50 micron region of each spectrum. Reflecting-layer model calculations are fitted by a nonlinear least-squares routine to the spectra at wavelengths beyond the absorption band to determine far-infrared optical parameters. It is found that, with a similar wavelength dependence, the absorption coefficient of a new black coating (ECP 2200) made by the Minnesota Mining and Manufacturing Company (3M) is less than 1/10 that of the discontinued coating, 3M Black Velvet Nextel. It is also shown that the extinction effect of a very rough substrate is negligible compared with the effects of increased upper surface roughness and coating thickness attendant upon fully covering a very rough substrate.

Smith, S. M.↗

Simulation and study of small numbers of random events

Random events were simulated by computer and subjected to various statistical methods to extract important parameters. Various forms of curve fitting were explored, such as least squares, least distance from a line, maximum likelihood. Problems considered were dead time, exponential decay, and spectrum extraction from cosmic ray data using binned data and data from individual events. Computer programs, mostly of an iterative nature, were developed to do these simulations and extractions and are partially listed as appendices. The mathematical basis for the compuer programs is given.

Shelton, R. D.↗

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

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

Absorption intensities and complex refractive indices of crystalline HCN, HC3N, and C4N2 in the infrared region

IR absorption intensities are presented for thin crystalline films of HCN, HC3N, and C4N2, together with n and k complex refractive indices determined on the basis of an iterative program for the Kramers-Konig integral via a least-squares, point-by-point fitting of the experimental transmission data. It is established that the transmission spectra generated by means of these n and k values can reproduce the experimental transmission observation values to within + or - 2 percent.

Masterson, C. M.↗

Photometry of the 4686 A emission line of gamma(2) Velorum from the South Pole

An automated optical telescope located at the Amundsen-Scott South Pole station on Antarctica, has been used to obtain more than 78 h of photometry of the He II emission line (4686 A) of the spectroscopic binary gamma(2) Velorum. These data were obtained on seven different days during the 1987 austral winter; the longest continuous run spans 19 h. Two independent period search techniques have been used to search for periodic behavior in the strength of the He II emission line of this Wolf-Rayet star. They are: (1) power spectrum analysis and (2) a first-order sine function fit to the data using least squares. Various multiplicities of a period on the order of 1.3 h with amplitudes of a few percent are found in most of these data. According to recent theoretical models of Wolf-Rayet stars, fluctuations in the He II emission line may indicate vibrational instability in gamma(2) Vel. These pulsations may, in turn, give rise to shocks which propagate outward and which may provide the necessary conditions for periodic changes in the state of a given region of the atmosphere to occur.

Taylor, Maryjane↗

Venusian highlands - Geoid to topography ratios and their implications

Geoid-to-topography ratios (GTRs) are estimated for 12 Venusian highland features to allow comparison with convection calculations and with terrestrial data of oceanic hot spots, swells, and plateaus. The geoid is estimated in the wavenumber domain from the isostatic formula, using the topography and the apparent depths isostatic compensation (ADC) for each region. In the space domain, the GTR is equal to the least squares slope of the linear fit of the geoid to the topograpy. The resulting GTR range is 7-31 m/km, which is much higher than terrestrial oceanic values (-1 to 5 m/km). The features fall into two distinct groups, one with a GTR range of 7-13 m/km, and one with a range of 19-25 m/km. A model for thermal thinning of a 100-km-thick lithosphere fits all values in the lower GTR group to within one standard deviation.

Smrekar, Suzanne E.↗

A combined surface/volume scattering retracking algorithm for ice sheet satellite altimetry

An algorithm that is based upon a combined surface-volume scattering model is developed. It can be used to retrack individual altimeter waveforms over ice sheets. An iterative least-squares procedure is used to fit the combined model to the return waveforms. The retracking algorithm comprises two distinct sections. The first generates initial model parameter estimates from a filtered altimeter waveform. The second uses the initial estimates, the theoretical model, and the waveform data to generate corrected parameter estimates. This retracking algorithm can be used to assess the accuracy of elevations produced from current retracking algorithms when subsurface volume scattering is present. This is extremely important so that repeated altimeter elevation measurements can be used to accurately detect changes in the mass balance of the ice sheets. By analyzing the distribution of the model parameters over large portions of the ice sheet, regional and seasonal variations in the near-surface properties of the snowpack can be quantified.

Davis, Curt H.↗

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

A Method of Estimating Residuals in Orbital Theory

The degree of approximation used in determining the orbits of earth satellites is reflected in the residuals (differences between calculated and observed positions). The least-squares procedure generally used to fit theory to observation tends to obscure the significance of theoretical parameters, so that the physical sources of residuals cease to be apparent. A method is outlined herein for estimating the magnitude of the residuals to be expected from an approximate theory presumed to have one missing or incorrect term.

Lecar, Myron↗

Air- and Self-Broadened Half Widths, Pressure-Induced Shifts, and Line Mixing in the Nu(sub 2) Band of (12)CH4

Lorentz self- and air-broadened half width and pressure-induced shift coefficients and their dependences on temperature have been measured from laboratory absorption spectra for nearly 130 transitions in the nu(sub 2) band of (12)CH4. In addition line mixing coefficients (using the relaxation matrix element formalism) for both self- and airbroadening were experimentally determined for the first time for a small number of transitions in this band. Accurate line positions and absolute line intensities were also determined. These parameters were obtained by analyzing high-resolution (approx. 0.003 to 0.01 per cm) laboratory spectra of high-purity natural CH4 and air-broadened CH4 recorded at temperatures between 226 and 297 K using the McMath-Pierce Fourier transform spectrometer (FTS) located at the National Solar Observatory on Kitt Peak, Arizona. A multispectrum nonlinear least squares technique was used to fit short (5-15 per cm) spectral intervals in 24-29 spectra simultaneously. Parameters were determined for nu(sub 2) transitions up to J" = 16. The variations of the measured broadening and shift parameters with the rotational quantum number index and tetrahedral symmetry species are examined. The present results are also compared with previous measurements available in the literature.

Smith, M. A. H.↗

Self- and Air-Broadened Line Shapes in the 2v3 P and R Branches of 12CH4

In this paper we report line shape parameters of 12CH4 for several hundred 2V(sub 3) transitions in the spectral regions 5891-5996 cm( exp -1) (P branch) and 6015-6115 cm(exp -1) (R branch). Air- and self-broadening coefficients were measured as a function of temperature; line mixing via off-diagonal relaxation matrix element coefficients was also obtained for 47 transition pairs. In total, nearly 1517 positions and intensities were retrieved, but many transitions were too weak for the line shape study. For this analysis, we used 25 high-resolution (0.0056 and 0.0067 cm(ex[ -1) and high signal-to-noise (S/N) spectra of high-purity 12CH4 and the same high-purity 12CH4 broadened by dry air recorded at different sample temperatures between 130 K and 295 K with the Bruker IFS 125HR Fourier transform spectrometer at JPL. Three different absorption cells were used (1) a White cell set to a path length of 13.09 m for room temperature data, (2) a single-pass 0.2038 m long coolable cell (for self-broadening) and (3) a multipass cell with 20.941 m total path coolable Herriott cell (for air-broadening). In total there were 13 spectra with pure 12CH4 (0.27-599 Torr) and 12 air-broadened spectra with total sample pressures of 80-805 Torr and volume mixing ratios (VMR) of methane between 0.18 and 1.0. An interactive multispectrum nonlinear least-squares technique was employed to fit the individual P10-P1 and R0-R10 manifolds in all the spectra simultaneously. Results obtained from the present analysis are compared to other recent measurements.

Devi, V. Malathy↗

Model Choice Metrics to Optimize Profile-QSAR Performance

Predicting molecular activity against protein targets is difficult because of the paucity of experimental data. Approaches like multitask modeling and collaborative filtering seek to improve model accuracy by leveraging results from multiple targets, but are limited because different compounds are measured with different assays, leading to sparse data matrices. Profile-QSAR (pQSAR) 2.0 addresses this problem by fitting a series of partial least squares models for each target, using as features the predictions from single-task models on the remaining targets. Here, this method has been shown to produce better results than single task and multitask models. However, the factors determining the success of pQSAR 2.0 have as yet not been characterized. In this paper we examine the experimental conditions that lead to better pQSAR models. We limit the amount of data available to the method by retraining with decreasing amounts of data and explore the model’s ability to generalize to compounds that have never been assayed. Finally, we look at the properties of training data needed to demonstrate pQSAR improvement.

Biological and medical sciences, Computer science↗

Least-Squares Self-Calibration of Imaging Array Data

When arrays are used to collect multiple appropriately-dithered images of the same region of sky, the resulting data set can be calibrated using a least-squares minimization procedure that determines the optimal fit between the data and a model of that data. The model parameters include the desired sky intensities as well as instrument parameters such as pixel-to-pixel gains and offsets. The least-squares solution simultaneously provides the formal error estimates for the model parameters. With a suitable observing strategy, the need for separate calibration observations is reduced or eliminated. We show examples of this calibration technique applied to HST NICMOS observations of the Hubble Deep Fields and simulated SIRTF IRAC observations.

Arendt, R. G.↗

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic (LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖, 𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, there are some instances where the fitting needs to be done programmatically. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

PolyFit: A C++ code for Polynomial Curve Fit with Calculation of Error Bars

In radiobiology, many dose-response results are modeled using the so-called linear-quadratic(LQ)model, which means that results are modeled as a function of dose Das 𝑅(𝐷)=𝛽0+𝛽1𝐷+𝛽2𝐷2. The coefficients 𝛽0, 𝛽1and 𝛽2are obtained from fitting a series of data points (𝑥𝑖,𝑦𝑖), which is usually done using a least-square method. The LQ and more generally the polynomial fit capability is implemented in many software that analyzes data. However, it is often convenient to do the fitting programmatically, especially when a large number of datasets should be analyzed. Furthermore, depending on the software used, some features may not be implemented. In this mini-review, I discuss the basis of polynomial fitting, including the calculation of errors on the coefficients and results, use of weighting and fixing the intercept value (the coefficient 𝛽0). A simple C++ code to perform the polynomial curve fitting is also provided. This code should be useful not only in radiobiology but in other fields of science as well.

Ianik Plante↗

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↗