ALBERT

Notes: Abstract By means of the Laplace transform, the behavior of a simplified model of the cardiovascular system is mathematically formulated. This formulation allows mathematical expression of the periodicity of the cardiac output and the systemic response. With the cardiac output represented as half of a sine function cycle, the systolic aortic pressure becomes the sum of a sine term and exponential terms, while the sum of the exponential terms alone represents the diastolic pressure. The characteristics of the mathematical expressions for systole and diastole are analyzed, and some relationships of potentially practical value are derived. Variation in the parameters of the system yields mathematical results consistent with the expected physical ones.

Notes: Abstract A generalization of Landahl's approximation method (H. D. Landahl,Bull. Math. Biophysics,15, 49–61, 1953) for non-linear diffusion problems is suggested. The method is applied to sorption, desorption, and free diffusion problems involving concentration-dependent diffusion coefficients. With some limitations, the results compare favorably with those obtained by numerical methods.

Notes: Abstract Recently a theorem for representing current generators in a volume conductor by the superposition of a central dipole, quadrupole, octopole, etc., has been established by G. C. K. Yeh, J. Martinek, and H. de Beaumont (Bull. Math. Biophysics,20, 203–16, 1958). This theorem makes possible the representation of any discrete or line, surface- or volume-distributed current source by a unique model which can be determined for each given case by surface potential measurements and closed form analysis. In this paper the multipole representations of an eccentric dipole and an eccentric double-layer are obtained in terms of the various parameters of the assumed singularities, and the contributions to surface potentials due to each of the multipoles are compared. Certain numerical results corresponding to those of E. Frank (Amer. Heart J.,46, 364–78, 1953) are carried out and compared. Furthermore, the multipole representation of a partially damaged double-layer is also determined and compared with that of an undamaged one. It is concluded that within the range of parameters corresponding to human subjects the higher-order multipoles can contribute significantly to the surface potentials compared with the dipole.

Notes: Abstract In line with a recent suggestion by the author (Bull. Math. Biophysics,20, 267–73, September, 1958) that not only does the organism as a whole map on the primordial, but that each organ can also be thus mapped, it is shown that the previously introduced abstract spaces, which represent an organism, contain subspaces which map continuously on the space of the primordial. Several theorems about those subspaces are proven.

Notes: Abstract The DNA-protein coding problem is given a general algebraic formulation, the consequences of which are then explored by standard mathematical methods. To keep the treatment self-contained, the mathematical techniques to be used are explained in detail. It is demonstrated that there exista priori a countably infinite number of different abstract DNA-protein codes, thereby showing that inductive attempts to construct such a code will most likely be fruitless. A notion of ergodicity is then introduced, which imposes a number of restrictions on the admissible codes, and, in fact, these considerations enable us toderive a small portion of a code which, if our hypothesis of ergodicity is correct, must occur in nature. Finally, we discuss briefly the problem as to whether there can exist more than one DNA-protein code in nature.

Notes: Abstract The present-day practices of electrocardiography and vectorardiography are based upon the theory that the surface potential differences can be assumed to be due to a single dipole inside the body. It is shown in this paper that a dipole cannot account for all the surface potentials due to realistic current generators, and hence the determination of the current generator from surface potential measurements based upon such a theory will lead to inconsistent representations of the heart for one and the same subject. To demonstrate this point two eccentric dipoles of different strengths and locations representing two muscle fibers are taken to be the current generator in a homogeneous spherical conductor. The exact surface potentials are then expressed by means of the “interior sphere theorem” of the authors. With these expressions the magnitude, direction, and location of the resultant dipole are determined by the method of D. Gabor and C. V. Nelson (J. App. Physics,25, 413–16, 1954). The surface potentials due to this resultant dipole are again exactly expressed by means of the “interior sphere theorem” and compared with those due to the eccentric dipoles assumed. It can be seen that the differences can be considerable. It is suggested that the multipole model of the authors (Bull. Math. Biophysics,20, 203–16, 1958) be used as a more accurate and the only unique representation of the heart.

Notes: Abstract In a preceding paper (Bull. Math. Biophysics 20, 71–93, 1958) the principle of biotopological mapping was formulated in terms of a continuous mapping of an abstract space, made from the set of biological properties which characterize the organism, by an appropriate definition of neighborhoods. In this paper it is shown that we may consider directly the mappings of the different sets of properties which characterize different organisms without taking recourse to abstract spaces. All the verificable conclusions made in the preceding paper remain valid. A serious difficulty mentioned previously is, however, avoided and the possibility of more general predictions is established.

Notes: Abstract The general Theory of Categories is applied to the study of the (M, R)-systems previously defined. A set of axioms is provided which characterize “abstract (M, R)-systems”, defined in terms of the Theory of Categories. It is shown that the replication of the repair components of these systems may be accounted for in a natural way within this framework, thereby obviating the need for anad hoc postulation of a replication mechanism. A time-lag structure is introduced into these abstract (M, R)-systems. In order to apply this structure to a discussion of the “morphology” of these systems, it is necessary to make certain assumptions which relate the morphology to the time lags. By so doing, a system of abstract biology is in effect constructed. In particular, a formulation of a general Principle of Optimal Design is proposed for these systems. It is shown under what conditions the repair mechanism of the system will be localized into a spherical region, suggestive of the nuclear arrangements in cells. The possibility of placing an abstract (M, R)-system into optimal form in more than one way is then investigated, and a necessary and sufficient condition for this occurrence is obtained. Some further implications of the above assumptions are then discussed.

Notes: Abstract The transient stage of the random dispersal of logistic populations is investigated, using a Sturm-Liouville series leading to an infinite system of non-linear integral equations. These equations are then solved via a successive approximation scheme. R. A. Fisher's (steady-state) velocity of advance paradox is discussed. An illustrative example is worked to the second order of approximation.

Notes: Abstract An approximation method using a sine function is used to solve the second degree growth equation for the case in which an organism may simultaneously become dispersed throughout a uniform region. The resulting expression for a special case is compared with the expression obtained by R. Barakat (1959,Bull. Math. Biophysics,21, 141–51), giving the first two terms, by an iterative, procedure. The agreement is satisfactory.

Notes: Abstract Diffusion through a flat pore into a large open region is proportional to the linear dimension of the pore and not to its area. This was first explained by Brown and Escombe (1900) for a circular pore and is here generalized, by means of a dimensional argument, to include any type of regular opening. The problem is further generalized to include diffusion through pores of finite thickness, finite distance apart, and into finite regions. Since this problem cannot be solved exactly, an approximation method is introduced. Reasons for the credibility of the approximation are presented. It is then shown, by means of the approximation method, that the diffusive flow through a pore is equal to the total concentration difference divided by the resistance of the system. The resistance, in turn, is the sum of the resistances of all portions of the system, each of which is calculated. The result is compared with results which have been calculated exactly for limiting cases and found to agree very well. The results are then applied to a standard method of computing pore size in membranes, and it is shown that the correction factor is negligible.

Notes: Abstract V. S. Ivlev [Experimental Ecology of Nutrition of Fishes, 1955, Moscow (in Russia)] has shown that the food uptake by fishes during a fixed interval of time is an exponential function of the concentration of food. Ivlev's equation is derived here, and it is shown that it can hold only for non-stationary conditions, such as prevailed in Ivlev's experiments. For a stationary state, the rate of food uptake should tend asymptotically to a limiting value as the concentration increases, but the variation is not exponential. Different other aspects of the problem are investigated, and definite new experimental procedures suggested. The implications of Ivlev's findings on the effect of non-uniformity of food distribution upon the rate of food consumption are studied from a mathematical point of view. The conclusion is reached that whereas a fish does not, in the process of eating, move directly to an individual food particle which it perceives, it does move more or less directly to large aggregates of particles, if the latter are distributed nonuniformly.

