Optimal mean-variance control for discrete-time linear systems with Markovian jumps and multiplicative noises

In this paper, we consider the stochastic optimal control problem of discrete-time linear systems subject to Markov jumps and multiplicative noises under two criteria. The first one is an unconstrained mean-variance trade-off performance criterion along the time, and the second one is a minimum variance criterion along the time with constraints on the expected output. We present explicit conditions for the existence of an optimal control strategy for the problems, generalizing previous results in the literature. We conclude the paper by presenting a numerical example of a multi-period portfolio selection problem with regime switching in which it is desired to minimize the sum of the variances of the portfolio along the time under the restriction of keeping the expected value of the portfolio greater than some minimum values specified by the investor. (C) 2011 Elsevier Ltd. All rights reserved.

Dispersive Estimates of Solutions to the Wave Equation with a Potential in Dimensions Two and Three

2000 Mathematics Subject Classification: 35L15, 35B40, 47F05.

A planar reconfigurable linear rigid-body motion linkage with two operation modes

A planar reconfigurable linear (also rectilinear) rigid-body motion linkage (RLRBML) with two operation modes, that is, linear rigid-body motion mode and lockup mode, is presented using only R (revolute) joints. The RLRBML does not require disassembly and external intervention to implement multi-task requirements. It is created via combining a Robert’s linkage and a double parallelogram linkage (with equal lengths of rocker links) arranged in parallel, which can convert a limited circular motion to a linear rigid-body motion without any reference guide way. This linear rigid-body motion is achieved since the double parallelogram linkage can guarantee the translation of the motion stage, and Robert’s linkage ensures the approximate straight line motion of its pivot joint connecting to the double parallelogram linkage. This novel RLRBML is under the linear rigid-body motion mode if the four rocker links in the double parallelogram linkage are not parallel. The motion stage is in the lockup mode if all of the four rocker links in the double parallelogram linkage are kept parallel in a tilted position (but the inner/outer two rocker links are still parallel). In the lockup mode, the motion stage of the RLRBML is prohibited from moving even under power off, but the double parallelogram linkage is still moveable for its own rotation application. It is noted that further RLRBMLs can be obtained from the above RLRBML by replacing Robert’s linkage with any other straight line motion linkage (such as Watt’s linkage). Additionally, a compact RLRBML and two single-mode linear rigid-body motion linkages are presented.

Existence Results for Abstract Partial Neutral Integro-differential Equation with Unbounded Delay

In this paper we study the existence and regularity of mild solutions for a class of abstract partial neutral integro-differential equations with unbounded delay.

