A computational model of the error detector of human visual accommodation

A mathematical model is proposed for the error detector of the human visual accommodative system. The model supposes that the accommodative error detector derives both the direction and the magnitude of the accommodative error from naturally-occuring oscillations of the lens and their effects on retinal-image contrast. Differential operators take the first derivatives of two time varying functions: lens power and retinal-image contrast. Directional information is obtained by comparing the signs of these two derivatives and magnitude information is obtained by comparing their amplitudes.Research conducted at the School of Optometry, University of California, BerkeleySupported by National Eye Institute grant EYO-3532-04(C.S.) and National Institutes of Health core grant # 1-445420-32011     

A computational model for dynamic analysis of the human gait

Biomechanical models are important tools in the study of human motion. This work proposes a computational model to analyse the dynamics of lower limb motion using a kinematic chain to represent the body segments and rotational joints linked by viscoelastic elements. The model uses anthropometric parameters, ground reaction forces and joint Cardan angles from subjects to analyse lower limb motion during the gait. The model allows evaluating these data in each body plane. Six healthy subjects walked on a treadmill to record the kinematic and kinetic data. In addition, anthropometric parameters were recorded to construct the model. The viscoelastic parameter values were fitted for the model joints (hip, knee and ankle). The proposed model demonstrated that manipulating the viscoelastic parameters between the body segments could fit the amplitudes and frequencies of motion. The data collected in this work have viscoelastic parameter values that follow a normal distribution, indicating that these values are directly related to the gait pattern. To validate the model, we used the values of the joint angles to perform a comparison between the model results and previously published data. The model results show a same pattern and range of values found in the literature for the human gait motion.     

A mathematical model of the human thermal system

This paper describes a mathematical model developed to simulate the physical characteristics of the human thermal system in the transient state. Physiological parameters, such as local metabolic heat generation rates, local blood flow rates, and rates of sweating, must be specified as input data. Automatic computation of these parameters will be built into the model at a later date when it is used to study thermal regulation in the human. Finite-difference techniques have been used to solve the heat conduction equation on a Control Data Corporation 1604 computer. Since numerical techniques were used, it was possible to include many more factors in this model than in previous ones. The body was divided into 15 geometric regions, which were the head, the thorax, the abdomen, and the proximal, medial, and distal segments of the arms and legs. Axial gradients in a given segment were neglected. In each segment, the large arteries and veins were approximated by an arterial pool and a venous pool which were distributed radially throughout the segment. Accumulation of heat in the blood of the large arteries and veins, and heat transfer from the large arteries and veins to the surrounding tissue were taken into account. The venous streams were collected together at the heart before flowing into the capillaries of the lungs. Each of the segments was subdivided into 15 radial sections, thereby allowing considerable freedom in the assignment of physical properties such as thermal conductivity and rate of blood flow to the capillaries. The program has been carefully checked for errors, and it is now being used to analyze some problems of current interest. This study was supported by the office of the Surgeon General, U.S. Army, under contract no. DA 49-193-MD-2005.     

A normal ocular motor system model that simulates the dual-mode fast phases of latent/manifest latent nystagmus

 The fast phases of latent/manifest latent nystagmus (LMLN) may either cause the target image to fall within (foveating) or outside (defoveating) the foveal area. We previously verified that both types are generated by the same mechanism as voluntary saccades and propose a hypothetical, dual-mode mechanism (computer model) for LMLN that utilizes normal ocular-motor control functions. Fixation data recorded during the past 30 years from 97 subjects with LMLN using both infrared and magnetic search coil oculography were used as a basis for our simulations. The MATLAB/Simulink software was used to construct a robust, modular, ocular motor system model, capable of simulating LMLN. Fast-phase amplitude versus both peak velocity and duration of simulated saccades were equivalent to those of saccades in normal subjects. Based on our LMLN studies, we constructed a hypothetical model in which the slow-phase velocity acted to trigger the change between foveating and defoveating LMLN fast phases. Foveating fast phases were generated during lower slow-phase velocities whereas defoveating fast phases occurred during higher slow-phase velocities. The bidirectional model simulated Alexander's law behavior under all viewing and fixation conditions. Our ocular-motor model accurately simulates LMLN patient ocular motility data and provides a hypothetical explanation for the conditions that result in both foveating and defoveating fast phases. As is the case for normal physiological saccades, the position error determined the saccadic amplitudes for foveating fast phases. However, the final slow-phase velocity determined the amplitudes of defoveating fast phases. In addition, we suggest that individuals with LMLN use their fixation subsystem to further decrease the slow-phase velocity as the target image approaches the foveal center. Received: 16 June 2000 / Accepted in revised form: 20 May 2001     

A two-scale network dynamic model for human mobility process