Notes: Abstract Some relational aspects of the property of self-reproduction of biological systems are studied. If in addition to the requirement of the property of self-reproduction we add also the requirement of adaptability of the organism to changing environment, this imposes certain conditions on the topology of the graphs which represent such systems. A further study of the relational properties of such systems seems to offer the possibility of deriving the principle of biological mapping from the requirement of self-reproduction and adaptability. An examination of the problem of the original formation of such self-reproducing systems in connection with the established fact of impossibility of spontaneous generation leads to the conclusion that an organism must inhibit such processes which, in the absence of organisms, would lead to spontaneous generation.

Notes: Abstract In the human, the antagonistic, extensor-flexor system of the leg is an example of a common type of neurophysiological feedback system. After a brief introduction to the neuroanatomy and physiology of this feedback system, the paper formulates transfer functions from temporal response data available in the literature. A feedback stability analysis, based on the extension of Nyquist's stability criteria to multiple-loop systems and utilizing flow-graph techniques, demonstrates the stable behavior of the system. Expressions are given relating the sensitivity of the system to variations in muscle response and Golgi tendon organ (tension receptor) response. By considering the events taking place at synapses and end-plates during “isometric tension-small knee angle excursion” conditions as stationary stochastic processes, an external “noise” input to the system is given, whose spectrum is derived from the statistics of a shot-process representation of these events. The paper concludes with some correlations between the analytical results and clinical syndromes.

Notes: Abstract In this continuation of a previous report it is shown how the Volterra population dynamics, which underlies the statistical theory, can be based on a variational principle; how the dynamics can be generalized as regards both the behavior of total populations and migration phenomena; and how many directly observable data, such as amplitudes and frequencies of oscillation of a population, fit into the statistical theory and can test it. Such a test is carried out in some detail using the fox-catch data of Elton, with a clear indication that the theory is capable of comprehending the major statistical properties of population-time curves. A final section sketches an extension of the theory to cover secular variations of external conditions such as temperature of the environment.

Notes: Abstract The parameter domain for which the quasi-steady state assumption is valid can be considerably extended merely by a simple change of variable. This is demonstrated for a variety of biologically significant examples taken from enzyme kinetics, immunology and ecology.

Notes: Abstract The ability of random fluctuations in selection to maintain genetic diversity is greatly increased when generations overlap. This result has been derived previously using genetic models with very special assumptions about the population age structure. Here we explore its robustness in more realistic population models, with very general age structure or physiological structure. For a range of genetic models (haploid, diploid, single and multilocus) we find that the condition for maintaining genetic diversity generalizes almost without change. Genetic diversity is maintained by selection if a product of the form (generation overlap)×(selection intensity)×(variability in the selection regime) is sufficiently large, where the generation overlap is measured in units of Fisher's reproductive value. This conclusion is based on a local evolutionary stability analysis, which differs from the standard “protected polymorphism” criterion for the maintenance of genetic diversity. Simulation results match the predictions from the local stability analysis, but not those from the protected polymorphism criterion. The condition obtained here for maintaining genetic diversity requires fitness fluctuations that are substantial but well within the range observed in many studies of natural populations.

Notes: Abstract Premised on relatively simple assumptions, mathematical models like those of Monod, Pirt or Droop inadequately explain the complex transient behavior of microbial populations. In particular, these models fail to explain many aspects of the dynamics of aTetrahymena pyriformis-Escherichia coli community. In this study an alternative approach, an individual-based model, is employed to investigate the growth and interactions ofTetrahymena pyriformis andE. coli in a batch culture. Due to improved representation of physiological processes, the model provides a better agreement with experimental data of bacterial density and ciliate biomass than previous modeling studies. It predicts a much larger coexistence domain than rudimentary models, dependence of biomass dynamics on initial conditions (bacteria to ciliate biomasses ratio) and appropriate timing of minimal bacteria density. Moreover, it is found that accumulation ofE. coli sized particles andE. coli toxic metabolites has a stabilizing effect on the system.

Notes: Abstract At the core of contemporarymorphometrics—the quantitative study of biological shape variation—is a synthesis of two originally divergent methodological styles. One contributory tradition is the multivariate analysis of covariance matrices originally developed as biometrics and now dominant across a broad expanse of applied statistics. This approach, couched solely in the linear geometry of covariance structures, ignores biomathematical aspects of the original measurements. The other tributary emphasizes the direct visualization of changes in biological form. However, making objective the biological meaning of the features seen in those diagrams was always problematical; also, the representation of variation, as distinct from pairwise difference, proved infeasible. To combine these two variants of biomathematical modeling into a valid praxis for quantitative studies of biological shape was a goal earnestly sought though most of this century. That goal was finally achieved in the 1980s when techniques from mathematical statistics, multivariate biometrics, non-Euclidean geometry and computer graphics were combined in a coherent new system of tools for the complete regionalized quantitative analysis oflandmark points together with the biomedical images in which they are seen. In this morphometric synthesis, correspondence of landmarks (biologically labeled geometric points, like “bridge of the nose”) across specimens is taken as a biomathematical primitive. The shapes of configurations of landmarks are defined as equivalence classes with respect to the Euclidean similarity group and then represented as single points in David Kendall'sshape space, a Riemannian manifold with Procrustes distance as metric. All conventional multivariate strategies carry over to the study of shape variation and covariation when shapes are interpreted in the tangent space to the shape manifold at an average shape. For biomathematical interpretation of such analyses, one needs a basis for the tangent space compatible with the reality of local biotheoretical processes and explanations at many different geometric scales, and one needs graphics for visualizing average shape differences and other statistical contrasts there. Both of these needs are managed by thethin-plate spline, a deformation function that has an unusually helpful linear algebra. The spline also links the biometrics of landmarks to deformation analysis of the images from which the landmarks originally arose. This article reviews the history and principal tools of this synthesis in their biomathematical and biometrical context and demonstrates their usefulness in a study of focal neuroanatomical anomalies in schizophrenia.

Notes: Abstract A competition model describing tumor-normal cell interaction with the added effects of periodically pulsed chemotherapy is discussed. The model describes parameter conditions needed to prevent relapse following attempts to remove the tumor or tumor metastasis. The effects of resistant tumor subpopulations are also investigated and recurrence prevention strategies are explored.

Notes: Abstract Increasing attention is being paid to the configuration and development of vascular structures and their possible correlations with physiological events. The study of angiogenesis in normal and pathological states as well as in the embryo and adult has provided new insights into the mechanism of vessel growth and organization of the vasculature. Various mathematical branching models have been developed. These constructions are mainly geometrical and only involve a branching phenomenon. We propose the use of a deterministic non-linear model based on physiological laws and hydrodynamics. Growth, branching and anastomosis, the three actual main events occurring in vascular growth, are included in this model. Space growth, including cells and vessels, is defined by a decreasing transformation. Space density and the length of new sprouts are controlled by a set of parameters. The conditions on these parameters are well established, which allows the production of realistic patterns.

Notes: Abstract The quasi-stationary distribution of a population within a system of interacting populations is approximated by a stochastic logistic process. The parameters of this process can be expressed in the parameters of the full system. Using the diffusion approximation, an expression for the expected extinction time is derived from this logistic process. Since the expected extinction time is expressed in the parameters of the full system, the effect of these parameters on the extinction risk can be easily evaluated, which may be of use for studies in ecology, conservation biology and epidemiology. The outcome is compared with simulation results for the case of a prey-predator system.

