ALBERT

Notizen: Abstract The partial differential equation of the random walk problem with persistence of direction and external bias is derived. By persistence of direction or internal bias we mean that the probability a particle will travel in a given direction need not be the same for all directions, but depends solely upon the particle's previous direction of motion. The external bias arises from an anisotropy of the medium or an external force on the particle. The problem is treated by considering that the net displacement of a particle arises from two factors, namely, that neither the probability of the particle traveling in any direction after turning nor the distance the particle travels in a given direction need be the same for all directions. A modified Fokker-Planck equation is first obtained using the assumptions that the particles have a distribution of travel times and speeds and that the average time of travel between turns need not be zero. The fional equation incopporating the assumption of a persistence of direction and an external bias is then derived. Applications to the study of diffusion and to long-chain polymers are then made.

Notizen: Abstract The transmission of some information or behavior pattern is treated as a flow of “particles” which execute random motions over a population of individuals and which may multiply or disappear. Equations are derived for the number density of these “particles” and from this is calculated the number of individuals through which the “particles” have passed. The results are applied to a number of situations such as 1) uniform spatial distribution with multiplication factor decreasing with time because of loss of interest or confusion of the information, 2) multiplication factor constant but the rate of spreal decreasing with multiple hearings, 3) one-dimensional region with a small starting region with or without an absorbing barrier 4) two-dimensional region with absorbing barrier, 5) continous sources of information within a small region in one dimension, 6) uniform spatial distribution in which individuals do not respond to more than one hearing.

Notizen: Abstract A situation is considered in which a fluid containing a substance flows through a vessel at a constant rate, the substance being permeable to the vessel wall. In the region outside the vessel there is supposed to be rapid mixing in the direction perpendicular to the axis of the vessel but no mixing longitudinally. The solution for the spatial distribution at any time is given for the case of an arbitrary initial distribution along the vessel length in the absence of an input. The solution is also given for the case of a single impulsive input, the concentration being initially zero everywhere.

Notizen: Abstract Organisms orient themselves to a stimulus by two general methods. One method is by directed orientation (taxis); the other is by undirected locomotory reaction (kinesis). An equation, and the methods for finding the necessary parameters of this equation, is derived for the distribution of organisms within a container, with the following limitations: (1) the organisms have no accommodation, (2) they are always active, and (3) the stimulus changes slowly with position. Necessary modifications of the equation are then derived, so that the last two limitations may be eliminated. The equation cannot be solved excatly because of its complexity; hence an approximation method must be used. This method is discussed, an approximate solution is found, and a time constant for equilibrium to be established is derived. Applications tovarious experiments in the literature are then made with fairly satisfactory results. A new interpretation of the theory of klino-kinesis with accommodation is found upon application of the equations developed to experimental work. Further limitations and uses of these equations are then discussed.

Notizen: Abstract Certain parameters are defined which roughly characterize the internal structure of networks. A given network structure uniquely determines the values of the parameters, but the reverse is not true. The parameters therefore define certain classes of networks. One of the parameters, thedispersion D(S) gives an indication of the “compactness” of the internal structure. Addition theorems and inequalities are derived relating the dispersions of sub-systems to the dispersion of the complete structure.

Notizen: Abstract A mathematical theory is developed which permits the determination of certain parameters of an inhomogenous tissue, such as a nerve trunk without its epineurium. The parameters are the permeability coefficients for entrance into an exit of a substance from the nerve fibers, and the diffusion coefficient of the interstitial material. The experimental data required are the dimensions of the cross-section, the average diameter of the fibers, and the ratio of the cross-sectional are of the fibers to the total cross-section, as well as the time course of the decrease of the fraction of the substance left in the nerve trunk, when the trunk is immersed in a bathing solution containing none of it.

Notizen: Abstract A model based on enzyme localization is developed which gives rise to an apparent active transport of a metabolite into or out of cells. The model is applied to three simple situations, using Fick's equation and the Rashevsky approximation. It is shown that the apparent efficiency can be made as large as desired if, for constant reaction, the outer cell region is made sufficiently small, or, for autocatalytic reaction, if the metabolite concentration in the outer region is sufficiently small. The physical limitations imposed by this mechanism are developed for all three situations.

Notizen: Abstract A previously derived iteration formula for a random net was applied to some data on the spread of information through a population. It was found that if the axon density (the only free parameter in the formula) is determined by the first pair of experimental values, the predicted spread is much more rapid than the observed one. If the successive values of the “apparent axon density” are calculated from the successive experimental values, it is noticed that this quantity at first suffers a sharp drop from an initial high value to its lowest value and then gradually “recovers”. An attempt is made to account for this behavior of the apparent axon density in terms of the “assumption of transitivity”, based on a certain socio-structural bias, namely, that the likely contacts of two individuals who themselves have been in contact are expected to be strongly overlapping. The assumption of transitivity leads to a drop in the apparent axon density from an arbitrary initial value to the vicinity of unity (if the actual axon density is not too small). However, the “recovery” is not accounted for, and thus the predicted spread turns out to beslower than the observed.

Notizen: Abstract The assumption of transitivity treated in part I is modified in various ways to describe an information-diffusion process, in which a certain amount of randomness of contact does occur. In one model a parameter is introduced which is indicative of a tendency to go beyond one's immediate vicinity to spread the information as the vicinity becomes saturated with knowers. In another model the randomness appears in the assumption that new knowers are uniformly distributed among the knowers. Two of the equations thus derived, each with two free parameters are in good agreement with experimental results.

Notizen: Abstract In the framework of the neural network theory effects similar to hypnotic displays are constructed. They are based on the associative paradigm involving non-linear interaction of excitatory and inhibitory channels with synaptic memory. The non-linearity of long-term memorizing processes may cause effects exhibited by blind spots, which are interpreted as the first stage of hypnosis. More complicated phenomena are discussed in terms of a two-layer network.

Notizen: Abstract Mutation is introduced into autocatalytic reaction networks. The differential equations obtained are neither of repliator-type nor can they be transformed straightway into a linear equation. Examples of low dimensional dynamical systems —n=2, 3 and 4 — are discussed and complete qualitative analysis is presented. Error thresholds known from simple replication-mutation kinetics with frequency independent replication rates occur here as well. Instead of cooperative transitions or higher order phase transitions the thresholds appear here as supercritical or subcritical bifurcations being analogous to first-order phase transitions.

Notizen: Abstract The non-linear behavior of a differential equations-based predator-prey model, incorporating a spatial refuge protecting a consant proportion of prey and with temperature-dependent parameters chosen appropriately for a mite interaction on fruit trees, is examined using the numerical bifurcation code AUTO 86. The most significant result of this analysis is the existence of a temperature interval in which increasing the amount of refuge dynamically destabilizes the system; and on part of this interval the interaction is less likely to persist in that predator and prey minimum population densities are lower than when no refuge is available. It is also shown that increasing the amount of refuge can lead to population outbreaks due to the presence of multiple stable states. The ecological implications of a refuge are discussed with respect to the biological control of mite pests.

