Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stochastic Processes”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1999

Scaling in nature: from DNA through heartbeats to weather

The purpose of this report is to describe some recent progress in applying scaling concepts to various systems in nature. We review several systems characterized by scaling laws such as DNA sequences, heartbeat rates and weather variations. We discuss the finding that the exponent alpha quantifying the scaling in DNA in smaller for coding than for noncoding sequences. We also discuss the application of fractal scaling analysis to the dynamics of heartbeat regulation, and report the recent finding that the scaling exponent alpha is smaller during sleep periods compared to wake periods. We also discuss the recent findings that suggest a universal scaling exponent characterizing the weather fluctuations.

Review↗
NASA NTRS1999

Statistical physics and physiology: monofractal and multifractal approaches

Even under healthy, basal conditions, physiologic systems show erratic fluctuations resembling those found in dynamical systems driven away from a single equilibrium state. Do such "nonequilibrium" fluctuations simply reflect the fact that physiologic systems are being constantly perturbed by external and intrinsic noise? Or, do these fluctuations actually, contain useful, "hidden" information about the underlying nonequilibrium control mechanisms? We report some recent attempts to understand the dynamics of complex physiologic fluctuations by adapting and extending concepts and methods developed very recently in statistical physics. Specifically, we focus on interbeat interval variability as an important quantity to help elucidate possibly non-homeostatic physiologic variability because (i) the heart rate is under direct neuroautonomic control, (ii) interbeat interval variability is readily measured by noninvasive means, and (iii) analysis of these heart rate dynamics may provide important practical diagnostic and prognostic information not obtainable with current approaches. The analytic tools we discuss may be used on a wider range of physiologic signals. We first review recent progress using two analysis methods--detrended fluctuation analysis and wavelets--sufficient for quantifying monofractual structures. We then describe recent work that quantifies multifractal features of interbeat interval series, and the discovery that the multifractal structure of healthy subjects is different than that of diseased subjects.

NASA Discipline Cardiopulmonary↗
NASA NTRS1998

Stochastic feedback and the regulation of biological rhythms

We propose a general approach to the question of how biological rhythms spontaneously self-regulate, based on the concept of "stochastic feedback". We illustrate this approach by considering at a coarse-grained level the neuroautonomic regulation of the heart rate. The model generates complex dynamics and successfully acounts for key characteristics of cardiac variability, including the l/f power spectrum, the functional form and scaling of the distribution of variations, and correlations in the Fourier phases indicating nonlinear dynamics.

NASA Discipline Cardiopulmonary↗
NASA NTRS1999

Modeling heart rate variability by stochastic feedback

We consider the question of how the cardiac rhythm spontaneously self-regulates and propose a new mechanism as a possible answer. We model the neuroautonomic regulation of the heart rate as a stochastic feedback system and find that the model successfully accounts for key characteristics of cardiac variability, including the 1/f power spectrum, the functional form and scaling of the distribution of variations of the interbeat intervals, and the correlations in the Fourier phases which indicate nonlinear dynamics.

Non-NASA Center↗
NASA NTRS1997

Implicit Kalman filtering

For an implicitly defined discrete system, a new algorithm for Kalman filtering is developed and an efficient numerical implementation scheme is proposed. Unlike the traditional explicit approach, the implicit filter can be readily applied to ill-conditioned systems and allows for generalization to descriptor systems. The implementation of the implicit filter depends on the solution of the congruence matrix equation (A1)(Px)(AT1) = Py. We develop a general iterative method for the solution of this equation, and prove necessary and sufficient conditions for convergence. It is shown that when the system matrices of an implicit system are sparse, the implicit Kalman filter requires significantly less computer time and storage to implement as compared to the traditional explicit Kalman filter. Simulation results are presented to illustrate and substantiate the theoretical developments.

Non-NASA Center↗
NASA NTRS1997

Deviations from uniform power law scaling in nonstationary time series

A classic problem in physics is the analysis of highly nonstationary time series that typically exhibit long-range correlations. Here we test the hypothesis that the scaling properties of the dynamics of healthy physiological systems are more stable than those of pathological systems by studying beat-to-beat fluctuations in the human heart rate. We develop techniques based on the Fano factor and Allan factor functions, as well as on detrended fluctuation analysis, for quantifying deviations from uniform power-law scaling in nonstationary time series. By analyzing extremely long data sets of up to N = 10(5) beats for 11 healthy subjects, we find that the fluctuations in the heart rate scale approximately uniformly over several temporal orders of magnitude. By contrast, we find that in data sets of comparable length for 14 subjects with heart disease, the fluctuations grow erratically, indicating a loss of scaling stability.

Non-NASA Center↗
NASA NTRS1997

Detection of chaotic determinism in time series from randomly forced maps