Notes: Abstract The cytokines are the information superhighway of the immune system. They are an important component of the integrated behavior of the system. In order to be able to have a good understanding of the immune system, we must be able to model the effect of cytokines and their combined effect. This work is a step in that direction. We study the combined effect of two cytokines: interleukin-2 (IL-2) and interleukin-4 (IL-4) on some cells of the immune system. Interleukin-2 and interleukin-4 are important growth and differentiation factors for B and T cells. Interleukin-4 antagonizes the effect of interleukin-2 on B cells and some T cells while it synergizes with interleukin-2 on other T cells. We build a mathematical model of the interaction of both cytokines on T and B cells as a building block toward a model of the Th1/Th2 cross-regulation. The response of a given cell to the combination of interleukin-2 and interleukin-4 is shown to involve competing dynamical effects which can lead to either antagnostic or synergistic combined effect.

Notes: Abstract We propose a mathematical approach to the modelling of self-organizing hierarchies in animal societies. This approach relies on a basic positive feedback mechanism that reinforces the ability of a given individual to win or to lose in a hierarchical interaction, depending on how many times it won or lost in previous interactions. Motivated by experiments carried out on primitively eusocial waspsPolistes, the model, is based on coupled differential equations supplemented with a small stochastic term. Numerical integrations allow many different hierarchical profiles to be obtained depending on the model parameters: (1) the particular form of the probability for an individual to win or lose a fight given its history, (2) the probability of interaction between two individuals, (3) the forgetting strength, which determines the rate at which events in the past are forgotten and no longer influence the force of an individual and (4) two individual recognition parameters, which set the contribution of individual recognition in the process of hierarchical genesis. We compare the results, expressed in terms of a hierarchical index or of the Landau number that describes the degree of linearity of the hierarchy, with various experimental results.

Notes: Abstract The normal process of dermal wound healing fails in some cases, due to fibro-proliferative disorders such as keloid and hypertrophic scars. These types of abnormal healing may be regarded as pathologically excessive responses to wounding in terms of fibroblastic cell profiles and their inflammatory growth-factor mediators. Biologically, these conditions are poorly understood and current medical treatments are thus unreliable. In this paper, the authors apply an existing deterministic mathematical model for fibroplasia and wound contraction in adult mammalian dermis (Olsenet al., J. theor. Biol. 177, 113–128, 1995) to investigate key clinical problems concerning these healing disorders. A caricature model is proposed which retains the fundamental cellular and chemical components of the full model, in order to analyse the spatiotemporal dynamics of the initiation, progression, cessation and regression of fibro-contractive diseases in relation to normal healing. This model accounts for fibroblastic cell migration, proliferation and death and growth-factor diffusion, production by cells and tissue removal/decay. Explicit results are obtained in terms of the model processes and parameters. The rate of cellular production of the chemical is shown to be critical to the development of a stable pathological state. Further, cessation and/or regression of the disease depend on appropriate spatiotemporally varying forms for this production rate, which can be understood in terms of the bistability of the normal dermal and pathological steady states—a central property of the model, which is evident from stability and bifurcation analyses. The work predicts novel, biologically realistic and testable pathogenic and control mechanisms, the understanding of which will lead toward more effective strategies for clinical therapy of fibro-proliferative disorders.

Notes: Abstract Populations often exhibit abrupt changes in abundance associated with a smooth, continuous change in some component of their environment, with the abruptness usually attributed to inter-specific interactions or physical extremes. This paper presents a spatially explicit single-species population model in which intra-specific interactions alone are responsible for such an abrupt change. The essential mechanism involves cooperation in both colonization (through enhanced recruitment near other individuals) and mortality (protection through a “safety-in-numbers” interaction). Large fluctuations in population density would likely be observable near the transition region.

Notes: Abstract A mathematical model of the nitrogen transformation cycle in an aquatic environment is studied. Using Pontryagin's maximum principle, a preferential utilization of ammonium to nitrate by phytoplankton is explained and verified by experimental data. A multiparameter bifurcation is given. The model was found to have four types of equilibrium sets. It is shown that a Hopf bifurcation may occur.

Notes: Abstract Parallel computation employing a domain decomposition method was used to calculate precisely without approximations the spatio-temporal distribution of Ca2+ in nerve terminals. The results showed, contrary to expectations, that for equal admitted Ca2+ currents at low (one channel open) and high (four channels open) depolarization, the average Ca2+ concentration at the release area is higher at the low depolarization. These calculations provide additional support for the Ca2+-voltage hypothesis for neurotransmitter release.

Notes: Abstract Type I hypersensitivity, which functions to protect the organism from parasites, is caused by binding of antigen to IgE antibodies pre-attached to the cell surface of tissue mast cells and their circulating counterparts, the basophils. In “allergy,” type I hypersensitivity is inappropriately induced by protein-based foreign substances (such as pollen) or protein components of insect stings, which in the normal course of events would be cleared from the organism without causing any damage. Paradoxically, a successful clinical treatment of allergy involves repeated immunization of allergic persons with low doses of the allergen—immunotherapy. Investigation of the available experimental evidence leads to the conclusion that the phenomena of immunotherapy are best addressed in terms of the interplay among the mechanism(s) of immune memory—Th1/Th2 cross-regulation—and the physical compart-mentalization of the immune system. These conclusions are illustrated with a numerical simulation.

Notes: Abstract In the presence of seasonal forcing, predator-prey models with quadratic interaction terms and weak dissipation can exhibit infinite numbers of coexisting periodic attractors corresponding to cycles of different magnitude and frequency. These motions are best understood with reference to the conservative case, for which the degree of dissipation is, by definition, zero. Here one observes the familiar mix of “regular” (neutrally stable orbits and tori) and chaotic motion typical of non-integrable Hamiltonian systems. Perturbing away from the conservative limit, the chaos becomes transitory. In addition, the invariant tori are destroyed and the neutrally stable periodic orbits becomes stable limit cycles, the basins of attraction of which are intertwined in a complicated fashion. As a result, stochastic perturbations can bounce the system from one basin to another with consequent changes in system behavior. Biologically, weak dissipation corresponds to the case in which predators are able to regulate the density of their prey well below carrying capacity.

Notes: Abstract The standard method for measuringin vitro antibiotic efficacy is based on a point observation of bacterial activity 18 hours after inoculation. The method, while simple, forgoes significant information by ignoring the dynamics of the interations between antibiotic and bacteria. This paper proposes a simple dynamic model describing these interactions. The model consists of two non-linear differential equations of the S-system type. Its parameter values are estimated, through the minimization of residual errors, from data on the effect of the carbapenem antibiotic imipenem onPseudomonas aeruginosa. The model adequately describes the dynamic behavior of the bacterial populations in the presence of the antibiotic: beginning with drug administration, then through the decline of the bacterial population and possibly ending with bacterial resurgence.

Notes: Abstract We have studied an ecological system of two species, which we denotestrong andweak, respectively, that compete for a single food resource. This system is modelled as a three component reaction-diffusion process. In the presence of a solitary pulse of increased resources, the weaker competitor can diffuse toward this surplus, gaining a competitive advantage and hence persisting in contraposition with the classical Lotka-Volterra result. An exact analytical solution has been found through a quantum mechanical analogy. A stability analysis of this solution against changes in different parameters has been carried out.