Notizen: Abstract In many applications of control theory on plant growth models biomass maximization is postulated to avoid analytically unsolvable problems while fruit maximization is commonly considered to be a more realistic criterion. In a special case, we are able to compare these criteria. Iwasa and Roughgarden (1984,Theor. Pop. Biol. 25, 78–105) have investigated a certain class of plant growth models using a fruit maximization criterion. They proved that, in the vegetative growth period, the organs follow a certain path of balanced growth. We show that this path remains optimal when biomass maximization is postulated. This underlines the importance of the balanced growth path found by Iwasa and Roughgarden. Furthermore, our result suggests that in the vegetative growth period the biomass maximization criterion is a good approximation of fruit maximization. In another theoretical control investigation, Schultzeet al. (1983,Oecologia 58, 169–177) derived a different type of balanced growth path. We apply the theory of Iwasa and Roughgarden to an improved version of the model of Schulzeet al. This leads to a new description of balanced growth between root and shoot that reflects non-linearities in the water uptake process and constitutes an interesting hypothesis for further experimental testing.

Notizen: Abstract In this paper the effects of changing the ion concentration in and around a sample of soft tissue are investigated. The triphasic theory developed by Laiet al. (1990,Biomechanics of Diarthrodial Joints, Vol. 1, Berlin, Springer-Verlag) is reduced to two coupled partial differential equations involving fluid ion concentration and tissue solid deformation. These equations are given in general form for Cartesian, cylindrical and spherical geometries. After solving the two equations quantities such as fluid velocity, fluid pressure, chemical potentials and chemical expansion stress may be easily calculated. In the Cartesian geometry comparison is made with the experimental and theoretical work of Myerset al. (1984,ASME J. biomech. Engng,106, 151–158). This dealt with changing the ion concentration of a salt shower on a strip of bovine articular cartilage. Results were obtained in both free swelling and isometric tension states, using an empirical formula to acount for ion induced deformation. The present theory predicts lower ion concentrations inside the tissue than this earlier work. A spherical sample of tissue subjected to a change in salt bath ion concentration is also considered. Numerical results are obtained for both hypertonic and hypotonic bathing solutions. Of particular interest is the finding that tissue may contract internally before reaching a final swollen equilibrium state or swell internally before finally contracting. By considering the relative magnitude, and also variation throughout the time course of terms in the governing equations, an even simpler system is deduced. As well as being linear the concentration equation in the new system is uncoupled. Results obtained from the linear system compare well with those from the spherical section. Thus, biological swelling situations may be modelled by a simple system of equations with the possibility, of approximate analytic solutions in certain cases.

Notizen: Abstract Many models of immune networks have been proposed since the original work of Jerne [1974,Ann. Immun. (Inst. Pasteur) 125C, 373–389]. Recently, a limited class of models (Weisbuchet al., 1990,J. theor. Biol. 146, 483–499) have been shown to maintain immunological memory by idiotypic network interactions. We examine generalizations of these models when the networks are both large and highly connected to study their memory capacity, i.e. their ability to account for immunization to a large number of random antigens. Our calculations show that in these minimal models, random connectivities with continuously distributed affinities reduce the memory capacity to essentially nil.

Notizen: Abstract A kinetic model is proposed to delineate the factors that determine the coronary reactive hyperemic response (RHR) to transient ischemia. The model comprises of myocardial-interstitial (M) and vascular (V) compartments. Vasodilator metabolites (VM) are produced in the M compartment during the interval of coronary occlusion. The rate of VM production is dependent on the flow rate during the ischemic period, the ratio of excess flow above the control level (R) to the loss of flow during occlusion period (D), the amount of oxygen stored and the degree of vasodilation in the V compartment prior to occlusion. Following a complete release of occlusion, VM are transported from the M to V compartment and are washed out or degraded with time. The time course of RHR is determined by the coronary patency which is proportional to VM concentration in the V compartment. Based on a set of numerical constants, the model is tested by simulating RHR to the various occlusion manoeuvres: a pair of 10 sec occlusions separated by brief release, a 15 sec release followed by a second brief occlusion, a brief release of an occlusion followed by restriced inflow and a period of restricted inflow after occlusion. The simulated results fit the experimental R/D and RH durations data of canine hearts. Factors that determine the impairment of RH capacity in coronary stenosis are suggested in terms of the model scheme.

Notizen: Abstract In the present paper a kinetic study is made of the behaviour of a Michaelis-Menten enzyme-catalysed reaction in the presence of irreversible inhibitors rendered unstable in the medium by their reaction with the product of enzymatic catalysis. A general mechanism involving competitive, non-competitive, uncompetitive and mixed irreversible inhibition with one or two steps has been analysed. The differential equation that describes the kinetics of the reaction is non-linear and computer simulations of its dynamic behaviour are presented. The results obtained show that the systems studied here present kinetic co-operativity for a target enzyme that follows the simple Michaelis-Menten mechanism in its action on the substrate, except in the case of an uncompetitive-type inhibitor.

Notizen: Abstract The relative contributions of mitochondrial β-oxidation and peroxisomal β-oxidation and peroxisomal ω-oxidation to the oxidation of a given fatty acidin vivo can be quantitated by an isotopic method. The approach requires infusion of a fatty acid labelled on two specific carbon atoms (e.g. [1-14C] and [11-14C] palmitate) to an isotopic steady state, with subsequent isolation and degradation of an acetylated conjugate as a product of the liver cytosolic acetyl CoA pool and of ketone bodies as a product of the liver mitochondrial acetyl CoA pool.

Notizen: Abstract Pancreatic β-cells in intact islets of Langerhans perfused with various glucose concentrations exhibit periodic bursting electrical activity (BEA) consisting of active and silent phases. The fraction of the time spent in the active phase is called the plateau fraction and appears to be strongly correlated with the rate of release of insulin from islets as glucose concentration is varied. Here this correlation is quantified and a theoretical development is presented in detail. Experimental rates of insulin release are correlated with “effective” plateau fractions over a range of glucose concentrations. There are a number of different models for BEA in pancreatic β-cells and a method is developed here to quantify the dependence of a glucose dependent parameter on glucose concentration. As an example, the plateau fractions computed from the Sherman-Rinzel-Keizer model are matched with experimental plateau fractions to obtain a relationship between the model's glucose-dependent parameter, β, and glucose concentration. Knowledge of the relationships between β and glucose concentration and between experimental measurements of rates of insulin release and plateau fractions permits the determination of theoretical rates of insulin release from the model.

Notizen: Abstract When a suspension of bacterial cells of the speciesBacillus subtilis is placed in a chamber with its upper surface open to the atmosphere complex bioconvection patterns are observed. These arise because the cells: (1) are denser than water; and (2) usually swim upwards, so that the density of an initially uniform suspension becomes greater at the top than the bottom. When the vertical density gradient becomes large enough, an overturning instability occurs which ultimately evolves into the observed patterns. The reason that the cells swim upwards is that they are aerotactic, i.e. they swim up gradients of oxygen, and they consume oxygen. These properties are incorporated in conservation equations for the cell (N) and oxygen (C) concentrations, and these are solved in the pre-instability phase of development whenN andC depend only on the vertical coordinate and time. Numerical results are obtained for both shallow- and deep-layer chambers, which are intrinsically different and require different mathematical and numerical treatments. It is found that, for both shallow and deep chambers, a thin boundary layer, densely packed with cells, forms near the surface. Beneath this layer the suspension becomes severely depleted of cells. Furthermore, in the deep chamber cases, a discontinuity in the cell concentration arises between this cell-depleted region and a cell-rich region further below, where no significant oxygen concentration gradients develop before the oxygen is fully consumed. The results obtained from the model are in good qualitative agreement with the experimental observations.

