Some biomedical applications of a non-linear curve fit method
Nonlinear curve fit method for estimating parameters of equations for biological models
SEARCH · Search NASA
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.
Nonlinear curve fit method for estimating parameters of equations for biological models
Boundary-layer transition data for angles of attack from 2.5 to 47 deg from a flight experiment with a cone that reentered at angles of attack up to 75 deg were analyzed and their local flow conditions are presented. The transition data were obtained from both acoustic and electrostatic sensors. The data from the acoustic and electrostatic sensors were correlated by use of three different sets of correlating parameters. For each set of correlating parameters, the transitional and turbulent data from the acoustic and electrostatic sensors were plotted separately. Each of the correlations was compared with a linear curve fit of previous cone-flight data at near zero angle of attack. Results indicate that the data from acoustic sensors falls below the linear curve fits to the flight data and indicates the sensitivity of the sensors. It is shown that the transition data from the electrostatic sensors tends to scatter much more than the data from the acoustic sensors. However, the data from the electrostatic sensors tends to scatter about the linear curve fits based on previous flight data.
A system and method of measuring a residence time in a gas-turbine engine is provided, whereby the method includes placing pressure sensors at a combustor entrance and at a turbine exit of the gas-turbine engine and measuring a combustor pressure at the combustor entrance and a turbine exit pressure at the turbine exit. The method further includes computing cross-spectrum functions between a combustor pressure sensor signal from the measured combustor pressure and a turbine exit pressure sensor signal from the measured turbine exit pressure, applying a linear curve fit to the cross-spectrum functions, and computing a post-combustion residence time from the linear curve fit.
To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of: (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disproportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit = -10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1 %) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucteotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life.
To identify the energy source that drives the biosynthesis of amino acids, lipids, and nucleotides from glucose, we calculated the free energy change due to redox disproportionation of the substrate carbon of: (1) 26-carbon fermentation reactions and (2) the biosynthesis of amino acids and lipids of E. coli from glucose. The free energy (cal/mmol of carbon) of these reactions was plotted as a function of the degree of redox disproportionation of carbon (disproportionative electron transfers (mmol)/mmol of carbon). The zero intercept and proportionality between energy yield and degree of redox disproportionation exhibited by this plot demonstrate that redox disproportionation is the principal energy source of these redox reactions (slope of linear fit = -10.4 cal/mmol of disproportionative electron transfers). The energy and disproportionation values of E. coli amino acid and lipid biosynthesis from glucose lie near this linear curve fit with redox disproportionation accounting for 84% and 96% (and ATP only 6% and 1%) of the total energy of amino acid and lipid biosynthesis, respectively. These observations establish that redox disproportionation of carbon, and not ATP, is the primary energy source driving amino acid and lipid biosynthesis from glucose. In contrast, we found that nucleotide biosynthesis involves very little redox disproportionation, and consequently depends almost entirely on ATP for energy. The function of sugar redox disproportionation as the major source of free energy for the biosynthesis of amino acids and lipids suggests that sugar disproportionation played a central role in the origin of metabolism, and probably the origin of life.
This paper presents an improved solution for curve fitting data to an exponential equation (Y = AeBt + C). This improvement is in four areas ? speed, stability, determinant processing time, and the removal of limits. The solution presented in this paper avoids iterative techniques and their stability errors by using three mathematical ideas ? discrete calculus, a special relationship (between exponential curves and the Mean Value Theorem for Derivatives), and a simple linear curve fit algorithm. This method can also be applied to fitting data to the general power law equation Y = AxB + C and the general geometric growth equation Y = AkBt + C.
An improved solution for curve fitting data to an exponential equation (y = Ae(exp Bt) + C) has been developed. This improvement is in four areas -- speed, stability, determinant processing time, and the removal of limits. The solution presented avoids iterative techniques and their stability errors by using three mathematical ideas: discrete calculus, a special relationship (be tween exponential curves and the Mean Value Theorem for Derivatives), and a simple linear curve fit algorithm. This method can also be applied to fitting data to the general power law equation y = Ax(exp B) + C and the general geometric growth equation y = Ak(exp Bt) + C.
A computer-oriented technique is presented for performing a nonlinear exponential regression analysis on decay-type experimental data. The technique involves the least squares procedure wherein the nonlinear problem is linearized by expansion in a Taylor series. A linear curve fitting procedure for determining the initial nominal estimates for the unknown exponential model parameters is included as an integral part of the technique. A correction matrix was derived and then applied to the nominal estimate to produce an improved set of model parameters. The solution cycle is repeated until some predetermined criterion is satisfied.
Signal processing techniques used in Mesosphere-Stratosphere-Troposphere (MST) radars are reviewed. Techniques which produce good estimates of the total power, frequency shift, and spectral width of the radar power spectra are considered. Non-linear curve fitting, autocovariance, autocorrelation, covariance, and maximum likelihood estimators are discussed.
The structure and kinetic energy budget of turbulence in the convective marine atmospheric boundary layer as observed by aircraft during a cold air outbreak have been studied using mixed layer scaling. The results are significantly different from those of previous studies under conditions closer to free convection. The normalized turbulent kinetic energy and turbulent transport are about twice those found during the Air Mass Transformation Experiment (AMTEX). This implies that for a given surface heating the present case is dynamically more active. The difference is mainly due to the greater importance of wind shear in the present case. This case is closer to the roll vortex regime, whereas AMTEX observed mesoscale cellular convection which is closer to free convection. Shear generation is found to provide a significant energy source, in addition to buoyancy production, to maintain a larger normalized turbulent kinetic energy and to balance a larger normalized dissipation. The interaction between turbulent pressure and divergence (i.e., pressure scrambling) is also found to transfer energy from the vertical to the horizontal components, and is expected to be stronger in roll vortices than in m esoscale cells. The sensible heat flux is found to fit well with a linear vertical profile in a clear or subcloud planetary boundary layer (PBL), in good agreement with the results of Lenschow et al., (1980). The heat flux ratio between the PBL top and the surface, derived from the linear fitted curve, is approximately -0.14, in good agreement with that derived from the lidar data for the same case. Near the PBL top, the heat flux profiles are consistent with those of Deardoff (1979) and Deardorff et al. (1980).
Previously reported methods of forecasting lightning threat using fields of graupel flux from WRF simulations are extended to include the simulated field of vertically integrated ice within storms. Although the ice integral shows less temporal variability than graupel flux, it provides more areal coverage, and can thus be used to create a lightning forecast that better matches the areal coverage of the lightning threat found in observations of flash extent density. A blended lightning forecast threat can be constructed that retains much of the desirable temporal sensitivity of the graupel flux method, while also incorporating the coverage benefits of the ice integral method. The graupel flux and ice integral fields contributing to the blended forecast are calibrated against observed lightning flash origin density data, based on Lightning Mapping Array observations from a series of case studies chosen to cover a wide range of flash rate conditions. Linear curve fits that pass through the origin are found to be statistically robust for the calibration procedures.
Curve fitting for modified paired comparisons using linearized nonlinear regression
Approximation of noisy data in the plane by straight lines or elliptic or single-branch hyperbolic curve segments arises in pattern recognition, data compaction, and other problems. The efficient search for and approximation of data by such curves were examined. Recursive least-squares linear curve-fitting was used, and ellipses and hyperbolas are parameterized as quadratic functions in x and y. The error minimized by the algorithm is interpreted, and central processing unit (CPU) times for estimating parameters for fitting straight lines and quadratic curves were determined and compared. CPU time for data search was also determined for the case of straight line fitting. Quadratic curve fitting is shown to require about six times as much CPU time as does straight line fitting, and curves relating CPU time and fitting error were determined for straight line fitting. Results are derived on early sequential determination of whether or not the underlying curve is a straight line.
An empirical method is developed to predict certain parameters of future solar activity cycles. Sunspot cycle statistics are examined, and curve fitting and linear regression analysis techniques are utilized.
A thrust stand calibration technique for use in testing repetitively-pulsed electric thrusters for in-space propulsion has been developed and tested using a modified hanging pendulum thrust stand. In the implementation of this technique, current pulses are applied to a solenoidal coil to produce a pulsed magnetic field that acts against the magnetic field produced by a permanent magnet mounted to the thrust stand pendulum arm. The force on the magnet is applied in this non-contact manner, with the entire pulsed force transferred to the pendulum arm through a piezoelectric force transducer to provide a time-accurate force measurement. Modeling of the pendulum arm dynamics reveals that after an initial transient in thrust stand motion the quasisteady average deflection of the thrust stand arm away from the unforced or zero position can be related to the average applied force through a simple linear Hooke s law relationship. Modeling demonstrates that this technique is universally applicable except when the pulsing period is increased to the point where it approaches the period of natural thrust stand motion. Calibration data were obtained using a modified hanging pendulum thrust stand previously used for steady-state thrust measurements. Data were obtained for varying impulse bit at constant pulse frequency and for varying pulse frequency. The two data sets exhibit excellent quantitative agreement with each other as the constant relating average deflection and average thrust match within the errors on the linear regression curve fit of the data. Quantitatively, the error on the calibration coefficient is roughly 1% of the coefficient value.
Experimental results are presented for the elastic and strength properties of T300/5208 graphite-epoxy at room temperature, 116K (-250 F), and 394K (+250 F). Results are presented for unidirectional 0, 90, and 45 degree laminates, and + or - 30, + or - 45, and + or - 60 degree angle-ply laminates. The stress-strain behavior of the 0 and 90 degree laminates is essentially linear for all three temperatures and that the stress-strain behavior of all other laminates is linear at 116K. A second-order curve provides the best fit for the temperature is linear at 116K. A second-order curve provides the best fit for the temperature dependence of the elastic modulus of all laminates and for the principal shear modulus. Poisson's ratio appears to vary linearly with temperature. all moduli decrease with increasing temperature except for E (sub 1) which exhibits a small increase. The strength temperature dependence is also quadratic for all laminates except the 0 degree - laminate which exhibits linear temperature dependence. In many cases the temperature dependence of properties is nearly linear.
Simple relations for determining the enthalpy and temperature of hydrogen-helium gas mixtures were developed for hydrogen volumetric compositions from 1.0 to 0.7. These relations are expressed as a function of pressure and density and are valid for a range of temperatures from 7,000 to 35,000 K and pressures from 0.10 to 3.14 MPa. The proportionality constant and exponents in the correlation equations were determined for each gas composition by applying a linear least squares curve fit to a large number of thermodynamic calculations obtained from a detailed computer code. Although these simple relations yielded thermodynamic properties suitable for many engineering applications, their accuracy was improved significantly by evaluating the proportionality constants at postshock conditions and correlating these values as a function of the gas composition and the product of freestream velocity and shock angle. The resulting equations for the proportionality constants in terms of velocity and gas composition and the corresponding simple realtions for enthalpy and temperature were incorporated into a flow field computational scheme. Comparison was good between the thermodynamic properties determined from these relations and those obtained by using a detailed computer code to determine the properties. Thus, an appreciable savings in computer time was realized with no significant loss in accuracy.
Landsat radiance data were used to test mathematical models relating diffuse reflectance to aquatic suspended solids concentration. Digital CCT data for Landsat passes over the Bay of Fundy, Nova Scotia were analyzed on a General Electric Co. Image 100 multispectral analysis system. Three data sets were studied separately and together in all combinations with and without solar angle correction. Statistical analysis and chromaticity analysis show that a nonlinear relationship between Landsat radiance and suspended solids concentration is better at curve-fitting than a linear relationship. In particular, the quasi-single-scattering diffuse reflectance model developed by Gordon and coworkers is corroborated. The Gordon model applied to 33 points of MSS 5 data combined from three dates produced r = 0.98.