Notes: Abstract Collapsible-tube flow with self-excited oscillations has been extensively investigated. Though physiologically relevant, forced oscillation coupled with self-excited oscillation has received little attention in this context. Based on an ODE model of collapsible-tube flow, the present study applies modern dynamics methods to investigate numerically the responses of forced oscillation to a limit-cycle oscillation which has topological characteristics discovered in previous unforced experiments. A devil's staircase and period-doubling cascades are presented with forcing frequency and amplitude as control parameters. In both cases, details are provided in a bifurcation diagram. Poincaré sections, a frequency spectrum and the largest Lyapunov exponents verify the existence of chaos in some circumstances. The thin fractal structure found in the strange attractors is believed to be a result of high damping and low stiffness in such systems.

Notes: Abstract In this work, we show that a one-dimensional model of the blood flow across the lungs can reproduce the evolution of a bolus versus the time. Solving the differential equation governing the bolus concentration in the framework of this model, we determine the solution which fulfills Gaussian initial boundary conditions. An effective parameter related to the ratio of a diffusion coefficient to the square of the mean speed of the flow is defined. The determination of its numerical values following a semi-empirical approach enables us to know accurately the mean transit time and the cardiac output. The results have been compared to other methods, and were found in good agreement. Such an approach could be of interest in all studies where the knowledge of flow—including micro-circulation—is needed.

Notes: Abstract How two species interact during and after colonization influences which of them will be present in each stage of succession. In the tolerance model of ecological succession in a patchy environment, empty patches can be colonized by any species, but the ability to tolerate reduced resource levels determines which species will exclude the other. Here, we analyze a meta-population model of the possible roles of competition in colonization and succession, using non-linear Markov chains as a mathematical framework. Different kinds of competition affect the final equilibrial, abundances of the species involved in qualitatively different ways. An explicit criterion is given to determine which interactions have stronger effects on the final equilibrial levels of the weaker, species. Precise conditions are stated for the co-existence of both species. Both species are more likely to co-exist in the presence of an intermediate disturbance frequency.

Notes: Abstract This work addresses the design of LoG filters in the frequency domain within a structure formed by the cascade of quasi-Gaussian and discrete Laplacian filters. The main feature of such a structure is that it requires half the number of convolutions of the classical structure in which the LoG transfer function is expressed as the sum of two separable transfer functions of 1-D Gaussian and LoG type. Such a perspective allows one to rephrase the design of IIR and FIR filters for edge detection as a frequency domain approximation problem solvable by standard digital filter design tools. The zero-phase IIR solutions have a good performance at low orders and approximation errors practically independent of the aperture parameter. The characteristics of the nearly linear-phase IIR filters solving the problem suggest the consideration of linear-phase FIR filters with zeros constrained on the unit circle. The use of such filters leads to remarkable computational savings with respect to the filters designed by impulse response sampling. The agreement between the edge values obtained by the filters designed according to the scheme proposed in this work and those obtained by standard techniques is very good.

Notes: Abstract We study boundary value problems for the time-harmonic form of the Maxwell equations, as well as for other related systems of equations, on arbitrary Lipschitz domains in the three-dimensional Euclidean space. The main goal is to develop the corresponding theory for Lp-integrable bounday data for optimal values of p's. We also discuss a number of relevant applications in electromagnetic scattering.

Notes: Abstract Beurling's algebra $A^*=\{f:\sum_{k=0}^{\infty} \sup_{k\le |m|} |\hat f (m)| 〈 \infty \}$ is considered. A* arises quite naturally in problems of summability of the Fourier series at Lebesgue points, whereas Wiener's algebra A of functions with absolutely convergent Fourier series arises when studying the norm convergence of linear means. Certainly, both algebras are used in some other areas. A* has many properties similar to those of A, but there are certain essential distinctions. A* is a regular Banach algebra, its space of maximal ideals coincides with $[-\pi,\pi],$ and its dual space is indicated. Analogs of Herz's and Wiener-Ditkin's theorems hold. Quantitative parameters in an analog of the Beurling-Pollard theorem differ from those for A. Several inclusion results comparing the algebra A* with certain Banach spaces of smooth functions are given. Some special properties of the analogous space for Fourier transforms on the real axis are presented. The paper ends with a summary of some open problems.

Notes: Abstract We construct an algebra of left-invariant pseudodifferential operators on SU(2). We require only that the symbols be homogeneous and C2. For Fourier-bandlimited symbols, we derive the expected formulae for composition and commutators and construct an orthonormal basis of common approximate eigenvectors that could be used to study spectral theory. Some remarks on applications to matrices of operators are made.

Notes: Abstract A novel and complete method is presented for constructing the generalized Bezout identity in polynomial form of the transfer matrix from its state-space representation of a reachable and observable system. As a result, one can easily apply well-developed synthesis and analysis theories to systems described in the frequency domain.

Notes: Abstract The aim of this paper is to obtain the uncertain value set in the complex plane for systems with real and complex parameters that are known to lie inside a ball in a weightedl p-norm. It generalizes previously available results and may be used to test the robust stability of polynomials whose coefficients lie in a weighted lp-ball.

Notes: Abstract This paper describes a new approach—a Fixed-Level Least Squares (FLLS) method for linear-phase FIR filter design. It is intended for rejection of the Gibbs phenomenon through the introduction of a set of equally spaced fixed levels in the transition band and subsequent redefinition of the approximated and weighted functions. Detailed mathematical solutions of the problem as well as many examples are given. The results in graphical form are shown as an output of the FLLS software model.

Notes: Abstract In this paper we examine order reduction of parabolic systems using modal truncation. The parabolic distributed system is first approximated using the Galerkin method. The system matrices have a special structure that allows us to find the approximate spectrum of the parabolic system. To do this we compute approximate inverses of tridiagonal, diagonally dominant symmetric matrices. The approximation leads to algorithms of orderO(n), as opposed to traditional algorithms of orderO(n), wheren is the order of the system. Finally, an example is presented to illustrate the proposed algorithm.

Notes: Abstract We give results concerning the problem of approximating the input-output maps of nonlinear discrete-time approximately finite-memory systems. Here the focus is on the linear dynamical parts of the approximating structures, and we give examples showing that these linear parts can be derived from asingle prespecified function that meets certain conditions. This is done in the context of an approximation theorem in which attention is focused on what we call “basic sets.” We also consider the related but very different problem of approximating nonlinear functionals using lattice operations or the usual linear ring operations. For this problem we give criteria, not just sufficient conditions, for approximation on compact subsets of reflexive Banach spaces (any Hilbert space is a reflexive Banach space).

Notes: Abstract New calibration methods are presented for bearing estimation with sensor location uncertainties, which are based on eigenvalue decomposition of the sample covariance matrix and three or more nondisjoint sources in known bearings. The methods can be applied to arbitrary arrays, including linear arrays, require fewer computations, and are suitable for low signal-to-noise ratio (SNR) cases. Computer simulations are given to illustrate the performance of the proposed methods.

Notes: Abstract A doubly recursive algorithm for time domain convolution with a piecewise linear weighting function is presented that combines the speed of a recursive (IIR) digital filter with the flexibility and ease of design of a nonrecursive (FIR) digital filter. The approach approximates the desired FIR weighting function by a sum-of-triangles weighting function. ForL triangles (or triangle pairs for a linear phase filter) the algorithm is of orderLN. The approximation improves with the number of triangles. A significant advantage of the algorithm compared to FFT filtering or direct convolution is that there is no necessity of a tradeoff between frequency response accuracy and computation time per output point as the data spacing decreases in the filtered signal. The computational complexity is dependent on the number of triangles chosen, not the width of the weighting function, so the algorithm is especially effective for filters with an inherently wide FIR weighting function.