Hyperchaos in a Chua`s circuit with two new added branches

In this work a fourth-order Chua`s circuit, capable of generating hyperchaotic oscillations in a wide range of parameters, is presented. The circuit is obtained by adding two new branches to the original topology of the Chua`s double scroll circuit. One of the added branches is a linear inductor-resistor series connection, and the other one is a nonlinear voltage-controlled current source. A theoretical analysis of the circuit equations is presented, along with numerical and experimental results.

On the Filtering Problem for Continuous-Time Markov Jump Linear Systems with no Observation of the Markov Chain

We consider in this paper the optimal stationary dynamic linear filtering problem for continuous-time linear systems subject to Markovian jumps in the parameters (LSMJP) and additive noise (Wiener process). It is assumed that only an output of the system is available and therefore the values of the jump parameter are not accessible. It is a well known fact that in this setting the optimal nonlinear filter is infinite dimensional, which makes the linear filtering a natural numerically, treatable choice. The goal is to design a dynamic linear filter such that the closed loop system is mean square stable and minimizes the stationary expected value of the mean square estimation error. It is shown that an explicit analytical solution to this optimal filtering problem is obtained from the stationary solution associated to a certain Riccati equation. It is also shown that the problem can be formulated using a linear matrix inequalities (LMI) approach, which can be extended to consider convex polytopic uncertainties on the parameters of the possible modes of operation of the system and on the transition rate matrix of the Markov process. As far as the authors are aware of this is the first time that this stationary filtering problem (exact and robust versions) for LSMJP with no knowledge of the Markov jump parameters is considered in the literature. Finally, we illustrate the results with an example.

Linear minimum mean square filter for discrete-time linear systems with Markov jumps and multiplicative noises

In this paper we obtain the linear minimum mean square estimator (LMMSE) for discrete-time linear systems subject to state and measurement multiplicative noises and Markov jumps on the parameters. It is assumed that the Markov chain is not available. By using geometric arguments we obtain a Kalman type filter conveniently implementable in a recurrence form. The stationary case is also studied and a proof for the convergence of the error covariance matrix of the LMMSE to a stationary value under the assumption of mean square stability of the system and ergodicity of the associated Markov chain is obtained. It is shown that there exists a unique positive semi-definite solution for the stationary Riccati-like filter equation and, moreover, this solution is the limit of the error covariance matrix of the LMMSE. The advantage of this scheme is that it is very easy to implement and all calculations can be performed offline. (c) 2011 Elsevier Ltd. All rights reserved.

Value of real time three-dimensional echocardiography in patients with hypertrophic cardiomyopathy: Comparison with two-dimensional echocardiography and magnetic resonance imaging

Real time three-dimensional echocardiography (RT3DE) has been demonstrated to be an accurate technique to quantify left ventricular (LV) volumes and function in different patient populations. We sought to determine the value of RT3DE for evaluating patients with hypertrophic cardiomyopathy (HCM), in comparison with cardiac magnetic resonance imaging (MRI). Methods: We studied 20 consecutive patients with HCM who underwent two-dimensional echocardiography (2DE), RT3DE, and MRI. Parameters analyzed by echocardiography and MRI included: wall thickness, LV volumes, ejection fraction (LVEF), mass, geometric index, and dyssynchrony index. Statistical analysis was performed by Lin agreement coefficient, Pearson linear correlation and Bland-Altman model. Results: There was excellent agreement between 2DE and RT3DE (Rc = 0.92), 2DE and MRI (Rc = 0.85), and RT3DE and MRI (Rc = 0.90) for linear measurements. Agreement indexes for LV end-diastolic and end-systolic volumes were Rc = 0.91 and Rc = 0.91 between 2DE and RT3DE, Rc = 0.94 and Rc = 0.95 between RT3DE and MRI, and Rc = 0.89 and Rc = 0.88 between 2DE and MRI, respectively. Satisfactory agreement was observed between 2DE and RT3DE (Rc = 0.75), RT3DE and MRI (Rc = 0.83), and 2DE and MRI (Rc = 0.73) for determining LVEF, with a mild underestimation of LVEF by 2DE, and smaller variability between RT3DE and MRI. Regarding LV mass, excellent agreement was observed between RT3DE and MRI (Rc = 0.96), with bias of -6.3 g (limits of concordance = 42.22 to -54.73 g). Conclusion: In patients with HCM, RT3DE demonstrated superior performance than 2DE for the evaluation of myocardial hypertrophy, LV volumes, LVEF, and LV mass.

Acoustic analysis of speech through a hearing aid: perceptual effects of changes with two-channel compression

Compression amplification significantly alters the acoustic speech signal in comparison to linear amplification. The central hypothesis of the present study was that the compression settings of a two-channel aid that best preserved the acoustic properties of speech compared to linear amplification would yield the best perceptual results, and that the compression settings that most altered the acoustic properties of speech compared to linear would yield significantly poorer speech perception. On the basis of initial acoustic analysis of the test stimuli recorded through a hearing aid, two different compression amplification settings were chosen for the perceptual study. Participants were 74 adults with mild to moderate sensorineural hearing impairment. Overall, the speech perception results supported the hypothesis. A further aim of the study was to determine if variation in participants' speech perception with compression amplification (compared to linear amplification) could be explained by the individual characteristics of age, degree of loss, dynamic range, temporal resolution, and frequency selectivity; however, no significant relationships were found.

Comparison of different numerical methods for the solution of the time-fractional reaction-diffusion equation with variable diffusion coefficient

In this work we perform a comparison of two different numerical schemes for the solution of the time-fractional diffusion equation with variable diffusion coefficient and a nonlinear source term. The two methods are the implicit numerical scheme presented in [M.L. Morgado, M. Rebelo, Numerical approximation of distributed order reaction- diffusion equations, Journal of Computational and Applied Mathematics 275 (2015) 216-227] that is adapted to our type of equation, and a colocation method where Chebyshev polynomials are used to reduce the fractional differential equation to a system of ordinary differential equations

A one-dimensional Keller-Segel equation with a drift issued from the boundary

We investigate in this note the dynamics of a one-dimensional Keller-Segel type model on the half-line. On the contrary to the classical configuration, the chemical production term is located on the boundary. We prove, under suitable assumptions, the following dichotomy which is reminiscent of the two-dimensional Keller-Segel system. Solutions are global if the mass is below the critical mass, they blow-up in finite time above the critical mass, and they converge to some equilibrium at the critical mass. Entropy techniques are presented which aim at providing quantitative convergence results for the subcritical case. This note is completed with a brief introduction to a more realistic model (still one-dimensional).

Gaussian estimates for the density of the non-linear stochastic heat equation in any space dimension

In this paper, we establish lower and upper Gaussian bounds for the probability density of the mild solution to the stochastic heat equation with multiplicative noise and in any space dimension. The driving perturbation is a Gaussian noise which is white in time with some spatially homogeneous covariance. These estimates are obtained using tools of the Malliavin calculus. The most challenging part is the lower bound, which is obtained by adapting a general method developed by Kohatsu-Higa to the underlying spatially homogeneous Gaussian setting. Both lower and upper estimates have the same form: a Gaussian density with a variance which is equal to that of the mild solution of the corresponding linear equation with additive noise.

Assembling genes from predicted exons in linear time with dynamic programming

In a number of programs for gene structure prediction in higher eukaryotic genomic sequences, exon prediction is decoupled from gene assembly: a large pool of candidate exons is predicted and scored from features located in the query DNA sequence, and candidate genes are assembled from such a pool as sequences of nonoverlapping frame-compatible exons. Genes are scored as a function of the scores of the assembled exons, and the highest scoring candidate gene is assumed to be the most likely gene encoded by the query DNA sequence. Considering additive gene scoring functions, currently available algorithms to determine such a highest scoring candidate gene run in time proportional to the square of the number of predicted exons. Here, we present an algorithm whose running time grows only linearly with the size of the set of predicted exons. Polynomial algorithms rely on the fact that, while scanning the set of predicted exons, the highest scoring gene ending in a given exon can be obtained by appending the exon to the highest scoring among the highest scoring genes ending at each compatible preceding exon. The algorithm here relies on the simple fact that such highest scoring gene can be stored and updated. This requires scanning the set of predicted exons simultaneously by increasing acceptor and donor position. On the other hand, the algorithm described here does not assume an underlying gene structure model. Indeed, the definition of valid gene structures is externally defined in the so-called Gene Model. The Gene Model specifies simply which gene features are allowed immediately upstream which other gene features in valid gene structures. This allows for great flexibility in formulating the gene identification problem. In particular it allows for multiple-gene two-strand predictions and for considering gene features other than coding exons (such as promoter elements) in valid gene structures.

Investigation of the composition of linear alkylbenzenes with emphasis on the identification and quantitation of some trace compounds using GS/MS system in both electron impact and chemical ionization modes

Linear alkylbenzenes, LAB, formed by the Alel3 or HF catalyzed alkylation of benzene are common raw materials for surfactant manufacture. Normally they are sulphonated using S03 or oleum to give the corresponding linear alkylbenzene sulphonates In >95 % yield. As concern has grown about the environmental impact of surfactants,' questions have been raised about the trace levels of unreacted raw materials, linear alkylbenzenes and minor impurities present in them. With the advent of modem analytical instruments and techniques, namely GCIMS, the opportunity has arisen to identify the exact nature of these impurities and to determine the actual levels of them present in the commercial linear ,alkylbenzenes. The object of the proposed study was to separate, identify and quantify major and minor components (1-10%) in commercial linear alkylbenzenes. The focus of this study was on the structure elucidation and determination of impurities and on the qualitative determination of them in all analyzed linear alkylbenzene samples. A gas chromatography/mass spectrometry, (GCIMS) study was performed o~ five samples from the same manufacturer (different production dates) and then it was followed by the analyses of ten commercial linear alkylbenzenes from four different suppliers. All the major components, namely linear alkylbenzene isomers, followed the same elution pattern with the 2-phenyl isomer eluting last. The individual isomers were identified by interpretation of their electron impact and chemical ionization mass spectra. The percent isomer distribution was found to be different from sample to sample. Average molecular weights were calculated using two methods, GC and GCIMS, and compared with the results reported on the Certificate of Analyses (C.O.A.) provided by the manufacturers of commercial linear alkylbenzenes. The GC results in most cases agreed with the reported values, whereas GC/MS results were significantly lower, between 0.41 and 3.29 amu. The minor components, impurities such as branched alkylbenzenes and dialkyltetralins eluted according to their molecular weights. Their fragmentation patterns were studied using electron impact ionization mode and their molecular weight ions confirmed by a 'soft ionization technique', chemical ionization. The level of impurities present i~ the analyzed commercial linear alkylbenzenes was expressed as the percent of the total sample weight, as well as, in mg/g. The percent of impurities was observed to vary between 4.5 % and 16.8 % with the highest being in sample "I". Quantitation (mg/g) of impurities such as branched alkylbenzenes and dialkyltetralins was done using cis/trans-l,4,6,7-tetramethyltetralin as an internal standard. Samples were analyzed using .GC/MS system operating under full scan and single ion monitoring data acquisition modes. The latter data acquisition mode, which offers higher sensitivity, was used to analyze all samples under investigation for presence of linear dialkyltetralins. Dialkyltetralins were reported quantitatively, whereas branched alkylbenzenes were reported semi-qualitatively. The GC/MS method that was developed during the course of this study allowed identification of some other trace impurities present in commercial LABs. Compounds such as non-linear dialkyltetralins, dialkylindanes, diphenylalkanes and alkylnaphthalenes were identified but their detailed structure elucidation and the quantitation was beyond the scope of this study. However, further investigation of these compounds will be the subject of a future study.

Oscillation for a second-order neutral differential equation with impulses

We consider a certain type of second-order neutral delay differential systems and we establish two results concerning the oscillation of solutions after the system undergoes controlled abrupt perturbations (called impulses). As a matter of fact, some particular non-impulsive cases of the system are oscillatory already. Thus, we are interested in finding adequate impulse controls under which our system remains oscillatory. (C) 2009 Elsevier Inc. All rights reserved.