The technological changes and educational expansion have created the heterogeneity in the human species. Clearly, this heterogeneity generates a structure in the population dynamics, namely: citizen, permanent resident, visitor, and etc. Furthermore, as the heterogeneity in the population increases, the human mobility between meta-populations patches also increases. Depending on spatial scales, a meta-population patch can be decomposed into sub-patches, for examples: homes, neighborhoods, towns, etc. Members of the population can move between the sub-patches. The dynamics of human mobility in a heterogeneous and scaled structured population is still its infancy level. In this work, an attempt is made to investigate the human mobility dynamics of heterogeneous and scaled structured population. We present a two scaled human mobility model for a meta-population. The sub regions and regions are interlinked via intra-and inter regional transport network systems. Under various types of growth order assumptions on the intra and interregional residence times of the residents of a sub region, different patterns of static behavior of the mobility process is studied. In addition, the results reveal that the system has a natural tendency to quarantine itself without total breaking a link in the transportation network system. Moreover, there is a threshold point for the largest intra regional visiting time of residents of a given sub region that leads to either a total isolation of the residents from other sub regions within the region or a partial isolation of residents from some of the sub regions within the region.     

A model of the neuro-musculo-skeletal system for human locomotion

The generation of human locomotion was examined by linking computational neuroscience with biomechanics from the perspective of nonlinear dynamical theory. We constructed a model of human locomotion, which includes a musculo-skeletal system with 8 segments and 20 muscles, a neural rhythm generator composed of 7 pairs of neural oscillators, and mechanisms for processing and transporting sensory and motor signals. Using a computer simulation, we found that locomotion emerged as a stable limit cycle that was generated by the global entrainment between the musculo-skeletal system, the neural system, and the environment. Moreover, the walking movements of the model could be compared quantitatively with those of experimental studies in humans.Part of this paper was presented to IVth International Symposium on Computer Simulation in Biomechanics, Paris, France, July 1, 1993     

A dynamic root system growth model based on L-Systems

Understanding the impact of roots and rhizosphere traits on plant resource efficiency is important, in particular in the light of upcoming shortages of mineral fertilizers and climate change with increasing frequency of droughts. We developed a modular approach to root growth and architecture modelling with a special focus on soil root interactions. The dynamic three-dimensional model is based on L-Systems, rewriting systems well-known in plant architecture modelling. We implemented the model in Matlab in a way that simplifies introducing new features as required. Different kinds of tropisms were implemented as stochastic processes that determine the position of the different roots in space. A simulation study was presented for phosphate uptake by a maize root system in a pot experiment. Different sink terms were derived from the root architecture, and the effects of gravitropism and chemotropism were demonstrated. This root system model is an open and flexible tool which can easily be coupled to different kinds of soil models.     

干扰条件下捕食者与被捕食者系统动态模型

建立了在干扰条件下捕食者与被捕食者系统动态模型,得出系统的演化方程和演化曲线及演化相图,分别讨论了无干扰和有干扰影响情况下物种的演化情况,说明考虑了干扰后的模型更为合理,也更有利于物种的进化和新物种的产生.     

A dynamic model for bilirubin binding to human serum albumin

Site-directed mutagenesis of human serum albumin was used to study the role of various amino acid residues in bilirubin binding. A comparison of thermodynamic, proteolytic, and x-ray crystallographic data from previous studies allowed a small number of amino acid residues in subdomain 2A to be selected as targets for substitution. The following recombinant human serum albumin species were synthesized in the yeast species Pichia pastoris: K195M, K199M, F211V, W214L, R218M, R222M, H242V, R257M, and wild type human serum albumin. The affinity of bilirubin was measured by two independent methods and found to be similar for all human serum albumin species. Examination of the absorption and circular dichroism spectra of bilirubin bound to its high affinity site revealed dramatic differences between the conformations of bilirubin bound to the above human serum albumin species. The absorption and circular dichroism spectra of bilirubin bound to the above human serum albumin species in aqueous solutions saturated with chloroform were also examined. The effect of certain amino acid substitutions on the conformation of bound bilirubin was altered by the addition of chloroform. In total, the present study suggests a dynamic, unusually flexible high affinity binding site for bilirubin on human serum albumin.     

A human cell model for dynamic testing of MR contrast agents