Notes: Abstract A solution of the problem of correlation analysis of stochastic amplitude-angle modulated signals is presented. The solution is obtained on the basis of the theory of linear stochastic processes for the general case of nonstationary correlated Gaussian noises in both modulation channels. These noises are regarded as the results of a linear filtration of Gaussian white noise processes. As an example of the analysis a case of band-pass modulation noises is considered.

Notes: Abstract Two-dimensional (2-D) periodically shift-variant (PSV) digital filters are considered. These filters have applications in processing video signals with cyclostationary noise, scrambling digital images, and 2-D multirate signal processing. The filters are formulated in the form of a Givone-Roesser (GR) state-space model with periodic coefficients. This PSV model is then presented in block form as a shift-invariant system that also has the same GR state-space form. This block form has reduced computations and ease of analysis. An algorithm is developed that transforms any given 2-D PSV GR system to its equivalent shift-invariant model. Invertibility of this model is an important consideration, especially in image scrambling applications. A condition is established for the invertibility of the shift-invariant model of the 2-D PSV system. Also, the inverse system can be easily computed from the original. The established results are illustrated with an example.

Notes: Abstract In a recent paper a method is described for constructing certain approximations to a general element in the closure of the convex hull of a subset of an inner product space. This is of interest in connection with neural networks. Here we give an algorithm that generates simpler approximants with somewhat less computational cost.

Notes: Abstract In this paper, we gave an exponential estimate for the solutions of singular systems of differential equations with a delay and establish a comparison inequality for the systems. Based on this we estabish the criteria for asymptotic stability for the composite singular systems.

Notes: Abstract This paper presents the VLSI architectures for three-level correlator design based on 1-μm CMOS technology. The architecture performs very high speed, real-time, three-level cross-correlation of signals. Two architectures, one for serial incoming samples of signals (serial data) and the other for stored signal samples (parallel data), are described in the paper.

Notes: Abstract This paper gives the existence of and bounds for the branch voltages and currents in a given infinite network, which on each branch are related by a maximal monotone resistance function. We give a simpler proof than that of Minty, with smaller bounds, in the finite network case.

Notes: Abstract An efficient procedure is presented for computing the frequency-response templates of uncertain rational transfer functions depending on two independent interval polynomials. It consists of a modified Cohen-Sutherland algorithm for the intersection test of two rectangles and a pivoting algorithm. Because the proposed procedure traces out the boundary of a frequency-response template directly, it is, in general, more efficient than the previously proposed algorithms.

Notes: Abstract In this paper, a novel VLSI algorithm for the computation of a two-dimensional discrete cosine transform is proposed. The 2D-DCT equation can be expressed by the sum of high order cosine functions, and the algorithm can be realized by combining a highly efficient first order recursive structure with some simplified matrix multiplications, which results in highly regular hardware architecture and simple routing. The algorithm has temporal and spatial locality of connection and can be segmentized for pipeline operations, so the computation time is greatly reduced. Owing to the simplicity in hardware structure, it is especially good for VLSI implementation.

Notes: Abstract We examine a generalized matched-filter problem in which the interference is a nonstationary process generated by passing white noise through a general linear time-varying filter. First a matched filter is constructed by transforming the problem into an equivalent formulation involving stationary interference and a time-varying propagation channel. Whereas the response of a time-invariant matched filter is sampled at its peak, the response of this time-varying matched filter is normalized before sampling to account for variations in the signal power. Next a matched filter is constructed using a spectral characterization of the nonstationary interference. This construction is then used to formulate a simplified solution for the case where the rate of variation in the nonstationary interference is sufficiently small. The different solutions are illustrated by a numerical example.

Notes: Abstract The purpose of this paper is to investigate the approximation capability of elliptic basis function (EBF) neural networks. The main results are: (1) A necessary and sufficient condition for a functionS′ (R 1) to be qualified as an activation function in the hidden layer of an EBF neural network is given. (2) Every nonzero functionG ε L 2(R n ) is qualified to be an activation function in the elliptic neural network to approximate any function in L2(Rn). (3) As applications, we give new proofs of the theorems concerning the approximation capability of affine basis function (ABF) neural networks and generalized radial basis function (GRBF) neural networks (including radial basis function neural networks) with arbitrary activation functions. In particular, we solve the problem of the approximation capability of sigma-pi neural networks.

Notes: Abstract The frequency response of linear interval systems under PID controllers is analyzed. A two-stage elimination algorithm is presented to describe the boundary of the response in the complex plane. Illustrative examples are included.

Notes: Abstract This article deals with a new intelligent analog circuit design system. Unlike previous design methods, this approach introduces a formal circuit representation for both the numerical and heuristic knowledge of the design system. The predicate logic circuit representation is proposed as a new formal analog circuit description language. The language's syntax and semantics provide a precise symbolic description of analog circuit functionality at higher levels of hierarchy and network component connectivities, together with CMOS transistor sizes as the lowest level of hierarchy. It is shown how sentence conversion rules of language grammer can be used to derive transistor level circuits from input performance specifications. Language sentences have the form of clauses. In addition to clause representation, a frame representation to reflect VLSI design hierarchy is also introduced. An original bidirectional inference mechanism with elements of hypothesis has been introduced to infer designs from the knowledge in clause and frame representations. The unique feature of this circuit knowledge representation is its ability to automate the analog CMOS circuit design process. The design methodology is described in detail. The proposed iterative closed loop design system adopts an expert system approach to provide topological synthesis, uses a SPICE circuit simulator to evaluate the circuit performance, and uses a new diagnostic expert system to provide advice on how to improve the design. The implementation of the methodology and associated experimental results for simple CMOS operational amplifier designs are also presented.

Notes: Abstract A global kinetic analysis of a model consisting of an autocatalytic zymogen-activation process, in which an irreversible inhibitor competes with the zymogen for the active site of the proteinase, and a monitoring coupled reaction, in which the enzyme acts upon one of its substrates, is presented. This analysis is based on the progress curves of any of the two products released in the monitoring reaction. The general solution is applied to an important particular case in which rapid equilibrium conditions prevail. Finally, we suggest a procedure to predict whether the inhibition or activation route dominates in the steady state of the system. These results generalize our previous analysis of simpler mechanisms.

Notes: Abstract In this work, treating the artery as a thick-walled cylindrical shell made of an incompressible, isotropic and elastic solid, utilizing the large deformation theory and the stress-strain relation proposed by Demiray (1976b,Trans. ASME Ser. E, J. Appl. Mech.,98, 194–197), an explicit expression for the pulse speed is obtained and the effect of lumen pressure and the axial stretch on wave speed is discussed. Numerical results indicate that the wave speed increases with lumen pressure but decreases with the axial stretch. The results of the present model are compared with our previous work (Demiray, 1988,J. Biomech. 21, 55–58) on the same subject.

Notes: Abstract For a one-locus selection model, Svirezhev introduced an integral variational principle by defining a Lagrangian which remained stationary on the trajectory followed by the population undergoing selection. It is shown here (i) that this principle can be extended to multiple loci in some simple cases and (ii) that the Lagrangian is defined by a straightforward generalization of the one-locus case, but (iii) that in two-locus or more general models there is no straightforward extension of this principle if linkage and epistasis are present. The population trajectories can be constructed as trajectories of steepest ascent in a Riemannian metric space. A general method is formulated to find the metric tensor and the surface in the metric space on which the trajectories, which characterize the variations in the gene structure of the population, lie. The local optimality principle holds good in such a space. In the special case when all possible linkage disequilibria are zero, the phase point of then-locus genetic system moves on the surface of the product space ofn higher dimensional unit spheres in a certain Riemannian metric space of gene frequencies so that the rate of change of mean fitness is maximum along the trajectory. In the two-locus case the corresponding surface is a hyper-torus.