Notizen: Abstract We describe a classification scheme for bursting oscillations which encompasses many of those found in the literature on bursting in excitable media. This is an extension of the scheme of Rinzel (inMathematical Topics in Population Biology, Springer, Berlin, 1987), put in the context of a sequence of horizontal cuts through a two-parameter bifurcation diagram. We use this to describe the phenomenological character of different types of bursting, addressing the issue of how well the bursting can be characterized given the limited amount of information often available in experimental settings.

Notizen: Abstract To ensure its sustained growth, a tumour may secrete chemical compounds which cause neighbouring capillaries to form sprouts which then migrate towards it, furnishing the tumour with an increased supply of nutrients. In this paper a mathematical model is presented which describes the migration of capillary sprouts in response to a chemoattractant field set up by a tumour-released angiogenic factor, sometimes termed a tumour angiogenesis factor (TAF). The resulting model admits travelling wave solutions which correspond either to successful neovascularization of the tumour or failure of the tumour to secure a vascular network, and which exhibit many of the characteristic features of angiogenesis. For example, the increasing speed of the vascular front, and the evolution of an increasingly developed vascular network behind the leading capillary tip front (the brush-border effect) are both discernible from the numerical simulations. Through the development and analysis of a simplified caricature model, valuable insight is gained into how the balance between chemotaxis, tip proliferation and tip death affects the tumour's ability to induce a vascular response from neighbouring blood vessels. In particular, it is possible to define the success of angiogenesis in terms of known parameters, thereby providing a potential framework for assessing the viability of tumour neovascularization in terms of measurable quantities.

Notizen: Abstract The paper presents a simple method for solving the optimal (LQG) control problem on an infinite time horizon for linear plants described by proper rational transfer function matrices. Since the class of proper plants discussed in the paper is more general than the one commonly used, additional properties of solutions are presented. Two numerical examples illustrate the theoretical considerations. The examples have been performed on an IBM PC using the proposed algorithm. This algorithm is a part of a public ASMCS package developed by the author for analysis, synthesis, and simulations of multi-input, multi-output (MIMO) control systems.

Notizen: Abstract We study the behavior of a Hilbert network (i.e., a finite or countably infinite network whose variables are in a Hilbert space and in which the associated total energy is finite) whose elements are affected by perturbations. More specifically, we will give estimates for a change of the current distribution caused by (a) perturbations of the elements of the nominal network when the voltage sources are fixed, and (b) a change of voltage sources in a network whose elements are perturbed. The conditions given in our theorems imply insensitivity and robust stability of the nominal network. The applications of the results are illustrated by an example of an infinite network.

Notizen: Abstract A new lattice filter structure to model two-dimensional (2-D) autoregressive (AR) fields is proposed. The proposed structure utilizes and extracts the information contained in the backward prediction error fields and their delayed versions. The main idea is to use two sets of reflection coefficients corresponding to two quadrant filters and to increase the number of reflection coefficients with the order of the lattice filter. Increasing the number of reflection coefficients at each stage produces a sufficient number of independent parameters to model AR fields up to order three, which is an improvement over the existing 2-D lattice filter structures. The improvement is confirmed by computer simulations. In addition, a relationship between the reflection coefficients and the AR coefficients is derived. It is also shown that the entropy contained in the backward prediction error field vector of the proposed structure is closer to the input entropy when compared to those contained in existing 2-D lattice filters.

Notizen: Abstract This paper addresses the problem of global asymptotic stability in recursive first hyperquadrant causalm-D digital filters. A set of simple-to-check 1-D conditions necessary for stability of them-D system is given. Generalizations tov-dimensional (v〈m) subsystems are also provided. These derived conditions are useful in determining asymptotic convergence ofm-D digital filters implemented in fixed or floating point arithmetic. The set of 1-D conditions can be considered the analogous result to practical bounded input bounded output (BIBO) stability for linearm-D systems.

Notizen: Abstract A new efficient algorithm for calculating the weighting coefficients of maximally linear, FIR digital differentiators is presented. Simple closed-form explicit and recursive formulas are derived in a very straightforward manner. Moreover, a simple recursive equation is established, relating coefficients of two digital differentiators of adjacent ranks.

Notizen: Abstract This paper presents a sequence of pseudorandom arrays with triangular symmetry. The sequence of arrays is also pseudorandom in nature, so it can be called a pseudorandom sequence of arrays, or PRSA. A circuit for the PRSA generator is given and some interesting properties of this type of PRSA are discussed. Also, some of the concepts presented in this paper are clarified with an example.

Notizen: Abstract We investigate the $L_p$ -error of approximation to a function $f\in L_p({\Bbb T}^d)$ by a linear combination $\sum_{k}c_ke_k$ of $n$ exponentials $e_k(x):= e^{i\langle k,x\rangle}=e^{i(k_1x_1+\cdots+k_dx_d)}$ on ${\Bbb T}^d,$ where the frequencies $k\in {\Bbb Z}^d$ are allowed to depend on $f.$ We bound this error in terms of the smoothness and other properties of $f$ and show that our bounds are best possible in the sense of approximation of certain classes of functions.

Notizen: Abstract Let $L[\,\cdot\,]$ be a nondivergent linear second-order uniformly elliptic partial differential operator defined on functions with domain $\Omega.$ Consider the question, "When is a function u a solution of $L[u] = 0$ on $\Omega$ ?" The naive answer, "u is a solution of $L[u] = 0$ on $\Omega$ if $u\in C^2(\Omega)$ and $L[u](x) = 0$ for all $x\in\Omega,$ " is clearly too limited. Indeed, if the coefficients of L are in $W^{1,2}\cap L^{\infty},$ then L can be rewritten in divergence form for which the notion of a "weak" solution can be applied. In this case there could be infinitely many functions that are "weak" but not classical solutions. More importantly, even if the coefficients of L are just bounded and measurable, the recent results of Krylov permit us to construct "solutions" of $L[u] = 0$ on $\Omega,$ and these "solutions" are generally no better than continuous; the "weak" solutions previously mentioned can be obtained by this construction, too. The preceding discussion provides us with an adequate extrinsic definition of solution (i.e., given a function u we either prove that it is or is not the result of such a construction) that has been used by several authors, but one that is not particularly satisfying or illuminating. Our major contribution in this paper is to show the following. I. There is an intrinsic definition of solution that is equivalent to the extrinsic one. II. Furthermore, the intrinsic definition is just the (now) well-known Crandall-Lions viscosity solution, modified in a natural way to accommodate measurable coefficients.

Notizen: Abstract We prove a Tauberian theorem of the form $\phi * g (x)\sim p(x)w(x)$ as $x \to \infty,$ where p(x) is a bounded periodic function and w(x) is a weighted function of power growth. It can be used to study the weighted average of the form $(T^\alpha (\hbox {ln }T)^\beta)^{-1}\int _0^T h(t) \, dt.$

