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
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
Iterative solution of large linear systems
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
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
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.)
Focal points and principal solutions of linear Hamiltonian systems revisited
Šepitka, Peter; Šimon Hilscher, Roman
2018-05-01
In this paper we present a novel view on the principal (and antiprincipal) solutions of linear Hamiltonian systems, as well as on the focal points of their conjoined bases. We present a new and unified theory of principal (and antiprincipal) solutions at a finite point and at infinity, and apply it to obtain new representation of the multiplicities of right and left proper focal points of conjoined bases. We show that these multiplicities can be characterized by the abnormality of the system in a neighborhood of the given point and by the rank of the associated T-matrix from the theory of principal (and antiprincipal) solutions. We also derive some additional important results concerning the representation of T-matrices and associated normalized conjoined bases. The results in this paper are new even for completely controllable linear Hamiltonian systems. We also discuss other potential applications of our main results, in particular in the singular Sturmian theory.
Solution methods for large systems of linear equations in BACCHUS
International Nuclear Information System (INIS)
Homann, C.; Dorr, B.
1993-05-01
The computer programme BACCHUS is used to describe steady state and transient thermal-hydraulic behaviour of a coolant in a fuel element with intact geometry in a fast breeder reactor. In such computer programmes generally large systems of linear equations with sparse matrices of coefficients, resulting from discretization of coolant conservation equations, must be solved thousands of times giving rise to large demands of main storage and CPU time. Direct and iterative solution methods of the systems of linear equations, available in BACCHUS, are described, giving theoretical details and experience with their use in the programme. Besides use of a method of lines, a Runge-Kutta-method, for solution of the partial differential equation is outlined. (orig.) [de
Solution of the fully fuzzy linear systems using iterative techniques
International Nuclear Information System (INIS)
Dehghan, Mehdi; Hashemi, Behnam; Ghatee, Mehdi
2007-01-01
This paper mainly intends to discuss the iterative solution of fully fuzzy linear systems which we call FFLS. We employ Dubois and Prade's approximate arithmetic operators on LR fuzzy numbers for finding a positive fuzzy vector x-tilde which satisfies A-tildex-tilde=b, where A-tilde and b-tilde are a fuzzy matrix and a fuzzy vector, respectively. Please note that the positivity assumption is not so restrictive in applied problems. We transform FFLS and propose iterative techniques such as Richardson, Jacobi, Jacobi overrelaxation (JOR), Gauss-Seidel, successive overrelaxation (SOR), accelerated overrelaxation (AOR), symmetric and unsymmetric SOR (SSOR and USSOR) and extrapolated modified Aitken (EMA) for solving FFLS. In addition, the methods of Newton, quasi-Newton and conjugate gradient are proposed from nonlinear programming for solving a fully fuzzy linear system. Various numerical examples are also given to show the efficiency of the proposed schemes
Riccati transformations and principal solutions of discrete linear systems
International Nuclear Information System (INIS)
Ahlbrandt, C.D.; Hooker, J.W.
1984-01-01
Consider a second-order linear matrix difference equation. A definition of principal and anti-principal, or recessive and dominant, solutions of the equation are given and the existence of principal and anti-principal solutions and the essential uniqueness of principal solutions is proven
Linear homotopy solution of nonlinear systems of equations in geodesy
Paláncz, Béla; Awange, Joseph L.; Zaletnyik, Piroska; Lewis, Robert H.
2010-01-01
A fundamental task in geodesy is solving systems of equations. Many geodetic problems are represented as systems of multivariate polynomials. A common problem in solving such systems is improper initial starting values for iterative methods, leading to convergence to solutions with no physical meaning, or to convergence that requires global methods. Though symbolic methods such as Groebner bases or resultants have been shown to be very efficient, i.e., providing solutions for determined systems such as 3-point problem of 3D affine transformation, the symbolic algebra can be very time consuming, even with special Computer Algebra Systems (CAS). This study proposes the Linear Homotopy method that can be implemented easily in high-level computer languages like C++ and Fortran that are faster than CAS by at least two orders of magnitude. Using Mathematica, the power of Homotopy is demonstrated in solving three nonlinear geodetic problems: resection, GPS positioning, and affine transformation. The method enlarging the domain of convergence is found to be efficient, less sensitive to rounding of numbers, and has lower complexity compared to other local methods like Newton-Raphson.
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 note on the time decay of solutions for the linearized Wigner-Poisson system
Gamba, Irene; Gualdani, Maria; Sparber, Christof
2009-01-01
We consider the one-dimensional Wigner-Poisson system of plasma physics, linearized around a (spatially homogeneous) Lorentzian distribution and prove that the solution of the corresponding linearized problem decays to zero in time. We also give
Solution of a System of Linear Equations with Fuzzy Numbers
Czech Academy of Sciences Publication Activity Database
Horčík, Rostislav
2008-01-01
Roč. 159, č. 14 (2008), s. 1788-1810 ISSN 0165-0114 R&D Projects: GA AV ČR KJB100300502 Institutional research plan: CEZ:AV0Z10300504 Keywords : fuzzy number * fuzzy interval * interval analysis * fuzzy arithmetic * fuzzy class theory * united solution set Subject RIV: BA - General Mathematics Impact factor: 1.833, year: 2008
DEFF Research Database (Denmark)
Ahmadi, Mohamadreza; Mojallali, Hamed; Wisniewski, Rafal
2012-01-01
This paper addresses the robust stability and control problem of uncertain piecewise linear switched systems where, instead of the conventional Carathe ́odory solutions, we allow for Filippov solutions. In other words, in contrast to the previous studies, solutions with infinite switching in fini...... algorithm is proposed to surmount the aforementioned matrix inequality conditions....... time along the facets and on faces of arbitrary dimensions are also taken into account. Firstly, based on earlier results, the stability problem of piecewise linear systems with Filippov solutions is translated into a number of linear matrix inequality feasibility tests. Subsequently, a set of matrix...... inequalities are brought forward, which determines the asymptotic stability of the Filippov solutions of a given uncertain piecewise linear system. Afterwards, bilinear matrix inequality conditions for synthesizing a robust controller with a guaranteed H∞ per- formance are formulated. Finally, a V-K iteration...
Solution of systems of linear algebraic equations by the method of summation of divergent series
International Nuclear Information System (INIS)
Kirichenko, G.A.; Korovin, Ya.S.; Khisamutdinov, M.V.; Shmojlov, V.I.
2015-01-01
A method for solving systems of linear algebraic equations has been proposed on the basis on the summation of the corresponding continued fractions. The proposed algorithm for solving systems of linear algebraic equations is classified as direct algorithms providing an exact solution in a finite number of operations. Examples of solving systems of linear algebraic equations have been presented and the effectiveness of the algorithm has been estimated [ru
A note on the time decay of solutions for the linearized Wigner-Poisson system
Gamba, Irene
2009-01-01
We consider the one-dimensional Wigner-Poisson system of plasma physics, linearized around a (spatially homogeneous) Lorentzian distribution and prove that the solution of the corresponding linearized problem decays to zero in time. We also give an explicit algebraic decay rate.
Asymptotic behavior of solutions of linear multi-order fractional differential equation systems
Diethelm, Kai; Siegmund, Stefan; Tuan, H. T.
2017-01-01
In this paper, we investigate some aspects of the qualitative theory for multi-order fractional differential equation systems. First, we obtain a fundamental result on the existence and uniqueness for multi-order fractional differential equation systems. Next, a representation of solutions of homogeneous linear multi-order fractional differential equation systems in series form is provided. Finally, we give characteristics regarding the asymptotic behavior of solutions to some classes of line...
Computer programs for the solution of systems of linear algebraic equations
Sequi, W. T.
1973-01-01
FORTRAN subprograms for the solution of systems of linear algebraic equations are described, listed, and evaluated in this report. Procedures considered are direct solution, iteration, and matrix inversion. Both incore methods and those which utilize auxiliary data storage devices are considered. Some of the subroutines evaluated require the entire coefficient matrix to be in core, whereas others account for banding or sparceness of the system. General recommendations relative to equation solving are made, and on the basis of tests, specific subprograms are recommended.
Preprocessing in Matlab Inconsistent Linear System for a Meaningful Least Squares Solution
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.
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
International Nuclear Information System (INIS)
Li Tatsien
1994-01-01
By means of the concept of the weak linear degeneracy, one gets the global existence and the sharp estimate of the lifespan of C 1 solutions to the Cauchy problem for general first order quasilinear hyperbolic systems with small initial data with compact support. (author). 23 refs, 1 fig
On the solution of a class of fuzzy system of linear equations
Indian Academy of Sciences (India)
J. Mathematics and Comput. Sci. 1: 1–5. Salkuyeh D K 2011 On the solution of the fuzzy Sylvester matrix equation. Soft Computing 15: 953–961. Senthilkumar P and Rajendran G 2011 New approach to solve symmetric fully fuzzy linear systems. S¯adhan¯a 36: 933–940. Wang K and Zheng B 2007 Block iterative methods ...
Parallelized preconditioned BiCGStab solution of sparse linear system equations in F-COBRA-TF
International Nuclear Information System (INIS)
Geemert, Rene van; Glück, Markus; Riedmann, Michael; Gabriel, Harry
2011-01-01
Recently, the in-house development of a preconditioned and parallelized BiCGStab solver has been pursued successfully in AREVA’s advanced sub-channel code F-COBRA-TF. This solver can be run either in a sequential computation mode on a single CPU, or in a parallel computation mode on multiple parallel CPUs. The developed procedure enables the computation of several thousands of successive sparse linear system solutions in F-COBRA-TF with acceptable wall clock run times. The current paper provides general information about F-COBRA-TF in terms of modeling capabilities and application areas, and points out where the relevance arises for the efficient iterative solution of sparse linear systems. Furthermore, the preconditioning and parallelization strategies in the developed BiCGStab iterative solution approach are discussed. The paper is concluded with a number of verification examples. (author)
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.
TOEPLITZ, Solution of Linear Equation System with Toeplitz or Circulant Matrix
International Nuclear Information System (INIS)
Garbow, B.
1984-01-01
Description of program or function: TOEPLITZ is a collection of FORTRAN subroutines for solving linear systems Ax=b, where A is a Toeplitz matrix, a Circulant matrix, or has one or several block structures based on Toeplitz or Circulant matrices. Such systems arise in problems of electrodynamics, acoustics, mathematical statistics, algebra, in the numerical solution of integral equations with a difference kernel, and in the theory of stationary time series and signals
An efficient parallel algorithm for the solution of a tridiagonal linear system of equations
Stone, H. S.
1971-01-01
Tridiagonal linear systems of equations are solved on conventional serial machines in a time proportional to N, where N is the number of equations. The conventional algorithms do not lend themselves directly to parallel computations on computers of the ILLIAC IV class, in the sense that they appear to be inherently serial. An efficient parallel algorithm is presented in which computation time grows as log sub 2 N. The algorithm is based on recursive doubling solutions of linear recurrence relations, and can be used to solve recurrence relations of all orders.
The numerical solution of linear multi-term fractional differential equations: systems of equations
Edwards, John T.; Ford, Neville J.; Simpson, A. Charles
2002-11-01
In this paper, we show how the numerical approximation of the solution of a linear multi-term fractional differential equation can be calculated by reduction of the problem to a system of ordinary and fractional differential equations each of order at most unity. We begin by showing how our method applies to a simple class of problems and we give a convergence result. We solve the Bagley Torvik equation as an example. We show how the method can be applied to a general linear multi-term equation and give two further examples.
On the economical solution method for a system of linear algebraic equations
Directory of Open Access Journals (Sweden)
Jan Awrejcewicz
2004-01-01
Full Text Available The present work proposes a novel optimal and exact method of solving large systems of linear algebraic equations. In the approach under consideration, the solution of a system of algebraic linear equations is found as a point of intersection of hyperplanes, which needs a minimal amount of computer operating storage. Two examples are given. In the first example, the boundary value problem for a three-dimensional stationary heat transfer equation in a parallelepiped in ℝ3 is considered, where boundary value problems of first, second, or third order, or their combinations, are taken into account. The governing differential equations are reduced to algebraic ones with the help of the finite element and boundary element methods for different meshes applied. The obtained results are compared with known analytical solutions. The second example concerns computation of a nonhomogeneous shallow physically and geometrically nonlinear shell subject to transversal uniformly distributed load. The partial differential equations are reduced to a system of nonlinear algebraic equations with the error of O(hx12+hx22. The linearization process is realized through either Newton method or differentiation with respect to a parameter. In consequence, the relations of the boundary condition variations along the shell side and the conditions for the solution matching are reported.
Non linear Euler-Poisson system. Part 1: global existence of low entropy solutions
International Nuclear Information System (INIS)
Cordier, S.
1995-05-01
In this work a 1-D model of electrons and ions plasma is considered. Electrons are supposed to be in Maxwell-Boltzmann thermodynamic equilibrium while ions are described with an isothermal flow model of charged particles submitted to a self-consistent electric field. A collision term between neutral particles and ions simulates the presence of neutral particles. This work demonstrates the existence of low entropy solutions for this simple model with arbitrary initial conditions. Most of the paper is devoted to the demonstration of this theorem and follows the successive steps: construction of a numerical scheme, recall of the classical properties of Riemann problem solutions using Glimm method, uniform estimations for the whole variation norm, and finally, convergence of the constructed solutions towards a low entropy solution for the non-linear Euler/Poisson system. Domains of application for this type of model are listed in the conclusion. (J.S.). 18 refs
Economic planning for electric energy systems: a multi objective linearized approach for solution
International Nuclear Information System (INIS)
Mata Medeiros Branco, T. da.
1986-01-01
The economic planning problem associated to the expansion and operation of electrical power systems is considered in this study, represented for a vectorial objective function in which the minimization of resources involved and maximization of attended demand constitute goals to be satisfied. Supposing all the variables involved with linear characteristic and considering the conflict existing among the objectives to be achieved, in order to find a solution, a multi objective linearized approach is proposed. This approximation utilizes the compromise programming technique and linear programming methods. Generation and transmission are simultaneously considered into the optimization process in which associated losses and the capacity of each line are included. Illustrated examples are also presented with results discussed. (author)
Asymptotical Behavior of the Solution of a SDOF Linear Fractionally Damped Vibration System
Directory of Open Access Journals (Sweden)
Z.H. Wang
2011-01-01
Full Text Available Fractional-order derivative has been shown an adequate tool to the study of so-called "anomalous" social and physical behaviors, in reflecting their non-local, frequency- and history-dependent properties, and it has been used to model practical systems in engineering successfully, including the famous Bagley-Torvik equation modeling forced motion of a rigid plate immersed in Newtonian fluid. The solutions of the initial value problems of linear fractional differential equations are usually expressed in terms of Mittag-Leffler functions or some other kind of power series. Such forms of solutions are not good for engineers not only in understanding the solutions but also in investigation. This paper proves that for the linear SDOF oscillator with a damping described by fractional-order derivative whose order is between 1 and 2, the solution of its initial value problem free of external excitation consists of two parts, the first one is the 'eigenfunction expansion' that is similar to the case without fractional-order derivative, and the second one is a definite integral that is independent of the eigenvalues (or characteristic roots. The integral disappears in the classical linear oscillator and it can be neglected from the solution when stationary solution is addressed. Moreover, the response of the fractionally damped oscillator under harmonic excitation is calculated in a similar way, and it is found that the fractional damping with order between 1 and 2 can be used to produce oscillation with large amplitude as well as to suppress oscillation, depending on the ratio of the excitation frequency and the natural frequency.
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.
Directory of Open Access Journals (Sweden)
T. H. S. Abdelaziz
2005-01-01
Full Text Available In this paper we introduce a complete parametric approach for solving the problem of eigenstructure assignment via state-derivative feedback for linear systems. This problem is always solvable for any controllable systems iff the open-loop system matrix is nonsingular. In this work, two parametric solutions to the feedback gain matrix are introduced that describe the available degrees of freedom offered by the state-derivative feedback in selecting the associated eigenvectors from an admissible class. These freedoms can be utilized to improve robustness of the closed-loop system. Accordingly, the sensitivity of the assigned eigenvalues to perturbations in the system and gain matrix is minimized. Numerical examples are included to show the effectiveness of the proposed approach.
International Nuclear Information System (INIS)
Dehghan, Mehdi; Shakourifar, Mohammad; Hamidi, Asgar
2009-01-01
The purpose of this study is to implement Adomian-Pade (Modified Adomian-Pade) technique, which is a combination of Adomian decomposition method (Modified Adomian decomposition method) and Pade approximation, for solving linear and nonlinear systems of Volterra functional equations. The results obtained by using Adomian-Pade (Modified Adomian-Pade) technique, are compared to those obtained by using Adomian decomposition method (Modified Adomian decomposition method) alone. The numerical results, demonstrate that ADM-PADE (MADM-PADE) technique, gives the approximate solution with faster convergence rate and higher accuracy than using the standard ADM (MADM).
A comparative study of iterative solutions to linear systems arising in quantum mechanics
International Nuclear Information System (INIS)
Jing Yanfei; Huang Tingzhu; Duan Yong; Carpentieri, Bruno
2010-01-01
This study is mainly focused on iterative solutions with simple diagonal preconditioning to two complex-valued nonsymmetric systems of linear equations arising from a computational chemistry model problem proposed by Sherry Li of NERSC. Numerical experiments show the feasibility of iterative methods to some extent when applied to the problems and reveal the competitiveness of our recently proposed Lanczos biconjugate A-orthonormalization methods to other classic and popular iterative methods. By the way, experiment results also indicate that application specific preconditioners may be mandatory and required for accelerating convergence.
Solution of single linear tridiagonal systems and vectorization of the ICCG algorithm on the Cray 1
International Nuclear Information System (INIS)
Kershaw, D.S.
1981-01-01
The numerical algorithms used to solve the physics equation in codes which model laser fusion are examined, it is found that a large number of subroutines require the solution of tridiagonal linear systems of equations. One dimensional radiation transport, thermal and suprathermal electron transport, ion thermal conduction, charged particle and neutron transport, all require the solution of tridiagonal systems of equations. The standard algorithm that has been used in the past on CDC 7600's will not vectorize and so cannot take advantage of the large speed increases possible on the Cray-1 through vectorization. There is however, an alternate algorithm for solving tridiagonal systems, called cyclic reduction, which allows for vectorization, and which is optimal for the Cray-1. Software based on this algorithm is now being used in LASNEX to solve tridiagonal linear systems in the subroutines mentioned above. The new algorithm runs as much as five times faster than the standard algorithm on the Cray-1. The ICCG method is being used to solve the diffusion equation with a nine-point coupling scheme on the CDC 7600. In going from the CDC 7600 to the Cray-1, a large part of the algorithm consists of solving tridiagonal linear systems on each L line of the Lagrangian mesh in a manner which is not vectorizable. An alternate ICCG algorithm for the Cray-1 was developed which utilizes a block form of the cyclic reduction algorithm. This new algorithm allows full vectorization and runs as much as five times faster than the old algorithm on the Cray-1. It is now being used in Cray LASNEX to solve the two-dimensional diffusion equation in all the physics subroutines mentioned above
On the completeness of the set of Bethe-Hulthen solutions of the linear Heisenberg system
International Nuclear Information System (INIS)
Caspers, W J; Labuz, M; Wal, A
2006-01-01
In this work we formulate the standard form of the solutions of the Heisenberg chain with periodic boundary conditions and show that these solutions can be transformed into the well-known Bethe-Hulthen solutions. The standard form is found by solving the secular problem, separated according to the irreducible representations of the translation group. The relevant parameters exp(ik j ) of the Bethe-Hulthen solutions are found from a set of linear equations with coefficients derived from the standard solutions. This correspondence between standard and Bethe-Hulthen solutions realizes the completeness of the Bethe-Hulthen method
DEFF Research Database (Denmark)
Köyluoglu, H.U.; Nielsen, Søren R.K.; Cakmak, A.S.
1994-01-01
perturbation method using stochastic differential equations. The joint statistical moments entering the perturbation solution are determined by considering an augmented dynamic system with state variables made up of the displacement and velocity vector and their first and second derivatives with respect......The paper deals with the first and second order statistical moments of the response of linear systems with random parameters subject to random excitation modelled as white-noise multiplied by an envelope function with random parameters. The method of analysis is basically a second order...... to the random parameters of the problem. Equations for partial derivatives are obtained from the partial differentiation of the equations of motion. The zero time-lag joint statistical moment equations for the augmented state vector are derived from the Itô differential formula. General formulation is given...
Some examples of non-linear systems and characteristics of their solutions
CSIR Research Space (South Africa)
Greben, JM
2006-07-01
Full Text Available . In contrast to certain other applications in complexity theory, these non-linear solutions are characterized by great stability. To go beyond the dominant non-perturbative solution one has to consider the source term as well. The parameter freedom...
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.
Iterative solution of linear systems in the 20th century
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
Direct interaction between linear electron transfer chains and solute transport systems in bacteria
Elferink, Marieke G.L.; Hellingwerf, Klaas J.; Belkum, Marco J. van; Poolman, Bert; Konings, Wil N.
1984-01-01
In studies on alanine and lactose transport in Rhodopseudomonas sphaeroides we have demonstrated that the rate of solute uptake in this phototrophic bacterium is regulated by the rate of light-induced cyclic electron transfer. In the present paper the interaction between linear electron transfer
Frank, T D; Beek, P J
2001-08-01
Recently, Küchler and Mensch [Stochastics Stochastics Rep. 40, 23 (1992)] derived exact stationary probability densities for linear stochastic delay differential equations. This paper presents an alternative derivation of these solutions by means of the Fokker-Planck approach introduced by Guillouzic [Phys. Rev. E 59, 3970 (1999); 61, 4906 (2000)]. Applications of this approach, which is argued to have greater generality, are discussed in the context of stochastic models for population growth and tracking movements.
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
LINPACK, Subroutine Library for Linear Equation System Solution and Matrix Calculation
International Nuclear Information System (INIS)
Dongarra, J.J.
1979-01-01
1 - Description of problem or function: LINPACK is a collection of FORTRAN subroutines which analyze and solve various classes of systems of simultaneous linear algebraic equations. The collection deals with general, banded, symmetric indefinite, symmetric positive definite, triangular, and tridiagonal square matrices, as well as with least squares problems and the QR and singular value decompositions of rectangular matrices. A subroutine-naming convention is employed in which each subroutine name consists of five letters which represent a coded specification (TXXYY) of the computation done by that subroutine. The first letter, T, indicates the matrix data type. Standard FORTRAN allows the use of three such types: S REAL, D DOUBLE PRECISION, and C COMPLEX. In addition, some FORTRAN systems allow a double-precision complex type: Z COMPLEX*16. The second and third letters of the subroutine name, XX, indicate the form of the matrix or its decomposition: GE: General, GB: General band, PO: Positive definite, PP: Positive definite packed, PB: Positive definite band, SI: Symmetric indefinite, SP: Symmetric indefinite packed, HI: Hermitian indefinite, HP: Hermitian indefinite packed, TR: Triangular, GT: General tridiagonal, PT: Positive definite tridiagonal, CH: Cholesky decomposition, QR: Orthogonal-triangular decomposition, SV: Singular value decomposition. The final two letters, YY, indicate the computation done by the particular subroutine: FA: Factor, CO: Factor and estimate condition, SL: Solve, DI: Determinant and/or inverse and/or inertia, DC: Decompose, UD: Update, DD: Down-date, EX Exchange. The following chart shows all the LINPACK subroutines. The initial 'S' in the names may be replaced by D, C or Z and the initial 'C' in the complex-only names may be replaced by a Z. SGE: FA, CO, SL, DI; SGB: FA, CO, SL, DI; SPO: FA, CO, SL, DI; SPP: FA, CO, SL, DI; SPB: FA, CO, SL, DI; SSI: FA, CO, SL, DI; SSP: FA, CO, SL, DI; CHI: FA, CO, SL, DI; CHP: FA, CO, SL, DI; STR
Directory of Open Access Journals (Sweden)
Dauda GuliburYAKUBU
2012-12-01
Full Text Available Accurate solutions to initial value systems of ordinary differential equations may be approximated efficiently by Runge-Kutta methods or linear multistep methods. Each of these has limitations of one sort or another. In this paper we consider, as a middle ground, the derivation of continuous general linear methods for solution of stiff systems of initial value problems in ordinary differential equations. These methods are designed to combine the advantages of both Runge-Kutta and linear multistep methods. Particularly, methods possessing the property of A-stability are identified as promising methods within this large class of general linear methods. We show that the continuous general linear methods are self-starting and have more ability to solve the stiff systems of ordinary differential equations, than the discrete ones. The initial value systems of ordinary differential equations are solved, for instance, without looking for any other method to start the integration process. This desirable feature of the proposed approach leads to obtaining very high accuracy of the solution of the given problem. Illustrative examples are given to demonstrate the novelty and reliability of the methods.
A block Krylov subspace time-exact solution method for linear ordinary differential equation systems
Bochev, Mikhail A.
2013-01-01
We propose a time-exact Krylov-subspace-based method for solving linear ordinary differential equation systems of the form $y'=-Ay+g(t)$ and $y"=-Ay+g(t)$, where $y(t)$ is the unknown function. The method consists of two stages. The first stage is an accurate piecewise polynomial approximation of
International Nuclear Information System (INIS)
Kalmykov, Mikhail Yu.; Kniehl, Bernd A.
2012-05-01
We argue that the Mellin-Barnes representations of Feynman diagrams can be used for obtaining linear systems of homogeneous differential equations for the original Feynman diagrams with arbitrary powers of propagators without recourse to the integration-by-parts technique. These systems of differential equation can be used (i) for the differential reductions to sets of basic functions and (ii) for counting the numbers of master-integrals.
Energy Technology Data Exchange (ETDEWEB)
Nygaard, K
1967-12-15
The numerical deconvolution of spectra is equivalent to the solution of a (large) system of linear equations with a matrix which is not necessarily a square matrix. The demand that the square sum of the residual errors shall be minimum is not in general sufficient to ensure a unique or 'sound' solution. Therefore other demands which may include the demand for minimum square errors are introduced which lead to 'sound' and 'non-oscillatory' solutions irrespective of the shape of the original matrix and of the determinant of the matrix of the normal equations.
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.
The solutions of second-order linear differential systems with constant delays
Diblík, Josef; Svoboda, Zdeněk
2017-07-01
The representations of solutions to initial problems for non-homogenous n-dimensional second-order differential equations with delays x″(t )-2 A x'(t -τ )+(A2+B2)x (t -2 τ )=f (t ) by means of special matrix delayed functions are derived. Square matrices A and B are commuting and τ > 0. Derived representations use what is called a delayed exponential of a matrix and results generalize some of known results previously derived for homogenous systems.
Chen, Y.-M.; Koniges, A. E.; Anderson, D. V.
1989-10-01
The biconjugate gradient method (BCG) provides an attractive alternative to the usual conjugate gradient algorithms for the solution of sparse systems of linear equations with nonsymmetric and indefinite matrix operators. A preconditioned algorithm is given, whose form resembles the incomplete L-U conjugate gradient scheme (ILUCG2) previously presented. Although the BCG scheme requires the storage of two additional vectors, it converges in a significantly lesser number of iterations (often half), while the number of calculations per iteration remains essentially the same.
Iterative solution of general sparse linear systems on clusters of workstations
Energy Technology Data Exchange (ETDEWEB)
Lo, Gen-Ching; Saad, Y. [Univ. of Minnesota, Minneapolis, MN (United States)
1996-12-31
Solving sparse irregularly structured linear systems on parallel platforms poses several challenges. First, sparsity makes it difficult to exploit data locality, whether in a distributed or shared memory environment. A second, perhaps more serious challenge, is to find efficient ways to precondition the system. Preconditioning techniques which have a large degree of parallelism, such as multicolor SSOR, often have a slower rate of convergence than their sequential counterparts. Finally, a number of other computational kernels such as inner products could ruin any gains gained from parallel speed-ups, and this is especially true on workstation clusters where start-up times may be high. In this paper we discuss these issues and report on our experience with PSPARSLIB, an on-going project for building a library of parallel iterative sparse matrix solvers.
Linear superposition solutions to nonlinear wave equations
International Nuclear Information System (INIS)
Liu Yu
2012-01-01
The solutions to a linear wave equation can satisfy the principle of superposition, i.e., the linear superposition of two or more known solutions is still a solution of the linear wave equation. We show in this article that many nonlinear wave equations possess exact traveling wave solutions involving hyperbolic, triangle, and exponential functions, and the suitable linear combinations of these known solutions can also constitute linear superposition solutions to some nonlinear wave equations with special structural characteristics. The linear superposition solutions to the generalized KdV equation K(2,2,1), the Oliver water wave equation, and the k(n, n) equation are given. The structure characteristic of the nonlinear wave equations having linear superposition solutions is analyzed, and the reason why the solutions with the forms of hyperbolic, triangle, and exponential functions can form the linear superposition solutions is also discussed
Energy Technology Data Exchange (ETDEWEB)
Clemens, M.; Weiland, T. [Technische Hochschule Darmstadt (Germany)
1996-12-31
In the field of computational electrodynamics the discretization of Maxwell`s equations using the Finite Integration Theory (FIT) yields very large, sparse, complex symmetric linear systems of equations. For this class of complex non-Hermitian systems a number of conjugate gradient-type algorithms is considered. The complex version of the biconjugate gradient (BiCG) method by Jacobs can be extended to a whole class of methods for complex-symmetric algorithms SCBiCG(T, n), which only require one matrix vector multiplication per iteration step. In this class the well-known conjugate orthogonal conjugate gradient (COCG) method for complex-symmetric systems corresponds to the case n = 0. The case n = 1 yields the BiCGCR method which corresponds to the conjugate residual algorithm for the real-valued case. These methods in combination with a minimal residual smoothing process are applied separately to practical 3D electro-quasistatical and eddy-current problems in electrodynamics. The practical performance of the SCBiCG methods is compared with other methods such as QMR and TFQMR.
Communications oriented programming of parallel iterative solutions of sparse linear systems
Patrick, M. L.; Pratt, T. W.
1986-01-01
Parallel algorithms are developed for a class of scientific computational problems by partitioning the problems into smaller problems which may be solved concurrently. The effectiveness of the resulting parallel solutions is determined by the amount and frequency of communication and synchronization and the extent to which communication can be overlapped with computation. Three different parallel algorithms for solving the same class of problems are presented, and their effectiveness is analyzed from this point of view. The algorithms are programmed using a new programming environment. Run-time statistics and experience obtained from the execution of these programs assist in measuring the effectiveness of these algorithms.
Energy Technology Data Exchange (ETDEWEB)
Nygaard, K
1968-09-15
From the point of view that no mathematical method can ever minimise or alter errors already made in a physical measurement, the classical least squares method has severe limitations which makes it unsuitable for the statistical analysis of many physical measurements. Based on the assumptions that the experimental errors are characteristic for each single experiment and that the errors must be properly estimated rather than minimised, a new method for solving large systems of linear equations is developed. The new method exposes the entire range of possible solutions before the decision is taken which of the possible solutions should be chosen as a representative one. The choice is based on physical considerations which (in two examples, curve fitting and unfolding of a spectrum) are presented in such a form that a computer is able to make the decision, A description of the computation is given. The method described is a tool for removing uncertainties due to conventional mathematical formulations (zero determinant, linear dependence) and which are not inherent in the physical problem as such. The method is therefore especially well fitted for unfolding of spectra.
International Nuclear Information System (INIS)
Nygaard, K.
1968-09-01
From the point of view that no mathematical method can ever minimise or alter errors already made in a physical measurement, the classical least squares method has severe limitations which makes it unsuitable for the statistical analysis of many physical measurements. Based on the assumptions that the experimental errors are characteristic for each single experiment and that the errors must be properly estimated rather than minimised, a new method for solving large systems of linear equations is developed. The new method exposes the entire range of possible solutions before the decision is taken which of the possible solutions should be chosen as a representative one. The choice is based on physical considerations which (in two examples, curve fitting and unfolding of a spectrum) are presented in such a form that a computer is able to make the decision, A description of the computation is given. The method described is a tool for removing uncertainties due to conventional mathematical formulations (zero determinant, linear dependence) and which are not inherent in the physical problem as such. The method is therefore especially well fitted for unfolding of spectra
Energy Technology Data Exchange (ETDEWEB)
Alleon, G. [EADS-CCR, 31 - Blagnac (France); Carpentieri, B.; Du, I.S.; Giraud, L.; Langou, J.; Martin, E. [Cerfacs, 31 - Toulouse (France)
2003-07-01
The boundary element method has become a popular tool for the solution of Maxwell's equations in electromagnetism. It discretizes only the surface of the radiating object and gives rise to linear systems that are smaller in size compared to those arising from finite element or finite difference discretizations. However, these systems are prohibitively demanding in terms of memory for direct methods and challenging to solve by iterative methods. In this paper we address the iterative solution via preconditioned Krylov methods of electromagnetic scattering problems expressed in an integral formulation, with main focus on the design of the pre-conditioner. We consider an approximate inverse method based on the Frobenius-norm minimization with a pattern prescribed in advance. The pre-conditioner is constructed from a sparse approximation of the dense coefficient matrix, and the patterns both for the pre-conditioner and for the coefficient matrix are computed a priori using geometric information from the mesh. We describe the implementation of the approximate inverse in an out-of-core parallel code that uses multipole techniques for the matrix-vector products, and show results on the numerical scalability of our method on systems of size up to one million unknowns. We propose an embedded iterative scheme based on the GMRES method and combined with multipole techniques, aimed at improving the robustness of the approximate inverse for large problems. We prove by numerical experiments that the proposed scheme enables the solution of very large and difficult problems efficiently at reduced computational and memory cost. Finally we perform a preliminary study on a spectral two-level pre-conditioner to enhance the robustness of our method. This numerical technique exploits spectral information of the preconditioned systems to build a low rank-update of the pre-conditioner. (authors)
International Nuclear Information System (INIS)
Alleon, G.; Carpentieri, B.; Du, I.S.; Giraud, L.; Langou, J.; Martin, E.
2003-01-01
The boundary element method has become a popular tool for the solution of Maxwell's equations in electromagnetism. It discretizes only the surface of the radiating object and gives rise to linear systems that are smaller in size compared to those arising from finite element or finite difference discretizations. However, these systems are prohibitively demanding in terms of memory for direct methods and challenging to solve by iterative methods. In this paper we address the iterative solution via preconditioned Krylov methods of electromagnetic scattering problems expressed in an integral formulation, with main focus on the design of the pre-conditioner. We consider an approximate inverse method based on the Frobenius-norm minimization with a pattern prescribed in advance. The pre-conditioner is constructed from a sparse approximation of the dense coefficient matrix, and the patterns both for the pre-conditioner and for the coefficient matrix are computed a priori using geometric information from the mesh. We describe the implementation of the approximate inverse in an out-of-core parallel code that uses multipole techniques for the matrix-vector products, and show results on the numerical scalability of our method on systems of size up to one million unknowns. We propose an embedded iterative scheme based on the GMRES method and combined with multipole techniques, aimed at improving the robustness of the approximate inverse for large problems. We prove by numerical experiments that the proposed scheme enables the solution of very large and difficult problems efficiently at reduced computational and memory cost. Finally we perform a preliminary study on a spectral two-level pre-conditioner to enhance the robustness of our method. This numerical technique exploits spectral information of the preconditioned systems to build a low rank-update of the pre-conditioner. (authors)
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.
The solution of linear systems of equations with a structural analysis code on the NAS CRAY-2
Poole, Eugene L.; Overman, Andrea L.
1988-01-01
Two methods for solving linear systems of equations on the NAS Cray-2 are described. One is a direct method; the other is an iterative method. Both methods exploit the architecture of the Cray-2, particularly the vectorization, and are aimed at structural analysis applications. To demonstrate and evaluate the methods, they were installed in a finite element structural analysis code denoted the Computational Structural Mechanics (CSM) Testbed. A description of the techniques used to integrate the two solvers into the Testbed is given. Storage schemes, memory requirements, operation counts, and reformatting procedures are discussed. Finally, results from the new methods are compared with results from the initial Testbed sparse Choleski equation solver for three structural analysis problems. The new direct solvers described achieve the highest computational rates of the methods compared. The new iterative methods are not able to achieve as high computation rates as the vectorized direct solvers but are best for well conditioned problems which require fewer iterations to converge to the solution.
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.
Czech Academy of Sciences Publication Activity Database
Dosoudilová, M.; Lomtatidze, Alexander; Šremr, Jiří
2015-01-01
Roč. 120, June (2015), s. 57-75 ISSN 0362-546X Institutional support: RVO:67985840 Keywords : half-linear system * Hartman -Wintner theorem * Kamenev theorem Subject RIV: BA - General Mathematics Impact factor: 1.125, year: 2015 http://www.sciencedirect.com/science/article/pii/S0362546X15000620
Subroutine for series solutions of linear differential equations
International Nuclear Information System (INIS)
Tasso, H.; Steuerwald, J.
1976-02-01
A subroutine for Taylor series solutions of systems of ordinary linear differential equations is descriebed. It uses the old idea of Lie series but allows simple implementation and is time-saving for symbolic manipulations. (orig.) [de
International Nuclear Information System (INIS)
Valat, J.
1960-12-01
Universal stability diagrams have been calculated and experimentally checked for Hill-Meissner type equations with square-wave coefficients. The study of these equations in the phase-plane has then made it possible to extend the periodic solution calculations to the case of non-linear differential equations with periodic square-wave coefficients. This theory has been checked experimentally. For non-linear coupled systems with constant coefficients, a search was first made for solutions giving an algebraic motion. The elliptical and Fuchs's functions solve such motions. The study of non-algebraic motions is more delicate, apart from the study of nonlinear Lissajous's motions. A functional analysis shows that it is possible however in certain cases to decouple the system and to find general solutions. For non-linear coupled systems with periodic square-wave coefficients it is then possible to calculate the conditions leading to periodic solutions, if the two non-linear associated systems with constant coefficients fall into one of the categories of the above paragraph. (author) [fr
General solution of linear vector supersymmetry
International Nuclear Information System (INIS)
Blasi, Alberto; Maggiore, Nicola
2007-01-01
We give the general solution of the Ward identity for the linear vector supersymmetry which characterizes all topological models. Such a solution, whose expression is quite compact and simple, greatly simplifies the study of theories displaying a supersymmetric algebraic structure, reducing to a few lines the proof of their possible finiteness. In particular, the cohomology technology, usually involved for the quantum extension of these theories, is completely bypassed. The case of Chern-Simons theory is taken as an example
Energy Technology Data Exchange (ETDEWEB)
Joubert, W. [Los Alamos National Lab., NM (United States); Carey, G.F. [Univ. of Texas, Austin, TX (United States)
1994-12-31
A great need exists for high performance numerical software libraries transportable across parallel machines. This talk concerns the PCG package, which solves systems of linear equations by iterative methods on parallel computers. The features of the package are discussed, as well as techniques used to obtain high performance as well as transportability across architectures. Representative numerical results are presented for several machines including the Connection Machine CM-5, Intel Paragon and Cray T3D parallel computers.
DEFF Research Database (Denmark)
Micaletti, R. C.; Cakmak, A. S.; Nielsen, Søren R. K.
structural properties. The resulting state-space formulation is a system of ordinary stochastic differential equations with random coefficient and deterministic initial conditions which are subsequently transformed into ordinary stochastic differential equations with deterministic coefficients and random......A method for computing the lower-order moments of randomly-excited multi-degree-of-freedom (MDOF) systems with random structural properties is proposed. The method is grounded in the techniques of stochastic calculus, utilizing a Markov diffusion process to model the structural system with random...... initial conditions. This transformation facilitates the derivation of differential equations which govern the evolution of the unconditional statistical moments of response. Primary consideration is given to linear systems and systems with odd polynomial nonlinearities, for in these cases...
Analytical exact solution of the non-linear Schroedinger equation
International Nuclear Information System (INIS)
Martins, Alisson Xavier; Rocha Filho, Tarcisio Marciano da
2011-01-01
Full text: In this work we present how to classify and obtain analytical solutions of the Schroedinger equation with a generic non-linearity in 1+1 dimensions. Our approach is based on the determination of Lie symmetry transformation mapping solutions into solutions, and non-classical symmetry transformations, mapping a given solution into itself. From these symmetries it is then possible to reduce the equation to a system of ordinary differential equations which can then be solved using standard methods. The generic non-linearity is handled by considering it as an additional unknown in the determining equations for the symmetry transformations. This results in an over-determined system of non-linear partial differential equations. Its solution can then be determined in some cases by reducing it to the so called involutive (triangular) form, and then solved. This reduction is very tedious and can only performed using a computer algebra system. Once the determining system is solved, we obtain the explicit form for the non-linearity admitting a Lie or non-classical symmetry. The analytical solutions are then derived by solving the reduced ordinary differential equations. The non-linear determining system for the non-classical symmetry transformations and Lie symmetry generators are obtaining using the computer algebra package SADE (symmetry analysis of differential equations), developed at our group. (author)
Linear operator inequalities for strongly stable weakly regular linear systems
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
Mercier, Sylvain; Gratton, Serge; Tardieu, Nicolas; Vasseur, Xavier
2017-12-01
Many applications in structural mechanics require the numerical solution of sequences of linear systems typically issued from a finite element discretization of the governing equations on fine meshes. The method of Lagrange multipliers is often used to take into account mechanical constraints. The resulting matrices then exhibit a saddle point structure and the iterative solution of such preconditioned linear systems is considered as challenging. A popular strategy is then to combine preconditioning and deflation to yield an efficient method. We propose an alternative that is applicable to the general case and not only to matrices with a saddle point structure. In this approach, we consider to update an existing algebraic or application-based preconditioner, using specific available information exploiting the knowledge of an approximate invariant subspace or of matrix-vector products. The resulting preconditioner has the form of a limited memory quasi-Newton matrix and requires a small number of linearly independent vectors. Numerical experiments performed on three large-scale applications in elasticity highlight the relevance of the new approach. We show that the proposed method outperforms the deflation method when considering sequences of linear systems with varying matrices.
International Nuclear Information System (INIS)
Mathews, M.D.; Ambekar, B.R.; Tyagi, A.K.
2005-01-01
Cell parameters and linear thermal expansion studies of the Th-M oxide systems with general compositions Th 1-x M x O 2-x/2 (M=Eu 3+ , Gd 3+ and Dy 3+ , 0.0= 1.5 in the ThO 2 lattice. The upper solid solubility limits of EuO 1.5 , GdO 1.5 and DyO 1.5 in the ThO 2 lattice under conditions of slow cooling from 1673K are represented as Th 0.50 Eu 0.50 O 1.75 , Th 0.60 Gd 0.40 O 1.80 and Th 0.85 Dy 0.15 O 1.925 , respectively. The linear thermal expansion (293-1123K) of MO 1.5 and their single-phase solid solutions with thoria were investigated by dilatometery. The average linear thermal expansion coefficients (α-bar ) of the compounds decrease on going from EuO 1.5 to DyO 1.5 . The values of α-bar for EuO 1.5 , GdO 1.5 and DyO 1.5 containing solid solutions showed a downward trend as a function of the dopant concentration. The linear thermal expansion (293-1473K) of the solid solutions investigated by high-temperature XRD also showed a similar trend
Numerical Solution of Heun Equation Via Linear Stochastic Differential Equation
Directory of Open Access Journals (Sweden)
Hamidreza Rezazadeh
2014-05-01
Full Text Available In this paper, we intend to solve special kind of ordinary differential equations which is called Heun equations, by converting to a corresponding stochastic differential equation(S.D.E.. So, we construct a stochastic linear equation system from this equation which its solution is based on computing fundamental matrix of this system and then, this S.D.E. is solved by numerically methods. Moreover, its asymptotic stability and statistical concepts like expectation and variance of solutions are discussed. Finally, the attained solutions of these S.D.E.s compared with exact solution of corresponding differential equations.
International Nuclear Information System (INIS)
Walrand, Stephan; Jamar, François; Pauwels, Stanislas
2009-01-01
Ill-posed linear systems occur in many different fields. A class of regularization methods, called constrained optimization, aims to determine the extremum of a penalty function whilst constraining an objective function to a likely value. We propose here a novel heuristic way to screen the local extrema satisfying the discrepancy principle. A modified version of the Landweber algorithm is used for the iteration process. After finding a local extremum, a bound is performed to the 'farthest' estimate in the data space still satisfying the discrepancy principle. Afterwards, the modified Landweber algorithm is again applied to find a new local extremum. This bound-iteration process is repeated until a satisfying solution is reached. For evaluation in nuclear medicine tomography, a novel penalty function that preserves the edge steps in the reconstructed solution was evaluated on Monte Carlo simulations and using real SPECT acquisitions as well. Surprisingly, the first bound always provided a significantly better solution in a wide range of statistics
DEFF Research Database (Denmark)
Gersborg-Hansen, Allan; Dammann, Bernd; Aage, Niels
concerned with developing a proper (COMSOL) model rather than developing efficient linear algebra solvers which motivates this investigation of the efficiency of the coupling COMSOL + SPL. The technicalities of making such a coupling is described in detail along with a measure of the speedup...
Energy Technology Data Exchange (ETDEWEB)
Valat, J [Commissariat a l' Energie Atomique, Saclay (France). Centre d' Etudes Nucleaires
1960-12-15
Universal stability diagrams have been calculated and experimentally checked for Hill-Meissner type equations with square-wave coefficients. The study of these equations in the phase-plane has then made it possible to extend the periodic solution calculations to the case of non-linear differential equations with periodic square-wave coefficients. This theory has been checked experimentally. For non-linear coupled systems with constant coefficients, a search was first made for solutions giving an algebraic motion. The elliptical and Fuchs's functions solve such motions. The study of non-algebraic motions is more delicate, apart from the study of nonlinear Lissajous's motions. A functional analysis shows that it is possible however in certain cases to decouple the system and to find general solutions. For non-linear coupled systems with periodic square-wave coefficients it is then possible to calculate the conditions leading to periodic solutions, if the two non-linear associated systems with constant coefficients fall into one of the categories of the above paragraph. (author) [French] Pour les equations du genre de Hill-Meissner a coefficients creneles, on a calcule des diagrammes universels de stabilite et ceux-ci ont ete verifies experimentalement. L'etude de ces equations dans le plan de phase a permis ensuite d'etendre le calcul des solutions periodiques au cas des equations differentielles non lineaires a coefficients periodiques creneles. Cette theorie a ete verifiee experimentalement. Pour Jes systemes couples non lineaires a coefficients constants, on a d'abord cherche les solutions menant a des mouvements algebriques. Les fonctions elliptiques et fuchsiennes uniformisent de tels mouvements. L'etude de mouvements non algebriques est plus delicate, a part l'etude des mouvements de Lissajous non lineaires. Une analyse fonctionnelle montre qu'il est toutefois possible dans certains cas de decoupler le systeme et de trouver des solutions generales. Pour les
On deformations of linear differential systems
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
Numerical solution of two-dimensional non-linear partial differential ...
African Journals Online (AJOL)
linear partial differential equations using a hybrid method. The solution technique involves discritizing the non-linear system of partial differential equations (PDEs) to obtain a corresponding nonlinear system of algebraic difference equations to be ...
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.
Samba, A. S.
1985-01-01
The problem of solving banded linear systems by direct (non-iterative) techniques on the Vector Processor System (VPS) 32 supercomputer is considered. Two efficient direct methods for solving banded linear systems on the VPS 32 are described. The vector cyclic reduction (VCR) algorithm is discussed in detail. The performance of the VCR on a three parameter model problem is also illustrated. The VCR is an adaptation of the conventional point cyclic reduction algorithm. The second direct method is the Customized Reduction of Augmented Triangles' (CRAT). CRAT has the dominant characteristics of an efficient VPS 32 algorithm. CRAT is tailored to the pipeline architecture of the VPS 32 and as a consequence the algorithm is implicitly vectorizable.
Dynamical systems and linear algebra
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)
Sartabanov, Zhaishylyk A.
2017-09-01
A new approach to the study of periodic by all independent variables system of equations with a differentiation operator solutions along the direction of the main diagonal and with delayed arguments is proposed. The essence of the approach is to reduce the study of the multi-periodic solution of a linear inhomogeneous system to the construction of a solution of a simpler linear differential-difference system on the basis of the method of variating arbitrary constants of the complete integral of a homogeneous system. An integral representation of the unique multiperiodic solution of an inhomogeneous system is presented, expressed by a functional series of terms given by multiple repeated integrals. An estimate is given for the norm of a multi-periodic solution.
Asymptotic solutions and spectral theory of linear wave equations
International Nuclear Information System (INIS)
Adam, J.A.
1982-01-01
This review contains two closely related strands. Firstly the asymptotic solution of systems of linear partial differential equations is discussed, with particular reference to Lighthill's method for obtaining the asymptotic functional form of the solution of a scalar wave equation with constant coefficients. Many of the applications of this technique are highlighted. Secondly, the methods and applications of the theory of the reduced (one-dimensional) wave equation - particularly spectral theory - are discussed. While the breadth of application and power of the techniques is emphasised throughout, the opportunity is taken to present to a wider readership, developments of the methods which have occured in some aspects of astrophysical (particularly solar) and geophysical fluid dynamics. It is believed that the topics contained herein may be of relevance to the applied mathematician or theoretical physicist interest in problems of linear wave propagation in these areas. (orig./HSI)
The linear ordering problem: an algorithm for the optimal solution ...
African Journals Online (AJOL)
In this paper we describe and implement an algorithm for the exact solution of the Linear Ordering problem. Linear Ordering is the problem of finding a linear order of the nodes of a graph such that the sum of the weights which are consistent with this order is as large as possible. It is an NP - Hard combinatorial optimisation ...
Linear perturbations of a self-similar solution of hydrodynamics with non-linear heat conduction
International Nuclear Information System (INIS)
Dubois-Boudesocque, Carine
2000-01-01
The stability of an ablative flow, where a shock wave is located upstream a thermal front, is of importance in inertial confinement fusion. The present model considers an exact self-similar solution to the hydrodynamic equations with non-linear heat conduction for a semi-infinite slab. For lack of an analytical solution, a high resolution numerical procedure is devised, which couples a finite difference method with a relaxation algorithm using a two-domain pseudo-spectral method. Stability of this solution is studied by introducing linear perturbation method within a Lagrangian-Eulerian framework. The initial and boundary value problem is solved by a splitting of the equations between a hyperbolic system and a parabolic equation. The boundary conditions of the hyperbolic system are treated, in the case of spectral methods, according to Thompson's approach. The parabolic equation is solved by an influence matrix method. These numerical procedures have been tested versus exact solutions. Considering a boundary heat flux perturbation, the space-time evolution of density, velocity and temperature are shown. (author) [fr
Solution of linear ill-posed problems using overcomplete dictionaries
Pensky, Marianna
2016-01-01
In the present paper we consider application of overcomplete dictionaries to solution of general ill-posed linear inverse problems. Construction of an adaptive optimal solution for such problems usually relies either on a singular value decomposition or representation of the solution via an orthonormal basis. The shortcoming of both approaches lies in the fact that, in many situations, neither the eigenbasis of the linear operator nor a standard orthonormal basis constitutes an appropriate co...
Introduction to linear systems of differential equations
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...
Feedback systems for linear colliders
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...
Solving Fully Fuzzy Linear System of Equations in General Form
Directory of Open Access Journals (Sweden)
A. Yousefzadeh
2012-06-01
Full Text Available In this work, we propose an approach for computing the positive solution of a fully fuzzy linear system where the coefficient matrix is a fuzzy $nimes n$ matrix. To do this, we use arithmetic operations on fuzzy numbers that introduced by Kaffman in and convert the fully fuzzy linear system into two $nimes n$ and $2nimes 2n$ crisp linear systems. If the solutions of these linear systems don't satisfy in positive fuzzy solution condition, we introduce the constrained least squares problem to obtain optimal fuzzy vector solution by applying the ranking function in given fully fuzzy linear system. Using our proposed method, the fully fuzzy linear system of equations always has a solution. Finally, we illustrate the efficiency of proposed method by solving some numerical examples.
Rational approximations to solutions of linear differential equations.
Chudnovsky, D V; Chudnovsky, G V
1983-08-01
Rational approximations of Padé and Padé type to solutions of differential equations are considered. One of the main results is a theorem stating that a simultaneous approximation to arbitrary solutions of linear differential equations over C(x) cannot be "better" than trivial ones implied by the Dirichlet box principle. This constitutes, in particular, the solution in the linear case of Kolchin's problem that the "Roth's theorem" holds for arbitrary solutions of algebraic differential equations. Complete effective proofs for several valuations are presented based on the Wronskian methods and graded subrings of Picard-Vessiot extensions.
Kibler, K. S.; Mcdaniel, G. A.
1981-01-01
A digital local linearization technique was used to solve a system of stiff differential equations which simulate a magnetic bearing assembly. The results prove the technique to be accurate, stable, and efficient when compared to a general purpose variable order Adams method with a stiff option.
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...
Approximate solution to neutron transport equation with linear anisotropic scattering
International Nuclear Information System (INIS)
Coppa, G.; Ravetto, P.; Sumini, M.
1983-01-01
A method to obtain an approximate solution to the transport equation, when both sources and collisions show a linearly anisotropic behavior, is outlined and the possible implications for numerical calculations in applied neutronics as well as shielding evaluations are investigated. The form of the differential system of equations taken by the method is quite handy and looks simpler and more manageable than any other today available technique. To go deeper into the efficiency of the method, some typical calculations concerning critical dimension of multiplying systems are then performed and the results are compared with the ones coming from the classical Ssub(N) approximations. The outcome of such calculations leads us to think of interesting developments of the method which could be quite useful in alternative to other today widespread approximate procedures, for any geometry, but especially for curved ones. (author)
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
Oscillatory behaviour of solutions of linear neutral differential ...
African Journals Online (AJOL)
The paper considers the contribution of space-time noise to the oscillatory behaviour of solutions of a linear neutral stochastic delay differential equation. It was established that under certain conditions on the time lags and their speed of adjustments, the presence of noise generates oscillation in the solution of the equation ...
Exponential estimates for solutions of half-linear differential equations
Czech Academy of Sciences Publication Activity Database
Řehák, Pavel
2015-01-01
Roč. 147, č. 1 (2015), s. 158-171 ISSN 0236-5294 Institutional support: RVO:67985840 Keywords : half-linear differential equation * decreasing solution * increasing solution * asymptotic behavior Subject RIV: BA - General Mathematics Impact factor: 0.469, year: 2015 http://link.springer.com/article/10.1007%2Fs10474-015-0522-9
An algorithm for the solution of dynamic linear programs
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
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
Oscillatory solutions of the Cauchy problem for linear differential equations
Directory of Open Access Journals (Sweden)
Gro Hovhannisyan
2015-06-01
Full Text Available We consider the Cauchy problem for second and third order linear differential equations with constant complex coefficients. We describe necessary and sufficient conditions on the data for the existence of oscillatory solutions. It is known that in the case of real coefficients the oscillatory behavior of solutions does not depend on initial values, but we show that this is no longer true in the complex case: hence in practice it is possible to control oscillatory behavior by varying the initial conditions. Our Proofs are based on asymptotic analysis of the zeros of solutions, represented as linear combinations of exponential functions.
Johnson, Thomas
2018-01-01
In a recent seminal paper \\cite{D--H--R} of Dafermos, Holzegel and Rodnianski the linear stability of the Schwarzschild family of black hole solutions to the Einstein vacuum equations was established by imposing a double null gauge. In this paper we shall prove that the Schwarzschild family is linearly stable as solutions to the Einstein vacuum equations by imposing instead a generalised wave gauge: all sufficiently regular solutions to the system of equations that result from linearising the...
Optimal Control of Switching Linear Systems
Directory of Open Access Journals (Sweden)
Ali Benmerzouga
2004-06-01
Full Text Available A solution to the control of switching linear systems with input constraints was given in Benmerzouga (1997 for both the conventional enumeration approach and the new approach. The solution given there turned out to be not unique. The main objective in this work is to determine the optimal control sequences {Ui(k , i = 1,..., M ; k = 0, 1, ..., N -1} which transfer the system from a given initial state X0 to a specific target state XT (or to be as close as possible by using the same discrete time solution obtained in Benmerzouga (1997 and minimizing a running cost-to-go function. By using the dynamic programming technique, the optimal solution is found for both approaches given in Benmerzouga (1997. The computational complexity of the modified algorithm is also given.
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
Systems of Inhomogeneous Linear Equations
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.
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
Iterative solution of linear equations in ODE codes. [Krylov subspaces
Energy Technology Data Exchange (ETDEWEB)
Gear, C. W.; Saad, Y.
1981-01-01
Each integration step of a stiff equation involves the solution of a nonlinear equation, usually by a quasi-Newton method that leads to a set of linear problems. Iterative methods for these linear equations are studied. Of particular interest are methods that do not require an explicit Jacobian, but can work directly with differences of function values using J congruent to f(x + delta) - f(x). Some numerical experiments using a modification of LSODE are reported. 1 figure, 2 tables.
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.)
Exact solution of some linear matrix equations using algebraic methods
Djaferis, T. E.; Mitter, S. K.
1977-01-01
A study is done of solution methods for Linear Matrix Equations including Lyapunov's equation, using methods of modern algebra. The emphasis is on the use of finite algebraic procedures which are easily implemented on a digital computer and which lead to an explicit solution to the problem. The action f sub BA is introduced a Basic Lemma is proven. The equation PA + BP = -C as well as the Lyapunov equation are analyzed. Algorithms are given for the solution of the Lyapunov and comment is given on its arithmetic complexity. The equation P - A'PA = Q is studied and numerical examples are given.
Numerical solution of non-linear diffusion problems
International Nuclear Information System (INIS)
Carmen, A. del; Ferreri, J.C.
1998-01-01
This paper presents a method for the numerical solution of non-linear diffusion problems using finite-differences in moving grids. Due to the presence of steep fronts in the solution domain and to the presence of advective terms originating in the grid movement, an implicit TVD scheme, first order in time and second order in space has been developed. Some algebraic details of the derivation are given. Results are shown for the pure advection of a scalar as a test case and an example dealing with the slow spreading of viscous fluids over plane surfaces. The agreement between numerical and analytical solutions is excellent. (author). 8 refs., 3 figs
Asymptotic formulae for solutions of half-linear differential equations
Czech Academy of Sciences Publication Activity Database
Řehák, Pavel
2017-01-01
Roč. 292, January (2017), s. 165-177 ISSN 0096-3003 Institutional support: RVO:67985840 Keywords : half-linear differential equation * nonoscillatory solution * regular variation Subject RIV: BA - General Mathematics OBOR OECD: Applied mathematics Impact factor: 1.738, year: 2016 http://www.sciencedirect.com/science/article/pii/S0096300316304581
On nonnegative solutions of second order linear functional differential equations
Czech Academy of Sciences Publication Activity Database
Lomtatidze, Alexander; Vodstrčil, Petr
2004-01-01
Roč. 32, č. 1 (2004), s. 59-88 ISSN 1512-0015 Institutional research plan: CEZ:AV0Z1019905 Keywords : second order linear functional differential equations * nonnegative solution * two-point boundary value problem Subject RIV: BA - General Mathematics
Radial solutions to semilinear elliptic equations via linearized operators
Directory of Open Access Journals (Sweden)
Phuong Le
2017-04-01
Full Text Available Let $u$ be a classical solution of semilinear elliptic equations in a ball or an annulus in $\\mathbb{R}^N$ with zero Dirichlet boundary condition where the nonlinearity has a convex first derivative. In this note, we prove that if the $N$-th eigenvalue of the linearized operator at $u$ is positive, then $u$ must be radially symmetric.
Minimal solution of linear formed fuzzy matrix equations
Directory of Open Access Journals (Sweden)
Maryam Mosleh
2012-10-01
Full Text Available In this paper according to the structured element method, the $mimes n$ inconsistent fuzzy matrix equation $Ailde{X}=ilde{B},$ which are linear formed by fuzzy structured element, is investigated. The necessary and sufficient condition for the existence of a fuzzy solution is also discussed. some examples are presented to illustrate the proposed method.
Linear systems optimal and robust control
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...
Identification problems in linear transformation system
International Nuclear Information System (INIS)
Delforge, Jacques.
1975-01-01
An attempt was made to solve the theoretical and numerical difficulties involved in the identification problem relative to the linear part of P. Delattre's theory of transformation systems. The theoretical difficulties are due to the very important problem of the uniqueness of the solution, which must be demonstrated in order to justify the value of the solution found. Simple criteria have been found when measurements are possible on all the equivalence classes, but the problem remains imperfectly solved when certain evolution curves are unknown. The numerical difficulties are of two kinds: a slow convergence of iterative methods and a strong repercussion of numerical and experimental errors on the solution. In the former case a fast convergence was obtained by transformation of the parametric space, while in the latter it was possible, from sensitivity functions, to estimate the errors, to define and measure the conditioning of the identification problem then to minimize this conditioning as a function of the experimental conditions [fr
A Solution to the Fundamental Linear Fractional Order Differential Equation
Hartley, Tom T.; Lorenzo, Carl F.
1998-01-01
This paper provides a solution to the fundamental linear fractional order differential equation, namely, (sub c)d(sup q, sub t) + ax(t) = bu(t). The impulse response solution is shown to be a series, named the F-function, which generalizes the normal exponential function. The F-function provides the basis for a qth order "fractional pole". Complex plane behavior is elucidated and a simple example, the inductor terminated semi- infinite lossy line, is used to demonstrate the theory.
A general method for enclosing solutions of interval linear equations
Czech Academy of Sciences Publication Activity Database
Rohn, Jiří
2012-01-01
Roč. 6, č. 4 (2012), s. 709-717 ISSN 1862-4472 R&D Projects: GA ČR GA201/09/1957; GA ČR GC201/08/J020 Institutional research plan: CEZ:AV0Z10300504 Keywords : interval linear equations * solution set * enclosure * absolute value inequality Subject RIV: BA - General Mathematics Impact factor: 1.654, year: 2012
Georgatzis, Konstantinos; Lal, Partha; Hawthorne, Christopher; Shaw, Martin; Piper, Ian; Tarbert, Claire; Donald, Rob; Williams, Christopher K I
2016-01-01
High-resolution, artefact-free and accurately annotated physiological data are desirable in patients with brain injury both to inform clinical decision-making and for intelligent analysis of the data in applications such as predictive modelling. We have quantified the quality of annotation surrounding artefactual events and propose a factorial switching linear dynamical systems (FSLDS) approach to automatically detect artefact in physiological data collected in the neurological intensive care unit (NICU). Retrospective analysis of the BrainIT data set to discover potential hypotensive events corrupted by artefact and identify the annotation of associated clinical interventions. Training of an FSLDS model on clinician-annotated artefactual events in five patients with severe traumatic brain injury. In a subset of 187 patients in the BrainIT database, 26.5 % of potential hypotensive events were abandoned because of artefactual data. Only 30 % of these episodes could be attributed to an annotated clinical intervention. As assessed by the area under the receiver operating characteristic curve metric, FSLDS model performance in automatically identifying the events of blood sampling, arterial line damping and patient handling was 0.978, 0.987 and 0.765, respectively. The influence of artefact on physiological data collected in the NICU is a significant problem. This pilot study using an FSLDS approach shows real promise and is under further development.
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.
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.
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.
An investigation on the solutions for the linear inverse problem in gamma ray tomography
International Nuclear Information System (INIS)
Araujo, Bruna G.M.; Dantas, Carlos C.; Santos, Valdemir A. dos; Finkler, Christine L.L.; Oliveira, Eric F. de; Melo, Silvio B.; Santos, M. Graca dos
2009-01-01
This paper the results obtained in single beam gamma ray tomography are investigated according to direct problem formulation and the applied solution for the linear system of equations. By image reconstruction based algebraic computational algorithms are used. The sparse under and over-determined linear system of equations was analyzed. Build in functions of Matlab software were applied and optimal solutions were investigate. Experimentally a section of the tube is scanned from various positions and at different angles. The solution, to find the vector of coefficients μ, from the vector of measured p values through the W matrix inversion, constitutes an inverse problem. A industrial tomography process requires a numerical solution of the system of equations. The definition of inverse problem according to Hadmard's is considered and as well the requirement of a well posed problem to find stable solutions. The formulation of the basis function and the computational algorithm to structure the weight matrix W were analyzed. For W full rank matrix the obtained solution is unique as expected. Total Least Squares was implemented which theory and computation algorithm gives adequate treatment for the problems due to non-unique solutions of the system of equations. Stability of the solution was investigating by means of a regularization technique and the comparison shows that it improves the results. An optimal solution as a function of the image quality, computation time and minimum residuals were quantified. The corresponding reconstructed images are shown in 3D graphics in order to compare with the solution. (author)
Oscillation of solutions of some higher order linear differential equations
Directory of Open Access Journals (Sweden)
Hong-Yan Xu
2009-11-01
Full Text Available In this paper, we deal with the order of growth and the hyper order of solutions of higher order linear differential equations $$f^{(k}+B_{k-1}f^{(k-1}+\\cdots+B_1f'+B_0f=F$$ where $B_j(z (j=0,1,\\ldots,k-1$ and $F$ are entire functions or polynomials. Some results are obtained which improve and extend previous results given by Z.-X. Chen, J. Wang, T.-B. Cao and C.-H. Li.
Description of All Solutions of a Linear Complementarity Problem
Czech Academy of Sciences Publication Activity Database
Rohn, Jiří
2009-01-01
Roč. 18, - (2009), s. 246-252 E-ISSN 1081-3810 R&D Projects: GA ČR GA201/09/1957; GA ČR GC201/08/J020 Institutional research plan: CEZ:AV0Z10300504 Keywords : linear complementarity problem * Moore-Penrose inverse * verified solution * absolute value equation Subject RIV: BA - General Mathematics Impact factor: 0.892, year: 2009 http://www.math.technion.ac.il/iic/ ela / ela -articles/articles/vol18_pp246-252.pdf
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
Directory of Open Access Journals (Sweden)
M. De la Sen
2009-01-01
Full Text Available This paper investigates the relations between the particular eigensolutions of a limiting functional differential equation of any order, which is the nominal (unperturbed linear autonomous differential equations, and the associate ones of the corresponding perturbed functional differential equation. Both differential equations involve point and distributed delayed dynamics including Volterra class dynamics. The proofs are based on a Perron-type theorem for functional equations so that the comparison is governed by the real part of a dominant zero of the characteristic equation of the nominal differential equation. The obtained results are also applied to investigate the global stability of the perturbed equation based on that of its corresponding limiting equation.
Linear quadratic optimization for positive LTI system
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.
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....
Application of Nearly Linear Solvers to Electric Power System Computation
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.
International Nuclear Information System (INIS)
Edery, D.
1983-11-01
The reduced system of the non linear resistive MHD equations is used in the 2-D one helicity approximation in the numerical computations of stationary tearing modes. The critical magnetic Raynolds number S (S=tausub(r)/tausub(H) where tausub(R) and tausub(H) are respectively the characteristic resistive and hydro magnetic times) and the corresponding linear solution are computed as a starting approximation for the full non linear equations. These equations are then treated numerically by an iterative procedure which is shown to be rapidly convergent. A numerical application is given in the last part of this paper
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
Perfect commuting-operator strategies for linear system games
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.
Linear quadratic Gaussian balancing for discrete-time infinite-dimensional linear systems
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
"Real-Time Optical Laboratory Linear Algebra Solution Of Partial Differential Equations"
Casasent, David; Jackson, James
1986-03-01
A Space Integrating (SI) Optical Linear Algebra Processor (OLAP) employing space and frequency-multiplexing, new partitioning and data flow, and achieving high accuracy performance with a non base-2 number system is described. Laboratory data on the performance of this system and the solution of parabolic Partial Differential Equations (PDEs) is provided. A multi-processor OLAP system is also described for the first time. It use in the solution of multiple banded matrices that frequently arise is then discussed. The utility and flexibility of this processor compared to digital systolic architectures should be apparent.
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)
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
Thermodynamics of (1-alkanol + linear monoether) systems
International Nuclear Information System (INIS)
Gonzalez, Juan Antonio; Mozo, Ismael; Garcia de la Fuente, Isaias; Cobos, Jose Carlos; Riesco, Nicolas
2008-01-01
Densities, ρ, and speeds of sound, u, of systems formed by 1-heptanol, or 1-octanol, or 1-decanol and dibutylether have been measured at a temperature of (293.15, 298.15, and 303.15) K and atmospheric pressure using a vibrating tube densimeter and sound analyser Anton Paar model DSA-5000. The ρ and u values were used to calculate excess molar volumes, V E , and deviations from the ideal behaviour of the thermal expansion coefficient, Δα p and of the isentropic compressibilities, Δκ S . The available database on molar excess enthalpies, H E , and V E for (1-alkanol + linear monoether) systems was used to investigate interactional and structural effects in such mixtures. The enthalpy of the OH...O bonds is lower for methanol solutions, and for the remainder systems, it is practically independent of the mixture compounds. The V E variation with the chain length of the 1-alkanol points out the existence of structural effects for systems including longer 1-alkanols. The ERAS model is applied to the studied mixtures. ERAS represents quite accurately H E and V E data using parameters which consistently depend on the molecular structure
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...
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 coefﬁcient matrix. The symmetric coefﬁcient 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.
Mallet, D. G.; McCue, S. W.
2009-01-01
The solution of linear ordinary differential equations (ODEs) is commonly taught in first-year undergraduate mathematics classrooms, but the understanding of the concept of a solution is not always grasped by students until much later. Recognizing what it is to be a solution of a linear ODE and how to postulate such solutions, without resorting to…
Observability of linear systems with saturated outputs
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.
General solutions of second-order linear difference equations of Euler type
Directory of Open Access Journals (Sweden)
Akane Hongyo
2017-01-01
Full Text Available The purpose of this paper is to give general solutions of linear difference equations which are related to the Euler-Cauchy differential equation \\(y^{\\prime\\prime}+(\\lambda/t^2y=0\\ or more general linear differential equations. We also show that the asymptotic behavior of solutions of the linear difference equations are similar to solutions of the linear differential equations.
Linear regression crash prediction models : issues and proposed solutions.
2010-05-01
The paper develops a linear regression model approach that can be applied to : crash data to predict vehicle crashes. The proposed approach involves novice data aggregation : to satisfy linear regression assumptions; namely error structure normality ...
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.
CEFR information management system solution
International Nuclear Information System (INIS)
Lu Fei; Zhao Jia'ning
2011-01-01
Based on finished information resources planning scheme for China sodium cooled experimental fast breeder reactor and the advanced information resources management solution concepts were applied, we got the building solution of CEFR information management systems. At the same time, the technical solutions of systems structures, logic structures, physical structures, development platforms and operation platforms for information resources management system in fast breeder reactors were developed, which provided programmatic introductions for development works in future. (authors)
Isolators Including Main Spring Linear Guide Systems
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.
Maximization of energy in the output of a linear system
International Nuclear Information System (INIS)
Dudley, D.G.
1976-01-01
A time-limited signal which, when passed through a linear system, maximizes the total output energy is considered. Previous work has shown that the solution is given by the eigenfunction associated with the maximum eigenvalue in a Hilbert-Schmidt integral equation. Analytical results are available for the case where the transfer function is a low-pass filter. This work is extended by obtaining a numerical solution to the integral equation which allows results for reasonably general transfer functions
Linear systems a measurement based approach
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.
Solution of second order linear fuzzy difference equation by Lagrange's multiplier method
Directory of Open Access Journals (Sweden)
Sankar Prasad Mondal
2016-06-01
Full Text Available In this paper we execute the solution procedure for second order linear fuzzy difference equation by Lagrange's multiplier method. In crisp sense the difference equation are easy to solve, but when we take in fuzzy sense it forms a system of difference equation which is not so easy to solve. By the help of Lagrange's multiplier we can solved it easily. The results are illustrated by two different numerical examples and followed by two applications.
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
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.
Localized solutions of non-linear Klein--Gordon equations
International Nuclear Information System (INIS)
Werle, J.
1977-05-01
Nondissipative, stationary solutions for a class of nonlinear Klein-Gordon equations for a scalar field were found explicitly. Since the field is different from zero only inside a sphere of definite radius, the solutions are called quantum droplets
Fundamental solution of the problem of linear programming and method of its determination
Petrunin, S. V.
1978-01-01
The idea of a fundamental solution to a problem in linear programming is introduced. A method of determining the fundamental solution and of applying this method to the solution of a problem in linear programming is proposed. Numerical examples are cited.
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.
Signals and transforms in linear systems analysis
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...
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...
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.)
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.
Energy Technology Data Exchange (ETDEWEB)
Lorber, A.A.; Carey, G.F.; Bova, S.W.; Harle, C.H. [Univ. of Texas, Austin, TX (United States)
1996-12-31
The connection between the solution of linear systems of equations by iterative methods and explicit time stepping techniques is used to accelerate to steady state the solution of ODE systems arising from discretized PDEs which may involve either physical or artificial transient terms. Specifically, a class of Runge-Kutta (RK) time integration schemes with extended stability domains has been used to develop recursion formulas which lead to accelerated iterative performance. The coefficients for the RK schemes are chosen based on the theory of Chebyshev iteration polynomials in conjunction with a local linear stability analysis. We refer to these schemes as Chebyshev Parameterized Runge Kutta (CPRK) methods. CPRK methods of one to four stages are derived as functions of the parameters which describe an ellipse {Epsilon} which the stability domain of the methods is known to contain. Of particular interest are two-stage, first-order CPRK and four-stage, first-order methods. It is found that the former method can be identified with any two-stage RK method through the correct choice of parameters. The latter method is found to have a wide range of stability domains, with a maximum extension of 32 along the real axis. Recursion performance results are presented below for a model linear convection-diffusion problem as well as non-linear fluid flow problems discretized by both finite-difference and finite-element methods.
International Nuclear Information System (INIS)
Shimizu, Yoshiaki
1988-01-01
Due to the simplicity and effectiveness, linear program has been popular in the actual optimization in various fields. In the previous study, the uncertainty involved in the model at the different stage of optimization was dealt with by post-optimizing analysis. But it often becomes insufficient to make a decision how to deal with an uncertain system especially suffering large parameter deviation. Recently in the field of processing systems, it is desired to obtain a flexible solution which can present the counterplan to a deviating system from a practical viewpoint. The scope of this preliminary note presents how to apply a methodology development to obtain the flexible solution of a linear program. For this purpose, a simple example associated with nuclear reactor decommissioning is shown. The problem to maximize a system performance given as an objective function under the constraint of the static behavior of the system is considered, and the flexible solution is determined. In Japan, the decommissioning of commercial nuclear power plants will being in near future, and the study using the retired research reactor JPDR is in progress. The planning of decontamination and the reuse of wastes is taken as the example. (Kako, I.)
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)
Unbounded solutions of quasi-linear difference equations
Czech Academy of Sciences Publication Activity Database
Cecchi, M.; Došlá, Zuzana; Marini, M.
2003-01-01
Roč. 45, 10-11 (2003), s. 1113-1123 ISSN 0898-1221 Institutional research plan: CEZ:AV0Z1019905 Keywords : nonlinear difference equation * possitive increasing solution * strongly increasing solution Subject RIV: BA - General Mathematics Impact factor: 0.498, year: 2003
Non linear photons: a non singular cosmological solution
International Nuclear Information System (INIS)
Alves, G.A.
1986-01-01
The validity of equivalence principle as principle of minimum coupling between field interactions, is discussed. The non minimum coupling between vector field and gravitational field, and some consequences of this coupling are analysed. Starting from spherical symmetry metric, the coupled field equations, obtaining exact solutions and interpreting these solutions, are solved. (M.C.K.) [pt
Unsteady Solution of Non-Linear Differential Equations Using Walsh Function Series
Gnoffo, Peter A.
2015-01-01
Walsh functions form an orthonormal basis set consisting of square waves. The discontinuous nature of square waves make the system well suited for representing functions with discontinuities. The product of any two Walsh functions is another Walsh function - a feature that can radically change an algorithm for solving non-linear partial differential equations (PDEs). The solution algorithm of non-linear differential equations using Walsh function series is unique in that integrals and derivatives may be computed using simple matrix multiplication of series representations of functions. Solutions to PDEs are derived as functions of wave component amplitude. Three sample problems are presented to illustrate the Walsh function series approach to solving unsteady PDEs. These include an advection equation, a Burgers equation, and a Riemann problem. The sample problems demonstrate the use of the Walsh function solution algorithms, exploiting Fast Walsh Transforms in multi-dimensions (O(Nlog(N))). Details of a Fast Walsh Reciprocal, defined here for the first time, enable inversion of aWalsh Symmetric Matrix in O(Nlog(N)) operations. Walsh functions have been derived using a fractal recursion algorithm and these fractal patterns are observed in the progression of pairs of wave number amplitudes in the solutions. These patterns are most easily observed in a remapping defined as a fractal fingerprint (FFP). A prolongation of existing solutions to the next highest order exploits these patterns. The algorithms presented here are considered a work in progress that provide new alternatives and new insights into the solution of non-linear PDEs.
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
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...
Linear dynamic coupling in geared rotor systems
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.
International Nuclear Information System (INIS)
Alvarez-Estrada, R.F.
1979-01-01
A comprehensive review of the inverse scattering solution of certain non-linear evolution equations of physical interest in one space dimension is presented. We explain in some detail the interrelated techniques which allow to linearize exactly the following equations: (1) the Korteweg and de Vries equation; (2) the non-linear Schrodinger equation; (3) the modified Korteweg and de Vries equation; (4) the Sine-Gordon equation. We concentrate in discussing the pairs of linear operators which accomplish such an exact linearization and the solution of the associated initial value problem. The application of the method to other non-linear evolution equations is reviewed very briefly
Some problems on non-linear semigroups and the blow-up of integral solutions
International Nuclear Information System (INIS)
Pavel, N.H.
1983-07-01
After some introductory remarks, this highly mathematical document considers a unifying approach in the theory of non-linear semigroups. Then a brief survey is given on blow-up of mild solutions from the semilinear case. Finally, the global behavior of solutions to non-linear evolution equations is addressed; it is found that classical results on the behavior of the maximal solution u as t up-arrow tsub(max) hold also for integral solutions
On output regulation for linear systems
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
Linear response theory for quantum open systems
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.
When to call a linear system nonnegative
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
Tikhonov theorem for linear hyperbolic systems
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...
Comparison of different methods for the solution of sets of linear equations
International Nuclear Information System (INIS)
Bilfinger, T.; Schmidt, F.
1978-06-01
The application of the conjugate-gradient methods as novel general iterative methods for the solution of sets of linear equations with symmetrical systems matrices led to this paper, where a comparison of these methods with the conventional differently accelerated Gauss-Seidel iteration was carried out. In additon, the direct Cholesky method was also included in the comparison. The studies referred mainly to memory requirement, computing time, speed of convergence, and accuracy of different conditions of the systems matrices, by which also the sensibility of the methods with respect to the influence of truncation errors may be recognized. (orig.) 891 RW [de
Conduction cooling systems for linear accelerator cavities
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.
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
Li, Yanning
2013-10-01
This article presents a new robust control framework for transportation problems in which the state is modeled by a first order scalar conservation law. Using an equivalent formulation based on a Hamilton-Jacobi equation, we pose the problem of controlling the state of the system on a network link, using boundary flow control, as a Linear Program. Unlike many previously investigated transportation control schemes, this method yields a globally optimal solution and is capable of handling shocks (i.e. discontinuities in the state of the system). We also demonstrate that the same framework can handle robust control problems, in which the uncontrollable components of the initial and boundary conditions are encoded in intervals on the right hand side of inequalities in the linear program. The lower bound of the interval which defines the smallest feasible solution set is used to solve the robust LP (or MILP if the objective function depends on boolean variables). Since this framework leverages the intrinsic properties of the Hamilton-Jacobi equation used to model the state of the system, it is extremely fast. Several examples are given to demonstrate the performance of the robust control solution and the trade-off between the robustness and the optimality. © 2013 IEEE.
Li, Yanning; Canepa, Edward S.; Claudel, Christian G.
2013-01-01
This article presents a new robust control framework for transportation problems in which the state is modeled by a first order scalar conservation law. Using an equivalent formulation based on a Hamilton-Jacobi equation, we pose the problem of controlling the state of the system on a network link, using boundary flow control, as a Linear Program. Unlike many previously investigated transportation control schemes, this method yields a globally optimal solution and is capable of handling shocks (i.e. discontinuities in the state of the system). We also demonstrate that the same framework can handle robust control problems, in which the uncontrollable components of the initial and boundary conditions are encoded in intervals on the right hand side of inequalities in the linear program. The lower bound of the interval which defines the smallest feasible solution set is used to solve the robust LP (or MILP if the objective function depends on boolean variables). Since this framework leverages the intrinsic properties of the Hamilton-Jacobi equation used to model the state of the system, it is extremely fast. Several examples are given to demonstrate the performance of the robust control solution and the trade-off between the robustness and the optimality. © 2013 IEEE.
General guidelines solution for linear programming with fuzzy coefficients
Directory of Open Access Journals (Sweden)
Sergio Gerardo de los Cobos Silva
2013-08-01
Full Text Available This work introduce to the Possibilistic Programming and the Fuzzy Programming as paradigms that allow to resolve problems of linear programming when the coefficients of the model or the restrictions on the same are presented as fuzzy numbers, rather than exact numbers (crisp. This work presents some examples based on [1].
Chaos as an intermittently forced linear system.
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.
POSITIVE SOLUTIONS TO SEMI-LINEAR SECOND-ORDER ORDINARY DIFFERENTIAL EQUATIONS IN BANACH SPACE
Institute of Scientific and Technical Information of China (English)
无
2008-01-01
In this paper,we study the existence of positive periodic solution to some second- order semi-linear differential equation in Banach space.By the fixed point index theory, we prove that the semi-linear differential equation has two positive periodic solutions.
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
Dual-range linearized transimpedance amplifier system
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).
On modulated complex non-linear dynamical systems
International Nuclear Information System (INIS)
Mahmoud, G.M.; Mohamed, A.A.; Rauh, A.
1999-01-01
This paper is concerned with the development of an approximate analytical method to investigate periodic solutions and their stability in the case of modulated non-linear dynamical systems whose equation of motion is describe. Such differential equations appear, for example, in problems of colliding particle beams in high-energy accelerators or one-mass systems with two or more degrees of freedom, e.g. rotors. The significance of periodic solutions lies on the fact that all non-periodic responses, if convergent, would approach to periodic solutions at the steady-state conditions. The example shows a good agreement between numerical and analytical results for small values of ε. The effect of the periodic modulation on the stability of the 2π-periodic solutions is discussed
No-signaling quantum key distribution: solution by linear programming
Hwang, Won-Young; Bae, Joonwoo; Killoran, Nathan
2015-02-01
We outline a straightforward approach for obtaining a secret key rate using only no-signaling constraints and linear programming. Assuming an individual attack, we consider all possible joint probabilities. Initially, we study only the case where Eve has binary outcomes, and we impose constraints due to the no-signaling principle and given measurement outcomes. Within the remaining space of joint probabilities, by using linear programming, we get bound on the probability of Eve correctly guessing Bob's bit. We then make use of an inequality that relates this guessing probability to the mutual information between Bob and a more general Eve, who is not binary-restricted. Putting our computed bound together with the Csiszár-Körner formula, we obtain a positive key generation rate. The optimal value of this rate agrees with known results, but was calculated in a more straightforward way, offering the potential of generalization to different scenarios.
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
Disturbance Decoupling of Switched Linear Systems
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
Uzawa method for fuzzy linear system
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.
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
Standard diffusive systems are well-posed linear systems
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.
Linear facility location in three dimensions - Models and solution methods
DEFF Research Database (Denmark)
Brimberg, Jack; Juel, Henrik; Schöbel, Anita
2002-01-01
We consider the problem of locating a line or a line segment in three-dimensional space, such that the sum of distances from the facility represented by the line (segment) to a given set of points is minimized. An example is planning the drilling of a mine shaft, with access to ore deposits through...... horizontal tunnels connecting the deposits and the shaft. Various models of the problem are developed and analyzed, and efficient solution methods are given....
Tunç, Cemil; Tunç, Osman
2016-01-01
In this paper, certain system of linear homogeneous differential equations of second-order is considered. By using integral inequalities, some new criteria for bounded and [Formula: see text]-solutions, upper bounds for values of improper integrals of the solutions and their derivatives are established to the considered system. The obtained results in this paper are considered as extension to the results obtained by Kroopnick (2014) [1]. An example is given to illustrate the obtained results.
Parameter identifiability of linear dynamical systems
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.
Role of statistical linearization in the solution of nonlinear stochastic equations
International Nuclear Information System (INIS)
Budgor, A.B.
1977-01-01
The solution of a generalized Langevin equation is referred to as a stochastic process. If the external forcing function is Gaussian white noise, the forward Kolmogarov equation yields the transition probability density function. Nonlinear problems must be handled by approximation procedures e.g., perturbation theories, eigenfunction expansions, and nonlinear optimization procedures. After some comments on the first two of these, attention is directed to the third, and the method of statistical linearization is used to demonstrate a relation to the former two. Nonlinear stochastic systems exhibiting sustained or forced oscillations and the centered nonlinear Schroedinger equation in the presence of Gaussian white noise excitation are considered as examples. 5 figures, 2 tables
Identification of general linear mechanical systems
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.
Tisdell, Christopher C.
2017-01-01
For over 50 years, the learning of teaching of "a priori" bounds on solutions to linear differential equations has involved a Euclidean approach to measuring the size of a solution. While the Euclidean approach to "a priori" bounds on solutions is somewhat manageable in the learning and teaching of the proofs involving…
A Fast Condensing Method for Solution of Linear-Quadratic Control Problems
DEFF Research Database (Denmark)
Frison, Gianluca; Jørgensen, John Bagterp
2013-01-01
consider a condensing (or state elimination) method to solve an extended version of the LQ control problem, and we show how to exploit the structure of this problem to both factorize the dense Hessian matrix and solve the system. Furthermore, we present two efficient implementations. The first......In both Active-Set (AS) and Interior-Point (IP) algorithms for Model Predictive Control (MPC), sub-problems in the form of linear-quadratic (LQ) control problems need to be solved at each iteration. The solution of these sub-problems is usually the main computational effort. In this paper we...... implementation is formally identical to the Riccati recursion based solver and has a computational complexity that is linear in the control horizon length and cubic in the number of states. The second implementation has a computational complexity that is quadratic in the control horizon length as well...
Large-time asymptotic behaviour of solutions of non-linear Sobolev-type equations
International Nuclear Information System (INIS)
Kaikina, Elena I; Naumkin, Pavel I; Shishmarev, Il'ya A
2009-01-01
The large-time asymptotic behaviour of solutions of the Cauchy problem is investigated for a non-linear Sobolev-type equation with dissipation. For small initial data the approach taken is based on a detailed analysis of the Green's function of the linear problem and the use of the contraction mapping method. The case of large initial data is also closely considered. In the supercritical case the asymptotic formulae are quasi-linear. The asymptotic behaviour of solutions of a non-linear Sobolev-type equation with a critical non-linearity of the non-convective kind differs by a logarithmic correction term from the behaviour of solutions of the corresponding linear equation. For a critical convective non-linearity, as well as for a subcritical non-convective non-linearity it is proved that the leading term of the asymptotic expression for large times is a self-similar solution. For Sobolev equations with convective non-linearity the asymptotic behaviour of solutions in the subcritical case is the product of a rarefaction wave and a shock wave. Bibliography: 84 titles.
Naganthran, Kohilavani; Nazar, Roslinda; Pop, Ioan
2018-05-01
This study investigated the influence of the non-linearly stretching/shrinking sheet on the boundary layer flow and heat transfer. A proper similarity transformation simplified the system of partial differential equations into a system of ordinary differential equations. This system of similarity equations is then solved numerically by using the bvp4c function in the MATLAB software. The generated numerical results presented graphically and discussed in the relevance of the governing parameters. Dual solutions found as the sheet stretched and shrunk in the horizontal direction. Stability analysis showed that the first solution is physically realizable whereas the second solution is not practicable.
Analytical Solutions to Non-linear Mechanical Oscillation Problems
DEFF Research Database (Denmark)
Kaliji, H. D.; Ghadimi, M.; Barari, Amin
2011-01-01
In this paper, the Max-Min Method is utilized for solving the nonlinear oscillation problems. The proposed approach is applied to three systems with complex nonlinear terms in their motion equations. By means of this method, the dynamic behavior of oscillation systems can be easily approximated u...
Lectures on algebraic system theory: Linear systems over rings
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.
Fundamental Matrix for a Class of Point Delay Linear Systems
International Nuclear Information System (INIS)
Sen, M. de la; Alastruey, C. F.
1998-01-01
It is difficult to establish explicit analytic forms for fundamental matrices of delayed linear systems. In this paper, an explicit form of exponential type is given for such a matrix in the case of punctual delays. The existence of real and complex fundamental matrices, for the case of real parameterizations of the differential system, is studied and discussed. Some additional commutativity properties involving the matrices parameters and the fundamental matrices as well as explicit expressions for the solution of the delayed differential system are also given. (Author)
Directory of Open Access Journals (Sweden)
Mervan Pašić
2016-10-01
Full Text Available We study non-monotone positive solutions of the second-order linear differential equations: $(p(tx'' + q(t x = e(t$, with positive $p(t$ and $q(t$. For the first time, some criteria as well as the existence and nonexistence of non-monotone positive solutions are proved in the framework of some properties of solutions $\\theta (t$ of the corresponding integrable linear equation: $(p(t\\theta''=e(t$. The main results are illustrated by many examples dealing with equations which allow exact non-monotone positive solutions not necessarily periodic. Finally, we pose some open questions.
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 ...
DEFF Research Database (Denmark)
Mejlbro, Leif
1997-01-01
An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians.......An alternative formula for the solution of linear differential equations of order n is suggested. When applicable, the suggested method requires fewer and simpler computations than the well-known method using Wronskians....
An injection system for a linear accelerator
International Nuclear Information System (INIS)
Santos, A.C.R.
1978-03-01
An injection system for the Linear Accelerator is developed using the parameters of machines at the Centro Brasileiro de Pesquisas Fisicas and the Instituto Militar de Engenharia. The proposed system consists basically of a prebuncher and a chopper. The pre-buncher is used to improve the energy resolution and also to increase the accelerator target current. The chopper is used to remove from the beam the electrons that have no possibility of attaining the desired energy and that are usually lost in the walls and the cavity tube, thus producing undesirable background. Theoretical development of the chopper is performed in order to obtain its dimensions for future construction. The complete design the pre-buncher and its feed supply system and the experimental verication of its performance are also presented. It is intended to give the necessary information for the design and construction of the complete injection system proposed. (Author) [pt
Operator approach to linear control systems
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.
Energy Technology Data Exchange (ETDEWEB)
Garcia Velarde, M
1977-07-01
Thermo convective instabilities in horizontal fluid layers are discussed with emphasis on the Rayleigh-Bernard model problem. Steady solutions and time-dependent phenomena (relaxation oscillations and transition to turbulence) are studied within the nonlinear Boussinesq-Oberbeck approximation. Homogeneous steady solutions, limit cycles, and inhomogeneous (ordered) spatial structures are also studied in simple reaction-diffusion systems. Lastly, the non-periodic attractor that appears at large Rayleigh numbers in the truncated Boussinesq-Oberbeck model of Lorenz, is constructed, and a discussion of turbulent behavior is given. (Author) 105 refs.
International Nuclear Information System (INIS)
Garcia Velarde, M.
1977-01-01
Thermoconvective instabilities in horizontal fluid layers are discussed with emphasis on the Rayleigh-Benard model problem. Steady solutions and time-dependent phenomena (relaxation oscillations and transition to turbulence) are studied within the nonlinear Boussinesq-Oberbeck approximation. Homogeneous steady solutions, limit cycles, and inhomogeneous (ordered) spatial structures are also studied in simple reaction-diffusion systems. Lastly, the non-periodic attractor that appears at large Rayleigh numbers in the truncated Boussinesq-Oberbeck model of Lorenz, is constructed, and a discussion of turbulent behavior is given. (author) [es
International Nuclear Information System (INIS)
Garcia Velarde, M.
1977-01-01
Thermo convective instabilities in horizontal fluid layers are discussed with emphasis on the Rayleigh-Bernard model problem. Steady solutions and time-dependent phenomena (relaxation oscillations and transition to turbulence) are studied within the nonlinear Boussinesq-Oberbeck approximation. Homogeneous steady solutions, limit cycles, and inhomogeneous (ordered) spatial structures are also studied in simple reaction-diffusion systems. Lastly, the non-periodic attractor that appears at large Rayleigh numbers in the truncated Boussinesq-Oberbeck model of Lorenz, is constructed, and a discussion of turbulent behavior is given. (Author) 105 refs
Nonautonomous linear system of the terrestrial carbon cycle
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
A parallel solver for huge dense linear systems
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
Directory of Open Access Journals (Sweden)
Hongchun Sun
2012-01-01
Full Text Available For the extended mixed linear complementarity problem (EML CP, we first present the characterization of the solution set for the EMLCP. Based on this, its global error bound is also established under milder conditions. The results obtained in this paper can be taken as an extension for the classical linear complementarity problems.
International Nuclear Information System (INIS)
Man, Yiu-Kwong
2010-01-01
In this communication, we present a method for computing the Liouvillian solution of second-order linear differential equations via algebraic invariant curves. The main idea is to integrate Kovacic's results on second-order linear differential equations with the Prelle-Singer method for computing first integrals of differential equations. Some examples on using this approach are provided. (fast track communication)
Dichotomy and almost automorphic solution of difference system
Directory of Open Access Journals (Sweden)
Samuel Castillo
2013-06-01
Full Text Available We study almost automorphic solutions of recurrence relations with values in a Banach space $V$ for quasilinear almost automorphic difference systems. Its linear part is a constant bounded linear operator $\\varLambda$ defined on $V$ satisfying an exponential dichotomy. We study the existence of almost automorphic solutions of the non-homogeneous linear difference equation and to quasilinear difference equation. Assuming global Lipschitz type conditions, we obtain Massera type results for these abstract systems. The case where the eigenvalues $\\lambda$ verify $\\left|\\lambda\\right|=1$ is also treated. An application to differential equations with piecewise constant argument is given.
Fast solution of elliptic partial differential equations using linear combinations of plane waves.
Pérez-Jordá, José M
2016-02-01
Given an arbitrary elliptic partial differential equation (PDE), a procedure for obtaining its solution is proposed based on the method of Ritz: the solution is written as a linear combination of plane waves and the coefficients are obtained by variational minimization. The PDE to be solved is cast as a system of linear equations Ax=b, where the matrix A is not sparse, which prevents the straightforward application of standard iterative methods in order to solve it. This sparseness problem can be circumvented by means of a recursive bisection approach based on the fast Fourier transform, which makes it possible to implement fast versions of some stationary iterative methods (such as Gauss-Seidel) consuming O(NlogN) memory and executing an iteration in O(Nlog(2)N) time, N being the number of plane waves used. In a similar way, fast versions of Krylov subspace methods and multigrid methods can also be implemented. These procedures are tested on Poisson's equation expressed in adaptive coordinates. It is found that the best results are obtained with the GMRES method using a multigrid preconditioner with Gauss-Seidel relaxation steps.
Considering system non-linearity in transmission pricing
International Nuclear Information System (INIS)
Oloomi-Buygi, M.; Salehizadeh, M. Reza
2008-01-01
In this paper a new approach for transmission pricing is presented. The contribution of a contract on power flow of a transmission line is used as extent-of-use criterion for transmission pricing. In order to determine the contribution of each contract on power flow of each transmission line, first the contribution of each contract on each voltage angle is determined, which is called voltage angle decomposition. To this end, DC power flow is used to compute a primary solution for voltage angle decomposition. To consider the impacts of system non-linearity on voltage angle decomposition, a method is presented to determine the share of different terms of sine argument in sine value. Then the primary solution is corrected in different iterations of decoupled Newton-Raphson power flow using the presented sharing method. The presented approach is applied to a 4-bus test system and IEEE 30-bus test system and the results are analyzed. (author)
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.
Non-linear dynamics of the passivity breakdown of iron in acidic solutions
Sazou, D
2003-01-01
Breakdown of the iron passivity in acid solutions accompanied by current oscillations was investigated by using electrochemical techniques, which reveal the non-linear dynamical response of the system in the current-potential (I-E) and current-time (I-t) planes. Current oscillations of the Fe-electrolyte electrochemical system were studied in the (a) absence and (b) presence of chlorides. In case (a) two oscillatory regions were distinguished; one at low potentials associated with the formation-dissolution of a ferrous salt and another at higher potentials associated with the formation-breakdown of the oxide film. Chaotic oscillations appear in the former region whereas periodic oscillations of a relaxation type appear in the latter region. In case (b), complex periodic and aperiodic oscillations are induced by small amounts of chlorides due to pitting corrosion. Pitting corrosion is a multistage localized process of a great technological importance. It consists of a local breakdown of the passive oxide film ...
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.
Conjugate gradient type methods for linear systems with complex symmetric coefficient matrices
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.
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, ...
Oscillation and nonoscillation results for solutions of half-linear equations with deviated argument
Czech Academy of Sciences Publication Activity Database
Drábek, P.; Kufner, Alois; Kuliev, K.
2017-01-01
Roč. 447, č. 1 (2017), s. 371-382 ISSN 0022-247X Institutional support: RVO:67985840 Keywords : half-linear equation * oscillatory solution * nonoscillatory solution Subject RIV: BA - General Mathematics OBOR OECD: Pure mathematics Impact factor: 1.064, year: 2016 http://www.sciencedirect.com/science/article/pii/S0022247X16306059
Directory of Open Access Journals (Sweden)
Yoshitsugu Takei
2015-01-01
Full Text Available Using two concrete examples, we discuss the multisummability of WKB solutions of singularly perturbed linear ordinary differential equations. Integral representations of solutions and a criterion for the multisummability based on the Cauchy-Heine transform play an important role in the proof.
Optimal linear precoding for indoor visible light communication system
Sifaou, Houssem
2017-07-31
Visible light communication (VLC) is an emerging technique that uses light-emitting diodes (LED) to combine communication and illumination. It is considered as a promising scheme for indoor wireless communication that can be deployed at reduced costs while offering high data rate performance. In this paper, we focus on the design of the downlink of a multi-user VLC system. Inherent to multi-user systems is the interference caused by the broadcast nature of the medium. Linear precoding based schemes are among the most popular solutions that have recently been proposed to mitigate inter-user interference. This paper focuses on the design of the optimal linear precoding scheme that solves the max-min signal-to-interference-plus-noise ratio (SINR) problem. The performance of the proposed precoding scheme is studied under different working conditions and compared with the classical zero-forcing precoding. Simulations have been provided to illustrate the high gain of the proposed scheme.
Nonautonomous linear Hamiltonian systems oscillation, spectral theory and control
Johnson, Russell; Novo, Sylvia; Núñez, Carmen; Fabbri, Roberta
2016-01-01
This monograph contains an in-depth analysis of the dynamics given by a linear Hamiltonian system of general dimension with nonautonomous bounded and uniformly continuous coefficients, without other initial assumptions on time-recurrence. Particular attention is given to the oscillation properties of the solutions as well as to a spectral theory appropriate for such systems. The book contains extensions of results which are well known when the coefficients are autonomous or periodic, as well as in the nonautonomous two-dimensional case. However, a substantial part of the theory presented here is new even in those much simpler situations. The authors make systematic use of basic facts concerning Lagrange planes and symplectic matrices, and apply some fundamental methods of topological dynamics and ergodic theory. Among the tools used in the analysis, which include Lyapunov exponents, Weyl matrices, exponential dichotomy, and weak disconjugacy, a fundamental role is played by the rotation number for linear Hami...
Solutions of the linearized Bach-Einstein equation in the static spherically symmetric case
International Nuclear Information System (INIS)
Schmidt, H.J.
1985-01-01
The Bach-Einstein equation linearized around Minkowski space-time is completely solved. The set of solutions depends on three parameters; a two-parameter subset of it becomes asymptotically flat. In that region the gravitational potential is of the type phi = -m/r + epsilon exp (-r/l). Because of the different asymptotic behaviour of both terms, it became necessary to linearize also around the Schwarzschild solution phi = -m/r. The linearized equation resulting in this case is discussed using qualitative methods. The result is that for m = 2l phi = -m/r + epsilon r -2 exp (-r/l) u, where u is some bounded function; m is arbitrary and epsilon again small. Further, the relation between the solution of the linearized and the full equation is discussed. (author)
On the summability of divergent power series solutions for certain first-order linear PDEs
Directory of Open Access Journals (Sweden)
Masaki Hibino
2015-01-01
Full Text Available This article is concerned with the study of the Borel summability of divergent power series solutions for certain singular first-order linear partial differential equations of nilpotent type. Our main purpose is to obtain conditions which coefficients of equations should satisfy in order to ensure the Borel summability of divergent solutions. We will see that there is a close affinity between the Borel summability of divergent solutions and global analytic continuation properties for coefficients of equations.
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.
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)
Linear concentration system; Sistema de concentracion lineal
Energy Technology Data Exchange (ETDEWEB)
Gonzalez Lugo, J.I; Leon Rovira, N; Aguayo Tellez, H [Instituto Tecnologico y de Estudios Superiores de Monterrey, Monterrey, Nuevo Leon (Mexico)]. E-mails: a00812662@itesm.mx; noel.leon@itesm.mx; haguayo@itesm.mx
2013-03-15
Solar linear concentration technologies to generate high temperatures are limited to the ranges of 200 to 500 degrees Celsius. While its performance has been tested through prototypes and pilot plants around the world, there are still areas of opportunity that can be exploited to obtain a linear concentration that achieves temperatures above this range in order to have a better use of the available solar energy. Because of this: It is possible to develop a linear concentration system that can track the sun with minimal movement of the absorber-receiver while maintaining temperatures above 850 degrees Celsius sufficient for industrial processes that require that temperature. The methodology consists of a series of stages (conceptual design, simulation, evaluation, development concept, results and validation) through which concepts are generated that allow design and evaluation of solar concentrator configurations with the help of simulation software. We have designed a linear parabolic concentrating system which comprises a set of mirrors segments with different focal lengths that works within the range of 600 degrees Celsius; however, it is advancing in the development of a double concentration to reach 850 degrees Celsius. [Spanish] Las tecnologias de concentracion lineal solar para generar altas temperaturas se ven limitadas a los rangos de 200 a 500 grados centigrados. Si bien su funcionamiento ha sido probado a traves de prototipos y plantas piloto alrededor del mundo, aun existen areas de oportunidad que pueden ser aprovechadas para obtener un sistema de concentracion lineal que permita alcanzar temperaturas mayores a este rango para asi tener un mejor aprovechamiento de la energia solar disponible. Debido a esto: Es posible desarrollar un sistema de concentracion lineal capaz de seguir la trayectoria del Sol con minimo movimiento del absorbedor-recibidor al mismo tiempo que mantiene temperaturas superiores a los 850 grados centigrados suficientes para
International Nuclear Information System (INIS)
Dubrovsky, V. G.; Topovsky, A. V.
2013-01-01
New exact solutions, nonstationary and stationary, of Veselov-Novikov (VN) equation in the forms of simple nonlinear and linear superpositions of arbitrary number N of exact special solutions u (n) , n= 1, …, N are constructed via Zakharov and Manakov ∂-dressing method. Simple nonlinear superpositions are represented up to a constant by the sums of solutions u (n) and calculated by ∂-dressing on nonzero energy level of the first auxiliary linear problem, i.e., 2D stationary Schrödinger equation. It is remarkable that in the zero energy limit simple nonlinear superpositions convert to linear ones in the form of the sums of special solutions u (n) . It is shown that the sums u=u (k 1 ) +...+u (k m ) , 1 ⩽k 1 2 m ⩽N of arbitrary subsets of these solutions are also exact solutions of VN equation. The presented exact solutions include as superpositions of special line solitons and also superpositions of plane wave type singular periodic solutions. By construction these exact solutions represent also new exact transparent potentials of 2D stationary Schrödinger equation and can serve as model potentials for electrons in planar structures of modern electronics.
Energy Technology Data Exchange (ETDEWEB)
Dubrovsky, V. G.; Topovsky, A. V. [Novosibirsk State Technical University, Karl Marx prosp. 20, Novosibirsk 630092 (Russian Federation)
2013-03-15
New exact solutions, nonstationary and stationary, of Veselov-Novikov (VN) equation in the forms of simple nonlinear and linear superpositions of arbitrary number N of exact special solutions u{sup (n)}, n= 1, Horizontal-Ellipsis , N are constructed via Zakharov and Manakov {partial_derivative}-dressing method. Simple nonlinear superpositions are represented up to a constant by the sums of solutions u{sup (n)} and calculated by {partial_derivative}-dressing on nonzero energy level of the first auxiliary linear problem, i.e., 2D stationary Schroedinger equation. It is remarkable that in the zero energy limit simple nonlinear superpositions convert to linear ones in the form of the sums of special solutions u{sup (n)}. It is shown that the sums u=u{sup (k{sub 1})}+...+u{sup (k{sub m})}, 1 Less-Than-Or-Slanted-Equal-To k{sub 1} < k{sub 2} < Horizontal-Ellipsis < k{sub m} Less-Than-Or-Slanted-Equal-To N of arbitrary subsets of these solutions are also exact solutions of VN equation. The presented exact solutions include as superpositions of special line solitons and also superpositions of plane wave type singular periodic solutions. By construction these exact solutions represent also new exact transparent potentials of 2D stationary Schroedinger equation and can serve as model potentials for electrons in planar structures of modern electronics.
Linear Actuator System for the NASA Docking System
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.
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 ...
International Nuclear Information System (INIS)
Paris, R.B.; Wood, A.D.
1984-11-01
The asymptotic expansions of solutions of a class of linear ordinary differential equations of arbitrary order n, containing a factor zsup(m) multiplying the lower order derivatives, are investigated for large values of z in the complex plane. Four classes of solutions are considered which exhibit the following behaviour as /z/ → infinity in certain sectors: (i) solutions whose behaviour is either exponentially large or algebraic (involving p ( < n) algebraic expansions), (ii) solutions which are exponentially small (iii) solutions with a single algebraic expansion and (iv) solutions which are even and odd functions of z whenever n+m is even. The asymptotic expansions of these solutions in a full neigbourhood of the point at infinity are obtained by means of the theory of the solutions in the case m=O developed in a previous paper
Control system analysis for the perturbed linear accelerator rf system
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
An implicit iterative scheme for solving large systems of linear equations
International Nuclear Information System (INIS)
Barry, J.M.; Pollard, J.P.
1986-12-01
An implicit iterative scheme for the solution of large systems of linear equations arising from neutron diffusion studies is presented. The method is applied to three-dimensional reactor studies and its performance is compared with alternative iterative approaches
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
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.
A linear complementarity method for the solution of vertical vehicle-track interaction
Zhang, Jian; Gao, Qiang; Wu, Feng; Zhong, Wan-Xie
2018-02-01
A new method is proposed for the solution of the vertical vehicle-track interaction including a separation between wheel and rail. The vehicle is modelled as a multi-body system using rigid bodies, and the track is treated as a three-layer beam model in which the rail is considered as an Euler-Bernoulli beam and both the sleepers and the ballast are represented by lumped masses. A linear complementarity formulation is directly established using a combination of the wheel-rail normal contact condition and the generalised-α method. This linear complementarity problem is solved using the Lemke algorithm, and the wheel-rail contact force can be obtained. Then the dynamic responses of the vehicle and the track are solved without iteration based on the generalised-α method. The same equations of motion for the vehicle and track are adopted at the different wheel-rail contact situations. This method can remove some restrictions, that is, time-dependent mass, damping and stiffness matrices of the coupled system, multiple equations of motion for the different contact situations and the effect of the contact stiffness. Numerical results demonstrate that the proposed method is effective for simulating the vehicle-track interaction including a separation between wheel and rail.
An extended GS method for dense linear systems
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.
Interpolation problem for the solutions of linear elasticity equations based on monogenic functions
Grigor'ev, Yuri; Gürlebeck, Klaus; Legatiuk, Dmitrii
2017-11-01
Interpolation is an important tool for many practical applications, and very often it is beneficial to interpolate not only with a simple basis system, but rather with solutions of a certain differential equation, e.g. elasticity equation. A typical example for such type of interpolation are collocation methods widely used in practice. It is known, that interpolation theory is fully developed in the framework of the classical complex analysis. However, in quaternionic analysis, which shows a lot of analogies to complex analysis, the situation is more complicated due to the non-commutative multiplication. Thus, a fundamental theorem of algebra is not available, and standard tools from linear algebra cannot be applied in the usual way. To overcome these problems, a special system of monogenic polynomials the so-called Pseudo Complex Polynomials, sharing some properties of complex powers, is used. In this paper, we present an approach to deal with the interpolation problem, where solutions of elasticity equations in three dimensions are used as an interpolation basis.
Cazzulani, Gabriele; Resta, Ferruccio; Ripamonti, Francesco
2012-04-01
During the last years, more and more mechanical applications saw the introduction of active control strategies. In particular, the need of improving the performances and/or the system health is very often associated to vibration suppression. This goal can be achieved considering both passive and active solutions. In this sense, many active control strategies have been developed, such as the Independent Modal Space Control (IMSC) or the resonant controllers (PPF, IRC, . . .). In all these cases, in order to tune and optimize the control strategy, the knowledge of the system dynamic behaviour is very important and it can be achieved both considering a numerical model of the system or through an experimental identification process. Anyway, dealing with non-linear or time-varying systems, a tool able to online identify the system parameters becomes a key-point for the control logic synthesis. The aim of the present work is the definition of a real-time technique, based on ARMAX models, that estimates the system parameters starting from the measurements of piezoelectric sensors. These parameters are returned to the control logic, that automatically adapts itself to the system dynamics. The problem is numerically investigated considering a carbon-fiber plate model forced through a piezoelectric patch.
Knotted solutions for linear and nonlinear theories: Electromagnetism and fluid dynamics
Directory of Open Access Journals (Sweden)
Daniel W.F. Alves
2017-10-01
Full Text Available We examine knotted solutions, the most simple of which is the “Hopfion”, from the point of view of relations between electromagnetism and ideal fluid dynamics. A map between fluid dynamics and electromagnetism works for initial conditions or for linear perturbations, allowing us to find new knotted fluid solutions. Knotted solutions are also found to be solutions of nonlinear generalizations of electromagnetism, and of quantum-corrected actions for electromagnetism coupled to other modes. For null configurations, electromagnetism can be described as a null pressureless fluid, for which we can find solutions from the knotted solutions of electromagnetism. We also map them to solutions of Euler's equations, obtained from a type of nonrelativistic reduction of the relativistic fluid equations.
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.
Classical solutions of non-linear sigma-models and their quantum fluctuations
International Nuclear Information System (INIS)
Din, A.M.
1980-05-01
I study the properties of O(N) and CPsup(n-1) non-linear sigma-models in the two dimensional Euclidean space. All classical solutions of the equations of motion can be characterized and in the CPsup(n-1) model they can be expressed in a simple and explicit way in terms of holomorphic vectors. The topological winding number and the action of the general CPsup(n-1) solution can be evaluated and the latter turns out always to be a integer multiple of 2π. I further discuss the stability of the solutions and the problem of one-loop calculations of quantum fluctuations around classical solutions
Lattice cluster theory of associating polymers. I. Solutions of linear telechelic polymer chains.
Dudowicz, Jacek; Freed, Karl F
2012-02-14
The lattice cluster theory (LCT) for the thermodynamics of a wide array of polymer systems has been developed by using an analogy to Mayer's virial expansions for non-ideal gases. However, the high-temperature expansion inherent to the LCT has heretofore precluded its application to systems exhibiting strong, specific "sticky" interactions. The present paper describes a reformulation of the LCT necessary to treat systems with both weak and strong, "sticky" interactions. This initial study concerns solutions of linear telechelic chains (with stickers at the chain ends) as the self-assembling system. The main idea behind this extension of the LCT lies in the extraction of terms associated with the strong interactions from the cluster expansion. The generalized LCT for sticky systems reduces to the quasi-chemical theory of hydrogen bonding of Panyioutou and Sanchez when correlation corrections are neglected in the LCT. A diagrammatic representation is employed to facilitate the evaluation of the corrections to the zeroth-order approximation from short range correlations. © 2012 American Institute of Physics
A solution approach for non-linear analysis of concrete members
International Nuclear Information System (INIS)
Hadi, N. M.; Das, S.
1999-01-01
Non-linear solution of reinforced concrete structural members, at and beyond its maximum strength poses complex numerical problems. This is due to the fact that concrete exhibits strain softening behaviour once it reaches its maximum strength. This paper introduces an improved non-linear solution capable to overcome the numerical problems efficiently. The paper also presents a new concept of modeling discrete cracks in concrete members by using gap elements. Gap elements are placed in between two adjacent concrete elements in tensile zone. The magnitude of elongation of gap elements, which represents the width of the crack in concrete, increases edith the increase of tensile stress in those elements. As a result, transfer of local from one concrete element to adjacent elements reduces. Results of non-linear finite element analysis of three concrete beams using this new solution strategy are compared with those obtained by other researchers, and a good agreement is achieved. (authors). 13 refs. 9 figs.,
International Nuclear Information System (INIS)
Grigoriu, Mircea; Samorodnitsky, Gennady
2004-01-01
Two methods are considered for assessing the asymptotic stability of the trivial solution of linear stochastic differential equations driven by Poisson white noise, interpreted as the formal derivative of a compound Poisson process. The first method attempts to extend a result for diffusion processes satisfying linear stochastic differential equations to the case of linear equations with Poisson white noise. The developments for the method are based on Ito's formula for semimartingales and Lyapunov exponents. The second method is based on a geometric ergodic theorem for Markov chains providing a criterion for the asymptotic stability of the solution of linear stochastic differential equations with Poisson white noise. Two examples are presented to illustrate the use and evaluate the potential of the two methods. The examples demonstrate limitations of the first method and the generality of the second method
System theory as applied differential geometry. [linear system
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.
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
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)
Improved harmonic balance approach to periodic solutions of non-linear jerk equations
International Nuclear Information System (INIS)
Wu, B.S.; Lim, C.W.; Sun, W.P.
2006-01-01
An analytical approximate approach for determining periodic solutions of non-linear jerk equations involving third-order time-derivative is presented. This approach incorporates salient features of both Newton's method and the method of harmonic balance. By appropriately imposing the method of harmonic balance to the linearized equation, the approach requires only one or two iterations to predict very accurate analytical approximate solutions for a large range of initial velocity amplitude. One typical example is used to verify and illustrate the usefulness and effectiveness of the proposed approach
Zeb, Salman; Yousaf, Muhammad
2017-01-01
In this article, we present a QR updating procedure as a solution approach for linear least squares problem with equality constraints. We reduce the constrained problem to unconstrained linear least squares and partition it into a small subproblem. The QR factorization of the subproblem is calculated and then we apply updating techniques to its upper triangular factor R to obtain its solution. We carry out the error analysis of the proposed algorithm to show that it is backward stable. We also illustrate the implementation and accuracy of the proposed algorithm by providing some numerical experiments with particular emphasis on dense problems.
Tisdell, Christopher C.
2017-11-01
For over 50 years, the learning of teaching of a priori bounds on solutions to linear differential equations has involved a Euclidean approach to measuring the size of a solution. While the Euclidean approach to a priori bounds on solutions is somewhat manageable in the learning and teaching of the proofs involving second-order, linear problems with constant co-efficients, we believe it is not pedagogically optimal. Moreover, the Euclidean method becomes pedagogically unwieldy in the proofs involving higher-order cases. The purpose of this work is to propose a simpler pedagogical approach to establish a priori bounds on solutions by considering a different way of measuring the size of a solution to linear problems, which we refer to as the Uber size. The Uber form enables a simplification of pedagogy from the literature and the ideas are accessible to learners who have an understanding of the Fundamental Theorem of Calculus and the exponential function, both usually seen in a first course in calculus. We believe that this work will be of mathematical and pedagogical interest to those who are learning and teaching in the area of differential equations or in any of the numerous disciplines where linear differential equations are used.
Linear and quadratic exponential modulation of the solutions of the paraxial wave equation
International Nuclear Information System (INIS)
Torre, A
2010-01-01
A review of well-known transformations, which allow us to pass from one solution of the paraxial wave equation (PWE) (in one transverse space variable) to another, is presented. Such transformations are framed within the unifying context of the Lie algebra formalism, being related indeed to symmetries of the PWE. Due to the closure property of the symmetry group of the PWE we are led to consider as not trivial only the linear and the quadratic exponential modulation (accordingly, accompanied by a suitable shift or scaling of the space variables) of the original solutions of the PWE, which are seen to be just conveyed by a linear and a quadratic exponential modulation of the relevant 'source' functions. We will see that recently introduced solutions of the 1D PWE in both rectangular and polar coordinates can be deduced from already known solutions through the resulting symmetry transformation related schemes
Solution to the Diffusion equation for multi groups in X Y geometry using Linear Perturbation theory
International Nuclear Information System (INIS)
Mugica R, C.A.
2004-01-01
Diverse methods exist to solve numerically the neutron diffusion equation for several energy groups in stationary state among those that highlight those of finite elements. In this work the numerical solution of this equation is presented using Raviart-Thomas nodal methods type finite element, the RT0 and RT1, in combination with iterative techniques that allow to obtain the approached solution in a quick form. Nevertheless the above mentioned, the precision of a method is intimately bound to the dimension of the approach space by cell, 5 for the case RT0 and 12 for the RT1, and/or to the mesh refinement, that makes the order of the problem of own value to solve to grow considerably. By this way if it wants to know an acceptable approach to the value of the effective multiplication factor of the system when this it has experimented a small perturbation it was appeal to the Linear perturbation theory with which is possible to determine it starting from the neutron flow and of the effective multiplication factor of the not perturbed case. Results are presented for a reference problem in which a perturbation is introduced in an assemble that simulates changes in the control bar. (Author)
Ladiges, Daniel R.; Sader, John E.
2018-05-01
Nanomechanical resonators and sensors, operated in ambient conditions, often generate low-Mach-number oscillating rarefied gas flows. Cercignani [C. Cercignani, J. Stat. Phys. 1, 297 (1969), 10.1007/BF01007482] proposed a variational principle for the linearized Boltzmann equation, which can be used to derive approximate analytical solutions of steady (time-independent) flows. Here we extend and generalize this principle to unsteady oscillatory rarefied flows and thus accommodate resonating nanomechanical devices. This includes a mathematical approach that facilitates its general use and allows for systematic improvements in accuracy. This formulation is demonstrated for two canonical flow problems: oscillatory Couette flow and Stokes' second problem. Approximate analytical formulas giving the bulk velocity and shear stress, valid for arbitrary oscillation frequency, are obtained for Couette flow. For Stokes' second problem, a simple system of ordinary differential equations is derived which may be solved to obtain the desired flow fields. Using this framework, a simple and accurate formula is provided for the shear stress at the oscillating boundary, again for arbitrary frequency, which may prove useful in application. These solutions are easily implemented on any symbolic or numerical package, such as Mathematica or matlab, facilitating the characterization of flows produced by nanomechanical devices and providing insight into the underlying flow physics.
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.
Analytical solution for a linearly graded-index-profile planar waveguide.
Touam, T; Yergeau, F
1993-01-20
An analytical solution is presented for the TE modes of a planar waveguide structure comprising a high-index guiding layer and a buried layer with a profile such that the square of the index varies linearly and matches the substrate and high-index guiding layer. The electric-field profiles and the dispersion relation are obtained and discussed, and a solution by the WKB method is compared.
Fuzzy solution of the linear programming problem with interval coefficients in the constraints
Dorota Kuchta
2005-01-01
A fuzzy concept of solving the linear programming problem with interval coefficients is proposed. For each optimism level of the decision maker (where the optimism concerns the certainty that no errors have been committed in the estimation of the interval coefficients and the belief that optimistic realisations of the interval coefficients will occur) another interval solution of the problem will be generated and the decision maker will be able to choose the final solution having a complete v...
Superconducting linear accelerator system for NSC
Indian Academy of Sciences (India)
This paper reports the construction of a superconducting linear accelerator as a booster to the 15 UD Pelletron accelerator at Nuclear Science Centre, New Delhi. The LINAC will use superconducting niobium quarter wave resonators as the accelerating element. Construction of the linear accelerator has progressed ...
Solutions of First-Order Volterra Type Linear Integrodifferential Equations by Collocation Method
Directory of Open Access Journals (Sweden)
Olumuyiwa A. Agbolade
2017-01-01
Full Text Available The numerical solutions of linear integrodifferential equations of Volterra type have been considered. Power series is used as the basis polynomial to approximate the solution of the problem. Furthermore, standard and Chebyshev-Gauss-Lobatto collocation points were, respectively, chosen to collocate the approximate solution. Numerical experiments are performed on some sample problems already solved by homotopy analysis method and finite difference methods. Comparison of the absolute error is obtained from the present method and those from aforementioned methods. It is also observed that the absolute errors obtained are very low establishing convergence and computational efficiency.
Fall with linear drag and Wien's displacement law: approximate solution and Lambert function
International Nuclear Information System (INIS)
Vial, Alexandre
2012-01-01
We present an approximate solution for the downward time of travel in the case of a mass falling with a linear drag force. We show how a quasi-analytical solution implying the Lambert function can be found. We also show that solving the previous problem is equivalent to the search for Wien's displacement law. These results can be of interest for undergraduate students, as they show that some transcendental equations found in physics may be solved without purely numerical methods. Moreover, as will be seen in the case of Wien's displacement law, solutions based on series expansion can be very accurate even with few terms. (paper)
Minimization of Linear Functionals Defined on| Solutions of Large-Scale Discrete Ill-Posed Problems
DEFF Research Database (Denmark)
Elden, Lars; Hansen, Per Christian; Rojas, Marielba
2003-01-01
The minimization of linear functionals de ned on the solutions of discrete ill-posed problems arises, e.g., in the computation of con dence intervals for these solutions. In 1990, Elden proposed an algorithm for this minimization problem based on a parametric-programming reformulation involving...... the solution of a sequence of trust-region problems, and using matrix factorizations. In this paper, we describe MLFIP, a large-scale version of this algorithm where a limited-memory trust-region solver is used on the subproblems. We illustrate the use of our algorithm in connection with an inverse heat...
International Nuclear Information System (INIS)
Datta, Dhurjati Prasad; Bose, Manoj Kumar
2004-01-01
We present a new one parameter family of second derivative discontinuous solutions to the simplest scale invariant linear ordinary differential equation. We also point out how the construction could be extended to generate families of higher derivative discontinuous solutions as well. The discontinuity can occur only for a subset of even order derivatives, viz., 2nd, 4th, 8th, 16th,.... The solutions are shown to break the discrete parity (reflection) symmetry of the underlying equation. These results are expected to gain significance in the contemporary search of a new dynamical principle for understanding complex phenomena in nature
Solutions of half-linear differential equations in the classes Gamma and Pi
Czech Academy of Sciences Publication Activity Database
Řehák, Pavel; Taddei, V.
2016-01-01
Roč. 29, 7-8 (2016), s. 683-714 ISSN 0893-4983 Institutional support: RVO:67985840 Keywords : half-linear differential equation * positive solution * asymptotic formula Subject RIV: BA - General Mathematics Impact factor: 0.565, year: 2016 http://projecteuclid.org/euclid.die/1462298681
Exact solutions of linearized Schwinger endash Dyson equation of fermion self-energy
International Nuclear Information System (INIS)
Zhou, B.
1997-01-01
The Schwinger endash Dyson equation of fermion self-energy in the linearization approximation is solved exactly in a theory with gauge and effective four-fermion interactions. Different expressions for the independent solutions, which, respectively, submit to irregular and regular ultraviolet boundary condition are derived and expounded. copyright 1997 American Institute of Physics
Bounded solutions of self-adjoint second order linear difference equations with periodic coeffients
Directory of Open Access Journals (Sweden)
Encinas A.M.
2018-02-01
Full Text Available In this work we obtain easy characterizations for the boundedness of the solutions of the discrete, self–adjoint, second order and linear unidimensional equations with periodic coefficients, including the analysis of the so-called discrete Mathieu equations as particular cases.
Gasyna, Zbigniew L.
2008-01-01
Computational experiment is proposed in which a linear algebra method is applied to the solution of the Schrodinger equation for a diatomic oscillator. Calculations of the vibration-rotation spectrum for the HCl molecule are presented and the results show excellent agreement with experimental data. (Contains 1 table and 1 figure.)
Symmetric linear systems - An application of algebraic systems theory
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.
Exact solution to the Coulomb wave using the linearized phase-amplitude method
Directory of Open Access Journals (Sweden)
Shuji Kiyokawa
2015-08-01
Full Text Available The author shows that the amplitude equation from the phase-amplitude method of calculating continuum wave functions can be linearized into a 3rd-order differential equation. Using this linearized equation, in the case of the Coulomb potential, the author also shows that the amplitude function has an analytically exact solution represented by means of an irregular confluent hypergeometric function. Furthermore, it is shown that the exact solution for the Coulomb potential reproduces the wave function for free space expressed by the spherical Bessel function. The amplitude equation for the large component of the Dirac spinor is also shown to be the linearized 3rd-order differential equation.
Heinkenschloss, Matthias
2005-01-01
We study a class of time-domain decomposition-based methods for the numerical solution of large-scale linear quadratic optimal control problems. Our methods are based on a multiple shooting reformulation of the linear quadratic optimal control problem as a discrete-time optimal control (DTOC) problem. The optimality conditions for this DTOC problem lead to a linear block tridiagonal system. The diagonal blocks are invertible and are related to the original linear quadratic optimal control problem restricted to smaller time-subintervals. This motivates the application of block Gauss-Seidel (GS)-type methods for the solution of the block tridiagonal systems. Numerical experiments show that the spectral radii of the block GS iteration matrices are larger than one for typical applications, but that the eigenvalues of the iteration matrices decay to zero fast. Hence, while the GS method is not expected to convergence for typical applications, it can be effective as a preconditioner for Krylov-subspace methods. This is confirmed by our numerical tests.A byproduct of this research is the insight that certain instantaneous control techniques can be viewed as the application of one step of the forward block GS method applied to the DTOC optimality system.
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.
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.
THE HIERARCHY OF COMPONENT FOR TASKS SOLUTION IN THE COURSE OF “LINEAR ALGEBRA”.
Directory of Open Access Journals (Sweden)
V. S. Kruglyk
2009-06-01
Full Text Available In the present article the case in point is an application of new informational technologies in the process of teaching Linear Algebra in Kherson State University. The component-oriented approach to the teaching Linear Algebra is examined as well as hierarchy of components which is used in the system and advantages of such approach in comparison with traditional approach.
Adding salt to a surfactant solution: Linear rheological response of the resulting morphologies
Energy Technology Data Exchange (ETDEWEB)
Gaudino, Danila; Pasquino, Rossana, E-mail: r.pasquino@unina.it; Grizzuti, Nino [DICMaPI, Università degli Studi di Napoli Federico II, P.le Tecchio 80, 80125 Napoli (Italy)
2015-11-15
The micellar system composed of Cetylpyridinium Chloride-Sodium Salicylate (CPyCl-NaSal) in brine aqueous solutions has been studied by systematically changing the salt concentration, in order to investigate the rheology of the arising morphologies. In particular, the zero-shear viscosity and the linear viscoelastic response have been measured as a function of the NaSal concentration (with [CPyCl] = 100 mM). The Newtonian viscosity shows a nonmonotonic dependence upon concentration, passing through a maximum at NaSal/CPyCl ≈ 0.6, and eventually dropping at higher salt concentrations. The progressive addition of salt determines first a transition from a Newtonian to a purely Maxwell-like behavior as the length of the micelles significantly increases. Beyond the peak viscosity, the viscoelastic data show two distinct features. On the one hand, the main relaxation time of the system strongly decreases, while the plateau modulus remains essentially constant. Calculations based on the rheological data show that, as the binding salt concentration increases, there is a decrease in micelles breaking rate and a decrease in their average length. On the other hand, in the same concentration region, a low-frequency elastic plateau is measured. Such a plateau is considered as the signature of a tenuous, but persistent branched network, whose existence is confirmed by cryo-transmission electron microscopy images.
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.
Practical application of equivalent linearization approaches to nonlinear piping systems
International Nuclear Information System (INIS)
Park, Y.J.; Hofmayer, C.H.
1995-01-01
The use of mechanical energy absorbers as an alternative to conventional hydraulic and mechanical snubbers for piping supports has attracted a wide interest among researchers and practitioners in the nuclear industry. The basic design concept of energy absorbers (EA) is to dissipate the vibration energy of piping systems through nonlinear hysteretic actions of EA exclamation point s under design seismic loads. Therefore, some type of nonlinear analysis needs to be performed in the seismic design of piping systems with EA supports. The equivalent linearization approach (ELA) can be a practical analysis tool for this purpose, particularly when the response approach (RSA) is also incorporated in the analysis formulations. In this paper, the following ELA/RSA methods are presented and compared to each other regarding their practice and numerical accuracy: Response approach using the square root of sum of squares (SRSS) approximation (denoted RS in this paper). Classical ELA based on modal combinations and linear random vibration theory (denoted CELA in this paper). Stochastic ELA based on direct solution of response covariance matrix (denoted SELA in this paper). New algorithms to convert response spectra to the equivalent power spectral density (PSD) functions are presented for both the above CELA and SELA methods. The numerical accuracy of the three EL are studied through a parametric error analysis. Finally, the practicality of the presented analysis is demonstrated in two application examples for piping systems with EA supports
Numerical solution of non-linear dual-phase-lag bioheat transfer equation within skin tissues.
Kumar, Dinesh; Kumar, P; Rai, K N
2017-11-01
This paper deals with numerical modeling and simulation of heat transfer in skin tissues using non-linear dual-phase-lag (DPL) bioheat transfer model under periodic heat flux boundary condition. The blood perfusion is assumed temperature-dependent which results in non-linear DPL bioheat transfer model in order to predict more accurate results. A numerical method of line which is based on finite difference and Runge-Kutta (4,5) schemes, is used to solve the present non-linear problem. Under specific case, the exact solution has been obtained and compared with the present numerical scheme, and we found that those are in good agreement. A comparison based on model selection criterion (AIC) has been made among non-linear DPL models when the variation of blood perfusion rate with temperature is of constant, linear and exponential type with the experimental data and it has been found that non-linear DPL model with exponential variation of blood perfusion rate is closest to the experimental data. In addition, it is found that due to absence of phase-lag phenomena in Pennes bioheat transfer model, it achieves steady state more quickly and always predict higher temperature than thermal and DPL non-linear models. The effect of coefficient of blood perfusion rate, dimensionless heating frequency and Kirchoff number on dimensionless temperature distribution has also been analyzed. The whole analysis is presented in dimensionless form. Copyright © 2017 Elsevier Inc. All rights reserved.
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
High density linear systems for fusion power
International Nuclear Information System (INIS)
Ellis, W.R.; Krakowski, R.A.
1975-01-01
The physics and technological limitations and uncertainties associated with the linear theta pinch are discussed in terms of a generalized energy balance, which has as its basis the ratio (Q/sub E/) of total electrical energy generated to net electrical energy consumed. Included in this total is the virtual energy of bred fissile fuel, if a hybrid blanket is used, as well as the actual of real energy deposited in the blanket by the fusion neutron. The advantages and disadvantages of the pulsed operation demanded by the linear theta pinch are also discussed
2013-01-01
This book consists of twenty seven chapters, which can be divided into three large categories: articles with the focus on the mathematical treatment of non-linear problems, including the methodologies, algorithms and properties of analytical and numerical solutions to particular non-linear problems; theoretical and computational studies dedicated to the physics and chemistry of non-linear micro-and nano-scale systems, including molecular clusters, nano-particles and nano-composites; and, papers focused on non-linear processes in medico-biological systems, including mathematical models of ferments, amino acids, blood fluids and polynucleic chains.
Analysis of Linear Hybrid Systems in CLP
DEFF Research Database (Denmark)
Banda, Gourinath; Gallagher, John Patrick
2009-01-01
In this paper we present a procedure for representing the semantics of linear hybrid automata (LHAs) as constraint logic programs (CLP); flexible and accurate analysis and verification of LHAs can then be performed using generic CLP analysis and transformation tools. LHAs provide an expressive...
Linear System of Equations, Matrix Inversion, and Linear Programming Using MS Excel
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…
Linear Scaling Solution of the Time-Dependent Self-Consistent-Field Equations
Directory of Open Access Journals (Sweden)
Matt Challacombe
2014-03-01
Full Text Available A new approach to solving the Time-Dependent Self-Consistent-Field equations is developed based on the double quotient formulation of Tsiper 2001 (J. Phys. B. Dual channel, quasi-independent non-linear optimization of these quotients is found to yield convergence rates approaching those of the best case (single channel Tamm-Dancoff approximation. This formulation is variational with respect to matrix truncation, admitting linear scaling solution of the matrix-eigenvalue problem, which is demonstrated for bulk excitons in the polyphenylene vinylene oligomer and the (4,3 carbon nanotube segment.
Modulation equations for spatially periodic systems: derivation and solutions
Schielen, R.; Doelman, A.
1996-01-01
We study a class of partial dierential equations in one spatial dimension, which can be seen as model equations for the analysis of pattern formation in physical systems dened on unbounded, weakly oscillating domains. We perform a linear and weakly nonlinear stability analysis for solutions that
Almost periodic solutions to systems of parabolic equations
Directory of Open Access Journals (Sweden)
Janpou Nee
1994-01-01
Full Text Available In this paper we show that the second-order differential solution is 2-almost periodic, provided it is 2-bounded, and the growth of the components of a non-linear function of a system of parabolic equation is bounded by any pair of con-secutive eigenvalues of the associated Dirichlet boundary value problems.
A SYSTEMIC VISION OF BIOLOGY: OVERCOMING LINEARITY
Directory of Open Access Journals (Sweden)
M. Mayer
2005-07-01
were used to build a hipermedia material. This technology permit overcomes a linear communication, improving the comprehension of the network perspective. The teachers speeches revealed their conceptual con- structions along the course, showed the development of the competences in identify interconnection points in the flow and chemical cycling of energy, compatible with a systemic view of life.
International Nuclear Information System (INIS)
LaChapelle, J.
2004-01-01
A path integral is presented that solves a general class of linear second order partial differential equations with Dirichlet/Neumann boundary conditions. Elementary kernels are constructed for both Dirichlet and Neumann boundary conditions. The general solution can be specialized to solve elliptic, parabolic, and hyperbolic partial differential equations with boundary conditions. This extends the well-known path integral solution of the Schroedinger/diffusion equation in unbounded space. The construction is based on a framework for functional integration introduced by Cartier/DeWitt-Morette
q-analogue of summability of formal solutions of some linear q-difference-differential equations
Directory of Open Access Journals (Sweden)
Hidetoshi Tahara
2015-01-01
Full Text Available Let \\(q\\gt 1\\. The paper considers a linear \\(q\\-difference-differential equation: it is a \\(q\\-difference equation in the time variable \\(t\\, and a partial differential equation in the space variable \\(z\\. Under suitable conditions and by using \\(q\\-Borel and \\(q\\-Laplace transforms (introduced by J.-P. Ramis and C. Zhang, the authors show that if it has a formal power series solution \\(\\hat{X}(t,z\\ one can construct an actual holomorphic solution which admits \\(\\hat{X}(t,z\\ as a \\(q\\-Gevrey asymptotic expansion of order \\(1\\.
ZHU, C. S.; ROBB, D. A.; EWINS, D. J.
2002-05-01
The multiple-solution response of rotors supported on squeeze film dampers is a typical non-linear phenomenon. The behaviour of the multiple-solution response in a flexible rotor supported on two identical squeeze film dampers with centralizing springs is studied by three methods: synchronous circular centred-orbit motion solution, numerical integration method and slow acceleration method using the assumption of a short bearing and cavitated oil film; the differences of computational results obtained by the three different methods are compared in this paper. It is shown that there are three basic forms for the multiple-solution response in the flexible rotor system supported on the squeeze film dampers, which are the resonant, isolated bifurcation and swallowtail bifurcation multiple solutions. In the multiple-solution speed regions, the rotor motion may be subsynchronous, super-subsynchronous, almost-periodic and even chaotic, besides synchronous circular centred, even if the gravity effect is not considered. The assumption of synchronous circular centred-orbit motion for the journal and rotor around the static deflection line can be used only in some special cases; the steady state numerical integration method is very useful, but time consuming. Using the slow acceleration method, not only can the multiple-solution speed regions be detected, but also the non-synchronous response regions.
Dynamics of unsymmetric piecewise-linear/non-linear systems using finite elements in time
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.
Solutions to the linearized Navier-Stokes equations for channel flow via the WKB approximation
Leonard, Anthony
2017-11-01
Progress on determining semi-analytical solutions to the linearized Navier-Stokes equations for incompressible channel flow, laminar and turbulent, is reported. Use of the WKB approximation yields, e.g., solutions to initial-value problem for the inviscid Orr-Sommerfeld equation in terms of the Bessel functions J+ 1 / 3 ,J- 1 / 3 ,J1 , and Y1 and their modified counterparts for any given wave speed c = ω /kx and k⊥ ,(k⊥2 =kx2 +kz2) . Of particular note to be discussed is a sequence i = 1 , 2 , . . . of homogeneous inviscid solutions with complex k⊥ i for each speed c, (0 < c <=Umax), in the downstream direction. These solutions for the velocity component normal to the wall v are localized in the plane parallel to the wall. In addition, for limited range of negative c, (- c * <= c <= 0) , we have found upstream-traveling homogeneous solutions with real k⊥(c) . In both cases the solutions for v serve as a source for corresponding solutions to the inviscid Squire equation for the vorticity component normal to the wall ωy.
Multiobjective Optimization of Linear Cooperative Spectrum Sensing: Pareto Solutions and Refinement.
Yuan, Wei; You, Xinge; Xu, Jing; Leung, Henry; Zhang, Tianhang; Chen, Chun Lung Philip
2016-01-01
In linear cooperative spectrum sensing, the weights of secondary users and detection threshold should be optimally chosen to minimize missed detection probability and to maximize secondary network throughput. Since these two objectives are not completely compatible, we study this problem from the viewpoint of multiple-objective optimization. We aim to obtain a set of evenly distributed Pareto solutions. To this end, here, we introduce the normal constraint (NC) method to transform the problem into a set of single-objective optimization (SOO) problems. Each SOO problem usually results in a Pareto solution. However, NC does not provide any solution method to these SOO problems, nor any indication on the optimal number of Pareto solutions. Furthermore, NC has no preference over all Pareto solutions, while a designer may be only interested in some of them. In this paper, we employ a stochastic global optimization algorithm to solve the SOO problems, and then propose a simple method to determine the optimal number of Pareto solutions under a computational complexity constraint. In addition, we extend NC to refine the Pareto solutions and select the ones of interest. Finally, we verify the effectiveness and efficiency of the proposed methods through computer simulations.
Directory of Open Access Journals (Sweden)
Qiong Liu
2012-01-01
Full Text Available We study the following fourth-order elliptic equations: Δ2+Δ=(,,∈Ω,=Δ=0,∈Ω, where Ω⊂ℝ is a bounded domain with smooth boundary Ω and (, is asymptotically linear with respect to at infinity. Using an equivalent version of Cerami's condition and the symmetric mountain pass lemma, we obtain the existence of multiple solutions for the equations.
On oscillations of solutions to second-order linear delay differential equations
Czech Academy of Sciences Publication Activity Database
Opluštil, Z.; Šremr, Jiří
2013-01-01
Roč. 20, č. 1 (2013), s. 65-94 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillatory solution Subject RIV: BA - General Mathematics Impact factor: 0.340, year: 2013 http://www.degruyter.com/view/j/gmj.2013.20.issue-1/gmj-2013-0001/gmj-2013-0001.xml?format=INT
Bounds and maximum principles for the solution of the linear transport equation
International Nuclear Information System (INIS)
Larsen, E.W.
1981-01-01
Pointwise bounds are derived for the solution of time-independent linear transport problems with surface sources in convex spatial domains. Under specified conditions, upper bounds are derived which, as a function of position, decrease with distance from the boundary. Also, sufficient conditions are obtained for the existence of maximum and minimum principles, and a counterexample is given which shows that such principles do not always exist
On oscillations of solutions to second-order linear delay differential equations
Czech Academy of Sciences Publication Activity Database
Opluštil, Z.; Šremr, Jiří
2013-01-01
Roč. 20, č. 1 (2013), s. 65-94 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : linear second-order delay differential equation * oscillatory solution Subject RIV: BA - General Mathematics Impact factor: 0.340, year: 2013 http://www.degruyter.com/view/j/gmj.2013.20.issue-1/gmj-2013-0001/gmj-2013-0001. xml ?format=INT
Remark on zeros of solutions of second-order linear ordinary differential equations
Czech Academy of Sciences Publication Activity Database
Dosoudilová, M.; Lomtatidze, Alexander
2016-01-01
Roč. 23, č. 4 (2016), s. 571-577 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : second-order linear equation * zeros of solutions * periodic boundary value problem Subject RIV: BA - General Mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2016.23.issue-4/gmj-2016-0052/gmj-2016-0052. xml
Exact Solution of Mutator Model with Linear Fitness and Finite Genome Length
Saakian, David B.
2017-08-01
We considered the infinite population version of the mutator phenomenon in evolutionary dynamics, looking at the uni-directional mutations in the mutator-specific genes and linear selection. We solved exactly the model for the finite genome length case, looking at the quasispecies version of the phenomenon. We calculated the mutator probability both in the statics and dynamics. The exact solution is important for us because the mutator probability depends on the genome length in a highly non-trivial way.
A General Construction of Linear Differential Equations with Solutions of Prescribed Properties
Czech Academy of Sciences Publication Activity Database
Neuman, František
2004-01-01
Roč. 17, č. 1 (2004), s. 71-76 ISSN 0893-9659 R&D Projects: GA AV ČR IAA1019902; GA ČR GA201/99/0295 Institutional research plan: CEZ:AV0Z1019905 Keywords : construction of linear differential equations * prescribed qualitative properties of solutions Subject RIV: BA - General Mathematics Impact factor: 0.414, year: 2004
Remark on zeros of solutions of second-order linear ordinary differential equations
Czech Academy of Sciences Publication Activity Database
Dosoudilová, M.; Lomtatidze, Alexander
2016-01-01
Roč. 23, č. 4 (2016), s. 571-577 ISSN 1072-947X Institutional support: RVO:67985840 Keywords : second-order linear equation * zero s of solutions * periodic boundary value problem Subject RIV: BA - General Mathematics Impact factor: 0.290, year: 2016 https://www.degruyter.com/view/j/gmj.2016.23.issue-4/gmj-2016-0052/gmj-2016-0052.xml
An algorithm for computing the hull of the solution set of interval linear equations
Czech Academy of Sciences Publication Activity Database
Rohn, Jiří
2011-01-01
Roč. 435, č. 2 (2011), s. 193-201 ISSN 0024-3795 R&D Projects: GA ČR GA201/09/1957; GA ČR GC201/08/J020 Institutional research plan: CEZ:AV0Z10300504 Keywords : interval linear equations * solution set * interval hull * algorithm * absolute value inequality Subject RIV: BA - General Mathematics Impact factor: 0.974, year: 2011
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.
Directory of Open Access Journals (Sweden)
Eusebio Eduardo Hernández Martinez
2013-01-01
Full Text Available In robotics, solving the direct kinematics problem (DKP for parallel robots is very often more difficult and time consuming than for their serial counterparts. The problem is stated as follows: given the joint variables, the Cartesian variables should be computed, namely the pose of the mobile platform. Most of the time, the DKP requires solving a non-linear system of equations. In addition, given that the system could be non-convex, Newton or Quasi-Newton (Dogleg based solvers get trapped on local minima. The capacity of such kinds of solvers to find an adequate solution strongly depends on the starting point. A well-known problem is the selection of such a starting point, which requires a priori information about the neighbouring region of the solution. In order to circumvent this issue, this article proposes an efficient method to select and to generate the starting point based on probabilistic learning. Experiments and discussion are presented to show the method performance. The method successfully avoids getting trapped on local minima without the need for human intervention, which increases its robustness when compared with a single Dogleg approach. This proposal can be extended to other structures, to any non-linear system of equations, and of course, to non-linear optimization problems.
Impact of quadratic non-linearity on the dynamics of periodic solutions of a wave equation
International Nuclear Information System (INIS)
Kolesov, Andrei Yu; Rozov, Nikolai Kh
2002-01-01
For the non-linear telegraph equation with homogeneous Dirichlet or Neumann conditions at the end-points of a finite interval the question of the existence and the stability of time-periodic solutions bifurcating from the zero equilibrium state is considered. The dynamics of these solutions under a change of the diffusion coefficient (that is, the coefficient of the second derivative with respect to the space variable) is investigated. For the Dirichlet boundary conditions it is shown that this dynamics substantially depends on the presence - or the absence - of quadratic terms in the non-linearity. More precisely, it is shown that a quadratic non-linearity results in the occurrence, under an unbounded decrease of diffusion, of an infinite sequence of bifurcations of each periodic solution. En route, the related issue of the limits of applicability of Yu.S. Kolesov's method of quasinormal forms to the construction of self-oscillations in singularly perturbed hyperbolic boundary value problems is studied
Energy Technology Data Exchange (ETDEWEB)
Setoguchi, T.; Manchu, Y.; Katsumata, M. [Toshiba Corp., Tokyo (Japan)
2000-04-01
Toshiba provides a range of information technology (IT) solutions called SmartEC Solution, which includes business-to-business electronic commerce systems and services based on international standards and industrial know-how, especially our electronic data interchange (EDI) know-how as a manufacturer. These IT solutions are supplied as services covering strategy planning, system integration, and application service provider based on five types of business-to-business electronic commerce. (author)
A parametric LTR solution for discrete-time systems
DEFF Research Database (Denmark)
Niemann, Hans Henrik; Jannerup, Ole Erik
1989-01-01
A parametric LTR (loop transfer recovery) solution for discrete-time compensators incorporating filtering observers which achieve exact recovery is presented for both minimum- and non-minimum-phase systems. First the recovery error, which defines the difference between the target loop transfer...... and the full loop transfer function, is manipulated into a general form involving the target loop transfer matrix and the fundamental recovery matrix. A parametric LTR solution based on the recovery matrix is developed. It is shown that the LQR/LTR (linear quadratic Gaussian/loop transfer recovery) solution...
Krylov Subspace Methods for Complex Non-Hermitian Linear Systems. Thesis
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.
On Attainability of Optimal Solutions for Linear Elliptic Equations with Unbounded Coefficients
Directory of Open Access Journals (Sweden)
P. I. Kogut
2011-12-01
Full Text Available We study an optimal boundary control problem (OCP associated to a linear elliptic equation —div (Vj/ + A(xVy = f describing diffusion in a turbulent flow. The characteristic feature of this equation is the fact that, in applications, the stream matrix A(x = [a,ij(x]i,j=i,...,N is skew-symmetric, ац(х = —a,ji(x, measurable, and belongs to L -space (rather than L°°. An optimal solution to such problem can inherit a singular character of the original stream matrix A. We show that optimal solutions can be attainable by solutions of special optimal boundary control problems.
The structure of solutions of the matrix linear unilateral polynomial equation with two variables
Directory of Open Access Journals (Sweden)
N. S. Dzhaliuk
2017-07-01
Full Text Available We investigate the structure of solutions of the matrix linear polynomial equation $A(\\lambdaX(\\lambda+B(\\lambdaY(\\lambda=C(\\lambda,$ in particular, possible degrees of the solutions. The solving of this equation is reduced to the solving of the equivalent matrix polynomial equation with matrix coefficients in triangular forms with invariant factors on the main diagonals, to which the matrices $A (\\lambda, B(\\lambda$ \\ and \\ $C(\\lambda$ are reduced by means of semiscalar equivalent transformations. On the basis of it, we have pointed out the bounds of the degrees of the matrix polynomial equation solutions. Necessary and sufficient conditions for the uniqueness of a solution with a minimal degree are established. An effective method for constructing minimal degree solutions of the equations is suggested. In this article, unlike well-known results about the estimations of the degrees of the solutions of the matrix polynomial equations in which both matrix coefficients are regular or at least one of them is regular, we have considered the case when the matrix polynomial equation has arbitrary matrix coefficients $A(\\lambda$ and $B(\\lambda.$
Spherically symmetric analysis on open FLRW solution in non-linear massive gravity
Energy Technology Data Exchange (ETDEWEB)
Chiang, Chien-I; Izumi, Keisuke; Chen, Pisin, E-mail: chienichiang@berkeley.edu, E-mail: izumi@phys.ntu.edu.tw, E-mail: chen@slac.stanford.edu [Leung Center for Cosmology and Particle Astrophysics, National Taiwan University, Taipei 10617, Taiwan (China)
2012-12-01
We study non-linear massive gravity in the spherically symmetric context. Our main motivation is to investigate the effect of helicity-0 mode which remains elusive after analysis of cosmological perturbation around an open Friedmann-Lemaitre-Robertson-Walker (FLRW) universe. The non-linear form of the effective energy-momentum tensor stemming from the mass term is derived for the spherically symmetric case. Only in the special case where the area of the two sphere is not deviated away from the FLRW universe, the effective energy momentum tensor becomes completely the same as that of cosmological constant. This opens a window for discriminating the non-linear massive gravity from general relativity (GR). Indeed, by further solving these spherically symmetric gravitational equations of motion in vacuum to the linear order, we obtain a solution which has an arbitrary time-dependent parameter. In GR, this parameter is a constant and corresponds to the mass of a star. Our result means that Birkhoff's theorem no longer holds in the non-linear massive gravity and suggests that energy can probably be emitted superluminously (with infinite speed) on the self-accelerating background by the helicity-0 mode, which could be a potential plague of this theory.
Krishnan, M.
2017-05-01
We present a model for calculating the net and effective electrical charge of globular macromolecules and linear polyelectrolytes such as proteins and DNA, given the concentration of monovalent salt and pH in solution. The calculation is based on a numerical solution of the non-linear Poisson-Boltzmann equation using a finite element discretized continuum approach. The model simultaneously addresses the phenomena of charge regulation and renormalization, both of which underpin the electrostatics of biomolecules in solution. We show that while charge regulation addresses the true electrical charge of a molecule arising from the acid-base equilibria of its ionizable groups, charge renormalization finds relevance in the context of a molecule's interaction with another charged entity. Writing this electrostatic interaction free energy in terms of a local electrical potential, we obtain an "interaction charge" for the molecule which we demonstrate agrees closely with the "effective charge" discussed in charge renormalization and counterion-condensation theories. The predictions of this model agree well with direct high-precision measurements of effective electrical charge of polyelectrolytes such as nucleic acids and disordered proteins in solution, without tunable parameters. Including the effective interior dielectric constant for compactly folded molecules as a tunable parameter, the model captures measurements of effective charge as well as published trends of pKa shifts in globular proteins. Our results suggest a straightforward general framework to model electrostatics in biomolecules in solution. In offering a platform that directly links theory and experiment, these calculations could foster a systematic understanding of the interrelationship between molecular 3D structure and conformation, electrical charge and electrostatic interactions in solution. The model could find particular relevance in situations where molecular crystal structures are not available or
International Nuclear Information System (INIS)
Qin, Hong; Davidson, Ronald C.
2011-01-01
In a linear trap confining a one-component nonneutral plasma, the external focusing force is a linear function of the configuration coordinates and/or the velocity coordinates. Linear traps include the classical Paul trap and the Penning trap, as well as the newly proposed rotating-radio- frequency traps and the Mobius accelerator. This paper describes a class of self-similar nonlinear solutions of nonneutral plasma in general time-dependent linear focusing devices, with self-consistent electrostatic field. This class of nonlinear solutions includes many known solutions as special cases.
A Linear Active Disturbance Rejection Control for a Ball and Rigid Triangle System
Directory of Open Access Journals (Sweden)
Carlos Aguilar-Ibanez
2016-01-01
Full Text Available This paper proposes an application of linear flatness control along with active disturbance rejection control (ADRC for the local stabilization and trajectory tracking problems in the underactuated ball and rigid triangle system. To this end, an observer-based linear controller of the ADRC type is designed based on the flat tangent linearization of the system around its corresponding unstable equilibrium rest position. It was accomplished through two decoupled linear extended observers and a single linear output feedback controller, with disturbance cancelation features. The controller guarantees locally exponentially asymptotic stability for the stabilization problem and practical local stability in the solution of the tracking error. An advantage of combining the flatness and the ADRC methods is that it possible to perform online estimates and cancels the undesirable effects of the higher-order nonlinearities discarded by the linearization approximation. Simulation indicates that the proposed controller behaves remarkably well, having an acceptable domain of attraction.
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.
Periodic inventory system in cafeteria using linear programming
Usop, Mohd Fais; Ishak, Ruzana; Hamdan, Ahmad Ridhuan
2017-11-01
Inventory management is an important factor in running a business. It plays a big role of managing the stock in cafeteria. If the inventories are failed to be managed wisely, it will affect the profit of the cafeteria. Therefore, the purpose of this study is to find the solution of the inventory management in cafeteria. Most of the cafeteria in Malaysia did not manage their stock well. Therefore, this study is to propose a database system of inventory management and to develop the inventory model in cafeteria management. In this study, new database system to improve the management of the stock in a weekly basis will be provided using Linear Programming Model to get the optimal range of the inventory needed for selected categories. Data that were collected by using the Periodic Inventory System at the end of the week within three months period being analyzed by using the Food Stock-take Database. The inventory model was developed from the collected data according to the category of the inventory in the cafeteria. Results showed the effectiveness of using the Periodic Inventory System and will be very helpful to the cafeteria management in organizing the inventory. Moreover, the findings in this study can reduce the cost of operation and increased the profit.
Partial Linearization of Mechanical Systems with Application to Observer Design
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
Simultaneous Balancing and Model Reduction of Switched Linear Systems
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
International Nuclear Information System (INIS)
Eckstein, U.; Harte, R.; Kraetzig, W.B.; Wittek, U.
1983-01-01
In order to describe nonlinear response and instability behaviour the paper starts with the total potential energy considering the basic kinematic equations of a consistent nonlinear shell theory for large displacements and moderate rotations. The material behaviour is assumed to be hyperelastic and isotropic. The incrementation and discretization of the total potential energy leads to the tangent stiffness relation, which is the central equation of computational algorithms based on combined incremental and iterative techniques. Here a symmetrized form of the RIKS/WEMPNER-algorithm for positive and negative load incrementation represents the basis of the nonlinear solution technique. To detect secondary equilibrium branches at points of neutral equilibrium within nonlinear primary paths a quadratic eigenvalue-problem has to be solved. In order to follow those complicated nonlinear response phenomena the RIKS/WEMPNER incrementation/iteration process is combined with a simultaneous solution of the linearized quadratic eigenvalue-problem. Additionally the essentials of a recently derived family of arbitrarily curved shell elements for linear (LACS) and geometrically nonlinear (NACS) shell problems are presented. The main advantage of these elements is the exact description of all geometric properties as well as the energy-equivalent representation of the applied loads in combination with an efficient algorithm to form the stiffness submatrices. Especially the NACS-elements are designed to improve the accuracy of the solution in the deep postbuckling range including moderate rotations. The derived finite elements and solution strategies are applied to a certain number of typical shell problems to prove the precision of the shell elements and to demonstrate the possibilities of tracing linear and nonlinear bifurcation problems as well as snap-through phenomena with and without secondary bifurcation branches. (orig.)
Linear System Control Using Stochastic Learning Automata
Ziyad, Nigel; Cox, E. Lucien; Chouikha, Mohamed F.
1998-01-01
This paper explains the use of a Stochastic Learning Automata (SLA) to control switching between three systems to produce the desired output response. The SLA learns the optimal choice of the damping ratio for each system to achieve a desired result. We show that the SLA can learn these states for the control of an unknown system with the proper choice of the error criteria. The results of using a single automaton are compared to using multiple automata.
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
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-.
Robust Hinf control of uncertain switched systems defined on polyhedral sets with Filippov solutions
DEFF Research Database (Denmark)
Ahmadi, Mohamadreza; Mojallali, Hamed; Wisniewski, Rafal
2012-01-01
This paper considers the control problem of a class of uncertain switched systems defined on polyhedral sets known as piecewise linear systems where, instead of the conventional Carathe ́odory solutions, Filippov solutions are studied. In other words, in contrast to the previous studies, solutions...
Analytical Solution for Optimum Design of Furrow Irrigation Systems
Kiwan, M. E.
1996-05-01
An analytical solution for the optimum design of furrow irrigation systems is derived. The non-linear calculus optimization method is used to formulate a general form for designing the optimum system elements under circumstances of maximizing the water application efficiency of the system during irrigation. Different system bases and constraints are considered in the solution. A full irrigation water depth is considered to be achieved at the tail of the furrow line. The solution is based on neglecting the recession and depletion times after off-irrigation. This assumption is valid in the case of open-end (free gradient) furrow systems rather than closed-end (closed dike) systems. Illustrative examples for different systems are presented and the results are compared with the output obtained using an iterative numerical solution method. The final derived solution is expressed as a function of the furrow length ratio (the furrow length to the water travelling distance). The function of water travelling developed by Reddy et al. is considered for reaching the optimum solution. As practical results from the study, the optimum furrow elements for free gradient systems can be estimated to achieve the maximum application efficiency, i.e. furrow length, water inflow rate and cutoff irrigation time.
A solution to the varying response of the linear power monitor induced by xenon poisoning
Energy Technology Data Exchange (ETDEWEB)
Godsey, T A; Randall, J D [Texas A and M University (United States)
1974-07-01
After conversion to FLIP fuel at Texas A and M, the fuel temperatures were examined very carefully. It was observed that the fuel temperature at 1 Mw varied over a wide range during the week. This variation was shown to be due to the variation in response of the linear CIC which was used to establish reactor power level. A modification of the linear power monitor was designed and installed. The response of this system was verified by using cobalt wires, fuel temperature, and a fission chamber located at 6 feet from the reactor core. The system has proven to be operationally satisfactory. (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...
Virtual Estimator for Piecewise Linear Systems Based on Observability Analysis
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
Linear wave systems on n-D spatial domains
Kurula, Mikael; Zwart, Heiko J.
2015-01-01
In this paper, we study the linear wave equation on an n-dimensional spatial domain.We show that there is a boundary triplet associated to the undamped wave equation. This enables us to characterise all boundary conditions for which the undamped wave equation possesses a unique solution
Stability of Bifurcating Stationary Solutions of the Artificial Compressible System
Teramoto, Yuka
2018-02-01
The artificial compressible system gives a compressible approximation of the incompressible Navier-Stokes system. The latter system is obtained from the former one in the zero limit of the artificial Mach number ɛ which is a singular limit. The sets of stationary solutions of both systems coincide with each other. It is known that if a stationary solution of the incompressible system is asymptotically stable and the velocity field of the stationary solution satisfies an energy-type stability criterion, then it is also stable as a solution of the artificial compressible one for sufficiently small ɛ . In general, the range of ɛ shrinks when the spectrum of the linearized operator for the incompressible system approaches to the imaginary axis. This can happen when a stationary bifurcation occurs. It is proved that when a stationary bifurcation from a simple eigenvalue occurs, the range of ɛ can be taken uniformly near the bifurcation point to conclude the stability of the bifurcating solution as a solution of the artificial compressible system.
Fu, Wei; Nijhoff, Frank W
2017-07-01
A unified framework is presented for the solution structure of three-dimensional discrete integrable systems, including the lattice AKP, BKP and CKP equations. This is done through the so-called direct linearizing transform, which establishes a general class of integral transforms between solutions. As a particular application, novel soliton-type solutions for the lattice CKP equation are obtained.
An overview of solution methods for multi-objective mixed integer linear programming programs
DEFF Research Database (Denmark)
Andersen, Kim Allan; Stidsen, Thomas Riis
Multiple objective mixed integer linear programming (MOMIP) problems are notoriously hard to solve to optimality, i.e. finding the complete set of non-dominated solutions. We will give an overview of existing methods. Among those are interactive methods, the two phases method and enumeration...... methods. In particular we will discuss the existing branch and bound approaches for solving multiple objective integer programming problems. Despite the fact that branch and bound methods has been applied successfully to integer programming problems with one criterion only a few attempts has been made...
On the existence of tunneling bounce solutions in piecewise linear potentials
International Nuclear Information System (INIS)
Dutta, Koushik; Hector, Cecelie; Konstandin, Thomas; Vaudrevange, Pascal M.; Westphal, Alexander
2012-02-01
Coleman tunneling in a general scalar potential with two non-degenerate minima is known to have an approximation in terms of a piecewise linear triangular-shaped potential with sharp 'kinks' at the place of the local minima. This approximate potential has a regime where the existence of the bounce solution needs the scalar field to 'wait' for some amount of Euclidean time at one of the 'kinks'. We discuss under which circumstances the correct bounce action can be consistently obtained as the limiting case of a regular scalar potential where 'kinks' are resolved as locally smooth 'cap' regions. (orig.)
Gradient remediability in linear distributed parabolic systems ...
African Journals Online (AJOL)
The aim of this paper is the introduction of a new concept that concerned the analysis of a large class of distributed parabolic systems. It is the general concept of gradient remediability. More precisely, we study with respect to the gradient observation, the existence of an input operator (gradient efficient actuators) ensuring ...
Exact solutions for a system of nonlinear plasma fluid equations
International Nuclear Information System (INIS)
Prahovic, M.G.; Hazeltine, R.D.; Morrison, P.J.
1991-04-01
A method is presented for constructing exact solutions to a system of nonlinear plasma fluid equations that combines the physics of reduced magnetohydrodynamics and the electrostatic drift-wave description of the Charney-Hasegawa-Mima equation. The system has nonlinearities that take the form of Poisson brackets involving the fluid field variables. The method relies on modifying a class of simple equilibrium solutions, but no approximations are made. A distinguishing feature is that the original nonlinear problem is reduced to the solution of two linear partial differential equations, one fourth-order and the other first-order. The first-order equation has Hamiltonian characteristics and is easily integrated, supplying information about the general structure of solutions. 6 refs
Growth and Zeros of Meromorphic Solutions to Second-Order Linear Differential Equations
Directory of Open Access Journals (Sweden)
Maamar Andasmas
2016-04-01
Full Text Available The main purpose of this article is to investigate the growth of meromorphic solutions to homogeneous and non-homogeneous second order linear differential equations f00+Af0+Bf = F, where A(z, B (z and F (z are meromorphic functions with finite order having only finitely many poles. We show that, if there exist a positive constants σ > 0, α > 0 such that |A(z| ≥ eα|z|σ as |z| → +∞, z ∈ H, where dens{|z| : z ∈ H} > 0 and ρ = max{ρ(B, ρ(F} < σ, then every transcendental meromorphic solution f has an infinite order. Further, we give some estimates of their hyper-order, exponent and hyper-exponent of convergence of distinct zeros.
Phantom solution in a non-linear Israel-Stewart theory
Cruz, Miguel; Cruz, Norman; Lepe, Samuel
2017-06-01
In this paper we present a phantom solution with a big rip singularity in a non-linear regime of the Israel-Stewart formalism. In this framework it is possible to extend this causal formalism in order to describe accelerated expansion, where assumption of near equilibrium is no longer valid. We assume a flat universe filled with a single viscous fluid ruled by a barotropic EoS, p = ωρ, which can represent a late time accelerated phase of the cosmic evolution. The solution allows to cross the phantom divide without evoking an exotic matter fluid and the effective EoS parameter is always lesser than -1 and constant in time.
ORACLS: A system for linear-quadratic-Gaussian control law design
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.
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.
Clemens, Joshua William
Game theory has application across multiple fields, spanning from economic strategy to optimal control of an aircraft and missile on an intercept trajectory. The idea of game theory is fascinating in that we can actually mathematically model real-world scenarios and determine optimal decision making. It may not always be easy to mathematically model certain real-world scenarios, nonetheless, game theory gives us an appreciation for the complexity involved in decision making. This complexity is especially apparent when the players involved have access to different information upon which to base their decision making (a nonclassical information pattern). Here we will focus on the class of adversarial two-player games (sometimes referred to as pursuit-evasion games) with nonclassical information pattern. We present a two-sided (simultaneous) optimization solution method for the two-player linear quadratic Gaussian (LQG) multistage game. This direct solution method allows for further interpretation of each player's decision making (strategy) as compared to previously used formal solution methods. In addition to the optimal control strategies, we present a saddle point proof and we derive an expression for the optimal performance index value. We provide some numerical results in order to further interpret the optimal control strategies and to highlight real-world application of this game-theoretic optimal solution.
On the discretization of linear fractional representations of LPV systems
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
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
DEFF Research Database (Denmark)
Ommen, Torben Schmidt; Markussen, Wiebke Brix; Elmegaard, Brian
2014-01-01
In the paper, three frequently used operation optimisation methods are examined with respect to their impact on operation management of the combined utility technologies for electric power and DH (district heating) of eastern Denmark. The investigation focusses on individual plant operation...... differences and differences between the solution found by each optimisation method. One of the investigated approaches utilises LP (linear programming) for optimisation, one uses LP with binary operation constraints, while the third approach uses NLP (non-linear programming). The LP model is used...... as a benchmark, as this type is frequently used, and has the lowest amount of constraints of the three. A comparison of the optimised operation of a number of units shows significant differences between the three methods. Compared to the reference, the use of binary integer variables, increases operation...
Feedback linearizing control of a MIMO power system
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.
Directory of Open Access Journals (Sweden)
Nahed S. Hussein
2014-01-01
Full Text Available A numerical boundary integral scheme is proposed for the solution to the system of eld equations of plane. The stresses are prescribed on one-half of the circle, while the displacements are given. The considered problem with mixed boundary conditions in the circle is replaced by two problems with homogeneous boundary conditions, one of each type, having a common solution. The equations are reduced to a system of boundary integral equations, which is then discretized in the usual way, and the problem at this stage is reduced to the solution to a rectangular linear system of algebraic equations. The unknowns in this system of equations are the boundary values of four harmonic functions which define the full elastic solution and the unknown boundary values of stresses or displacements on proper parts of the boundary. On the basis of the obtained results, it is inferred that a stress component has a singularity at each of the two separation points, thought to be of logarithmic type. The results are discussed and boundary plots are given. We have also calculated the unknown functions in the bulk directly from the given boundary conditions using the boundary collocation method. The obtained results in the bulk are discussed and three-dimensional plots are given. A tentative form for the singular solution is proposed and the corresponding singular stresses and displacements are plotted in the bulk. The form of the singular tangential stress is seen to be compatible with the boundary values obtained earlier. The efficiency of the used numerical schemes is discussed.
Portable, x-band, linear accelerator systems
International Nuclear Information System (INIS)
Schonberg, R.G.; Deruyter, H.; Fowkes, W.R.; Johnson, W.A.; Miller, R.H.; Potter, J.M.; Weaver, J.N.
1985-01-01
Three light-weight, x-band, electron accelerators have been developed to provide a series of highly portable sources of x-rays and neutrons for nondestructive testing. The 1.5 MeV x-ray unit has a 200 kW magnetron for an RF source and an air-cooled, traveling wave accelerating structure to minimize its weight. The 4 and 6 MeV units share the same drive system which contains a 1.2 MW magnetron. The 4 MeV unit uses a traveling-wave guide to produce x-rays and the 6MeV unit uses a standing-wave guide to produce x-rays or neutrons. The choice of 9.3 GHz was dictated by the availability of a high power coaxial magnetron and by the obvious dimensional and weight advantages of a higher frequency over the more common S-band frequencies around 3 GHz
Stability margin of linear systems with parameters described by fuzzy numbers.
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.
An implicit meshless scheme for the solution of transient non-linear Poisson-type equations
Bourantas, Georgios
2013-07-01
A meshfree point collocation method is used for the numerical simulation of both transient and steady state non-linear Poisson-type partial differential equations. Particular emphasis is placed on the application of the linearization method with special attention to the lagging of coefficients method and the Newton linearization method. The localized form of the Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are performed for regular nodal distributions, stressing the positivity conditions that make the resulting system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through representative and well-established benchmark problems. © 2013 Elsevier Ltd.
An implicit meshless scheme for the solution of transient non-linear Poisson-type equations
Bourantas, Georgios; Burganos, Vasilis N.
2013-01-01
A meshfree point collocation method is used for the numerical simulation of both transient and steady state non-linear Poisson-type partial differential equations. Particular emphasis is placed on the application of the linearization method with special attention to the lagging of coefficients method and the Newton linearization method. The localized form of the Moving Least Squares (MLS) approximation is employed for the construction of the shape functions, in conjunction with the general framework of the point collocation method. Computations are performed for regular nodal distributions, stressing the positivity conditions that make the resulting system stable and convergent. The accuracy and the stability of the proposed scheme are demonstrated through representative and well-established benchmark problems. © 2013 Elsevier Ltd.
Structure Learning in Stochastic Non-linear Dynamical Systems
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.
Portable, x-band, linear accelerator systems
International Nuclear Information System (INIS)
Schonberg, R.G.; Deruyter, H.; Fowkes, W.R.; Johnson, W.A.; Miller, R.H.; Potter, J.M.; Weaver, J.N.
1985-01-01
Three light-weight, x-band, electron accelerators have been developed to provide a series of highly portable sources of x-rays and neutrons for non-destructive testing. The 1.5 MeV x-ray unit has a 200 kW magnetron for an RF source and an air-cooled, traveling wave accelerating structure to minimize its weight. The 4 and 6 MeV units share the same drive system which contains a 1.2 MW magnetron. The 4 MeV unit uses a traveling-wave guide to produce x-rays and the 6MeV unit uses a standing-wave guide to produce x-rays or neutrons. The choice of 9.3 GHz was dictated by the availability of a high power coaxial magnetron and by the obvious dimensional and weight advantages of a higher frequency over the more common S-band frequencies around 3 GHz
A neural network method for solving a system of linear variational inequalities
International Nuclear Information System (INIS)
Lan Hengyou; Cui Yishun
2009-01-01
In this paper, we transmute the solution for a new system of linear variational inequalities to an equilibrium point of neural networks, and by using analytic technique, some sufficient conditions are presented. Further, the estimation of the exponential convergence rates of the neural networks is investigated. The new and useful results obtained in this paper generalize and improve the corresponding results of recent works.
Exponential stability of switched linear systems with time-varying delay
Directory of Open Access Journals (Sweden)
Satiracoo Pairote
2007-11-01
Full Text Available We use a Lyapunov-Krasovskii functional approach to establish the exponential stability of linear systems with time-varying delay. Our delay-dependent condition allows to compute simultaneously the two bounds that characterize the exponential stability rate of the solution. A simple procedure for constructing switching rule is also presented.
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
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.
Iterative algorithms for large sparse linear systems on parallel computers
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.
Simultaneous Balancing and Model Reduction of Switched Linear Systems
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
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...
Uniqueness of solutions of relay systems, Special Issue on Hybrid Systems
Lootsma, Y.J.; van der Schaft, Arjan; Camlıbel, M.K.
1999-01-01
Conditions are given for uniqueness of solutions of linear time-invariant systems under relay feedback. From a hybrid dynamical point of view this entails the deterministic specification of the discrete transition rules. The results are based on the formulation of relay systems as complementarity
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
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.
Quantitative Pointwise Estimate of the Solution of the Linearized Boltzmann Equation
Lin, Yu-Chu; Wang, Haitao; Wu, Kung-Chien
2018-04-01
We study the quantitative pointwise behavior of the solutions of the linearized Boltzmann equation for hard potentials, Maxwellian molecules and soft potentials, with Grad's angular cutoff assumption. More precisely, for solutions inside the finite Mach number region (time like region), we obtain the pointwise fluid structure for hard potentials and Maxwellian molecules, and optimal time decay in the fluid part and sub-exponential time decay in the non-fluid part for soft potentials. For solutions outside the finite Mach number region (space like region), we obtain sub-exponential decay in the space variable. The singular wave estimate, regularization estimate and refined weighted energy estimate play important roles in this paper. Our results extend the classical results of Liu and Yu (Commun Pure Appl Math 57:1543-1608, 2004), (Bull Inst Math Acad Sin 1:1-78, 2006), (Bull Inst Math Acad Sin 6:151-243, 2011) and Lee et al. (Commun Math Phys 269:17-37, 2007) to hard and soft potentials by imposing suitable exponential velocity weight on the initial condition.
Quantitative Pointwise Estimate of the Solution of the Linearized Boltzmann Equation
Lin, Yu-Chu; Wang, Haitao; Wu, Kung-Chien
2018-06-01
We study the quantitative pointwise behavior of the solutions of the linearized Boltzmann equation for hard potentials, Maxwellian molecules and soft potentials, with Grad's angular cutoff assumption. More precisely, for solutions inside the finite Mach number region (time like region), we obtain the pointwise fluid structure for hard potentials and Maxwellian molecules, and optimal time decay in the fluid part and sub-exponential time decay in the non-fluid part for soft potentials. For solutions outside the finite Mach number region (space like region), we obtain sub-exponential decay in the space variable. The singular wave estimate, regularization estimate and refined weighted energy estimate play important roles in this paper. Our results extend the classical results of Liu and Yu (Commun Pure Appl Math 57:1543-1608, 2004), (Bull Inst Math Acad Sin 1:1-78, 2006), (Bull Inst Math Acad Sin 6:151-243, 2011) and Lee et al. (Commun Math Phys 269:17-37, 2007) to hard and soft potentials by imposing suitable exponential velocity weight on the initial condition.
Error compensation for hybrid-computer solution of linear differential equations
Kemp, N. H.
1970-01-01
Z-transform technique compensates for digital transport delay and digital-to-analog hold. Method determines best values for compensation constants in multi-step and Taylor series projections. Technique also provides hybrid-calculation error compared to continuous exact solution, plus system stability properties.
International Nuclear Information System (INIS)
Bouaziz, M.N.; Aziz, Abdul
2010-01-01
A novel concept of double optimal linearization is introduced and used to obtain a simple and accurate solution for the temperature distribution in a straight rectangular convective-radiative fin with temperature dependent thermal conductivity. The solution is built from the classical solution for a pure convection fin of constant thermal conductivity which appears in terms of hyperbolic functions. When compared with the direct numerical solution, the double optimally linearized solution is found to be accurate within 4% for a range of radiation-conduction and thermal conductivity parameters that are likely to be encountered in practice. The present solution is simple and offers superior accuracy compared with the fairly complex approximate solutions based on the homotopy perturbation method, variational iteration method, and the double series regular perturbation method. The fin efficiency expression resembles the classical result for the constant thermal conductivity convecting fin. The present results are easily usable by the practicing engineers in their thermal design and analysis work involving fins.
International Nuclear Information System (INIS)
Belmonte-Beitia, Juan; Konotop, Vladimir V.; Perez-Garcia, Victor M.; Vekslerchik, Vadym E.
2009-01-01
Using similarity transformations we construct explicit solutions of the nonlinear Schroedinger equation with linear and nonlinear periodic potentials. We present explicit forms of spatially localized and periodic solutions, and study their properties. We put our results in the framework of the exploited perturbation techniques and discuss their implications on the properties of associated linear periodic potentials and on the possibilities of stabilization of gap solitons using polychromatic lattices.
Global stability and exact solution of an arbitrary-solute nonlinear cellular mass transport system.
Benson, James D
2014-12-01
The prediction of the cellular state as a function of extracellular concentrations and temperatures has been of interest to physiologists for nearly a century. One of the most widely used models in the field is one where mass flux is linearly proportional to the concentration difference across the membrane. These fluxes define a nonlinear differential equation system for the intracellular state, which when coupled with appropriate initial conditions, define the intracellular state as a function of the extracellular concentrations of both permeating and nonpermeating solutes. Here we take advantage of a reparametrization scheme to extend existing stability results to a more general setting and to a develop analytical solutions to this model for an arbitrary number of extracellular solutes. Copyright © 2014 Elsevier Inc. All rights reserved.
Solutions to estimation problems for scalar hamilton-jacobi equations using linear programming
Claudel, Christian G.; Chamoin, Timothee; Bayen, Alexandre M.
2014-01-01
This brief presents new convex formulations for solving estimation problems in systems modeled by scalar Hamilton-Jacobi (HJ) equations. Using a semi-analytic formula, we show that the constraints resulting from a HJ equation are convex, and can be written as a set of linear inequalities. We use this fact to pose various (and seemingly unrelated) estimation problems related to traffic flow-engineering as a set of linear programs. In particular, we solve data assimilation and data reconciliation problems for estimating the state of a system when the model and measurement constraints are incompatible. We also solve traffic estimation problems, such as travel time estimation or density estimation. For all these problems, a numerical implementation is performed using experimental data from the Mobile Century experiment. In the context of reproducible research, the code and data used to compute the results presented in this brief have been posted online and are accessible to regenerate the results. © 2013 IEEE.
Stability analysis of switched linear systems defined by graphs
Athanasopoulos, N.; Lazar, M.
2014-01-01
We present necessary and sufficient conditions for global exponential stability for switched discrete-time linear systems, under arbitrary switching, which is constrained within a set of admissible transitions. The class of systems studied includes the family of systems under arbitrary switching,
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 ...
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...
Linear Optimization of Frequency Spectrum Assignments Across System
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
International Nuclear Information System (INIS)
Gene Golub; Kwok Ko
2009-01-01
The solutions of sparse eigenvalue problems and linear systems constitute one of the key computational kernels in the discretization of partial differential equations for the modeling of linear accelerators. The computational challenges faced by existing techniques for solving those sparse eigenvalue problems and linear systems call for continuing research to improve on the algorithms so that ever increasing problem size as required by the physics application can be tackled. Under the support of this award, the filter algorithm for solving large sparse eigenvalue problems was developed at Stanford to address the computational difficulties in the previous methods with the goal to enable accelerator simulations on then the world largest unclassified supercomputer at NERSC for this class of problems. Specifically, a new method, the Hemitian skew-Hemitian splitting method, was proposed and researched as an improved method for solving linear systems with non-Hermitian positive definite and semidefinite matrices.
Dynamic stability of a vertically excited non-linear continuous system
Czech Academy of Sciences Publication Activity Database
Náprstek, Jiří; Fischer, Cyril
2015-01-01
Roč. 155, July (2015), s. 106-114 ISSN 0045-7949 R&D Projects: GA ČR(CZ) GA15-01035S Institutional support: RVO:68378297 Keywords : non-linear systems * auto-parametric systems * semi-trivial solution * dynamic stability * system recovery * post- critical response Subject RIV: JM - Building Engineering Impact factor: 2.425, year: 2015 http://www.sciencedirect.com/science/article/pii/S0045794915000024
Stochastic Least-Squares Petrov--Galerkin Method for Parameterized Linear Systems
Energy Technology Data Exchange (ETDEWEB)
Lee, Kookjin [Univ. of Maryland, College Park, MD (United States). Dept. of Computer Science; Carlberg, Kevin [Sandia National Lab. (SNL-CA), Livermore, CA (United States); Elman, Howard C. [Univ. of Maryland, College Park, MD (United States). Dept. of Computer Science and Inst. for Advanced Computer Studies
2018-03-29
Here, we consider the numerical solution of parameterized linear systems where the system matrix, the solution, and the right-hand side are parameterized by a set of uncertain input parameters. We explore spectral methods in which the solutions are approximated in a chosen finite-dimensional subspace. It has been shown that the stochastic Galerkin projection technique fails to minimize any measure of the solution error. As a remedy for this, we propose a novel stochatic least-squares Petrov--Galerkin (LSPG) method. The proposed method is optimal in the sense that it produces the solution that minimizes a weighted $\\ell^2$-norm of the residual over all solutions in a given finite-dimensional subspace. Moreover, the method can be adapted to minimize the solution error in different weighted $\\ell^2$-norms by simply applying a weighting function within the least-squares formulation. In addition, a goal-oriented seminorm induced by an output quantity of interest can be minimized by defining a weighting function as a linear functional of the solution. We establish optimality and error bounds for the proposed method, and extensive numerical experiments show that the weighted LSPG method outperforms other spectral methods in minimizing corresponding target weighted norms.
Fractional order differentiation by integration: An application to fractional linear systems
Liu, Dayan
2013-02-04
In this article, we propose a robust method to compute the output of a fractional linear system defined through a linear fractional differential equation (FDE) with time-varying coefficients, where the input can be noisy. We firstly introduce an estimator of the fractional derivative of an unknown signal, which is defined by an integral formula obtained by calculating the fractional derivative of a truncated Jacobi polynomial series expansion. We then approximate the FDE by applying to each fractional derivative this formal algebraic integral estimator. Consequently, the fractional derivatives of the solution are applied on the used Jacobi polynomials and then we need to identify the unknown coefficients of the truncated series expansion of the solution. Modulating functions method is used to estimate these coefficients by solving a linear system issued from the approximated FDE and some initial conditions. A numerical result is given to confirm the reliability of the proposed method. © 2013 IFAC.
Stability analysis of switched linear systems defined by graphs
Athanasopoulos, Nikolaos; Lazar, Mircea
2015-01-01
We present necessary and sufficient conditions for global exponential stability for switched discrete-time linear systems, under arbitrary switching, which is constrained within a set of admissible transitions. The class of systems studied includes the family of systems under arbitrary switching, periodic systems, and systems with minimum and maximum dwell time specifications. To reach the result, we describe the set of rules that define the admissible transitions with a weighted directed gra...
Modeling and analysis of linear hyperbolic systems of balance laws
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...
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.
On the existence of tunneling bounce solutions in piecewise linear potentials
Energy Technology Data Exchange (ETDEWEB)
Dutta, Koushik; Hector, Cecelie; Konstandin, Thomas; Vaudrevange, Pascal M.; Westphal, Alexander
2012-02-15
Coleman tunneling in a general scalar potential with two non-degenerate minima is known to have an approximation in terms of a piecewise linear triangular-shaped potential with sharp 'kinks' at the place of the local minima. This approximate potential has a regime where the existence of the bounce solution needs the scalar field to 'wait' for some amount of Euclidean time at one of the 'kinks'. We discuss under which circumstances the correct bounce action can be consistently obtained as the limiting case of a regular scalar potential where 'kinks' are resolved as locally smooth 'cap' regions. (orig.)
International Nuclear Information System (INIS)
Shimizu, Yoshiaki
1991-01-01
In recent complicated nuclear systems, there are increasing demands for developing highly advanced procedures for various problems-solvings. Among them keen interests have been paid on man-machine communications to improve both safety and economy factors. Many optimization methods have been good enough to elaborate on these points. In this preliminary note, we will concern with application of linear programming (LP) for this purpose. First we will present a new superior version of the generalized PAPA method (GEPAPA) to solve LP problems. We will then examine its effectiveness when applied to derive dynamic matrix control (DMC) as the LP solution. The approach is to aim at the above goal through a quality control of process that will appear in the system. (author)
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)
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.)
Non-linear dynamics and alternating 'flip' solutions in ferrofluidic Taylor-Couette flow
Altmeyer, Sebastian
2018-04-01
This study treats with the influence of a symmetry-breaking transversal magnetic field on the nonlinear dynamics of ferrofluidic Taylor-Couette flow - flow confined between two concentric independently rotating cylinders. We detected alternating 'flip' solutions which are flow states featuring typical characteristics of slow-fast-dynamics in dynamical systems. The flip corresponds to a temporal change in the axial wavenumber and we find them to appear either as pure 2-fold axisymmetric (due to the symmetry-breaking nature of the applied transversal magnetic field) or involving non-axisymmetric, helical modes in its interim solution. The latter ones show features of typical ribbon solutions. In any case the flip solutions have a preferential first axial wavenumber which corresponds to the more stable state (slow dynamics) and second axial wavenumber, corresponding to the short appearing more unstable state (fast dynamics). However, in both cases the flip time grows exponential with increasing the magnetic field strength before the flip solutions, living on 2-tori invariant manifolds, cease to exist, with lifetime going to infinity. Further we show that ferrofluidic flow turbulence differ from the classical, ordinary (usually at high Reynolds number) turbulence. The applied magnetic field hinders the free motion of ferrofluid partials and therefore smoothen typical turbulent quantities and features so that speaking of mildly chaotic dynamics seems to be a more appropriate expression for the observed motion.
Evaluation of Security Solutions for Android Systems
Shabtai, Asaf; Mimran, Dudu; Elovici, Yuval
2015-01-01
With the increasing usage of smartphones a plethora of security solutions are being designed and developed. Many of the security solutions fail to cope with advanced attacks and are not aways properly designed for smartphone platforms. Therefore, there is a need for a methodology to evaluate their effectiveness. Since the Android operating system has the highest market share today, we decided to focus on it in this study in which we review some of the state-of-the-art security solutions for A...
Existence of infinitely many periodic solutions for second-order nonautonomous Hamiltonian systems
Directory of Open Access Journals (Sweden)
Wen Guan
2015-04-01
Full Text Available By using minimax methods and critical point theory, we obtain infinitely many periodic solutions for a second-order nonautonomous Hamiltonian systems, when the gradient of potential energy does not exceed linear growth.
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 afﬁne linear parametervarying systems (A LPV). The class of control structures includes decentralized of any order, ﬁxed 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....
Soliton solutions in a diatomic lattice system
International Nuclear Information System (INIS)
Yajima, Nobuo; Satsuma, Junkichi.
1979-04-01
A continuum limit is considered for a diatomic lattice system with a cubic nonlinearity. A long wave equation describing the interaction of acoustic and optical modes is obtained. It reduces, in certain approximations, to equations having coupled wave solutions. The solutions exhibit trapping of an optical mode by an acoustic soliton. The form of the trapped optical wave depends on the mass ratio of adjacent particles in the diatomic lattice. (author)
Forced vibration of nonlinear system with symmetrical piecewise-linear characteristics
International Nuclear Information System (INIS)
Watanabe, Takeshi
1983-01-01
It is fairly difficult to treat exactly the analysis of a vibrating system including some play because it is accompanied by a strong nonlinear phenomenon of collision. The author attempted the theoretical analysis by the exact solution using series solution and the approximate solution, treating the forced vibration of a system having some play as the forced vibration of a continuous system with nonlinear boundary condition or the colliding vibration of a continuum. In this report, the problem of such system with play is treated as a nonlinear system having the symmetrical, piecewise linear characteristics of one degree of freedom. That is, it is considered that at the time of collision due to play, the collided body causes the deformation accompanied by triangular hystersis elastically and plastically, and the spring characteristics of restitution force change piecewise by the collision. The exact solution using series solution and the approximate solution are performed, and the effectiveness of these theoretical solutions is confirmed by comparing with the solution using an analog computer. The relation between the accuracy of two analysis methods and nonlinear parameters is shown by the examples of numerical calculation. (Kako, I.)
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
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
State space and input-output linear systems
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...
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...
Directory of Open Access Journals (Sweden)
Luo Li-Qin
2016-01-01
Full Text Available In this paper, we investigate the value distribution of meromorphic solutions of homogeneous and non-homogeneous complex linear differential-difference equations, and obtain the results on the relations between the order of the solutions and the convergence exponents of the zeros, poles, a-points and small function value points of the solutions, which show the relations in the case of non-homogeneous equations are sharper than the ones in the case of homogeneous equations.
Unification of three linear models for the transient visual system
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
Punctuated equilibrium in a non-linear system of action
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
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.
INPUT-OUTPUT STRUCTURE OF LINEAR-DIFFERENTIAL ALGEBRAIC SYSTEMS
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
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...
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
Criteria for stability of linear dynamical systems with multiple delays ...
African Journals Online (AJOL)
In this study we considered a linear Dynamical system with multiple delays and find suitable conditions on the systems parameters such that for a given initial function, we can define a mapping in a carefully chosen complete metric space on which the mapping has a unique fixed point. An asymptotic stability theory for the ...
A data-acquisition system for high speed linear CCD
International Nuclear Information System (INIS)
Liu Zhiyan; Chen Xiangcai; Jiang Xiaoshan; Zhang Hongyu; Liang Zhongwang; Xiang Haisheng; Hu Jun
2010-01-01
A data-acquisition system for high speed linear CCD (Charge Coupled device) is mainly introduced. The optical fiber transmission technology is used. The data is sent to PC through USB or PCI interface. The construction of the system, the design of the PCI interface hardware, software design and the design of the control program running on host computer are also introduced. (authors)
Partial Synchronization Manifolds for Linearly Time-Delay Coupled Systems
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...
The linear sizes tolerances and fits system modernization
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.
Mean Transit Time and Mean Residence Time for Linear Diffusion–Convection–Reaction Transport System
Directory of Open Access Journals (Sweden)
Jacek Waniewski
2007-01-01
Full Text Available Characteristic times for transport processes in biological systems may be evaluated as mean transit times (MTTs (for transit states or mean residence times (MRT (for steady states. It is shown in a general framework of a (linear reaction–diffusion–convection equation that these two times are related. Analytical formulas are also derived to calculate moments of exit time distribution using solutions for a stationary state of the system.
On Energy Efficient Mobile Hydraulic Systems : with Focus on Linear Actuation
Heybroek, Kim
2017-01-01
In this dissertation, energy efficient hydraulic systems are studied. The research focuses on solutions for linear actuators in mobile applications, with emphasis on construction machines. Alongside the aspect of energy efficiency, the thesis deals with competing aspects in hydraulic system design found in the development of construction machines. Simulation models and controls for different concepts are developed, taking the whole machine into account. In line with this work, several proof o...
Damped oscillations of linear systems a mathematical introduction
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 ...
Ultra-high Frequency Linear Fiber Optic Systems
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...
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...
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
Linear dynamical quantum systems analysis, synthesis, and control
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 ...
Exact solutions to chaotic and stochastic systems
González, J. A.; Reyes, L. I.; Guerrero, L. E.
2001-03-01
We investigate functions that are exact solutions to chaotic dynamical systems. A generalization of these functions can produce truly random numbers. For the first time, we present solutions to random maps. This allows us to check, analytically, some recent results about the complexity of random dynamical systems. We confirm the result that a negative Lyapunov exponent does not imply predictability in random systems. We test the effectiveness of forecasting methods in distinguishing between chaotic and random time series. Using the explicit random functions, we can give explicit analytical formulas for the output signal in some systems with stochastic resonance. We study the influence of chaos on the stochastic resonance. We show, theoretically, the existence of a new type of solitonic stochastic resonance, where the shape of the kink is crucial. Using our models we can predict specific patterns in the output signal of stochastic resonance systems.
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
The graphics software of the Saclay linear accelerator control system
International Nuclear Information System (INIS)
Gournay, J.F.
1987-06-01
The Control system of the Saclay Linear Accelerator is based upon modern technology hardware. In the graphic software, pictures are created in exactly the same manner for all the graphic devices supported by the system. The informations used to draw a picture are stored in an array called a graphic segment. Three output primitives are used to add graphic material in a segment. Three coordinate systems are defined
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.
Chaos synchronization of a unified chaotic system via partial linearization
International Nuclear Information System (INIS)
Yu Yongguang; Li Hanxiong; Duan Jian
2009-01-01
A partial linearization method is proposed for realizing the chaos synchronization of an unified chaotic system. Through synchronizing partial state of the chaotic systems can result in the synchronization of their entire states, and the resulting controller is singularity free. The results can be easily extended to the synchronization of other similar chaotic systems. Simulation results are conducted to show the effectiveness of the method.
Hasegawa, Chihiro; Duffull, Stephen B
2018-02-01
Pharmacokinetic-pharmacodynamic systems are often expressed with nonlinear ordinary differential equations (ODEs). While there are numerous methods to solve such ODEs these methods generally rely on time-stepping solutions (e.g. Runge-Kutta) which need to be matched to the characteristics of the problem at hand. The primary aim of this study was to explore the performance of an inductive approximation which iteratively converts nonlinear ODEs to linear time-varying systems which can then be solved algebraically or numerically. The inductive approximation is applied to three examples, a simple nonlinear pharmacokinetic model with Michaelis-Menten elimination (E1), an integrated glucose-insulin model and an HIV viral load model with recursive feedback systems (E2 and E3, respectively). The secondary aim of this study was to explore the potential advantages of analytically solving linearized ODEs with two examples, again E3 with stiff differential equations and a turnover model of luteinizing hormone with a surge function (E4). The inductive linearization coupled with a matrix exponential solution provided accurate predictions for all examples with comparable solution time to the matched time-stepping solutions for nonlinear ODEs. The time-stepping solutions however did not perform well for E4, particularly when the surge was approximated by a square wave. In circumstances when either a linear ODE is particularly desirable or the uncertainty in matching the integrator to the ODE system is of potential risk, then the inductive approximation method coupled with an analytical integration method would be an appropriate alternative.
SNR Estimation in Linear Systems with Gaussian Matrices
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
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.
Experimental quantum computing to solve systems of linear equations.
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.
Should Pruning be a Pre-Processor of any Linear System?
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
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....
Energy Technology Data Exchange (ETDEWEB)
Fuentes, Carlos [Instituto Mexicano de Tecnologia del Agua, Jiutepec, Morelos (Mexico); Parlangue, Jean-Yves [Departamento de Agricultura e Ingenieria Biologica (United States); Haverkamp, Randel; Vauclin, Michael [Laboratorio de Estudio de las Transferencias en Hidrologia y Medio ambiente (France)
2001-12-01
The exact solution of the one-dimensional vertical infiltration equation is deducted, when the hydraulic diffusivity is considered constant and the hydraulic conductivity is a combination of both a linear and quadratic functions of the soil water content. This quasi-linear solution includes as particular cases, both the classical solution known as linear soil and the Knight solution. The cumulative infiltrated water as a function of time provided by the quasi-linear solution has been compared with the cumulative infiltrated water obtained from the numerical solution of the Richards equation on three different soils of contrasting hydrodynamic properties. The good agreement between the two solutions has shown that the quasi-linear solution can be used on soils where the accepted hypothesis, on hydraulic diffusivity and hydraulic conductivity, for its deduction is not satisfied. [Spanish] Se deduce la solucion exacta de la ecuacion de la infiltracion unidimensional vertical cuando la difusividad hidraulica es considerada constante y la conductividad hidraulica es una combinacion de una funcion lineal y una cuadratica del contenido volumetrico de agua. Esta solucion cuasi-lineal de la infiltracion contiene, como casos particulares, la solucion clasica conocida como suelo lineal y la solucion de Knight. La lamina infiltrada acumulada en funcion del tiempo proporcionada por la solucion cuasi-lineal se ha comparado con la lamina infiltrada proporcionada por la solucion numerica de la ecuacion de Richards en tres suelos de propiedades hidrodinamicas contrastantes. El buen acuerdo entre las laminas infiltradas ha mostrado que la solucion cuasi-lineal puede utilizarse en suelos donde la difusividad y la conductividad hidraulicas no satisfacen los supuestos de la deduccion.
DEFF Research Database (Denmark)
Fitzek, Frank; Toth, Tamas; Szabados, Áron
2014-01-01
This paper advocates the use of random linear network coding for storage in distributed clouds in order to reduce storage and traffic costs in dynamic settings, i.e. when adding and removing numerous storage devices/clouds on-the-fly and when the number of reachable clouds is limited. We introduce...... various network coding approaches that trade-off reliability, storage and traffic costs, and system complexity relying on probabilistic recoding for cloud regeneration. We compare these approaches with other approaches based on data replication and Reed-Solomon codes. A simulator has been developed...... to carry out a thorough performance evaluation of the various approaches when relying on different system settings, e.g., finite fields, and network/storage conditions, e.g., storage space used per cloud, limited network use, and limited recoding capabilities. In contrast to standard coding approaches, our...
Linear-constraint wavefront control for exoplanet coronagraphic imaging systems
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.
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
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
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 ...
Force analysis of linear induction motor for magnetic levitation system
Kuijpers, A.A.; Nemlioglu, C.; Sahin, F.; Verdel, A.J.D.; Compter, J.C.; Lomonova, E.
2010-01-01
This paper presents the analyses of thrust and normal forces of linear induction motor (LIM) segments which are implemented in a rotating ring system. To obtain magnetic levitation in a cost effective and sustainable way, decoupled control of thrust and normal forces is required. This study includes
Input design for linear dynamic systems using maxmin criteria
DEFF Research Database (Denmark)
Sadegh, Payman; Hansen, Lars H.; Madsen, Henrik
1998-01-01
This paper considers the problem of input design for maximizing the smallest eigenvalue of the information matrix for linear dynamic systems. The optimization of the smallest eigenvalue is of interest in parameter estimation and parameter change detection problems. We describe a simple cutting...
Generating Nice Linear Systems for Matrix Gaussian Elimination
Homewood, L. James
2004-01-01
In this article an augmented matrix that represents a system of linear equations is called nice if a sequence of elementary row operations that reduces the matrix to row-echelon form, through matrix Gaussian elimination, does so by restricting all entries to integers in every step. Many instructors wish to use the example of matrix Gaussian…
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%.
Robust self-triggered MPC for constrained linear systems
Brunner, F.D.; Heemels, W.P.M.H.; Allgöwer, F.
2014-01-01
In this paper we propose a robust self-triggered model predictive control algorithm for linear systems with additive bounded disturbances and hard constraints on the inputs and state. In self-triggered control, at every sampling instant the time until the next sampling instant is computed online
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...
Relative controllability and null controllability of linear delay systems ...
African Journals Online (AJOL)
Necessary and sufficient conditions are established for the relative, absolute controllability and null controllability of the generalized linear delay system and its discrete prototype. The paper presents illuminating examples on previous controllability results by Manitius and Olbrot [7] and carries over the results of Onwuatu [8] ...
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)
Optimal Linear Responses for Markov Chains and Stochastically Perturbed Dynamical Systems
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.
Information Systems Solutions for Environmental Sustainability
DEFF Research Database (Denmark)
Gholami, Roya; Watson, Richard T.; Hasan, Helen
2016-01-01
We contend that too few information systems (IS) academics engage in impactful research that offers solutions to global warming despite the fact that climate change is one of the most critical challenges facing this generation. Climate change is a major threat to global sustainability in the 21st...... themselves in creating solutions for environmental problems. Moreover, information is a perquisite for assessing the state of the environment and making appropriate decisions to ameliorate identified problems. Indeed, the IS scholarly community needs to help create a sustainable society. While...
Periodic solutions of dissipative systems revisited
Directory of Open Access Journals (Sweden)
Górniewicz Lech
2006-01-01
Full Text Available We reprove in an extremely simple way the classical theorem that time periodic dissipative systems imply the existence of harmonic periodic solutions, in the case of uniqueness. We will also show that, in the lack of uniqueness, the existence of harmonics is implied by uniform dissipativity. The localization of starting points and multiplicity of periodic solutions will be established, under suitable additional assumptions, as well. The arguments are based on the application of various asymptotic fixed point theorems of the Lefschetz and Nielsen type.
Periodic solutions of dissipative systems revisited
Directory of Open Access Journals (Sweden)
Lech Górniewicz
2006-05-01
Full Text Available We reprove in an extremely simple way the classical theorem that time periodic dissipative systems imply the existence of harmonic periodic solutions, in the case of uniqueness. We will also show that, in the lack of uniqueness, the existence of harmonics is implied by uniform dissipativity. The localization of starting points and multiplicity of periodic solutions will be established, under suitable additional assumptions, as well. The arguments are based on the application of various asymptotic fixed point theorems of the Lefschetz and Nielsen type.
Liu, Tao; Huang, Jie
2017-04-17
This paper presents a discrete-time recurrent neural network approach to solving systems of linear equations with two features. First, the system of linear equations may not have a unique solution. Second, the system matrix is not known precisely, but a sequence of matrices that converges to the unknown system matrix exponentially is known. The problem is motivated from solving the output regulation problem for linear systems. Thus, an application of our main result leads to an online solution to the output regulation problem for linear systems.
Directory of Open Access Journals (Sweden)
Chun-Fu Chen
2014-03-01
Full Text Available Linear analytical study on the mechanical sensitivity in large deflection of unsymmetrically layered and laterally loaded piezoelectric plate under pretension is conducted. von Karman plate theory for large deflection is utilized but extended to the case of an unsymmetrically layered plate embedded with a piezoelectric layer. The governing equations thus obtained are simplified by omitting the arising nonlinear terms, yielding a Bessel or modified Bessel equation for the lateral slope. Depending on the relative magnitude of the piezoelectric effect, for both cases, analytical solutions of various geometrical responses are developed and formulated via Bessel and modified Bessel functions. The associated ultimate radial stresses are further derived following lamina constitutive law to evaluate the mechanical sensitivity of the considered plate. For a nearly monolithic plate under a very low applied voltage, the results are in good agreement with those for a single-layered case due to pure mechanical load available in literature, and thus the present approach is checked. For a two-layered unsymmetric plate made of typical silicon-based materials, a sound piezoelectric effect is illustrated particularly in a low pretension condition.
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
Flexible solution for interoperable cloud healthcare systems.
Vida, Mihaela Marcella; Lupşe, Oana Sorina; Stoicu-Tivadar, Lăcrămioara; Bernad, Elena
2012-01-01
It is extremely important for the healthcare domain to have a standardized communication because will improve the quality of information and in the end the resulting benefits will improve the quality of patients' life. The standards proposed to be used are: HL7 CDA and CCD. For a better access to the medical data a solution based on cloud computing (CC) is investigated. CC is a technology that supports flexibility, seamless care, and reduced costs of the medical act. To ensure interoperability between healthcare information systems a solution creating a Web Custom Control is presented. The control shows the database tables and fields used to configure the two standards. This control will facilitate the work of the medical staff and hospital administrators, because they can configure the local system easily and prepare it for communication with other systems. The resulted information will have a higher quality and will provide knowledge that will support better patient management and diagnosis.
Monitoring and control system of the Saclay electron linear accelerator
International Nuclear Information System (INIS)
Lafontaine, Antoine
1974-01-01
A description is given of the automatic monitoring and control system of the 60MeV electron linear accelerator of the Centre d'Etudes Nucleaires de Saclay. The paper is mostly concerned with the programmation of the system. However, in a real time device, there is a very close association between computer and electronics, the latter are therefore described in details and make up most of the paper. [fr
A new timing system for the Stanford Linear Collider
International Nuclear Information System (INIS)
Paffrath, L.; Bernstein, D.; Kang, H.; Koontz, R.; Leger, G.; Pierce, W.; Ross, M.; Wilmunder, A.
1985-01-01
In order to be able to meet the goals of the Stanford Linear Collider, a much more precise timing system had to be implemented. This paper describes the specification and design of this system, and the results obtained from its use on 1/3 of the SLAC linac. The functions of various elements are described, and a programmable delay unit (PDU) is described in detail
Hyperchaotic encryption based on multi-scroll piecewise linear Systems
Czech Academy of Sciences Publication Activity Database
García-Martínez, M.; Ontanon-García, L.J.; Campos-Cantón, E.; Čelikovský, Sergej
2015-01-01
Roč. 270, č. 1 (2015), s. 413-424 ISSN 0096-3003 R&D Projects: GA ČR GA13-20433S Institutional support: RVO:67985556 Keywords : Hyperchaotic encryption * Piecewise linear systems * Stream cipher * Pseudo-random bit generator * Chaos theory * Multi-scrollattractors Subject RIV: BC - Control Systems Theory Impact factor: 1.345, year: 2015 http://library.utia.cas.cz/separaty/2015/TR/celikovsky-0446895.pdf
Global Linear Representations of Nonlinear Systems and the Adjoint Map
Banks, S.P.
1988-01-01
In this paper we shall study the global linearization of nonlinear systems on a manifold by two methods. The first consists of an expansion of the vector field in the space of square integrable vector fields. In the second method we use the adjoint representation of the Lie algebra vector fields to obtain an infinite-dimensional matrix representation of the system. A connection between the two approaches will be developed.
Comments on new iterative methods for solving linear systems
Directory of Open Access Journals (Sweden)
Wang Ke
2017-06-01
Full Text Available Some new iterative methods were presented by Du, Zheng and Wang for solving linear systems in [3], where it is shown that the new methods, comparing to the classical Jacobi or Gauss-Seidel method, can be applied to more systems and have faster convergence. This note shows that their methods are suitable for more matrices than positive matrices which the authors suggested through further analysis and numerical examples.
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.
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
Design and performance of the Stanford Linear Collider Control System
International Nuclear Information System (INIS)
Melen, R.E.
1984-10-01
The success of the Stanford Linear Collider (SLC) will be dependent upon the implementation of a very large advanced computer-based instrumentation and control system. This paper describes the architectural design of this system as well as a critique of its performance. This critique is based on experience obtained from its use in the control and monitoring of 1/3 of the SLAC linac and in support of an expensive experimental machine physics experimental program. 11 references, 3 figures
Geon-type solutions of the non-linear Heisenberg-Klein-Gordon equation
International Nuclear Information System (INIS)
Mielke, E.W.; Scherzer, R.
1980-10-01
As a model for a ''unitary'' field theory of extended particles we consider the non-linear Klein-Gordon equation - associated with a ''squared'' Heisenberg-Pauli-Weyl non-linear spinor equation - coupled to strong gravity. Using a stationary spherical ansatz for the complex scalar field as well as for the background metric generated via Einstein's field equation, we are able to study the effects of the scalar self-interaction as well as of the classical tensor forces. By numerical integration we obtain a continuous spectrum of localized, gravitational solitons resembling the geons previously constructed for the Einstein-Maxwell system by Wheeler. A self-generated curvature potential originating from the curved background partially confines the Schroedinger type wave functions within the ''scalar geon''. For zero angular momentum states and normalized scalar charge the spectrum for the total gravitational energy of these solitons exhibits a branching with respect to the number of nodes appearing in the radial part of the scalar field. Preliminary studies for higher values of the corresponding ''principal quantum number'' reveal that a kind of fine splitting of the energy levels occurs, which may indicate a rich, particle-like structure of these ''quantized geons''. (author)
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.
Projective-Dual Method for Solving Systems of Linear Equations with Nonnegative Variables
Ganin, B. V.; Golikov, A. I.; Evtushenko, Yu. G.
2018-02-01
In order to solve an underdetermined system of linear equations with nonnegative variables, the projection of a given point onto its solutions set is sought. The dual of this problem—the problem of unconstrained maximization of a piecewise-quadratic function—is solved by Newton's method. The problem of unconstrained optimization dual of the regularized problem of finding the projection onto the solution set of the system is considered. A connection of duality theory and Newton's method with some known algorithms of projecting onto a standard simplex is shown. On the example of taking into account the specifics of the constraints of the transport linear programming problem, the possibility to increase the efficiency of calculating the generalized Hessian matrix is demonstrated. Some examples of numerical calculations using MATLAB are presented.
Self-Tuning Control of Linear Systems Followed by Deadzones
Directory of Open Access Journals (Sweden)
K. Kazlauskas
2014-02-01
Full Text Available The aim of the present paper is to increase the efficiency of self-tuning generalized minimum variance (GMV control of linear time-invariant (LTI systems followed by deadzone nonlinearities. An approach, based on reordering of observations to be processed for the reconstruction of an unknown internal signal that acts between LTI system and a static nonlinear block of the closed-loop Wiener system, has been developed. The results of GMV self-tuning control of the second order LTI system with an ordinary deadzone are given.
Algorithmic Approach to Abstracting Linear Systems by Timed Automata
DEFF Research Database (Denmark)
Sloth, Christoffer; Wisniewski, Rafael
2011-01-01
This paper proposes an LMI-based algorithm for abstracting dynamical systems by timed automata, which enables automatic formal verification of linear systems. The proposed abstraction is based on partitioning the state space of the system using positive invariant sets, generated by Lyapunov...... functions. This partitioning ensures that the vector field of the dynamical system is transversal to all facets of the cells, which induces some desirable properties of the abstraction. The algorithm is based on identifying intersections of level sets of quadratic Lyapunov functions, and determining...
Directory of Open Access Journals (Sweden)
Sukhpreet Kaur Sidhu
2014-01-01
Full Text Available The drawbacks of the existing methods to obtain the fuzzy optimal solution of such linear programming problems, in which coefficients of the constraints are represented by real numbers and all the other parameters as well as variables are represented by symmetric trapezoidal fuzzy numbers, are pointed out, and to resolve these drawbacks, a new method (named as Mehar method is proposed for the same linear programming problems. Also, with the help of proposed Mehar method, a new method, much easy as compared to the existing methods, is proposed to deal with the sensitivity analysis of the same type of linear programming problems.
International Nuclear Information System (INIS)
Fernandes, L.; Friedlander, A.; Guedes, M.; Judice, J.
2001-01-01
This paper addresses a General Linear Complementarity Problem (GLCP) that has found applications in global optimization. It is shown that a solution of the GLCP can be computed by finding a stationary point of a differentiable function over a set defined by simple bounds on the variables. The application of this result to the solution of bilinear programs and LCPs is discussed. Some computational evidence of its usefulness is included in the last part of the paper
International Nuclear Information System (INIS)
Anastassi, Z. A.; Simos, T. E.
2010-01-01
We develop a new family of explicit symmetric linear multistep methods for the efficient numerical solution of the Schroedinger equation and related problems with oscillatory solution. The new methods are trigonometrically fitted and have improved intervals of periodicity as compared to the corresponding classical method with constant coefficients and other methods from the literature. We also apply the methods along with other known methods to real periodic problems, in order to measure their efficiency.
Fundamentals of linear systems for physical scientists and engineers
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...
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)
Mixed problems for linear symmetric hyperbolic systems with characteristic boundary conditions
International Nuclear Information System (INIS)
Secchi, P.
1994-01-01
We consider the initial-boundary value problem for symmetric hyperbolic systems with characteristic boundary of constant multiplicity. In the linear case we give some results about the existence of regular solutions in suitable functions spaces which take in account the loss of regularity in the normal direction to the characteristic boundary. We also consider the equations of ideal magneto-hydrodynamics under perfectly conducting wall boundary conditions and give some results about the solvability of such mixed problem. (author). 16 refs
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.
Challenging problems and solutions in intelligent systems
Grzegorzewski, Przemysław; Kacprzyk, Janusz; Owsiński, Jan; Penczek, Wojciech; Zadrożny, Sławomir
2016-01-01
This volume presents recent research, challenging problems and solutions in Intelligent Systems– covering the following disciplines: artificial and computational intelligence, fuzzy logic and other non-classic logics, intelligent database systems, information retrieval, information fusion, intelligent search (engines), data mining, cluster analysis, unsupervised learning, machine learning, intelligent data analysis, (group) decision support systems, intelligent agents and multi-agent systems, knowledge-based systems, imprecision and uncertainty handling, electronic commerce, distributed systems, etc. The book defines a common ground for sometimes seemingly disparate problems and addresses them by using the paradigm of broadly perceived intelligent systems. It presents a broad panorama of a multitude of theoretical and practical problems which have been successfully dealt with using the paradigm of intelligent computing.
On preconditioner updates for sequences of saddle-point linear systems
Directory of Open Access Journals (Sweden)
Simone Valentina De
2018-02-01
Full Text Available Updating preconditioners for the solution of sequences of large and sparse saddle- point linear systems via Krylov methods has received increasing attention in the last few years, because it allows to reduce the cost of preconditioning while keeping the efficiency of the overall solution process. This paper provides a short survey of the two approaches proposed in the literature for this problem: updating the factors of a preconditioner available in a block LDLT form, and updating a preconditioner via a limited-memory technique inspired by quasi-Newton methods.
Energy-Water System Solutions | Energy Analysis | NREL
System Solutions Energy-Water System Solutions NREL has been a pioneer in the development of energy -water system solutions that explicitly address and optimize energy-water tradeoffs. NREL has evaluated energy-water system solutions for Department of Defense bases, islands, communities recovering from
Essential uncontrollability of discrete linear, time-invariant, dynamical systems
Cliff, E. M.
1975-01-01
The concept of a 'best approximating m-dimensional subspace' for a given set of vectors in n-dimensional whole space is introduced. Such a subspace is easily described in terms of the eigenvectors of an associated Gram matrix. This technique is used to approximate an achievable set for a discrete linear time-invariant dynamical system. This approximation characterizes the part of the state space that may be reached using modest levels of control. If the achievable set can be closely approximated by a proper subspace of the whole space then the system is 'essentially uncontrollable'. The notion finds application in studies of failure-tolerant systems, and in decoupling.
The new control system of the Saclay linear accelerator
International Nuclear Information System (INIS)
Gournay, J.F.; Gourcy, G.; Garreau, F.; Giraud, A.; Rouault, J.
1985-05-01
A new control system for the Safety Linear Accelerator is now being designed. The computer control architecture is based on 3 dedicated VME crates with MC68000 micro-processors: one crate with a disk-based operating system will run the high level application programs and the data base management facilities, another one will manage the man-machine communications and the third one will interface the system to the linac equipments. Communications between the VME microcomputers will be done through 16 bit parallel links. The software is modular and organized in specific layers, the data base is fully distributed. About 90% of the code is written in Fortran
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.
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....
Coherent versus Measurement Feedback: Linear Systems Theory for Quantum Information
Directory of Open Access Journals (Sweden)
Naoki Yamamoto
2014-11-01
Full Text Available To control a quantum system via feedback, we generally have two options in choosing a control scheme. One is the coherent feedback, which feeds the output field of the system, through a fully quantum device, back to manipulate the system without involving any measurement process. The other one is measurement-based feedback, which measures the output field and performs a real-time manipulation on the system based on the measurement results. Both schemes have advantages and disadvantages, depending on the system and the control goal; hence, their comparison in several situations is important. This paper considers a general open linear quantum system with the following specific control goals: backaction evasion, generation of a quantum nondemolished variable, and generation of a decoherence-free subsystem, all of which have important roles in quantum information science. Some no-go theorems are proven, clarifying that those goals cannot be achieved by any measurement-based feedback control. On the other hand, it is shown that, for each control goal there exists a coherent feedback controller accomplishing the task. The key idea to obtain all the results is system theoretic characterizations of the above three notions in terms of controllability and observability properties or transfer functions of linear systems, which are consistent with their standard definitions.
Sparse symmetric preconditioners for dense linear systems in electromagnetism
Carpentieri, Bruno; Duff, Iain S.; Giraud, Luc; Monga Made, M. Magolu
2004-01-01
We consider symmetric preconditioning strategies for the iterative solution of dense complex symmetric non-Hermitian systems arising in computational electromagnetics. In particular, we report on the numerical behaviour of the classical incomplete Cholesky factorization as well as some of its recent
A Combined MPI-CUDA Parallel Solution of Linear and Nonlinear Poisson-Boltzmann Equation
Directory of Open Access Journals (Sweden)
José Colmenares
2014-01-01
Full Text Available The Poisson-Boltzmann equation models the electrostatic potential generated by fixed charges on a polarizable solute immersed in an ionic solution. This approach is often used in computational structural biology to estimate the electrostatic energetic component of the assembly of molecular biological systems. In the last decades, the amount of data concerning proteins and other biological macromolecules has remarkably increased. To fruitfully exploit these data, a huge computational power is needed as well as software tools capable of exploiting it. It is therefore necessary to move towards high performance computing and to develop proper parallel implementations of already existing and of novel algorithms. Nowadays, workstations can provide an amazing computational power: up to 10 TFLOPS on a single machine equipped with multiple CPUs and accelerators such as Intel Xeon Phi or GPU devices. The actual obstacle to the full exploitation of modern heterogeneous resources is efficient parallel coding and porting of software on such architectures. In this paper, we propose the implementation of a full Poisson-Boltzmann solver based on a finite-difference scheme using different and combined parallel schemes and in particular a mixed MPI-CUDA implementation. Results show great speedups when using the two schemes, achieving an 18.9x speedup using three GPUs.
Ryo, Ikehata
Uniform energy and L2 decay of solutions for linear wave equations with localized dissipation will be given. In order to derive the L2-decay property of the solution, a useful device whose idea comes from Ikehata-Matsuyama (Sci. Math. Japon. 55 (2002) 33) is used. In fact, we shall show that the L2-norm and the total energy of solutions, respectively, decay like O(1/ t) and O(1/ t2) as t→+∞ for a kind of the weighted initial data.
Optimal approximation of linear systems by artificial immune response
Institute of Scientific and Technical Information of China (English)
无
2006-01-01
This paper puts forward a novel artificial immune response algorithm for optimal approximation of linear systems. A quaternion model of artificial immune response is proposed for engineering computing. The model abstracts four elements, namely, antigen, antibody, reaction rules among antibodies, and driving algorithm describing how the rules are applied to antibodies, to simulate the process of immune response. Some reaction rules including clonal selection rules, immunological memory rules and immune regulation rules are introduced. Using the theorem of Markov chain, it is proofed that the new model is convergent. The experimental study on the optimal approximation of a stable linear system and an unstable one show that the approximate models searched by the new model have better performance indices than those obtained by some existing algorithms including the differential evolution algorithm and the multi-agent genetic algorithm.
Large linear magnetoresistivity in strongly inhomogeneous planar and layered systems
International Nuclear Information System (INIS)
Bulgadaev, S.A.; Kusmartsev, F.V.
2005-01-01
Explicit expressions for magnetoresistance R of planar and layered strongly inhomogeneous two-phase systems are obtained, using exact dual transformation, connecting effective conductivities of in-plane isotropic two-phase systems with and without magnetic field. These expressions allow to describe the magnetoresistance of various inhomogeneous media at arbitrary concentrations x and magnetic fields H. All expressions show large linear magnetoresistance effect with different dependencies on the phase concentrations. The corresponding plots of the x- and H-dependencies of R(x,H) are represented for various values, respectively, of magnetic field and concentrations at some values of inhomogeneity parameter. The obtained results show a remarkable similarity with the existing experimental data on linear magnetoresistance in silver chalcogenides Ag 2+δ Se. A possible physical explanation of this similarity is proposed. It is shown that the random, stripe type, structures of inhomogeneities are the most suitable for a fabrication of magnetic sensors and a storage of information at room temperatures
International Nuclear Information System (INIS)
Phan Thanh An; Phan Le Na; Ngo Quoc Chung
2004-05-01
We describe a practical implementation for finding parametric domain for asymptotic stability with probability one of zero solution of linear Ito stochastic differential equations based on Korenevskij and Mitropolskij's sufficient condition and our sufficient conditions. Numerical results show that all of these sufficient conditions are crucial in the implementation. (author)
Directory of Open Access Journals (Sweden)
Mehmet Tarik Atay
2013-01-01
Full Text Available The Variational Iteration Method (VIM and Modified Variational Iteration Method (MVIM are used to find solutions of systems of stiff ordinary differential equations for both linear and nonlinear problems. Some examples are given to illustrate the accuracy and effectiveness of these methods. We compare our results with exact results. In some studies related to stiff ordinary differential equations, problems were solved by Adomian Decomposition Method and VIM and Homotopy Perturbation Method. Comparisons with exact solutions reveal that the Variational Iteration Method (VIM and the Modified Variational Iteration Method (MVIM are easier to implement. In fact, these methods are promising methods for various systems of linear and nonlinear stiff ordinary differential equations. Furthermore, VIM, or in some cases MVIM, is giving exact solutions in linear cases and very satisfactory solutions when compared to exact solutions for nonlinear cases depending on the stiffness ratio of the stiff system to be solved.
Linear Quantum Systems: Non-Classical States and Robust Stability
2016-06-29
modulation and entanglement in a compound gradient echo memory, Physical Review A 93(2) 023809 2016. We present a theoretical model for a Kerr...Carvalho, M. Hedges and M R James, Analysis of the operation of gradient echo memories using a quantum input-output model, New Journal of Physics , 15...new structured uncertainty methods that ensure robust stability of quantum systems based on nominal linear models, and (v) physical realizability
Dynamic logic architecture based on piecewise-linear systems
International Nuclear Information System (INIS)
Peng Haipeng; Liu Fei; Li Lixiang; Yang Yixian; Wang Xue
2010-01-01
This Letter explores piecewise-linear systems to construct dynamic logic architecture. The proposed schemes can discriminate the two input signals and obtain 16 kinds of logic operations by different combinations of parameters and conditions for determining the output. Each logic cell performs more flexibly, that makes it possible to achieve complex logic operations more simply and construct computing architecture with less logic cells. We also analyze the various performances of our schemes under different conditions and the characteristics of these schemes.
CHEBYSHEV ACCELERATION TECHNIQUE FOR SOLVING FUZZY LINEAR SYSTEM
Directory of Open Access Journals (Sweden)
S.H. Nasseri
2011-07-01
Full Text Available In this paper, Chebyshev acceleration technique is used to solve the fuzzy linear system (FLS. This method is discussed in details and followed by summary of some other acceleration techniques. Moreover, we show that in some situations that the methods such as Jacobi, Gauss-Sidel, SOR and conjugate gradient is divergent, our proposed method is applicable and the acquired results are illustrated by some numerical examples.
CHEBYSHEV ACCELERATION TECHNIQUE FOR SOLVING FUZZY LINEAR SYSTEM
Directory of Open Access Journals (Sweden)
S.H. Nasseri
2009-10-01
Full Text Available In this paper, Chebyshev acceleration technique is used to solve the fuzzy linear system (FLS. This method is discussed in details and followed by summary of some other acceleration techniques. Moreover, we show that in some situations that the methods such as Jacobi, Gauss-Sidel, SOR and conjugate gradient is divergent, our proposed method is applicable and the acquired results are illustrated by some numerical examples.
Efficient Preconditioning of Sequences of Nonsymmetric Linear Systems
Czech Academy of Sciences Publication Activity Database
Duintjer Tebbens, Jurjen; Tůma, Miroslav
2007-01-01
Roč. 29, č. 5 (2007), s. 1918-1941 ISSN 1064-8275 R&D Projects: GA AV ČR 1ET400300415; GA AV ČR KJB100300703 Institutional research plan: CEZ:AV0Z10300504 Keywords : preconditioned iterative methods * sparse matrices * sequences of linear algebraic systems * incomplete factorizations * factorization updates * Gauss–Jordan transformations * minimum spanning tree Subject RIV: BA - General Mathematics Impact factor: 1.784, year: 2007
AZTEC: A parallel iterative package for the solving linear systems
Energy Technology Data Exchange (ETDEWEB)
Hutchinson, S.A.; Shadid, J.N.; Tuminaro, R.S. [Sandia National Labs., Albuquerque, NM (United States)
1996-12-31
We describe a parallel linear system package, AZTEC. The package incorporates a number of parallel iterative methods (e.g. GMRES, biCGSTAB, CGS, TFQMR) and preconditioners (e.g. Jacobi, Gauss-Seidel, polynomial, domain decomposition with LU or ILU within subdomains). Additionally, AZTEC allows for the reuse of previous preconditioning factorizations within Newton schemes for nonlinear methods. Currently, a number of different users are using this package to solve a variety of PDE applications.
Thach, Trung Thanh; Shin, Donghyuk; Han, Seungsu; Lee, Sangho
2016-04-01
The conformational flexibility of linkage-specific polyubiquitin chains enables ubiquitylated proteins and their receptors to be involved in a variety of cellular processes. Linear or Met1-linked polyubiquitin chains, associated with nondegradational cellular signalling pathways, have been known to adopt multiple conformations from compact to extended conformations. However, the extent of such conformational flexibility remains open. Here, the crystal structure of linear Ub2 was determined in a more compact conformation than that of the previously known structure (PDB entry 3axc). The two structures differ significantly from each other, as shown by an r.m.s.d. between C(α) atoms of 3.1 Å. The compactness of the linear Ub2 structure in comparison with PDB entry 3axc is supported by smaller values of the radius of gyration (Rg; 18 versus 18.9 Å) and the maximum interatomic distance (Dmax; 55.5 versus 57.8 Å). Extra intramolecular hydrogen bonds formed among polar residues between the distal and proximal ubiquitin moieties seem to contribute to stabilization of the compact conformation of linear Ub2. An ensemble of three semi-extended and extended conformations of linear Ub2 was also observed by small-angle X-ray scattering (SAXS) analysis in solution. In addition, the conformational heterogeneity in linear polyubiquitin chains is clearly manifested by SAXS analyses of linear Ub3 and Ub4: at least three distinct solution conformations are observed in each chain, with the linear Ub3 conformations being compact. The results expand the extent of conformational space of linear polyubiquitin chains and suggest that changes in the conformational ensemble may be pivotal in mediating multiple signalling pathways.
Feedback Linearization Controller for a Wind Energy Power System
Directory of Open Access Journals (Sweden)
Muthana Alrifai
2016-09-01
Full Text Available This paper deals with the control of a doubly-fed induction generator (DFIG-based variable speed wind turbine power system. A system of eight ordinary differential equations is used to model the wind energy conversion system. The generator has a wound rotor type with back-to-back three-phase power converter bridges between its rotor and the grid; it is modeled using the direct-quadrature rotating reference frame with aligned stator flux. An input-state feedback linearization controller is proposed for the wind energy power system. The controller guarantees that the states of the system track the desired states. Simulation results are presented to validate the proposed control scheme. Moreover, further simulation results are shown to investigate the robustness of the proposed control scheme to changes in some of the parameters of the system.
International Nuclear Information System (INIS)
Theodorakis, Stavros
2003-01-01
We emulate the cubic term Ψ 3 in the nonlinear Schroedinger equation by a piecewise linear term, thus reducing the problem to a set of uncoupled linear inhomogeneous differential equations. The resulting analytic expressions constitute an excellent approximation to the exact solutions, as is explicitly shown in the case of the kink, the vortex, and a δ function trap. Such a piecewise linear emulation can be used for any differential equation where the only nonlinearity is a Ψ 3 one. In particular, it can be used for the nonlinear Schroedinger equation in the presence of harmonic traps, giving analytic Bose-Einstein condensate solutions that reproduce very accurately the numerically calculated ones in one, two, and three dimensions
International Nuclear Information System (INIS)
Goncalves, Glenio A.; Bodmann, Bardo; Bogado, Sergio; Vilhena, Marco T.
2008-01-01
Analytical solutions for neutron transport in cylindrical geometry is available for isotropic problems, but to the best of our knowledge for anisotropic problems are not available, yet. In this work, an analytical solution for the neutron transport equation in an infinite cylinder assuming anisotropic scattering is reported. Here we specialize the solution, without loss of generality, for the linearly anisotropic problem using the combined decomposition and HTS N methods. The key feature of this method consists in the application of the decomposition method to the anisotropic problem by virtue of the fact that the inverse of the operator associated to isotropic problem is well know and determined by the HTS N approach. So far, following the idea of the decomposition method, we apply this operator to the integral term, assuming that the angular flux appearing in the integrand is considered to be equal to the HTS N solution interpolated by polynomial considering only even powers. This leads to the first approximation for an anisotropic solution. Proceeding further, we replace this solution for the angular flux in the integral and apply again the inverse operator for the isotropic problem in the integral term and obtain a new approximation for the angular flux. This iterative procedure yields a closed form solution for the angular flux. This methodology can be generalized, in a straightforward manner, for transport problems with any degree of anisotropy. For the sake of illustration, we report numerical simulations for linearly anisotropic transport problems. (author)
Stochastic Reformulations of Linear Systems: Algorithms and Convergence Theory
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.
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
Stochastic Reformulations of Linear Systems: Algorithms and Convergence Theory
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.
International Nuclear Information System (INIS)
Pilipchuk, L. A.; Pilipchuk, A. S.
2015-01-01
In this paper we propose the theory of decomposition, methods, technologies, applications and implementation in Wol-fram Mathematica for the constructing the solutions of the sparse linear systems. One of the applications is the Sensor Location Problem for the symmetric graph in the case when split ratios of some arc flows can be zeros. The objective of that application is to minimize the number of sensors that are assigned to the nodes. We obtain a sparse system of linear algebraic equations and research its matrix rank. Sparse systems of these types appear in generalized network flow programming problems in the form of restrictions and can be characterized as systems with a large sparse sub-matrix representing the embedded network structure
Fast computation of the Maslov index for hyperbolic linear systems with periodic coefficients
International Nuclear Information System (INIS)
Chardard, F; Dias, F; Bridges, T J
2006-01-01
The Maslov index is a topological property of periodic orbits of finite-dimensional Hamiltonian systems that is widely used in semiclassical quantization, quantum chaology, stability of waves and classical mechanics. The Maslov index is determined from the analysis of a linear Hamiltonian system with periodic coefficients. In this paper, a numerical scheme is devised to compute the Maslov index for hyperbolic linear systems when the phase space has a low dimension. The idea is to compute on the exterior algebra of the ambient vector space, where the Lagrangian subspace representing the unstable subspace is reduced to a line. When the exterior algebra is projectified the Lagrangian subspace always forms a closed loop. The idea is illustrated by application to Hamiltonian systems on a phase space of dimension 4. The theory is used to compute the Maslov index for the spectral problem associated with periodic solutions of the fifth-order Korteweg de Vries equation