Notes: Abstract Most recent models of the immune network are based upon a phenomenological log bell-shaped interaction function. This function depends on a single parameter, the “field”, which is the sum of all ligand concentrations weighted by their respective affinities. The typical behavior of these models is dominated by percolation, a phenomenon in which a local stimulus spreads globally throughout the network. The usual reason for employing a log bell-shaped interaction function is that B cells are activated by cross-linking of their surface immunoglobulin receptors. Here we formally derive a new phenomenological log bell-shaped function from the chemistry of receptor cross-linking by bivalent ligand. Specifying how this new function depends on the ligand concentrations requires two fields: a binding field and a cross-linking field. When we compare the activation functions for ligand-receptor pairs with different affinities, the one-field and the two-field functions differ markedly. In the case of the one-field activation function, its graph is shifted to increasingly higher concentration as the affinity decreases but keeps its width and height. In the case of the two-field activation function, the graph of a low-affinity interaction is nested within the graphs of all higher-affinity interactions. We show that this difference in the relations among activation functions for different affinities radically changes the network behavior. In models that described B cell proliferation using the one-field activation function, network behavior was dominated by low-affinity interactions. Conversely, in our new model, the high-affinity interactions are the most significant. As a consequence, percolation is no longer the only typical network behavior.

Notes: Abstract The dynamics of a microbial community consisting of a eucaryotic ciliateTetrahymena pyriformis and procaryoticEscherichia coli in a batch culture is explored by employing an individual-based approach. In this portion of the article, Part I, population models are presented. Because both models are individual-based, models of individual organisms are developed prior to construction of the population models. The individual models use an energy budget method in which growth depends on energy gain from feeding and energy sinks such as maintenance and reproduction. These models are not limited by simplifying assumptions about constant yield, constant energy sinks and Monod growth kinetics as are traditional models of microbal organisms. Population models are generated from individual models by creating distinct individual types and assigning to each type the number of real individuals they represent. A population is a compilation of individual types that vary in a phase of cell cycle and physiological parameters such as filtering rate for ciliates and maximum anabolic rate for bacteria. An advantage of the developed models is that they realistically describe the growth of the individual cells feeding on resource which varies in density and composition. Part II, the core of the project, integrates models into a dynamic microbial community and provides model analysis based upon available data.

Notes: Abstract Mathematical models are developed for the chemotherapy of AIDS. The models are systems of differential equations describing the interaction of the HIV infected immune system with AZT chemotherapy. The models produce the three types of qualitative clinical behavior: anuninfected steady state, aninfected steady state (latency) and aprogression to AIDS state. The effect of treatment is to perturb the system from progression to AIDS back to latency. Simulation of treatment schedules is provided for the consideration of treatment regimes. The following issues of chemotherapy are addressed: (i) daily frequency of treatment, (ii) early versus late initiation of treatment and (iii) intermittent treatment with intervals of no treatment. The simulations suggest the following properties of AZT chemotherapy: (i) the daily period of treatment does not affect the outcome of the treatment, (ii) treatment should not begin until after the final decline of T cells begins (not until the T cell population falls below approximately 300 mm−3) and then, it should be administered immediately and (iii) a possible strategy for treatment which may cope with side effects and/or resistance, is to treat intermittently with chemotherapy followed by interruptions in the treatment during which either a different drug or no treatment is administered. These properties are revealed in the simulations, as the model equations incorporate AZT chemotherapy as a weakly effective treatment process. We incorporate into the model the fact that AZT treatment does not eliminate HIV, but only restrains its progress. The mathematical model, although greatly simplified as a description of an extremely complex process, offers a means to pose hypotheses concerning treatment protocols, simulate alternative strategies and guide the qualitative understanding of AIDS chemotherapy.

Notes: Abstract A stochastic cellular automata model for the population dynamics of the army antEciton burchelli on Barro Colorado Island in Panama is set up. It is simulated on the computer and shown to give good agreement with biological data. It is analysed using two approximations akin to the mean field approximation in statistical mechanics, and good agreement with the simulations is obtained. Finally, the role of distance between successive statary phase bivouacs is discussed with regard to the rate of colony growth. There are two aspects of the biological system studied here that make it of general importance. First, the population is structured, since the size of each colony of army ants is crucial. Second, the spatial behaviour of the population, as in many others, is not diffusion-like, although it is random. This has implications for the kind of model that is chosen.

Notes: Abstract An analytical model for thermal damage of retinal tissue due to absorption of laser energy by finite-sized melanin granules is developed. Since melanin is the primary absorber of visible and near-IR light in the skin and in the retina, bulk heating of tissue can be determined by superposition of individual melanin granule effects. Granules are modeled as absorbing spheres surrounded by an infinite medium of water. Analytical solutions to the heat equation result in computations that are quick and accurate. Moreover, the model does not rely on symmetric beam profiles, and so arbitrary images can be studied. The important contribution of this model is to provide a more accurate biological description of submillisecond pulse exposures than previous retinal models, while achieving agreement for longer pulses. This model can also be naturally extended into the sub-microsecond domain by including vaporization as a damage mechanism. It therefore represents the beginning of a model which can be applied across the entire pulse duration domain.

Notes: Abstract Most of the elastic tubes found in the mammalian body will collapse from a distended circular cross section and when collapsed may undergo flow-induced oscillations. A mathematical model describing fluid flow in a collapsible tube is analysed using the software package AUTO-86. AUTO-86 is used for continuation and bifurcation problems in systems of non-linear ordinary differential equations. The model is a third-order lumped-parameter type and is based on the classical “Starling resistor”; it describes the unsteady flow behaviour and, in particular, the experimentally observed self-excited oscillations, in a way which is simple enough to give physical understanding, yet still firmly based on fluid mechanical principles. Some of the bifurcation types found in this model bear close resemblance to the types suggested by experimental observations of self-excited oscillations in collapsible tubes; they thus shed some light on the various topological changes which occur in practice, particularly in view of the fact that some of the points found numerically are diffcult to achieve experimentally, while the existence of others can only be inferred indirectly and uncertainly from experiment.

Notes: Abstract A solution algorithm yielding the pressure and flow-rate distributions for steady flow in an arbitrary, tree-like network is provided. Given the tree topology, the conductance of each segment and the pressure distribution at the boundary nodes, the solution is obtained from a simple recursion based on perfect Gauss elimination. An iterative solution method using this algorithm is suggested to solve for the pressure and flow-rate distributions in an arbitrary diverging-converging (arterial-venous) network consisting of two tree-like networks which are connected to each other at the capillary nodes. A number of special solutions for tree-like networks are obtained for which the general algorithm is either simplified or can be replaced by closed form solutions of the pressure and flow-rate distributions. These special solutions can also be obtained for each tree of diverging-converging networks having particular topologies and conductance distributions. Sample numerical results are provided.

Notes: Abstract Previous game-theoretic models of reciprocity have assumed that populations are large and organisms effectively sessile. This paper analyzes an iterated prisoner's dilemma among non-sessile organisms in a finite population, on the assumption that an individual's chance of remaining in one place is not influenced by a partner's behavior. This mode of interaction is suitable for analyzing potentially cooperative behaviors that are secondary to the advantage of group formation, e.g. allogrooming among social mammals. The analysis yields necessary conditions for stable reciprocity in terms of three parameters, namely, a benefit/cost ratio, the probability of further interaction and the probability of partner retention. The results suggest that, in highly mobile organisms such as fish, birds and mammals, reciprocity may be stable only if the population is small and the relative benefit and future interaction probability are both large.