Notizen: Abstract In this paper, we present a new class of two-dimensional FIR-median hybrid (FMH) filters, which we call separable FMH filters, for image noise smoothing. In a separable FMH filter, a one-dimensional FMH filter is applied to the rows and the columns of an image successively. The deterministic properties of separable and cross window FMH filters are discussed. Under certain assumptions, it is proved that a root of a separable FMH filter is a root of the corresponding cross window FMH filter. The noise attenuating properties of a separable FMH filter are studied and compared with those of the separable median filter, the two-dimensional median filter with a square window, and the cross window FMH filter.

Notizen: Abstract The Chinese remainder theorem is a fundamental technique widely employed in digital signal processing for designing fast algorithms for computing convolutions. Classically, it has two versions. One is over a ring of integers and the second is over a ring of polynomials with coefficients defined over a field. In our previous papers, we developed an extension to this well-known theorem for the case of a ring of polynomials with coefficients defined over a finite ring of integers. The objective was to generalize number-theoretictransforms, which turn out to be a special case of this extension. This paper focuses on the extension of the Chinese remainder theorem for processing complex-valued integer sequences. Once again, the present work generalizes the complex-number-theoretic transforms. The impetus for this work is provided by the occurrence of complex integer sequences in digital signal processing and the desire to process them using exact arithmetic.

Notizen: Abstract A method to improve the stability of the forced response of a recursive digital filter with nonlinearities due to finite wordlength is presented. An analysis of applying a stability criterion on the nonlinear filter is performed by transforming the problem into a mathematical problem of finding a zero set. As a result of the mathematical analysis, one can find a maximal sector size that bounds an equivalent nonlinearity, so that stability is ensured. This sector size is influenced by the choice of the desired nonlinearity characteristic appropriate to the problem, and the amount of the input scaling needed. Another improvement of the filter's stability is obtained by adding a feedback gain, the value of which is determined by using the zero set-finding method. Simulation results showing the improved stability properties of a filter designed by this method are presented.

Notizen: Abstract In this paper a fully systolic, bit level, three-port unconstrained adaptor, which constitutes the main nontrivial building block of ladder wave digital filters, is generated. The one-dimensional binary convolution is used as the underlying algorithm for the implementation of a multiplication. The Isb-first input data organization approach is adopted and thecanonical mapping methodology is used to fully systolize the unconstrained parallel three-port adaptor at the bit level with piplining period a=1. The technique is based on a transformation of the adaptor's signal-flow graph, so that unidirectional data flow takes place. A ring-systolic scheme is proposed for implementing communications among adaptors.

Notizen: Abstract We consider the problem of robust stochastic adaptive control of not necessarily minimum phase systems in the presence of unmodelled dynamics. Stochastic gradient algorithms with parameter projection and modified gain sequence are used for the estimation of the unknown controller parameters. Global stability of the adaptive system is achieved without requiring the strictly positive real condition and the persistency exciting condition to be satisfied.

Notizen: Abstract A theoretical framework for the investigation of the qualitative behavior of differential-algebraic equations (DAEs) near an equilibrium point is established. The key notion of our approach is the notion of regularity. A DAE is called regular locally around an equilibrium point if there is a unique vector field such that the solutions of the DAE and the vector field are in one-to-one correspondence in a neighborhood of this equilibrium point. Sufficient conditions for the regularity of an equilibrium point are stated. This in turn allows us to translate several local results, as formulated for vector fields, to DAEs that are regular locally around a given equilibrium point (e.g. Local Stable and Unstable Manifold Theorem, Hopf theorem). It is important that these theorems are stated in terms of the given problem and not in terms of the corresponding vector field.

