Linear-constraint wavefront control for exoplanet coronagraphic imaging systems

Science.gov (United States)

Sun, He; Eldorado Riggs, A. J.; Kasdin, N. Jeremy; Vanderbei, Robert J.; Groff, Tyler Dean

2017-01-01

A coronagraph is a leading technology for achieving high-contrast imaging of exoplanets in a space telescope. It uses a system of several masks to modify the diffraction and achieve extremely high contrast in the image plane around target stars. However, coronagraphic imaging systems are very sensitive to optical aberrations, so wavefront correction using deformable mirrors (DMs) is necessary to avoid contrast degradation in the image plane. Electric field conjugation (EFC) and Stroke minimization (SM) are two primary high-contrast wavefront controllers explored in the past decade. EFC minimizes the average contrast in the search areas while regularizing the strength of the control inputs. Stroke minimization calculates the minimum DM commands under the constraint that a target average contrast is achieved. Recently in the High Contrast Imaging Lab at Princeton University (HCIL), a new linear-constraint wavefront controller based on stroke minimization was developed and demonstrated using numerical simulation. Instead of only constraining the average contrast over the entire search area, the new controller constrains the electric field of each single pixel using linear programming, which could led to significant increases in speed of the wavefront correction and also create more uniform dark holes. As a follow-up of this work, another linear-constraint controller modified from EFC is demonstrated theoretically and numerically and the lab verification of the linear-constraint controllers is reported. Based on the simulation and lab results, the pros and cons of linear-constraint controllers are carefully compared with EFC and stroke minimization.

Linear versus non-linear supersymmetry, in general

Energy Technology Data Exchange (ETDEWEB)

Ferrara, Sergio [Theoretical Physics Department, CERN,CH-1211 Geneva 23 (Switzerland); INFN - Laboratori Nazionali di Frascati,Via Enrico Fermi 40, I-00044 Frascati (Italy); Department of Physics and Astronomy, UniversityC.L.A.,Los Angeles, CA 90095-1547 (United States); Kallosh, Renata [SITP and Department of Physics, Stanford University,Stanford, California 94305 (United States); Proeyen, Antoine Van [Institute for Theoretical Physics, Katholieke Universiteit Leuven,Celestijnenlaan 200D, B-3001 Leuven (Belgium); Wrase, Timm [Institute for Theoretical Physics, Technische Universität Wien,Wiedner Hauptstr. 8-10, A-1040 Vienna (Austria)

2016-04-12

We study superconformal and supergravity models with constrained superfields. The underlying version of such models with all unconstrained superfields and linearly realized supersymmetry is presented here, in addition to the physical multiplets there are Lagrange multiplier (LM) superfields. Once the equations of motion for the LM superfields are solved, some of the physical superfields become constrained. The linear supersymmetry of the original models becomes non-linearly realized, its exact form can be deduced from the original linear supersymmetry. Known examples of constrained superfields are shown to require the following LM’s: chiral superfields, linear superfields, general complex superfields, some of them are multiplets with a spin.

Linear versus non-linear supersymmetry, in general

International Nuclear Information System (INIS)

Ferrara, Sergio; Kallosh, Renata; Proeyen, Antoine Van; Wrase, Timm

2016-01-01

We study superconformal and supergravity models with constrained superfields. The underlying version of such models with all unconstrained superfields and linearly realized supersymmetry is presented here, in addition to the physical multiplets there are Lagrange multiplier (LM) superfields. Once the equations of motion for the LM superfields are solved, some of the physical superfields become constrained. The linear supersymmetry of the original models becomes non-linearly realized, its exact form can be deduced from the original linear supersymmetry. Known examples of constrained superfields are shown to require the following LM’s: chiral superfields, linear superfields, general complex superfields, some of them are multiplets with a spin.

Novel methods for Solving Economic Dispatch of Security-Constrained Unit Commitment Based on Linear Programming

Science.gov (United States)

Guo, Sangang

2017-09-01

There are two stages in solving security-constrained unit commitment problems (SCUC) within Lagrangian framework: one is to obtain feasible units’ states (UC), the other is power economic dispatch (ED) for each unit. The accurate solution of ED is more important for enhancing the efficiency of the solution to SCUC for the fixed feasible units’ statues. Two novel methods named after Convex Combinatorial Coefficient Method and Power Increment Method respectively based on linear programming problem for solving ED are proposed by the piecewise linear approximation to the nonlinear convex fuel cost functions. Numerical testing results show that the methods are effective and efficient.

A New Interpolation Approach for Linearly Constrained Convex Optimization

KAUST Repository

Espinoza, Francisco

2012-08-01

In this thesis we propose a new class of Linearly Constrained Convex Optimization methods based on the use of a generalization of Shepard\\'s interpolation formula. We prove the properties of the surface such as the interpolation property at the boundary of the feasible region and the convergence of the gradient to the null space of the constraints at the boundary. We explore several descent techniques such as steepest descent, two quasi-Newton methods and the Newton\\'s method. Moreover, we implement in the Matlab language several versions of the method, particularly for the case of Quadratic Programming with bounded variables. Finally, we carry out performance tests against Matab Optimization Toolbox methods for convex optimization and implementations of the standard log-barrier and active-set methods. We conclude that the steepest descent technique seems to be the best choice so far for our method and that it is competitive with other standard methods both in performance and empirical growth order.

Q-deformed systems and constrained dynamics

International Nuclear Information System (INIS)

Shabanov, S.V.

1993-01-01

It is shown that quantum theories of the q-deformed harmonic oscillator and one-dimensional free q-particle (a free particle on the 'quantum' line) can be obtained by the canonical quantization of classical Hamiltonian systems with commutative phase-space variables and a non-trivial symplectic structure. In the framework of this approach, classical dynamics of a particle on the q-line coincides with the one of a free particle with friction. It is argued that q-deformed systems can be treated as ordinary mechanical systems with the second-class constraints. In particular, second-class constrained systems corresponding to the q-oscillator and q-particle are given. A possibility of formulating q-deformed systems via gauge theories (first-class constrained systems) is briefly discussed. (orig.)

Model Predictive Control Based on Kalman Filter for Constrained Hammerstein-Wiener Systems

Directory of Open Access Journals (Sweden)

Man Hong

2013-01-01

Full Text Available To precisely track the reactor temperature in the entire working condition, the constrained Hammerstein-Wiener model describing nonlinear chemical processes such as in the continuous stirred tank reactor (CSTR is proposed. A predictive control algorithm based on the Kalman filter for constrained Hammerstein-Wiener systems is designed. An output feedback control law regarding the linear subsystem is derived by state observation. The size of reaction heat produced and its influence on the output are evaluated by the Kalman filter. The observation and evaluation results are calculated by the multistep predictive approach. Actual control variables are computed while considering the constraints of the optimal control problem in a finite horizon through the receding horizon. The simulation example of the CSTR tester shows the effectiveness and feasibility of the proposed algorithm.

Robust model predictive control for constrained continuous-time nonlinear systems

Science.gov (United States)

Sun, Tairen; Pan, Yongping; Zhang, Jun; Yu, Haoyong

2018-02-01

In this paper, a robust model predictive control (MPC) is designed for a class of constrained continuous-time nonlinear systems with bounded additive disturbances. The robust MPC consists of a nonlinear feedback control and a continuous-time model-based dual-mode MPC. The nonlinear feedback control guarantees the actual trajectory being contained in a tube centred at the nominal trajectory. The dual-mode MPC is designed to ensure asymptotic convergence of the nominal trajectory to zero. This paper extends current results on discrete-time model-based tube MPC and linear system model-based tube MPC to continuous-time nonlinear model-based tube MPC. The feasibility and robustness of the proposed robust MPC have been demonstrated by theoretical analysis and applications to a cart-damper springer system and a one-link robot manipulator.

Mathematical Modeling of Constrained Hamiltonian Systems

NARCIS (Netherlands)

Schaft, A.J. van der; Maschke, B.M.

1995-01-01

Network modelling of unconstrained energy conserving physical systems leads to an intrinsic generalized Hamiltonian formulation of the dynamics. Constrained energy conserving physical systems are directly modelled as implicit Hamiltonian systems with regard to a generalized Dirac structure on the

An improved partial bundle method for linearly constrained minimax problems

Directory of Open Access Journals (Sweden)

Chunming Tang

2016-02-01

Full Text Available In this paper, we propose an improved partial bundle method for solving linearly constrained minimax problems. In order to reduce the number of component function evaluations, we utilize a partial cutting-planes model to substitute for the traditional one. At each iteration, only one quadratic programming subproblem needs to be solved to obtain a new trial point. An improved descent test criterion is introduced to simplify the algorithm. The method produces a sequence of feasible trial points, and ensures that the objective function is monotonically decreasing on the sequence of stability centers. Global convergence of the algorithm is established. Moreover, we utilize the subgradient aggregation strategy to control the size of the bundle and therefore overcome the difficulty of computation and storage. Finally, some preliminary numerical results show that the proposed method is effective.

Application of pattern search method to power system security constrained economic dispatch with non-smooth cost function

International Nuclear Information System (INIS)

Al-Othman, A.K.; El-Naggar, K.M.

2008-01-01

Direct search methods are evolutionary algorithms used to solve optimization problems. (DS) methods do not require any information about the gradient of the objective function at hand while searching for an optimum solution. One of such methods is Pattern Search (PS) algorithm. This paper presents a new approach based on a constrained pattern search algorithm to solve a security constrained power system economic dispatch problem (SCED) with non-smooth cost function. Operation of power systems demands a high degree of security to keep the system satisfactorily operating when subjected to disturbances, while and at the same time it is required to pay attention to the economic aspects. Pattern recognition technique is used first to assess dynamic security. Linear classifiers that determine the stability of electric power system are presented and added to other system stability and operational constraints. The problem is formulated as a constrained optimization problem in a way that insures a secure-economic system operation. Pattern search method is then applied to solve the constrained optimization formulation. In particular, the method is tested using three different test systems. Simulation results of the proposed approach are compared with those reported in literature. The outcome is very encouraging and proves that pattern search (PS) is very applicable for solving security constrained power system economic dispatch problem (SCED). In addition, valve-point effect loading and total system losses are considered to further investigate the potential of the PS technique. Based on the results, it can be concluded that the PS has demonstrated ability in handling highly nonlinear discontinuous non-smooth cost function of the SCED. (author)

Reachable Distance Space: Efficient Sampling-Based Planning for Spatially Constrained Systems

KAUST Repository

Xinyu Tang,; Thomas, S.; Coleman, P.; Amato, N. M.

2010-01-01

reachable distance space (RD-space), in which all configurations lie in the set of constraint-satisfying subspaces. This enables us to directly sample the constrained subspaces with complexity linear in the number of the robot's degrees of freedom

Chance-constrained/stochastic linear programming model for acid rain abatement. I. Complete colinearity and noncolinearity

Energy Technology Data Exchange (ETDEWEB)

Ellis, J H; McBean, E A; Farquhar, G J

1985-01-01

A Linear Programming model is presented for development of acid rain abatement strategies in eastern North America. For a system comprised of 235 large controllable point sources and 83 uncontrolled area sources, it determines the least-cost method of reducing SO/sub 2/ emissions to satisfy maximum wet sulfur deposition limits at 20 sensitive receptor locations. In this paper, the purely deterministic model is extended to a probabilistic form by incorporating the effects of meteorologic variability on the long-range pollutant transport processes. These processes are represented by source-receptor-specific transfer coefficients. Experiments for quantifying the spatial variability of transfer coefficients showed their distributions to be approximately lognormal with logarithmic standard deviations consistently about unity. Three methods of incorporating second-moment random variable uncertainty into the deterministic LP framework are described: Two-Stage Programming Under Uncertainty, Chance-Constrained Programming and Stochastic Linear Programming. A composite CCP-SLP model is developed which embodies the two-dimensional characteristics of transfer coefficient uncertainty. Two probabilistic formulations are described involving complete colinearity and complete noncolinearity for the transfer coefficient covariance-correlation structure. The completely colinear and noncolinear formulations are considered extreme bounds in a meteorologic sense and yield abatement strategies of largely didactic value. Such strategies can be characterized as having excessive costs and undesirable deposition results in the completely colinear case and absence of a clearly defined system risk level (other than expected-value) in the noncolinear formulation.

A Chance-Constrained Economic Dispatch Model in Wind-Thermal-Energy Storage System

Directory of Open Access Journals (Sweden)

Yanzhe Hu

2017-03-01

Full Text Available As a type of renewable energy, wind energy is integrated into the power system with more and more penetration levels. It is challenging for the power system operators (PSOs to cope with the uncertainty and variation of the wind power and its forecasts. A chance-constrained economic dispatch (ED model for the wind-thermal-energy storage system (WTESS is developed in this paper. An optimization model with the wind power and the energy storage system (ESS is first established with the consideration of both the economic benefits of the system and less wind curtailments. The original wind power generation is processed by the ESS to obtain the final wind power output generation (FWPG. A Gaussian mixture model (GMM distribution is adopted to characterize the probabilistic and cumulative distribution functions with an analytical expression. Then, a chance-constrained ED model integrated by the wind-energy storage system (W-ESS is developed by considering both the overestimation costs and the underestimation costs of the system and solved by the sequential linear programming method. Numerical simulation results using the wind power data in four wind farms are performed on the developed ED model with the IEEE 30-bus system. It is verified that the developed ED model is effective to integrate the uncertain and variable wind power. The GMM distribution could accurately fit the actual distribution of the final wind power output, and the ESS could help effectively decrease the operation costs.

Reachable Distance Space: Efficient Sampling-Based Planning for Spatially Constrained Systems

KAUST Repository

Xinyu Tang,

2010-01-25

Motion planning for spatially constrained robots is difficult due to additional constraints placed on the robot, such as closure constraints for closed chains or requirements on end-effector placement for articulated linkages. It is usually computationally too expensive to apply sampling-based planners to these problems since it is difficult to generate valid configurations. We overcome this challenge by redefining the robot\\'s degrees of freedom and constraints into a new set of parameters, called reachable distance space (RD-space), in which all configurations lie in the set of constraint-satisfying subspaces. This enables us to directly sample the constrained subspaces with complexity linear in the number of the robot\\'s degrees of freedom. In addition to supporting efficient sampling of configurations, we show that the RD-space formulation naturally supports planning and, in particular, we design a local planner suitable for use by sampling-based planners. We demonstrate the effectiveness and efficiency of our approach for several systems including closed chain planning with multiple loops, restricted end-effector sampling, and on-line planning for drawing/sculpting. We can sample single-loop closed chain systems with 1,000 links in time comparable to open chain sampling, and we can generate samples for 1,000-link multi-loop systems of varying topologies in less than a second. © 2010 The Author(s).

SU-G-BRA-08: Diaphragm Motion Tracking Based On KV CBCT Projections with a Constrained Linear Regression Optimization

Energy Technology Data Exchange (ETDEWEB)

Wei, J [City College of New York, New York, NY (United States); Chao, M [The Mount Sinai Medical Center, New York, NY (United States)

2016-06-15

Purpose: To develop a novel strategy to extract the respiratory motion of the thoracic diaphragm from kilovoltage cone beam computed tomography (CBCT) projections by a constrained linear regression optimization technique. Methods: A parabolic function was identified as the geometric model and was employed to fit the shape of the diaphragm on the CBCT projections. The search was initialized by five manually placed seeds on a pre-selected projection image. Temporal redundancies, the enabling phenomenology in video compression and encoding techniques, inherent in the dynamic properties of the diaphragm motion together with the geometrical shape of the diaphragm boundary and the associated algebraic constraint that significantly reduced the searching space of viable parabolic parameters was integrated, which can be effectively optimized by a constrained linear regression approach on the subsequent projections. The innovative algebraic constraints stipulating the kinetic range of the motion and the spatial constraint preventing any unphysical deviations was able to obtain the optimal contour of the diaphragm with minimal initialization. The algorithm was assessed by a fluoroscopic movie acquired at anteriorposterior fixed direction and kilovoltage CBCT projection image sets from four lung and two liver patients. The automatic tracing by the proposed algorithm and manual tracking by a human operator were compared in both space and frequency domains. Results: The error between the estimated and manual detections for the fluoroscopic movie was 0.54mm with standard deviation (SD) of 0.45mm, while the average error for the CBCT projections was 0.79mm with SD of 0.64mm for all enrolled patients. The submillimeter accuracy outcome exhibits the promise of the proposed constrained linear regression approach to track the diaphragm motion on rotational projection images. Conclusion: The new algorithm will provide a potential solution to rendering diaphragm motion and ultimately

SU-G-BRA-08: Diaphragm Motion Tracking Based On KV CBCT Projections with a Constrained Linear Regression Optimization

International Nuclear Information System (INIS)

Wei, J; Chao, M

2016-01-01

Purpose: To develop a novel strategy to extract the respiratory motion of the thoracic diaphragm from kilovoltage cone beam computed tomography (CBCT) projections by a constrained linear regression optimization technique. Methods: A parabolic function was identified as the geometric model and was employed to fit the shape of the diaphragm on the CBCT projections. The search was initialized by five manually placed seeds on a pre-selected projection image. Temporal redundancies, the enabling phenomenology in video compression and encoding techniques, inherent in the dynamic properties of the diaphragm motion together with the geometrical shape of the diaphragm boundary and the associated algebraic constraint that significantly reduced the searching space of viable parabolic parameters was integrated, which can be effectively optimized by a constrained linear regression approach on the subsequent projections. The innovative algebraic constraints stipulating the kinetic range of the motion and the spatial constraint preventing any unphysical deviations was able to obtain the optimal contour of the diaphragm with minimal initialization. The algorithm was assessed by a fluoroscopic movie acquired at anteriorposterior fixed direction and kilovoltage CBCT projection image sets from four lung and two liver patients. The automatic tracing by the proposed algorithm and manual tracking by a human operator were compared in both space and frequency domains. Results: The error between the estimated and manual detections for the fluoroscopic movie was 0.54mm with standard deviation (SD) of 0.45mm, while the average error for the CBCT projections was 0.79mm with SD of 0.64mm for all enrolled patients. The submillimeter accuracy outcome exhibits the promise of the proposed constrained linear regression approach to track the diaphragm motion on rotational projection images. Conclusion: The new algorithm will provide a potential solution to rendering diaphragm motion and ultimately

A Reduced Dantzig-Wolfe Decomposition for a Suboptimal Linear MPC

DEFF Research Database (Denmark)

Standardi, Laura; Poulsen, Niels Kjølstad; Jørgensen, John Bagterp

2014-01-01

Linear Model Predictive Control (MPC) is an efficient control technique that repeatedly solves online constrained linear programs. In this work we propose an economic linear MPC strategy for operation of energy systems consisting of multiple and independent power units. These systems cooperate...

Finite-time convergent recurrent neural network with a hard-limiting activation function for constrained optimization with piecewise-linear objective functions.

Science.gov (United States)

Liu, Qingshan; Wang, Jun

2011-04-01

This paper presents a one-layer recurrent neural network for solving a class of constrained nonsmooth optimization problems with piecewise-linear objective functions. The proposed neural network is guaranteed to be globally convergent in finite time to the optimal solutions under a mild condition on a derived lower bound of a single gain parameter in the model. The number of neurons in the neural network is the same as the number of decision variables of the optimization problem. Compared with existing neural networks for optimization, the proposed neural network has a couple of salient features such as finite-time convergence and a low model complexity. Specific models for two important special cases, namely, linear programming and nonsmooth optimization, are also presented. In addition, applications to the shortest path problem and constrained least absolute deviation problem are discussed with simulation results to demonstrate the effectiveness and characteristics of the proposed neural network.

Adaptive Constrained Optimal Control Design for Data-Based Nonlinear Discrete-Time Systems With Critic-Only Structure.

Science.gov (United States)

Luo, Biao; Liu, Derong; Wu, Huai-Ning

2018-06-01

Reinforcement learning has proved to be a powerful tool to solve optimal control problems over the past few years. However, the data-based constrained optimal control problem of nonaffine nonlinear discrete-time systems has rarely been studied yet. To solve this problem, an adaptive optimal control approach is developed by using the value iteration-based Q-learning (VIQL) with the critic-only structure. Most of the existing constrained control methods require the use of a certain performance index and only suit for linear or affine nonlinear systems, which is unreasonable in practice. To overcome this problem, the system transformation is first introduced with the general performance index. Then, the constrained optimal control problem is converted to an unconstrained optimal control problem. By introducing the action-state value function, i.e., Q-function, the VIQL algorithm is proposed to learn the optimal Q-function of the data-based unconstrained optimal control problem. The convergence results of the VIQL algorithm are established with an easy-to-realize initial condition . To implement the VIQL algorithm, the critic-only structure is developed, where only one neural network is required to approximate the Q-function. The converged Q-function obtained from the critic-only VIQL method is employed to design the adaptive constrained optimal controller based on the gradient descent scheme. Finally, the effectiveness of the developed adaptive control method is tested on three examples with computer simulation.

Non linear system become linear system

Directory of Open Access Journals (Sweden)

Petre Bucur

2007-01-01

Full Text Available The present paper refers to the theory and the practice of the systems regarding non-linear systems and their applications. We aimed the integration of these systems to elaborate their response as well as to highlight some outstanding features.

A Primal-Dual Interior Point-Linear Programming Algorithm for MPC

DEFF Research Database (Denmark)

Edlund, Kristian; Sokoler, Leo Emil; Jørgensen, John Bagterp

2009-01-01

Constrained optimal control problems for linear systems with linear constraints and an objective function consisting of linear and l1-norm terms can be expressed as linear programs. We develop an efficient primal-dual interior point algorithm for solution of such linear programs. The algorithm...

A Globally Convergent Parallel SSLE Algorithm for Inequality Constrained Optimization

Directory of Open Access Journals (Sweden)

Zhijun Luo

2014-01-01

Full Text Available A new parallel variable distribution algorithm based on interior point SSLE algorithm is proposed for solving inequality constrained optimization problems under the condition that the constraints are block-separable by the technology of sequential system of linear equation. Each iteration of this algorithm only needs to solve three systems of linear equations with the same coefficient matrix to obtain the descent direction. Furthermore, under certain conditions, the global convergence is achieved.

Event-triggered decentralized robust model predictive control for constrained large-scale interconnected systems

Directory of Open Access Journals (Sweden)

Ling Lu

2016-12-01

Full Text Available This paper considers the problem of event-triggered decentralized model predictive control (MPC for constrained large-scale linear systems subject to additive bounded disturbances. The constraint tightening method is utilized to formulate the MPC optimization problem. The local predictive control law for each subsystem is determined aperiodically by relevant triggering rule which allows a considerable reduction of the computational load. And then, the robust feasibility and closed-loop stability are proved and it is shown that every subsystem state will be driven into a robust invariant set. Finally, the effectiveness of the proposed approach is illustrated via numerical simulations.

A Linear Programming Approach to Routing Control in Networks of Constrained Nonlinear Positive Systems with Concave Flow Rates

Science.gov (United States)

Arneson, Heather M.; Dousse, Nicholas; Langbort, Cedric

2014-01-01

We consider control design for positive compartmental systems in which each compartment's outflow rate is described by a concave function of the amount of material in the compartment.We address the problem of determining the routing of material between compartments to satisfy time-varying state constraints while ensuring that material reaches its intended destination over a finite time horizon. We give sufficient conditions for the existence of a time-varying state-dependent routing strategy which ensures that the closed-loop system satisfies basic network properties of positivity, conservation and interconnection while ensuring that capacity constraints are satisfied, when possible, or adjusted if a solution cannot be found. These conditions are formulated as a linear programming problem. Instances of this linear programming problem can be solved iteratively to generate a solution to the finite horizon routing problem. Results are given for the application of this control design method to an example problem. Key words: linear programming; control of networks; positive systems; controller constraints and structure.

Affine Lie algebraic origin of constrained KP hierarchies

International Nuclear Information System (INIS)

Aratyn, H.; Gomes, J.F.; Zimerman, A.H.

1994-07-01

It is presented an affine sl(n+1) algebraic construction of the basic constrained KP hierarchy. This hierarchy is analyzed using two approaches, namely linear matrix eigenvalue problem on hermitian symmetric space and constrained KP Lax formulation and we show that these approaches are equivalent. The model is recognized to be generalized non-linear Schroedinger (GNLS) hierarchy and it is used as a building block for a new class of constrained KP hierarchies. These constrained KP hierarchies are connected via similarity-Backlund transformations and interpolate between GNLS and multi-boson KP-Toda hierarchies. The construction uncovers origin of the Toda lattice structure behind the latter hierarchy. (author). 23 refs

ROBUST MPC FOR STABLE LINEAR SYSTEMS

Directory of Open Access Journals (Sweden)

M.A. Rodrigues

2002-03-01

Full Text Available In this paper, a new model predictive controller (MPC, which is robust for a class of model uncertainties, is developed. Systems with stable dynamics and time-invariant model uncertainty are treated. The development herein proposed is focused on real industrial systems where the controller is part of an on-line optimization scheme and works in the output-tracking mode. In addition, the system has a time-varying number of degrees of freedom since some of the manipulated inputs may become constrained. Moreover, the number of controlled outputs may also vary during system operation. Consequently, the actual system may show operating conditions with a number of controlled outputs larger than the number of available manipulated inputs. The proposed controller uses a state-space model, which is aimed at the representation of the output-predicted trajectory. Based on this model, a cost function is proposed whereby the output error is integrated along an infinite prediction horizon. It is considered the case of multiple operating points, where the controller stabilizes a set of models corresponding to different operating conditions for the system. It is shown that closed-loop stability is guaranteed by the feasibility of a linear matrix optimization problem.

Constrained State Estimation for Individual Localization in Wireless Body Sensor Networks

Directory of Open Access Journals (Sweden)

Xiaoxue Feng

2014-11-01

Full Text Available Wireless body sensor networks based on ultra-wideband radio have recently received much research attention due to its wide applications in health-care, security, sports and entertainment. Accurate localization is a fundamental problem to realize the development of effective location-aware applications above. In this paper the problem of constrained state estimation for individual localization in wireless body sensor networks is addressed. Priori knowledge about geometry among the on-body nodes as additional constraint is incorporated into the traditional filtering system. The analytical expression of state estimation with linear constraint to exploit the additional information is derived. Furthermore, for nonlinear constraint, first-order and second-order linearizations via Taylor series expansion are proposed to transform the nonlinear constraint to the linear case. Examples between the first-order and second-order nonlinear constrained filters based on interacting multiple model extended kalman filter (IMM-EKF show that the second-order solution for higher order nonlinearity as present in this paper outperforms the first-order solution, and constrained IMM-EKF obtains superior estimation than IMM-EKF without constraint. Another brownian motion individual localization example also illustrates the effectiveness of constrained nonlinear iterative least square (NILS, which gets better filtering performance than NILS without constraint.

Constrained State Estimation for Individual Localization in Wireless Body Sensor Networks

Science.gov (United States)

Feng, Xiaoxue; Snoussi, Hichem; Liang, Yan; Jiao, Lianmeng

2014-01-01

Wireless body sensor networks based on ultra-wideband radio have recently received much research attention due to its wide applications in health-care, security, sports and entertainment. Accurate localization is a fundamental problem to realize the development of effective location-aware applications above. In this paper the problem of constrained state estimation for individual localization in wireless body sensor networks is addressed. Priori knowledge about geometry among the on-body nodes as additional constraint is incorporated into the traditional filtering system. The analytical expression of state estimation with linear constraint to exploit the additional information is derived. Furthermore, for nonlinear constraint, first-order and second-order linearizations via Taylor series expansion are proposed to transform the nonlinear constraint to the linear case. Examples between the first-order and second-order nonlinear constrained filters based on interacting multiple model extended kalman filter (IMM-EKF) show that the second-order solution for higher order nonlinearity as present in this paper outperforms the first-order solution, and constrained IMM-EKF obtains superior estimation than IMM-EKF without constraint. Another brownian motion individual localization example also illustrates the effectiveness of constrained nonlinear iterative least square (NILS), which gets better filtering performance than NILS without constraint. PMID:25390408

Constrained state estimation for individual localization in wireless body sensor networks.

Science.gov (United States)

Feng, Xiaoxue; Snoussi, Hichem; Liang, Yan; Jiao, Lianmeng

2014-11-10

Wireless body sensor networks based on ultra-wideband radio have recently received much research attention due to its wide applications in health-care, security, sports and entertainment. Accurate localization is a fundamental problem to realize the development of effective location-aware applications above. In this paper the problem of constrained state estimation for individual localization in wireless body sensor networks is addressed. Priori knowledge about geometry among the on-body nodes as additional constraint is incorporated into the traditional filtering system. The analytical expression of state estimation with linear constraint to exploit the additional information is derived. Furthermore, for nonlinear constraint, first-order and second-order linearizations via Taylor series expansion are proposed to transform the nonlinear constraint to the linear case. Examples between the first-order and second-order nonlinear constrained filters based on interacting multiple model extended kalman filter (IMM-EKF) show that the second-order solution for higher order nonlinearity as present in this paper outperforms the first-order solution, and constrained IMM-EKF obtains superior estimation than IMM-EKF without constraint. Another brownian motion individual localization example also illustrates the effectiveness of constrained nonlinear iterative least square (NILS), which gets better filtering performance than NILS without constraint.

Robust Management of Combined Heat and Power Systems via Linear Decision Rules

DEFF Research Database (Denmark)

Zugno, Marco; Morales González, Juan Miguel; Madsen, Henrik

2014-01-01

The heat and power outputs of Combined Heat and Power (CHP) units are jointly constrained. Hence, the optimal management of systems including CHP units is a multicommodity optimization problem. Problems of this type are stochastic, owing to the uncertainty inherent both in the demand for heat and...... linear decision rules to guarantee both tractability and a correct representation of the dynamic aspects of the problem. Numerical results from an illustrative example confirm the value of the proposed approach....

Secure Fusion Estimation for Bandwidth Constrained Cyber-Physical Systems Under Replay Attacks.

Science.gov (United States)

Chen, Bo; Ho, Daniel W C; Hu, Guoqiang; Yu, Li; Bo Chen; Ho, Daniel W C; Guoqiang Hu; Li Yu; Chen, Bo; Ho, Daniel W C; Hu, Guoqiang; Yu, Li

2018-06-01

State estimation plays an essential role in the monitoring and supervision of cyber-physical systems (CPSs), and its importance has made the security and estimation performance a major concern. In this case, multisensor information fusion estimation (MIFE) provides an attractive alternative to study secure estimation problems because MIFE can potentially improve estimation accuracy and enhance reliability and robustness against attacks. From the perspective of the defender, the secure distributed Kalman fusion estimation problem is investigated in this paper for a class of CPSs under replay attacks, where each local estimate obtained by the sink node is transmitted to a remote fusion center through bandwidth constrained communication channels. A new mathematical model with compensation strategy is proposed to characterize the replay attacks and bandwidth constrains, and then a recursive distributed Kalman fusion estimator (DKFE) is designed in the linear minimum variance sense. According to different communication frameworks, two classes of data compression and compensation algorithms are developed such that the DKFEs can achieve the desired performance. Several attack-dependent and bandwidth-dependent conditions are derived such that the DKFEs are secure under replay attacks. An illustrative example is given to demonstrate the effectiveness of the proposed methods.

Minimal constrained supergravity

Energy Technology Data Exchange (ETDEWEB)

Cribiori, N. [Dipartimento di Fisica e Astronomia “Galileo Galilei”, Università di Padova, Via Marzolo 8, 35131 Padova (Italy); INFN, Sezione di Padova, Via Marzolo 8, 35131 Padova (Italy); Dall' Agata, G., E-mail: dallagat@pd.infn.it [Dipartimento di Fisica e Astronomia “Galileo Galilei”, Università di Padova, Via Marzolo 8, 35131 Padova (Italy); INFN, Sezione di Padova, Via Marzolo 8, 35131 Padova (Italy); Farakos, F. [Dipartimento di Fisica e Astronomia “Galileo Galilei”, Università di Padova, Via Marzolo 8, 35131 Padova (Italy); INFN, Sezione di Padova, Via Marzolo 8, 35131 Padova (Italy); Porrati, M. [Center for Cosmology and Particle Physics, Department of Physics, New York University, 4 Washington Place, New York, NY 10003 (United States)

2017-01-10

We describe minimal supergravity models where supersymmetry is non-linearly realized via constrained superfields. We show that the resulting actions differ from the so called “de Sitter” supergravities because we consider constraints eliminating directly the auxiliary fields of the gravity multiplet.

Minimal constrained supergravity

International Nuclear Information System (INIS)

Cribiori, N.; Dall'Agata, G.; Farakos, F.; Porrati, M.

2017-01-01

We describe minimal supergravity models where supersymmetry is non-linearly realized via constrained superfields. We show that the resulting actions differ from the so called “de Sitter” supergravities because we consider constraints eliminating directly the auxiliary fields of the gravity multiplet.

Minimal constrained supergravity

Directory of Open Access Journals (Sweden)

N. Cribiori

2017-01-01

Full Text Available We describe minimal supergravity models where supersymmetry is non-linearly realized via constrained superfields. We show that the resulting actions differ from the so called “de Sitter” supergravities because we consider constraints eliminating directly the auxiliary fields of the gravity multiplet.

Constrained Optimal Stochastic Control of Non-Linear Wave Energy Point Absorbers

DEFF Research Database (Denmark)

Sichani, Mahdi Teimouri; Chen, Jian-Bing; Kramer, Morten

2014-01-01

to extract energy. Constrains are enforced on the control force to prevent large structural stresses in the floater at specific hot spots with the risk of inducing fatigue damage, or because the demanded control force cannot be supplied by the actuator system due to saturation. Further, constraints...... are enforced on the motion of the floater to prevent it from hitting the bottom of the sea or to make unacceptable jumps out of the water. The applied control law, which is of the feedback type with feedback from the displacement, velocity, and acceleration of the floater, contains two unprovided gain...

Linearization of the Lorenz system

International Nuclear Information System (INIS)

Li, Chunbiao; Sprott, Julien Clinton; Thio, Wesley

2015-01-01

A partial and complete piecewise linearized version of the Lorenz system is proposed. The linearized versions have an independent total amplitude control parameter. Additional further linearization leads naturally to a piecewise linear version of the diffusionless Lorenz system. A chaotic circuit with a single amplitude controller is then implemented using a new switch element, producing a chaotic oscillation that agrees with the numerical calculation for the piecewise linear diffusionless Lorenz system. - Highlights: • A partial and complete piecewise linearized version of the Lorenz system are addressed. • The linearized versions have an independent total amplitude control parameter. • A piecewise linear version of the diffusionless Lorenz system is derived by further linearization. • A corresponding chaotic circuit without any multiplier is implemented for the chaotic oscillation

Linearization of the Lorenz system

Energy Technology Data Exchange (ETDEWEB)

Li, Chunbiao, E-mail: goontry@126.com [School of Electronic & Information Engineering, Nanjing University of Information Science & Technology, Nanjing 210044 (China); Engineering Technology Research and Development Center of Jiangsu Circulation Modernization Sensor Network, Jiangsu Institute of Commerce, Nanjing 211168 (China); Sprott, Julien Clinton [Department of Physics, University of Wisconsin–Madison, Madison, WI 53706 (United States); Thio, Wesley [Department of Electrical and Computer Engineering, The Ohio State University, Columbus, OH 43210 (United States)

2015-05-08

A partial and complete piecewise linearized version of the Lorenz system is proposed. The linearized versions have an independent total amplitude control parameter. Additional further linearization leads naturally to a piecewise linear version of the diffusionless Lorenz system. A chaotic circuit with a single amplitude controller is then implemented using a new switch element, producing a chaotic oscillation that agrees with the numerical calculation for the piecewise linear diffusionless Lorenz system. - Highlights: • A partial and complete piecewise linearized version of the Lorenz system are addressed. • The linearized versions have an independent total amplitude control parameter. • A piecewise linear version of the diffusionless Lorenz system is derived by further linearization. • A corresponding chaotic circuit without any multiplier is implemented for the chaotic oscillation.

Multiplicative algorithms for constrained non-negative matrix factorization

KAUST Repository

Peng, Chengbin

2012-12-01

Non-negative matrix factorization (NMF) provides the advantage of parts-based data representation through additive only combinations. It has been widely adopted in areas like item recommending, text mining, data clustering, speech denoising, etc. In this paper, we provide an algorithm that allows the factorization to have linear or approximatly linear constraints with respect to each factor. We prove that if the constraint function is linear, algorithms within our multiplicative framework will converge. This theory supports a large variety of equality and inequality constraints, and can facilitate application of NMF to a much larger domain. Taking the recommender system as an example, we demonstrate how a specialized weighted and constrained NMF algorithm can be developed to fit exactly for the problem, and the tests justify that our constraints improve the performance for both weighted and unweighted NMF algorithms under several different metrics. In particular, on the Movielens data with 94% of items, the Constrained NMF improves recall rate 3% compared to SVD50 and 45% compared to SVD150, which were reported as the best two in the top-N metric. © 2012 IEEE.

Stability Criterion of Linear Stochastic Systems Subject to Mixed H2/Passivity Performance

Directory of Open Access Journals (Sweden)

Cheung-Chieh Ku

2015-01-01

Full Text Available The H2 control scheme and passivity theory are applied to investigate the stability criterion of continuous-time linear stochastic system subject to mixed performance. Based on the stochastic differential equation, the stochastic behaviors can be described as multiplicative noise terms. For the considered system, the H2 control scheme is applied to deal with the problem on minimizing output energy. And the asymptotical stability of the system can be guaranteed under desired initial conditions. Besides, the passivity theory is employed to constrain the effect of external disturbance on the system. Moreover, the Itô formula and Lyapunov function are used to derive the sufficient conditions which are converted into linear matrix inequality (LMI form for applying convex optimization algorithm. Via solving the sufficient conditions, the state feedback controller can be established such that the asymptotical stability and mixed performance of the system are achieved in the mean square. Finally, the synchronous generator system is used to verify the effectiveness and applicability of the proposed design method.

Quantization of soluble classical constrained systems

International Nuclear Information System (INIS)

Belhadi, Z.; Menas, F.; Bérard, A.; Mohrbach, H.

2014-01-01

The derivation of the brackets among coordinates and momenta for classical constrained systems is a necessary step toward their quantization. Here we present a new approach for the determination of the classical brackets which does neither require Dirac’s formalism nor the symplectic method of Faddeev and Jackiw. This approach is based on the computation of the brackets between the constants of integration of the exact solutions of the equations of motion. From them all brackets of the dynamical variables of the system can be deduced in a straightforward way

Quantization of soluble classical constrained systems

Energy Technology Data Exchange (ETDEWEB)

Belhadi, Z. [Laboratoire de physique et chimie quantique, Faculté des sciences, Université Mouloud Mammeri, BP 17, 15000 Tizi Ouzou (Algeria); Laboratoire de physique théorique, Faculté des sciences exactes, Université de Bejaia, 06000 Bejaia (Algeria); Menas, F. [Laboratoire de physique et chimie quantique, Faculté des sciences, Université Mouloud Mammeri, BP 17, 15000 Tizi Ouzou (Algeria); Ecole Nationale Préparatoire aux Etudes d’ingéniorat, Laboratoire de physique, RN 5 Rouiba, Alger (Algeria); Bérard, A. [Equipe BioPhysStat, Laboratoire LCP-A2MC, ICPMB, IF CNRS No 2843, Université de Lorraine, 1 Bd Arago, 57078 Metz Cedex (France); Mohrbach, H., E-mail: herve.mohrbach@univ-lorraine.fr [Equipe BioPhysStat, Laboratoire LCP-A2MC, ICPMB, IF CNRS No 2843, Université de Lorraine, 1 Bd Arago, 57078 Metz Cedex (France)

2014-12-15

The derivation of the brackets among coordinates and momenta for classical constrained systems is a necessary step toward their quantization. Here we present a new approach for the determination of the classical brackets which does neither require Dirac’s formalism nor the symplectic method of Faddeev and Jackiw. This approach is based on the computation of the brackets between the constants of integration of the exact solutions of the equations of motion. From them all brackets of the dynamical variables of the system can be deduced in a straightforward way.

Lagrangian formalism for constrained systems. 2. Gauge symmetries

International Nuclear Information System (INIS)

Pyatov, P.N.

1990-01-01

Using the Lagrangian formalism for constrained systems all gauge symmetries peculiar for a given Lagrangian system and in establishing the relation between them and the constraints are constructed. Besides, the question about the possible dependence of gauge transformations on accelerations and other higher order time derivatives of coordinates is clarified. 14 refs

Incomplete Dirac reduction of constrained Hamiltonian systems

Energy Technology Data Exchange (ETDEWEB)

Chandre, C., E-mail: chandre@cpt.univ-mrs.fr

2015-10-15

First-class constraints constitute a potential obstacle to the computation of a Poisson bracket in Dirac’s theory of constrained Hamiltonian systems. Using the pseudoinverse instead of the inverse of the matrix defined by the Poisson brackets between the constraints, we show that a Dirac–Poisson bracket can be constructed, even if it corresponds to an incomplete reduction of the original Hamiltonian system. The uniqueness of Dirac brackets is discussed. The relevance of this procedure for infinite dimensional Hamiltonian systems is exemplified.

Linear system theory

Science.gov (United States)

Callier, Frank M.; Desoer, Charles A.

1991-01-01

The aim of this book is to provide a systematic and rigorous access to the main topics of linear state-space system theory in both the continuous-time case and the discrete-time case; and the I/O description of linear systems. The main thrusts of the work are the analysis of system descriptions and derivations of their properties, LQ-optimal control, state feedback and state estimation, and MIMO unity-feedback systems.

Coherent states in constrained systems

International Nuclear Information System (INIS)

Nakamura, M.; Kojima, K.

2001-01-01

When quantizing the constrained systems, there often arise the quantum corrections due to the non-commutativity in the re-ordering of constraint operators in the products of operators. In the bosonic second-class constraints, furthermore, the quantum corrections caused by the uncertainty principle should be taken into account. In order to treat these corrections simultaneously, the alternative projection technique of operators is proposed by introducing the available minimal uncertainty states of the constraint operators. Using this projection technique together with the projection operator method (POM), these two kinds of quantum corrections were investigated

Modelos lineares e não lineares inteiros para problemas da mochila bidimensional restrita a 2 estágios Linear and nonlinear integer models for constrained two-stage two-dimensional knapsack problems

Directory of Open Access Journals (Sweden)

Horacio Hideki Yanasse

2013-01-01

Full Text Available Neste trabalho revemos alguns modelos lineares e não lineares inteiros para gerar padrões de corte bidimensionais guilhotinados de 2 estágios, incluindo os casos exato e não exato e restrito e irrestrito. Esses problemas são casos particulares do problema da mochila bidimensional. Apresentamos também novos modelos para gerar esses padrões de corte, baseados em adaptações ou extensões de modelos para gerar padrões de corte bidimensionais restritos 1-grupo. Padrões 2 estágios aparecem em diferentes processos de corte, como, por exemplo, em indústrias de móveis e de chapas de madeira. Os modelos são úteis para a pesquisa e o desenvolvimento de métodos de solução mais eficientes, explorando estruturas particulares, a decomposição do modelo, relaxações do modelo etc. Eles também são úteis para a avaliação do desempenho de heurísticas, já que permitem (pelo menos para problemas de tamanho moderado uma estimativa do gap de otimalidade de soluções obtidas por heurísticas. Para ilustrar a aplicação dos modelos, analisamos os resultados de alguns experimentos computacionais com exemplos da literatura e outros gerados aleatoriamente. Os resultados foram produzidos usando um software comercial conhecido e mostram que o esforço computacional necessário para resolver os modelos pode ser bastante diferente.In this work we review some linear and nonlinear integer models to generate two stage two-dimensional guillotine cutting patterns, including the constrained, non constrained, exact and non exact cases. These problems are particular cases of the two dimensional knapsack problems. We also present new models to generate these cutting patterns, based on adaptations and extensions of models that generate one-group constrained two dimensional cutting patterns. Two stage patterns arise in different cutting processes like, for instance, in the furniture industry and wooden hardboards. The models are useful for the research and

Linear optical response of finite systems using multishift linear system solvers

Energy Technology Data Exchange (ETDEWEB)

Hübener, Hannes; Giustino, Feliciano [Department of Materials, University of Oxford, Oxford OX1 3PH (United Kingdom)

2014-07-28

We discuss the application of multishift linear system solvers to linear-response time-dependent density functional theory. Using this technique the complete frequency-dependent electronic density response of finite systems to an external perturbation can be calculated at the cost of a single solution of a linear system via conjugate gradients. We show that multishift time-dependent density functional theory yields excitation energies and oscillator strengths in perfect agreement with the standard diagonalization of the response matrix (Casida's method), while being computationally advantageous. We present test calculations for benzene, porphin, and chlorophyll molecules. We argue that multishift solvers may find broad applicability in the context of excited-state calculations within density-functional theory and beyond.

Positive solution of non-square fully Fuzzy linear system of equation in general form using least square method

Directory of Open Access Journals (Sweden)

Reza Ezzati

2014-08-01

Full Text Available In this paper, we propose the least square method for computing the positive solution of a non-square fully fuzzy linear system. To this end, we use Kaffman' arithmetic operations on fuzzy numbers \\cite{17}. Here, considered existence of exact solution using pseudoinverse, if they are not satisfy in positive solution condition, we will compute fuzzy vector core and then we will obtain right and left spreads of positive fuzzy vector by introducing constrained least squares problem. Using our proposed method, non-square fully fuzzy linear system of equations always has a solution. Finally, we illustrate the efficiency of proposed method by solving some numerical examples.

Pole shifting with constrained output feedback

International Nuclear Information System (INIS)

Hamel, D.; Mensah, S.; Boisvert, J.

1984-03-01

The concept of pole placement plays an important role in linear, multi-variable, control theory. It has received much attention since its introduction, and several pole shifting algorithms are now available. This work presents a new method which allows practical and engineering constraints such as gain limitation and controller structure to be introduced right into the pole shifting design strategy. This is achieved by formulating the pole placement problem as a constrained optimization problem. Explicit constraints (controller structure and gain limits) are defined to identify an admissible region for the feedback gain matrix. The desired pole configuration is translated into an appropriate cost function which must be closed-loop minimized. The resulting constrained optimization problem can thus be solved with optimization algorithms. The method has been implemented as an algorithmic interactive module in a computer-aided control system design package, MVPACK. The application of the method is illustrated to design controllers for an aircraft and an evaporator. The results illustrate the importance of controller structure on overall performance of a control system

Quasicanonical structure of optimal control in constrained discrete systems

Science.gov (United States)

Sieniutycz, S.

2003-06-01

This paper considers discrete processes governed by difference rather than differential equations for the state transformation. The basic question asked is if and when Hamiltonian canonical structures are possible in optimal discrete systems. Considering constrained discrete control, general optimization algorithms are derived that constitute suitable theoretical and computational tools when evaluating extremum properties of constrained physical models. The mathematical basis of the general theory is the Bellman method of dynamic programming (DP) and its extension in the form of the so-called Carathéodory-Boltyanski (CB) stage criterion which allows a variation of the terminal state that is otherwise fixed in the Bellman's method. Two relatively unknown, powerful optimization algorithms are obtained: an unconventional discrete formalism of optimization based on a Hamiltonian for multistage systems with unconstrained intervals of holdup time, and the time interval constrained extension of the formalism. These results are general; namely, one arrives at: the discrete canonical Hamilton equations, maximum principles, and (at the continuous limit of processes with free intervals of time) the classical Hamilton-Jacobi theory along with all basic results of variational calculus. Vast spectrum of applications of the theory is briefly discussed.

Model Predictive Control for Linear Complementarity and Extended Linear Complementarity Systems

Directory of Open Access Journals (Sweden)

Bambang Riyanto

2005-11-01

Full Text Available In this paper, we propose model predictive control method for linear complementarity and extended linear complementarity systems by formulating optimization along prediction horizon as mixed integer quadratic program. Such systems contain interaction between continuous dynamics and discrete event systems, and therefore, can be categorized as hybrid systems. As linear complementarity and extended linear complementarity systems finds applications in different research areas, such as impact mechanical systems, traffic control and process control, this work will contribute to the development of control design method for those areas as well, as shown by three given examples.

Small-kernel constrained-least-squares restoration of sampled image data

Science.gov (United States)

Hazra, Rajeeb; Park, Stephen K.

1992-10-01

Constrained least-squares image restoration, first proposed by Hunt twenty years ago, is a linear image restoration technique in which the restoration filter is derived by maximizing the smoothness of the restored image while satisfying a fidelity constraint related to how well the restored image matches the actual data. The traditional derivation and implementation of the constrained least-squares restoration filter is based on an incomplete discrete/discrete system model which does not account for the effects of spatial sampling and image reconstruction. For many imaging systems, these effects are significant and should not be ignored. In a recent paper Park demonstrated that a derivation of the Wiener filter based on the incomplete discrete/discrete model can be extended to a more comprehensive end-to-end, continuous/discrete/continuous model. In a similar way, in this paper, we show that a derivation of the constrained least-squares filter based on the discrete/discrete model can also be extended to this more comprehensive continuous/discrete/continuous model and, by so doing, an improved restoration filter is derived. Building on previous work by Reichenbach and Park for the Wiener filter, we also show that this improved constrained least-squares restoration filter can be efficiently implemented as a small-kernel convolution in the spatial domain.

Closed-form solutions for linear regulator-design of mechanical systems including optimal weighting matrix selection

Science.gov (United States)

Hanks, Brantley R.; Skelton, Robert E.

1991-01-01

This paper addresses the restriction of Linear Quadratic Regulator (LQR) solutions to the algebraic Riccati Equation to design spaces which can be implemented as passive structural members and/or dampers. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical systems. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist. Some examples of simple spring mass systems are shown to illustrate key points.

Null Space Integration Method for Constrained Multibody Systems with No Constraint Violation

International Nuclear Information System (INIS)

Terze, Zdravko; Lefeber, Dirk; Muftic, Osman

2001-01-01

A method for integrating equations of motion of constrained multibody systems with no constraint violation is presented. A mathematical model, shaped as a differential-algebraic system of index 1, is transformed into a system of ordinary differential equations using the null-space projection method. Equations of motion are set in a non-minimal form. During integration, violations of constraints are corrected by solving constraint equations at the position and velocity level, utilizing the metric of the system's configuration space, and projective criterion to the coordinate partitioning method. The method is applied to dynamic simulation of 3D constrained biomechanical system. The simulation results are evaluated by comparing them to the values of characteristic parameters obtained by kinematics analysis of analyzed motion based unmeasured kinematics data

Systems analysis on the condition of market penetration for hydrogen technologies using linear programming model

International Nuclear Information System (INIS)

Kato, K.; Ihara, S.

1993-01-01

Hydrogen is expected to be an important energy carrier, especially in the frame of global warming problem solution. The purpose of this study is to examine the condition of market penetration of hydrogen technologies in reducing CO 2 emissions. A multi-time-period linear programming model (MARKAL, Market Allocation)) is used to explore technology options and cost for meeting the energy demands while reducing CO 2 emissions from energy systems. The results show that hydrogen technologies become economical when CO 2 emissions are stringently constrained. 9 figs., 2 refs

Splines and polynomial tools for flatness-based constrained motion planning

Science.gov (United States)

Suryawan, Fajar; De Doná, José; Seron, María

2012-08-01

This article addresses the problem of trajectory planning for flat systems with constraints. Flat systems have the useful property that the input and the state can be completely characterised by the so-called flat output. We propose a spline parametrisation for the flat output, the performance output, the states and the inputs. Using this parametrisation the problem of constrained trajectory planning can be cast into a simple quadratic programming problem. An important result is that the B-spline parametrisation used gives exact results for constrained linear continuous-time system. The result is exact in the sense that the constrained signal can be made arbitrarily close to the boundary without having intersampling issues (as one would have in sampled-data systems). Simulation examples are presented, involving the generation of rest-to-rest trajectories. In addition, an experimental result of the method is also presented, where two methods to generate trajectories for a magnetic-levitation (maglev) system in the presence of constraints are compared and each method's performance is discussed. The first method uses the nonlinear model of the plant, which turns out to belong to the class of flat systems. The second method uses a linearised version of the plant model around an operating point. In every case, a continuous-time description is used. The experimental results on a real maglev system reported here show that, in most scenarios, the nonlinear and linearised models produce almost similar, indistinguishable trajectories.

Constraining new physics models with isotope shift spectroscopy

Science.gov (United States)

Frugiuele, Claudia; Fuchs, Elina; Perez, Gilad; Schlaffer, Matthias

2017-07-01

Isotope shifts of transition frequencies in atoms constrain generic long- and intermediate-range interactions. We focus on new physics scenarios that can be most strongly constrained by King linearity violation such as models with B -L vector bosons, the Higgs portal, and chameleon models. With the anticipated precision, King linearity violation has the potential to set the strongest laboratory bounds on these models in some regions of parameter space. Furthermore, we show that this method can probe the couplings relevant for the protophobic interpretation of the recently reported Be anomaly. We extend the formalism to include an arbitrary number of transitions and isotope pairs and fit the new physics coupling to the currently available isotope shift measurements.

A Simply Constrained Optimization Reformulation of KKT Systems Arising from Variational Inequalities

International Nuclear Information System (INIS)

Facchinei, F.; Fischer, A.; Kanzow, C.; Peng, J.-M.

1999-01-01

The Karush-Kuhn-Tucker (KKT) conditions can be regarded as optimality conditions for both variational inequalities and constrained optimization problems. In order to overcome some drawbacks of recently proposed reformulations of KKT systems, we propose casting KKT systems as a minimization problem with nonnegativity constraints on some of the variables. We prove that, under fairly mild assumptions, every stationary point of this constrained minimization problem is a solution of the KKT conditions. Based on this reformulation, a new algorithm for the solution of the KKT conditions is suggested and shown to have some strong global and local convergence properties

Feedback systems for linear colliders

CERN Document Server

Hendrickson, L; Himel, Thomas M; Minty, Michiko G; Phinney, N; Raimondi, Pantaleo; Raubenheimer, T O; Shoaee, H; Tenenbaum, P G

1999-01-01

Feedback systems are essential for stable operation of a linear collider, providing a cost-effective method for relaxing tight tolerances. In the Stanford Linear Collider (SLC), feedback controls beam parameters such as trajectory, energy, and intensity throughout the accelerator. A novel dithering optimization system which adjusts final focus parameters to maximize luminosity contributed to achieving record performance in the 1997-98 run. Performance limitations of the steering feedback have been investigated, and improvements have been made. For the Next Linear Collider (NLC), extensive feedback systems are planned as an intregal part of the design. Feedback requiremetns for JLC (the Japanese Linear Collider) are essentially identical to NLC; some of the TESLA requirements are similar but there are significant differences. For NLC, algorithms which incorporate improvements upon the SLC implementation are being prototyped. Specialized systems for the damping rings, rf and interaction point will operate at hi...

On gauge fixing and quantization of constrained Hamiltonian systems

International Nuclear Information System (INIS)

Dayi, O.F.

1989-06-01

In constrained Hamiltonian systems which possess first class constraints some subsidiary conditions should be imposed for detecting physical observables. This issue and quantization of the system are clarified. It is argued that the reduced phase space and Dirac method of quantization, generally, differ only in the definition of the Hilbert space one should use. For the dynamical systems possessing second class constraints the definition of physical Hilbert space in the BFV-BRST operator quantization method is different from the usual definition. (author). 18 refs

Feedback Systems for Linear Colliders

International Nuclear Information System (INIS)

1999-01-01

Feedback systems are essential for stable operation of a linear collider, providing a cost-effective method for relaxing tight tolerances. In the Stanford Linear Collider (SLC), feedback controls beam parameters such as trajectory, energy, and intensity throughout the accelerator. A novel dithering optimization system which adjusts final focus parameters to maximize luminosity contributed to achieving record performance in the 1997-98 run. Performance limitations of the steering feedback have been investigated, and improvements have been made. For the Next Linear Collider (NLC), extensive feedback systems are planned as an integral part of the design. Feedback requirements for JLC (the Japanese Linear Collider) are essentially identical to NLC; some of the TESLA requirements are similar but there are significant differences. For NLC, algorithms which incorporate improvements upon the SLC implementation are being prototyped. Specialized systems for the damping rings, rf and interaction point will operate at high bandwidth and fast response. To correct for the motion of individual bunches within a train, both feedforward and feedback systems are planned. SLC experience has shown that feedback systems are an invaluable operational tool for decoupling systems, allowing precision tuning, and providing pulse-to-pulse diagnostics. Feedback systems for the NLC will incorporate the key SLC features and the benefits of advancing technologies

Application of a Double-Sided Chance-Constrained Integer Linear Program for Optimization of the Incremental Value of Ecosystem Services in Jilin Province, China

Directory of Open Access Journals (Sweden)

Baofeng Cai

2017-08-01

Full Text Available The Interconnected River System Network Project (IRSNP is a significant water supply engineering project, which is capable of effectively utilizing flood resources to generate ecological value, by connecting 198 lakes and ponds in western Jilin, northeast China. In this article, an optimization research approach has been proposed to maximize the incremental value of IRSNP ecosystem services. A double-sided chance-constrained integer linear program (DCCILP method has been proposed to support the optimization, which can deal with uncertainties presented as integers or random parameters that appear on both sides of the decision variable at the same time. The optimal scheme indicates that after rational optimization, the total incremental value of ecosystem services from the interconnected river system network project increased 22.25%, providing an increase in benefits of 3.26 × 109 ¥ compared to the original scheme. Most of the functional area is swamp wetland, which provides the greatest ecological benefits. Adjustment services increased obviously, implying that the optimization scheme prioritizes ecological benefits rather than supply and production services.

Standard diffusive systems are well-posed linear systems

NARCIS (Netherlands)

Matignon, Denis; Zwart, Heiko J.

2004-01-01

The class of well-posed linear systems as introduced by Salamon has become a well-understood class of systems, see e.g. the work of Weiss and the book of Staffans. Many partial partial differential equations with boundary control and point observation can be formulated as a well-posed linear system.

Constrained quadratic stabilization of discrete-time uncertain nonlinear multi-model systems using piecewise affine state-feedback

Directory of Open Access Journals (Sweden)

Olav Slupphaug

1999-07-01

Full Text Available In this paper a method for nonlinear robust stabilization based on solving a bilinear matrix inequality (BMI feasibility problem is developed. Robustness against model uncertainty is handled. In different non-overlapping regions of the state-space called clusters the plant is assumed to be an element in a polytope which vertices (local models are affine systems. In the clusters containing the origin in their closure, the local models are restricted to be linear systems. The clusters cover the region of interest in the state-space. An affine state-feedback is associated with each cluster. By utilizing the affinity of the local models and the state-feedback, a set of linear matrix inequalities (LMIs combined with a single nonconvex BMI are obtained which, if feasible, guarantee quadratic stability of the origin of the closed-loop. The feasibility problem is attacked by a branch-and-bound based global approach. If the feasibility check is successful, the Liapunov matrix and the piecewise affine state-feedback are given directly by the feasible solution. Control constraints are shown to be representable by LMIs or BMIs, and an application of the control design method to robustify constrained nonlinear model predictive control is presented. Also, the control design method is applied to a simple example.

A feasible DY conjugate gradient method for linear equality constraints

Science.gov (United States)

LI, Can

2017-09-01

In this paper, we propose a feasible conjugate gradient method for solving linear equality constrained optimization problem. The method is an extension of the Dai-Yuan conjugate gradient method proposed by Dai and Yuan to linear equality constrained optimization problem. It can be applied to solve large linear equality constrained problem due to lower storage requirement. An attractive property of the method is that the generated direction is always feasible and descent direction. Under mild conditions, the global convergence of the proposed method with exact line search is established. Numerical experiments are also given which show the efficiency of the method.

A New Continuous-Time Equality-Constrained Optimization to Avoid Singularity.

Science.gov (United States)

Quan, Quan; Cai, Kai-Yuan

2016-02-01

In equality-constrained optimization, a standard regularity assumption is often associated with feasible point methods, namely, that the gradients of constraints are linearly independent. In practice, the regularity assumption may be violated. In order to avoid such a singularity, a new projection matrix is proposed based on which a feasible point method to continuous-time, equality-constrained optimization is developed. First, the equality constraint is transformed into a continuous-time dynamical system with solutions that always satisfy the equality constraint. Second, a new projection matrix without singularity is proposed to realize the transformation. An update (or say a controller) is subsequently designed to decrease the objective function along the solutions of the transformed continuous-time dynamical system. The invariance principle is then applied to analyze the behavior of the solution. Furthermore, the proposed method is modified to address cases in which solutions do not satisfy the equality constraint. Finally, the proposed optimization approach is applied to three examples to demonstrate its effectiveness.

PWR control system design using advanced linear and non-linear methodologies

International Nuclear Information System (INIS)

Rabindran, N.; Whitmarsh-Everiss, M.J.

2004-01-01

Consideration is here given to the methodology deployed for non-linear heuristic analysis in the time domain supported by multi-variable linear control system design methods for the purposes of operational dynamics and control system analysis. This methodology is illustrated by the application of structural singular value μ analysis to Pressurised Water Reactor control system design. (author)

Dynamic linearization system for a radiation gauge

International Nuclear Information System (INIS)

Panarello, J.A.

1977-01-01

The linearization system and process converts a high resolution non-linear analog input signal, representative of the thickness of an object, into a high resolution linear analog output signal suitable for use in driving a variety of output devices. The system requires only a small amount of memory for storing pre-calculated non-linear correction coefficients. The system channels the input signal to separate circuit paths so that it may be used directly to; locate an appropriate correction coefficient; develop a correction term after an appropriate correction coefficient is located; and develop a linearized signal having the same high resolution inherent in the input signal. The system processes the linearized signal to compensate for the possible errors introduced by radiation source noise. The processed linearized signal is the high resolution linear analog output signal which accurately represents the thickness of the object being gauged

Coherent Structures and Entropy in Constrained, Modulationally Unstable, Nonintegrable Systems

International Nuclear Information System (INIS)

Rumpf, Benno; Newell, Alan C.

2001-01-01

Many studies have shown that nonintegrable systems with modulational instabilities constrained by more than one conservation law exhibit universal long time behavior involving large coherent structures in a sea of small fluctuations. We show how this behavior can be explained in detail by simple thermodynamic arguments

A constrained approach to multiscale stochastic simulation of chemically reacting systems

KAUST Repository

Cotter, Simon L.

2011-01-01

Stochastic simulation of coupled chemical reactions is often computationally intensive, especially if a chemical system contains reactions occurring on different time scales. In this paper, we introduce a multiscale methodology suitable to address this problem, assuming that the evolution of the slow species in the system is well approximated by a Langevin process. It is based on the conditional stochastic simulation algorithm (CSSA) which samples from the conditional distribution of the suitably defined fast variables, given values for the slow variables. In the constrained multiscale algorithm (CMA) a single realization of the CSSA is then used for each value of the slow variable to approximate the effective drift and diffusion terms, in a similar manner to the constrained mean-force computations in other applications such as molecular dynamics. We then show how using the ensuing Fokker-Planck equation approximation, we can in turn approximate average switching times in stochastic chemical systems. © 2011 American Institute of Physics.

Reliability evaluation of linear multi-state consecutively-connected systems constrained by m consecutive and n total gaps

International Nuclear Information System (INIS)

Yu, Huan; Yang, Jun; Peng, Rui; Zhao, Yu

2016-01-01

This paper extends the linear multi-state consecutively-connected system (LMCCS) to the case of LMCCS-MN, where MN denotes the dual constraints of m consecutive gaps and n total gaps. All the nodes are distributed along a line and form a sequence. The distances between the adjacent nodes are usually non-uniform. The nodes except the last one can contain statistically independent multi-state connection elements (MCEs). Each MCE can provide a connection between the node at which it is located and the next nodes along the sequence. The LMCCS-MN fails if it meets either of the two constraints. The universal generating function technique is adopted to evaluate the system reliability. The optimal allocations of LMCCS-MN with two different types of failures are solved by genetic algorithm. Finally, two examples are given for the demonstration of the proposed model. - Highlights: • A new model of multi-state consecutively-connected system (LMCCS-MN) is proposed. • The non-uniform distributed nodes are involved in the proposed LMCCS-MN model. • An algorithm for system reliability evaluation is provided by the UGF method. • The computational complexity of the proposed algorithm is discussed in detail. • Optimal element allocation problem is formulated and solved.

Cosmicflows Constrained Local UniversE Simulations

Science.gov (United States)

Sorce, Jenny G.; Gottlöber, Stefan; Yepes, Gustavo; Hoffman, Yehuda; Courtois, Helene M.; Steinmetz, Matthias; Tully, R. Brent; Pomarède, Daniel; Carlesi, Edoardo

2016-01-01

This paper combines observational data sets and cosmological simulations to generate realistic numerical replicas of the nearby Universe. The latter are excellent laboratories for studies of the non-linear process of structure formation in our neighbourhood. With measurements of radial peculiar velocities in the local Universe (cosmicflows-2) and a newly developed technique, we produce Constrained Local UniversE Simulations (CLUES). To assess the quality of these constrained simulations, we compare them with random simulations as well as with local observations. The cosmic variance, defined as the mean one-sigma scatter of cell-to-cell comparison between two fields, is significantly smaller for the constrained simulations than for the random simulations. Within the inner part of the box where most of the constraints are, the scatter is smaller by a factor of 2 to 3 on a 5 h-1 Mpc scale with respect to that found for random simulations. This one-sigma scatter obtained when comparing the simulated and the observation-reconstructed velocity fields is only 104 ± 4 km s-1, I.e. the linear theory threshold. These two results demonstrate that these simulations are in agreement with each other and with the observations of our neighbourhood. For the first time, simulations constrained with observational radial peculiar velocities resemble the local Universe up to a distance of 150 h-1 Mpc on a scale of a few tens of megaparsecs. When focusing on the inner part of the box, the resemblance with our cosmic neighbourhood extends to a few megaparsecs (<5 h-1 Mpc). The simulations provide a proper large-scale environment for studies of the formation of nearby objects.

Constrained non-linear optimization in 3D reflexion tomography; Problemes d'optimisation non-lineaire avec contraintes en tomographie de reflexion 3D

Energy Technology Data Exchange (ETDEWEB)

Delbos, F

2004-11-01

Reflexion tomography allows the determination of a subsurface velocity model from the travel times of seismic waves. The introduction of a priori information in this inverse problem can lead to the resolution of a constrained non-linear least-squares problem. The goal of the thesis is to improve the resolution techniques of this optimization problem, whose main difficulties are its ill-conditioning, its large scale and an expensive cost function in terms of CPU time. Thanks to a detailed study of the problem and to numerous numerical experiments, we justify the use of a sequential quadratic programming method, in which the tangential quadratic programs are solved by an original augmented Lagrangian method. We show the global linear convergence of the latter. The efficiency and robustness of the approach are demonstrated on several synthetic examples and on two real data cases. (author)

Indirect synthesis of multi-degree of freedom transient systems. [linear programming for a kinematically linear system

Science.gov (United States)

Pilkey, W. D.; Chen, Y. H.

1974-01-01

An indirect synthesis method is used in the efficient optimal design of multi-degree of freedom, multi-design element, nonlinear, transient systems. A limiting performance analysis which requires linear programming for a kinematically linear system is presented. The system is selected using system identification methods such that the designed system responds as closely as possible to the limiting performance. The efficiency is a result of the method avoiding the repetitive systems analyses accompanying other numerical optimization methods.

A Homogeneous and Self-Dual Interior-Point Linear Programming Algorithm for Economic Model Predictive Control

DEFF Research Database (Denmark)

Sokoler, Leo Emil; Frison, Gianluca; Skajaa, Anders

2015-01-01

We develop an efficient homogeneous and self-dual interior-point method (IPM) for the linear programs arising in economic model predictive control of constrained linear systems with linear objective functions. The algorithm is based on a Riccati iteration procedure, which is adapted to the linear...... system of equations solved in homogeneous and self-dual IPMs. Fast convergence is further achieved using a warm-start strategy. We implement the algorithm in MATLAB and C. Its performance is tested using a conceptual power management case study. Closed loop simulations show that 1) the proposed algorithm...

The bounds of feasible space on constrained nonconvex quadratic programming

Science.gov (United States)

Zhu, Jinghao

2008-03-01

This paper presents a method to estimate the bounds of the radius of the feasible space for a class of constrained nonconvex quadratic programmingsE Results show that one may compute a bound of the radius of the feasible space by a linear programming which is known to be a P-problem [N. Karmarkar, A new polynomial-time algorithm for linear programming, Combinatorica 4 (1984) 373-395]. It is proposed that one applies this method for using the canonical dual transformation [D.Y. Gao, Canonical duality theory and solutions to constrained nonconvex quadratic programming, J. Global Optimization 29 (2004) 377-399] for solving a standard quadratic programming problem.

Dynamical systems and linear algebra

OpenAIRE

Colonius, Fritz (Prof.)

2007-01-01

Dynamical systems and linear algebra / F. Colonius, W. Kliemann. - In: Handbook of linear algebra / ed. by Leslie Hogben. - Boca Raton : Chapman & Hall/CRC, 2007. - S. 56,1-56,22. - (Discrete mathematics and its applications)

A constrained regularization method for inverting data represented by linear algebraic or integral equations

Science.gov (United States)

Provencher, Stephen W.

1982-09-01

CONTIN is a portable Fortran IV package for inverting noisy linear operator equations. These problems occur in the analysis of data from a wide variety experiments. They are generally ill-posed problems, which means that errors in an unregularized inversion are unbounded. Instead, CONTIN seeks the optimal solution by incorporating parsimony and any statistical prior knowledge into the regularizor and absolute prior knowledge into equallity and inequality constraints. This can be greatly increase the resolution and accuracyh of the solution. CONTIN is very flexible, consisting of a core of about 50 subprograms plus 13 small "USER" subprograms, which the user can easily modify to specify special-purpose constraints, regularizors, operator equations, simulations, statistical weighting, etc. Specjial collections of USER subprograms are available for photon correlation spectroscopy, multicomponent spectra, and Fourier-Bessel, Fourier and Laplace transforms. Numerically stable algorithms are used throughout CONTIN. A fairly precise definition of information content in terms of degrees of freedom is given. The regularization parameter can be automatically chosen on the basis of an F-test and confidence region. The interpretation of the latter and of error estimates based on the covariance matrix of the constrained regularized solution are discussed. The strategies, methods and options in CONTIN are outlined. The program itself is described in the following paper.

Final focus systems for linear colliders

International Nuclear Information System (INIS)

Erickson, R.A.

1987-11-01

The final focus system of a linear collider must perform two primary functions, it must focus the two opposing beams so that their transverse dimensions at the interaction point are small enough to yield acceptable luminosity, and it must steer the beams together to maintain collisions. In addition, the final focus system must transport the outgoing beams to a location where they can be recycled or safely dumped. Elementary optical considerations for linear collider final focus systems are discussed, followed by chromatic aberrations. The design of the final focus system of the SLAC Linear Collider (SLC) is described. Tuning and diagnostics and steering to collision are discussed. Most of the examples illustrating the concepts covered are drawn from the SLC, but the principles and conclusions are said to be generally applicable to other linear collider designs as well. 26 refs., 17 figs

Uniqueness conditions for constrained three-way factor decompositions with linearly dependent loadings

NARCIS (Netherlands)

Stegeman, Alwin; De Almeida, Andre L. F.

2009-01-01

In this paper, we derive uniqueness conditions for a constrained version of the parallel factor (Parafac) decomposition, also known as canonical decomposition (Candecomp). Candecomp/Parafac (CP) decomposes a three-way array into a prespecified number of outer product arrays. The constraint is that

Linear operator inequalities for strongly stable weakly regular linear systems

NARCIS (Netherlands)

Curtain, RF

2001-01-01

We consider the question of the existence of solutions to certain linear operator inequalities (Lur'e equations) for strongly stable, weakly regular linear systems with generating operators A, B, C, 0. These operator inequalities are related to the spectral factorization of an associated Popov

Reduction of Linear Functional Systems using Fuhrmann's Equivalence

Directory of Open Access Journals (Sweden)

Mohamed S. Boudellioua

2016-11-01

Full Text Available Functional systems arise in the treatment of systems of partial differential equations, delay-differential equations, multidimensional equations, etc. The problem of reducing a linear functional system to a system containing fewer equations and unknowns was first studied by Serre. Finding an equivalent presentation of a linear functional system containing fewer equations and fewer unknowns can generally simplify both the study of the structural properties of the linear functional system and of different numerical analysis issues, and it can sometimes help in solving the linear functional system. In this paper, Fuhrmann's equivalence is used to present a constructive result on the reduction of under-determined linear functional systems to a single equation involving a single unknown. This equivalence transformation has been studied by a number of authors and has been shown to play an important role in the theory of linear functional systems.

Adaptive Fuzzy Output-Constrained Fault-Tolerant Control of Nonlinear Stochastic Large-Scale Systems With Actuator Faults.

Science.gov (United States)

Li, Yongming; Ma, Zhiyao; Tong, Shaocheng

2017-09-01

The problem of adaptive fuzzy output-constrained tracking fault-tolerant control (FTC) is investigated for the large-scale stochastic nonlinear systems of pure-feedback form. The nonlinear systems considered in this paper possess the unstructured uncertainties, unknown interconnected terms and unknown nonaffine nonlinear faults. The fuzzy logic systems are employed to identify the unknown lumped nonlinear functions so that the problems of structured uncertainties can be solved. An adaptive fuzzy state observer is designed to solve the nonmeasurable state problem. By combining the barrier Lyapunov function theory, adaptive decentralized and stochastic control principles, a novel fuzzy adaptive output-constrained FTC approach is constructed. All the signals in the closed-loop system are proved to be bounded in probability and the system outputs are constrained in a given compact set. Finally, the applicability of the proposed controller is well carried out by a simulation example.

Assessment of oscillatory stability constrained available transfer capability

International Nuclear Information System (INIS)

Jain, T.; Singh, S.N.; Srivastava, S.C.

2009-01-01

This paper utilizes a bifurcation approach to compute oscillatory stability constrained available transfer capability (ATC) in an electricity market having bilateral as well as multilateral transactions. Oscillatory instability in non-linear systems can be related to Hopf bifurcation. At the Hopf bifurcation, one pair of the critical eigenvalues of the system Jacobian reaches imaginary axis. A new optimization formulation, including Hopf bifurcation conditions, has been developed in this paper to obtain the dynamic ATC. An oscillatory stability based contingency screening index, which takes into account the impact of transactions on severity of contingency, has been utilized to identify critical contingencies to be considered in determining ATC. The proposed method has been applied for dynamic ATC determination on a 39-bus New England system and a practical 75-bus Indian system considering composite static load as well as dynamic load models. (author)

Introduction to linear systems of differential equations

CERN Document Server

Adrianova, L Ya

1995-01-01

The theory of linear systems of differential equations is one of the cornerstones of the whole theory of differential equations. At its root is the concept of the Lyapunov characteristic exponent. In this book, Adrianova presents introductory material and further detailed discussions of Lyapunov exponents. She also discusses the structure of the space of solutions of linear systems. Classes of linear systems examined are from the narrowest to widest: 1)�autonomous, 2)�periodic, 3)�reducible to autonomous, 4)�nearly reducible to autonomous, 5)�regular. In addition, Adrianova considers the following: stability of linear systems and the influence of perturbations of the coefficients on the stability the criteria of uniform stability and of uniform asymptotic stability in terms of properties of the solutions several estimates of the growth rate of solutions of a linear system in terms of its coefficients How perturbations of the coefficients change all the elements of the spectrum of the system is defin...

Window observers for linear systems

Directory of Open Access Journals (Sweden)

Utkin Vadim

2000-01-01

Full Text Available Given a linear system x ˙ = A x + B u with output y = C x and a window function ω ( t , i.e., ∀ t , ω ( t ∈ {0,1 }, and assuming that the window function is Lebesgue measurable, we refer to the following observer, x ˆ = A x + B u + ω ( t L C ( x − x ˆ as a window observer. The stability issue is treated in this paper. It is proven that for linear time-invariant systems, the window observer can be stabilized by an appropriate design under a very mild condition on the window functions, albeit for linear time-varying system, some regularity of the window functions is required to achieve observer designs with the asymptotic stability. The corresponding design methods are developed. An example is included to illustrate the possible applications

Balanced truncation for linear switched systems

DEFF Research Database (Denmark)

Petreczky, Mihaly; Wisniewski, Rafal; Leth, John-Josef

2013-01-01

In this paper, we present a theoretical analysis of the model reduction algorithm for linear switched systems from Shaker and Wisniewski (2011, 2009) and . This algorithm is a reminiscence of the balanced truncation method for linear parameter varying systems (Wood et al., 1996) [3]. Specifically...

System Entropy Measurement of Stochastic Partial Differential Systems

Directory of Open Access Journals (Sweden)

Bor-Sen Chen

2016-03-01

Full Text Available System entropy describes the dispersal of a system’s energy and is an indication of the disorder of a physical system. Several system entropy measurement methods have been developed for dynamic systems. However, most real physical systems are always modeled using stochastic partial differential dynamic equations in the spatio-temporal domain. No efficient method currently exists that can calculate the system entropy of stochastic partial differential systems (SPDSs in consideration of the effects of intrinsic random fluctuation and compartment diffusion. In this study, a novel indirect measurement method is proposed for calculating of system entropy of SPDSs using a Hamilton–Jacobi integral inequality (HJII-constrained optimization method. In other words, we solve a nonlinear HJII-constrained optimization problem for measuring the system entropy of nonlinear stochastic partial differential systems (NSPDSs. To simplify the system entropy measurement of NSPDSs, the global linearization technique and finite difference scheme were employed to approximate the nonlinear stochastic spatial state space system. This allows the nonlinear HJII-constrained optimization problem for the system entropy measurement to be transformed to an equivalent linear matrix inequalities (LMIs-constrained optimization problem, which can be easily solved using the MATLAB LMI-toolbox (MATLAB R2014a, version 8.3. Finally, several examples are presented to illustrate the system entropy measurement of SPDSs.

Constrained non-linear optimization in 3D reflexion tomography; Problemes d'optimisation non-lineaire avec contraintes en tomographie de reflexion 3D

Energy Technology Data Exchange (ETDEWEB)

Delbos, F.

2004-11-01

Reflexion tomography allows the determination of a subsurface velocity model from the travel times of seismic waves. The introduction of a priori information in this inverse problem can lead to the resolution of a constrained non-linear least-squares problem. The goal of the thesis is to improve the resolution techniques of this optimization problem, whose main difficulties are its ill-conditioning, its large scale and an expensive cost function in terms of CPU time. Thanks to a detailed study of the problem and to numerous numerical experiments, we justify the use of a sequential quadratic programming method, in which the tangential quadratic programs are solved by an original augmented Lagrangian method. We show the global linear convergence of the latter. The efficiency and robustness of the approach are demonstrated on several synthetic examples and on two real data cases. (author)

On pole structure assignment in linear systems

Czech Academy of Sciences Publication Activity Database

Loiseau, J.-J.; Zagalak, Petr

2009-01-01

Roč. 82, č. 7 (2009), s. 1179-1192 ISSN 0020-7179 R&D Projects: GA ČR(CZ) GA102/07/1596 Institutional research plan: CEZ:AV0Z10750506 Keywords : linear systems * linear state feedback * pole structure assignment Subject RIV: BC - Control Systems Theory Impact factor: 1.124, year: 2009 http://library.utia.cas.cz/separaty/2009/AS/zagalak-on pole structure assignment in linear systems.pdf

Linear systems a measurement based approach

CERN Document Server

Bhattacharyya, S P; Mohsenizadeh, D N

2014-01-01

This brief presents recent results obtained on the analysis, synthesis and design of systems described by linear equations. It is well known that linear equations arise in most branches of science and engineering as well as social, biological and economic systems. The novelty of this approach is that no models of the system are assumed to be available, nor are they required. Instead, a few measurements made on the system can be processed strategically to directly extract design values that meet specifications without constructing a model of the system, implicitly or explicitly. These new concepts are illustrated by applying them to linear DC and AC circuits, mechanical, civil and hydraulic systems, signal flow block diagrams and control systems. These applications are preliminary and suggest many open problems. The results presented in this brief are the latest effort in this direction and the authors hope these will lead to attractive alternatives to model-based design of engineering and other systems.

Stochastic risk-averse coordinated scheduling of grid integrated energy storage units in transmission constrained wind-thermal systems within a conditional value-at-risk framework

International Nuclear Information System (INIS)

Hemmati, Reza; Saboori, Hedayat; Saboori, Saeid

2016-01-01

In recent decades, wind power resources have been integrated in the power systems increasingly. Besides confirmed benefits, utilization of large share of this volatile source in power generation portfolio has been faced system operators with new challenges in terms of uncertainty management. It is proved that energy storage systems are capable to handle projected uncertainty concerns. Risk-neutral methods have been proposed in the previous literature to schedule storage units considering wind resources uncertainty. Ignoring risk of the cost distributions with non-desirable properties may result in experiencing high costs in some unfavorable scenarios with high probability. In order to control the risk of the operator decisions, this paper proposes a new risk-constrained two-stage stochastic programming model to make optimal decisions on energy storage and thermal units in a transmission constrained hybrid wind-thermal power system. Risk-aversion procedure is explicitly formulated using the conditional value-at-risk measure, because of possessing distinguished features compared to the other risk measures. The proposed model is a mixed integer linear programming considering transmission network, thermal unit dynamics, and storage devices constraints. The simulations results demonstrate that taking the risk of the problem into account will affect scheduling decisions considerably depend on the level of the risk-aversion. - Highlights: • Risk of the operation decisions is handled by using risk-averse programming. • Conditional value-at-risk is used as risk measure. • Optimal risk level is obtained based on the cost/benefit analysis. • The proposed model is a two-stage stochastic mixed integer linear programming. • The unit commitment is integrated with ESSs and wind power penetration.

Final Focus Systems in Linear Colliders

International Nuclear Information System (INIS)

Raubenheimer, Tor

1998-01-01

In colliding beam facilities, the ''final focus system'' must demagnify the beams to attain the very small spot sizes required at the interaction points. The first final focus system with local chromatic correction was developed for the Stanford Linear Collider where very large demagnifications were desired. This same conceptual design has been adopted by all the future linear collider designs as well as the SuperConducting Supercollider, the Stanford and KEK B-Factories, and the proposed Muon Collider. In this paper, the over-all layout, physics constraints, and optimization techniques relevant to the design of final focus systems for high-energy electron-positron linear colliders are reviewed. Finally, advanced concepts to avoid some of the limitations of these systems are discussed

An inexact fuzzy-chance-constrained air quality management model.

Science.gov (United States)

Xu, Ye; Huang, Guohe; Qin, Xiaosheng

2010-07-01

Regional air pollution is a major concern for almost every country because it not only directly relates to economic development, but also poses significant threats to environment and public health. In this study, an inexact fuzzy-chance-constrained air quality management model (IFAMM) was developed for regional air quality management under uncertainty. IFAMM was formulated through integrating interval linear programming (ILP) within a fuzzy-chance-constrained programming (FCCP) framework and could deal with uncertainties expressed as not only possibilistic distributions but also discrete intervals in air quality management systems. Moreover, the constraints with fuzzy variables could be satisfied at different confidence levels such that various solutions with different risk and cost considerations could be obtained. The developed model was applied to a hypothetical case of regional air quality management. Six abatement technologies and sulfur dioxide (SO2) emission trading under uncertainty were taken into consideration. The results demonstrated that IFAMM could help decision-makers generate cost-effective air quality management patterns, gain in-depth insights into effects of the uncertainties, and analyze tradeoffs between system economy and reliability. The results also implied that the trading scheme could achieve lower total abatement cost than a nontrading one.

Systems of Inhomogeneous Linear Equations

Science.gov (United States)

Scherer, Philipp O. J.

Many problems in physics and especially computational physics involve systems of linear equations which arise e.g. from linearization of a general nonlinear problem or from discretization of differential equations. If the dimension of the system is not too large standard methods like Gaussian elimination or QR decomposition are sufficient. Systems with a tridiagonal matrix are important for cubic spline interpolation and numerical second derivatives. They can be solved very efficiently with a specialized Gaussian elimination method. Practical applications often involve very large dimensions and require iterative methods. Convergence of Jacobi and Gauss-Seidel methods is slow and can be improved by relaxation or over-relaxation. An alternative for large systems is the method of conjugate gradients.

Application of Constrained Linear MPC to a Spray Dryer

DEFF Research Database (Denmark)

Petersen, Lars Norbert; Poulsen, Niels Kjølstad; Niemann, Hans Henrik

2014-01-01

In this paper we develop a linear model predictive control (MPC) algorithm for control of a two stage spray dryer. The states are estimated by a stationary Kalman filter. A non-linear first-principle engineering model is developed to simulate the spray drying process. The model is validated against...... experimental data and able to precisely predict the temperatures, the air humidity and the residual moisture in the dryer. The MPC controls these variables to the target and reject disturbances. Spray drying is a cost-effective method to evaporate water from liquid foods and produces a free flowing powder...

CLFs-based optimization control for a class of constrained visual servoing systems.

Science.gov (United States)

Song, Xiulan; Miaomiao, Fu

2017-03-01

In this paper, we use the control Lyapunov function (CLF) technique to present an optimized visual servo control method for constrained eye-in-hand robot visual servoing systems. With the knowledge of camera intrinsic parameters and depth of target changes, visual servo control laws (i.e. translation speed) with adjustable parameters are derived by image point features and some known CLF of the visual servoing system. The Fibonacci method is employed to online compute the optimal value of those adjustable parameters, which yields an optimized control law to satisfy constraints of the visual servoing system. The Lyapunov's theorem and the properties of CLF are used to establish stability of the constrained visual servoing system in the closed-loop with the optimized control law. One merit of the presented method is that there is no requirement of online calculating the pseudo-inverse of the image Jacobian's matrix and the homography matrix. Simulation and experimental results illustrated the effectiveness of the method proposed here. Copyright © 2016 ISA. Published by Elsevier Ltd. All rights reserved.

Autonomous Navigation with Constrained Consistency for C-Ranger

Directory of Open Access Journals (Sweden)

Shujing Zhang

2014-06-01

Full Text Available Autonomous underwater vehicles (AUVs have become the most widely used tools for undertaking complex exploration tasks in marine environments. Their synthetic ability to carry out localization autonomously and build an environmental map concurrently, in other words, simultaneous localization and mapping (SLAM, are considered to be pivotal requirements for AUVs to have truly autonomous navigation. However, the consistency problem of the SLAM system has been greatly ignored during the past decades. In this paper, a consistency constrained extended Kalman filter (EKF SLAM algorithm, applying the idea of local consistency, is proposed and applied to the autonomous navigation of the C-Ranger AUV, which is developed as our experimental platform. The concept of local consistency (LC is introduced after an explicit theoretical derivation of the EKF-SLAM system. Then, we present a locally consistency-constrained EKF-SLAM design, LC-EKF, in which the landmark estimates used for linearization are fixed at the beginning of each local time period, rather than evaluated at the latest landmark estimates. Finally, our proposed LC-EKF algorithm is experimentally verified, both in simulations and sea trials. The experimental results show that the LC-EKF performs well with regard to consistency, accuracy and computational efficiency.

Low-lying excited states by constrained DFT

Science.gov (United States)

Ramos, Pablo; Pavanello, Michele

2018-04-01

Exploiting the machinery of Constrained Density Functional Theory (CDFT), we propose a variational method for calculating low-lying excited states of molecular systems. We dub this method eXcited CDFT (XCDFT). Excited states are obtained by self-consistently constraining a user-defined population of electrons, Nc, in the virtual space of a reference set of occupied orbitals. By imposing this population to be Nc = 1.0, we computed the first excited state of 15 molecules from a test set. Our results show that XCDFT achieves an accuracy in the predicted excitation energy only slightly worse than linear-response time-dependent DFT (TDDFT), but without incurring into problems of variational collapse typical of the more commonly adopted ΔSCF method. In addition, we selected a few challenging processes to test the limits of applicability of XCDFT. We find that in contrast to TDDFT, XCDFT is capable of reproducing energy surfaces featuring conical intersections (azobenzene and H3) with correct topology and correct overall energetics also away from the intersection. Venturing to condensed-phase systems, XCDFT reproduces the TDDFT solvatochromic shift of benzaldehyde when it is embedded by a cluster of water molecules. Thus, we find XCDFT to be a competitive method among single-reference methods for computations of excited states in terms of time to solution, rate of convergence, and accuracy of the result.

Fast Solvers for Dense Linear Systems

Energy Technology Data Exchange (ETDEWEB)

Kauers, Manuel [Research Institute for Symbolic Computation (RISC), Altenbergerstrasse 69, A4040 Linz (Austria)

2008-10-15

It appears that large scale calculations in particle physics often require to solve systems of linear equations with rational number coefficients exactly. If classical Gaussian elimination is applied to a dense system, the time needed to solve such a system grows exponentially in the size of the system. In this tutorial paper, we present a standard technique from computer algebra that avoids this exponential growth: homomorphic images. Using this technique, big dense linear systems can be solved in a much more reasonable time than using Gaussian elimination over the rationals.

Linear quadratic optimization for positive LTI system

Science.gov (United States)

Muhafzan, Yenti, Syafrida Wirma; Zulakmal

2017-05-01

Nowaday the linear quadratic optimization subject to positive linear time invariant (LTI) system constitute an interesting study considering it can become a mathematical model of variety of real problem whose variables have to nonnegative and trajectories generated by these variables must be nonnegative. In this paper we propose a method to generate an optimal control of linear quadratic optimization subject to positive linear time invariant (LTI) system. A sufficient condition that guarantee the existence of such optimal control is discussed.

State space and input-output linear systems

CERN Document Server

Delchamps, David F

1988-01-01

It is difficult for me to forget the mild sense of betrayal I felt some ten years ago when I discovered, with considerable dismay, that my two favorite books on linear system theory - Desoer's Notes for a Second Course on Linear Systems and Brockett's Finite Dimensional Linear Systems - were both out of print. Since that time, of course, linear system theory has undergone a transformation of the sort which always attends the maturation of a theory whose range of applicability is expanding in a fashion governed by technological developments and by the rate at which such advances become a part of engineering practice. The growth of the field has inspired the publication of some excellent books; the encyclopedic treatises by Kailath and Chen, in particular, come immediately to mind. Nonetheless, I was inspired to write this book primarily by my practical needs as a teacher and researcher in the field. For the past five years, I have taught a one semester first year gradu­ ate level linear system theory course i...

Linear collider systems and costs

International Nuclear Information System (INIS)

Loew, G.A.

1993-05-01

The purpose of this paper is to examine some of the systems and sub-systems involved in so-called ''conventional'' e + e - linear colliders and to study how their design affects the overall cost of these machines. There are presently a total of at least six 500 GeV c. of m. linear collider projects under study in the world. Aside from TESLA (superconducting linac at 1.3 GHz) and CLIC (two-beam accelerator with main linac at 30GHz), the other four proposed e + e - linear colliders can be considered ''conventional'' in that their main linacs use the proven technique of driving room temperature accelerator sections with pulsed klystrons and modulators. The centrally distinguishing feature between these projects is their main linac rf frequency: 3 GHz for the DESY machine, 11.424 GHz for the SLAC and JLC machines, and 14 GHz for the VLEPP machine. The other systems, namely the electron and positron sources, preaccelerators, compressors, damping rings and final foci, are fairly similar from project to project. Probably more than 80% of the cost of these linear colliders will be incurred in the two main linacs facing each other and it is therefore in their design and construction that major savings or extra costs may be found

New approach to solve symmetric fully fuzzy linear systems

Indian Academy of Sciences (India)

concepts of fuzzy set theory and then define a fully fuzzy linear system of equations. .... To represent the above problem as fully fuzzy linear system, we represent x .... Fully fuzzy linear systems can be solved by Linear programming approach, ...

A receding horizon scheme for discrete-time polytopic linear parameter varying systems in networked architectures

International Nuclear Information System (INIS)

Franzè, Giuseppe; Lucia, Walter; Tedesco, Francesco

2014-01-01

This paper proposes a Model Predictive Control (MPC) strategy to address regulation problems for constrained polytopic Linear Parameter Varying (LPV) systems subject to input and state constraints in which both plant measurements and command signals in the loop are sent through communication channels subject to time-varying delays (Networked Control System (NCS)). The results here proposed represent a significant extension to the LPV framework of a recent Receding Horizon Control (RHC) scheme developed for the so-called robust case. By exploiting the parameter availability, the pre-computed sequences of one- step controllable sets inner approximations are less conservative than the robust counterpart. The resulting framework guarantees asymptotic stability and constraints fulfilment regardless of plant uncertainties and time-delay occurrences. Finally, experimental results on a laboratory two-tank test-bed show the effectiveness of the proposed approach

Iterated non-linear model predictive control based on tubes and contractive constraints.

Science.gov (United States)

Murillo, M; Sánchez, G; Giovanini, L

2016-05-01

This paper presents a predictive control algorithm for non-linear systems based on successive linearizations of the non-linear dynamic around a given trajectory. A linear time varying model is obtained and the non-convex constrained optimization problem is transformed into a sequence of locally convex ones. The robustness of the proposed algorithm is addressed adding a convex contractive constraint. To account for linearization errors and to obtain more accurate results an inner iteration loop is added to the algorithm. A simple methodology to obtain an outer bounding-tube for state trajectories is also presented. The convergence of the iterative process and the stability of the closed-loop system are analyzed. The simulation results show the effectiveness of the proposed algorithm in controlling a quadcopter type unmanned aerial vehicle. Copyright © 2016 ISA. Published by Elsevier Ltd. All rights reserved.

STABILITY OF LINEAR SYSTEMS WITH MARKOVIAN JUMPS

Directory of Open Access Journals (Sweden)

Jorge Enrique Mayta Guillermo

2016-12-01

Full Text Available In this work we will analyze the stability of linear systems governed by a Markov chain, this family is known in the specialized literature as linear systems with Markov jumps or by its acronyms in English MJLS as it is denoted in [1]. Linear systems governed by a Markov chain are dynamic systems with abrupt changes. We give some denitions of stability for the MJLS system, where these types of stability are equivalent as long as the state space of the Markov chain is nite. Finally we present a theorem that characterizes the stochastic stability by means of an equation of the Lyapunov type. The result is a generalization of a theorem in classical theory.

Using Real and Simulated TNOs to Constrain the Outer Solar System

Science.gov (United States)

Kaib, Nathan

2018-04-01

Over the past 2-3 decades our understanding of the outer solar system’s history and current state has evolved dramatically. An explosion in the number of detected trans-Neptunian objects (TNOs) coupled with simultaneous advances in numerical models of orbital dynamics has driven this rapid evolution. However, successfully constraining the orbital architecture and evolution of the outer solar system requires accurately comparing simulation results with observational datasets. This process is challenging because observed datasets are influenced by orbital discovery biases as well as TNO size and albedo distributions. Meanwhile, such influences are generally absent from numerical results. Here I will review recent work I and others have undertaken using numerical simulations in concert with catalogs of observed TNOs to constrain the outer solar system’s current orbital architecture and past evolution.

Solving discretely-constrained MPEC problems with applications in electric power markets

International Nuclear Information System (INIS)

Gabriel, Steven A.; Leuthold, Florian U.

2010-01-01

Many of the European energy markets are characterized by dominant players that own a large share of their respective countries' generation capacities. In addition to that, there is a significant lack of cross-border transmission capacity. Combining both facts justifies the assumption that these dominant players are able to influence the market outcome of an internal European energy market due to strategic behavior. In this paper, we present a mathematical formulation in order to solve a Stackelberg game for a network-constrained energy market using integer programming. The strategic player is the Stackelberg leader and the independent system operator (including the decisions of the competitive fringe firms) acts as follower. We assume that there is one strategic player which results in a mathematical program with equilibrium constraints (MPEC). This MPEC is reformulated as mixed-integer linear program (MILP) by using disjunctive constraints and linearization. The MILP formulation gives the opportunity to solve the problems reliably and paves the way to add discrete constraints to the original MPEC formulation which can be used in order to solve discretely-constrained mathematical programs with equilibrium constraints (DC-MPECs). We report computational results for a small illustrative network as well as a stylized Western European grid with realistic data. (author)

Solving discretely-constrained MPEC problems with applications in electric power markets

Energy Technology Data Exchange (ETDEWEB)

Gabriel, Steven A. [1143 Glenn L. Martin Hall, Department of Civil and Environmental Engineering, University of Maryland, College Park, MD 20742-3021 (United States); Leuthold, Florian U. [Chair of Energy Economics and Public Sector Management, Dresden University of Technology, 01069 Dresden (Germany)

2010-01-15

Many of the European energy markets are characterized by dominant players that own a large share of their respective countries' generation capacities. In addition to that, there is a significant lack of cross-border transmission capacity. Combining both facts justifies the assumption that these dominant players are able to influence the market outcome of an internal European energy market due to strategic behavior. In this paper, we present a mathematical formulation in order to solve a Stackelberg game for a network-constrained energy market using integer programming. The strategic player is the Stackelberg leader and the independent system operator (including the decisions of the competitive fringe firms) acts as follower. We assume that there is one strategic player which results in a mathematical program with equilibrium constraints (MPEC). This MPEC is reformulated as mixed-integer linear program (MILP) by using disjunctive constraints and linearization. The MILP formulation gives the opportunity to solve the problems reliably and paves the way to add discrete constraints to the original MPEC formulation which can be used in order to solve discretely-constrained mathematical programs with equilibrium constraints (DC-MPECs). We report computational results for a small illustrative network as well as a stylized Western European grid with realistic data. (author)

ITMETH, Iterative Routines for Linear System

International Nuclear Information System (INIS)

Greenbaum, A.

1989-01-01

1 - Description of program or function: ITMETH is a collection of iterative routines for solving large, sparse linear systems. 2 - Method of solution: ITMETH solves general linear systems of the form AX=B using a variety of methods: Jacobi iteration; Gauss-Seidel iteration; incomplete LU decomposition or matrix splitting with iterative refinement; diagonal scaling, matrix splitting, or incomplete LU decomposition with the conjugate gradient method for the problem AA'Y=B, X=A'Y; bi-conjugate gradient method with diagonal scaling, matrix splitting, or incomplete LU decomposition; and ortho-min method with diagonal scaling, matrix splitting, or incomplete LU decomposition. ITMETH also solves symmetric positive definite linear systems AX=B using the conjugate gradient method with diagonal scaling or matrix splitting, or the incomplete Cholesky conjugate gradient method

Cross-training workers in dual resource constrained systems with heterogeneous processing times

NARCIS (Netherlands)

Bokhorst, J. A. C.; Gaalman, G. J. C.

2009-01-01

In this paper, we explore the effect of cross-training workers in Dual Resource Constrained (DRC) systems with machines having different mean processing times. By means of queuing and simulation analysis, we show that the detrimental effects of pooling (cross-training) previously found in single

Numerical Methods for Solution of the Extended Linear Quadratic Control Problem

DEFF Research Database (Denmark)

Jørgensen, John Bagterp; Frison, Gianluca; Gade-Nielsen, Nicolai Fog

2012-01-01

In this paper we present the extended linear quadratic control problem, its efficient solution, and a discussion of how it arises in the numerical solution of nonlinear model predictive control problems. The extended linear quadratic control problem is the optimal control problem corresponding...... to the Karush-Kuhn-Tucker system that constitute the majority of computational work in constrained nonlinear and linear model predictive control problems solved by efficient MPC-tailored interior-point and active-set algorithms. We state various methods of solving the extended linear quadratic control problem...... and discuss instances in which it arises. The methods discussed in the paper have been implemented in efficient C code for both CPUs and GPUs for a number of test examples....

Dynamics of unsymmetric piecewise-linear/non-linear systems using finite elements in time

Science.gov (United States)

Wang, Yu

1995-08-01

The dynamic response and stability of a single-degree-of-freedom system with unsymmetric piecewise-linear/non-linear stiffness are analyzed using the finite element method in the time domain. Based on a Hamilton's weak principle, this method provides a simple and efficient approach for predicting all possible fundamental and sub-periodic responses. The stability of the steady state response is determined by using Floquet's theory without any special effort for calculating transition matrices. This method is applied to a number of examples, demonstrating its effectiveness even for a strongly non-linear problem involving both clearance and continuous stiffness non-linearities. Close agreement is found between available published findings and the predictions of the finite element in time approach, which appears to be an efficient and reliable alternative technique for non-linear dynamic response and stability analysis of periodic systems.

Isolators Including Main Spring Linear Guide Systems

Science.gov (United States)

Goold, Ryan (Inventor); Buchele, Paul (Inventor); Hindle, Timothy (Inventor); Ruebsamen, Dale Thomas (Inventor)

2017-01-01

Embodiments of isolators, such as three parameter isolators, including a main spring linear guide system are provided. In one embodiment, the isolator includes first and second opposing end portions, a main spring mechanically coupled between the first and second end portions, and a linear guide system extending from the first end portion, across the main spring, and toward the second end portion. The linear guide system expands and contracts in conjunction with deflection of the main spring along the working axis, while restricting displacement and rotation of the main spring along first and second axes orthogonal to the working axis.

Displacement measurement system for linear array detector

International Nuclear Information System (INIS)

Zhang Pengchong; Chen Ziyu; Shen Ji

2011-01-01

It presents a set of linear displacement measurement system based on encoder. The system includes displacement encoders, optical lens and read out circuit. Displacement read out unit includes linear CCD and its drive circuit, two amplifier circuits, second order Butterworth low-pass filter and the binarization circuit. The coding way is introduced, and various parts of the experimental signal waveforms are given, and finally a linear experimental test results are given. The experimental results are satisfactory. (authors)

Generalized Cross-Gramian for Linear Systems

DEFF Research Database (Denmark)

Shaker, Hamid Reza

2012-01-01

The cross-gramian is a well-known matrix with embedded controllability and observability information. The cross-gramian is related to the Hankel operator and the Hankel singular values of a linear square system and it has several interesting properties. These properties make the cross...... square symmetric systems, the ordinary cross-gramian does not exist. To cope with this problem, a new generalized cross-gramian is introduced in this paper. In contrast to the ordinary cross-gramian, the generalized cross-gramian can be easily obtained for general linear systems and therefore can be used...

NP-Hardness of optimizing the sum of Rational Linear Functions over an Asymptotic-Linear-Program

OpenAIRE

Chermakani, Deepak Ponvel

2012-01-01

We convert, within polynomial-time and sequential processing, an NP-Complete Problem into a real-variable problem of minimizing a sum of Rational Linear Functions constrained by an Asymptotic-Linear-Program. The coefficients and constants in the real-variable problem are 0, 1, -1, K, or -K, where K is the time parameter that tends to positive infinity. The number of variables, constraints, and rational linear functions in the objective, of the real-variable problem is bounded by a polynomial ...

On Optimal Feedback Control for Stationary Linear Systems

International Nuclear Information System (INIS)

Russell, David L.

2010-01-01

We study linear-quadratic optimal control problems for finite dimensional stationary linear systems AX+BU=Z with output Y=CX+DU from the viewpoint of linear feedback solution. We interpret solutions in relation to system robustness with respect to disturbances Z and relate them to nonlinear matrix equations of Riccati type and eigenvalue-eigenvector problems for the corresponding Hamiltonian system. Examples are included along with an indication of extensions to continuous, i.e., infinite dimensional, systems, primarily of elliptic type.

Breaking down the barriers of using strong authentication and encryption in resource constrained embedded systems

Science.gov (United States)

Knobler, Ron; Scheffel, Peter; Jackson, Scott; Gaj, Kris; Kaps, Jens Peter

2013-05-01

Various embedded systems, such as unattended ground sensors (UGS), are deployed in dangerous areas, where they are subject to compromise. Since numerous systems contain a network of devices that communicate with each other (often times with commercial off the shelf [COTS] radios), an adversary is able to intercept messages between system devices, which jeopardizes sensitive information transmitted by the system (e.g. location of system devices). Secret key algorithms such as AES are a very common means to encrypt all system messages to a sufficient security level, for which lightweight implementations exist for even very resource constrained devices. However, all system devices must use the appropriate key to encrypt and decrypt messages from each other. While traditional public key algorithms (PKAs), such as RSA and Elliptic Curve Cryptography (ECC), provide a sufficiently secure means to provide authentication and a means to exchange keys, these traditional PKAs are not suitable for very resource constrained embedded systems or systems which contain low reliability communication links (e.g. mesh networks), especially as the size of the network increases. Therefore, most UGS and other embedded systems resort to pre-placed keys (PPKs) or other naïve schemes which greatly reduce the security and effectiveness of the overall cryptographic approach. McQ has teamed with the Cryptographic Engineering Research Group (CERG) at George Mason University (GMU) to develop an approach using revolutionary cryptographic techniques that provides both authentication and encryption, but on resource constrained embedded devices, without the burden of large amounts of key distribution or storage.

Convergence Guaranteed Nonlinear Constraint Model Predictive Control via I/O Linearization

Directory of Open Access Journals (Sweden)

Xiaobing Kong

2013-01-01

Full Text Available Constituting reliable optimal solution is a key issue for the nonlinear constrained model predictive control. Input-output feedback linearization is a popular method in nonlinear control. By using an input-output feedback linearizing controller, the original linear input constraints will change to nonlinear constraints and sometimes the constraints are state dependent. This paper presents an iterative quadratic program (IQP routine on the continuous-time system. To guarantee its convergence, another iterative approach is incorporated. The proposed algorithm can reach a feasible solution over the entire prediction horizon. Simulation results on both a numerical example and the continuous stirred tank reactors (CSTR demonstrate the effectiveness of the proposed method.

Perfect commuting-operator strategies for linear system games

Science.gov (United States)

Cleve, Richard; Liu, Li; Slofstra, William

2017-01-01

Linear system games are a generalization of Mermin's magic square game introduced by Cleve and Mittal. They show that perfect strategies for linear system games in the tensor-product model of entanglement correspond to finite-dimensional operator solutions of a certain set of non-commutative equations. We investigate linear system games in the commuting-operator model of entanglement, where Alice and Bob's measurement operators act on a joint Hilbert space, and Alice's operators must commute with Bob's operators. We show that perfect strategies in this model correspond to possibly infinite-dimensional operator solutions of the non-commutative equations. The proof is based around a finitely presented group associated with the linear system which arises from the non-commutative equations.

Constrained systems described by Nambu mechanics

International Nuclear Information System (INIS)

Lassig, C.C.; Joshi, G.C.

1996-01-01

Using the framework of Nambu's generalised mechanics, we obtain a new description of constrained Hamiltonian dynamics, involving the introduction of another degree of freedom in phase space, and the necessity of defining the action integral on a world sheet. We also discuss the problem of quantizing Nambu mechanics. (authors). 5 refs

Adaptive Multi-Agent Systems for Constrained Optimization

Science.gov (United States)

Macready, William; Bieniawski, Stefan; Wolpert, David H.

2004-01-01

Product Distribution (PD) theory is a new framework for analyzing and controlling distributed systems. Here we demonstrate its use for distributed stochastic optimization. First we review one motivation of PD theory, as the information-theoretic extension of conventional full-rationality game theory to the case of bounded rational agents. In this extension the equilibrium of the game is the optimizer of a Lagrangian of the (probability distribution of) the joint state of the agents. When the game in question is a team game with constraints, that equilibrium optimizes the expected value of the team game utility, subject to those constraints. The updating of the Lagrange parameters in the Lagrangian can be viewed as a form of automated annealing, that focuses the MAS more and more on the optimal pure strategy. This provides a simple way to map the solution of any constrained optimization problem onto the equilibrium of a Multi-Agent System (MAS). We present computer experiments involving both the Queen s problem and K-SAT validating the predictions of PD theory and its use for off-the-shelf distributed adaptive optimization.

Control system analysis for the perturbed linear accelerator rf system

CERN Document Server

Sung Il Kwon

2002-01-01

This paper addresses the modeling problem of the linear accelerator RF system in SNS. Klystrons are modeled as linear parameter varying systems. The effect of the high voltage power supply ripple on the klystron output voltage and the output phase is modeled as an additive disturbance. The cavity is modeled as a linear system and the beam current is modeled as the exogenous disturbance. The output uncertainty of the low level RF system which results from the uncertainties in the RF components and cabling is modeled as multiplicative uncertainty. Also, the feedback loop uncertainty and digital signal processing signal conditioning subsystem uncertainties are lumped together and are modeled as multiplicative uncertainty. Finally, the time delays in the loop are modeled as a lumped time delay. For the perturbed open loop system, the closed loop system performance, and stability are analyzed with the PI feedback controller.

CONTROL SYSTEM ANALYSIS FOR THE PERTURBED LINEAR ACCELERATOR RF SYSTEM

International Nuclear Information System (INIS)

SUNG-IL KWON; AMY H. REGAN

2002-01-01

This paper addresses the modeling problem of the linear accelerator RF system in SNS. Klystrons are modeled as linear parameter varying systems. The effect of the high voltage power supply ripple on the klystron output voltage and the output phase is modeled as an additive disturbance. The cavity is modeled as a linear system and the beam current is modeled as the exogenous disturbance. The output uncertainty of the low level RF system which results from the uncertainties in the RF components and cabling is modeled as multiplicative uncertainty. Also, the feedback loop uncertainty and digital signal processing signal conditioning subsystem uncertainties are lumped together and are modeled as multiplicative uncertainty. Finally, the time delays in the loop are modeled as a lumped time delay. For the perturbed open loop system, the closed loop system performance, and stability are analyzed with the PI feedback controller

Sensitivity theory for general non-linear algebraic equations with constraints

International Nuclear Information System (INIS)

Oblow, E.M.

1977-04-01

Sensitivity theory has been developed to a high state of sophistication for applications involving solutions of the linear Boltzmann equation or approximations to it. The success of this theory in the field of radiation transport has prompted study of possible extensions of the method to more general systems of non-linear equations. Initial work in the U.S. and in Europe on the reactor fuel cycle shows that the sensitivity methodology works equally well for those non-linear problems studied to date. The general non-linear theory for algebraic equations is summarized and applied to a class of problems whose solutions are characterized by constrained extrema. Such equations form the basis of much work on energy systems modelling and the econometrics of power production and distribution. It is valuable to have a sensitivity theory available for these problem areas since it is difficult to repeatedly solve complex non-linear equations to find out the effects of alternative input assumptions or the uncertainties associated with predictions of system behavior. The sensitivity theory for a linear system of algebraic equations with constraints which can be solved using linear programming techniques is discussed. The role of the constraints in simplifying the problem so that sensitivity methodology can be applied is highlighted. The general non-linear method is summarized and applied to a non-linear programming problem in particular. Conclusions are drawn in about the applicability of the method for practical problems

Linear and non-linear Modified Gravity forecasts with future surveys

Science.gov (United States)

Casas, Santiago; Kunz, Martin; Martinelli, Matteo; Pettorino, Valeria

2017-12-01

Modified Gravity theories generally affect the Poisson equation and the gravitational slip in an observable way, that can be parameterized by two generic functions (η and μ) of time and space. We bin their time dependence in redshift and present forecasts on each bin for future surveys like Euclid. We consider both Galaxy Clustering and Weak Lensing surveys, showing the impact of the non-linear regime, with two different semi-analytical approximations. In addition to these future observables, we use a prior covariance matrix derived from the Planck observations of the Cosmic Microwave Background. In this work we neglect the information from the cross correlation of these observables, and treat them as independent. Our results show that η and μ in different redshift bins are significantly correlated, but including non-linear scales reduces or even eliminates the correlation, breaking the degeneracy between Modified Gravity parameters and the overall amplitude of the matter power spectrum. We further apply a Zero-phase Component Analysis and identify which combinations of the Modified Gravity parameter amplitudes, in different redshift bins, are best constrained by future surveys. We extend the analysis to two particular parameterizations of μ and η and consider, in addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future surveys will be able to constrain the current values of η and μ at the 2-5% level when using only linear scales (wavevector k < 0 . 15 h/Mpc), depending on the specific time parameterization; sensitivity improves to about 1% when non-linearities are included.

Closed-form solutions for linear regulator design of mechanical systems including optimal weighting matrix selection

Science.gov (United States)

Hanks, Brantley R.; Skelton, Robert E.

1991-01-01

Vibration in modern structural and mechanical systems can be reduced in amplitude by increasing stiffness, redistributing stiffness and mass, and/or adding damping if design techniques are available to do so. Linear Quadratic Regulator (LQR) theory in modern multivariable control design, attacks the general dissipative elastic system design problem in a global formulation. The optimal design, however, allows electronic connections and phase relations which are not physically practical or possible in passive structural-mechanical devices. The restriction of LQR solutions (to the Algebraic Riccati Equation) to design spaces which can be implemented as passive structural members and/or dampers is addressed. A general closed-form solution to the optimal free-decay control problem is presented which is tailored for structural-mechanical system. The solution includes, as subsets, special cases such as the Rayleigh Dissipation Function and total energy. Weighting matrix selection is a constrained choice among several parameters to obtain desired physical relationships. The closed-form solution is also applicable to active control design for systems where perfect, collocated actuator-sensor pairs exist.

A CLASS OF NONMONOTONE TRUST REGION ALGORITHMS FOR LINEARLY CONSTRAINED OPTIMIZATION%线性约束优化的一类非单调信赖域算法

Institute of Scientific and Technical Information of China (English)

葛恒武; 陈中文

2002-01-01

We present a class of nonmonotone trust region algorithms for linearly constrained optimization in this paper.The algorithm may adjust automatically the scope of the monotonicity by the degree that the quadratic model is "trusted".Under the suitable conditions,it is proved that any limit point of the infinite sequence generated by the algorithm is the Kuhn-Tucker point of the primal problem.Finally,some numerical results show that the new algorithm is very effective.

A constrained Delaunay discretization method for adaptively meshing highly discontinuous geological media

Science.gov (United States)

Wang, Yang; Ma, Guowei; Ren, Feng; Li, Tuo

2017-12-01

A constrained Delaunay discretization method is developed to generate high-quality doubly adaptive meshes of highly discontinuous geological media. Complex features such as three-dimensional discrete fracture networks (DFNs), tunnels, shafts, slopes, boreholes, water curtains, and drainage systems are taken into account in the mesh generation. The constrained Delaunay triangulation method is used to create adaptive triangular elements on planar fractures. Persson's algorithm (Persson, 2005), based on an analogy between triangular elements and spring networks, is enriched to automatically discretize a planar fracture into mesh points with varying density and smooth-quality gradient. The triangulated planar fractures are treated as planar straight-line graphs (PSLGs) to construct piecewise-linear complex (PLC) for constrained Delaunay tetrahedralization. This guarantees the doubly adaptive characteristic of the resulted mesh: the mesh is adaptive not only along fractures but also in space. The quality of elements is compared with the results from an existing method. It is verified that the present method can generate smoother elements and a better distribution of element aspect ratios. Two numerical simulations are implemented to demonstrate that the present method can be applied to various simulations of complex geological media that contain a large number of discontinuities.

Energy balance in a system with quasispherical linear compression

International Nuclear Information System (INIS)

Es'kov, A.G.; Kozlov, N.P.; Kurtmullaev, R.K.; Semenov, V.N.; Khvesyuk, V.I.; Yaminskii, A.V.

1983-01-01

This letter reports the resists of some experimental studies and a numerical simulation of the Tor-linear fusion system, 1 in which a heavy plasma shell with a closed magnetic structure is compressed in a quasispherical manner. The parameters of the Tor-Linear, at the Kurchatov Institute of Atomic Energy in Moscow are as follows: The energy stored in the system which accelerates the linear is E = 0.5 MJ; the linear mass is m = 0.2 kg; the working volume of the linear module is 1.5 x 10 -3 m 3 ; the linear velocity is approx.10 3 m/s; the guiding field in the toriod in the linear is 1--10 x 10 21 m -3 ; and the intial volume of the plasma in the linear chamber is 2.5 x 10 -4 m 3 . In this series of experiments, new solutions were developed for all the systems of the plasma--linear complex of the Tor-Linear: to produce a plasma toroid, to transport it, and to trap it in the linear cavity

Time-optimal feedback control for linear systems

International Nuclear Information System (INIS)

Mirica, S.

1976-01-01

The paper deals with the results of qualitative investigations of the time-optimal feedback control for linear systems with constant coefficients. In the first section, after some definitions and notations, two examples are given and it is shown that even the time-optimal control problem for linear systems with constant coefficients which looked like ''completely solved'' requires a further qualitative investigation of the stability to ''permanent perturbations'' of optimal feedback control. In the second section some basic results of the linear time-optimal control problem are reviewed. The third section deals with the definition of Boltyanskii's ''regular synthesis'' and its connection to Filippov's theory of right-hand side discontinuous differential equations. In the fourth section a theorem is proved concerning the stability to perturbations of time-optimal feedback control for linear systems with scalar control. In the last two sections it is proved that, if the matrix which defines the system has only real eigenvalues or is three-dimensional, the time-optimal feedback control defines a regular synthesis and therefore is stable to perturbations. (author)

Robust media processing on programmable power-constrained systems

Science.gov (United States)

McVeigh, Jeff

2005-03-01

To achieve consumer-level quality, media systems must process continuous streams of audio and video data while maintaining exacting tolerances on sampling rate, jitter, synchronization, and latency. While it is relatively straightforward to design fixed-function hardware implementations to satisfy worst-case conditions, there is a growing trend to utilize programmable multi-tasking solutions for media applications. The flexibility of these systems enables support for multiple current and future media formats, which can reduce design costs and time-to-market. This paper provides practical engineering solutions to achieve robust media processing on such systems, with specific attention given to power-constrained platforms. The techniques covered in this article utilize the fundamental concepts of algorithm and software optimization, software/hardware partitioning, stream buffering, hierarchical prioritization, and system resource and power management. A novel enhancement to dynamically adjust processor voltage and frequency based on buffer fullness to reduce system power consumption is examined in detail. The application of these techniques is provided in a case study of a portable video player implementation based on a general-purpose processor running a non real-time operating system that achieves robust playback of synchronized H.264 video and MP3 audio from local storage and streaming over 802.11.

Useful tools for non-linear systems: Several non-linear integral inequalities

Czech Academy of Sciences Publication Activity Database

Agahi, H.; Mohammadpour, A.; Mesiar, Radko; Vaezpour, M. S.

2013-01-01

Roč. 49, č. 1 (2013), s. 73-80 ISSN 0950-7051 R&D Projects: GA ČR GAP402/11/0378 Institutional support: RVO:67985556 Keywords : Monotone measure * Comonotone functions * Integral inequalities * Universal integral Subject RIV: BA - General Mathematics Impact factor: 3.058, year: 2013 http://library.utia.cas.cz/separaty/2013/E/mesiar-useful tools for non-linear systems several non-linear integral inequalities.pdf

Signals and transforms in linear systems analysis

CERN Document Server

Wasylkiwskyj, Wasyl

2013-01-01

Signals and Transforms in Linear Systems Analysis covers the subject of signals and transforms, particularly in the context of linear systems theory. Chapter 2 provides the theoretical background for the remainder of the text. Chapter 3 treats Fourier series and integrals. Particular attention is paid to convergence properties at step discontinuities. This includes the Gibbs phenomenon and its amelioration via the Fejer summation techniques. Special topics include modulation and analytic signal representation, Fourier transforms and analytic function theory, time-frequency analysis and frequency dispersion. Fundamentals of linear system theory for LTI analogue systems, with a brief account of time-varying systems, are covered in Chapter 4 . Discrete systems are covered in Chapters 6 and 7.  The Laplace transform treatment in Chapter 5 relies heavily on analytic function theory as does Chapter 8 on Z -transforms. The necessary background on complex variables is provided in Appendix A. This book is intended to...

Stability Constrained Efficiency Optimization for Droop Controlled DC-DC Conversion System

DEFF Research Database (Denmark)

Meng, Lexuan; Dragicevic, Tomislav; Guerrero, Josep M.

2013-01-01

implementing tertiary regulation. Moreover, system dynamic is affected when shifting VRs. Therefore, the stability is considered in optimization by constraining the eigenvalues arising from dynamic state space model of the system. Genetic algorithm is used in searching for global efficiency optimum while....... As the efficiency of each converter changes with output power, virtual resistances (VRs) are set as decision variables for adjusting power sharing proportion among converters. It is noteworthy that apart from restoring the voltage deviation, secondary control plays an important role to stabilize dc bus voltage when...

Chaos as an intermittently forced linear system.

Science.gov (United States)

Brunton, Steven L; Brunton, Bingni W; Proctor, Joshua L; Kaiser, Eurika; Kutz, J Nathan

2017-05-30

Understanding the interplay of order and disorder in chaos is a central challenge in modern quantitative science. Approximate linear representations of nonlinear dynamics have long been sought, driving considerable interest in Koopman theory. We present a universal, data-driven decomposition of chaos as an intermittently forced linear system. This work combines delay embedding and Koopman theory to decompose chaotic dynamics into a linear model in the leading delay coordinates with forcing by low-energy delay coordinates; this is called the Hankel alternative view of Koopman (HAVOK) analysis. This analysis is applied to the Lorenz system and real-world examples including Earth's magnetic field reversal and measles outbreaks. In each case, forcing statistics are non-Gaussian, with long tails corresponding to rare intermittent forcing that precedes switching and bursting phenomena. The forcing activity demarcates coherent phase space regions where the dynamics are approximately linear from those that are strongly nonlinear.The huge amount of data generated in fields like neuroscience or finance calls for effective strategies that mine data to reveal underlying dynamics. Here Brunton et al.develop a data-driven technique to analyze chaotic systems and predict their dynamics in terms of a forced linear model.

The theory of a general quantum system interacting with a linear dissipative system

International Nuclear Information System (INIS)

Feynman, R.P.; Vernon, F.L.

2000-01-01

A formalism has been developed, using Feynman's space-time formulation of nonrelativistic quantum mechanics whereby the behavior of a system of interest, which is coupled to other external quantum systems, may be calculated in terms of its own variables only. It is shown that the effect of the external systems in such a formalism can always be included in a general class of functionals (influence functionals) of the coordinates of the system only. The properties of influence functionals for general systems are examined. Then, specific forms of influence functionals representing the effect of definite and random classical forces, linear dissipative systems at finite temperatures, and combinations of these are analyzed in detail. The linear system analysis is first done for perfectly linear systems composed of combinations of harmonic oscillators, loss being introduced by continuous distributions of oscillators. Then approximately linear systems and restrictions necessary for the linear behavior are considered. Influence functionals for all linear systems are shown to have the same form in terms of their classical response functions. In addition, a fluctuation-dissipation theorem is derived relating temperature and dissipation of the linear system to a fluctuating classical potential acting on the system of interest which reduces to the Nyquist-Johnson relation for noise in the case of electric circuits. Sample calculations of transition probabilities for the spontaneous emission of an atom in free space and in a cavity are made. Finally, a theorem is proved showing that within the requirements of linearity all sources of noise or quantum fluctuation introduced by maser-type amplification devices are accounted for by a classical calculation of the characteristics of the maser

Sensitive Constrained Optimal PMU Allocation with Complete Observability for State Estimation Solution

Directory of Open Access Journals (Sweden)

R. Manam

2017-12-01

Full Text Available In this paper, a sensitive constrained integer linear programming approach is formulated for the optimal allocation of Phasor Measurement Units (PMUs in a power system network to obtain state estimation. In this approach, sensitive buses along with zero injection buses (ZIB are considered for optimal allocation of PMUs in the network to generate state estimation solutions. Sensitive buses are evolved from the mean of bus voltages subjected to increase of load consistently up to 50%. Sensitive buses are ranked in order to place PMUs. Sensitive constrained optimal PMU allocation in case of single line and no line contingency are considered in observability analysis to ensure protection and control of power system from abnormal conditions. Modeling of ZIB constraints is included to minimize the number of PMU network allocations. This paper presents optimal allocation of PMU at sensitive buses with zero injection modeling, considering cost criteria and redundancy to increase the accuracy of state estimation solution without losing observability of the whole system. Simulations are carried out on IEEE 14, 30 and 57 bus systems and results obtained are compared with traditional and other state estimation methods available in the literature, to demonstrate the effectiveness of the proposed method.

Asymptotic Likelihood Distribution for Correlated & Constrained Systems

CERN Document Server

Agarwal, Ujjwal

2016-01-01

It describes my work as summer student at CERN. The report discusses the asymptotic distribution of the likelihood ratio for total no. of parameters being h and 2 out of these being are constrained and correlated.

Input-constrained model predictive control via the alternating direction method of multipliers

DEFF Research Database (Denmark)

Sokoler, Leo Emil; Frison, Gianluca; Andersen, Martin S.

2014-01-01

This paper presents an algorithm, based on the alternating direction method of multipliers, for the convex optimal control problem arising in input-constrained model predictive control. We develop an efficient implementation of the algorithm for the extended linear quadratic control problem (LQCP......) with input and input-rate limits. The algorithm alternates between solving an extended LQCP and a highly structured quadratic program. These quadratic programs are solved using a Riccati iteration procedure, and a structure-exploiting interior-point method, respectively. The computational cost per iteration...... is quadratic in the dimensions of the controlled system, and linear in the length of the prediction horizon. Simulations show that the approach proposed in this paper is more than an order of magnitude faster than several state-of-the-art quadratic programming algorithms, and that the difference in computation...

Linear and non-linear energy barriers in systems of interacting single-domain ferromagnetic particles

International Nuclear Information System (INIS)

Petrila, Iulian; Bodale, Ilie; Rotarescu, Cristian; Stancu, Alexandru

2011-01-01

A comparative analysis between linear and non-linear energy barriers used for modeling statistical thermally-excited ferromagnetic systems is presented. The linear energy barrier is obtained by new symmetry considerations about the anisotropy energy and the link with the non-linear energy barrier is also presented. For a relevant analysis we compare the effects of linear and non-linear energy barriers implemented in two different models: Preisach-Neel and Ising-Metropolis. The differences between energy barriers which are reflected in different coercive field dependence of the temperature are also presented. -- Highlights: → The linear energy barrier is obtained from symmetry considerations. → The linear and non-linear energy barriers are calibrated and implemented in Preisach-Neel and Ising-Metropolis models. → The temperature and time effects of the linear and non-linear energy barriers are analyzed.

Krylov Subspace Methods for Complex Non-Hermitian Linear Systems. Thesis

Science.gov (United States)

Freund, Roland W.

1991-01-01

We consider Krylov subspace methods for the solution of large sparse linear systems Ax = b with complex non-Hermitian coefficient matrices. Such linear systems arise in important applications, such as inverse scattering, numerical solution of time-dependent Schrodinger equations, underwater acoustics, eddy current computations, numerical computations in quantum chromodynamics, and numerical conformal mapping. Typically, the resulting coefficient matrices A exhibit special structures, such as complex symmetry, or they are shifted Hermitian matrices. In this paper, we first describe a Krylov subspace approach with iterates defined by a quasi-minimal residual property, the QMR method, for solving general complex non-Hermitian linear systems. Then, we study special Krylov subspace methods designed for the two families of complex symmetric respectively shifted Hermitian linear systems. We also include some results concerning the obvious approach to general complex linear systems by solving equivalent real linear systems for the real and imaginary parts of x. Finally, numerical experiments for linear systems arising from the complex Helmholtz equation are reported.

Normal form of linear systems depending on parameters

International Nuclear Information System (INIS)

Nguyen Huynh Phan.

1995-12-01

In this paper we resolve completely the problem to find normal forms of linear systems depending on parameters for the feedback action that we have studied for the special case of controllable linear systems. (author). 24 refs

Numerical solution of large sparse linear systems

International Nuclear Information System (INIS)

Meurant, Gerard; Golub, Gene.

1982-02-01

This note is based on one of the lectures given at the 1980 CEA-EDF-INRIA Numerical Analysis Summer School whose aim is the study of large sparse linear systems. The main topics are solving least squares problems by orthogonal transformation, fast Poisson solvers and solution of sparse linear system by iterative methods with a special emphasis on preconditioned conjuguate gradient method [fr

Fast Combinatorial Algorithm for the Solution of Linearly Constrained Least Squares Problems

Science.gov (United States)

Van Benthem, Mark H.; Keenan, Michael R.

2008-11-11

A fast combinatorial algorithm can significantly reduce the computational burden when solving general equality and inequality constrained least squares problems with large numbers of observation vectors. The combinatorial algorithm provides a mathematically rigorous solution and operates at great speed by reorganizing the calculations to take advantage of the combinatorial nature of the problems to be solved. The combinatorial algorithm exploits the structure that exists in large-scale problems in order to minimize the number of arithmetic operations required to obtain a solution.

Optimization of an implicit constrained multi-physics system for motor wheels of electric vehicle

International Nuclear Information System (INIS)

Lei, Fei; Du, Bin; Liu, Xin; Xie, Xiaoping; Chai, Tian

2016-01-01

In this paper, implicit constrained multi-physics model of a motor wheel for an electric vehicle is built and then optimized. A novel optimization approach is proposed to solve the compliance problem between implicit constraints and stochastic global optimization. Firstly, multi-physics model of motor wheel is built from the theories of structural mechanics, electromagnetism and thermal physics. Then, implicit constraints are applied from the vehicle performances and magnetic characteristics. Implicit constrained optimization is carried out by a series of unconstrained optimization and verifications. In practice, sequentially updated subspaces are designed to completely substitute the original design space in local areas. In each subspace, a solution is obtained and is then verified by the implicit constraints. Optimal solutions which satisfy the implicit constraints are accepted as final candidates. The final global optimal solution is optimized from those candidates. Discussions are carried out to discover the differences between optimal solutions with unconstrained problem and different implicit constrained problems. Results show that the implicit constraints have significant influences on the optimal solution and the proposed approach is effective in finding the optimals. - Highlights: • An implicit constrained multi-physics model is built for sizing a motor wheel. • Vehicle dynamic performances are applied as implicit constraints for nonlinear system. • An efficient novel optimization is proposed to explore the constrained design space. • The motor wheel is optimized to achieve maximum efficiency on vehicle dynamics. • Influences of implicit constraints on vehicle performances are compared and analyzed.

Linear System of Equations, Matrix Inversion, and Linear Programming Using MS Excel

Science.gov (United States)

El-Gebeily, M.; Yushau, B.

2008-01-01

In this note, we demonstrate with illustrations two different ways that MS Excel can be used to solve Linear Systems of Equation, Linear Programming Problems, and Matrix Inversion Problems. The advantage of using MS Excel is its availability and transparency (the user is responsible for most of the details of how a problem is solved). Further, we…

Lectures on algebraic system theory: Linear systems over rings

Science.gov (United States)

Kamen, E. W.

1978-01-01

The presentation centers on four classes of systems that can be treated as linear systems over a ring. These are: (1) discrete-time systems over a ring of scalars such as the integers; (2) continuous-time systems containing time delays; (3) large-scale discrete-time systems; and (4) time-varying discrete-time systems.

Final focus systems for linear colliders

International Nuclear Information System (INIS)

Helm, R.; Irwin, J.

1992-08-01

Final focus systems for linear colliders present many exacting challenges in beam optics, component design, and beam quality. Efforts to resolve these problems as they relate to a new generation of linear colliders are under way at several laboratories around the world. We will outline criteria for final focus systems and discuss the current state of understanding and resolution of the outstanding problems. We will discuss tolerances on alignment, field quality and stability for optical elements, and the implications for beam parameters such as emittance, energy spread, bunch length, and stability in position and energy. Beam-based correction procedures, which in principle can alleviate many of the tolerances, will be described. Preliminary results from the Final Focus Test Beam (FFTB) under construction at SLAC will be given. Finally, we mention conclusions from operating experience at the Stanford Linear Collider (SLC)

Final focus systems for linear colliders

International Nuclear Information System (INIS)

Helm, R.; Irwing, J.

1992-01-01

Final focus systems for linear colliders present many exacting challenges in beam optics, component design, and beam quality. Efforts to resolve these problems as they relate to a new generation of linear colliders are under way at several laboratories around the world. We outline criteria for final focus systems and discuss the current state of understanding and resolution of the outstanding problems. We discuss tolerances on alignment, field quality and stability for optical elements, and the implications for beam parameters such as emittance, energy spread , bunch length, and stability in position and energy. Beam-based correction procedures, which in principle can alleviate many of the tolerances, are described. Preliminary results from the Final Focus Test Beam (FFTB) under construction at SLAC are given. Finally, we mention conclusions from operating experience at the Stanford Linear Collider (SLC). (Author) 16 refs., 4 tabs., 6 figs

SLFP: a stochastic linear fractional programming approach for sustainable waste management.

Science.gov (United States)

Zhu, H; Huang, G H

2011-12-01

A stochastic linear fractional programming (SLFP) approach is developed for supporting sustainable municipal solid waste management under uncertainty. The SLFP method can solve ratio optimization problems associated with random information, where chance-constrained programming is integrated into a linear fractional programming framework. It has advantages in: (1) comparing objectives of two aspects, (2) reflecting system efficiency, (3) dealing with uncertainty expressed as probability distributions, and (4) providing optimal-ratio solutions under different system-reliability conditions. The method is applied to a case study of waste flow allocation within a municipal solid waste (MSW) management system. The obtained solutions are useful for identifying sustainable MSW management schemes with maximized system efficiency under various constraint-violation risks. The results indicate that SLFP can support in-depth analysis of the interrelationships among system efficiency, system cost and system-failure risk. Copyright © 2011 Elsevier Ltd. All rights reserved.

Constrained dansyl derivatives reveal bacterial specificity of highly conserved thymidylate synthases.

Science.gov (United States)

Calò, Sanuele; Tondi, Donatella; Ferrari, Stefania; Venturelli, Alberto; Ghelli, Stefano; Costi, Maria Paola

2008-03-25

The elucidation of the structural/functional specificities of highly conserved enzymes remains a challenging area of investigation, and enzymes involved in cellular replication are important targets for functional studies and drug discovery. Thymidylate synthase (TS, ThyA) governs the synthesis of thymidylate for use in DNA synthesis. The present study focused on Lactobacillus casei TS (LcTS) and Escherichia coli TS (EcTS), which exhibit 50 % sequence identity and strong folding similarity. We have successfully designed and validated a chemical model in which linear, but not constrained, dansyl derivatives specifically complement the LcTS active site. Conversely, chemically constrained dansyl derivatives showed up to 1000-fold improved affinity for EcTS relative to the inhibitory activity of linear derivatives. This study demonstrates that the accurate design of small ligands can uncover functional features of highly conserved enzymes.

Dynamical computation of constrained flexible systems using a modal Udwadia-Kalaba formulation: Application to musical instruments.

Science.gov (United States)

Antunes, J; Debut, V

2017-02-01

Most musical instruments consist of dynamical subsystems connected at a number of constraining points through which energy flows. For physical sound synthesis, one important difficulty deals with enforcing these coupling constraints. While standard techniques include the use of Lagrange multipliers or penalty methods, in this paper, a different approach is explored, the Udwadia-Kalaba (U-K) formulation, which is rooted on analytical dynamics but avoids the use of Lagrange multipliers. This general and elegant formulation has been nearly exclusively used for conceptual systems of discrete masses or articulated rigid bodies, namely, in robotics. However its natural extension to deal with continuous flexible systems is surprisingly absent from the literature. Here, such a modeling strategy is developed and the potential of combining the U-K equation for constrained systems with the modal description is shown, in particular, to simulate musical instruments. Objectives are twofold: (1) Develop the U-K equation for constrained flexible systems with subsystems modelled through unconstrained modes; and (2) apply this framework to compute string/body coupled dynamics. This example complements previous work [Debut, Antunes, Marques, and Carvalho, Appl. Acoust. 108, 3-18 (2016)] on guitar modeling using penalty methods. Simulations show that the proposed technique provides similar results with a significant improvement in computational efficiency.

Online constrained model-based reinforcement learning

CSIR Research Space (South Africa)

Van Niekerk, B

2017-08-01

Full Text Available Constrained Model-based Reinforcement Learning Benjamin van Niekerk School of Computer Science University of the Witwatersrand South Africa Andreas Damianou∗ Amazon.com Cambridge, UK Benjamin Rosman Council for Scientific and Industrial Research, and School... MULTIPLE SHOOTING Using direct multiple shooting (Bock and Plitt, 1984), problem (1) can be transformed into a structured non- linear program (NLP). First, the time horizon [t0, t0 + T ] is partitioned into N equal subintervals [tk, tk+1] for k = 0...

Linear Actuator System for the NASA Docking System

Science.gov (United States)

Dick, Brandon N.; Oesch, Christopher; Rupp, Timothy W.

2017-01-01

The Linear Actuator System (LAS) is a major sub-system within the NASA Docking System (NDS). The NDS Block 1 will be used on the Boeing Crew Space Transportation (CST-100) system to achieve docking with the International Space Station. Critical functions in the Soft Capture aspect of docking are performed by the LAS. This paper describes the general function of the LAS, the system's key requirements and technical challenges, and the development and qualification approach for the system.

A simple two stage optimization algorithm for constrained power economic dispatch

International Nuclear Information System (INIS)

Huang, G.; Song, K.

1994-01-01

A simple two stage optimization algorithm is proposed and investigated for fast computation of constrained power economic dispatch control problems. The method is a simple demonstration of the hierarchical aggregation-disaggregation (HAD) concept. The algorithm first solves an aggregated problem to obtain an initial solution. This aggregated problem turns out to be classical economic dispatch formulation, and it can be solved in 1% of overall computation time. In the second stage, linear programming method finds optimal solution which satisfies power balance constraints, generation and transmission inequality constraints and security constraints. Implementation of the algorithm for IEEE systems and EPRI Scenario systems shows that the two stage method obtains average speedup ratio 10.64 as compared to classical LP-based method

Linear integral equations and soliton systems

International Nuclear Information System (INIS)

Quispel, G.R.W.

1983-01-01

A study is presented of classical integrable dynamical systems in one temporal and one spatial dimension. The direct linearizations are given of several nonlinear partial differential equations, for example the Korteweg-de Vries equation, the modified Korteweg-de Vries equation, the sine-Gordon equation, the nonlinear Schroedinger equation, and the equation of motion for the isotropic Heisenberg spin chain; the author also discusses several relations between these equations. The Baecklund transformations of these partial differential equations are treated on the basis of a singular transformation of the measure (or equivalently of the plane-wave factor) occurring in the corresponding linear integral equations, and the Baecklund transformations are used to derive the direct linearization of a chain of so-called modified partial differential equations. Finally it is shown that the singular linear integral equations lead in a natural way to the direct linearizations of various nonlinear difference-difference equations. (Auth.)

Dual-range linearized transimpedance amplifier system

Science.gov (United States)

Wessendorf, Kurt O.

2010-11-02

A transimpedance amplifier system is disclosed which simultaneously generates a low-gain output signal and a high-gain output signal from an input current signal using a single transimpedance amplifier having two different feedback loops with different amplification factors to generate two different output voltage signals. One of the feedback loops includes a resistor, and the other feedback loop includes another resistor in series with one or more diodes. The transimpedance amplifier system includes a signal linearizer to linearize one or both of the low- and high-gain output signals by scaling and adding the two output voltage signals from the transimpedance amplifier. The signal linearizer can be formed either as an analog device using one or two summing amplifiers, or alternately can be formed as a digital device using two analog-to-digital converters and a digital signal processor (e.g. a microprocessor or a computer).

Structure Learning in Stochastic Non-linear Dynamical Systems

Science.gov (United States)

Morris, R. D.; Smelyanskiy, V. N.; Luchinsky, D. G.

2005-12-01

A great many systems can be modeled in the non-linear dynamical systems framework, as x˙ = f(x) + ξ(t), where f(x) is the potential function for the system, and ξ(t) is the driving noise. Modeling the potential using a set of basis functions, we derive the posterior for the basis coefficients. A more challenging problem is to determine the set of basis functions that are required to model a particular system. We show that using the Bayesian Information Criteria (BIC) to rank models, and the beam search technique, that we can accurately determine the structure of simple non-linear dynamical system models, and the structure of the coupling between non-linear dynamical systems where the individual systems are known. This last case has important ecological applications, for example in predator-prey systems, where the very structure of the coupling between predator-prey pairs can have great ecological significance.

Choosing health, constrained choices.

Science.gov (United States)

Chee Khoon Chan

2009-12-01

In parallel with the neo-liberal retrenchment of the welfarist state, an increasing emphasis on the responsibility of individuals in managing their own affairs and their well-being has been evident. In the health arena for instance, this was a major theme permeating the UK government's White Paper Choosing Health: Making Healthy Choices Easier (2004), which appealed to an ethos of autonomy and self-actualization through activity and consumption which merited esteem. As a counterpoint to this growing trend of informed responsibilization, constrained choices (constrained agency) provides a useful framework for a judicious balance and sense of proportion between an individual behavioural focus and a focus on societal, systemic, and structural determinants of health and well-being. Constrained choices is also a conceptual bridge between responsibilization and population health which could be further developed within an integrative biosocial perspective one might refer to as the social ecology of health and disease.

Seismic analysis of equipment system with non-linearities such as gap and friction using equivalent linearization method

International Nuclear Information System (INIS)

Murakami, H.; Hirai, T.; Nakata, M.; Kobori, T.; Mizukoshi, K.; Takenaka, Y.; Miyagawa, N.

1989-01-01

Many of the equipment systems of nuclear power plants contain a number of non-linearities, such as gap and friction, due to their mechanical functions. It is desirable to take such non-linearities into account appropriately for the evaluation of the aseismic soundness. However, in usual design works, linear analysis method with rough assumptions is applied from engineering point of view. An equivalent linearization method is considered to be one of the effective analytical techniques to evaluate non-linear responses, provided that errors to a certain extent are tolerated, because it has greater simplicity in analysis and economization in computing time than non-linear analysis. The objective of this paper is to investigate the applicability of the equivalent linearization method to evaluate the maximum earthquake response of equipment systems such as the CANDU Fuelling Machine which has multiple non- linearities

Identification of general linear mechanical systems

Science.gov (United States)

Sirlin, S. W.; Longman, R. W.; Juang, J. N.

1983-01-01

Previous work in identification theory has been concerned with the general first order time derivative form. Linear mechanical systems, a large and important class, naturally have a second order form. This paper utilizes this additional structural information for the purpose of identification. A realization is obtained from input-output data, and then knowledge of the system input, output, and inertia matrices is used to determine a set of linear equations whereby we identify the remaining unknown system matrices. Necessary and sufficient conditions on the number, type and placement of sensors and actuators are given which guarantee identificability, and less stringent conditions are given which guarantee generic identifiability. Both a priori identifiability and a posteriori identifiability are considered, i.e., identifiability being insured prior to obtaining data, and identifiability being assured with a given data set.

System theory as applied differential geometry. [linear system

Science.gov (United States)

Hermann, R.

1979-01-01

The invariants of input-output systems under the action of the feedback group was examined. The approach used the theory of Lie groups and concepts of modern differential geometry, and illustrated how the latter provides a basis for the discussion of the analytic structure of systems. Finite dimensional linear systems in a single independent variable are considered. Lessons of more general situations (e.g., distributed parameter and multidimensional systems) which are increasingly encountered as technology advances are presented.

Flexure Based Linear and Rotary Bearings

Science.gov (United States)

Voellmer, George M. (Inventor)

2016-01-01

A flexure based linear bearing includes top and bottom parallel rigid plates; first and second flexures connecting the top and bottom plates and constraining exactly four degrees of freedom of relative motion of the plates, the four degrees of freedom being X and Y axis translation and rotation about the X and Y axes; and a strut connecting the top and bottom plates and further constraining exactly one degree of freedom of the plates, the one degree of freedom being one of Z axis translation and rotation about the Z axis.

Numerical Estimation of Balanced and Falling States for Constrained Legged Systems

Science.gov (United States)

Mummolo, Carlotta; Mangialardi, Luigi; Kim, Joo H.

2017-08-01

Instability and risk of fall during standing and walking are common challenges for biped robots. While existing criteria from state-space dynamical systems approach or ground reference points are useful in some applications, complete system models and constraints have not been taken into account for prediction and indication of fall for general legged robots. In this study, a general numerical framework that estimates the balanced and falling states of legged systems is introduced. The overall approach is based on the integration of joint-space and Cartesian-space dynamics of a legged system model. The full-body constrained joint-space dynamics includes the contact forces and moments term due to current foot (or feet) support and another term due to altered contact configuration. According to the refined notions of balanced, falling, and fallen, the system parameters, physical constraints, and initial/final/boundary conditions for balancing are incorporated into constrained nonlinear optimization problems to solve for the velocity extrema (representing the maximum perturbation allowed to maintain balance without changing contacts) in the Cartesian space at each center-of-mass (COM) position within its workspace. The iterative algorithm constructs the stability boundary as a COM state-space partition between balanced and falling states. Inclusion in the resulting six-dimensional manifold is a necessary condition for a state of the given system to be balanced under the given contact configuration, while exclusion is a sufficient condition for falling. The framework is used to analyze the balance stability of example systems with various degrees of complexities. The manifold for a 1-degree-of-freedom (DOF) legged system is consistent with the experimental and simulation results in the existing studies for specific controller designs. The results for a 2-DOF system demonstrate the dependency of the COM state-space partition upon joint-space configuration (elbow-up vs

Application of Nearly Linear Solvers to Electric Power System Computation

Science.gov (United States)

Grant, Lisa L.

To meet the future needs of the electric power system, improvements need to be made in the areas of power system algorithms, simulation, and modeling, specifically to achieve a time frame that is useful to industry. If power system time-domain simulations could run in real-time, then system operators would have situational awareness to implement online control and avoid cascading failures, significantly improving power system reliability. Several power system applications rely on the solution of a very large linear system. As the demands on power systems continue to grow, there is a greater computational complexity involved in solving these large linear systems within reasonable time. This project expands on the current work in fast linear solvers, developed for solving symmetric and diagonally dominant linear systems, in order to produce power system specific methods that can be solved in nearly-linear run times. The work explores a new theoretical method that is based on ideas in graph theory and combinatorics. The technique builds a chain of progressively smaller approximate systems with preconditioners based on the system's low stretch spanning tree. The method is compared to traditional linear solvers and shown to reduce the time and iterations required for an accurate solution, especially as the system size increases. A simulation validation is performed, comparing the solution capabilities of the chain method to LU factorization, which is the standard linear solver for power flow. The chain method was successfully demonstrated to produce accurate solutions for power flow simulation on a number of IEEE test cases, and a discussion on how to further improve the method's speed and accuracy is included.

On output regulation for linear systems

NARCIS (Netherlands)

Saberi, Ali; Stoorvogel, Antonie Arij; Sannuti, Peddapullaiah

For both continuous- and discrete-time systems, we revisit the output regulation problem for linear systems. We generalize the problem formulation in order • to expand the class of reference or disturbance signals, • to utilize the derivative or feedforward information of reference signals whenever

Ultra-high Frequency Linear Fiber Optic Systems

CERN Document Server

Lau, Kam

2011-01-01

This book provides an in-depth treatment of both linear fiber-optic systems and their key enabling devices. It presents a concise but rigorous treatment of the theory and practice of analog (linear) fiber-optics links and systems that constitute the foundation of Hybrid Fiber Coax infrastructure in present-day CATV distribution and cable modem Internet access. Emerging applications in remote fiber-optic feed for free-space millimeter wave enterprise campus networks are also described. Issues such as dispersion and interferometric noise are treated quantitatively, and means for mitigating them are explained. This broad but concise text will thus be invaluable not only to students of fiber-optics communication but also to practicing engineers. To the second edition of this book important new aspects of linear fiber-optic transmission technologies are added, such as high level system architectural issues, algorithms for deriving the optimal frequency assignment, directly modulated or externally modulated laser t...

Advanced Computational Methods for Security Constrained Financial Transmission Rights: Structure and Parallelism

Energy Technology Data Exchange (ETDEWEB)

Elbert, Stephen T.; Kalsi, Karanjit; Vlachopoulou, Maria; Rice, Mark J.; Glaesemann, Kurt R.; Zhou, Ning

2012-07-26

Financial Transmission Rights (FTRs) help power market participants reduce price risks associated with transmission congestion. FTRs are issued based on a process of solving a constrained optimization problem with the objective to maximize the FTR social welfare under power flow security constraints. Security constraints for different FTR categories (monthly, seasonal or annual) are usually coupled and the number of constraints increases exponentially with the number of categories. Commercial software for FTR calculation can only provide limited categories of FTRs due to the inherent computational challenges mentioned above. In this paper, a novel non-linear dynamical system (NDS) approach is proposed to solve the optimization problem. The new formulation and performance of the NDS solver is benchmarked against widely used linear programming (LP) solvers like CPLEX™ and tested on large-scale systems using data from the Western Electricity Coordinating Council (WECC). The NDS is demonstrated to outperform the widely used CPLEX algorithms while exhibiting superior scalability. Furthermore, the NDS based solver can be easily parallelized which results in significant computational improvement.

Risk-constrained self-scheduling of a fuel and emission constrained power producer using rolling window procedure

International Nuclear Information System (INIS)

Kazempour, S. Jalal; Moghaddam, Mohsen Parsa

2011-01-01

This work addresses a relevant methodology for self-scheduling of a price-taker fuel and emission constrained power producer in day-ahead correlated energy, spinning reserve and fuel markets to achieve a trade-off between the expected profit and the risk versus different risk levels based on Markowitz's seminal work in the area of portfolio selection. Here, a set of uncertainties including price forecasting errors and available fuel uncertainty are considered. The latter uncertainty arises because of uncertainties in being called for reserve deployment in the spinning reserve market and availability of power plant. To tackle the price forecasting errors, variances of energy, spinning reserve and fuel prices along with their covariances which are due to markets correlation are taken into account using relevant historical data. In order to tackle available fuel uncertainty, a framework for self-scheduling referred to as rolling window is proposed. This risk-constrained self-scheduling framework is therefore formulated and solved as a mixed-integer non-linear programming problem. Furthermore, numerical results for a case study are discussed. (author)

Minimal solution of general dual fuzzy linear systems

International Nuclear Information System (INIS)

Abbasbandy, S.; Otadi, M.; Mosleh, M.

2008-01-01

Fuzzy linear systems of equations, play a major role in several applications in various area such as engineering, physics and economics. In this paper, we investigate the existence of a minimal solution of general dual fuzzy linear equation systems. Two necessary and sufficient conditions for the minimal solution existence are given. Also, some examples in engineering and economic are considered

On the stability of non-linear systems

International Nuclear Information System (INIS)

Guelman, M.

1968-09-01

A study is made of the absolute stability of nonlinear systems, using Liapounov's second method and taking into account the results obtained from V.M. Popov's work. The results already established are first presented, in particular concerning the frequency domain criterions for absolute stability of automatic control systems containing one single non linearity. The results have been extended to show the existence of a limiting parabola. New use is then made of the methods studied for deriving absolute stability criterions for a system containing a different type of non linearity. Finally, the results obtained are considered from the point of view of Aizerman's conjecture. (author) [fr

Linear dynamic coupling in geared rotor systems

Science.gov (United States)

David, J. W.; Mitchell, L. D.

1986-01-01

The effects of high frequency oscillations caused by the gear mesh, on components of a geared system that can be modeled as rigid discs are analyzed using linear dynamic coupling terms. The coupled, nonlinear equations of motion for a disc attached to a rotating shaft are presented. The results of a trial problem analysis show that the inclusion of the linear dynamic coupling terms can produce significant changes in the predicted response of geared rotor systems, and that the produced sideband responses are greater than the unbalanced response. The method is useful in designing gear drives for heavy-lift helicopters, industrial speed reducers, naval propulsion systems, and heavy off-road equipment.

ORACLS: A system for linear-quadratic-Gaussian control law design

Science.gov (United States)

Armstrong, E. S.

1978-01-01

A modern control theory design package (ORACLS) for constructing controllers and optimal filters for systems modeled by linear time-invariant differential or difference equations is described. Numerical linear-algebra procedures are used to implement the linear-quadratic-Gaussian (LQG) methodology of modern control theory. Algorithms are included for computing eigensystems of real matrices, the relative stability of a matrix, factored forms for nonnegative definite matrices, the solutions and least squares approximations to the solutions of certain linear matrix algebraic equations, the controllability properties of a linear time-invariant system, and the steady state covariance matrix of an open-loop stable system forced by white noise. Subroutines are provided for solving both the continuous and discrete optimal linear regulator problems with noise free measurements and the sampled-data optimal linear regulator problem. For measurement noise, duality theory and the optimal regulator algorithms are used to solve the continuous and discrete Kalman-Bucy filter problems. Subroutines are also included which give control laws causing the output of a system to track the output of a prescribed model.

System identication of a linearized hysteretic system using covariance driven stochastic subspace identication

DEFF Research Database (Denmark)

Bajric, Anela

A single mass Bouc-Wen oscillator with linear static restoring force contribution is approximated by an equivalent linear system. The aim of the linearized model is to emulate the correct force-displacement response of the Bouc-Wenmodel with characteristic hysteretic behaviour. The linearized mod...

Maximum Constrained Directivity of Oversteered End-Fire Sensor Arrays

Directory of Open Access Journals (Sweden)

Andrea Trucco

2015-06-01

Full Text Available For linear arrays with fixed steering and an inter-element spacing smaller than one half of the wavelength, end-fire steering of a data-independent beamformer offers better directivity than broadside steering. The introduction of a lower bound on the white noise gain ensures the necessary robustness against random array errors and sensor mismatches. However, the optimum broadside performance can be obtained using a simple processing architecture, whereas the optimum end-fire performance requires a more complicated system (because complex weight coefficients are needed. In this paper, we reconsider the oversteering technique as a possible way to simplify the processing architecture of equally spaced end-fire arrays. We propose a method for computing the amount of oversteering and the related real-valued weight vector that allows the constrained directivity to be maximized for a given inter-element spacing. Moreover, we verify that the maximized oversteering performance is very close to the optimum end-fire performance. We conclude that optimized oversteering is a viable method for designing end-fire arrays that have better constrained directivity than broadside arrays but with a similar implementation complexity. A numerical simulation is used to perform a statistical analysis, which confirms that the maximized oversteering performance is robust against sensor mismatches.

A convex optimization approach for solving large scale linear systems

Directory of Open Access Journals (Sweden)

Debora Cores

2017-01-01

Full Text Available The well-known Conjugate Gradient (CG method minimizes a strictly convex quadratic function for solving large-scale linear system of equations when the coefficient matrix is symmetric and positive definite. In this work we present and analyze a non-quadratic convex function for solving any large-scale linear system of equations regardless of the characteristics of the coefficient matrix. For finding the global minimizers, of this new convex function, any low-cost iterative optimization technique could be applied. In particular, we propose to use the low-cost globally convergent Spectral Projected Gradient (SPG method, which allow us to extend this optimization approach for solving consistent square and rectangular linear system, as well as linear feasibility problem, with and without convex constraints and with and without preconditioning strategies. Our numerical results indicate that the new scheme outperforms state-of-the-art iterative techniques for solving linear systems when the symmetric part of the coefficient matrix is indefinite, and also for solving linear feasibility problems.

Observability of linear systems with saturated outputs

NARCIS (Netherlands)

Koplon, R.; Sontag, E.D.; Hautus, M.L.J.

1994-01-01

We present necessary and sufficient conditions for observability of the class of output-saturated systems. These are linear systems whose output passes through a saturation function before it can be measured.

On non-linear dynamics of a coupled electro-mechanical system

DEFF Research Database (Denmark)

Darula, Radoslav; Sorokin, Sergey

2012-01-01

Electro-mechanical devices are an example of coupled multi-disciplinary weakly non-linear systems. Dynamics of such systems is described in this paper by means of two mutually coupled differential equations. The first one, describing an electrical system, is of the first order and the second one...... excitation. The results are verified using a numerical model created in MATLAB Simulink environment. Effect of non-linear terms on dynamical response of the coupled system is investigated; the backbone and envelope curves are analyzed. The two phenomena, which exist in the electro-mechanical system: (a......, for mechanical system, is of the second order. The governing equations are coupled via linear and weakly non-linear terms. A classical perturbation method, a method of multiple scales, is used to find a steadystate response of the electro-mechanical system exposed to a harmonic close-resonance mechanical...

Simultaneous Balancing and Model Reduction of Switched Linear Systems

NARCIS (Netherlands)

Monshizadeh, Nima; Trentelman, Hendrikus; Camlibel, M.K.

2011-01-01

In this paper, first, balanced truncation of linear systems is revisited. Then, simultaneous balancing of multiple linear systems is investigated. Necessary and sufficient conditions are introduced to identify the case where simultaneous balancing is possible. The validity of these conditions is not

The linear sizes tolerances and fits system modernization

Science.gov (United States)

Glukhov, V. I.; Grinevich, V. A.; Shalay, V. V.

2018-04-01

The study is carried out on the urgent topic for technical products quality providing in the tolerancing process of the component parts. The aim of the paper is to develop alternatives for improving the system linear sizes tolerances and dimensional fits in the international standard ISO 286-1. The tasks of the work are, firstly, to classify as linear sizes the elements additionally linear coordinating sizes that determine the detail elements location and, secondly, to justify the basic deviation of the tolerance interval for the element's linear size. The geometrical modeling method of real details elements, the analytical and experimental methods are used in the research. It is shown that the linear coordinates are the dimensional basis of the elements linear sizes. To standardize the accuracy of linear coordinating sizes in all accuracy classes, it is sufficient to select in the standardized tolerance system only one tolerance interval with symmetrical deviations: Js for internal dimensional elements (holes) and js for external elements (shafts). The main deviation of this coordinating tolerance is the average zero deviation, which coincides with the nominal value of the coordinating size. Other intervals of the tolerance system are remained for normalizing the accuracy of the elements linear sizes with a fundamental change in the basic deviation of all tolerance intervals is the maximum deviation corresponding to the limit of the element material: EI is the lower tolerance for the of the internal elements (holes) sizes and es is the upper tolerance deviation for the outer elements (shafts) sizes. It is the sizes of the material maximum that are involved in the of the dimensional elements mating of the shafts and holes and determine the fits type.

Rf system specifications for a linear accelerator

International Nuclear Information System (INIS)

Young, A.; Eaton, L.E.

1992-01-01

A linear accelerator contains many systems; however, the most complex and costly is the RF system. The goal of an RF system is usually simply stated as maintaining the phase and amplitude of the RF signal within a given tolerance to accelerate the charged particle beam. An RF system that drives a linear accelerator needs a complete system specification, which should contain specifications for all the subsystems (i.e., high-power RF, low-level RF, RF generation/distribution, and automation control). This paper defines a format for the specifications of these subsystems and discusses each RF subsystem independently to provide a comprehensive understanding of the function of each subsystem. This paper concludes with an example of a specification spreadsheet allowing one to input the specifications of a subsystem. Thus, some fundamental parameters (i.e., the cost and size) of the RF system can be determined

Relative null controllability of linear systems with multiple delays in ...

African Journals Online (AJOL)

varying multiple delays in state and control are developed. If the uncontrolled system is uniformly asymptotically stable, and if the linear system is controllable, then the linear system is null controllable. Journal of the Nigerian Association of ...

Solution of a Complex Least Squares Problem with Constrained Phase.

Science.gov (United States)

Bydder, Mark

2010-12-30

The least squares solution of a complex linear equation is in general a complex vector with independent real and imaginary parts. In certain applications in magnetic resonance imaging, a solution is desired such that each element has the same phase. A direct method for obtaining the least squares solution to the phase constrained problem is described.

Superconducting linear accelerator system for NSC

Indian Academy of Sciences (India)

59, No. 5. — journal of. November 2002 physics pp. 849–858. Superconducting linear accelerator system for NSC ... cryogenics facility, RF electronics development, facilities for fabricating niobium resonators indige- ... Prototype resonator was.

Linear quadratic Gaussian balancing for discrete-time infinite-dimensional linear systems

NARCIS (Netherlands)

Opmeer, MR; Curtain, RF

2004-01-01

In this paper, we study the existence of linear quadratic Gaussian (LQG)-balanced realizations for discrete-time infinite-dimensional systems. LQG-balanced realizations are those for which the smallest nonnegative self-adjoint solutions of the control and filter Riccati equations are equal. We show

A separation theorem for the stochastic sampled-data LQG problem. [control of continuous linear plant disturbed by white noise

Science.gov (United States)

Halyo, N.; Caglayan, A. K.

1976-01-01

This paper considers the control of a continuous linear plant disturbed by white plant noise when the control is constrained to be a piecewise constant function of time; i.e. a stochastic sampled-data system. The cost function is the integral of quadratic error terms in the state and control, thus penalizing errors at every instant of time while the plant noise disturbs the system continuously. The problem is solved by reducing the constrained continuous problem to an unconstrained discrete one. It is shown that the separation principle for estimation and control still holds for this problem when the plant disturbance and measurement noise are Gaussian.

Tikhonov theorem for linear hyperbolic systems

OpenAIRE

Tang , Ying; Prieur , Christophe; Girard , Antoine

2015-01-01

International audience; A class of linear systems of conservation laws with a small perturbation parameter is introduced. By setting the perturbation parameter to zero, two subsystems, the reduced system standing for the slow dynamics and the boundary-layer system representing the fast dynamics, are computed. It is first proved that the exponential stability of the full system implies the stability of both subsystems. Secondly, a counter example is given to indicate that the converse is not t...

Conduction cooling systems for linear accelerator cavities

Science.gov (United States)

Kephart, Robert

2017-05-02

A conduction cooling system for linear accelerator cavities. The system conducts heat from the cavities to a refrigeration unit using at least one cavity cooler interconnected with a cooling connector. The cavity cooler and cooling connector are both made from solid material having a very high thermal conductivity of approximately 1.times.10.sup.4 W m.sup.-1 K.sup.-1 at temperatures of approximately 4 degrees K. This allows for very simple and effective conduction of waste heat from the linear accelerator cavities to the cavity cooler, along the cooling connector, and thence to the refrigeration unit.

Theoretical analysis of balanced truncation for linear switched systems

DEFF Research Database (Denmark)

Petreczky, Mihaly; Wisniewski, Rafal; Leth, John-Josef

2012-01-01

In this paper we present theoretical analysis of model reduction of linear switched systems based on balanced truncation, presented in [1,2]. More precisely, (1) we provide a bound on the estimation error using L2 gain, (2) we provide a system theoretic interpretation of grammians and their singu......In this paper we present theoretical analysis of model reduction of linear switched systems based on balanced truncation, presented in [1,2]. More precisely, (1) we provide a bound on the estimation error using L2 gain, (2) we provide a system theoretic interpretation of grammians...... for showing this independence is realization theory of linear switched systems. [1] H. R. Shaker and R. Wisniewski, "Generalized gramian framework for model/controller order reduction of switched systems", International Journal of Systems Science, Vol. 42, Issue 8, 2011, 1277-1291. [2] H. R. Shaker and R....... Wisniewski, "Switched Systems Reduction Framework Based on Convex Combination of Generalized Gramians", Journal of Control Science and Engineering, 2009....

A comparison between linear and toroidal Extrap systems

International Nuclear Information System (INIS)

Lehnert, B.

1988-09-01

The Extrap scheme consists of a Z-pinch immersed in an octupole field generated by currents in a set of external conductors. A comparison between linear and toroidal Extrap geometry is made in this paper. As compared to toroidal systems, linear geometry has the advantages of relative simplicity and of a current drive by means of electrodes. Linear devices are convenient for basic studies of Extrap, at moderately high pinch currents and plasma temperatures. Within the parameter ranges of experiments at high pinch currents and plasma temperatures, linear systems have on the other hand some substantial disadvantages, on account of the plasma interaction with the end regions. This results in a limitation of the energy confinement time, and leads in the case of an ohmically heated plasma to excessively high plasma densities and small pinch radii which also complicate the introduction of the external conductors. (author)

Correlated Levy Noise in Linear Dynamical Systems

International Nuclear Information System (INIS)

Srokowski, T.

2011-01-01

Linear dynamical systems, driven by a non-white noise which has the Levy distribution, are analysed. Noise is modelled by a specific stochastic process which is defined by the Langevin equation with a linear force and the Levy distributed symmetric white noise. Correlation properties of the process are discussed. The Fokker-Planck equation driven by that noise is solved. Distributions have the Levy shape and their width, for a given time, is smaller than for processes in the white noise limit. Applicability of the adiabatic approximation in the case of the linear force is discussed. (author)

Reliability modelling and simulation of switched linear system ...

African Journals Online (AJOL)

Reliability modelling and simulation of switched linear system control using temporal databases. ... design of fault-tolerant real-time switching systems control and modelling embedded micro-schedulers for complex systems maintenance.

Constrained Vapor Bubble Experiment

Science.gov (United States)

Gokhale, Shripad; Plawsky, Joel; Wayner, Peter C., Jr.; Zheng, Ling; Wang, Ying-Xi

2002-11-01

Microgravity experiments on the Constrained Vapor Bubble Heat Exchanger, CVB, are being developed for the International Space Station. In particular, we present results of a precursory experimental and theoretical study of the vertical Constrained Vapor Bubble in the Earth's environment. A novel non-isothermal experimental setup was designed and built to study the transport processes in an ethanol/quartz vertical CVB system. Temperature profiles were measured using an in situ PC (personal computer)-based LabView data acquisition system via thermocouples. Film thickness profiles were measured using interferometry. A theoretical model was developed to predict the curvature profile of the stable film in the evaporator. The concept of the total amount of evaporation, which can be obtained directly by integrating the experimental temperature profile, was introduced. Experimentally measured curvature profiles are in good agreement with modeling results. For microgravity conditions, an analytical expression, which reveals an inherent relation between temperature and curvature profiles, was derived.

A supply function model for representing the strategic bidding of the producers in constrained electricity markets

International Nuclear Information System (INIS)

Bompard, Ettore; Napoli, Roberto; Lu, Wene; Jiang, Xiuchen

2010-01-01

The modeling of the bidding behaviour of the producer is a key-point in the modeling and simulation of the competitive electricity markets. In our paper, the linear supply function model is applied so as to find the Supply Function Equilibrium analytically. It also proposed a new and efficient approach to find SFEs for the network constrained electricity markets by finding the best slope of the supply function with the help of changing the intercept, and the method can be applied on the large systems. The approach proposed is applied to study IEEE-118 bus test systems and the comparison between bidding slope and bidding intercept is presented, as well, with reference to the test system. (author)

Damped oscillations of linear systems a mathematical introduction

CERN Document Server

Veselić, Krešimir

2011-01-01

The theory of linear damped oscillations was originally developed more than hundred years ago and is still of vital research interest to engineers, mathematicians and physicists alike. This theory plays a central role in explaining the stability of mechanical structures in civil engineering, but it also has applications in other fields such as electrical network systems and quantum mechanics. This volume gives an introduction to linear finite dimensional damped systems as they are viewed by an applied mathematician. After a short overview of the physical principles leading to the linear system model, a largely self-contained mathematical theory for this model is presented. This includes the geometry of the underlying indefinite metric space, spectral theory of J-symmetric matrices and the associated quadratic eigenvalue problem. Particular attention is paid to the sensitivity issues which influence numerical computations. Finally, several recent research developments are included, e.g. Lyapunov stability and ...

A parallel solver for huge dense linear systems

Science.gov (United States)

Badia, J. M.; Movilla, J. L.; Climente, J. I.; Castillo, M.; Marqués, M.; Mayo, R.; Quintana-Ortí, E. S.; Planelles, J.

2011-11-01

HDSS (Huge Dense Linear System Solver) is a Fortran Application Programming Interface (API) to facilitate the parallel solution of very large dense systems to scientists and engineers. The API makes use of parallelism to yield an efficient solution of the systems on a wide range of parallel platforms, from clusters of processors to massively parallel multiprocessors. It exploits out-of-core strategies to leverage the secondary memory in order to solve huge linear systems O(100.000). The API is based on the parallel linear algebra library PLAPACK, and on its Out-Of-Core (OOC) extension POOCLAPACK. Both PLAPACK and POOCLAPACK use the Message Passing Interface (MPI) as the communication layer and BLAS to perform the local matrix operations. The API provides a friendly interface to the users, hiding almost all the technical aspects related to the parallel execution of the code and the use of the secondary memory to solve the systems. In particular, the API can automatically select the best way to store and solve the systems, depending of the dimension of the system, the number of processes and the main memory of the platform. Experimental results on several parallel platforms report high performance, reaching more than 1 TFLOP with 64 cores to solve a system with more than 200 000 equations and more than 10 000 right-hand side vectors. New version program summaryProgram title: Huge Dense System Solver (HDSS) Catalogue identifier: AEHU_v1_1 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEHU_v1_1.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: Standard CPC licence, http://cpc.cs.qub.ac.uk/licence/licence.html No. of lines in distributed program, including test data, etc.: 87 062 No. of bytes in distributed program, including test data, etc.: 1 069 110 Distribution format: tar.gz Programming language: Fortran90, C Computer: Parallel architectures: multiprocessors, computer clusters Operating system

SNR Estimation in Linear Systems with Gaussian Matrices

KAUST Repository

Suliman, Mohamed Abdalla Elhag; Alrashdi, Ayed; Ballal, Tarig; Al-Naffouri, Tareq Y.

2017-01-01

This letter proposes a highly accurate algorithm to estimate the signal-to-noise ratio (SNR) for a linear system from a single realization of the received signal. We assume that the linear system has a Gaussian matrix with one sided left correlation. The unknown entries of the signal and the noise are assumed to be independent and identically distributed with zero mean and can be drawn from any distribution. We use the ridge regression function of this linear model in company with tools and techniques adapted from random matrix theory to achieve, in closed form, accurate estimation of the SNR without prior statistical knowledge on the signal or the noise. Simulation results show that the proposed method is very accurate.

SNR Estimation in Linear Systems with Gaussian Matrices

KAUST Repository

Suliman, Mohamed Abdalla Elhag

2017-09-27

This letter proposes a highly accurate algorithm to estimate the signal-to-noise ratio (SNR) for a linear system from a single realization of the received signal. We assume that the linear system has a Gaussian matrix with one sided left correlation. The unknown entries of the signal and the noise are assumed to be independent and identically distributed with zero mean and can be drawn from any distribution. We use the ridge regression function of this linear model in company with tools and techniques adapted from random matrix theory to achieve, in closed form, accurate estimation of the SNR without prior statistical knowledge on the signal or the noise. Simulation results show that the proposed method is very accurate.

An extended GS method for dense linear systems

Science.gov (United States)

Niki, Hiroshi; Kohno, Toshiyuki; Abe, Kuniyoshi

2009-09-01

Davey and Rosindale [K. Davey, I. Rosindale, An iterative solution scheme for systems of boundary element equations, Internat. J. Numer. Methods Engrg. 37 (1994) 1399-1411] derived the GSOR method, which uses an upper triangular matrix [Omega] in order to solve dense linear systems. By applying functional analysis, the authors presented an expression for the optimum [Omega]. Moreover, Davey and Bounds [K. Davey, S. Bounds, A generalized SOR method for dense linear systems of boundary element equations, SIAM J. Comput. 19 (1998) 953-967] also introduced further interesting results. In this note, we employ a matrix analysis approach to investigate these schemes, and derive theorems that compare these schemes with existing preconditioners for dense linear systems. We show that the convergence rate of the Gauss-Seidel method with preconditioner PG is superior to that of the GSOR method. Moreover, we define some splittings associated with the iterative schemes. Some numerical examples are reported to confirm the theoretical analysis. We show that the EGS method with preconditioner produces an extremely small spectral radius in comparison with the other schemes considered.

Optimal Linear Responses for Markov Chains and Stochastically Perturbed Dynamical Systems

Science.gov (United States)

Antown, Fadi; Dragičević, Davor; Froyland, Gary

2018-03-01

The linear response of a dynamical system refers to changes to properties of the system when small external perturbations are applied. We consider the little-studied question of selecting an optimal perturbation so as to (i) maximise the linear response of the equilibrium distribution of the system, (ii) maximise the linear response of the expectation of a specified observable, and (iii) maximise the linear response of the rate of convergence of the system to the equilibrium distribution. We also consider the inhomogeneous, sequential, or time-dependent situation where the governing dynamics is not stationary and one wishes to select a sequence of small perturbations so as to maximise the overall linear response at some terminal time. We develop the theory for finite-state Markov chains, provide explicit solutions for some illustrative examples, and numerically apply our theory to stochastically perturbed dynamical systems, where the Markov chain is replaced by a matrix representation of an approximate annealed transfer operator for the random dynamical system.

Anti-D3 branes and moduli in non-linear supergravity

Science.gov (United States)

Garcia del Moral, Maria P.; Parameswaran, Susha; Quiroz, Norma; Zavala, Ivonne

2017-10-01

Anti-D3 branes and non-perturbative effects in flux compactifications spontaneously break supersymmetry and stabilise moduli in a metastable de Sitter vacua. The low energy 4D effective field theory description for such models would be a supergravity theory with non-linearly realised supersymmetry. Guided by string theory modular symmetry, we compute this non-linear supergravity theory, including dependence on all bulk moduli. Using either a constrained chiral superfield or a constrained vector field, the uplifting contribution to the scalar potential from the anti-D3 brane can be parameterised either as an F-term or Fayet-Iliopoulos D-term. Using again the modular symmetry, we show that 4D non-linear supergravities that descend from string theory have an enhanced protection from quantum corrections by non-renormalisation theorems. The superpotential giving rise to metastable de Sitter vacua is robust against perturbative string-loop and α' corrections.

On deformations of linear differential systems

NARCIS (Netherlands)

Gontsov, R.R.; Poberezhnyi, V.A.; Helminck, G.F.

2011-01-01

This article concerns deformations of meromorphic linear differential systems. Problems relating to their existence and classification are reviewed, and the global and local behaviour of solutions to deformation equations in a neighbourhood of their singular set is analysed. Certain classical

Popov–Belevitch–Hautus type tests for the controllability of linear complementarity systems

NARCIS (Netherlands)

Camlibel, M. Kanat

2007-01-01

It is well-known that checking certain controllability properties of very simple piecewise linear systems are undecidable problems. This paper deals with the controllability problem of a class of piecewise linear systems, known as linear complementarity systems. By exploiting the underlying

A method for fitting regression splines with varying polynomial order in the linear mixed model.

Science.gov (United States)

Edwards, Lloyd J; Stewart, Paul W; MacDougall, James E; Helms, Ronald W

2006-02-15

The linear mixed model has become a widely used tool for longitudinal analysis of continuous variables. The use of regression splines in these models offers the analyst additional flexibility in the formulation of descriptive analyses, exploratory analyses and hypothesis-driven confirmatory analyses. We propose a method for fitting piecewise polynomial regression splines with varying polynomial order in the fixed effects and/or random effects of the linear mixed model. The polynomial segments are explicitly constrained by side conditions for continuity and some smoothness at the points where they join. By using a reparameterization of this explicitly constrained linear mixed model, an implicitly constrained linear mixed model is constructed that simplifies implementation of fixed-knot regression splines. The proposed approach is relatively simple, handles splines in one variable or multiple variables, and can be easily programmed using existing commercial software such as SAS or S-plus. The method is illustrated using two examples: an analysis of longitudinal viral load data from a study of subjects with acute HIV-1 infection and an analysis of 24-hour ambulatory blood pressure profiles.

Uniform discretizations: a quantization procedure for totally constrained systems including gravity

Energy Technology Data Exchange (ETDEWEB)

Campiglia, Miguel [Instituto de Fisica, Facultad de Ciencias, Igua 4225, esq. Mataojo, Montevideo (Uruguay); Di Bartolo, Cayetano [Departamento de Fisica, Universidad Simon BolIvar, Aptdo. 89000, Caracas 1080-A (Venezuela); Gambini, Rodolfo [Instituto de Fisica, Facultad de Ciencias, Igua 4225, esq. Mataojo, Montevideo (Uruguay); Pullin, Jorge [Department of Physics and Astronomy, Louisiana State University, Baton Rouge, LA 70803-4001 (United States)

2007-05-15

We present a new method for the quantization of totally constrained systems including general relativity. The method consists in constructing discretized theories that have a well defined and controlled continuum limit. The discrete theories are constraint-free and can be readily quantized. This provides a framework where one can introduce a relational notion of time and that nevertheless approximates in a well defined fashion the theory of interest. The method is equivalent to the group averaging procedure for many systems where the latter makes sense and provides a generalization otherwise. In the continuum limit it can be shown to contain, under certain assumptions, the 'master constraint' of the 'Phoenix project'. It also provides a correspondence principle with the classical theory that does not require to consider the semiclassical limit.

A study on switched linear system identification using game ...

African Journals Online (AJOL)

A study on switched linear system identification using game-theoretic strategies and neural computing. ... This study deals with application of game-theoretic strategies and neural computing to switched linear ... AJOL African Journals Online.

Conjugate gradient type methods for linear systems with complex symmetric coefficient matrices

Science.gov (United States)

Freund, Roland

1989-01-01

We consider conjugate gradient type methods for the solution of large sparse linear system Ax equals b with complex symmetric coefficient matrices A equals A(T). Such linear systems arise in important applications, such as the numerical solution of the complex Helmholtz equation. Furthermore, most complex non-Hermitian linear systems which occur in practice are actually complex symmetric. We investigate conjugate gradient type iterations which are based on a variant of the nonsymmetric Lanczos algorithm for complex symmetric matrices. We propose a new approach with iterates defined by a quasi-minimal residual property. The resulting algorithm presents several advantages over the standard biconjugate gradient method. We also include some remarks on the obvious approach to general complex linear systems by solving equivalent real linear systems for the real and imaginary parts of x. Finally, numerical experiments for linear systems arising from the complex Helmholtz equation are reported.

Dynamically constrained ensemble perturbations – application to tides on the West Florida Shelf

Directory of Open Access Journals (Sweden)

F. Lenartz

2009-07-01

Full Text Available A method is presented to create an ensemble of perturbations that satisfies linear dynamical constraints. A cost function is formulated defining the probability of each perturbation. It is shown that the perturbations created with this approach take the land-sea mask into account in a similar way as variational analysis techniques. The impact of the land-sea mask is illustrated with an idealized configuration of a barrier island. Perturbations with a spatially variable correlation length can be also created by this approach. The method is applied to a realistic configuration of the West Florida Shelf to create perturbations of the M2 tidal parameters for elevation and depth-averaged currents. The perturbations are weakly constrained to satisfy the linear shallow-water equations. Despite that the constraint is derived from an idealized assumption, it is shown that this approach is applicable to a non-linear and baroclinic model. The amplitude of spurious transient motions created by constrained perturbations of initial and boundary conditions is significantly lower compared to perturbing the variables independently or to using only the momentum equation to compute the velocity perturbations from the elevation.

Disturbance Decoupling of Switched Linear Systems

NARCIS (Netherlands)

Yurtseven, E.; Heemels, W.P.M.H.; Camlibel, M.K.

2010-01-01

In this paper we consider disturbance decoupling problems for switched linear systems. We will provide necessary and sufficient conditions for three different versions of disturbance decoupling, which differ based on which signals are considered to be the disturbance. In the first version the

Feedback linearizing control of a MIMO power system

Science.gov (United States)

Ilyes, Laszlo

Prior research has demonstrated that either the mechanical or electrical subsystem of a synchronous electric generator may be controlled using single-input single-output (SISO) nonlinear feedback linearization. This research suggests a new approach which applies nonlinear feedback linearization to a multi-input multi-output (MIMO) model of the synchronous electric generator connected to an infinite bus load model. In this way, the electrical and mechanical subsystems may be linearized and simultaneously decoupled through the introduction of a pair of auxiliary inputs. This allows well known, linear, SISO control methods to be effectively applied to the resulting systems. The derivation of the feedback linearizing control law is presented in detail, including a discussion on the use of symbolic math processing as a development tool. The linearizing and decoupling properties of the control law are validated through simulation. And finally, the robustness of the control law is demonstrated.

A first-order multigrid method for bound-constrained convex optimization

Czech Academy of Sciences Publication Activity Database

Kočvara, Michal; Mohammed, S.

2016-01-01

Roč. 31, č. 3 (2016), s. 622-644 ISSN 1055-6788 R&D Projects: GA ČR(CZ) GAP201/12/0671 Grant - others:European Commission - EC(XE) 313781 Institutional support: RVO:67985556 Keywords : bound-constrained optimization * multigrid methods * linear complementarity problems Subject RIV: BA - General Mathematics Impact factor: 1.023, year: 2016 http://library.utia.cas.cz/separaty/2016/MTR/kocvara-0460326.pdf

Nonautonomous linear system of the terrestrial carbon cycle

Science.gov (United States)

Luo, Y.

2012-12-01

Carbon cycle has been studied by uses of observation through various networks, field and laboratory experiments, and simulation models. Much less has been done on theoretical thinking and analysis to understand fundament properties of carbon cycle and then guide observatory, experimental, and modeling research. This presentation is to explore what would be the theoretical properties of terrestrial carbon cycle and how those properties can be used to make observatory, experimental, and modeling research more effective. Thousands of published data sets from litter decomposition and soil incubation studies almost all indicate that decay processes of litter and soil organic carbon can be well described by first order differential equations with one or more pools. Carbon pool dynamics in plants and soil after disturbances (e.g., wildfire, clear-cut of forests, and plows of soil for cropping) and during natural recovery or ecosystem restoration also exhibit characteristics of first-order linear systems. Thus, numerous lines of empirical evidence indicate that the terrestrial carbon cycle can be adequately described as a nonautonomous linear system. The linearity reflects the nature of the carbon cycle that carbon, once fixed by photosynthesis, is linearly transferred among pools within an ecosystem. The linear carbon transfer, however, is modified by nonlinear functions of external forcing variables. In addition, photosynthetic carbon influx is also nonlinearly influenced by external variables. This nonautonomous linear system can be mathematically expressed by a first-order linear ordinary matrix equation. We have recently used this theoretical property of terrestrial carbon cycle to develop a semi-analytic solution of spinup. The new methods have been applied to five global land models, including NCAR's CLM and CABLE models and can computationally accelerate spinup by two orders of magnitude. We also use this theoretical property to develop an analytic framework to

Linearization of Nonautonomous Impulsive System with Nonuniform Exponential Dichotomy

Directory of Open Access Journals (Sweden)

Yongfei Gao

2014-01-01

Full Text Available This paper gives a version of Hartman-Grobman theorem for the impulsive differential equations. We assume that the linear impulsive system has a nonuniform exponential dichotomy. Under some suitable conditions, we proved that the nonlinear impulsive system is topologically conjugated to its linear system. Indeed, we do construct the topologically equivalent function (the transformation. Moreover, the method to prove the topological conjugacy is quite different from those in previous works (e.g., see Barreira and Valls, 2006.

Partial Linearization of Mechanical Systems with Application to Observer Design

NARCIS (Netherlands)

Sarras, Ioannis; Venkatraman, Aneesh; Ortega, Romeo; Schaft, Arjan van der

2008-01-01

We consider general mechanical systems and establish a necessary and sufficient condition for the existence of a suitable change in the generalized momentum coordinates such that the new dynamics become linear in the transformed momenta. The class of systems which can be (partially) linearized by

FAST modularization framework for wind turbine simulation: full-system linearization

Science.gov (United States)

Jonkman, J. M.; Jonkman, B. J.

2016-09-01

The wind engineering community relies on multiphysics engineering software to run nonlinear time-domain simulations e.g. for design-standards-based loads analysis. Although most physics involved in wind energy are nonlinear, linearization of the underlying nonlinear system equations is often advantageous to understand the system response and exploit well- established methods and tools for analyzing linear systems. This paper presents the development and verification of the new linearization functionality of the open-source engineering tool FAST v8 for land-based wind turbines, as well as the concepts and mathematical background needed to understand and apply it correctly.

A Proposed Method for Solving Fuzzy System of Linear Equations

Directory of Open Access Journals (Sweden)

Reza Kargar

2014-01-01

Full Text Available This paper proposes a new method for solving fuzzy system of linear equations with crisp coefficients matrix and fuzzy or interval right hand side. Some conditions for the existence of a fuzzy or interval solution of m×n linear system are derived and also a practical algorithm is introduced in detail. The method is based on linear programming problem. Finally the applicability of the proposed method is illustrated by some numerical examples.

Virtual Estimator for Piecewise Linear Systems Based on Observability Analysis

Science.gov (United States)

Morales-Morales, Cornelio; Adam-Medina, Manuel; Cervantes, Ilse; Vela-Valdés and, Luis G.; García Beltrán, Carlos Daniel

2013-01-01

This article proposes a virtual sensor for piecewise linear systems based on observability analysis that is in function of a commutation law related with the system's outpu. This virtual sensor is also known as a state estimator. Besides, it presents a detector of active mode when the commutation sequences of each linear subsystem are arbitrary and unknown. For the previous, this article proposes a set of virtual estimators that discern the commutation paths of the system and allow estimating their output. In this work a methodology in order to test the observability for piecewise linear systems with discrete time is proposed. An academic example is presented to show the obtained results. PMID:23447007

Symmetric linear systems - An application of algebraic systems theory

Science.gov (United States)

Hazewinkel, M.; Martin, C.

1983-01-01

Dynamical systems which contain several identical subsystems occur in a variety of applications ranging from command and control systems and discretization of partial differential equations, to the stability augmentation of pairs of helicopters lifting a large mass. Linear models for such systems display certain obvious symmetries. In this paper, we discuss how these symmetries can be incorporated into a mathematical model that utilizes the modern theory of algebraic systems. Such systems are inherently related to the representation theory of algebras over fields. We will show that any control scheme which respects the dynamical structure either implicitly or explicitly uses the underlying algebra.

A new linearized equation for servo valve in hydraulic control systems

International Nuclear Information System (INIS)

Kim, Tae Hyung; Lee, Ill Yeong

2002-01-01

In the procedure of the hydraulic control system analysis, a linearized approximate equation described by the first order term of Taylor's series has been widely used. Such a linearized equation is effective just near the operating point. And, as of now, there are no general standards on how to determine the operating point of a servo valve in the process of applying the linearized equation. So, in this study, a new linearized equation for valve characteristics is proposed as a modified form of the existing linearized equation. And, a method for selecting an optimal operating point is proposed for the new linearized equation. The effectiveness of the new linearized equation is confirmed through numerical simulations and experiments for a model hydraulic control system

Geometrically constrained kinematic global navigation satellite systems positioning: Implementation and performance

Science.gov (United States)

Asgari, Jamal; Mohammadloo, Tannaz H.; Amiri-Simkooei, Ali Reza

2015-09-01

GNSS kinematic techniques are capable of providing precise coordinates in extremely short observation time-span. These methods usually determine the coordinates of an unknown station with respect to a reference one. To enhance the precision, accuracy, reliability and integrity of the estimated unknown parameters, GNSS kinematic equations are to be augmented by possible constraints. Such constraints could be derived from the geometric relation of the receiver positions in motion. This contribution presents the formulation of the constrained kinematic global navigation satellite systems positioning. Constraints effectively restrict the definition domain of the unknown parameters from the three-dimensional space to a subspace defined by the equation of motion. To test the concept of the constrained kinematic positioning method, the equation of a circle is employed as a constraint. A device capable of moving on a circle was made and the observations from 11 positions on the circle were analyzed. Relative positioning was conducted by considering the center of the circle as the reference station. The equation of the receiver's motion was rewritten in the ECEF coordinates system. A special attention is drawn onto how a constraint is applied to kinematic positioning. Implementing the constraint in the positioning process provides much more precise results compared to the unconstrained case. This has been verified based on the results obtained from the covariance matrix of the estimated parameters and the empirical results using kinematic positioning samples as well. The theoretical standard deviations of the horizontal components are reduced by a factor ranging from 1.24 to 2.64. The improvement on the empirical standard deviation of the horizontal components ranges from 1.08 to 2.2.

Free and constrained symplectic integrators for numerical general relativity

International Nuclear Information System (INIS)

Richter, Ronny; Lubich, Christian

2008-01-01

We consider symplectic time integrators in numerical general relativity and discuss both free and constrained evolution schemes. For free evolution of ADM-like equations we propose the use of the Stoermer-Verlet method, a standard symplectic integrator which here is explicit in the computationally expensive curvature terms. For the constrained evolution we give a formulation of the evolution equations that enforces the momentum constraints in a holonomically constrained Hamiltonian system and turns the Hamilton constraint function from a weak to a strong invariant of the system. This formulation permits the use of the constraint-preserving symplectic RATTLE integrator, a constrained version of the Stoermer-Verlet method. The behavior of the methods is illustrated on two effectively (1+1)-dimensional versions of Einstein's equations, which allow us to investigate a perturbed Minkowski problem and the Schwarzschild spacetime. We compare symplectic and non-symplectic integrators for free evolution, showing very different numerical behavior for nearly-conserved quantities in the perturbed Minkowski problem. Further we compare free and constrained evolution, demonstrating in our examples that enforcing the momentum constraints can turn an unstable free evolution into a stable constrained evolution. This is demonstrated in the stabilization of a perturbed Minkowski problem with Dirac gauge, and in the suppression of the propagation of boundary instabilities into the interior of the domain in Schwarzschild spacetime

When to call a linear system nonnegative

NARCIS (Netherlands)

Nieuwenhuis, J.W.

1998-01-01

In this paper we will consider discrete time invariant linear systems that allow for an input-state-output representation with a finite dimensional state space, and that have a finite number of inputs and outputs. The basic issue in this paper is when to call these systems nonnegative. An important

Parameter identifiability of linear dynamical systems

Science.gov (United States)

Glover, K.; Willems, J. C.

1974-01-01

It is assumed that the system matrices of a stationary linear dynamical system were parametrized by a set of unknown parameters. The question considered here is, when can such a set of unknown parameters be identified from the observed data? Conditions for the local identifiability of a parametrization are derived in three situations: (1) when input/output observations are made, (2) when there exists an unknown feedback matrix in the system and (3) when the system is assumed to be driven by white noise and only output observations are made. Also a sufficient condition for global identifiability is derived.

Stability problems for linear hyperbolic systems

International Nuclear Information System (INIS)

Eckhoff, K.S.

1975-05-01

The stability properties for the trivial solution of a general linear hyperbolic system of partial differential equations of the first order are studied. It is shown that results may be obtained by studying the stability properties of certain systems of ordinary differential equations which can be constructed from the hyperbolic system (the so-called transport equations). In some cases the associated stability problem for the transport equations can in fact be shown to be equivalent to the stability problem for the hyperbolic system, but in general the transport equations will only give the necessary conditions for stability. (Auth.)

Linear dynamical quantum systems analysis, synthesis, and control

CERN Document Server

Nurdin, Hendra I

2017-01-01

This monograph provides an in-depth treatment of the class of linear-dynamical quantum systems. The monograph presents a detailed account of the mathematical modeling of these systems using linear algebra and quantum stochastic calculus as the main tools for a treatment that emphasizes a system-theoretic point of view and the control-theoretic formulations of quantum versions of familiar problems from the classical (non-quantum) setting, including estimation and filtering, realization theory, and feedback control. Both measurement-based feedback control (i.e., feedback control by a classical system involving a continuous-time measurement process) and coherent feedback control (i.e., feedback control by another quantum system without the intervention of any measurements in the feedback loop) are treated. Researchers and graduates studying systems and control theory, quantum probability and stochastics or stochastic control whether from backgrounds in mechanical or electrical engineering or applied mathematics ...

Simultaneous Balancing and Model Reduction of Switched Linear Systems

OpenAIRE

Monshizadeh, Nima; Trentelman, Hendrikus; Camlibel, M.K.

2011-01-01

In this paper, first, balanced truncation of linear systems is revisited. Then, simultaneous balancing of multiple linear systems is investigated. Necessary and sufficient conditions are introduced to identify the case where simultaneous balancing is possible. The validity of these conditions is not limited to a certain type of balancing, and they are applicable for different types of balancing corresponding to different equations, like Lyapunov or Riccati equations. The results obtained are ...

Solar photovoltaic water pumping system using a new linear actuator

OpenAIRE

Andrada Gascón, Pedro; Castro, Javier

2007-01-01

In this paper a photovoltaic solar pumping system using a new linear actuator is presented. This linear actuator is a double-sided flat two-phase variable-reluctance linear stepper motor that moves a piston-type water pump with the help of a rope, a pulley and a counterweight. The entire actuator pump ensemble is controlled by a simple electronic unit that manages the electric power generated by a photovoltaic array. The proposed system is suitable for rural communities in developing...

Solution of generalized shifted linear systems with complex symmetric matrices

International Nuclear Information System (INIS)

Sogabe, Tomohiro; Hoshi, Takeo; Zhang, Shao-Liang; Fujiwara, Takeo

2012-01-01

We develop the shifted COCG method [R. Takayama, T. Hoshi, T. Sogabe, S.-L. Zhang, T. Fujiwara, Linear algebraic calculation of Green’s function for large-scale electronic structure theory, Phys. Rev. B 73 (165108) (2006) 1–9] and the shifted WQMR method [T. Sogabe, T. Hoshi, S.-L. Zhang, T. Fujiwara, On a weighted quasi-residual minimization strategy of the QMR method for solving complex symmetric shifted linear systems, Electron. Trans. Numer. Anal. 31 (2008) 126–140] for solving generalized shifted linear systems with complex symmetric matrices that arise from the electronic structure theory. The complex symmetric Lanczos process with a suitable bilinear form plays an important role in the development of the methods. The numerical examples indicate that the methods are highly attractive when the inner linear systems can efficiently be solved.

A methodology for controlling motion and constraint forces in holonomically constrained systems

International Nuclear Information System (INIS)

Sapio, Vincent De; Srinivasa, Narayan

2015-01-01

Holonomic constraints are ubiquitous in multibody systems. We present an approach to effectively address the control of holonomically constrained systems using a novel decomposition of task, constraint, and posture space. In addition to providing a natural approach for motion control in the presence of constraints, this scheme also allows for concurrent specification of desired constraint forces, given sufficient actuation. It does this by exposing both motion coordinates and constraint forces within the control formalism, allowing for substantial flexibility in control synthesis. Implementations are presented based on a partitioning of the constraint forces into controlled and uncontrolled subsets, as well as a specification of implicit conditions on the constraint forces. A number of examples demonstrate the practical efficacy of the approach. Finally, a system-level methodology for constraint management during robot interactions with the environment is presented

A methodology for controlling motion and constraint forces in holonomically constrained systems

Energy Technology Data Exchange (ETDEWEB)

Sapio, Vincent De, E-mail: vdesapio@hrl.com; Srinivasa, Narayan, E-mail: nsrinivasa@hrl.com [HRL Laboratories, LLC, Information and Systems Sciences Laboratory (United States)

2015-02-15

Holonomic constraints are ubiquitous in multibody systems. We present an approach to effectively address the control of holonomically constrained systems using a novel decomposition of task, constraint, and posture space. In addition to providing a natural approach for motion control in the presence of constraints, this scheme also allows for concurrent specification of desired constraint forces, given sufficient actuation. It does this by exposing both motion coordinates and constraint forces within the control formalism, allowing for substantial flexibility in control synthesis. Implementations are presented based on a partitioning of the constraint forces into controlled and uncontrolled subsets, as well as a specification of implicit conditions on the constraint forces. A number of examples demonstrate the practical efficacy of the approach. Finally, a system-level methodology for constraint management during robot interactions with the environment is presented.

Economic MPC for a linear stochastic system of energy units

DEFF Research Database (Denmark)

Jørgensen, John Bagterp; Sokoler, Leo Emil; Standardi, Laura

2016-01-01

This paper summarizes comprehensively the work in four recent PhD theses from the Technical University of Denmark related to Economic MPC of future power systems. Future power systems will consist of a large number of decentralized power producers and a large number of controllable power consumers...... in addition to stochastic power producers such as wind turbines and solar power plants. Control of such large scale systems requires new control algorithms. In this paper, we formulate the control of such a system as an Economic Model Predictive Control (MPC) problem. When the power producers and controllable...... power consumers have linear dynamics, the Economic MPC may be expressed as a linear program. We provide linear models for a number of energy units in an energy system, formulate an Economic MPC for coordination of such a system. We indicate how advances in computational MPC makes the solutions...

Experimental quantum computing to solve systems of linear equations.

Science.gov (United States)

Cai, X-D; Weedbrook, C; Su, Z-E; Chen, M-C; Gu, Mile; Zhu, M-J; Li, Li; Liu, Nai-Le; Lu, Chao-Yang; Pan, Jian-Wei

2013-06-07

Solving linear systems of equations is ubiquitous in all areas of science and engineering. With rapidly growing data sets, such a task can be intractable for classical computers, as the best known classical algorithms require a time proportional to the number of variables N. A recently proposed quantum algorithm shows that quantum computers could solve linear systems in a time scale of order log(N), giving an exponential speedup over classical computers. Here we realize the simplest instance of this algorithm, solving 2×2 linear equations for various input vectors on a quantum computer. We use four quantum bits and four controlled logic gates to implement every subroutine required, demonstrating the working principle of this algorithm.

Weak regularizability and pole assignment for non-square linear systems

Czech Academy of Sciences Publication Activity Database

Korotka, Tetiana; Loiseau, J. J.; Zagalak, Petr

2012-01-01

Roč. 48, č. 6 (2012), s. 1065-1088 ISSN 0023-5954 R&D Projects: GA ČR GAP103/12/2431 Keywords : linear systems * linear state feedback * pole assignment Subject RIV: BC - Control Systems Theory Impact factor: 0.619, year: 2012 http://library.utia.cas.cz/separaty/2013/AS/korotka-0386325.pdf

Communication Schemes with Constrained Reordering of Resources

DEFF Research Database (Denmark)

Popovski, Petar; Utkovski, Zoran; Trillingsgaard, Kasper Fløe

2013-01-01

This paper introduces a communication model inspired by two practical scenarios. The first scenario is related to the concept of protocol coding, where information is encoded in the actions taken by an existing communication protocol. We investigate strategies for protocol coding via combinatorial...... reordering of the labelled user resources (packets, channels) in an existing, primary system. However, the degrees of freedom of the reordering are constrained by the operation of the primary system. The second scenario is related to communication systems with energy harvesting, where the transmitted signals...... are constrained by the energy that is available through the harvesting process. We have introduced a communication model that covers both scenarios and elicits their key feature, namely the constraints of the primary system or the harvesting process. We have shown how to compute the capacity of the channels...

Experimental Approaches at Linear Colliders

International Nuclear Information System (INIS)

Jaros, John A

2002-01-01

Precision measurements have played a vital role in our understanding of elementary particle physics. Experiments performed using e + e - collisions have contributed an essential part. Recently, the precision measurements at LEP and SLC have probed the standard model at the quantum level and severely constrained the mass of the Higgs boson [1]. Coupled with the limits on the Higgs mass from direct searches [2], this enables the mass to be constrained to be in the range 115-205 GeV. Developments in accelerator R and D have matured to the point where one could contemplate construction of a linear collider with initial energy in the 500 GeV range and a credible upgrade path to ∼ 1 TeV. Now is therefore the correct time to critically evaluate the case for such a facility

Linearly Refined Session Types

Directory of Open Access Journals (Sweden)

Pedro Baltazar

2012-11-01

Full Text Available Session types capture precise protocol structure in concurrent programming, but do not specify properties of the exchanged values beyond their basic type. Refinement types are a form of dependent types that can address this limitation, combining types with logical formulae that may refer to program values and can constrain types using arbitrary predicates. We present a pi calculus with assume and assert operations, typed using a session discipline that incorporates refinement formulae written in a fragment of Multiplicative Linear Logic. Our original combination of session and refinement types, together with the well established benefits of linearity, allows very fine-grained specifications of communication protocols in which refinement formulae are treated as logical resources rather than persistent truths.

Constrained Perturbation Regularization Approach for Signal Estimation Using Random Matrix Theory

KAUST Repository

Suliman, Mohamed Abdalla Elhag

2016-10-06

In this work, we propose a new regularization approach for linear least-squares problems with random matrices. In the proposed constrained perturbation regularization approach, an artificial perturbation matrix with a bounded norm is forced into the system model matrix. This perturbation is introduced to improve the singular-value structure of the model matrix and, hence, the solution of the estimation problem. Relying on the randomness of the model matrix, a number of deterministic equivalents from random matrix theory are applied to derive the near-optimum regularizer that minimizes the mean-squared error of the estimator. Simulation results demonstrate that the proposed approach outperforms a set of benchmark regularization methods for various estimated signal characteristics. In addition, simulations show that our approach is robust in the presence of model uncertainty.

Linear Optimization of Frequency Spectrum Assignments Across System

Science.gov (United States)

2016-03-01

selection tools, frequency allocation, transmission optimization, electromagnetic maneuver warfare, electronic protection, assignment model 15. NUMBER ...Characteristics Modeled ...............................................................29 Table 10.   Antenna Systems Modeled , Number of Systems and...surveillance EW early warning GAMS general algebraic modeling system GHz gigahertz IDE integrated development environment ILP integer linear program

Decentralized linear quadratic power system stabilizers for multi ...

Indian Academy of Sciences (India)

Introduction. Modern excitation systems considerably enhance the overall transient stability of power systems ..... to the local bus rather than the angle δ measured with respect to the remote bus. ... With this in view, the linear and nonlinear per-.

Linear response theory for quantum open systems

OpenAIRE

Wei, J. H.; Yan, YiJing

2011-01-01

Basing on the theory of Feynman's influence functional and its hierarchical equations of motion, we develop a linear response theory for quantum open systems. Our theory provides an effective way to calculate dynamical observables of a quantum open system at its steady-state, which can be applied to various fields of non-equilibrium condensed matter physics.

Lowest-order constrained variational method for simple many-fermion systems

International Nuclear Information System (INIS)

Alexandrov, I.; Moszkowski, S.A.; Wong, C.W.

1975-01-01

The authors study the potential energy of many-fermion systems calculated by the lowest-order constrained variational (LOCV) method of Pandharipande. Two simple two-body interactions are used. For a simple hard-core potential in a dilute Fermi gas, they find that the Huang-Yang exclusion correction can be used to determine a healing distance. The result is close to the older Pandharipande prescription for the healing distance. For a hard core plus attractive exponential potential, the LOCV result agrees closely with the lowest-order separation method of Moszkowski and Scott. They find that the LOCV result has a shallow minimum as a function of the healing distance at the Moszkowski-Scott separation distance. The significance of the absence of a Brueckner dispersion correction in the LOCV result is discussed. (Auth.)

Constraining ALPs with linear and circular polarisation measurements of quasar light

International Nuclear Information System (INIS)

Payez, Alexandre

2013-09-01

We discuss the constraints derived on the mixing of photons with light pseudoscalars using the distributions of good-quality linear and circular polarisation measurements of light from the least polarised classes of quasars. We also provide the dependence of our limit on the average electron density in the local supercluster for nearly massless particles.

Constraining ALPs with linear and circular polarisation measurements of quasar light

Energy Technology Data Exchange (ETDEWEB)

Payez, Alexandre [Liege Univ. (Belgium). IFPA Group; Deutsches Elektronen-Synchrotron (DESY), Hamburg (Germany)

2013-09-15

We discuss the constraints derived on the mixing of photons with light pseudoscalars using the distributions of good-quality linear and circular polarisation measurements of light from the least polarised classes of quasars. We also provide the dependence of our limit on the average electron density in the local supercluster for nearly massless particles.

A general digital computer procedure for synthesizing linear automatic control systems

International Nuclear Information System (INIS)

Cummins, J.D.

1961-10-01

The fundamental concepts required for synthesizing a linear automatic control system are considered. A generalized procedure for synthesizing automatic control systems is demonstrated. This procedure has been programmed for the Ferranti Mercury and the IBM 7090 computers. Details of the programmes are given. The procedure uses the linearized set of equations which describe the plant to be controlled as the starting point. Subsequent computations determine the transfer functions between any desired variables. The programmes also compute the root and phase loci for any linear (and some non-linear) configurations in the complex plane, the open loop and closed loop frequency responses of a system, the residues of a function of the complex variable 's' and the time response corresponding to these residues. With these general programmes available the design of 'one point' automatic control systems becomes a routine scientific procedure. Also dynamic assessments of plant may be carried out. Certain classes of multipoint automatic control problems may also be solved with these procedures. Autonomous systems, invariant systems and orthogonal systems may also be studied. (author)

A new active absorption system and its performance to linear and non-linear waves

DEFF Research Database (Denmark)

Andersen, Thomas Lykke; Clavero, M.; Frigaard, Peter Bak

2016-01-01

Highlights •An active absorption system for wavemakers has been developed. •The theory for flush mounted gauges has been extended to cover also small gaps. •The new system has been validated in a wave flume with wavemakers in both ends. •A generation and absorption procedure for highly non-linear...

On stabilization of linear systems with stochastic disturbances and input saturation

NARCIS (Netherlands)

Stoorvogel, A.A.; Weiland, S.; Saberi, A.

2004-01-01

It is well-known that for linear systems internal asymptotic stability implies external stability in the sense that when the external input is in Lp then also the state will be in Lp. However, for the control of linear systems with saturation where the controlled system is nonlinear this implication

Pareto optimality in infinite horizon linear quadratic differential games

NARCIS (Netherlands)

Reddy, P.V.; Engwerda, J.C.

2013-01-01

In this article we derive conditions for the existence of Pareto optimal solutions for linear quadratic infinite horizon cooperative differential games. First, we present a necessary and sufficient characterization for Pareto optimality which translates to solving a set of constrained optimal

Linear programming foundations and extensions

CERN Document Server

Vanderbei, Robert J

2001-01-01

Linear Programming: Foundations and Extensions is an introduction to the field of optimization. The book emphasizes constrained optimization, beginning with a substantial treatment of linear programming, and proceeding to convex analysis, network flows, integer programming, quadratic programming, and convex optimization. The book is carefully written. Specific examples and concrete algorithms precede more abstract topics. Topics are clearly developed with a large number of numerical examples worked out in detail. Moreover, Linear Programming: Foundations and Extensions underscores the purpose of optimization: to solve practical problems on a computer. Accordingly, the book is coordinated with free efficient C programs that implement the major algorithms studied: -The two-phase simplex method; -The primal-dual simplex method; -The path-following interior-point method; -The homogeneous self-dual methods. In addition, there are online JAVA applets that illustrate various pivot rules and variants of the simplex m...

Canonical symmetry of a constrained Hamiltonian system and canonical Ward identity

International Nuclear Information System (INIS)

Li, Zi-ping

1995-01-01

An algorithm for the construction of the generators of the gauge transformation of a constrained Hamiltonian system is given. The relationships among the coefficients connecting the first constraints in the generator are made clear. Starting from the phase space generating function of the Green function, the Ward identity in canonical formalism is deduced. We point out that the quantum equations of motion in canonical form for a system with singular Lagrangian differ from the classical ones whether Dirac's conjecture holds true or not. Applications of the present formulation to the Abelian and non-Abelian gauge theories are given. The expressions for PCAC and generalized PCAC of the AVV vertex are derived exactly from another point of view. A new form of the Ward identity for gauge-ghost proper vertices is obtained which differs from the usual Ward-Takahashi identity arising from the BRS invariance

Evolutionary constrained optimization

CERN Document Server

Deb, Kalyanmoy

2015-01-01

This book makes available a self-contained collection of modern research addressing the general constrained optimization problems using evolutionary algorithms. Broadly the topics covered include constraint handling for single and multi-objective optimizations; penalty function based methodology; multi-objective based methodology; new constraint handling mechanism; hybrid methodology; scaling issues in constrained optimization; design of scalable test problems; parameter adaptation in constrained optimization; handling of integer, discrete and mix variables in addition to continuous variables; application of constraint handling techniques to real-world problems; and constrained optimization in dynamic environment. There is also a separate chapter on hybrid optimization, which is gaining lots of popularity nowadays due to its capability of bridging the gap between evolutionary and classical optimization. The material in the book is useful to researchers, novice, and experts alike. The book will also be useful...

Universal Linear Precoding for NBI-Proof Widely Linear Equalization in MC Systems

Directory of Open Access Journals (Sweden)

Donatella Darsena

2007-09-01

Full Text Available In multicarrier (MC systems, transmitter redundancy, which is introduced by means of finite-impulse response (FIR linear precoders, allows for perfect or zero-forcing (ZF equalization of FIR channels (in the absence of noise. Recently, it has been shown that the noncircular or improper nature of some symbol constellations offers an intrinsic source of redundancy, which can be exploited to design efficient FIR widely-linear (WL receiving structures for MC systems operating in the presence of narrowband interference (NBI. With regard to both cyclic-prefixed and zero-padded transmission techniques, it is shown in this paper that, with appropriately designed precoders, it is possible to synthesize in both cases WL-ZF universal equalizers, which guarantee perfect symbol recovery for any FIR channel. Furthermore, it is theoretically shown that the intrinsic redundancy of the improper symbol sequence also enables WL-ZF equalization, based on the minimum mean output-energy criterion, with improved NBI suppression capabilities. Finally, results of numerical simulations are presented, which assess the merits of the proposed precoding designs and validate the theoretical analysis carried out.

ON THE STABILIZATION OF THE LINEAR HYBRID SYSTEM STRUCTURE

Directory of Open Access Journals (Sweden)

Kirillov

2014-11-01

Full Text Available The linear control hybrid system, consisting of a fi- nite set of subsystems (modes having different dimensions, is considered. The moments of reset time are determined by some complementary function – evolutionary time. This function satisfies the special complementary ordinary differential equation. The mode stabilization problem is solved for some class of piecewise linear controls. The method of stabilization relies on the set of invariant planes, the existence of which is due to the special form of the hybrid system.

Event-Triggered Output-Feedback Control for Disturbed Linear Systems

Directory of Open Access Journals (Sweden)

Hao Jiang

2018-01-01

Full Text Available In the last few decades, event-triggered control received considerable attention, because of advantages in reducing the resource utilization, such as communication load and processor. In this paper, we propose an event-triggered output-feedback controller for disturbed linear systems, in order to achieve both better resource utilization and disturbance attenuation properties at the same time. Based on our prior work on state-feedback H∞ control for disturbed systems, we propose an approach to design an output-feedback H∞ controller for the system whose states are not completely observable, and a sufficient condition guaranteeing the asymptotic stability and robustness of the system is given in the form of LMIs (Linear Matrix Inequalities.

INPUT-OUTPUT STRUCTURE OF LINEAR-DIFFERENTIAL ALGEBRAIC SYSTEMS

NARCIS (Netherlands)

KUIJPER, M; SCHUMACHER, JM

Systems of linear differential and algebraic equations occur in various ways, for instance, as a result of automated modeling procedures and in problems involving algebraic constraints, such as zero dynamics and exact model matching. Differential/algebraic systems may represent an input-output

New approach to solve fully fuzzy system of linear equations using ...

Indian Academy of Sciences (India)

This paper proposes two new methods to solve fully fuzzy system of linear equations. The fuzzy system has been converted to a crisp system of linear equations by using single and double parametric form of fuzzy numbers to obtain the non-negative solution. Double parametric form of fuzzy numbers is defined and applied ...

Continuation of Sets of Constrained Orbit Segments

DEFF Research Database (Denmark)

Schilder, Frank; Brøns, Morten; Chamoun, George Chaouki

Sets of constrained orbit segments of time continuous flows are collections of trajectories that represent a whole or parts of an invariant set. A non-trivial but simple example is a homoclinic orbit. A typical representation of this set consists of an equilibrium point of the flow and a trajectory...... that starts close and returns close to this fixed point within finite time. More complicated examples are hybrid periodic orbits of piecewise smooth systems or quasi-periodic invariant tori. Even though it is possible to define generalised two-point boundary value problems for computing sets of constrained...... orbit segments, this is very disadvantageous in practice. In this talk we will present an algorithm that allows the efficient continuation of sets of constrained orbit segments together with the solution of the full variational problem....

Linear circuits, systems and signal processing: theory and application

International Nuclear Information System (INIS)

Byrnes, C.I.; Saeks, R.E.; Martin, C.F.

1988-01-01

In part because of its universal role as a first approximation of more complicated behaviour and in part because of the depth and breadth of its principle paradigms, the study of linear systems continues to play a central role in control theory and its applications. Enhancing more traditional applications to aerospace and electronics, application areas such as econometrics, finance, and speech and signal processing have contributed to a renaissance in areas such as realization theory and classical automatic feedback control. Thus, the last few years have witnessed a remarkable research effort expended in understanding both new algorithms and new paradigms for modeling and realization of linear processes and in the analysis and design of robust control strategies. The papers in this volume reflect these trends in both the theory and applications of linear systems and were selected from the invited and contributed papers presented at the 8th International Symposium on the Mathematical Theory of Networks and Systems held in Phoenix on June 15-19, 1987

Automatic frequency control system for driving a linear accelerator

International Nuclear Information System (INIS)

Helgesson, A.L.

1976-01-01

An automatic frequency control system is described for maintaining the drive frequency applied to a linear accelerator to produce maximum particle output from the accelerator. The particle output amplitude is measured and the frequency of the radio frequency source powering the linear accelerator is adjusted to maximize particle output amplitude

Stability and Linear Quadratic Differential Games of Discrete-Time Markovian Jump Linear Systems with State-Dependent Noise

Directory of Open Access Journals (Sweden)

Huiying Sun

2014-01-01

Full Text Available We mainly consider the stability of discrete-time Markovian jump linear systems with state-dependent noise as well as its linear quadratic (LQ differential games. A necessary and sufficient condition involved with the connection between stochastic Tn-stability of Markovian jump linear systems with state-dependent noise and Lyapunov equation is proposed. And using the theory of stochastic Tn-stability, we give the optimal strategies and the optimal cost values for infinite horizon LQ stochastic differential games. It is demonstrated that the solutions of infinite horizon LQ stochastic differential games are concerned with four coupled generalized algebraic Riccati equations (GAREs. Finally, an iterative algorithm is presented to solve the four coupled GAREs and a simulation example is given to illustrate the effectiveness of it.

Modeling and analysis of linear hyperbolic systems of balance laws

CERN Document Server

Bartecki, Krzysztof

2016-01-01

This monograph focuses on the mathematical modeling of distributed parameter systems in which mass/energy transport or wave propagation phenomena occur and which are described by partial differential equations of hyperbolic type. The case of linear (or linearized) 2 x 2 hyperbolic systems of balance laws is considered, i.e., systems described by two coupled linear partial differential equations with two variables representing physical quantities, depending on both time and one-dimensional spatial variable. Based on practical examples of a double-pipe heat exchanger and a transportation pipeline, two typical configurations of boundary input signals are analyzed: collocated, wherein both signals affect the system at the same spatial point, and anti-collocated, in which the input signals are applied to the two different end points of the system. The results of this book emerge from the practical experience of the author gained during his studies conducted in the experimental installation of a heat exchange cente...

On differential operators generating iterative systems of linear ODEs of maximal symmetry algebra

Science.gov (United States)

Ndogmo, J. C.

2017-06-01

Although every iterative scalar linear ordinary differential equation is of maximal symmetry algebra, the situation is different and far more complex for systems of linear ordinary differential equations, and an iterative system of linear equations need not be of maximal symmetry algebra. We illustrate these facts by examples and derive families of vector differential operators whose iterations are all linear systems of equations of maximal symmetry algebra. Some consequences of these results are also discussed.

FAST Modularization Framework for Wind Turbine Simulation: Full-System Linearization: Preprint

Energy Technology Data Exchange (ETDEWEB)

Jonkman, Jason; Jonkman, Bonnie

2016-11-01

The wind engineering community relies on multiphysics engineering software to run nonlinear time-domain simulations e.g. for design-standards-based loads analysis. Although most physics involved in wind energy are nonlinear, linearization of the underlying nonlinear system equations is often advantageous to understand the system response and exploit well-established methods and tools for analyzing linear systems. This paper presents the development and verification of the new linearization functionality of the open-source engineering tool FAST v8 for land-based wind turbines, as well as the concepts and mathematical background needed to understand and apply it correctly.

Collimation systems in the next linear collider

International Nuclear Information System (INIS)

Merminga, N.; Irwin, J.; Helm, R.; Ruth, R.D.

1991-02-01

Experience indicates that beam collimation will be an essential element of the next generation e + E - linear colliders. A proposal for using nonlinear lenses to drive beam tails to large amplitudes was presented in a previous paper. Here we study the optimization of such systems including effects of wakefields and optical aberrations. Protection and design of the scrapers in these systems are discussed. 9 refs., 7 figs

Linearly and nonlinearly bidirectionally coupled synchronization of hyperchaotic systems

International Nuclear Information System (INIS)

Zhou Jin; Lu Junan; Wu Xiaoqun

2007-01-01

To date, there have been many results about unidirectionally coupled synchronization of chaotic systems. However, much less work is reported on bidirectionally-coupled synchronization. In this paper, we investigate the synchronization of two bidirectionally coupled Chen hyperchaotic systems, which are coupled linearly and nonlinearly respectively. Firstly, linearly coupled synchronization of two hyperchaotic Chen systems is investigated, and a theorem on how to choose the coupling coefficients are developed to guarantee the global asymptotical synchronization of two coupled hyperchaotic systems. Analysis shows that the choice of the coupling coefficients relies on the bound of the chaotic system. Secondly, the nonlinearly coupled synchronization is studied; a sufficient condition for the locally asymptotical synchronization is derived, which is independent of the bound of the hyperchaotic system. Finally, numerical simulations are included to verify the effectiveness and feasibility of the developed theorems

Stability Tests of Positive Fractional Continuous-time Linear Systems with Delays

Directory of Open Access Journals (Sweden)

Tadeusz Kaczorek

2013-06-01

Full Text Available Necessary and sufficient conditions for the asymptotic stability of positive fractional continuous-time linear systems with many delays are established. It is shown that: 1 the asymptotic stability of the positive fractional system is independent of their delays, 2 the checking of the asymptotic stability of the positive fractional systems with delays can be reduced to checking of the asymptotic stability of positive standard linear systems without delays.

Time and Memory Efficient Online Piecewise Linear Approximation of Sensor Signals.

Science.gov (United States)

Grützmacher, Florian; Beichler, Benjamin; Hein, Albert; Kirste, Thomas; Haubelt, Christian

2018-05-23

Piecewise linear approximation of sensor signals is a well-known technique in the fields of Data Mining and Activity Recognition. In this context, several algorithms have been developed, some of them with the purpose to be performed on resource constrained microcontroller architectures of wireless sensor nodes. While microcontrollers are usually constrained in computational power and memory resources, all state-of-the-art piecewise linear approximation techniques either need to buffer sensor data or have an execution time depending on the segment’s length. In the paper at hand, we propose a novel piecewise linear approximation algorithm, with a constant computational complexity as well as a constant memory complexity. Our proposed algorithm’s worst-case execution time is one to three orders of magnitude smaller and its average execution time is three to seventy times smaller compared to the state-of-the-art Piecewise Linear Approximation (PLA) algorithms in our experiments. In our evaluations, we show that our algorithm is time and memory efficient without sacrificing the approximation quality compared to other state-of-the-art piecewise linear approximation techniques, while providing a maximum error guarantee per segment, a small parameter space of only one parameter, and a maximum latency of one sample period plus its worst-case execution time.

Iterative solution of linear systems in the 20­th century

NARCIS (Netherlands)

Saad, Y.; Vorst, H.A. van der

2000-01-01

This paper sketches the main research developments in the area of iterative methods for solving linear systems during the 20th century. Although iterative methods for solving linear systems find their origin in the early nineteenth century (work by Gauss), the field has seen an explosion of

Linear systems optimal and robust control

CERN Document Server

Sinha, Alok

2007-01-01

Introduction Overview Contents of the Book State Space Description of a Linear System Transfer Function of a Single Input/Single Output (SISO) System State Space Realizations of a SISO System SISO Transfer Function from a State Space Realization Solution of State Space Equations Observability and Controllability of a SISO System Some Important Similarity Transformations Simultaneous Controllability and Observability Multiinput/Multioutput (MIMO) Systems State Space Realizations of a Transfer Function Matrix Controllability and Observability of a MIMO System Matrix-Fraction Description (MFD) MFD of a Transfer Function Matrix for the Minimal Order of a State Space Realization Controller Form Realization from a Right MFD Poles and Zeros of a MIMO Transfer Function Matrix Stability Analysis State Feedback Control and Optimization State Variable Feedback for a Single Input System Computation of State Feedback Gain Matrix for a Multiinput System State Feedback Gain Matrix for a Multi...

Uzawa method for fuzzy linear system

OpenAIRE

Ke Wang

2013-01-01

An Uzawa method is presented for solving fuzzy linear systems whose coefficient matrix is crisp and the right-hand side column is arbitrary fuzzy number vector. The explicit iterative scheme is given. The convergence is analyzed with convergence theorems and the optimal parameter is obtained. Numerical examples are given to illustrate the procedure and show the effectiveness and efficiency of the method.

Control of Non-linear Marine Cooling System

DEFF Research Database (Denmark)

Hansen, Michael; Stoustrup, Jakob; Bendtsen, Jan Dimon

2011-01-01

We consider the problem of designing control laws for a marine cooling system used for cooling the main engine and auxiliary components aboard several classes of container vessels. We focus on achieving simple set point control for the system and do not consider compensation of the non-linearitie......-linearities, closed circuit flow dynamics or transport delays that are present in the system. Control laws are therefore designed using classical control theory and the performance of the design is illustrated through two simulation examples....

Inexact nonlinear improved fuzzy chance-constrained programming model for irrigation water management under uncertainty

Science.gov (United States)

Zhang, Chenglong; Zhang, Fan; Guo, Shanshan; Liu, Xiao; Guo, Ping

2018-01-01

An inexact nonlinear mλ-measure fuzzy chance-constrained programming (INMFCCP) model is developed for irrigation water allocation under uncertainty. Techniques of inexact quadratic programming (IQP), mλ-measure, and fuzzy chance-constrained programming (FCCP) are integrated into a general optimization framework. The INMFCCP model can deal with not only nonlinearities in the objective function, but also uncertainties presented as discrete intervals in the objective function, variables and left-hand side constraints and fuzziness in the right-hand side constraints. Moreover, this model improves upon the conventional fuzzy chance-constrained programming by introducing a linear combination of possibility measure and necessity measure with varying preference parameters. To demonstrate its applicability, the model is then applied to a case study in the middle reaches of Heihe River Basin, northwest China. An interval regression analysis method is used to obtain interval crop water production functions in the whole growth period under uncertainty. Therefore, more flexible solutions can be generated for optimal irrigation water allocation. The variation of results can be examined by giving different confidence levels and preference parameters. Besides, it can reflect interrelationships among system benefits, preference parameters, confidence levels and the corresponding risk levels. Comparison between interval crop water production functions and deterministic ones based on the developed INMFCCP model indicates that the former is capable of reflecting more complexities and uncertainties in practical application. These results can provide more reliable scientific basis for supporting irrigation water management in arid areas.

Structured Control of Affine Linear Parameter Varying Systems

DEFF Research Database (Denmark)

Adegas, Fabiano Daher; Stoustrup, Jakob

2011-01-01

This paper presents a new procedure to design structured controllers for discrete-time affine linear parametervarying systems (A LPV). The class of control structures includes decentralized of any order, fixed order output feedback, simultaneous plant-control design, among others. A parametervarying...... non-convex condition for an upper bound on the induced L2-norm performance is solved by an iterative linear matrix inequalities (LMI) optimization algorithm. Numerical examples demostrate the effectiveness of the proposed approach....

Phase and amplitude detection system for the Stanford Linear Accelerator

International Nuclear Information System (INIS)

Fox, J.D.; Schwarz, H.D.

1983-01-01

A computer controlled phase and amplitude detection system to measure and stabilize the rf power sources in the Stanford Linear Accelerator is described. This system measures the instantaneous phase and amplitude of a 1 microsecond 2856 MHz rf pulse and will be used for phase feedback control and for amplitude and phase jitter detection. This paper discusses the measurement system performance requirements for the operation of the Stanford Linear Collider, and the design and implementation of the phase and amplitude detection system. The fundamental software algorithms used in the measurement are described, as is the performance of the prototype phase and amplitude detector system

Feedback-linearization and feedback-feedforward decentralized control for multimachine power system

Energy Technology Data Exchange (ETDEWEB)

De Tuglie, Enrico [Dipartimento di Ingegneria dell' Ambiente, e per lo Sviluppo Sostenibile - DIASS, Politecnico di Bari, Viale del Turismo 8, 74100 Taranto (Italy); Iannone, Silvio Marcello; Torelli, Francesco [Dipartimento di Elettrotecnica, ed Elettronica - DEE, Politecnico di Bari, Via Re David 200, 70125 Bari (Italy)

2008-03-15

In this paper a decentralized nonlinear controller for large-scale power systems is investigated. The proposed controller design is based on the input-output feedback linearization methodology. In order to overcome computational difficulties in adopting such methodology, the overall interconnected nonlinear system, given as n-order, is analyzed as a cascade connection of an n{sub 1}-order nonlinear subsystem and an n{sub 2}-order linear subsystem. The controller design is obtained by applying input-output feedback linearization to the nonlinear subsystem and adopting a tracking control scheme, based on feedback-feedforward technique, for the linear subsystem. In the assumed system model, which is characterised by an interconnected structure between generating units, a decentralised adaptive controller is implemented by decentralizing these constraints. The use of a totally decentralised controller implies a system performance decay with respect to performance when the system is equipped with a centralised controller. Fortunately, the robustness of the proposed controller, based on input-output feedback procedure, guarantees good performance in terms of disturbance even when disturbances are caused by decentralization of interconnection constraints. Test results, provided on the IEEE 30 bus test system, demonstrate the effectiveness and practical applicability of proposed methodology. (author)

Koopman Invariant Subspaces and Finite Linear Representations of Nonlinear Dynamical Systems for Control.

Science.gov (United States)

Brunton, Steven L; Brunton, Bingni W; Proctor, Joshua L; Kutz, J Nathan

2016-01-01

In this wIn this work, we explore finite-dimensional linear representations of nonlinear dynamical systems by restricting the Koopman operator to an invariant subspace spanned by specially chosen observable functions. The Koopman operator is an infinite-dimensional linear operator that evolves functions of the state of a dynamical system. Dominant terms in the Koopman expansion are typically computed using dynamic mode decomposition (DMD). DMD uses linear measurements of the state variables, and it has recently been shown that this may be too restrictive for nonlinear systems. Choosing the right nonlinear observable functions to form an invariant subspace where it is possible to obtain linear reduced-order models, especially those that are useful for control, is an open challenge. Here, we investigate the choice of observable functions for Koopman analysis that enable the use of optimal linear control techniques on nonlinear problems. First, to include a cost on the state of the system, as in linear quadratic regulator (LQR) control, it is helpful to include these states in the observable subspace, as in DMD. However, we find that this is only possible when there is a single isolated fixed point, as systems with multiple fixed points or more complicated attractors are not globally topologically conjugate to a finite-dimensional linear system, and cannot be represented by a finite-dimensional linear Koopman subspace that includes the state. We then present a data-driven strategy to identify relevant observable functions for Koopman analysis by leveraging a new algorithm to determine relevant terms in a dynamical system by ℓ1-regularized regression of the data in a nonlinear function space; we also show how this algorithm is related to DMD. Finally, we demonstrate the usefulness of nonlinear observable subspaces in the design of Koopman operator optimal control laws for fully nonlinear systems using techniques from linear optimal control.ork, we explore finite

Constrained control allocation for systems with redundant control effectors

OpenAIRE

Bordignon, Kenneth A

1996-01-01

Control allocation is examined for linear time-invariant problems that have more controls than degrees of freedom. The controls are part of a physical system and are subject to limits on their maximum positions. A control allocation scheme commands control deflections in response to some desired output. The ability of a control allocation scheme to produce the desired output without violating the physical position constraints is used to compare allocation schemes. Methods are develop...

Linearized Alternating Direction Method of Multipliers for Constrained Nonconvex Regularized Optimization

Science.gov (United States)

2016-11-22

structure of the graph, we replace the ℓ1- norm by the nonconvex Capped -ℓ1 norm , and obtain the Generalized Capped -ℓ1 regularized logistic regression...X. M. Yuan. Linearized augmented lagrangian and alternating direction methods for nuclear norm minimization. Mathematics of Computation, 82(281):301...better approximations of ℓ0- norm theoretically and computationally beyond ℓ1- norm , for example, the compressive sensing (Xiao et al., 2011). The

An Equilibrium Chance-Constrained Multiobjective Programming Model with Birandom Parameters and Its Application to Inventory Problem

Directory of Open Access Journals (Sweden)

Zhimiao Tao

2013-01-01

Full Text Available An equilibrium chance-constrained multiobjective programming model with birandom parameters is proposed. A type of linear model is converted into its crisp equivalent model. Then a birandom simulation technique is developed to tackle the general birandom objective functions and birandom constraints. By embedding the birandom simulation technique, a modified genetic algorithm is designed to solve the equilibrium chance-constrained multiobjective programming model. We apply the proposed model and algorithm to a real-world inventory problem and show the effectiveness of the model and the solution method.

Fractional action-like variational problems in holonomic, non-holonomic and semi-holonomic constrained and dissipative dynamical systems

Energy Technology Data Exchange (ETDEWEB)

EL-Nabulsi, Ahmad Rami [Department of Nuclear and Energy Engineering, Cheju National University, Ara-dong 1, Jeju 690-756 (Korea, Republic of)], E-mail: nabulsiahmadrami@yahoo.fr

2009-10-15

We communicate through this work the fractional calculus of variations and its corresponding Euler-Lagrange equations in 1D constrained holonomic, non-holonomic, and semi-holonomic dissipative dynamical system. The extension of the laws obtained to the 2D space state is done. Some interesting consequences are revealed.

On the discretization of linear fractional representations of LPV systems

NARCIS (Netherlands)

Toth, R.; Lovera, M.; Heuberger, P.S.C.; Corno, M.; Hof, Van den P.M.J.

2012-01-01

Commonly, controllers for linear parameter-varying (LPV) systems are designed in continuous time using a linear fractional representation (LFR) of the plant. However, the resulting controllers are implemented on digital hardware. Furthermore, discrete-time LPV synthesis approaches require a

An algorithm for the solution of dynamic linear programs

Science.gov (United States)

Psiaki, Mark L.

1989-01-01

The algorithm's objective is to efficiently solve Dynamic Linear Programs (DLP) by taking advantage of their special staircase structure. This algorithm constitutes a stepping stone to an improved algorithm for solving Dynamic Quadratic Programs, which, in turn, would make the nonlinear programming method of Successive Quadratic Programs more practical for solving trajectory optimization problems. The ultimate goal is to being trajectory optimization solution speeds into the realm of real-time control. The algorithm exploits the staircase nature of the large constraint matrix of the equality-constrained DLPs encountered when solving inequality-constrained DLPs by an active set approach. A numerically-stable, staircase QL factorization of the staircase constraint matrix is carried out starting from its last rows and columns. The resulting recursion is like the time-varying Riccati equation from multi-stage LQR theory. The resulting factorization increases the efficiency of all of the typical LP solution operations over that of a dense matrix LP code. At the same time numerical stability is ensured. The algorithm also takes advantage of dynamic programming ideas about the cost-to-go by relaxing active pseudo constraints in a backwards sweeping process. This further decreases the cost per update of the LP rank-1 updating procedure, although it may result in more changes of the active set that if pseudo constraints were relaxed in a non-stagewise fashion. The usual stability of closed-loop Linear/Quadratic optimally-controlled systems, if it carries over to strictly linear cost functions, implies that the saving due to reduced factor update effort may outweigh the cost of an increased number of updates. An aerospace example is presented in which a ground-to-ground rocket's distance is maximized. This example demonstrates the applicability of this class of algorithms to aerospace guidance. It also sheds light on the efficacy of the proposed pseudo constraint relaxation

Model-Checking of Linear-Time Properties in Multi-Valued Systems

OpenAIRE

Li, Yongming; Droste, Manfred; Lei, Lihui

2012-01-01

In this paper, we study model-checking of linear-time properties in multi-valued systems. Safety property, invariant property, liveness property, persistence and dual-persistence properties in multi-valued logic systems are introduced. Some algorithms related to the above multi-valued linear-time properties are discussed. The verification of multi-valued regular safety properties and multi-valued $\\omega$-regular properties using lattice-valued automata are thoroughly studied. Since the law o...

Frequency Constrained ShiftCP Modeling of Neuroimaging Data

DEFF Research Database (Denmark)

Mørup, Morten; Hansen, Lars Kai; Madsen, Kristoffer H.

2011-01-01

The shift invariant multi-linear model based on the CandeComp/PARAFAC (CP) model denoted ShiftCP has proven useful for the modeling of latency changes in trial based neuroimaging data[17]. In order to facilitate component interpretation we presently extend the shiftCP model such that the extracted...... components can be constrained to pertain to predefined frequency ranges such as alpha, beta and gamma activity. To infer the number of components in the model we propose to apply automatic relevance determination by imposing priors that define the range of variation of each component of the shiftCP model...

Global stabilization of linear continuous time-varying systems with bounded controls

International Nuclear Information System (INIS)

Phat, V.N.

2004-08-01

This paper deals with the problem of global stabilization of a class of linear continuous time-varying systems with bounded controls. Based on the controllability of the nominal system, a sufficient condition for the global stabilizability is proposed without solving any Riccati differential equation. Moreover, we give sufficient conditions for the robust stabilizability of perturbation/uncertain linear time-varying systems with bounded controls. (author)

Should Pruning be a Pre-Processor of any Linear System?

Science.gov (United States)

Sen, Syamal K.; Ramakrishnan, Suja; Agarwal, Ravi P.; Shaykhian, Gholam Ali

2011-01-01

There are many real-world problems whose mathematical models turn out to be linear systems Ax = b, where A is an m x n matrix. Each equation of the linear system is an information. An information, in a physical problem, such as 4 mangoes, 6 bananas, and 5 oranges cost $10, is mathematically modeled as an equation 4x(sub 1) + 6x(sub 2) + 5x(sub 3) = 10 , where x(sub 1), x(sub 2), x(sub 3) are each cost of one mango, that of one banana, and that of one orange, respectively. All the information put together in a specified context, constitutes the physical problem and need not be all distinct. Some of these could be redundant, which cannot be readily identified by inspection. The resulting mathematical model will thus have equations corresponding to this redundant information and hence are linearly dependent and thus superfluous. Consequently, these equations once identified should be better pruned in the process of solving the system. The benefits are (i) less computation and hence less error and consequently a better quality of solution and (ii) reduced storage requirements. In literature, the pruning concept is not in vogue so far although it is most desirable. It is assumed that at least one information, i.e. one equation is known to be correct and which will be our first equation. In a numerical linear system, the system could be slightly inconsistent or inconsistent of varying degree. If the system is too inconsistent, then we should fall back on to the physical problem (PP), check the correctness of the PP derived from the material universe, modify it, if necessary, and then check the corresponding mathematical model (MM) and correct it. In nature/material universe, inconsistency is completely nonexistent. If the MM becomes inconsistent, it could be due to error introduced by the concerned measuring device and/or due to assumptions made on the PP to obtain an MM which is relatively easily solvable or simply due to human error. No measuring device can usually

A Decomposition Method for Security Constrained Economic Dispatch of a Three-Layer Power System

Science.gov (United States)

Yang, Junfeng; Luo, Zhiqiang; Dong, Cheng; Lai, Xiaowen; Wang, Yang

2018-01-01

This paper proposes a new decomposition method for the security-constrained economic dispatch in a three-layer large-scale power system. The decomposition is realized using two main techniques. The first is to use Ward equivalencing-based network reduction to reduce the number of variables and constraints in the high-layer model without sacrificing accuracy. The second is to develop a price response function to exchange signal information between neighboring layers, which significantly improves the information exchange efficiency of each iteration and results in less iterations and less computational time. The case studies based on the duplicated RTS-79 system demonstrate the effectiveness and robustness of the proposed method.

Applications of equivalent linearization approaches to nonlinear piping systems

International Nuclear Information System (INIS)

Park, Y.; Hofmayer, C.; Chokshi, N.

1997-01-01

The piping systems in nuclear power plants, even with conventional snubber supports, are highly complex nonlinear structures under severe earthquake loadings mainly due to various mechanical gaps in support structures. Some type of nonlinear analysis is necessary to accurately predict the piping responses under earthquake loadings. The application of equivalent linearization approaches (ELA) to seismic analyses of nonlinear piping systems is presented. Two types of ELA's are studied; i.e., one based on the response spectrum method and the other based on the linear random vibration theory. The test results of main steam and feedwater piping systems supported by snubbers and energy absorbers are used to evaluate the numerical accuracy and limitations

Design of offline measuring system for radiation damage effects on linear CCD

International Nuclear Information System (INIS)

Zhang Yong; Tang Benqi; Xiao Zhigang; Wang Zujun; Huang Fang; Huang Shaoyan

2004-01-01

The paper discusses the hardware design of offline measuring system for radiation damage effects on linear CCD. Some credible results were achieved by using this system. The test results indicate that the system is available for the study of the radiation damage effects on linear CCD. (authors)

Koopman Invariant Subspaces and Finite Linear Representations of Nonlinear Dynamical Systems for Control

Science.gov (United States)

Brunton, Steven L.; Brunton, Bingni W.; Proctor, Joshua L.; Kutz, J. Nathan

2016-01-01

In this work, we explore finite-dimensional linear representations of nonlinear dynamical systems by restricting the Koopman operator to an invariant subspace spanned by specially chosen observable functions. The Koopman operator is an infinite-dimensional linear operator that evolves functions of the state of a dynamical system. Dominant terms in the Koopman expansion are typically computed using dynamic mode decomposition (DMD). DMD uses linear measurements of the state variables, and it has recently been shown that this may be too restrictive for nonlinear systems. Choosing the right nonlinear observable functions to form an invariant subspace where it is possible to obtain linear reduced-order models, especially those that are useful for control, is an open challenge. Here, we investigate the choice of observable functions for Koopman analysis that enable the use of optimal linear control techniques on nonlinear problems. First, to include a cost on the state of the system, as in linear quadratic regulator (LQR) control, it is helpful to include these states in the observable subspace, as in DMD. However, we find that this is only possible when there is a single isolated fixed point, as systems with multiple fixed points or more complicated attractors are not globally topologically conjugate to a finite-dimensional linear system, and cannot be represented by a finite-dimensional linear Koopman subspace that includes the state. We then present a data-driven strategy to identify relevant observable functions for Koopman analysis by leveraging a new algorithm to determine relevant terms in a dynamical system by ℓ1-regularized regression of the data in a nonlinear function space; we also show how this algorithm is related to DMD. Finally, we demonstrate the usefulness of nonlinear observable subspaces in the design of Koopman operator optimal control laws for fully nonlinear systems using techniques from linear optimal control. PMID:26919740

Kalman filtering for time-delayed linear systems

Institute of Scientific and Technical Information of China (English)

LU Xiao; WANG Wei

2006-01-01

This paper is to study the linear minimum variance estimation for discrete- time systems. A simple approach to the problem is presented by developing re-organized innovation analysis for the systems with instantaneous and double time-delayed measurements. It is shown that the derived estimator involves solving three different standard Kalman filtering with the same dimension as the original system. The obtained results form the basis for solving some complicated problems such as H∞ fixed-lag smoothing, preview control, H∞ filtering and control with time delays.

A canonical form of the equation of motion of linear dynamical systems

Science.gov (United States)

Kawano, Daniel T.; Salsa, Rubens Goncalves; Ma, Fai; Morzfeld, Matthias

2018-03-01

The equation of motion of a discrete linear system has the form of a second-order ordinary differential equation with three real and square coefficient matrices. It is shown that, for almost all linear systems, such an equation can always be converted by an invertible transformation into a canonical form specified by two diagonal coefficient matrices associated with the generalized acceleration and displacement. This canonical form of the equation of motion is unique up to an equivalence class for non-defective systems. As an important by-product, a damped linear system that possesses three symmetric and positive definite coefficients can always be recast as an undamped and decoupled system.

Daylighting System Based on Novel Design of Linear Fresnel lens

Directory of Open Access Journals (Sweden)

Thanh Tuan Pham

2017-10-01

Full Text Available In this paper, we present a design and optical simulation of a daylighting system using a novel design of linear Fresnel lens, which is constructed based on the conservation of optical path length and edge ray theorem. The linear Fresnel lens can achieve a high uniformity by using a new idea of design in which each groove of the lens distributes sunlight uniformly over the receiver so that the whole lens also uniformly distributes sunlight over the receiver. In this daylighting system, the novel design of linear Fresnel lens significantly improves the uniformity of collector and distributor. Therefore, it can help to improve the performance of the daylighting system. The structure of the linear Fresnel lenses is designed by using Matlab. Then, the structure of lenses is appreciated by ray tracing in LightToolsTM to find out the optimum lens shape. In addition, the simulation is performed by using LightToolsTM to estimate the efficiency of the daylighting system. The results show that the designed collector can achieve the efficiency of ~80% with the tolerance of ~0.60 and the concentration ratio of 340 times, while the designed distributor can reach a high uniformity of >90%.

Linear and nonlinear dynamic systems in financial time series prediction

Directory of Open Access Journals (Sweden)

Salim Lahmiri

2012-10-01

Full Text Available Autoregressive moving average (ARMA process and dynamic neural networks namely the nonlinear autoregressive moving average with exogenous inputs (NARX are compared by evaluating their ability to predict financial time series; for instance the S&P500 returns. Two classes of ARMA are considered. The first one is the standard ARMA model which is a linear static system. The second one uses Kalman filter (KF to estimate and predict ARMA coefficients. This model is a linear dynamic system. The forecasting ability of each system is evaluated by means of mean absolute error (MAE and mean absolute deviation (MAD statistics. Simulation results indicate that the ARMA-KF system performs better than the standard ARMA alone. Thus, introducing dynamics into the ARMA process improves the forecasting accuracy. In addition, the ARMA-KF outperformed the NARX. This result may suggest that the linear component found in the S&P500 return series is more dominant than the nonlinear part. In sum, we conclude that introducing dynamics into the ARMA process provides an effective system for S&P500 time series prediction.

Scheduling Aircraft Landings under Constrained Position Shifting

Science.gov (United States)

Balakrishnan, Hamsa; Chandran, Bala

2006-01-01

Optimal scheduling of airport runway operations can play an important role in improving the safety and efficiency of the National Airspace System (NAS). Methods that compute the optimal landing sequence and landing times of aircraft must accommodate practical issues that affect the implementation of the schedule. One such practical consideration, known as Constrained Position Shifting (CPS), is the restriction that each aircraft must land within a pre-specified number of positions of its place in the First-Come-First-Served (FCFS) sequence. We consider the problem of scheduling landings of aircraft in a CPS environment in order to maximize runway throughput (minimize the completion time of the landing sequence), subject to operational constraints such as FAA-specified minimum inter-arrival spacing restrictions, precedence relationships among aircraft that arise either from airline preferences or air traffic control procedures that prevent overtaking, and time windows (representing possible control actions) during which each aircraft landing can occur. We present a Dynamic Programming-based approach that scales linearly in the number of aircraft, and describe our computational experience with a prototype implementation on realistic data for Denver International Airport.

H-Shaped Multiple Linear Motor Drive Platform Control System Design Based on an Inverse System Method

Directory of Open Access Journals (Sweden)

Caiyan Qin

2017-12-01

Full Text Available Due to its simple mechanical structure and high motion stability, the H-shaped platform has been increasingly widely used in precision measuring, numerical control machining and semiconductor packaging equipment, etc. The H-shaped platform is normally driven by multiple (three permanent magnet synchronous linear motors. The main challenges for H-shaped platform-control include synchronous control between the two linear motors in the Y direction as well as total positioning error of the platform mover, a combination of position deviation in X and Y directions. To deal with the above challenges, this paper proposes a control strategy based on the inverse system method through state feedback and dynamic decoupling of the thrust force. First, mechanical dynamics equations have been deduced through the analysis of system coupling based on the platform structure. Second, the mathematical model of the linear motors and the relevant coordinate transformation between dq-axis currents and ABC-phase currents are analyzed. Third, after the main concept of inverse system method being explained, the inverse system model of the platform control system has been designed after defining relevant system variables. Inverse system model compensates the original nonlinear coupled system into pseudo-linear decoupled linear system, for which typical linear control methods, like PID, can be adopted to control the system. The simulation model of the control system is built in MATLAB/Simulink and the simulation result shows that the designed control system has both small synchronous deviation and small total trajectory tracking error. Furthermore, the control program has been run on NI controller for both fixed-loop-time and free-loop-time modes, and the test result shows that the average loop computation time needed is rather small, which makes it suitable for real industrial applications. Overall, it proves that the proposed new control strategy can be used in

a Continuous-Time Positive Linear System

Directory of Open Access Journals (Sweden)

Kyungsup Kim

2013-01-01

Full Text Available This paper discusses a computational method to construct positive realizations with sparse matrices for continuous-time positive linear systems with multiple complex poles. To construct a positive realization of a continuous-time system, we use a Markov sequence similar to the impulse response sequence that is used in the discrete-time case. The existence of the proposed positive realization can be analyzed with the concept of a polyhedral convex cone. We provide a constructive algorithm to compute positive realizations with sparse matrices of some positive systems under certain conditions. A sufficient condition for the existence of a positive realization, under which the proposed constructive algorithm works well, is analyzed.

Chance-constrained programming approach to natural-gas curtailment decisions

Energy Technology Data Exchange (ETDEWEB)

Guldmann, J M

1981-10-01

This paper presents a modeling methodology for the determination of optimal-curtailment decisions by a gas-distribution utility during a chronic gas-shortage situation. Based on the end-use priority approach, a linear-programming model is formulated, that reallocates the available gas supply among the utility's customers while minimizing fuel switching, unemployment, and utility operating costs. This model is then transformed into a chance-constrained program in order to account for the weather-related variability of the gas requirements. The methodology is applied to the East Ohio Gas Company. 16 references, 2 figures, 3 tables.

Iterative algorithms for large sparse linear systems on parallel computers

Science.gov (United States)

Adams, L. M.

1982-01-01

Algorithms for assembling in parallel the sparse system of linear equations that result from finite difference or finite element discretizations of elliptic partial differential equations, such as those that arise in structural engineering are developed. Parallel linear stationary iterative algorithms and parallel preconditioned conjugate gradient algorithms are developed for solving these systems. In addition, a model for comparing parallel algorithms on array architectures is developed and results of this model for the algorithms are given.

Incremental Closed-loop Identification of Linear Parameter Varying Systems

DEFF Research Database (Denmark)

Bendtsen, Jan Dimon; Trangbæk, Klaus

2011-01-01

, closed-loop system identification is more difficult than open-loop identification. In this paper we prove that the so-called Hansen Scheme, a technique known from linear time-invariant systems theory for transforming closed-loop system identification problems into open-loop-like problems, can be extended...

Stability and response bounds of non-conservative linear systems

DEFF Research Database (Denmark)

Pommer, Christian

2003-01-01

For a linear system of second order differential equations the stability is studied by Lyapunov's direct method. The Lyapunov matrix equation is solved and a sufficient condition for stability is expressed by the system matrices. For a system which satisfies the condition for stability the Lyapunov...

Galerkin projection methods for solving multiple related linear systems

Energy Technology Data Exchange (ETDEWEB)

Chan, T.F.; Ng, M.; Wan, W.L.

1996-12-31

We consider using Galerkin projection methods for solving multiple related linear systems A{sup (i)}x{sup (i)} = b{sup (i)} for 1 {le} i {le} s, where A{sup (i)} and b{sup (i)} are different in general. We start with the special case where A{sup (i)} = A and A is symmetric positive definite. The method generates a Krylov subspace from a set of direction vectors obtained by solving one of the systems, called the seed system, by the CG method and then projects the residuals of other systems orthogonally onto the generated Krylov subspace to get the approximate solutions. The whole process is repeated with another unsolved system as a seed until all the systems are solved. We observe in practice a super-convergence behaviour of the CG process of the seed system when compared with the usual CG process. We also observe that only a small number of restarts is required to solve all the systems if the right-hand sides are close to each other. These two features together make the method particularly effective. In this talk, we give theoretical proof to justify these observations. Furthermore, we combine the advantages of this method and the block CG method and propose a block extension of this single seed method. The above procedure can actually be modified for solving multiple linear systems A{sup (i)}x{sup (i)} = b{sup (i)}, where A{sup (i)} are now different. We can also extend the previous analytical results to this more general case. Applications of this method to multiple related linear systems arising from image restoration and recursive least squares computations are considered as examples.

Operator approach to linear control systems

CERN Document Server

Cheremensky, A

1996-01-01

Within the framework of the optimization problem for linear control systems with quadratic performance index (LQP), the operator approach allows the construction of a systems theory including a number of particular infinite-dimensional optimization problems with hardly visible concreteness. This approach yields interesting interpretations of these problems and more effective feedback design methods. This book is unique in its emphasis on developing methods for solving a sufficiently general LQP. Although this is complex material, the theory developed here is built on transparent and relatively simple principles, and readers with less experience in the field of operator theory will find enough material to give them a good overview of the current state of LQP theory and its applications. Audience: Graduate students and researchers in the fields of mathematical systems theory, operator theory, cybernetics, and control systems.

Construction of exact invariants of time-dependent linear nonholonomic dynamical systems

International Nuclear Information System (INIS)

Fu Jingli; Jimenez, Salvador; Tang Yifa; Vazquez, Luis

2008-01-01

In this work, we build exact dynamical invariants for time-dependent, linear, nonholonomic Hamiltonian systems in two dimensions. Our aim is to obtain an additional insight into the theoretical understanding of generalized Hamilton canonical equations. In particular, we investigate systems represented by a quadratic Hamiltonian subject to linear nonholonomic constraints. We use a Lie algebraic method on the systems to build the invariants. The role and scope of these invariants is pointed out

Construction of exact invariants of time-dependent linear nonholonomic dynamical systems

Energy Technology Data Exchange (ETDEWEB)

Fu Jingli [Institute of Mathematical Physics, Zhejiang Sci-Tech University, Hangzhou 310018 (China)], E-mail: sqfujingli@163.com; Jimenez, Salvador [Departamento de Matematica Aplicada TTII, E.T.S.I. Telecomunicacion, Universidad Politecnica de Madrid, 28040 Madrid (Spain); Tang Yifa [State Key Laboratory of Scientific and Engineering Computing, ICMSEC, Academy of Mathematics and Systems Science, Chinese Academy of Sciences, PO Box 2719, Beijing 100080 (China); Vazquez, Luis [Departamento de Matematica Aplicada Facultad de Informatica, Universidad Complutense de Madrid, 28040 Madrid (Spain); Centro de Astrobiologia (CSIC-INTA), Torrejon de Ardoz, 28850 Madrid (Spain)

2008-03-03

In this work, we build exact dynamical invariants for time-dependent, linear, nonholonomic Hamiltonian systems in two dimensions. Our aim is to obtain an additional insight into the theoretical understanding of generalized Hamilton canonical equations. In particular, we investigate systems represented by a quadratic Hamiltonian subject to linear nonholonomic constraints. We use a Lie algebraic method on the systems to build the invariants. The role and scope of these invariants is pointed out.

H 2 guaranteed cost control of discrete linear systems

Directory of Open Access Journals (Sweden)

Colmenares W.

2000-01-01

Full Text Available This paper presents necessary and sufficient conditions for the existence of a quadratically stabilizing output feedback controller which also assures H 2 guaranteed cost performance on a discrete linear uncertain system where the uncertainty is of the norm bounded type. The conditions are presented as a collection of linear matrix inequalities.The solution, however requires a search over a scalar parameter space.

Linear problems and Baecklund transformations for the Hirota-Ohta system

International Nuclear Information System (INIS)

Adler, V.E.; Postnikov, V.V.

2011-01-01

The auxiliary linear problems are presented for all discretization levels of the Hirota-Ohta system. The structure of these linear problems coincides essentially with the structure of Nonlinear Schroedinger hierarchy. The squared eigenfunction constraints are found which relate Hirota-Ohta and Kulish-Sklyanin vectorial NLS hierarchies.

Computer Based Dose Control System on Linear Accelerator

International Nuclear Information System (INIS)

Taxwim; Djoko-SP; Widi-Setiawan; Agus-Budi Wiyatna

2000-01-01

The accelerator technology has been used for radio therapy. DokterKaryadi Hospital in Semarang use electron or X-ray linear accelerator (Linac)for cancer therapy. One of the control parameter of linear accelerator isdose rate. It is particle current or amount of photon rate to the target. Thecontrol of dose rate in linac have been done by adjusting repetition rate ofanode pulse train of electron source. Presently the control is stillproportional control. To enhance the quality of the control result (minimalstationer error, velocity and stability), the dose control system has beendesigned by using the PID (Proportional Integral Differential) controlalgorithm and the derivation of transfer function of control object.Implementation of PID algorithm control system is done by giving an input ofdose error (the different between output dose and dose rate set point). Theoutput of control system is used for correction of repetition rate set pointfrom pulse train of electron source anode. (author)

A Projected Non-linear Conjugate Gradient Method for Interactive Inverse Kinematics

DEFF Research Database (Denmark)

Engell-Nørregård, Morten; Erleben, Kenny

2009-01-01

Inverse kinematics is the problem of posing an articulated figure to obtain a wanted goal, without regarding inertia and forces. Joint limits are modeled as bounds on individual degrees of freedom, leading to a box-constrained optimization problem. We present A projected Non-linear Conjugate...... Gradient optimization method suitable for box-constrained optimization problems for inverse kinematics. We show application on inverse kinematics positioning of a human figure. Performance is measured and compared to a traditional Jacobian Transpose method. Visual quality of the developed method...

Switching control of linear systems for generating chaos

International Nuclear Information System (INIS)

Liu Xinzhi; Teo, Kok-Lay; Zhang Hongtao; Chen Guanrong

2006-01-01

In this paper, a new switching method is developed, which can be applied to generating different types of chaos or chaos-like dynamics from two or more linear systems. A numerical simulation is given to illustrate the generated chaotic dynamic behavior of the systems with some variable parameters. Finally, a circuit is built to realize various chaotic dynamical behaviors

Consys Linear Control System Design Software Package

International Nuclear Information System (INIS)

Diamantidis, Z.

1987-01-01

This package is created in order to help engineers, researchers, students and all who work on linear control systems. The software includes all time and frequency domain analysises, spectral analysises and networks, active filters and regulators design aids. The programmes are written on Hewlett Packard computer in Basic 4.0

First-order systems of linear partial differential equations: normal forms, canonical systems, transform methods

Directory of Open Access Journals (Sweden)

Heinz Toparkus

2014-04-01

Full Text Available In this paper we consider first-order systems with constant coefficients for two real-valued functions of two real variables. This is both a problem in itself, as well as an alternative view of the classical linear partial differential equations of second order with constant coefficients. The classification of the systems is done using elementary methods of linear algebra. Each type presents its special canonical form in the associated characteristic coordinate system. Then you can formulate initial value problems in appropriate basic areas, and you can try to achieve a solution of these problems by means of transform methods.

Constraining dark sector perturbations I: cosmic shear and CMB lensing

International Nuclear Information System (INIS)

Battye, Richard A.; Moss, Adam; Pearson, Jonathan A.

2015-01-01

We present current and future constraints on equations of state for dark sector perturbations. The equations of state considered are those corresponding to a generalized scalar field model and time-diffeomorphism invariant L(g) theories that are equivalent to models of a relativistic elastic medium and also Lorentz violating massive gravity. We develop a theoretical understanding of the observable impact of these models. In order to constrain these models we use CMB temperature data from Planck, BAO measurements, CMB lensing data from Planck and the South Pole Telescope, and weak galaxy lensing data from CFHTLenS. We find non-trivial exclusions on the range of parameters, although the data remains compatible with w=−1. We gauge how future experiments will help to constrain the parameters. This is done via a likelihood analysis for CMB experiments such as CoRE and PRISM, and tomographic galaxy weak lensing surveys, focussing in on the potential discriminatory power of Euclid on mildly non-linear scales

Constraining dark sector perturbations I: cosmic shear and CMB lensing

Science.gov (United States)

Battye, Richard A.; Moss, Adam; Pearson, Jonathan A.

2015-04-01

We present current and future constraints on equations of state for dark sector perturbations. The equations of state considered are those corresponding to a generalized scalar field model and time-diffeomorphism invariant Script L(g) theories that are equivalent to models of a relativistic elastic medium and also Lorentz violating massive gravity. We develop a theoretical understanding of the observable impact of these models. In order to constrain these models we use CMB temperature data from Planck, BAO measurements, CMB lensing data from Planck and the South Pole Telescope, and weak galaxy lensing data from CFHTLenS. We find non-trivial exclusions on the range of parameters, although the data remains compatible with w=-1. We gauge how future experiments will help to constrain the parameters. This is done via a likelihood analysis for CMB experiments such as CoRE and PRISM, and tomographic galaxy weak lensing surveys, focussing in on the potential discriminatory power of Euclid on mildly non-linear scales.

Chosen interval methods for solving linear interval systems with special type of matrix

Science.gov (United States)

Szyszka, Barbara

2013-10-01

The paper is devoted to chosen direct interval methods for solving linear interval systems with special type of matrix. This kind of matrix: band matrix with a parameter, from finite difference problem is obtained. Such linear systems occur while solving one dimensional wave equation (Partial Differential Equations of hyperbolic type) by using the central difference interval method of the second order. Interval methods are constructed so as the errors of method are enclosed in obtained results, therefore presented linear interval systems contain elements that determining the errors of difference method. The chosen direct algorithms have been applied for solving linear systems because they have no errors of method. All calculations were performed in floating-point interval arithmetic.

Punctuated equilibrium in a non-linear system of action

NARCIS (Netherlands)

J.S. Timmermans (Jos)

2008-01-01

textabstractColeman's equilibrium model of social development, the Linear System of Action, is extended to cover the dynamics of societal transitions. The model implemented has the characteristics of a dissipative system. A variation and selection algorithm favoring the retention of relatively

Periodic solutions of asymptotically linear Hamiltonian systems without twist conditions

Energy Technology Data Exchange (ETDEWEB)

Cheng Rong [Coll. of Mathematics and Physics, Nanjing Univ. of Information Science and Tech., Nanjing (China); Dept. of Mathematics, Southeast Univ., Nanjing (China); Zhang Dongfeng [Dept. of Mathematics, Southeast Univ., Nanjing (China)

2010-05-15

In dynamical system theory, especially in many fields of applications from mechanics, Hamiltonian systems play an important role, since many related equations in mechanics can be written in an Hamiltonian form. In this paper, we study the existence of periodic solutions for a class of Hamiltonian systems. By applying the Galerkin approximation method together with a result of critical point theory, we establish the existence of periodic solutions of asymptotically linear Hamiltonian systems without twist conditions. Twist conditions play crucial roles in the study of periodic solutions for asymptotically linear Hamiltonian systems. The lack of twist conditions brings some difficulty to the study. To the authors' knowledge, very little is known about the case, where twist conditions do not hold. (orig.)

Observability of linear control systems on Lie groups

International Nuclear Information System (INIS)

Ayala, V.; Hacibekiroglu, A.K.

1995-01-01

In this paper, we study the observability problem for a linear control system Σ on a Lie group G. The drift vector field of Σ is an infinitesimal automorphism of G and the control vectors are elements in the Lie algebra of G. We establish algebraic conditions to characterize locally and globally observability for Σ. As in the linear case on R n , these conditions are independent of the control vector. We give an algorithm on the co-tangent bundle of G to calculate the equivalence class of the neutral element. (author). 6 refs

Client's Constraining Factors to Construction Project Management

African Journals Online (AJOL)

factors as a significant system that constrains project management success of public and ... finance for the project and prompt payment for work executed; clients .... consideration of the loading patterns of these variables, the major factor is ...

A conceptual design of Final Focus Systems for linear colliders

International Nuclear Information System (INIS)

Brown, K.L.

1987-06-01

Linear colliders are a relatively recent development in the evolution of particle accelerators. This report discusses some of the approaches that have been considered for the design of Final Focus Systems to demagnify the beam exiting from a linac to the small size suitable for collisions at the interaction point. The system receiving the most attention is the one adopted for the SLAC Linear Collider. However, the theory and optical techniques discussed should be applicable to the design efforts for future machines

Chance Constrained Input Relaxation to Congestion in Stochastic DEA. An Application to Iranian Hospitals.

Science.gov (United States)

Kheirollahi, Hooshang; Matin, Behzad Karami; Mahboubi, Mohammad; Alavijeh, Mehdi Mirzaei

2015-01-01

This article developed an approached model of congestion, based on relaxed combination of inputs, in stochastic data envelopment analysis (SDEA) with chance constrained programming approaches. Classic data envelopment analysis models with deterministic data have been used by many authors to identify congestion and estimate its levels; however, data envelopment analysis with stochastic data were rarely used to identify congestion. This article used chance constrained programming approaches to replace stochastic models with "deterministic equivalents". This substitution leads us to non-linear problems that should be solved. Finally, the proposed method based on relaxed combination of inputs was used to identify congestion input in six Iranian hospital with one input and two outputs in the period of 2009 to 2012.

Robust output feedback H-infinity control and filtering for uncertain linear systems

CERN Document Server

Chang, Xiao-Heng

2014-01-01

"Robust Output Feedback H-infinity Control and Filtering for Uncertain Linear Systems" discusses new and meaningful findings on robust output feedback H-infinity control and filtering for uncertain linear systems, presenting a number of useful and less conservative design results based on the linear matrix inequality (LMI) technique. Though primarily intended for graduate students in control and filtering, the book can also serve as a valuable reference work for researchers wishing to explore the area of robust H-infinity control and filtering of uncertain systems. Dr. Xiao-Heng Chang is a Professor at the College of Engineering, Bohai University, China.

Lag synchronization of chaotic systems with time-delayed linear

Indian Academy of Sciences (India)

In this paper, the lag synchronization of chaotic systems with time-delayed linear terms via impulsive control is investigated. Based on the stability theory of impulsive delayed differential equations, some sufficient conditions are obtained guaranteeing the synchronized behaviours between two delayed chaotic systems.

Frequency Interval Cross Gramians for Linear and Bilinear Systems

DEFF Research Database (Denmark)

Jazlan, Ahmad; Sreeram, Victor; Shaker, Hamid Reza

2017-01-01

In many control engineering problems, it is desired to analyze the systems at particular frequency intervals of interest. This paper focuses on the development of frequency interval cross gramians for both linear and bilinear systems. New generalized Sylvester equations for calculating the freque...

CFORM- LINEAR CONTROL SYSTEM DESIGN AND ANALYSIS: CLOSED FORM SOLUTION AND TRANSIENT RESPONSE OF THE LINEAR DIFFERENTIAL EQUATION

Science.gov (United States)

Jamison, J. W.

1994-01-01

CFORM was developed by the Kennedy Space Center Robotics Lab to assist in linear control system design and analysis using closed form and transient response mechanisms. The program computes the closed form solution and transient response of a linear (constant coefficient) differential equation. CFORM allows a choice of three input functions: the Unit Step (a unit change in displacement); the Ramp function (step velocity); and the Parabolic function (step acceleration). It is only accurate in cases where the differential equation has distinct roots, and does not handle the case for roots at the origin (s=0). Initial conditions must be zero. Differential equations may be input to CFORM in two forms - polynomial and product of factors. In some linear control analyses, it may be more appropriate to use a related program, Linear Control System Design and Analysis (KSC-11376), which uses root locus and frequency response methods. CFORM was written in VAX FORTRAN for a VAX 11/780 under VAX VMS 4.7. It has a central memory requirement of 30K. CFORM was developed in 1987.

Finite Abstractions of Max-Plus-Linear Systems : Theory and Algorithms

NARCIS (Netherlands)

Adzkiya, D.

2014-01-01

Max-Plus-Linear (MPL) systems are a class of discrete-event systems with a continuous state space characterizing the timing of the underlying sequential discrete events. These systems are predisposed to describe the timing synchronization between interleaved processes. MPL systems are employed in

A representation theorem for linear discrete-space systems

Directory of Open Access Journals (Sweden)

Sandberg Irwin W.

1998-01-01

Full Text Available The cornerstone of the theory of discrete-time single-input single-output linear systems is the idea that every such system has an input–output map H that can be represented by a convolution or the familiar generalization of a convolution. This thinking involves an oversight which is corrected in this note by adding an additional term to the representation.

Fuzzy Constrained Predictive Optimal Control of High Speed Train with Actuator Dynamics

Directory of Open Access Journals (Sweden)

Xi Wang

2016-01-01

Full Text Available We investigate the problem of fuzzy constrained predictive optimal control of high speed train considering the effect of actuator dynamics. The dynamics feature of the high speed train is modeled as a cascade of cars connected by flexible couplers, and the formulation is mathematically transformed into a Takagi-Sugeno (T-S fuzzy model. The goal of this study is to design a state feedback control law at each decision step to enhance safety, comfort, and energy efficiency of high speed train subject to safety constraints on the control input. Based on Lyapunov stability theory, the problem of optimizing an upper bound on the cruise control cost function subject to input constraints is reduced to a convex optimization problem involving linear matrix inequalities (LMIs. Furthermore, we analyze the influences of second-order actuator dynamics on the fuzzy constrained predictive controller, which shows risk of potentially deteriorating the overall system. Employing backstepping method, an actuator compensator is proposed to accommodate for the influence of the actuator dynamics. The experimental results show that with the proposed approach high speed train can track the desired speed, the relative coupler displacement between the neighbouring cars is stable at the equilibrium state, and the influence of actuator dynamics is reduced, which demonstrate the validity and effectiveness of the proposed approaches.

An analysis of the electromagnetic field in multi-polar linear induction system

International Nuclear Information System (INIS)

Chervenkova, Todorka; Chervenkov, Atanas

2002-01-01

In this paper a new method for determination of the electromagnetic field vectors in a multi-polar linear induction system (LIS) is described. The analysis of the electromagnetic field has been done by four dimensional electromagnetic potentials in conjunction with theory of the magnetic loops . The electromagnetic field vectors are determined in the Minkovski's space as elements of the Maxwell's tensor. The results obtained are compared with those got from the analysis made by the finite elements method (FEM).With the method represented in this paper one can determine the electromagnetic field vectors in the multi-polar linear induction system using four-dimensional potential. A priority of this method is the obtaining of analytical results for the electromagnetic field vectors. These results are also valid for linear media. The dependencies are valid also at high speeds of movement. The results of the investigated linear induction system are comparable to those got by the finite elements method. The investigations may be continued in the determination of other characteristics such as drag force, levitation force, etc. The method proposed in this paper for an analysis of linear induction system can be used for optimization calculations. (Author)

New results for exponential synchronization of linearly coupled ordinary differential systems

International Nuclear Information System (INIS)

Tong Ping; Chen Shi-Hua

2017-01-01

This paper investigates the exponential synchronization of linearly coupled ordinary differential systems. The intrinsic nonlinear dynamics may not satisfy the QUAD condition or weak-QUAD condition. First, it gives a new method to analyze the exponential synchronization of the systems. Second, two theorems and their corollaries are proposed for the local or global exponential synchronization of the coupled systems. Finally, an application to the linearly coupled Hopfield neural networks and several simulations are provided for verifying the effectiveness of the theoretical results. (paper)

Dynamic Response of Non-Linear Inelsatic Systems to Poisson-Driven Stochastic Excitations

DEFF Research Database (Denmark)

Nielsen, Søren R. K.; Iwankiewicz, R.

of an equivalent linearization techni que and substituting the non-analytical non-linearity in the original system by the cubic form in the pertinent state variables. The response moments are evaluated for the equivalent systems with the help of a generalized Ito's differential rule. The analytical results...

Constraining omega and bias from the Stromlo-APM survey

International Nuclear Information System (INIS)

Loveday, J.

1995-05-01

Galaxy redshift surveys provide a distorted picture of the universe due to the non-Hubble component of galaxy motions. By measuring such distortions in the linear regime one can constrain the quantity β = Ω 0.6 where Ω is the cosmological density parameter and b is the (linear) bias factor for optically-selected galaxies. In this paper we estimate β from the Stromlo-APM redshift survey by comparing the amplitude of the direction-averaged redshift space correlation function to the real space correlation function. We find a 95% confidence upper limit of β = 0.75, with a 'best estimate' of β ∼ 0.48. A bias parameter b ∼ 2 is thus required if Ω ≡ 1. However, higher-order correlations measured from the APM galaxy survey indicate a low value for the bias parameter b ∼ 1, requiring that Q approx-lt 0.6

SLAP, Large Sparse Linear System Solution Package

International Nuclear Information System (INIS)

Greenbaum, A.

1987-01-01

1 - Description of program or function: SLAP is a set of routines for solving large sparse systems of linear equations. One need not store the entire matrix - only the nonzero elements and their row and column numbers. Any nonzero structure is acceptable, so the linear system solver need not be modified when the structure of the matrix changes. Auxiliary storage space is acquired and released within the routines themselves by use of the LRLTRAN POINTER statement. 2 - Method of solution: SLAP contains one direct solver, a band matrix factorization and solution routine, BAND, and several interactive solvers. The iterative routines are as follows: JACOBI, Jacobi iteration; GS, Gauss-Seidel Iteration; ILUIR, incomplete LU decomposition with iterative refinement; DSCG and ICCG, diagonal scaling and incomplete Cholesky decomposition with conjugate gradient iteration (for symmetric positive definite matrices only); DSCGN and ILUGGN, diagonal scaling and incomplete LU decomposition with conjugate gradient interaction on the normal equations; DSBCG and ILUBCG, diagonal scaling and incomplete LU decomposition with bi-conjugate gradient iteration; and DSOMN and ILUOMN, diagonal scaling and incomplete LU decomposition with ORTHOMIN iteration

New approach to solve fully fuzzy system of linear equations using ...

Indian Academy of Sciences (India)

Known example problems are solved to illustrate the efficacy and ... The concept of fuzzy set and fuzzy number were first introduced by Zadeh .... (iii) Fully fuzzy linear systems can be solved by linear programming approach, Gauss elim-.

On the convergence of the dynamic series solution of a constrained ...

African Journals Online (AJOL)

The one dimensional problem of analysing the dynamic behaviour of an elevated water tower with elastic deflection–control device and subjected to a dynamic load was examined in [2]. The constrained elastic system was modeled as a column carrying a concentrated mass at its top and elastically constrained at a point ...

Modeling of non-linear CHP efficiency curves in distributed energy systems

DEFF Research Database (Denmark)

Milan, Christian; Stadler, Michael; Cardoso, Gonçalo

2015-01-01

Distributed energy resources gain an increased importance in commercial and industrial building design. Combined heat and power (CHP) units are considered as one of the key technologies for cost and emission reduction in buildings. In order to make optimal decisions on investment and operation...... for these technologies, detailed system models are needed. These models are often formulated as linear programming problems to keep computational costs and complexity in a reasonable range. However, CHP systems involve variations of the efficiency for large nameplate capacity ranges and in case of part load operation......, which can be even of non-linear nature. Since considering these characteristics would turn the models into non-linear problems, in most cases only constant efficiencies are assumed. This paper proposes possible solutions to address this issue. For a mixed integer linear programming problem two...

Fragment approach to constrained density functional theory calculations using Daubechies wavelets

International Nuclear Information System (INIS)

Ratcliff, Laura E.; Genovese, Luigi; Mohr, Stephan; Deutsch, Thierry

2015-01-01

In a recent paper, we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions are optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper, we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e., without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments

Fragment approach to constrained density functional theory calculations using Daubechies wavelets

Energy Technology Data Exchange (ETDEWEB)

Ratcliff, Laura E., E-mail: lratcliff@anl.gov [Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439 (United States); Université de Grenoble Alpes, CEA, INAC-SP2M, L-Sim, F-38000 Grenoble (France); Genovese, Luigi; Mohr, Stephan; Deutsch, Thierry [Université de Grenoble Alpes, CEA, INAC-SP2M, L-Sim, F-38000 Grenoble (France)

2015-06-21

In a recent paper, we presented a linear scaling Kohn-Sham density functional theory (DFT) code based on Daubechies wavelets, where a minimal set of localized support functions are optimized in situ and therefore adapted to the chemical properties of the molecular system. Thanks to the systematically controllable accuracy of the underlying basis set, this approach is able to provide an optimal contracted basis for a given system: accuracies for ground state energies and atomic forces are of the same quality as an uncontracted, cubic scaling approach. This basis set offers, by construction, a natural subset where the density matrix of the system can be projected. In this paper, we demonstrate the flexibility of this minimal basis formalism in providing a basis set that can be reused as-is, i.e., without reoptimization, for charge-constrained DFT calculations within a fragment approach. Support functions, represented in the underlying wavelet grid, of the template fragments are roto-translated with high numerical precision to the required positions and used as projectors for the charge weight function. We demonstrate the interest of this approach to express highly precise and efficient calculations for preparing diabatic states and for the computational setup of systems in complex environments.

Stability margin of linear systems with parameters described by fuzzy numbers.

Science.gov (United States)

Husek, Petr

2011-10-01

This paper deals with the linear systems with uncertain parameters described by fuzzy numbers. The problem of determining the stability margin of those systems with linear affine dependence of the coefficients of a characteristic polynomial on system parameters is studied. Fuzzy numbers describing the system parameters are allowed to be characterized by arbitrary nonsymmetric membership functions. An elegant solution, graphical in nature, based on generalization of the Tsypkin-Polyak plot is presented. The advantage of the presented approach over the classical robust concept is demonstrated on a control of the Fiat Dedra engine model and a control of the quarter car suspension model.

A Type System for the Vectorial Aspect of the Linear-Algebraic Lambda-Calculus

Directory of Open Access Journals (Sweden)

Pablo Arrighi

2012-07-01

Full Text Available We describe a type system for the linear-algebraic lambda-calculus. The type system accounts for the part of the language emulating linear operators and vectors, i.e. it is able to statically describe the linear combinations of terms resulting from the reduction of programs. This gives rise to an original type theory where types, in the same way as terms, can be superposed into linear combinations. We show that the resulting typed lambda-calculus is strongly normalizing and features a weak subject-reduction.

Econometrics analysis of consumer behaviour: a linear expenditure system applied to energy

International Nuclear Information System (INIS)

Giansante, C.; Ferrari, V.

1996-12-01

In economics literature the expenditure system specification is a well known subject. The problem is to define a coherent representation of consumer behaviour through functional forms easy to calculate. In this work it is used the Stone-Geary Linear Expenditure System and its multi-level decision process version. The Linear Expenditure system is characterized by an easy calculating estimation procedure, and its multi-level specification allows substitution and complementary relations between goods. Moreover, the utility function separability condition on which the Utility Tree Approach is based, justifies to use an estimation procedure in two or more steps. This allows to use an high degree of expenditure categories disaggregation, impossible to reach the Linear Expediture System. The analysis is applied to energy sectors

Identification of Nonlinear Dynamic Systems Possessing Some Non-linearities

Directory of Open Access Journals (Sweden)

Y. N. Pavlov

2015-01-01

Full Text Available The subject of this work is the problem of identification of nonlinear dynamic systems based on the experimental data obtained by applying test signals to the system. The goal is to determinate coefficients of differential equations of systems by experimental frequency hodographs and separate similar, but different, in essence, forces: dissipative forces with the square of the first derivative in the motion equations and dissipative force from the action of dry friction. There was a proposal to use the harmonic linearization method to approximate each of the nonlinearity of "quadratic friction" and "dry friction" by linear friction with the appropriate harmonic linearization coefficient.Assume that a frequency transfer function of the identified system has a known form. Assume as well that there are disturbances while obtaining frequency characteristics of the realworld system. As a result, the points of experimentally obtained hodograph move randomly. Searching for solution of the identification problem was in the hodograph class, specified by the system model, which has the form of the frequency transfer function the same as the form of the frequency transfer function of the system identified. Minimizing a proximity criterion (measure of the experimentally obtained system hodograph and the system hodograph model for all the experimental points described and previously published by one of the authors allowed searching for the unknown coefficients of the frequenc ransfer function of the system model. The paper shows the possibility to identify a nonlinear dynamic system with multiple nonlinearities, obtained on the experimental samples of the frequency system hodograph. The proposed algorithm allows to select the nonlinearity of the type "quadratic friction" and "dry friction", i.e. also in the case where the nonlinearity is dependent on the same dynamic parameter, in particular, on the derivative of the system output value. For the dynamic

Linearity improvement on wide-range log signal of neutron measurement system for HANARO

International Nuclear Information System (INIS)

Kim, Young-Ki; Tuetken, Jeffrey S.

1998-01-01

This paper discusses engineering activities for improving the linearity characteristics of the Log Power signal from the neutron measurement system for HANARO. This neutron measurement system uses a fission chamber based detector which covers 10.3 decade-wide range from 10 -8 % full power(FP) up to 200%FP, The Log Power signal is designed to control the reactor at low power levels where most of the reactor physics tests are carried out. Therefore, the linearity characteristics of the Log Power signal is the major factor for accurate reactor power control. During the commissioning of the neutron measurement system, it was found that the linearity characteristics of the Log Power signal, especially near 10 -2 %FP, were not accurate enough for controlling the reactor during physics testing. Analysis of the system linearity data directly measured with reactor operating determined that the system was not operating per the design characteristics established from previous installations. The linearity data, which were taken as the reactor was increased in power, were sent to manufacturer's engineering group and a follow-up measures based on the analysis were then fed back to the field. Through step by step trouble-shooting activities, which included minor circuit modifications and alignment procedure changes, the linearity characteristics have been successfully improved and now exceed minimum performance requirements. This paper discusses the trouble-shooting techniques applied, the changes in the linearity characteristics, special circumstances in the HANARO application and the final resolution. (author)

Influence of capacity- and time-constrained intermediate storage in two-stage food production systems

DEFF Research Database (Denmark)

Akkerman, Renzo; van Donk, Dirk Pieter; Gaalman, Gerard

2007-01-01

In food processing, two-stage production systems with a batch processor in the first stage and packaging lines in the second stage are common and mostly separated by capacity- and time-constrained intermediate storage. This combination of constraints is common in practice, but the literature hardly...... of systems like this. Contrary to the common sense in operations management, the LPT rule is able to maximize the total production volume per day. Furthermore, we show that adding one tank has considerable effects. Finally, we conclude that the optimal setup frequency for batches in the first stage...... pays any attention to this. In this paper, we show how various capacity and time constraints influence the performance of a specific two-stage system. We study the effects of several basic scheduling and sequencing rules in the presence of these constraints in order to learn the characteristics...

Optimization of a constrained linear monochromator design for neutral atom beams.

Science.gov (United States)

Kaltenbacher, Thomas

2016-04-01

A focused ground state, neutral atom beam, exploiting its de Broglie wavelength by means of atom optics, is used for neutral atom microscopy imaging. Employing Fresnel zone plates as a lens for these beams is a well established microscopy technique. To date, even for favorable beam source conditions a minimal focus spot size of slightly below 1μm was reached. This limitation is essentially given by the intrinsic spectral purity of the beam in combination with the chromatic aberration of the diffraction based zone plate. Therefore, it is important to enhance the monochromaticity of the beam, enabling a higher spatial resolution, preferably below 100nm. We propose to increase the monochromaticity of a neutral atom beam by means of a so-called linear monochromator set-up - a Fresnel zone plate in combination with a pinhole aperture - in order to gain more than one order of magnitude in spatial resolution. This configuration is known in X-ray microscopy and has proven to be useful, but has not been applied to neutral atom beams. The main result of this work is optimal design parameters based on models for this linear monochromator set-up followed by a second zone plate for focusing. The optimization was performed for minimizing the focal spot size and maximizing the centre line intensity at the detector position for an atom beam simultaneously. The results presented in this work are for, but not limited to, a neutral helium atom beam. Copyright © 2016 Elsevier B.V. All rights reserved.

A real-time Java tool chain for resource constrained platforms

DEFF Research Database (Denmark)

Korsholm, Stephan Erbs; Søndergaard, Hans; Ravn, Anders P.

2013-01-01

The Java programming language was originally developed for embedded systems, but the resource requirements of previous and current Java implementations - especially memory consumption - tend to exclude them from being used on a significant class of resource constrained embedded platforms. The con......The Java programming language was originally developed for embedded systems, but the resource requirements of previous and current Java implementations - especially memory consumption - tend to exclude them from being used on a significant class of resource constrained embedded platforms...... by integrating: (1) a lean virtual machine (HVM) without any external dependencies on POSIX-like libraries or other OS functionalities, (2) a hardware abstraction layer, implemented almost entirely in Java through the use of hardware objects, first level interrupt handlers, and native variables, and (3....... An evaluation of the presented solution shows that the miniCDj benchmark gets reduced to a size where it can run on resource constrained platforms....

Stochastic Reformulations of Linear Systems: Algorithms and Convergence Theory

KAUST Repository

Richtarik, Peter; Taká č, Martin

2017-01-01

We develop a family of reformulations of an arbitrary consistent linear system into a stochastic problem. The reformulations are governed by two user-defined parameters: a positive definite matrix defining a norm, and an arbitrary discrete or continuous distribution over random matrices. Our reformulation has several equivalent interpretations, allowing for researchers from various communities to leverage their domain specific insights. In particular, our reformulation can be equivalently seen as a stochastic optimization problem, stochastic linear system, stochastic fixed point problem and a probabilistic intersection problem. We prove sufficient, and necessary and sufficient conditions for the reformulation to be exact. Further, we propose and analyze three stochastic algorithms for solving the reformulated problem---basic, parallel and accelerated methods---with global linear convergence rates. The rates can be interpreted as condition numbers of a matrix which depends on the system matrix and on the reformulation parameters. This gives rise to a new phenomenon which we call stochastic preconditioning, and which refers to the problem of finding parameters (matrix and distribution) leading to a sufficiently small condition number. Our basic method can be equivalently interpreted as stochastic gradient descent, stochastic Newton method, stochastic proximal point method, stochastic fixed point method, and stochastic projection method, with fixed stepsize (relaxation parameter), applied to the reformulations.

Stochastic Reformulations of Linear Systems: Algorithms and Convergence Theory

KAUST Repository

Richtarik, Peter

2017-06-04

We develop a family of reformulations of an arbitrary consistent linear system into a stochastic problem. The reformulations are governed by two user-defined parameters: a positive definite matrix defining a norm, and an arbitrary discrete or continuous distribution over random matrices. Our reformulation has several equivalent interpretations, allowing for researchers from various communities to leverage their domain specific insights. In particular, our reformulation can be equivalently seen as a stochastic optimization problem, stochastic linear system, stochastic fixed point problem and a probabilistic intersection problem. We prove sufficient, and necessary and sufficient conditions for the reformulation to be exact. Further, we propose and analyze three stochastic algorithms for solving the reformulated problem---basic, parallel and accelerated methods---with global linear convergence rates. The rates can be interpreted as condition numbers of a matrix which depends on the system matrix and on the reformulation parameters. This gives rise to a new phenomenon which we call stochastic preconditioning, and which refers to the problem of finding parameters (matrix and distribution) leading to a sufficiently small condition number. Our basic method can be equivalently interpreted as stochastic gradient descent, stochastic Newton method, stochastic proximal point method, stochastic fixed point method, and stochastic projection method, with fixed stepsize (relaxation parameter), applied to the reformulations.

New approach to solve symmetric fully fuzzy linear systems

Indian Academy of Sciences (India)

In this paper, we present a method to solve fully fuzzy linear systems with symmetric coefficient matrix. The symmetric coefficient matrix is decomposed into two systems of equations by using Cholesky method and then a solution can be obtained. Numerical examples are given to illustrate our method.

Evolution of the Argonne Tandem Linear Accelerator System (ATLAS) control system

International Nuclear Information System (INIS)

Power, M.; Munson, F.

2012-01-01

Given that the Argonne Tandem Linear Accelerator System (ATLAS) recently celebrated its 25. anniversary, this paper will explore the past, present, and future of the ATLAS Control System, and how it has evolved along with the accelerator and control system technology. ATLAS as we know it today, originated with a Tandem Van de Graff in the sixties. With the addition of the Booster section in the late seventies, came the first computerized control. ATLAS itself was placed into service on June 25, 1985, and was the world's first superconducting linear accelerator for ions. Since its dedication as a National User Facility, more than a thousand experiments by more than 2,000 users worldwide, have taken advantage of the unique capabilities it provides. Today, ATLAS continues to be a user facility for physicists who study the particles that form the heart of atoms. Its most recent addition, CARIBU (Californium Rare Isotope Breeder Upgrade), creates special beams that feed into ATLAS. ATLAS is similar to a living organism, changing and responding to new technological challenges and research needs. As it continues to evolve, so does the control system: from the original days using a DEC PDP-11/34 computer and two CAMAC crates, to a DEC Alpha computer running Vsystem software and more than twenty CAMAC crates, to distributed computers and VME systems. Future upgrades are also in the planning stages that will continue to evolve the control system. (authors)

Variance-Constrained Robust Estimation for Discrete-Time Systems with Communication Constraints

Directory of Open Access Journals (Sweden)

Baofeng Wang

2014-01-01

Full Text Available This paper is concerned with a new filtering problem in networked control systems (NCSs subject to limited communication capacity, which includes measurement quantization, random transmission delay, and packets loss. The measurements are first quantized via a logarithmic quantizer and then transmitted through a digital communication network with random delay and packet loss. The three communication constraints phenomena which can be seen as a class of uncertainties are formulated by a stochastic parameter uncertainty system. The purpose of the paper is to design a linear filter such that, for all the communication constraints, the error state of the filtering process is mean square bounded and the steady-state variance of the estimation error for each state is not more than the individual prescribed upper bound. It is shown that the desired filtering can effectively be solved if there are positive definite solutions to a couple of algebraic Riccati-like inequalities or linear matrix inequalities. Finally, an illustrative numerical example is presented to demonstrate the effectiveness and flexibility of the proposed design approach.

Structure-constrained sparse canonical correlation analysis with an application to microbiome data analysis.

Science.gov (United States)

Chen, Jun; Bushman, Frederic D; Lewis, James D; Wu, Gary D; Li, Hongzhe

2013-04-01

Motivated by studying the association between nutrient intake and human gut microbiome composition, we developed a method for structure-constrained sparse canonical correlation analysis (ssCCA) in a high-dimensional setting. ssCCA takes into account the phylogenetic relationships among bacteria, which provides important prior knowledge on evolutionary relationships among bacterial taxa. Our ssCCA formulation utilizes a phylogenetic structure-constrained penalty function to impose certain smoothness on the linear coefficients according to the phylogenetic relationships among the taxa. An efficient coordinate descent algorithm is developed for optimization. A human gut microbiome data set is used to illustrate this method. Both simulations and real data applications show that ssCCA performs better than the standard sparse CCA in identifying meaningful variables when there are structures in the data.

Partial Synchronization Manifolds for Linearly Time-Delay Coupled Systems

OpenAIRE

Steur, Erik; van Leeuwen, Cees; Michiels, Wim

2014-01-01

Sometimes a network of dynamical systems shows a form of incomplete synchronization characterized by synchronization of some but not all of its systems. This type of incomplete synchronization is called partial synchronization. Partial synchronization is associated with the existence of partial synchronization manifolds, which are linear invariant subspaces of C, the state space of the network of systems. We focus on partial synchronization manifolds in networks of system...

Genetic design of interpolated non-linear controllers for linear plants

International Nuclear Information System (INIS)

Ajlouni, N.

2000-01-01

The techniques of genetic algorithms are proposed as a means of designing non-linear PID control systems. It is shown that the use of genetic algorithms for this purpose results in highly effective non-linear PID control systems. These results are illustrated by using genetic algorithms to design a non-linear PID control system and contrasting the results with an optimally tuned linear PID controller. (author)

Constrained evolution in numerical relativity

Science.gov (United States)

Anderson, Matthew William

The strongest potential source of gravitational radiation for current and future detectors is the merger of binary black holes. Full numerical simulation of such mergers can provide realistic signal predictions and enhance the probability of detection. Numerical simulation of the Einstein equations, however, is fraught with difficulty. Stability even in static test cases of single black holes has proven elusive. Common to unstable simulations is the growth of constraint violations. This work examines the effect of controlling the growth of constraint violations by solving the constraints periodically during a simulation, an approach called constrained evolution. The effects of constrained evolution are contrasted with the results of unconstrained evolution, evolution where the constraints are not solved during the course of a simulation. Two different formulations of the Einstein equations are examined: the standard ADM formulation and the generalized Frittelli-Reula formulation. In most cases constrained evolution vastly improves the stability of a simulation at minimal computational cost when compared with unconstrained evolution. However, in the more demanding test cases examined, constrained evolution fails to produce simulations with long-term stability in spite of producing improvements in simulation lifetime when compared with unconstrained evolution. Constrained evolution is also examined in conjunction with a wide variety of promising numerical techniques, including mesh refinement and overlapping Cartesian and spherical computational grids. Constrained evolution in boosted black hole spacetimes is investigated using overlapping grids. Constrained evolution proves to be central to the host of innovations required in carrying out such intensive simulations.

Fitting and forecasting coupled dark energy in the non-linear regime

Energy Technology Data Exchange (ETDEWEB)

Casas, Santiago; Amendola, Luca; Pettorino, Valeria; Vollmer, Adrian [Institut für Theoretische Physik, Ruprecht-Karls-Universität Heidelberg, Philosophenweg 16, Heidelberg, 69120 Germany (Germany); Baldi, Marco, E-mail: casas@thphys.uni-heidelberg.de, E-mail: l.amendola@thphys.uni-heidelberg.de, E-mail: mail@marcobaldi.it, E-mail: v.pettorino@thphys.uni-heidelberg.de, E-mail: vollmer@thphys.uni-heidelberg.de [Dipartimento di Fisica e Astronomia, Alma Mater Studiorum Università di Bologna, viale Berti Pichat, 6/2, Bologna, I-40127 Italy (Italy)

2016-01-01

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range 0z=–1.6 and wave modes below 0k=1 h/Mpc. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, β{sup 2}, with precision smaller than 4% and all other cosmological parameters better than 1%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications.

Fitting and forecasting coupled dark energy in the non-linear regime

International Nuclear Information System (INIS)

Casas, Santiago; Amendola, Luca; Pettorino, Valeria; Vollmer, Adrian; Baldi, Marco

2016-01-01

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range 0z=–1.6 and wave modes below 0k=1 h/Mpc. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, β 2 , with precision smaller than 4% and all other cosmological parameters better than 1%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications

Closed-loop Identification for Control of Linear Parameter Varying Systems

DEFF Research Database (Denmark)

Bendtsen, Jan Dimon; Trangbæk, Klaus

2014-01-01

, closed- loop system identification is more difficult than open-loop identification. In this paper we prove that the so-called Hansen Scheme, a technique known from linear time-invariant systems theory for transforming closed-loop system identification problems into open-loop-like problems, can...

A system for monitoring the radiation effects of a proton linear accelerator

Energy Technology Data Exchange (ETDEWEB)

Skorkin, V. M., E-mail: skorkin@inr.ru; Belyanski, K. L.; Skorkin, A. V. [Russian Academy of Sciences, Institute for Nuclear Research (Russian Federation)

2016-12-15

The system for real-time monitoring of radioactivity of a high-current proton linear accelerator detects secondary neutron emission from proton beam losses in transport channels and measures the activity of radionuclides in gas and aerosol emissions and the radiation background in the environment affected by a linear accelerator. The data provided by gamma, beta, and neutron detectors are transferred over a computer network to the central server. The system allows one to monitor proton beam losses, the activity of gas and aerosol emissions, and the radiation emission level of a linear accelerator in operation.

Linear systems surviving the first breakdown; Systemes unbaires survivant a la premiere panne

Energy Technology Data Exchange (ETDEWEB)

Uberschlag, J [Commissariat a l' Energie Atomique, Saclay (France). Centre d' Etudes Nucleaires

1967-07-01

Various types of linear systems are described which are not affected by the first breakdown. They make it possible to operate continuously and are thus very reliable. This is because the first breakdown which occurs affects only very slightly the operation. These components can be replaced during working. The operation, the errors, and the detection are briefly considered in the case of three different designs of linear servo systems. An attempt at comparison is made, it could be developed in a particular case. (author) [French] On decrit divers types de systemes lineaires survivant a la premiere panne. Ils permettent un fonctionnement permanent et sont donc d'une tres grande fiabilite. En effet, une panne, qui peut etre signalee, perturbe peu le fonctionnement. Ces composants peuvent etre remplaces en marche. Les considerations de fonctionnement, d'erreur, de detection des pannes sont succinctement presentees, sur trois schemas de systemes lineaires asservis. Une tentative de comparaison est faite, qui pourrait etre developpee dans des cas precis. (auteur)

Sparse Linear Solver for Power System Analysis Using FPGA

National Research Council Canada - National Science Library

Johnson, J. R; Nagvajara, P; Nwankpa, C

2005-01-01

.... Numerical solution to load flow equations are typically computed using Newton-Raphson iteration, and the most time consuming component of the computation is the solution of a sparse linear system...

An inexact multistage fuzzy-stochastic programming for regional electric power system management constrained by environmental quality.

Science.gov (United States)

Fu, Zhenghui; Wang, Han; Lu, Wentao; Guo, Huaicheng; Li, Wei

2017-12-01

Electric power system involves different fields and disciplines which addressed the economic system, energy system, and environment system. Inner uncertainty of this compound system would be an inevitable problem. Therefore, an inexact multistage fuzzy-stochastic programming (IMFSP) was developed for regional electric power system management constrained by environmental quality. A model which concluded interval-parameter programming, multistage stochastic programming, and fuzzy probability distribution was built to reflect the uncertain information and dynamic variation in the case study, and the scenarios under different credibility degrees were considered. For all scenarios under consideration, corrective actions were allowed to be taken dynamically in accordance with the pre-regulated policies and the uncertainties in reality. The results suggest that the methodology is applicable to handle the uncertainty of regional electric power management systems and help the decision makers to establish an effective development plan.

Reinforcement learning solution for HJB equation arising in constrained optimal control problem.

Science.gov (United States)

Luo, Biao; Wu, Huai-Ning; Huang, Tingwen; Liu, Derong

2015-11-01

The constrained optimal control problem depends on the solution of the complicated Hamilton-Jacobi-Bellman equation (HJBE). In this paper, a data-based off-policy reinforcement learning (RL) method is proposed, which learns the solution of the HJBE and the optimal control policy from real system data. One important feature of the off-policy RL is that its policy evaluation can be realized with data generated by other behavior policies, not necessarily the target policy, which solves the insufficient exploration problem. The convergence of the off-policy RL is proved by demonstrating its equivalence to the successive approximation approach. Its implementation procedure is based on the actor-critic neural networks structure, where the function approximation is conducted with linearly independent basis functions. Subsequently, the convergence of the implementation procedure with function approximation is also proved. Finally, its effectiveness is verified through computer simulations. Copyright © 2015 Elsevier Ltd. All rights reserved.

Linear-array systems for aerospace NDE

International Nuclear Information System (INIS)

Smith, Robert A.; Willsher, Stephen J.; Bending, Jamie M.

1999-01-01

Rapid large-area inspection of composite structures for impact damage and multi-layered aluminum skins for corrosion has been a recognized priority for several years in both military and civil aerospace applications. Approaches to this requirement have followed two clearly different routes: the development of novel large-area inspection systems, and the enhancement of current ultrasonic or eddy-current methods to reduce inspection times. Ultrasonic inspection is possible with standard flaw detection equipment but the addition of a linear ultrasonic array could reduce inspection times considerably. In order to investigate their potential, 9-element and 17-element linear ultrasonic arrays for composites, and 64-element arrays for aluminum skins, have been developed to DERA specifications for use with the ANDSCAN area scanning system. A 5 m 2 composite wing surface has been scanned with a scan resolution of approximately 3 mm in 6 hours. With subsequent software and hardware improvements all four composite wing surfaces (top/bottom, left/right) of a military fighter aircraft can potentially be inspected in less than a day. Array technology has been very widely used in the medical ultrasound field although rarely above 10 MHz, whereas lap-joint inspection requires a pulse center-frequency of 12 to 20 MHz in order to resolve the separate interfaces in the lap joint. A 128 mm-long multi-element array of 5 mmx2 mm ultrasonic elements for use with the ANDSCAN scanning software was produced to a DERA specification by an NDT manufacturer with experience in the medical imaging field. This paper analyses the performance of the transducers that have been produced and evaluates their use in scanning systems of different configurations

Euclidean null controllability of linear systems with delays in state ...

African Journals Online (AJOL)

Sufficient conditions are developed for the Euclidean controllability of linear systems with delay in state and in control. Namely, if the uncontrolled system is uniformly asymptotically stable and the control equation proper, then the control system is Euclidean null controllable. Journal of the Nigerian Association of ...

Comparison of the nuclear code systems LINEAR-RECENT-NJOY and NJOY

International Nuclear Information System (INIS)

Seehusen, J.

1983-07-01

The reconstructed cross sections of the code systems LINEAR-RECENT-GROUPIE (Version 1982) and NJOY (Version 1982) have been compared for several materials. Some fuel cycle isotopes and structural materials of the ENDF/B-4 general purpose and ENDF/B-5 dosimetry files have been choosen. The reconstructed total, capture and fission cross sections calculated by LINEAR-RECENT and NJOY have been analized. The two sets of pointwise cross sections differ significantly. Another disagreement was found in the transformation of ENDF/B-4 and 5 files into data with a linear interpolation scheme. Unshielded multigroup constants at O 0 K (620 groups, SANDII) have been calculated by the three code systems LINEAR-RECENT-GROUPIE, NJOY and RESEND5-INTEND. The code system RESEND5-INTEND calculates wrong group constants and should not be used any more. The two sets of group constants obtained from ENDF/B-4 data using GROUPIE and NJOY differ for some group constants by more than 2%. Some disagreements at low energies (10 -3 -eV) of the total cross section of Na-23 and Al-27 are difficult to understand. For ENDF/B-5 dosimetry data the capture group constants differ significantly. (Author) [pt

Preprocessing in Matlab Inconsistent Linear System for a Meaningful Least Squares Solution

Science.gov (United States)

Sen, Symal K.; Shaykhian, Gholam Ali

2011-01-01

Mathematical models of many physical/statistical problems are systems of linear equations Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the . minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

A real-time Java tool chain for resource constrained platforms

DEFF Research Database (Denmark)

Korsholm, Stephan E.; Søndergaard, Hans; Ravn, Anders Peter

2014-01-01

The Java programming language was originally developed for embedded systems, but the resource requirements of previous and current Java implementations – especially memory consumption – tend to exclude them from being used on a significant class of resource constrained embedded platforms. The con......The Java programming language was originally developed for embedded systems, but the resource requirements of previous and current Java implementations – especially memory consumption – tend to exclude them from being used on a significant class of resource constrained embedded platforms...... by integrating the following: (1) a lean virtual machine without any external dependencies on POSIX-like libraries or other OS functionalities; (2) a hardware abstraction layer, implemented almost entirely in Java through the use of hardware objects, first level interrupt handlers, and native variables; and (3....... An evaluation of the presented solution shows that the miniCDj benchmark gets reduced to a size where it can run on resource constrained platforms....

Development of Kinematic 3D Laser Scanning System for Indoor Mapping and As-Built BIM Using Constrained SLAM

Directory of Open Access Journals (Sweden)

Jaehoon Jung

2015-10-01

Full Text Available The growing interest and use of indoor mapping is driving a demand for improved data-acquisition facility, efficiency and productivity in the era of the Building Information Model (BIM. The conventional static laser scanning method suffers from some limitations on its operability in complex indoor environments, due to the presence of occlusions. Full scanning of indoor spaces without loss of information requires that surveyors change the scanner position many times, which incurs extra work for registration of each scanned point cloud. Alternatively, a kinematic 3D laser scanning system, proposed herein, uses line-feature-based Simultaneous Localization and Mapping (SLAM technique for continuous mapping. Moreover, to reduce the uncertainty of line-feature extraction, we incorporated constrained adjustment based on an assumption made with respect to typical indoor environments: that the main structures are formed of parallel or orthogonal line features. The superiority of the proposed constrained adjustment is its reduction for uncertainties of the adjusted lines, leading to successful data association process. In the present study, kinematic scanning with and without constrained adjustment were comparatively evaluated in two test sites, and the results confirmed the effectiveness of the proposed system. The accuracy of the 3D mapping result was additionally evaluated by comparison with the reference points acquired by a total station: the Euclidean average distance error was 0.034 m for the seminar room and 0.043 m for the corridor, which satisfied the error tolerance for point cloud acquisition (0.051 m according to the guidelines of the General Services Administration for BIM accuracy.

Decentralized linear quadratic power system stabilizers for multi ...

Indian Academy of Sciences (India)

Linear quadratic stabilizers are well-known for their superior control capabilities when compared to the conventional lead–lag power system stabilizers. However, they have not seen much of practical importance as the state variables are generally not measurable; especially the generator rotor angle measurement is not ...

Stability Analysis for Multi-Parameter Linear Periodic Systems

DEFF Research Database (Denmark)

Seyranian, A.P.; Solem, Frederik; Pedersen, Pauli

1999-01-01

This paper is devoted to stability analysis of general linear periodic systems depending on real parameters. The Floquet method and perturbation technique are the basis of the development. We start out with the first and higher-order derivatives of the Floquet matrix with respect to problem...

Quantization of the Linearized Kepler Problem

OpenAIRE

Guerrero, Julio; Perez, Jose Miguel

2003-01-01

The linearized Kepler problem is considered, as obtained from the Kustaanheimo-Stiefel (K-S)transformation, both for negative and positive energies. The symmetry group for the Kepler problem turns out to be SU(2,2). For negative energies, the Hamiltonian of Kepler problem can be realized as the sum of the energies of four harmonic oscillator with the same frequency, with a certain constrain. For positive energies, it can be realized as the sum of the energies of four repulsive oscillator with...

An algorithm for mass matrix calculation of internally constrained molecular geometries

International Nuclear Information System (INIS)

Aryanpour, Masoud; Dhanda, Abhishek; Pitsch, Heinz

2008-01-01

Dynamic models for molecular systems require the determination of corresponding mass matrix. For constrained geometries, these computations are often not trivial but need special considerations. Here, assembling the mass matrix of internally constrained molecular structures is formulated as an optimization problem. Analytical expressions are derived for the solution of the different possible cases depending on the rank of the constraint matrix. Geometrical interpretations are further used to enhance the solution concept. As an application, we evaluate the mass matrix for a constrained molecule undergoing an electron-transfer reaction. The preexponential factor for this reaction is computed based on the harmonic model

An algorithm for mass matrix calculation of internally constrained molecular geometries.

Science.gov (United States)

Aryanpour, Masoud; Dhanda, Abhishek; Pitsch, Heinz

2008-01-28

Dynamic models for molecular systems require the determination of corresponding mass matrix. For constrained geometries, these computations are often not trivial but need special considerations. Here, assembling the mass matrix of internally constrained molecular structures is formulated as an optimization problem. Analytical expressions are derived for the solution of the different possible cases depending on the rank of the constraint matrix. Geometrical interpretations are further used to enhance the solution concept. As an application, we evaluate the mass matrix for a constrained molecule undergoing an electron-transfer reaction. The preexponential factor for this reaction is computed based on the harmonic model.

Non-linear feedback control of the p53 protein-mdm2 inhibitor system using the derivative-free non-linear Kalman filter.

Science.gov (United States)

Rigatos, Gerasimos G

2016-06-01

It is proven that the model of the p53-mdm2 protein synthesis loop is a differentially flat one and using a diffeomorphism (change of state variables) that is proposed by differential flatness theory it is shown that the protein synthesis model can be transformed into the canonical (Brunovsky) form. This enables the design of a feedback control law that maintains the concentration of the p53 protein at the desirable levels. To estimate the non-measurable elements of the state vector describing the p53-mdm2 system dynamics, the derivative-free non-linear Kalman filter is used. Moreover, to compensate for modelling uncertainties and external disturbances that affect the p53-mdm2 system, the derivative-free non-linear Kalman filter is re-designed as a disturbance observer. The derivative-free non-linear Kalman filter consists of the Kalman filter recursion applied on the linearised equivalent of the protein synthesis model together with an inverse transformation based on differential flatness theory that enables to retrieve estimates for the state variables of the initial non-linear model. The proposed non-linear feedback control and perturbations compensation method for the p53-mdm2 system can result in more efficient chemotherapy schemes where the infusion of medication will be better administered.

Design techniques for large scale linear measurement systems

International Nuclear Information System (INIS)

Candy, J.V.

1979-03-01

Techniques to design measurement schemes for systems modeled by large scale linear time invariant systems, i.e., physical systems modeled by a large number (> 5) of ordinary differential equations, are described. The techniques are based on transforming the physical system model to a coordinate system facilitating the design and then transforming back to the original coordinates. An example of a three-stage, four-species, extraction column used in the reprocessing of spent nuclear fuel elements is presented. The basic ideas are briefly discussed in the case of noisy measurements. An example using a plutonium nitrate storage vessel (reprocessing) with measurement uncertainty is also presented

Maximum principle for a stochastic delayed system involving terminal state constraints.

Science.gov (United States)

Wen, Jiaqiang; Shi, Yufeng

2017-01-01

We investigate a stochastic optimal control problem where the controlled system is depicted as a stochastic differential delayed equation; however, at the terminal time, the state is constrained in a convex set. We firstly introduce an equivalent backward delayed system depicted as a time-delayed backward stochastic differential equation. Then a stochastic maximum principle is obtained by virtue of Ekeland's variational principle. Finally, applications to a state constrained stochastic delayed linear-quadratic control model and a production-consumption choice problem are studied to illustrate the main obtained result.

Quadratic Plus Linear Operators which Preserve Pure States of Quantum Systems: Small Dimensions

International Nuclear Information System (INIS)

Saburov, Mansoor

2014-01-01

A mathematical formalism of quantum mechanics says that a pure state of a quantum system corresponds to a vector of norm 1 and an observable is a self-adjoint operator on the space of states. It is of interest to describe all linear or nonlinear operators which preserve the pure states of the system. In the linear case, it is nothing more than isometries of Hilbert spaces. In the nonlinear case, this problem was open. In this paper, in the small dimensional spaces, we shall describe all quadratic plus linear operators which preserve pure states of the quantum system

Development of linear flow rate control system for eccentric butter-fly valve

International Nuclear Information System (INIS)

Kwak, K. K.; Cho, S. W.; Park, J. S.; Cho, J. H.; Song, I. T.; Kim, J. G.; Kwon, S. J.; Kim, I. J.; Park, W. K.

1999-12-01

Butter-fly valves are advantageous over gate, globe, plug, and ball valves in a variety of installations, particularly in the large sizes. The purpose of this project development of linear flow rate control system for eccentric butter-fly valve (intelligent butter-fly valve system). The intelligent butter-fly valve system consist of a valve body, micro controller. The micro controller consist of torque control system, pressure censor, worm and worm gear and communication line etc. The characteristics of intelligent butter-fly valve system as follows: Linear flow rate control function. Digital remote control function. guard function. Self-checking function. (author)

Generation Method of Multipiecewise Linear Chaotic Systems Based on the Heteroclinic Shil’nikov Theorem and Switching Control

Directory of Open Access Journals (Sweden)

Chunyan Han

2015-01-01

Full Text Available Based on the heteroclinic Shil’nikov theorem and switching control, a kind of multipiecewise linear chaotic system is constructed in this paper. Firstly, two fundamental linear systems are constructed via linearization of a chaotic system at its two equilibrium points. Secondly, a two-piecewise linear chaotic system which satisfies the Shil’nikov theorem is generated by constructing heteroclinic loop between equilibrium points of the two fundamental systems by switching control. Finally, another multipiecewise linear chaotic system that also satisfies the Shil’nikov theorem is obtained via alternate translation of the two fundamental linear systems and heteroclinic loop construction of adjacent equilibria for the multipiecewise linear system. Some basic dynamical characteristics, including divergence, Lyapunov exponents, and bifurcation diagrams of the constructed systems, are analyzed. Meanwhile, computer simulation and circuit design are used for the proposed chaotic systems, and they are demonstrated to be effective for the method of chaos anticontrol.

Synchronization of linearly coupled unified chaotic systems based on linear balanced feedback scheme with constraints

International Nuclear Information System (INIS)

Chen, H.-H.; Chen, C.-S.; Lee, C.-I

2009-01-01

This paper investigates the synchronization of unidirectional and bidirectional coupled unified chaotic systems. A balanced coupling coefficient control method is presented for global asymptotic synchronization using the Lyapunov stability theorem and a minimum scheme with no constraints/constraints. By using the result of the above analysis, the balanced coupling coefficients are then designed to achieve the chaos synchronization of linearly coupled unified chaotic systems. The feasibility and effectiveness of the proposed chaos synchronization scheme are verified via numerical simulations.

Fundamentals of linear systems for physical scientists and engineers

CERN Document Server

Puri, N N

2009-01-01

Thanks to the advent of inexpensive computing, it is possible to analyze, compute, and develop results that were unthinkable in the '60s. Control systems, telecommunications, robotics, speech, vision, and digital signal processing are but a few examples of computing applications. While there are many excellent resources available that focus on one or two topics, few books cover most of the mathematical techniques required for a broader range of applications. Fundamentals of Linear Systems for Physical Scientists and Engineers is such a resource. The book draws from diverse areas of engineering and the physical sciences to cover the fundamentals of linear systems. Assuming no prior knowledge of complex mathematics on the part of the reader, the author uses his nearly 50 years of teaching experience to address all of the necessary mathematical techniques. Original proofs, hundreds of examples, and proven theorems illustrate and clarify the material. An extensive table provides Lyapunov functions for differentia...

ORACLS- OPTIMAL REGULATOR ALGORITHMS FOR THE CONTROL OF LINEAR SYSTEMS (CDC VERSION)

Science.gov (United States)

Armstrong, E. S.

1994-01-01

This control theory design package, called Optimal Regulator Algorithms for the Control of Linear Systems (ORACLS), was developed to aid in the design of controllers and optimal filters for systems which can be modeled by linear, time-invariant differential and difference equations. Optimal linear quadratic regulator theory, currently referred to as the Linear-Quadratic-Gaussian (LQG) problem, has become the most widely accepted method of determining optimal control policy. Within this theory, the infinite duration time-invariant problems, which lead to constant gain feedback control laws and constant Kalman-Bucy filter gains for reconstruction of the system state, exhibit high tractability and potential ease of implementation. A variety of new and efficient methods in the field of numerical linear algebra have been combined into the ORACLS program, which provides for the solution to time-invariant continuous or discrete LQG problems. The ORACLS package is particularly attractive to the control system designer because it provides a rigorous tool for dealing with multi-input and multi-output dynamic systems in both continuous and discrete form. The ORACLS programming system is a collection of subroutines which can be used to formulate, manipulate, and solve various LQG design problems. The ORACLS program is constructed in a manner which permits the user to maintain considerable flexibility at each operational state. This flexibility is accomplished by providing primary operations, analysis of linear time-invariant systems, and control synthesis based on LQG methodology. The input-output routines handle the reading and writing of numerical matrices, printing heading information, and accumulating output information. The basic vector-matrix operations include addition, subtraction, multiplication, equation, norm construction, tracing, transposition, scaling, juxtaposition, and construction of null and identity matrices. The analysis routines provide for the following

Stability analysis of linear switching systems with time delays

International Nuclear Information System (INIS)

Li Ping; Zhong Shouming; Cui Jinzhong

2009-01-01

The issue of stability analysis of linear switching system with discrete and distributed time delays is studied in this paper. An appropriate switching rule is applied to guarantee the stability of the whole switching system. Our results use a Riccati-type Lyapunov functional under a condition on the time delay. So, switching systems with mixed delays are developed. A numerical example is given to illustrate the effectiveness of our results.

A Sparse Approximate Inverse Preconditioner for Nonsymmetric Linear Systems

Czech Academy of Sciences Publication Activity Database

Benzi, M.; Tůma, Miroslav

1998-01-01

Roč. 19, č. 3 (1998), s. 968-994 ISSN 1064-8275 R&D Projects: GA ČR GA201/93/0067; GA AV ČR IAA230401 Keywords : large sparse systems * interative methods * preconditioning * approximate inverse * sparse linear systems * sparse matrices * incomplete factorizations * conjugate gradient -type methods Subject RIV: BA - General Mathematics Impact factor: 1.378, year: 1998

Iterative solution of large linear systems

CERN Document Server

Young, David Matheson

1971-01-01

This self-contained treatment offers a systematic development of the theory of iterative methods. Its focal point resides in an analysis of the convergence properties of the successive overrelaxation (SOR) method, as applied to a linear system with a consistently ordered matrix. The text explores the convergence properties of the SOR method and related techniques in terms of the spectral radii of the associated matrices as well as in terms of certain matrix norms. Contents include a review of matrix theory and general properties of iterative methods; SOR method and stationary modified SOR meth

Adaptive H∞ synchronization of chaotic systems via linear and nonlinear feedback control

International Nuclear Information System (INIS)

Fu Shi-Hui; Lu Qi-Shao; Du Ying

2012-01-01

Adaptive H ∞ synchronization of chaotic systems via linear and nonlinear feedback control is investigated. The chaotic systems are redesigned by using the generalized Hamiltonian systems and observer approach. Based on Lyapunov's stability theory, linear and nonlinear feedback control of adaptive H ∞ synchronization is established in order to not only guarantee stable synchronization of both master and slave systems but also reduce the effect of external disturbance on an H ∞ -norm constraint. Adaptive H ∞ synchronization of chaotic systems via three kinds of control is investigated with applications to Lorenz and Chen systems. Numerical simulations are also given to identify the effectiveness of the theoretical analysis. (general)

Relaxation Methods for Strictly Convex Regularizations of Piecewise Linear Programs

International Nuclear Information System (INIS)

Kiwiel, K. C.

1998-01-01

We give an algorithm for minimizing the sum of a strictly convex function and a convex piecewise linear function. It extends several dual coordinate ascent methods for large-scale linearly constrained problems that occur in entropy maximization, quadratic programming, and network flows. In particular, it may solve exact penalty versions of such (possibly inconsistent) problems, and subproblems of bundle methods for nondifferentiable optimization. It is simple, can exploit sparsity, and in certain cases is highly parallelizable. Its global convergence is established in the recent framework of B -functions (generalized Bregman functions)

Time-dependent switched discrete-time linear systems control and filtering

CERN Document Server

Zhang, Lixian; Shi, Peng; Lu, Qiugang

2016-01-01

This book focuses on the basic control and filtering synthesis problems for discrete-time switched linear systems under time-dependent switching signals. Chapter 1, as an introduction of the book, gives the backgrounds and motivations of switched systems, the definitions of the typical time-dependent switching signals, the differences and links to other types of systems with hybrid characteristics and a literature review mainly on the control and filtering for the underlying systems. By summarizing the multiple Lyapunov-like functions (MLFs) approach in which different requirements on comparisons of Lyapunov function values at switching instants, a series of methodologies are developed for the issues on stability and stabilization, and l2-gain performance or tube-based robustness for l∞ disturbance, respectively, in Chapters 2 and 3. Chapters 4 and 5 are devoted to the control and filtering problems for the time-dependent switched linear systems with either polytopic uncertainties or measurable time-varying...

Comparison of modal spectral and non-linear time history analysis of a piping system

International Nuclear Information System (INIS)

Gerard, R.; Aelbrecht, D.; Lafaille, J.P.

1987-01-01

A typical piping system of the discharge line of the chemical and volumetric control system, outside the containment, between the penetration and the heat exchanger, an operating power plant was analyzed using four different methods: Modal spectral analysis with 2% constant damping, modal spectral analysis using ASME Code Case N411 (PVRC damping), linear time history analysis, non-linear time history analysis. This paper presents an estimation of the conservatism of the linear methods compared to the non-linear analysis. (orig./HP)

A constrained robust least squares approach for contaminant release history identification

Science.gov (United States)

Sun, Alexander Y.; Painter, Scott L.; Wittmeyer, Gordon W.

2006-04-01

Contaminant source identification is an important type of inverse problem in groundwater modeling and is subject to both data and model uncertainty. Model uncertainty was rarely considered in the previous studies. In this work, a robust framework for solving contaminant source recovery problems is introduced. The contaminant source identification problem is first cast into one of solving uncertain linear equations, where the response matrix is constructed using a superposition technique. The formulation presented here is general and is applicable to any porous media flow and transport solvers. The robust least squares (RLS) estimator, which originated in the field of robust identification, directly accounts for errors arising from model uncertainty and has been shown to significantly reduce the sensitivity of the optimal solution to perturbations in model and data. In this work, a new variant of RLS, the constrained robust least squares (CRLS), is formulated for solving uncertain linear equations. CRLS allows for additional constraints, such as nonnegativity, to be imposed. The performance of CRLS is demonstrated through one- and two-dimensional test problems. When the system is ill-conditioned and uncertain, it is found that CRLS gave much better performance than its classical counterpart, the nonnegative least squares. The source identification framework developed in this work thus constitutes a reliable tool for recovering source release histories in real applications.

Robustness-tracking control based on sliding mode and H∞ theory for linear servo system

Institute of Scientific and Technical Information of China (English)

TIAN Yan-feng; GUO Qing-ding

2005-01-01

A robustness-tracking control scheme based on combining H∞ robust control and sliding mode control is proposed for a direct drive AC permanent-magnet linear motor servo system to solve the conflict between tracking and robustness of the linear servo system. The sliding mode tracking controller is designed to ensure the system has a fast tracking characteristic to the command, and the H∞ robustness controller suppresses the disturbances well within the close loop( including the load and the end effect force of linear motor etc. ) and effectively minimizes the chattering of sliding mode control which influences the steady state performance of the system. Simulation results show that this control scheme enhances the track-command-ability and the robustness of the linear servo system, and in addition, it has a strong robustness to parameter variations and resistance disturbances.

On energy efficient power allocation for power-constrained systems

KAUST Repository

Sboui, Lokman

2014-09-01

Recently, the energy efficiency (EE) has become an important factor when designing new wireless communication systems. Due to economic and environmental challenges, new trends and efforts are oriented toward “green” communication especially for energy-constrained applications such as wireless sensors network and cognitive radio. To this end, we analyze the power allocation scheme that maximizes the EE defined as rate over the total power including circuit power. We derive an explicit expression of the optimal power with instantaneous channel gain based on EE criterion. We show that the relation between the EE and the spectral efficiency (SE) when the optimal power is adopted is strictly increasing in contrast with the SE-EE trade-off discussed in the literature. We also solve a non-convex problem and compute explicitly the optimal power for ergodic EE under either a peak or an average power constraint. When the instantaneous channel is not available, we provide the optimal power equation and compute simple sub-optimal power. In the numerical results, we show that the sup-optimal solution is very close to the optimal solution. In addition, we show that the absence of the channel state information (CSI) only affects the EE and the SE performances at high power regime compared to the full CSI case.

Constrained Quantum Mechanics: Chaos in Non-Planar Billiards

Science.gov (United States)

Salazar, R.; Tellez, G.

2012-01-01

We illustrate some of the techniques to identify chaos signatures at the quantum level using as guiding examples some systems where a particle is constrained to move on a radial symmetric, but non-planar, surface. In particular, two systems are studied: the case of a cone with an arbitrary contour or "dunce hat billiard" and the rectangular…

Constrained non-linear multi-objective optimisation of preventive maintenance scheduling for offshore wind farms

Science.gov (United States)

Zhong, Shuya; Pantelous, Athanasios A.; Beer, Michael; Zhou, Jian

2018-05-01

Offshore wind farm is an emerging source of renewable energy, which has been shown to have tremendous potential in recent years. In this blooming area, a key challenge is that the preventive maintenance of offshore turbines should be scheduled reasonably to satisfy the power supply without failure. In this direction, two significant goals should be considered simultaneously as a trade-off. One is to maximise the system reliability and the other is to minimise the maintenance related cost. Thus, a non-linear multi-objective programming model is proposed including two newly defined objectives with thirteen families of constraints suitable for the preventive maintenance of offshore wind farms. In order to solve our model effectively, the nondominated sorting genetic algorithm II, especially for the multi-objective optimisation is utilised and Pareto-optimal solutions of schedules can be obtained to offer adequate support to decision-makers. Finally, an example is given to illustrate the performances of the devised model and algorithm, and explore the relationships of the two targets with the help of a contrast model.

Macroscopically constrained Wang-Landau method for systems with multiple order parameters and its application to drawing complex phase diagrams

Science.gov (United States)

Chan, C. H.; Brown, G.; Rikvold, P. A.

2017-05-01

A generalized approach to Wang-Landau simulations, macroscopically constrained Wang-Landau, is proposed to simulate the density of states of a system with multiple macroscopic order parameters. The method breaks a multidimensional random-walk process in phase space into many separate, one-dimensional random-walk processes in well-defined subspaces. Each of these random walks is constrained to a different set of values of the macroscopic order parameters. When the multivariable density of states is obtained for one set of values of fieldlike model parameters, the density of states for any other values of these parameters can be obtained by a simple transformation of the total system energy. All thermodynamic quantities of the system can then be rapidly calculated at any point in the phase diagram. We demonstrate how to use the multivariable density of states to draw the phase diagram, as well as order-parameter probability distributions at specific phase points, for a model spin-crossover material: an antiferromagnetic Ising model with ferromagnetic long-range interactions. The fieldlike parameters in this model are an effective magnetic field and the strength of the long-range interaction.

Synchronization and Control of Linearly Coupled Singular Systems

Directory of Open Access Journals (Sweden)

Fang Qingxiang

2013-01-01

Full Text Available The synchronization and control problem of linearly coupled singular systems is investigated. The uncoupled dynamical behavior at each node is general and can be chaotic or, otherwise the coupling matrix is not assumed to be symmetrical. Some sufficient conditions for globally exponential synchronization are derived based on Lyapunov stability theory. These criteria, which are in terms of linear matrix inequality (LMI, indicate that the left and right eigenvectors corresponding to eigenvalue zero of the coupling matrix play key roles in the stability analysis of the synchronization manifold. The controllers are designed for state feedback control and pinning control, respectively. Finally, a numerical example is provided to illustrate the effectiveness of the proposed conditions.

Unification of three linear models for the transient visual system

NARCIS (Netherlands)

Brinker, den A.C.

1989-01-01

Three different linear filters are considered as a model describing the experimentally determined triphasic impulse responses of discs. These impulse responses arc associated with the transient visual system. Each model reveals a different feature of the system. Unification of the models is

Solution of linear and nonlinear matrix systems. Application to a nonlinear diffusion equation

International Nuclear Information System (INIS)

Bonnet, M.; Meurant, G.

1978-01-01

Different methods of solution of linear and nonlinear algebraic systems are applied to the nonlinear system obtained by discretizing a nonlinear diffusion equation. For linear systems, methods in general use of alternating directions type or Gauss Seidel's methods are compared to more recent ones of the type of generalized conjugate gradient; the superiority of the latter is shown by numerical examples. For nonlinear systems, a method on nonlinear conjugate gradient is studied as also Newton's method and some of its variants. It should be noted, however that Newton's method is found to be more efficient when coupled with a good method for solution of the linear system. To conclude, such methods are used to solve a nonlinear diffusion problem and the numerical results obtained are to be compared [fr

Solution of linear and nonlinear matrix systems. Application to a nonlinear diffusion equation

International Nuclear Information System (INIS)

Bonnet, M.; Meurant, G.

1978-01-01

The object of this study is to compare different methods of solving linear and nonlinear algebraic systems and to apply them to the nonlinear system obtained by discretizing a nonlinear diffusion equation. For linear systems the conventional methods of alternating direction type or Gauss Seidel's methods are compared to more recent ones of the type of generalized conjugate gradient; the superiority of the latter is shown by numerical examples. For nonlinear systems, a method of nonlinear conjugate gradient is studied together with Newton's method and some of its variants. It should be noted, however, that Newton's method is found to be more efficient when coupled with a good method for solving the linear system. As a conclusion, these methods are used to solve a nonlinear diffusion problem and the numerical results obtained are compared [fr

Linear local stability of electrostatic drift modes in helical systems

International Nuclear Information System (INIS)

Yamagishi, O.; Nakajima, N.; Sugama, H.; Nakamura, Y.

2003-01-01

We investigate the stability of the drift wave in helical systems. For this purpose, we solve the linear local gyrokinetic-Poisson equation, in the electrostatic regime. As a model of helical plasmas, Large helical Device (LHD) is considered. The equation we apply is rather exact in the framework of linear gyrokinetic theory, where only the approximation is the ballooning representation. In this paper, we consider only collisionless cases. All the frequency regime can be naturally reated without any assumptions, and in such cases, ion temperature gradient modes (ITG), trapped electron modes (TEM), and electron temperature gradient modes (ETG) are expected to become unstable linearly independently. (orig.)

Recovery coefficients as a test of system linearity of response in PET

International Nuclear Information System (INIS)

Geworski, L.; Munz, D.L.; Knoop, B.; Hofmann, M.; Knapp, W.H.

2002-01-01

Aim: New imaging protocols have created an increasing demand for quantitation in dedicated PET. Besides attenuation and scatter correction the recovery correction, accounting for the instrument's limited spatial resolution, has gained importance. For clinical practicability these corrections should work independent from the object, i.e. from the actual distribution of emitter and absorber. Aim of the study was to test this object independency, i.e. system linearity of response, by comparing recovery coefficients (RC) determined for different object geometries. In fact, this comparison may serve as a final test on system linearity of response, as measured on the quantitative accuracy by which the activity concentration in small lesions can be recovered. Method: For hot and cold spot imaging situations spatial distribution of activity is different. Therefore, scatter correction algorithm has to deal with different scatter distributions. If all factors disturbing system linearity, specifically scatter and attenuation, are corrected to a sufficient degree of accuracy, the system behaves linearly resulting in the theoretical relationship. CSRC = (1-HSRC). Thus, this equation, applied hot and cold spot measurements, will serve as a test on the effectiveness of the corrections and, hence, as a test of system linearity of response. Following IEC standard procedures (IEC 61675-1) measurements were done with and without interplane septa (2D/3D) on an ECAT EXACT 922 using a cylindrical phantom containing six spheres of different diameters (10 mm - 40 mm). All data were corrected for attenuation (transmission scan) and scatter (2D: deconvolution, 3D: scatter model), as implemented in the scanner's standard software. Recovery coefficients were determined for cold (CSRC) and hot (HSRC) lesions using both 2D and 3D acquisition mode. Results: CSRC directly measured versus CSRC calculated according to eq. (1) from HSRC resulted in an excellent agreement for both 2D and 3D data

Shape space exploration of constrained meshes

KAUST Repository

Yang, Yongliang; Yang, Yijun; Pottmann, Helmut; Mitra, Niloy J.

2011-01-01

We present a general computational framework to locally characterize any shape space of meshes implicitly prescribed by a collection of non-linear constraints. We computationally access such manifolds, typically of high dimension and co-dimension, through first and second order approximants, namely tangent spaces and quadratically parameterized osculant surfaces. Exploration and navigation of desirable subspaces of the shape space with regard to application specific quality measures are enabled using approximants that are intrinsic to the underlying manifold and directly computable in the parameter space of the osculant surface. We demonstrate our framework on shape spaces of planar quad (PQ) meshes, where each mesh face is constrained to be (nearly) planar, and circular meshes, where each face has a circumcircle. We evaluate our framework for navigation and design exploration on a variety of inputs, while keeping context specific properties such as fairness, proximity to a reference surface, etc. © 2011 ACM.

Shape space exploration of constrained meshes

KAUST Repository

Yang, Yongliang

2011-12-12

We present a general computational framework to locally characterize any shape space of meshes implicitly prescribed by a collection of non-linear constraints. We computationally access such manifolds, typically of high dimension and co-dimension, through first and second order approximants, namely tangent spaces and quadratically parameterized osculant surfaces. Exploration and navigation of desirable subspaces of the shape space with regard to application specific quality measures are enabled using approximants that are intrinsic to the underlying manifold and directly computable in the parameter space of the osculant surface. We demonstrate our framework on shape spaces of planar quad (PQ) meshes, where each mesh face is constrained to be (nearly) planar, and circular meshes, where each face has a circumcircle. We evaluate our framework for navigation and design exploration on a variety of inputs, while keeping context specific properties such as fairness, proximity to a reference surface, etc. © 2011 ACM.

MIROS: a hybrid real-time energy-efficient operating system for the resource-constrained wireless sensor nodes.

Science.gov (United States)

Liu, Xing; Hou, Kun Mean; de Vaulx, Christophe; Shi, Hongling; El Gholami, Khalid

2014-09-22

Operating system (OS) technology is significant for the proliferation of the wireless sensor network (WSN). With an outstanding OS; the constrained WSN resources (processor; memory and energy) can be utilized efficiently. Moreover; the user application development can be served soundly. In this article; a new hybrid; real-time; memory-efficient; energy-efficient; user-friendly and fault-tolerant WSN OS MIROS is designed and implemented. MIROS implements the hybrid scheduler and the dynamic memory allocator. Real-time scheduling can thus be achieved with low memory consumption. In addition; it implements a mid-layer software EMIDE (Efficient Mid-layer Software for User-Friendly Application Development Environment) to decouple the WSN application from the low-level system. The application programming process can consequently be simplified and the application reprogramming performance improved. Moreover; it combines both the software and the multi-core hardware techniques to conserve the energy resources; improve the node reliability; as well as achieve a new debugging method. To evaluate the performance of MIROS; it is compared with the other WSN OSes (TinyOS; Contiki; SOS; openWSN and mantisOS) from different OS concerns. The final evaluation results prove that MIROS is suitable to be used even on the tight resource-constrained WSN nodes. It can support the real-time WSN applications. Furthermore; it is energy efficient; user friendly and fault tolerant.

MIROS: A Hybrid Real-Time Energy-Efficient Operating System for the Resource-Constrained Wireless Sensor Nodes

Science.gov (United States)

Liu, Xing; Hou, Kun Mean; de Vaulx, Christophe; Shi, Hongling; Gholami, Khalid El

2014-01-01

Operating system (OS) technology is significant for the proliferation of the wireless sensor network (WSN). With an outstanding OS; the constrained WSN resources (processor; memory and energy) can be utilized efficiently. Moreover; the user application development can be served soundly. In this article; a new hybrid; real-time; memory-efficient; energy-efficient; user-friendly and fault-tolerant WSN OS MIROS is designed and implemented. MIROS implements the hybrid scheduler and the dynamic memory allocator. Real-time scheduling can thus be achieved with low memory consumption. In addition; it implements a mid-layer software EMIDE (Efficient Mid-layer Software for User-Friendly Application Development Environment) to decouple the WSN application from the low-level system. The application programming process can consequently be simplified and the application reprogramming performance improved. Moreover; it combines both the software and the multi-core hardware techniques to conserve the energy resources; improve the node reliability; as well as achieve a new debugging method. To evaluate the performance of MIROS; it is compared with the other WSN OSes (TinyOS; Contiki; SOS; openWSN and mantisOS) from different OS concerns. The final evaluation results prove that MIROS is suitable to be used even on the tight resource-constrained WSN nodes. It can support the real-time WSN applications. Furthermore; it is energy efficient; user friendly and fault tolerant. PMID:25248069

Well logging system with linearity control

International Nuclear Information System (INIS)

Jones, J.M.

1973-01-01

Apparatus is described for controlling the gain of a nuclear well logging system comprising: (1) means for measuring the energy spectrum of gamma rays produced by earth formation materials surrounding a well borehole; (2) means for measuring the number of counts of a gamma rays having an energy falling within each of at least two predetermined energy band portions of the gamma ray energy spectrum; (3) means for generating a signal proportional to the ratio of the gamma ray counts and for comparing the ratio signal with at least one constant ratio calibration signal; (4) means for generating an error signal representative of the difference of the ratio signal and the constant ratio calibration signal; and (5) means for using the error signal to control the linearity of the well logging system. (author)

Refined Fuchs inequalities for systems of linear differential equations

International Nuclear Information System (INIS)

Gontsov, R R

2004-01-01

We refine the Fuchs inequalities obtained by Corel for systems of linear meromorphic differential equations given on the Riemann sphere. Fuchs inequalities enable one to estimate the sum of exponents of the system over all its singular points. We refine these well-known inequalities by considering the Jordan structure of the leading coefficient of the Laurent series for the matrix of the right-hand side of the system in the neighbourhood of a singular point

GPS-based ionospheric tomography with a constrained adaptive ...

Indian Academy of Sciences (India)

Gauss weighted function is introduced to constrain the tomography system in the new method. It can resolve the ... the research focus in the fields of space geodesy and ... ment of GNSS such as GPS, Glonass, Galileo and. Compass, as these ...

Principle and analysis of a linear motor driving system for HTS levitation applications

International Nuclear Information System (INIS)

Jin, Jian X.; Guo, You G.; Zhu, Jian G.

2007-01-01

High temperature superconductor (HTS) high levitation force density with passive and self-stabilizing features allows a number of special applications to be developed. Linear motor driving systems are commonly required for those applications such as levitated transport systems. In this paper a prototype linear motor driving system with HTS is analyzed with calculation details including its magnetic fields and driving forces presented in the paper

Integer Linear Programming for Constrained Multi-Aspect Committee Review Assignment

Science.gov (United States)

Karimzadehgan, Maryam; Zhai, ChengXiang

2011-01-01

Automatic review assignment can significantly improve the productivity of many people such as conference organizers, journal editors and grant administrators. A general setup of the review assignment problem involves assigning a set of reviewers on a committee to a set of documents to be reviewed under the constraint of review quota so that the reviewers assigned to a document can collectively cover multiple topic aspects of the document. No previous work has addressed such a setup of committee review assignments while also considering matching multiple aspects of topics and expertise. In this paper, we tackle the problem of committee review assignment with multi-aspect expertise matching by casting it as an integer linear programming problem. The proposed algorithm can naturally accommodate any probabilistic or deterministic method for modeling multiple aspects to automate committee review assignments. Evaluation using a multi-aspect review assignment test set constructed using ACM SIGIR publications shows that the proposed algorithm is effective and efficient for committee review assignments based on multi-aspect expertise matching. PMID:22711970

A Stationary North-Finding Scheme for an Azimuth Rotational IMU Utilizing a Linear State Equality Constraint

Directory of Open Access Journals (Sweden)

Huapeng Yu

2015-02-01

Full Text Available The Kalman filter (KF has always been used to improve north-finding performance under practical conditions. By analyzing the characteristics of the azimuth rotational inertial measurement unit (ARIMU on a stationary base, a linear state equality constraint for the conventional KF used in the fine north-finding filtering phase is derived. Then, a constrained KF using the state equality constraint is proposed and studied in depth. Estimation behaviors of the concerned navigation errors when implementing the conventional KF scheme and the constrained KF scheme during stationary north-finding are investigated analytically by the stochastic observability approach, which can provide explicit formulations of the navigation errors with influencing variables. Finally, multiple practical experimental tests at a fixed position are done on a postulate system to compare the stationary north-finding performance of the two filtering schemes. In conclusion, this study has successfully extended the utilization of the stochastic observability approach for analytic descriptions of estimation behaviors of the concerned navigation errors, and the constrained KF scheme has demonstrated its superiority over the conventional KF scheme for ARIMU stationary north-finding both theoretically and practically.