Time series from biological system often display fluctuations in the measured variables. Much effort has been directed at determining whether this variability reflects deterministic chaos, or whether it is merely "noise". Despite this effort, it has been difficult to establish the presence of chaos in time series from biological sytems. The output from a biological system is probably the result of both its internal dynamics, and the input to the system from the surroundings. This implies that the system should be viewed as a mixed system with both stochastic and deterministic components. We present a method that appears to be useful in deciding whether determinism is present in a time series, and if this determinism has chaotic attributes, i.e., a positive characteristic exponent that leads to sensitivity to initial conditions. The method relies on fitting a nonlinear autoregressive model to the time series followed by an estimation of the characteristic exponents of the model over the observed probability distribution of states for the system. The method is tested by computer simulations, and applied to heart rate variability data.

Non-NASA Center↗
NASA NTRS1996

Possible effects of protracted exposure on the additivity of risks from space radiations

Conventional radiation risk assessments are presently based on the additivity assumption. This assumption states that risks from individual components of a complex radiation field involving many different types of radiation can be added to yield the total risk of the complex radiation field. If the assumption is not correct, the summations and integrations performed to obtain the presently quoted risk estimates are not appropriate. This problem is particularly important in the area of space radiation risk evaluation because of the many different types of high- and low-LET radiation present in the galactic cosmic ray environment. For both low- and high-LET radiations at low enough dose rates, the present convention is that the addivity assumption holds. Mathematically, the total risk, Rtot is assumed to be Rtot = summation (i) Ri where the summation runs over the different types of radiation present. If the total dose (or fluence) from each component is such that the interaction between biological lesions caused by separate single track traversals is negligible within a given cell, it is presently considered to be reasonable to accept the additivity assumption. However, when the exposure is protracted over many cell doubling times (as will be the case for extended missions to the moon or Mars), the possibility exists that radiation effects that depend on multiple cellular events over a long time period, such as is probably the case in radiation-induced carcinogenesis, may not be additive in the above sense and the exposure interval may have to be included in the evaluation procedure. It is shown, however, that "inverse" dose-rate effects are not expected from intermediate LET radiations arising from the galactic cosmic ray environment due to the "sensitive-window-in-the-cell-cycle" hypothesis.

NASA Discipline Number 45-10↗
NASA NTRS1995

Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series

The healthy heartbeat is traditionally thought to be regulated according to the classical principle of homeostasis whereby physiologic systems operate to reduce variability and achieve an equilibrium-like state [Physiol. Rev. 9, 399-431 (1929)]. However, recent studies [Phys. Rev. Lett. 70, 1343-1346 (1993); Fractals in Biology and Medicine (Birkhauser-Verlag, Basel, 1994), pp. 55-65] reveal that under normal conditions, beat-to-beat fluctuations in heart rate display the kind of long-range correlations typically exhibited by dynamical systems far from equilibrium [Phys. Rev. Lett. 59, 381-384 (1987)]. In contrast, heart rate time series from patients with severe congestive heart failure show a breakdown of this long-range correlation behavior. We describe a new method--detrended fluctuation analysis (DFA)--for quantifying this correlation property in non-stationary physiological time series. Application of this technique shows evidence for a crossover phenomenon associated with a change in short and long-range scaling exponents. This method may be of use in distinguishing healthy from pathologic data sets based on differences in these scaling properties.

NASA Discipline Number 14-10↗
NASA NTRS2001

Noise effects on the complex patterns of abnormal heartbeats

Patients at high risk for sudden death often exhibit complex heart rhythms in which abnormal heartbeats are interspersed with normal heartbeats. We analyze such a complex rhythm in a single patient over a 12-h period and show that the rhythm can be described by a theoretical model consisting of two interacting oscillators with stochastic elements. By varying the magnitude of the noise, we show that for an intermediate level of noise, the model gives best agreement with key statistical features of the dynamics.

NASA Discipline Cardiopulmonary↗
NASA NTRS2000

Biological effects of cosmic radiation: deterministic and stochastic

Our basic understanding of the biological responses to cosmic radiations comes in large part from an international series of ground-based laboratory studies, where accelerators have provided the source of representative charged particle radiations. Most of the experimental studies have been performed using acute exposures to a single radiation type at relatively high doses and dose rates. However, most exposures in flight occur from low doses of mixed radiation fields at low-dose rates. This paper provides a brief overview of existing pertinent clinical and biological radiation data and the limitations associated with data available from specific components of the radiation fields in airflight and space travel.

Flight Experiment↗
NASA NTRS1993

Numerical experiments with model monophyletic and paraphyletic taxa

The problem of how accurately paraphyletic taxa versus monophyletic (i.e., holophyletic) groups (clades) capture underlying species patterns of diversity and extinction is explored with Monte Carlo simulations. Phylogenies are modeled as stochastic trees. Paraphyletic taxa are defined in an arbitrary manner by randomly choosing progenitors and clustering all descendants not belonging to other taxa. These taxa are then examined to determine which are clades, and the remaining paraphyletic groups are dissected to discover monophyletic subgroups. Comparisons of diversity patterns and extinction rates between modeled taxa and lineages indicate that paraphyletic groups can adequately capture lineage information under a variety of conditions of diversification and mass extinction. This suggests that these groups constitute more than mere "taxonomic noise" in this context. But, strictly monophyletic groups perform somewhat better, especially with regard to mass extinctions. However, when low levels of paleontologic sampling are simulated, the veracity of clades deteriorates, especially with respect to diversity, and modeled paraphyletic taxa often capture more information about underlying lineages. Thus, for studies of diversity and taxic evolution in the fossil record, traditional paleontologic genera and families need not be rejected in favor of cladistically-defined taxa.