Notes: Abstract We consider the kinetics of an autocatalytic reaction network in which replication and catalytic actions are separated by a translation step. We find that the behaviour of such a system is closely related to second-order replicator equations, which describe the kinetics of autocatalytic reaction networks in which the replicators act also as catalysts. In fact, the qualitative dynamics seems to be described almost entirely be the second-order reaction rates of the replication step. For two species we recover the qualitative dynamics of the replicator equations. Larger networks show some deviations, however. A hypercyclic system consisting of three interacting species can converge toward a stable limit cycle in contrast to the replicator equation case. A singular perturbation analysis shows that the replication-translation system reduces to a second-order replicator equation if translation is fast. The influence of mutations on replication-translation networks is also very similar to the behavior of selection-mutation equations.

Notes: Abstract In this work, we studied the propagation of non-linear waves in a pre-stressed thin elastic tube filled with an inviscid fluid. In the analysis, analogous to the physiological conditions of the arteries, the tube is assumed to be subject to a uniform pressureP 0 and a constant axial stretch ratio λz. In the course of blood flow it is assumed that a large dynamic displacement is superimposed on this static field. Furthermore, assuming that the displacement gradient in the axial direction is small, the non-linear equation of motion of the tube is obtained. Using the reductive perturbation technique, the propagation of weakly non-linear waves in the long-wave approximation is investigated. It is shown that the governing equations reduce to the Korteweg-deVries equation which admits a solitary wave solution. The result is discussed for some elastic materials existing in the literature.

Notes: Abstract The present paper deals with the bifurcation analysis of a simple food chain model consisting of components like detritus, nutrients, microorganisms, phytoplankton and zooplankton in an aquatic environment. The food chain model is described by a system of differential equations. If the length of the food chain (LFCH) is equal to 3 or 4, then an asymptotically stable equilibrium exists. For LFCH=5 or 6 the non-trivial equilibrium is unstable and the food-chain model has periodic orbits.

Notes: Abstract We develop a macroscopic model for delivering drug to brain tumors. The model accounts for bulk convective and diffusive transport across the blood-brain barrier and through the interstitial space. Through mathematical analysis and simulations, we assess the effects of changing parameters (within physiological bounds) on drug delivery. We find that there is an optimal treatment for convective drug delivery to the center of the tumor. We interpret this phenomenon in terms of traffic flow. The implications of our analyses on existing chemotherapeutic protocols are discussed.

Notes: Abstract The concept of shape space proposed by Perelson and Oster (1979,J. Theor. Biol. 81, 645–670) has been a useful tool for theoretical immunologists, who have invoked it to model idiotypic binding, which plays a significant role in mathematical models of immune networks. The actual construction of such a space from its definition requires specialized experimental information, which is not completely available. In this article, we discuss, with illustrative examples, how graphical representations similar to the idea of shape space can be derived by analyzing real affinity matrices, and the relative merits of such representations to approximations that might be obtained by the approach of Perelson and Oster. We also give directions for future research with a view toward applications.

Notes: Abstract A novel intact circular dsDNA supercoil is proposed as an alternative to the conventional DNA supercoil, so that the two complementary strands of ssDNA circles are separable without any covalent bond breakage. This new structure can be visualized by using two tubings: one black and one clear. Twist the black tubing a number of times and connect its two ends. Do the same for the clear tubing. Then wrap the two tubings together. This forms the separable or novel supercoil. On the other hand, the conventional supercoil can be modeled by twisting the black and clear tubings together and then connect their respective ends, so that the two tubings are not separable unless one of them is cut. Experimentally, in the absence of any enzyme, many intact plasmid dsDNA circles give two bands on agarose gel electrophoresis under a certain given condition, while the same plasmid molecules after cutting once by a restriction enzyme give only one band under the same, condition. In the case of intact pUC19 plasmids, these two bands can then be, recovered and sequenced separately, using two primers in opposite directions. Each band gives mostly one sequence which is complementary to that of the other band. The combination of the above theoretical model and experimental results strongly suggests that there is an alternative structure of DNA which does not have the usual difficulty of unwinding, rewinding and requiring numerous covalent bond breakages and ligations during semiconservative replication.

Notes: Abstract In this paper, a closed-form identification of possibly nonminimum phase multichannel moving average (MA) processes is derived by exploiting the eigenstructures of the observation cumulant matrices using the ESPRIT algorithm. The proposed approach allows the combination of statistics of different orders for better performance and offers reduced computation complexity when compared with existing iterative approaches. Simulations are also presented to demonstrate the performance of the proposed algorithm.

Notes: Abstract This paper presents an approach to estimate the directions-of-arrival (DOAs) of narrowband coherent (completely corrected) signals by an arbitrary geometry sensor array. The concept of focussing matrices is introduced to estimate the outputs of a virtual uniform linear array from the outputs of the real array in order to use spatial smoothing techniques. Unlike the work presented in [5], the proposed approach avoids multiple singular value decomposition (SVD) computations in the various sectors because it is not necessary to divide the field of view of the array into sectors. Simulation results are given to support the approach.

Notes: Abstract An optimal deconvolution filter design method is proposed in this paper for signal transmission systems with small perturbation of parameters. The perturbative parameters of the transmission channel and noise model are of probabilistic structures. A realizable filter is derived to minimize the mean square estimation error from the viewpoint of frequency domain. The calculus of variation technique and the spectral factorization method are used in the design procedure. The design method is suitable for the deconvolution of both minimum-phase and nonminimum-phase perturbative transmission systems. The minimum mean square error of the optimal deconvolution filter is also discussed. Finally, an example is given to illustrate the simulation results of the proposed optimal deconvolution filter.

Notes: Abstract In this paper we present a simple and practical algorithm for the estimation of uncertain parameters of linear systems. The uncertainty is twofold, involving random observation noise, and possible jumps in the parameter values. The jumps may occur at unknown points in time, and are of unknown magnitudes and directions. The algorithm is based on the Kalman filter, with a single-sample hypothesis test, which is used to employ a three-state decision rule (yes, no, maybe). The “maybe” choice invokes a fading memory Kalman filter. The overall algorithm contains the constant parameter filter, fading memory filter, and the set of tests and rules that enable it to switch back and forth between the two filters. Application examples are presented.

Notes: Abstract A novel approach that employs ordered binary decision diagrams (OBDDs) is contributed to factorize multi-level logic functions by requiring as few literals as possible. A logic function with PLA format is represented as an OBDD form first. A heuristic decision method of variable ordering, called theorder lookahead method, is derived for the construction of OBDDs. This method is based on the constant cofactor and the number of erasable logic terms for each input variable. The total execution time of the OBDD construction by the above ordering decisions is very fast for some MCNC benchmarks. With the above OBDDs, we introduce a simple yet effective graph manipulation, calledEXT, to obtain a minimal number of literals in the Boolean function. This greedyEXT algorithm consists mainly of two phases. The first phase, calledgraph analysis, identifies the similarities between nodes on the same level in the OBDD. The second phase, calledtree analysis, utilizes the above features to extract the common parts of the nodes. TheEXT procedure runs from the bottom level up to the top level of the OBDD. The computational complexity depends on the number of nodes in the OBDD. The results of simulations show thatEXT has a very fast CPU execution time and a competitive literal ratio with other methods for some MCNC benchmarks.EXT will produce the smallest literal number, especially for structured or symmetric circuits.

Notes: Abstract In this study we explore the use of nonlinear embedding maps to expand the dimension of the input space. The efficacy of such maps to speed training and to enhance performance is illustrated through several examples. A natural connection to nonlinear synaptic interconnects is also developed.

