一类分数阶微分方程的最优控制

主要研究了一类分数阶微分方程的最优控制问题.通过Oustaloup迭代逼近,可以将分数阶微分算子在频率域范围内进行近似.那么原先的分数阶微分方程就转换为一般的常微分方程.利用这个关系,可以得到关于分数阶微分方程的两个定理和一个引理.最后给出一个例子说明该方法的有效性.     

The dynamic systems approach to control and regulation of intracellular networks

Systems theory and cell biology have enjoyed a long relationship that has received renewed interest in recent years in the context of systems biology. The term 'systems' in systems biology comes from systems theory or dynamic systems theory: systems biology is defined through the application of systems- and signal-oriented approaches for an understanding of inter- and intra-cellular dynamic processes. The aim of the present text is to review the systems and control perspective of dynamic systems. The biologist's conceptual framework for representing the variables of a biochemical reaction network, and for describing their relationships, are pathway maps. A principal goal of systems biology is to turn these static maps into dynamic models, which can provide insight into the temporal evolution of biochemical reaction networks. Towards this end, we review the case for differential equation models as a 'natural' representation of causal entailment in pathways. Block-diagrams, commonly used in the engineering sciences, are introduced and compared to pathway maps. The stimulus-response representation of a molecular system is a necessary condition for an understanding of dynamic interactions among the components that make up a pathway. Using simple examples, we show how biochemical reactions are modelled in the dynamic systems framework and visualized using block-diagrams.     

Toward a classification of isodynamic feed-forward motifs

A preceding study analysed how the topology of network motifs affects the overall rate of the underlying biochemical processes. Surprisingly, it was shown that topologically non-isomorphic motifs can still be isodynamic in the sense that they exhibit the exact same performance rate. Because of the high prevalence of feed-forward functional modules in biological networks, one may hypothesize that evolution tends to favour motifs with faster dynamics. As a step towards ranking the efficiency of feed-forward network motifs, we use a linear flow model to prove theorems establishing that certain classes of motifs are isodynamic. In partitioning the class of all motifs on n nodes into equivalence classes based upon their dynamics, we establish a basis for comparing the efficiency/performance rates of different motifs. The potential biological importance of the theorems is briefly discussed and is the subject of an ongoing large-scale project.     

Repeated applications of a transdermal patch: analytical solution and optimal control of the delivery rate

The integral transform technique was implemented to solve a mathematical model developed for percutaneous drug absorption. The model included repeated application and removal of a patch from the skin. Fick's second law of diffusion was used to study the transport of a medicinal agent through the vehicle and subsequent penetration into the stratum corneum. Eigenmodes and eigenvalues were computed and introduced into an inversion formula to estimate the delivery rate and the amount of drug in the vehicle and the skin. A dynamic programming algorithm calculated the optimal doses necessary to achieve a desired transdermal flux. The analytical method predicted profiles that were in close agreement with published numerical solutions and provided an automated strategy to perform therapeutic drug monitoring and control.     

From molecular networks to qualitative cell behavior

Adaptation and behavior are characteristics of life which are fundamentally dynamic. If we want to model the living cell we have to describe it as a dynamic system. Typical dynamic models are based on quantitative differential equations requiring very detailed kinetic knowledge. Alternative modeling techniques for less fine-grained information are better suited to available functional genomics data. As such, constraint-based techniques and qualitative modeling have proven themselves to be valid approaches in cell biology. These approaches offer formal support to check the consistency of molecular networks against phenotypic observations in the light of dynamic systems.     

Engineering virtual cardiac tissue

The kinetics of proteins involved in ion transfer, sequestration and binding in cardiac cells can be modelled to construct a model of the electrical activity of isolated cardiac cells as a system of ordinary differential equations. These cell models may be incorporated into tissue models, which, when combined with histology and anatomy, form virtual tissues. The effects of changes in specific protein expression, or changes in protein kinetics, produced by mutations or pharmacological agents, can be simulated using these tissue models and used to account for the whole organ effects of changes in specific ion-transport protein activity.     

Bistability and oscillations in chemical reaction networks

Bifurcation theory is one of the most widely used approaches for analysis of dynamical behaviour of chemical and biochemical reaction networks. Some of the interesting qualitative behaviour that are analyzed are oscillations and bistability (a situation where a system has at least two coexisting stable equilibria). Both phenomena have been identified as central features of many biological and biochemical systems. This paper, using the theory of stoichiometric network analysis (SNA) and notions from algebraic geometry, presents sufficient conditions for a reaction network to display bifurcations associated with these phenomena. The advantage of these conditions is that they impose fewer algebraic conditions on model parameters than conditions associated with standard bifurcation theorems. To derive the new conditions, a coordinate transformation will be made that will guarantee the existence of branches of positive equilibria in the system. This is particularly useful in mathematical biology, where only positive variable values are considered to be meaningful. The first part of the paper will be an extended introduction to SNA and algebraic geometry-related methods which are used in the coordinate transformation and set up of the theorems. In the second part of the paper we will focus on the derivation of bifurcation conditions using SNA and algebraic geometry. Conditions will be derived for three bifurcations: the saddle-node bifurcation, a simple branching point, both linked to bistability, and a simple Hopf bifurcation. The latter is linked to oscillatory behaviour. The conditions derived are sufficient and they extend earlier results from stoichiometric network analysis as can be found in (Aguda and Clarke in J Chem Phys 87:3461–3470, 1987; Clarke and Jiang in J Chem Phys 99:4464–4476, 1993; Gatermann et al. in J Symb Comput 40:1361–1382, 2005). In these papers some necessary conditions for two of these bifurcations were given. A set of examples will illustrate that algebraic conditions arising from given sufficient bifurcation conditions are not more difficult to interpret nor harder to calculate than those arising from necessary bifurcation conditions. Hence an increasing amount of information is gained at no extra computational cost. The theory can also be used in a second step for a systematic bifurcation analysis of larger reaction networks. We have added a dedication of the paper to K. Gatermann.     