To determine the initial feasibility of using magnetic resonance (MR) imaging to detect early atherosclerosis, we investigated inflammatory cells labeled with a positive contrast agent in an endothelial cell-based testing system. The human monocytic cell line THP-1 was labeled by overnight incubation with a gadolinium colloid (Gado CELLTrack) prior to determination of the in vitro release profile from T1-weighted MR images. Next, MR signals arising from both a synthetic model of THP-1/human umbilical vein endothelial cell (HUVEC) accumulation and the dynamic adhesion of THP-1 cells to activated HUVECs under flow were obtained. THP-1 cells were found to be successfully--but not optimally--labeled with gadolinium colloid, and MR images demonstrated increased signal from labeled cells in both the synthetic and dynamic THP-1/HUVEC models. The observed THP-1 contrast release profile was rapid, suggesting the need for an agent that is optimized for retention in the target cells for use in further studies. Detection of labeled THP-1 cells was accomplished with no signal enhancement from unlabeled cells. These achievements demonstrate the feasibility of targeting early atherosclerosis with MR imaging, and suggest that using an in vitro system like the one described provides a rapid, efficient, and cost-effective way to support the development and evaluation of novel MR contrast agents.     

A joint model of the contractile system of striated muscle (dynamic version)

The Joint Model of the Contractile System of Muscle is a theoretical construction and the question is whether it may express adequately some aspects of the biological original. In addition to the previous results which presented preliminary calculations, simulation experiments have been performed with a dynamic version of the model. The relations found between the variables studied (length, shortening velocity, isotonic contraction load, shortening heat production rate and excess of effective activity), did not essentially differ from those established by preliminary calculations; the choice of appropriate parameter values and the introduction of some additional features in the system made it possible to obtain a better approximation of relations known from measurements on the biological object. The model discloses some relevant general aspects of the contractile system of "independent" generators of force and thus contributes to a specification of the relationships between molecular and macroscopic levels of contraction phenomena. Future development or correction of the concepts under consideration depend on more suitable empirical data and further simulation experiments.     

Analysis of human crystalline lens accommodation

The behavior of the human crystalline lens during accommodation is analytically studied. The lens is modeled as a closed axisymmetrical membrane shell of varying thickness enclosing an incompressible liquid. To simulate zonular tension associated with lenticular accommodation, an axisymmetrical radial force or displacement is imposed around the shell equator. Two second-order, simultaneous, nonlinear governing differential equations are derived. Numerical results, obtained from the investigation of human lens profiles of three independently published MRI images and a drawing of a microphotograph, demonstrate that when zonular traction within the physiological force range of the ciliary muscle is exerted, both central lens thickness and central optical power increase. Qualitatively, these increases are independent of lens shape. However, the magnitude of these changes is dependent on the initial profile of the lens and is enhanced by the "natural" variation in capsular thickness. Only when a pulling force significantly exceeds the force capacity of the ciliary muscle does the lens flatten and its central thickness and optical power decrease.     

A ratio model of perceived speed in the human visual system

The perceived speed of moving images changes over time. Prolonged viewing of a pattern (adaptation) leads to an exponential decrease in its perceived speed. Similarly, responses of neurones tuned to motion reduce exponentially over time. It is tempting to link these phenomena. However, under certain conditions, perceived speed increases after adaptation and the time course of these perceptual effects varies widely. We propose a model that comprises two temporally tuned mechanisms whose sensitivities reduce exponentially over time. Perceived speed is taken as the ratio of these filters' outputs. The model captures increases and decreases in perceived speed following adaptation and describes our data well with just four free parameters. Whilst the model captures perceptual time courses that vary widely, parameter estimates for the time constants of the underlying filters are in good agreement with estimates of the time course of adaptation of direction selective neurones in the mammalian visual system.     

A dynamic mathematical model of the chemostat

A number of experimental studies on the dynamic, behavior of the chemostat have shown that the specific growth rate does not, instantaneously adjust to changes in the concentration of limiting substrate in the chemostat following disturbances in the steady state input limiting substrate concentration or in the steady state dilution rate. Instead of an instantaneous response, as would be predicted by the Monod equation, experimental studies have shown that the specific growth rate experiences a dynamic lag in responding to the changes in the concentration of limiting substrate in the culture vessel. The observed dynamic lag has been recognized by researchers in such terms as an inertial phenomenon and as a hysteresis effect, but as yet a systems engineering approach has not been applied to the observed data. The present paper criticizes the use of the Monod equation as a dynamic relationship and offers as an alternative a dynamic equation relating specific growth rate to the limiting substrate concentration in the chemostat. Following the development of equations, experimental methods of evaluating parameters are discussed. Dynamic responses of analog simulations (incorporating the newly derived equations) are compared with the dynamic responses predicted by the Monod equation and with the dynamic responses of experimental chemostats.     

A dynamic model of the vertebrate retina

In this paper a mathematical model of the retina was proposed to clarify the spatio-temporal information processing mechanism in the retina of vertebrates. In order to explain spatio-temporal characteristics of an on-center receptive field of a ganglion cell, excitatory and inhibitory cell layers were introduced of which time lags increased with the lateral distance from a point of stimulation. The characteristics of this model were found to agree well with the physiological data: e.g., this model shows on-response to the input stimulus given on the center, off-response to the input on the surround, and on-off response to the input on the border between on- and off-response regions of the on-center field.