Notizen: Abstract The classical notion of the λ-generalized nullspace, defined on a matrixA εR n×n,where λ is an eigenvalue, is extended to the case of ordered pairs of matrices(F, G), F, G ε R m×nwhere the associated pencilsF − G is right regular. It is shown that for every α εC ∪ {∞} generalized eigenvalue of (F, G), an ascending nested sequence of spaces {P α i ,i=1, 2,...} and a descending nested sequence of spaces {ie495-02 i=1, 2,...} are defined from the α-Toeplitz matrices of (F, G); the first sequence has a maximal elementM α * , the α-generalized nullspace of (F, G), which is the element of the sequence corresponding to the index τα, the α-index of annihilation of (F, G), whereas the second sequence has the first elementP α * as its maximal element, the α-prime space of (F, G). The geometric properties of the {M α i ,i=1, 2,...,τα and {P α i ,i=1, 2,...sets, as well as their interrelations are investigated and are shown to be intimately related to the existence of nested basis matrices of the nullspaces of the α-Toeplitz matrices of (F, G). These nested basis matrices characterize completely the geometry ofM α * and provide a systematic procedure for the selection of maximal length linearly independent vector chains characterizing theα-Segre characteristic of (F, G).

Notizen: Abstract The robustness problem of stability for large-scale uncertain systems with a class of multiple time delays is addressed in this paper. By applying the complex Lyapunov stability theorem, the matrix measure techniques, and norm inequalities, a new approach for solving a general case of the above problem is proposed. Several robust stability conditions, delay-dependent or delay-independent, are derived to guarantee the asymptotic stability and exponential stability of the uncertain large-scale time-delay systems. Moreover, these obtained results can also be applied to the stabilization design.

Notizen: Abstract Successful speech recognition is highly dependent on appropriate speech segmentation. The poor efficiency of the sequential detection of abrupt changes in the signals with relatively short stationary intervals, as is the case with speech signals, can be improved by the off-line maximum likelihood segmentation algorithm. In this paper the new segmentation algorithm is presented. For the a priori known number of segments, the algorithm determines such signal partitions for which the sum of segment distortion is minimal. The generalized maximum likelihood distortion measure has been introduced, and has proven to be particularly efficient on short signal segments. In the case of an unknown number of segments, its estimate is obtained comparing the reduction of the distortion. The asymptotic properties of the distortion sequence have been analyzed, which led to the definition of the presented segmentation algorithm. The introduced measure can be applied both to the AR and ARMA models. The segmentation algorithm is verified on test signals as well as on the natural speech signal, for which the pitch synchronous framing scheme is applied. The experimental results also include a comparison of the AR and ARMA model-based segmentations. The first results show that ARMA model-based segmentation gives somewhat better results than the AR model algorithm.

Notizen: Abstract A new bit-serial architecture for implementation of high order FIR filters is introduced, as well as example FPGA and CMOS realizations. This structure exploits the simplicity of coefficients that consist of two power-of-two terms to yield efficient implementations. Quantization effects are discussed and a simple block scaling method for reducing rounding and truncation noise in high order filters is also presented.

Notizen: Abstract This paper describes an implementation of 2-D FIR and IIR linear digital filters via VLSI array processors. The underlying realization structures are based on the matrix decomposition approach. The 2-D concurrent processing is used in order to implement the row and column delays within the cycle time. A high degree of concurrency is achieved by exploiting the pipelining of the array processors with the inherent parallelism of the matrix decomposition structure. The resulting structures are modular, and regular, use only local communication and internal local feedback loops, and achieve high throughput and sampling rates.

Notizen: Abstract This article studies dynamical systems under perturbations. We prove an Equilibrium Equivalence Theorem that guarantees that the dynamics of the system remains unchanged under perturbations with certain fairly general assumptions. We also prove that a power system is robust with respect to parameter changes in generic situations. Concepts of equilibrium equivalence and equilibrium equivalence structural stability are developed and are applied to studies of bifurcations of vector fields on noncompact manifolds. A constructive approach to equilibrium equivalence structural stability verification is emphasized. General results on structural stability of vector fields on differential manifolds are established and important applications of this theory to stability analysis are considered.

Notizen: Abstract The problem of designing a stabilizing compensator for a control system to achieve prescribed initial value constraintsυ (i)(0+)=yi is considered. Indeed, modulo certain technical conditions, such a compensator exists if and only if yi=0;i= 0,1,...,r〈p〉 +r〈t〉 −2; wherer〈p〉 is the relative degree of the plant andr〈t〉 is the relative degree of the system input. This theorem is derived and a complete parameterization of the set of compensators that achieve the prescribed design constraints is formulated.

Notizen: Abstract Wideband source location in array signal processing has received much attention in the literature lately. Methods such as the Coherent Signal Subspace (CSS) method proposed by Wang and Kaveh [12], and the signal subspace method used by Cadzow [3], are typical of the approaches used to tackle the multiple wideband source location problem. Most of these methods are variations of the narrowband high-resolution methods. Grenier [5], on the other hand, has applied the idea of time-dependent Auto-Regressive (AR) modeling [7] for a nonstationary process to the frequency domain AR modeling of the sensor outputs in a linear array and has been able to produce good results for a wideband signal. The AR coefficients in the model are expanded in a set of frequency-dependent basis functions. The choice of the basis functions was deemed immaterial and the method works even when only one snapshot of the array output is available. In this paper, we re-examine this method and present an extension of the frequency-dependent AR modeling approach to a planar array. It is shown that the use of a set of sinc functions for representing the frequency-dependent AR coefficients accurately tracks their evolution in the frequency domain, and gives superior performance compared to that when power or Legendre functions are used. We also propose two methods for smoothing the spatial spectra, from which the source locations are determined. Comparison with the CSS method are also presented.

Notizen: Abstract Centralized methods for source location using sensor arrays have computational and communication burdens that increase significantly with the number of sensors in the array. Therefore, these methods may not be usable in the applications involving very large arrays. In such applications, the data processing may need to be decentralized. This paper introduces two methods for decentralized array processing, based on the recently proposed MODE algorithm. For prescribed nonoverlapping subarrays, both methods are shown to be statistically optimal in the sense that asymptotically they provide the most accurate decentralized estimates of source location parameters. The problem of subarray selection to further optimize the estimation accuracy is only briefly addressed. The two methods are intended for different types of applications: the first should be preferred when there exist significant possibilities for local processing or for parallel computation in the central processor; otherwise the second method should be preferred. The accuracy of the two decentralized methods is compared to the centralized Cramér-Rao bound, both analytically and numerically, in order to provide indications about the loss of accuracy associated with decentralized processing.

Notizen: Abstract We present in this paper a recursive-in-order least-squares (LS) algorithm to compute efficiently the parameters of a 2-D Gaussian Markov random field (GMRF) model. The algorithm is based on the fact that the least-squares estimation of the parameters of a 2-D noncausal GMRF model is consistent and the coefficient matrix in the normal equation has near-to-block-Toeplitz structure. Hence, it has a Levinson-like form for the updating of model parameters by introducing auxiliary variables. Moreover, this paper proposes the concept ofrecursive path for 2-D recursive-in-order algorithms, and points out that there exists a tradeoff between fast computation of the parameters and accurate choice of model support; a compromise recursive path is then suggested where the orders change alternately in two directions. The computational complexity of the developed algorithm is analyzed, and the results show that the algorithm is more efficient when either the image size or the model support is larger. It is found that the total number of multiplications (mps) involved in the new algorithm is only about 14% of that in the conventional LS method when the image size is 512 × 512 and the neighbor set of the model is a 17 × 17 window. Computer simulation results using the recursive-in-order algorithm developed in this paper and the conventional LS method are given to verify the correctness of the new algorithm.

Notizen: Abstract Cross terms are an inherent consequence of the second order nature of Cohen's class TFDs (Time-Frequency Distributions) [5], [6]. They are manifest in a TFD of multicomponent signals as spurious artifacts arising from interactions between the various signal components, and they can often appear at times and/or frequencies inconsistent with the underlying physical nature of the signal, causing misinterpretation [2], [3], [4]. There are many time frequency distributions that avoid the cross term effect; the best are the Choi-Williams ED (Exponential Distribution) [1] and Levin's IPS (Instantaneous Power Spectrum) [9]. In this paper we combine the cross term reducing philosophy of the ED and IPS to obtain a new TFD that most effectively reduces the cross term effect. Surprisingly, the new TFD also satisfies most desired TFD properties.

Notizen: Abstract The impressed field, “Window Field” (WF), due to a half-wave action potential on a muscle fiber, has been calculated on the basis of potential theory. It has been shown that in spite of the small intensity of the field, its integrated action can transfer the energy needed to induce, contraction from the membrane to the interior of the fiber. The energy of polarization has been found to be sufficient to exceed the energy of, thermal agitation on that length of fiber, which can be identified as the length of a sarcomere. The changes of ion concentration, caused by the WF, if calculated on the assumption of the semipermeability of theZ membranes, was found to be equal to the changes necessary to induce contraction of actomyosinin vitro.

Notizen: Abstract Some general properties of the solution of the diffusion equation are deduced for the steady-state, spherically symmetric system. On the basis of these developments some results of N. Rashevsky (Bull. Math. Biophysics,11, 15, 1949) are discussed and the results of a previous investigation (Hearon,Bull. Math. Biophysics,12, 135, 1950b) are extended to more general conditions. In particular these extensions apply to the flow of a soluteagainst its concentration gradient, the nonzero gradient of an inert metabolite, and theaccumulation or exclusion of an inert metabolite in a metabolic system.

Notizen: Abstract The approximation method of N. Rashevsky is discussed and reviewed. It is shown that in addition to theexplicit assumptions and approximations there is involved the assumption that the rate of metabolism is the same at every point in the cell and that theaverage rate of metabolism is different from zero. An expression is given for the error in the approximate method when the rate of metabolism is any function of the concentration. It is also shown that a solution in theform of that obtained by the approximate method is not possible if the generalized laws of diffusion are assumed to apply.

Notizen: Abstract Rate equations for the enzymatic oxidation of succinic acid are derived on the assumption that when a single molecule of substrate combines with an enzyme molecule, it can do so with either one or two sites on the enzyme, and that oxidation occurs only in the second case. In addition it is assumed that the product of the reaction, fumaric acid, combines reversibly with the enzyme. With certain enzyme preparations the data fitted such an equation satisfactorily. In others the rate was that of a first-order reaction, but addition of cytochrome changed it to the former type. It was concluded that the transfer of hydrogen to oxygen was a first-order reaction and dominated the whole rate when enzyme preparations were used which had been washed relatively free of cytochrome. When the limiting factor was succino-dehydrogenase the rates followed the new equation. Criteria for recognizing noncompetitive inhibition are given, and inhibition by di-tertiary butyl peroxide was shown to be of this type.

Notizen: Abstract This paper deals with the following question: Which distributions of radiosensitivity in a population can lead to an exponential survival curve? The problem is solved exactly, with statistical fluctuations in dose fully accounted for. It is shown that only an exponential distribution of sensitivities can give rise to an exponential survival curve.

Notizen: Abstract An approximation method is introduced which enables a number of diffusion-type problems to be solved in an approximate but simple manner. Many cases require only the solution of a simple first-order differential equation. The method is applied to a number of cases in which the exact solutions are available. A comparison shows that the method is quite satisfactory in these cases. The method is applied to diffusion problems with rate of consumption proportional to concentration or to the square of the concentration. In the latter case, the result obtained is essentially the same as that found by H. G. Landau (1950) after elaborate calculations.

Notizen: Abstract It is shown that a slight modification of a model of excitatory phenomena in irritable tissues, which has been treated before, exhibits spontaneous oscillations. The frequency of these oscillations and the time-course of the potential across the model membrane have been determined, together with the dependence of some of their characteristics on some important parameters, particularly (Ca++).

Notizen: Abstract It is demonstrated that, if the variations of viscosity throughout a cell are considered, swelling stresses may produce elongation and division. To do this it is necessary to generalize Betti's theorem to cover systems containing viscosity gradients and such a generalization is presented. On the basis of two special assumptions it is shown that most of the results of the diffusion drag theory of cell division may be duplicated by the present theory.

Notizen: Abstract The theory of imitative behavior, developed previously, is applied to the case of two social groups which are separated spatially. If the information of each group as to the behavior of the other is complete, the case reduces to that of a single group. When any information is lacking at all, the two groups are independent. If we have two mutually exclusive behaviorsA andB, all four combinationsAA, AB, BA, andBB are possible. If the mutual information gradually increases from zero, then for a certain value of it, the group which is more informed about the behavior of the other will change to that behavior if it did not already exhibit it. If for constant information the size of the group increases, then above a certain threshold value, the larger group imposes its behavior on the smaller.

Notizen: Abstract Local regulation of blood flow as determined by capillary diameter and the number of open capillaries in a region is considered. The local changes in capillary diameter and in the number of open capillaries are assumed to be due to concentration changes of a diffusible, nonspecified metabolite. This metabolite is produced in the tissue and carried away by the blood stream. Using these assumptions and applying pertinent data on capillaries, deductions are made concerning: (a) the law of blood flow as a function of temperature and capillary radius for the hyperemia of high temperature, (b) high flow as it depends on metabolism during strenuous exercise of muscle, and (c) a first approximation to the time duration of occlusion hyperemia.

Notizen: Abstract The necessary and sufficient condition is given forn integers to be the score structure of a society with a dominance relation. A proof is also given for a theorem showing that there are members who dominate every other member either directly or indirectly through a single intermediate member.

Notizen: Abstract On the basis of a previous general formulation (Bull. Math. Biophysics,15, 21–29, 1953a) a discussion is given of the error in the approximation method of N. Rashevsky. This error, inherent in the method when the metabolic rate is different at each point in the cell, is discussed in detail and numerical values are presented for two particular cases: the rate proportional to the concentration and the rate a prescribed function of the spatial coordinates. It is shown that the formulation for the first case also applies to several other cases, that the error is negligible provided the rate is sufficiently small, and that the error is fairly sensitive to the cell size. If the rate depends upon the coordinatesalone a small rate is not sufficient to insure a negligible error. The relations between the exact method, the standard approximate method, an earlier approximate method (Physics,7 260, 1936), and a more recent refinement (Bull. Math. Biophysics,10, 201, 1948) of the standard method are discussed.

Notizen: Abstract It is shown on the basis of (1) conservation of mass, (2) positive concentrations, and (3) the principle of detail balancing that periodic reactions cannot occur in a closed system described bylinear differential equations. The matrix,A, of the rate equations must be such that |A|=0,a ij〉0 fori≠j,a ii〈0, andVAV −1=B, whereV is diagonal andB is symmetric. These properties ofA imply that the latent roots are real and non-positive and that neither catalysis nor inhibition can be described bylinear equations. It is further shown that periodic reactions cannot occur in anopen system for which the matrix associated with the chemical reactions has the above properties and in which thesimple law of diffusion is obeyed. The relation of these results to Onsager's reciprocal relations and to previous work on periodic and cyclic chemical reactions is discussed. The utility of certain of these results for the treatment of isotope kinetics is indicated.

Notizen: Abstract It is demonstrated that an explanation of the small radius effect or the so-called sigma phenomenon may be obtained by noting that one of the effects of the presence of suspended particles in a flowing fluid is to increase the velocity of flow near the wall over that existing in the absence of particles. This effect may be considered equivalent to relaxing the boundary conditions at the wall. An expression for the viscosity is compared with data and fit is found to be good.

Notizen: Abstract The solution for the spatial distribution of ions in a Donnan equilibrium has been given by J. H. Bartlett and R. A. Kromhout (1952). The present note gives an explicit solution for the case in which the length of the region containing the membrane is large; in biological situations this requires only that the length considered should be greater than a few hundred Ångstrom units. The Donnan equilibrium may be considered to be a special case of a situation in which forces other than electrical act upon the ions; in particular, it represents the case in which only one ion is acted upon and the energy difference on the two sides of the membrane is infinite. An expression is given for the difference in energy of theith in terms of the electrical potential and of the ion concentrations. As an illustration, the results are applied to nerve membrane potentials.

Notizen: Abstract A mechanism is described which accounts for the active transport of Na+ ions through a membrane. It is assumed that at one side of the membrane the ion combines with a carrier ion, the resulting carrier compound then diffuses through the membrane and decomposes at the other side of the membrane. The free diffusion of the ions is also taken into account. The time rate of accumulation of the ion in question at the latter side of the membrane is calculated in terms of the concentrations of the ion at both sides of the membrane.

Notizen: Abstract The recent extension of the approximation method is applied to enable us to arrive at the time course of the concentrations at both sides of a membrane. From the differential equations which govern these, the steady-state solution is obtained in terms of the parameters, which are determined by the thickness of the diffusion layers, the chemical composition and reactions, and the diffusion constant of the membrane.

Notizen: Abstract An equation is derived from the spread of a “state” by contact through a thoroughly mixed population, in which the probability of transmission depends both on the over-all duration of the process and on the time an individual has been in the “state.” Cases in which this probability is a function of only one or the other of the two “times” are worked out. It is shown that in the case of dependence on “private time” alone the asymptotic value of the fraction of the population effected is the same as that derived by the random net approach.

Notizen: Abstract The reflection of pressure waves in a fluid enclosed within a tube with an elastic wall is studied for the case of a localized change in diameter of the tube. The concept of impedance is introduced. The relation of the reflection characteristics of the parts of the tube at either side of the change is derived on the basis of the continuity of pressure and mass flow at the site of the change. This relations is used to derive the expression for the ratio of the pressure oscillations measured in front of, and behind, the constriction in terms of the constants of the system. As a result, a method is indicated to locate the coarctation from measurements of the pressures in front of, and behind it.

Notizen: Abstract The question is raised concerning the possible causes of abnormally small standard deviations found in counting samples in which particles are distributed at random (e.g., blood cells, fat globules in milk, etc.). The effect of discarding abnormal samples is discounted inasmuch as small standard deviations occur even when all samples are counted. An approximation method is used to calculate the effect of finite particle size, of known repulsive forces between particles and of convection currents. This calculation shows that neither finite size nor the known repulsive forces are sufficient to account for the observed abnormality of standard deviation, but that convection currents can possibly account for it. The possible presence of long-range repulsive forces cannot, however, be excluded.

Notizen: Abstract The effect of finite particle size on the standard deviation in sample counts is computed for the one-dimensional case. To a first order of approximation the correction is found to be identifical with that found by H. de Vries (1953) using a general approximation method.

Notizen: Abstract A critical examination of the “classical” theories of photoreception in view of more recent experimental findings yields the result that these theories do not possess the property to describe all the more significant phenomena of photoreception correctly, and to some extent suffer the lack of more general applicability. The basis for a new and presumably more general theory of photoreception based on dynamical aspects is laid out. Emphasis is put on the time course of afferent and efferent excitation in the photoreception model, consisting of a receptor element, an afferent and an efferent neuron of the one-factor Rashevsky-type, and an effector organ.

Notizen: Abstract A mathematical model for the development of human society, beginning with the earliest stages of urban cultures, is outlined. In the early stages of history, behavior was characterized largely by adherence to a number of beliefs and prejudices of diffeirent kinds, which were accepted on faith and not subject to critical rational analysis. Due to psychobiological variability a very small number of individuals spontaneously appear at all times who challenge the accepted beliefs and prejudices and do not follow the accepted patterns of social behavior. The effect of these individuals upon the rest of the society, especially upon the younger generation, depends on the facilities with which information spreads in society. In earliest societies, when modern methods of mass communication were unknown, the channels of communication were practically identical with the channels of economic transport. The latter in its turn depended on the nature of the roads, and especially on the presence of waterway, which facilitated transportation. The sizes of the earliest cities and the distances between them were largely determined by relative ease of transporation. Expressions are derived for the average size of the earliest cities and for the average distance between them. The calculated average populations of the earliest cities are of the order of 103; the distance of the order of 102 km. Both are in agreement with some archaeological findings. An expression for the time spaon required for the development from the earliest stages of urban cultures to the present time is derived and shown to depend on the specific shoreline of the country, that is, the length of the shorline divided by the area of the country. It is pointed out that western Europe's specific shoreline, including land bordering both seas and rivers, is ten times as large as the shoreline area of other parts of the world. It is shown that this greater specific shoreline may account quantitatively for the faster social and technological development of western Europe in the last few centuries. The calculated total span of time of development from earliest urban cultures to our days is found to be of the order of magnitude of ten thousand years. It is shown that the model accounts for the existence at the present time of primitive cultures. A number of suggestions is made in regard to other possible applications of mathematics to history.

Notizen: Abstract The assumptions latent in the derivation of the integral equation of Branson are rendered explicit and discussed. It is shown that the equation is valid only for systems in which the substance disappears according to a linear rate law.

Notizen: Abstract It is shown that the validity of Branson's integral description of metabolizing systems is subject to severe limitations. The validity is insured only in cases where the reaction is of first order, or quasi of first order. In all other cases Branson's equation has to be modified to insure general applicability. The consequences of a different definition of the metabolizing functionF have also been investigated. With the new definitionF describes the pure effect of metabolization. It is found that in this case the integral equation is only capable of describing first-order reactions. With a slight modification of the integral equation it is possible to describe metabolites “with age”, which do not have reactions of definite order, but which satisfy the superposition principle.

Notizen: Abstract A problem in probability is stated with included the problem of the distribution of bacterial mutants as a special case. This problem is solved exactly but since the resulting expressions are too complicated for practical use, various approximate expressions for the distribution are considered, especially for the bacterial mutation case.

Notizen: Abstract Simple reaction and discrimination reaction, under the influence of imitation, are considered for the situation in which the stimulus or the stimuli vary slowly with time. The result is analogous to hysteresis under certain conditions. The calculations are facilitated by the solution of $$x = \int_{ - \infty }^{a + \beta x} {g\left( \xi \right)d} \xi ,$$ g(ξ) being the normal error function. Values ofx(α, β) are given in a table.

Notizen: Abstract On the basis of simple physical considerations the blood flow in a branching circulatory system is studied. The case of two groups of parallel vessels is treated. The vessels of the same group are supposed to be identical. The resistance of each group is determined by the resistance of each vessel in the group and by the number of vessels in the group. From the dependence of the resistance of each vessel on its radius an expression is obtained for the blood flow through each group of vessels in terms of the numbers and sizes of the vessels in each group. The number of open vessels in an organ and the radius of each of those vessels are assumed to depend on the metabolic rate of that organ. The relations so obtained, together with the expression above, are applied to derive the blood flow through an organ as a function of the metabolic rate of that organ. It is indicated that the relations obtained might describe the shifting of blood from one organ to another if the activity of one of them changes. A way is pointed out to treat neural regulation of this phenomenon.

Notizen: Abstract A stochastic model of population growth is treated using the Bellman-Harris theory of agedependent stochastic branching processes. The probability distribution for the population size at any time and the expectation are obtained when it is assumed that there is probability (1−σ), 0≤σ〈1, of the organism dividing into two at the end of its lifetime, and probability σ that division will not take place.

Notizen: Abstract In a previous paper, in which a possible mathematical approach to history was outlined, it was shown that urbanization plays an important part in the propagation of new ideas. The rate of such propagation influences the rate of historical developments. The present paper deals in more detail with possible mechanisms of formation of earliest cities. Equations are derived which give the limiting size of such cities and their rate of growth. Of particular importance for the spread of new ideas is the spread of information. The latter largely depends on the fraction of individuals who travel between city and country. Expressions for this quantity are derived. An approach is outlined to the mathematical study of the earliest social classes, which may have been formed as a result of military, religious, or economic stratifications.

Notizen: Abstract The velocity of propagation of a disturbance wave in a liquid flowing in a distensible tube is computed. The mathematical model is more general than those used in previous analyses: the tube wall properties are realistic; the convective part of the axial inertia forces is taken into account; radial inertia forces of both the fluid and tube wall are present; viscous stresses are present. Four parameters influencing the velocity of propagation are obtained and discussed. Curves are plotted illustrating the effects of the parameters. Contrary to the results of previous analyses, viscous effects are shown to be appreciable in blood flow. It is also shown that radial inertia effects can be important in laboratory set-ups.

Notizen: Abstract A general theory of the drying of frozen tissue is developed and applied to the measurement of the drying rate of frozen guinea pig liver. It is shown that for a given temperature of the subliming ice crystals the mininum drying time of a piece of guinea pig liver is greater than the minimum sublimation time of a piece of ice of the same size and shape by a factor of the order of one thousand. This fact has many implications in the design of freeze-dry apparatus which will be developed in a following paper.

Notizen: Abstract The equations governing the time course of the exchange of substances between the blood in the capillaries and the extracellular space are solved for the case of substances which do not penetrate the cells. The equations given relate the time course of the exchange process to the various tissue and circulation parameters such as the specific capillary wall area, the pore area, the inter-capillary distance, the size of the extra-vascular, extra-cellular space, the diffusion coefficient in this space, and the velocity of blood in the capillaries. Some experimental work on capillary exchange is discussed in relation to the theory and estimates are made of the relative importance of the various tissue and circulation parameters in the exchange of substances in different tissues.

Notizen: Abstract The formation of neuronal networks requires axonal growth towards target neutons. A simple set of grammar rules is introduced to describe axonal growth towards target cells situated both at short and long distances from the growing neuron. Growth for short distances is descrbed by growth following the highest gradient of a chemical compound (which is spread by diffusion from the targets). This approach fails to describe long-distance growth, which is addressed by adopting a graph grammar theory for growing trees. With these rules a flexible tool to draw network of neurons by computer can be developed.

Notizen: Abstract A pair of growth control triads are used to describe coincident tumor growth and liver regeneration after partial hepatectomy. The models are extensions of previous growth control models which describe tumor growth in an unperturbed host (Michelson and Leith, 1991,Bull. math. Biol. 53, 639–656; idem, 1992, Proceedings of the Third International Conference on Communications and Control, Vol. 2, pp. 481–490; idem, 1992,Bull. math. Biol. 55, 993–1011; idem,J. theor. Biol. 169, 327–338). The linkage between the two triads depends upon systemic signals carried by soluble factors, and mathematical descriptors based upon biological first principals are proposed. The sources of the growth factors, their targets and the processing of their signals are investigated. Analyses of equilibrium in the constant coefficients case and simulated growth curves for the dynamic system are presented, and the effects of growth factor-induced mitogenesis and angiogenesis are discussed in particular. A case is made for early and late responses in the coupled control system. The biology of the signal processing paradigm is placed within a new theoretical context and discussed with regard to tumor adaptation, liver differentiation and the development of a tumor hypoxic fraction.

Notizen: Abstract We consider the problem of optimal stabilization and control of populations which follow the Leslie model dynamics, within state space and control systems theory and methodology. Various types of culling strategies are formulated and introduced into the Leslie model as control inputs, and their effect on global asymptotic stability is investigated. Our new approach provides answers to several unexplored problems. We show that in general it is possible to achieve a desired stable equilibrium population level, through the design of a class ofshifted-proportional stabilizing culling policies. Further, we formulate general non-linear constrained opitmization problems, for obtaining the cost-optimal policy among this generally infinite class of such stabilizing policies. The theoretical findings are illustrated through the solution of the problem over an infinite planning horizon for a numerical example. A comparative study of the costs and dynamic effects of various culling strategies also supports the mathematical results.

Notizen: Abstract A new approach for data assimilation, which is based on the adjoint method, but allows the computer code for the adjoint to be constructed directly from the model computer code, is described. This technique is straightforward and reduces the chance of introducing errors in the construction of the adjoint code. Implementation of the technique is illustrated by applying it to a simple predator-prey model in a model fitting mode. A series of identical twin numerical experiments are used to show that this data assimilation approach can successfully recover model parameters as well as initial conditions. However, the ease with which these values are recovered is dependent on the form of the model equations as well as on the type and amount of data that are available. Additional numerical experiments show that sufficient coefficient and parameter recoveries are possible even when the assimilated data contain significant random noise. Thus, for biological systems that can be described by ecosystem models, the adjoint method represents a powerful approach for estimating values for little-known biological parameters, such as initial conditions, growth rates, and mortality rates.

Notizen: Abstract The effect of group size on behavioral parameters of the Oriental hornet,Vespa orientalis, was assessed experimentally under laboratory conditions. Hornet groups of various sizes (ranging from 1 to 100 individuals per group) comprised of young individuals (0–24 hr of age) devoid of a queen were placed in artificial breeding boxes (ABBs). The following three quantitative parameters were evaluated: the amount and rate of building as a function of the number of hornets in the group, the rate of oviposition as, related to group size and the longevity of hornets as a function of their group size. The probability for the occurrence of these events was similarly considered and additional behavioral parameters were only assessed qualitatively. Results of this investigation revealed a relation between the three mentioned quantitative behavioral parameters and the number of hornets per group. The number of hornets per group was positively related to the extent of building, the number of cells built by a group is $$2\pi \sqrt {group size} $$ , but negatively related to the rate of building. As for the delay of building, a non-monotone relation was found. The relation between number of hornets per group and the oviposition delay was found to be non-monotone; the number of hornets per group and their longevity were found to be inversely related. Discrepanices were recorded on the very small (1–2 individuals) or very large (100 individuals) hornet groups.

Notizen: Abstract A new method for reconstructing evolutionary relationship among bacteria by use of rRNA sequence data is proposed. The method is based on the concept of fuzzy classification of probabilitiesp(i), p(i/j) andp(i/j*) (i=A,G,C,U) of each sequence. The resulting partition tree shares common features of previous works but has some new peculiarities.

Notizen: Abstract We describe the behaviour of motile microorganisms (e.g. flagellates) attracted by “gyrotaxis” to a sinking, non-motile particle (e.g. an algal cell). The model is based on the application of Stokes' solution for the flow field around the settling cell. The volume within which the flagellate is attracted to the sinking particle is determined from the trajectories of the flagellate. The model of gyrotaxis has several applications; these include the colonization of sinking marine snow particles with motile microoganisms and suspension feeding by protozoa.

Notizen: Abstract The effect of varying habitat dimensionality on the dynamics of a model predator-prey system is examined using an individual-based simulation. The general results are that in one dimension fluctuations in abundance of prey and predators occur over a large range of spatial scales (extinctions occur over many spatial scales). In two dimensions (and low mobilities of prey and predators) the dynamics become more predictably periodic at local scales and constant at larger scales due to statistical stabilization. In three dimensions, the model can become “phase-locked” with prey and predators displaying oscillations in abundance over large spatial scales.

Notizen: Abstract A rather complete model of the gluconeogenic pathway was used, with the known separate pools of mitochondrial and cytosolic oxalacetate, malate and aspartate. The fumarase, malate dehydrogenase and glutamate oxalacetate transaminase reactions were assumed to be isotopically actively reversible, but none at isotopic equilibrium. Malate was assumed to exchange actively between the mitochondrial and cytosol, while aspartate exchange was more limited, in agreement with the known electrogenic nature of aspartate export from the mitochondria. This model was fit to14C data obtained in hepatocyte studies, and to the whole rat14C data obtained by Heath and Rose (Biochem J. 227, 851–876, 1985). The latter data were easily fit to our model, when a single mitochondrial oxalacetate pool was assumed. However, invoking two mitochondrial oxalacetate pools, as proposed by Heath and Rose, with the oxalacetate formed via pyruvate carboxylase preferentially channelled to gluconeogenesis, could not be fit with the known differences in scrambling in glucose and glutamate produced from L[3-14C]lactate.