Non-NASA Center↗
NASA NTRS1994

Limits to randomness in paleobiologic models: the case of Phanerozoic species diversity

The question of how random, or unconstrained, paleobiologic models should be is examined with a case study: Signor's (1982, 1985) inverse calculation of levels of marine species diversity through the Phanerozoic. His calculation involved an ingenious model that estimated species numbers and species abundances in the world oceans of the past by correcting known numbers of fossil species for variations in sedimentary rocks available for sampling and in effort paleontologists might devote to sampling. The model proves robust to changes in possible shapes of species-abundance distributions, but it is sensitive to alterations in the assumption that paleontologists collect fossils at random. If it is assumed that ease of collecting varies with age of sediment (with the Cenozoic offering easy sampling) or that paleontologists tend to seek out rarer fossils, results of the inverse calculation change. In particular, the magnitude of the calculated Cenozoic diversity increase always declines from the factor of about seven as originally reported to something considerably smaller. This leaves open the problem of the magnitude of Cenozoic increase in marine species diversity, awaiting better empirical data and, perhaps, more exacting models, random or otherwise.

Non-NASA Center↗
NASA NTRS1994

Physiological time-series analysis: what does regularity quantify?

Approximate entropy (ApEn) is a recently developed statistic quantifying regularity and complexity that appears to have potential application to a wide variety of physiological and clinical time-series data. The focus here is to provide a better understanding of ApEn to facilitate its proper utilization, application, and interpretation. After giving the formal mathematical description of ApEn, we provide a multistep description of the algorithm as applied to two contrasting clinical heart rate data sets. We discuss algorithm implementation and interpretation and introduce a general mathematical hypothesis of the dynamics of a wide class of diseases, indicating the utility of ApEn to test this hypothesis. We indicate the relationship of ApEn to variability measures, the Fourier spectrum, and algorithms motivated by study of chaotic dynamics. We discuss further mathematical properties of ApEn, including the choice of input parameters, statistical issues, and modeling considerations, and we conclude with a section on caveats to ensure correct ApEn utilization.

NASA Discipline Cardiopulmonary↗
NASA NTRS2004

cWINNOWER algorithm for finding fuzzy dna motifs

The cWINNOWER algorithm detects fuzzy motifs in DNA sequences rich in protein-binding signals. A signal is defined as any short nucleotide pattern having up to d mutations differing from a motif of length l. The algorithm finds such motifs if a clique consisting of a sufficiently large number of mutated copies of the motif (i.e., the signals) is present in the DNA sequence. The cWINNOWER algorithm substantially improves the sensitivity of the winnower method of Pevzner and Sze by imposing a consensus constraint, enabling it to detect much weaker signals. We studied the minimum detectable clique size qc as a function of sequence length N for random sequences. We found that qc increases linearly with N for a fast version of the algorithm based on counting three-member sub-cliques. Imposing consensus constraints reduces qc by a factor of three in this case, which makes the algorithm dramatically more sensitive. Our most sensitive algorithm, which counts four-member sub-cliques, needs a minimum of only 13 signals to detect motifs in a sequence of length N = 12,000 for (l, d) = (15, 4). Copyright Imperial College Press.

Evaluation Studies↗
NASA NTRS2000

Estimation of EOP From VLBI: Direct Approach

The currently adopted strategy of Earth Orientation Parameters (EOP) estimation from Very Long Baseline Interferometry (VLBI) is to estimate six parameters: Universal Time 1 (UT1), UT1 rate, pole positions, and nutation offsets for each 24-hour session independently. Then the resulting time series of raw EOP are filtered and a regression analysis is performed to obtain nutation coefficients, polhode of the pole, and other physical parameters. Thus, the latter parameters are obtained indirectly in two stages. An alternative approach of direct estimation of the final EOP is presented. Pole coordinates and UT1 are considered as a sum of three components: the low-period component that is modeled by a cubic spline, the harmonic component that includes forced nutation, precession and sub-daily variations of EOP, and the stochastic component that is modeled by a linear spline with segment length 1-2 hours. All parameters are obtained in a single LSQ solution using all available data.

POLAR MOTION↗
NASA NTRS2005

Stochastic Representation of Chaos using Terminal Attractors

A nonlinear version of the Liouville equation based upon terminal attractors is proposed for describing post-instability motions of dynamical systems with exponential divergence of trajectories such as those leading to chaos and turbulence. As a result, the post-instability motions are represented by expectations, variances, and higher moments of the state variables as functions of time. The proposed approach can be applied to conservative chaos, and in particular, to n-bodies problem, as well as to dissipative systems, and in particular, to chaotic attractors and turbulence.

stochastic processes↗