Optimization of harvesting age in an integral age-dependent model of population dynamics

The paper deals with optimal control in a linear integral age-dependent model of population dynamics. A problem for maximizing the harvesting return on a finite time horizon is formulated and analyzed. The optimal controls are the harvesting age and the rate of population removal by harvesting. The gradient and necessary condition for an extremum are derived. A qualitative analysis of the problem is provided. The model shows the presence of a zero-investment period. A preliminary asymptotic analysis indicates possible turnpike properties of the optimal harvesting age. Biological interpretation of all results is provided.     

Minimizing the dependency ratio in a population with below-replacement fertility through immigration

Many industrialized countries face fertility rates below replacement level, combined with declining mortality especially in older ages. Consequently, the populations of these countries have started to age. One important indicator of age structures is the dependency ratio which is the ratio of the nonworking age population to the working age population. In this work we find the age-specific immigration profile that minimizes the dependency ratio in a stationary population with below-replacement fertility. It is assumed that the number of immigrants per age is limited. We consider two alternative policies. In the first one, we fix the total number of people who annually immigrate to a country. In the second one, we prescribe the size of the receiving country’s population. For both cases we provide numerical results for the optimal immigration profile, for the resulting age structure of the population, as well as for the dependency ratio.     

Stability properties of pulse vaccination strategy in SEIR epidemic model

The problem of the applicability of the pulse vaccination strategy (PVS) for the stable eradication of some relevant general class of infectious diseases is analyzed in terms of study of local asymptotic stability (LAS) and global asymptotic stability (GAS) of the periodic eradication solution for the SEIR epidemic model in which is included the PVS. Demographic variations due or not to diseased-related fatalities are also considered. Due to the non-triviality of the Floquet's matrix associate to the studied model, the LAS is studied numerically and in this way it is found a simple approximate (but analytical) sufficient criterion which is an extension of the LAS constraint for the stability of the trivial equilibrium in SEIR model without vaccination. The numerical simulations also seem to suggest that the PVS is slightly more efficient than the continuous vaccination strategy. Analytically, the GAS of the eradication solutions is studied and it is demonstrated that the above criteria for the LAS guarantee also the GAS.     

Novel recurrent neural network for modelling biological networks: Oscillatory p53 interaction dynamics

Understanding the control of cellular networks consisting of gene and protein interactions and their emergent properties is a central activity of Systems Biology research. For this, continuous, discrete, hybrid, and stochastic methods have been proposed. Currently, the most common approach to modelling accurate temporal dynamics of networks is ordinary differential equations (ODE). However, critical limitations of ODE models are difficulty in kinetic parameter estimation and numerical solution of a large number of equations, making them more suited to smaller systems. In this article, we introduce a novel recurrent artificial neural network (RNN) that addresses above limitations and produces a continuous model that easily estimates parameters from data, can handle a large number of molecular interactions and quantifies temporal dynamics and emergent systems properties. This RNN is based on a system of ODEs representing molecular interactions in a signalling network. Each neuron represents concentration change of one molecule represented by an ODE. Weights of the RNN correspond to kinetic parameters in the system and can be adjusted incrementally during network training. The method is applied to the p53-Mdm2 oscillation system – a crucial component of the DNA damage response pathways activated by a damage signal. Simulation results indicate that the proposed RNN can successfully represent the behaviour of the p53-Mdm2 oscillation system and solve the parameter estimation problem with high accuracy. Furthermore, we presented a modified form of the RNN that estimates parameters and captures systems dynamics from sparse data collected over relatively large time steps. We also investigate the robustness of the p53-Mdm2 system using the trained RNN under various levels of parameter perturbation to gain a greater understanding of the control of the p53-Mdm2 system. Its outcomes on robustness are consistent with the current biological knowledge of this system. As more quantitative data become available on individual proteins, the RNN would be able to refine parameter estimation and mapping of temporal dynamics of individual signalling molecules as well as signalling networks as a system. Moreover, RNN can be used to modularise large signalling networks.     

Phase and amplitude instability in delay-diffusion population models

We analyze three simple population models that include generation times in their growth dynamics. In the presence of a slow one-dimensional diffusion process it is shown that for sufficiently small wave numbers, the amplitude of the homogeneous limit cycle solution is unstable and the phase of the bifurcating diffusion wave obeys a Burgers' type equation with a negative coefficient in the diffusion term. Numerical solutions for the phase and amplitude in the post-critical regime display a turbulent like behavior in space and time when the size of the system is larger than some critical value. This result follows from the coupling between the delay and diffusion terms.     

Optimal Reproduction Strategies in Two Species of Mound-Building Termites

We formulate a mathematical model for food collection and production of workers and nymphs in 2 species of mound building termites. We maximise the number of nymphs (reproductives) produced by each colony over its lifetime with respect to the proportion of eggs that hatch as nymphs as opposed to workers. The results predict that food storage has a very important influence on the pattern of nymph and worker production. Food storage affects the part of the year that nymph production dominates, whether nymphs and workers are produced at the same time or not, and the existence of a final phase in the colony’s life when a very large number of nymphs but no workers are produced.