Notes: Abstract We construct a general model of a nonlinear input-output system consisting of finitely many blocks whose variables are in extended spaces. This model encompasses all configurations used in control which have several feedback and feedforward loops. We show that such a system is essentially the traditional MI-MO feedback system. Moreover, assuming that one block of the system is a plant affected by perturbations, we derive conditions guaranteeing linear insensitivity and/or robust stability of the whole system. In particular, we consider systems whose variables are continuous, vector-valued functions on [0, ∞), and whose blocks are described by nonlinear Volterra operators. The applications of the results are illustrated by several examples of feedback-feedforward systems.

Notes: Abstract The design of a constant resistance multiplexer composed of three reactance ladder networks connected either in series or in parallel is a very practical problem. In this paper, an improvement on Norton's method and a set of necessary and sufficient conditions on the realizability of a Butterworth or elliptic constant resistance multiplexer are presented. Two illustrative examples are given to show the design procedure.

Notes: Abstract This paper extends Leslie's vector and matrix treatment of populations. A simple matrix is given for species in which adult mortality and fertility are independent of age, but in which the juvenile mortality rate differs from the adult. The population vector can be changed into a population matrix. This should allow equations using functions of the size of the population to be developed. Genetic variables such as sex or other polymorphisms can be introduced, and the notation allows different systems of selection or non-random mating to be specified.

Notes: Abstract This paper presents a new interval Pade approximation method to convert a continuous-time (discrete-time) uncertain linear system to an equivalent discrete-time (continuous-time) uncertain model via interval arithmetic operations. Based on the inclusion theorem related to the interval arithmetic, the interval Pade's approximants and their associated interval error matrices with interval arguments are obtained via the Pade's approximants and their associated error matrices with degenerate (real) arguments, respectively. Tighter error bounds of various approximate uncertain models with respect to their exact uncertain models are determined and used to modify the obtained Pade's approximants, so that the resulting approximate uncertain models are able to tightly enclose the original uncertain systems. Thus, the analysis and design of the original uncertain systems can be indirectly carried out using the converted uncertain models in either the continuous-time or the discrete-time domain.

Notes: Abstract Three new concepts — breakdown points, breakdown probabilities, and midpoint sensitivity curves — for stack filter analysis are introduced and analyzed in this paper. Breakdown points and probabilities can be used as measures of the robustness of stack filters. Midpoint sensitivity curves in turn give information on how sensitive the output of a stack filter is to the changes of a single value in the input window. The second major contribution of this paper is the extension of the current optimality theory of stack filters. This theory combines noise attenuation and different constraints on the filter's behavior. New constraints are introduced in this paper. A new optimization approach based on breakdown probability as a noise attenuation measure is also derived. In certain special cases it is shown that the optimal stack filter that achieves the best noise attenuation subject to given constraints can be obtained in closed form. An algorithm for finding this form is given in this paper, and its modification for finding a stack filter having (approximately) a required rank selection vector is presented.

Notes: Abstract This paper considers the estimation of the parameters of a linear discrete-time system from noise-perturbed input and output measurements. The conditions imposed on the system are fairly general. In particular, the input and output noises are allowed to be auto-correlated and they may be cross-correlated as well. The proposed method makes use of an instrumental variable (IV)-vector to produce a covariance matrix that is pre- and postmultiplied by some prechosen weights. The singular vectors of the so-obtained matrix possess complete information on the system parameters. A weighted subspace fitting (WSF) method is then applied to the aforementioned singular vectors to consistently estimate the parameters of the system. The IV-WSF technique suggested herein is noniterative and easy to implement, and has a small computational burden. The asymptotic distribution of its estimation errors is derived and the result is used to motivate the choice of the weighting matrix in the WSF step and also to predict the estimation accuracy. Numerical examples are included to illustrate the performance achievable by the method.

Notes: Abstract This paper considers a (J, J′)-spectral factorization of a signed spectral matrix associated with a discrete-time improper transfer matrix by using the descriptor system representation and the generalized algebraic Riccati equation (GARE). Under the assumption that the improper transfer matrix is stable, a numerical algorithm for (J, J′)-spectral factorization is developed based on solutions of the generalized eigenvalue problem (GEP) and a related matrix equation. For an unstable transfer matrix, aJ-unitary conjugation method is applied to obtain a stable transfer matrix generating the same signed spectral matrix. Then a Hankel norm model reduction problem is briefly discussed. A simple example is also included to illustrate the numerical procedure.

Notes: Abstract By means of norm,M-matrix, and matrix measure techniques, this paper estimates several restricted regions in the complex plane in which all eigenvalues of a class of discrete time-delay systems subjected to highly structured parametric perturbations are located. Both the stability and the instability conditions for these systems are also investigated via the proposed schemes. Two numerical examples are given to verify the correctness and demonstrate the applicability of the quantitative results.

Notes: Abstract In this paper, we discuss a unified approach and its simple algorithm for estimating and tracking the location of the desired talker in the far field and the near field by microphone arrays. The approach is unified for the far field and the near field and can enhance the effective power of the talker's signal. The simple algorithm further avoids singular-value decomposition (SVD) computations and is suitable for real-time processing. Computer simulations and experiments using actual array data support the techniques.

Notes: Abstract We study some explicit functions introduced by Riemann, Jordan, Levy, Kahane... . These functions share the property of having a dense set of discontinuities. We prove that they are examples of multifractal functions.

Notes: Abstract We show that uniqueness and existence for signal reconstruction from multiscale edges in the Mallat and Zhong algorithm become possible if we restrict our signals to Paley-Wiener space, band-limit our wavelets, and irregularly sample at the wavelet transform (absolute) maxima—the edges—while possibly including (enough) extra points at each level. We do this in a setting that closely resembles the numerical analysis setting of Mallat and Zhong and that seems to capture something of the essence of their (practical) reconstruction method. Our work builds on a uniqueness result for reconstructing an L2 signal from irregular sampling of its wavelet transform of Grochenig and the related work of Benedetto, Heller, Mallat, and Zhong. We show that the rate of convergence for this reconstruction algorithm is geometric and computable in advance. Finally, we consider the effect on the rate of convergence of not sampling enough local maxima.

Notes: Abstract The main purpose of this paper is to give a procedure to "mollify" the low-pass filters of a large number of Minimally Supported Frequency (MSF) wavelets so that the smoother functions obtained in this way are also low-pass filters for an MRA. Hence, we are able to approximate (in the L2-norm) MSF wavelets by wavelets with any desired degree of smoothness on the Fourier transform side. Although the MSF wavelets we consider are bandlimited, this may not be true for their smooth approximations. This phenomena is related to the invariant cycles under the transformation $x\mapsto 2x (\mbox{mod}2\pi).$ We also give a characterization of all low-pass filters for MSF wavelets. Throughout the paper new and interesting examples of wavelets are described.

Notes: Abstract Systematic methods are developed for simplified implementations of cascaded stack and WOS filters. For recursive stack and WOS filters, corresponding nonrecursive implementations are given, with linear complexity with respect to the number of iterations. Dynamic domino logic is proposed for VLSI hardware implementation of positive Boolean functions, and a pipelined stack filter architecture is described.

Notes: Abstract In this study we consider a multilayer perceptron network with sigmoidal activation and trained via the backpropagation algorithm. The output of all neurons is collected and a simple linear regression is performed. It is shown that untrained networks with randomly chosen coefficients perform comparably with fully trained networks. This result casts a new light on the role of activation functions, the impact of dimensionality, and the efficacy of training algorithms such as backpropagation.