Integrated vehicle dynamics control using State Dependent Riccati Equations

NARCIS (Netherlands)

Bonsen, B.; Mansvelders, R.; Vermeer, E.

2010-01-01

In this paper we discuss a State Dependent Riccati Equations (SDRE) solution for Integrated Vehicle Dynamics Control (IVDC). The SDRE approach is a nonlinear variant of the well known Linear Quadratic Regulator (LQR) and implements a quadratic cost function optimization. A modified version of this

Error Estimates for Approximate Solutions of the Riccati Equation with Real or Complex Potentials

Science.gov (United States)

Finster, Felix; Smoller, Joel

2010-09-01

A method is presented for obtaining rigorous error estimates for approximate solutions of the Riccati equation, with real or complex potentials. Our main tool is to derive invariant region estimates for complex solutions of the Riccati equation. We explain the general strategy for applying these estimates and illustrate the method in typical examples, where the approximate solutions are obtained by gluing together WKB and Airy solutions of corresponding one-dimensional Schrödinger equations. Our method is motivated by, and has applications to, the analysis of linear wave equations in the geometry of a rotating black hole.

Convergence of the standard RLS method and UDUT factorisation of covariance matrix for solving the algebraic Riccati equation of the DLQR via heuristic approximate dynamic programming

Science.gov (United States)

Moraes Rêgo, Patrícia Helena; Viana da Fonseca Neto, João; Ferreira, Ernesto M.

2015-08-01

The main focus of this article is to present a proposal to solve, via UDUT factorisation, the convergence and numerical stability problems that are related to the covariance matrix ill-conditioning of the recursive least squares (RLS) approach for online approximations of the algebraic Riccati equation (ARE) solution associated with the discrete linear quadratic regulator (DLQR) problem formulated in the actor-critic reinforcement learning and approximate dynamic programming context. The parameterisations of the Bellman equation, utility function and dynamic system as well as the algebra of Kronecker product assemble a framework for the solution of the DLQR problem. The condition number and the positivity parameter of the covariance matrix are associated with statistical metrics for evaluating the approximation performance of the ARE solution via RLS-based estimators. The performance of RLS approximators is also evaluated in terms of consistence and polarisation when associated with reinforcement learning methods. The used methodology contemplates realisations of online designs for DLQR controllers that is evaluated in a multivariable dynamic system model.

Solving Algebraic Riccati Equation Real Time for Integrated Vehicle Dynamics Control

NARCIS (Netherlands)

Kunnappillil Madhusudhanan, A; Corno, M.; Bonsen, B.; Holweg, E.

2012-01-01

In this paper we present a comparison study of different computational methods to implement State Dependent Riccati Equation (SDRE) based control in real time for a vehicle dynamics control application. Vehicles are mechatronic systems with nonlinear dynamics. One of the promising nonlinear control

PEAK COVARIANCE STABILITY OF A RANDOM RICCATI EQUATION ARISING FROM KALMAN FILTERING WITH OBSERVATION LOSSES

Institute of Scientific and Technical Information of China (English)

Li XIE; Lihua XIE

2007-01-01

We consider the stability of a random Riccati equation with a Markovian binary jump coefficient. More specifically, we are concerned with the boundedness of the solution of a random Riccati difference equation arising from Kalman filtering with measurement losses. A sufficient condition for the peak covariance stability is obtained which has a simpler form and is shown to be less conservative in some cases than a very recent result in existing literature. Furthermore, we show that a known sufficient condition is also necessary when the observability index equals one.

On Traveling Waves in Lattices: The Case of Riccati Lattices

Science.gov (United States)

Dimitrova, Zlatinka

2012-09-01

The method of simplest equation is applied for analysis of a class of lattices described by differential-difference equations that admit traveling-wave solutions constructed on the basis of the solution of the Riccati equation. We denote such lattices as Riccati lattices. We search for Riccati lattices within two classes of lattices: generalized Lotka-Volterra lattices and generalized Holling lattices. We show that from the class of generalized Lotka-Volterra lattices only the Wadati lattice belongs to the class of Riccati lattices. Opposite to this many lattices from the Holling class are Riccati lattices. We construct exact traveling wave solutions on the basis of the solution of Riccati equation for three members of the class of generalized Holling lattices.

Algorithm for solving polynomial algebraic Riccati equations and its application

Czech Academy of Sciences Publication Activity Database

Augusta, Petr; Augustová, Petra

2012-01-01

Roč. 1, č. 4 (2012), s. 237-242 ISSN 2223-7038 R&D Projects: GA ČR GPP103/12/P494 Institutional support: RVO:67985556 Keywords : Numerical algorithms * algebraic Riccati equation * spatially distributed systems * optimal control Subject RIV: BC - Control Systems Theory http://lib.physcon.ru/doc?id=8b4876d6a57d

On the Solution of the Rational Matrix Equation

Directory of Open Access Journals (Sweden)

Faßbender Heike

2007-01-01

Full Text Available We study numerical methods for finding the maximal symmetric positive definite solution of the nonlinear matrix equation , where is symmetric positive definite and is nonsingular. Such equations arise for instance in the analysis of stationary Gaussian reciprocal processes over a finite interval. Its unique largest positive definite solution coincides with the unique positive definite solution of a related discrete-time algebraic Riccati equation (DARE. We discuss how to use the butterfly algorithm to solve the DARE. This approach is compared to several fixed-point and doubling-type iterative methods suggested in the literature.

A new Riccati equation rational expansion method and its application to (2 + 1)-dimensional Burgers equation

International Nuclear Information System (INIS)

Wang Qi; Chen Yong; Zhang Hongqing

2005-01-01

In this paper, we present a new Riccati equation rational expansion method to uniformly construct a series of exact solutions for nonlinear evolution equations. Compared with most existing tanh methods and other sophisticated methods, the proposed method not only recover some known solutions, but also find some new and general solutions. The solutions obtained in this paper include rational triangular periodic wave solutions, rational solitary wave solutions and rational wave solutions. The efficiency of the method can be demonstrated on (2 + 1)-dimensional Burgers equation

A Riccati solution for the ideal MHD plasma response with applications to real-time stability control

Science.gov (United States)

Glasser, Alexander; Kolemen, Egemen; Glasser, A. H.

2018-03-01

Active feedback control of ideal MHD stability in a tokamak requires rapid plasma stability analysis. Toward this end, we reformulate the δW stability method with a Hamilton-Jacobi theory, elucidating analytical and numerical features of the generic tokamak ideal MHD stability problem. The plasma response matrix is demonstrated to be the solution of an ideal MHD matrix Riccati differential equation. Since Riccati equations are prevalent in the control theory literature, such a shift in perspective brings to bear a range of numerical methods that are well-suited to the robust, fast solution of control problems. We discuss the usefulness of Riccati techniques in solving the stiff ordinary differential equations often encountered in ideal MHD stability analyses—for example, in tokamak edge and stellarator physics. We demonstrate the applicability of such methods to an existing 2D ideal MHD stability code—DCON [A. H. Glasser, Phys. Plasmas 23, 072505 (2016)]—enabling its parallel operation in near real-time, with wall-clock time ≪1 s . Such speed may help enable active feedback ideal MHD stability control, especially in tokamak plasmas whose ideal MHD equilibria evolve with inductive timescale τ≳ 1s—as in ITER.

Improved Riccati Transfer Matrix Method for Free Vibration of Non-Cylindrical Helical Springs Including Warping

Directory of Open Access Journals (Sweden)

A.M. Yu

2012-01-01

Full Text Available Free vibration equations for non-cylindrical (conical, barrel, and hyperboloidal types helical springs with noncircular cross-sections, which consist of 14 first-order ordinary differential equations with variable coefficients, are theoretically derived using spatially curved beam theory. In the formulation, the warping effect upon natural frequencies and vibrating mode shapes is first studied in addition to including the rotary inertia, the shear and axial deformation influences. The natural frequencies of the springs are determined by the use of improved Riccati transfer matrix method. The element transfer matrix used in the solution is calculated using the Scaling and Squaring method and Pad'e approximations. Three examples are presented for three types of springs with different cross-sectional shapes under clamped-clamped boundary condition. The accuracy of the proposed method has been compared with the FEM results using three-dimensional solid elements (Solid 45 in ANSYS code. Numerical results reveal that the warping effect is more pronounced in the case of non-cylindrical helical springs than that of cylindrical helical springs, which should be taken into consideration in the free vibration analysis of such springs.

Solving Matrix Equations on Multi-Core and Many-Core Architectures

Directory of Open Access Journals (Sweden)

Peter Benner

2013-11-01

Full Text Available We address the numerical solution of Lyapunov, algebraic and differential Riccati equations, via the matrix sign function, on platforms equipped with general-purpose multicore processors and, optionally, one or more graphics processing units (GPUs. In particular, we review the solvers for these equations, as well as the underlying methods, analyze their concurrency and scalability and provide details on their parallel implementation. Our experimental results show that this class of hardware provides sufficient computational power to tackle large-scale problems, which only a few years ago would have required a cluster of computers.

Backward Stochastic Riccati Equations and Infinite Horizon L-Q Optimal Control with Infinite Dimensional State Space and Random Coefficients

International Nuclear Information System (INIS)

Guatteri, Giuseppina; Tessitore, Gianmario

2008-01-01

We study the Riccati equation arising in a class of quadratic optimal control problems with infinite dimensional stochastic differential state equation and infinite horizon cost functional. We allow the coefficients, both in the state equation and in the cost, to be random.In such a context backward stochastic Riccati equations are backward stochastic differential equations in the whole positive real axis that involve quadratic non-linearities and take values in a non-Hilbertian space. We prove existence of a minimal non-negative solution and, under additional assumptions, its uniqueness. We show that such a solution allows to perform the synthesis of the optimal control and investigate its attractivity properties. Finally the case where the coefficients are stationary is addressed and an example concerning a controlled wave equation in random media is proposed

Exact Solutions for Fractional Differential-Difference Equations by an Extended Riccati Sub-ODE Method

International Nuclear Information System (INIS)

Feng Qinghua

2013-01-01

In this paper, an extended Riccati sub-ODE method is proposed to establish new exact solutions for fractional differential-difference equations in the sense of modified Riemann—Liouville derivative. By a fractional complex transformation, a given fractional differential-difference equation can be turned into another differential-difference equation of integer order. The validity of the method is illustrated by applying it to solve the fractional Hybrid lattice equation and the fractional relativistic Toda lattice system. As a result, some new exact solutions including hyperbolic function solutions, trigonometric function solutions and rational solutions are established. (general)

Experimental evaluation of optimal Vehicle Dynamic Control based on the State Dependent Riccati Equation technique

NARCIS (Netherlands)

Alirezaei, M.; Kanarachos, S.A.; Scheepers, B.T.M.; Maurice, J.P.

2013-01-01

Development and experimentally evaluation of an optimal Vehicle Dynamic Control (VDC) strategy based on the State Dependent Riccati Equation (SDRE) control technique is presented. The proposed nonlinear controller is based on a nonlinear vehicle model with nonlinear tire characteristics. A novel

THE EXISTENCE OF THE STABILIZING SOLUTION OF THE RICCATI EQUATION ARISING IN DISCRETE-TIME STOCHASTIC ZERO SUM LQ DYNAMIC GAMES WITH PERIODIC COEFFICIENTS

Directory of Open Access Journals (Sweden)

Vasile Dr ̆agan

2017-06-01

Full Text Available We investigate the problem for solving a discrete-time periodic gen- eralized Riccati equation with an indefinite sign of the quadratic term. A necessary condition for the existence of bounded and stabilizing solution of the discrete-time Riccati equation with an indefinite quadratic term is derived. The stabilizing solution is positive semidefinite and satisfies the introduced sign conditions. The proposed condition is illustrated via a numerical example.

State-dependent differential Riccati equation to track control of time-varying systems with state and control nonlinearities.

Science.gov (United States)

Korayem, M H; Nekoo, S R

2015-07-01

This work studies an optimal control problem using the state-dependent Riccati equation (SDRE) in differential form to track for time-varying systems with state and control nonlinearities. The trajectory tracking structure provides two nonlinear differential equations: the state-dependent differential Riccati equation (SDDRE) and the feed-forward differential equation. The independence of the governing equations and stability of the controller are proven along the trajectory using the Lyapunov approach. Backward integration (BI) is capable of solving the equations as a numerical solution; however, the forward solution methods require the closed-form solution to fulfill the task. A closed-form solution is introduced for SDDRE, but the feed-forward differential equation has not yet been obtained. Different ways of solving the problem are expressed and analyzed. These include BI, closed-form solution with corrective assumption, approximate solution, and forward integration. Application of the tracking problem is investigated to control robotic manipulators possessing rigid or flexible joints. The intention is to release a general program for automatic implementation of an SDDRE controller for any manipulator that obeys the Denavit-Hartenberg (D-H) principle when only D-H parameters are received as input data. Copyright © 2015 ISA. Published by Elsevier Ltd. All rights reserved.

Complex Riccati equations as a link between different approaches for the description of dissipative and irreversible systems

International Nuclear Information System (INIS)

Schuch, Dieter

2012-01-01

Quantum mechanics is essentially described in terms of complex quantities like wave functions. The interesting point is that phase and amplitude of the complex wave function are not independent of each other, but coupled by some kind of conservation law. This coupling exists in time-independent quantum mechanics and has a counterpart in its time-dependent form. It can be traced back to a reformulation of quantum mechanics in terms of nonlinear real Ermakov equations or equivalent complex nonlinear Riccati equations, where the quadratic term in the latter equation explains the origin of the phase-amplitude coupling. Since realistic physical systems are always in contact with some kind of environment this aspect is also taken into account. In this context, different approaches for describing open quantum systems, particularly effective ones, are discussed and compared. Certain kinds of nonlinear modifications of the Schrödinger equation are discussed as well as their interrelations and their relations to linear approaches via non-unitary transformations. The modifications of the aforementioned Ermakov and Riccati equations when environmental effects are included can be determined in the time-dependent case. From formal similarities conclusions can be drawn how the equations of time-independent quantum mechanics can be modified to also incluce the enviromental aspects.

Periodic Sturm-Liouville problems related to two Riccati equations of constant coefficients

International Nuclear Information System (INIS)

Khmelnytskaya, K.V.; Rosu, H.C.; Gonzalez, A.

2010-01-01

We consider two closely related Riccati equations of constant parameters whose particular solutions are used to construct the corresponding class of supersymmetrically coupled second-order differential equations. We solve analytically these parametric periodic problems along the whole real axis. Next, the analytically solved model is used as a case study for a powerful numerical approach that is employed here for the first time in the investigation of the energy band structure of periodic not necessarily regular potentials. The approach is based on the well-known self-matching procedure of James (1949) and implements the spectral parameter power series solutions introduced by Kravchenko (2008). We obtain additionally an efficient series representation of the Hill discriminant based on Kravchenko's series.

On the Solution of the Rational Matrix Equation X=Q+LX−1LT

Directory of Open Access Journals (Sweden)

Heike Faßbender

2007-01-01

Full Text Available We study numerical methods for finding the maximal symmetric positive definite solution of the nonlinear matrix equation X=Q+LX−1LT, where Q is symmetric positive definite and L is nonsingular. Such equations arise for instance in the analysis of stationary Gaussian reciprocal processes over a finite interval. Its unique largest positive definite solution coincides with the unique positive definite solution of a related discrete-time algebraic Riccati equation (DARE. We discuss how to use the butterfly SZ algorithm to solve the DARE. This approach is compared to several fixed-point and doubling-type iterative methods suggested in the literature.

New Trace Bounds for the Product of Two Matrices and Their Applications in the Algebraic Riccati Equation

Directory of Open Access Journals (Sweden)

Liu Jianzhou

2009-01-01

Full Text Available By using singular value decomposition and majorization inequalities, we propose new inequalities for the trace of the product of two arbitrary real square matrices. These bounds improve and extend the recent results. Further, we give their application in the algebraic Riccati equation. Finally, numerical examples have illustrated that our results are effective and superior.

A numerical technique for solving fractional optimal control problems and fractional Riccati differential equations

Directory of Open Access Journals (Sweden)

F. Ghomanjani

2016-10-01

Full Text Available In the present paper, we apply the Bezier curves method for solving fractional optimal control problems (OCPs and fractional Riccati differential equations. The main advantage of this method is that it can reduce the error of the approximate solutions. Hence, the solutions obtained using the Bezier curve method give good approximations. Some numerical examples are provided to confirm the accuracy of the proposed method. All of the numerical computations have been performed on a PC using several programs written in MAPLE 13.

New exact solutions of the(2+1-dimensional Broer-Kaup equation by the consistent Riccati expansion method

Directory of Open Access Journals (Sweden)

Jiang Ying

2017-01-01

Full Text Available In this work, we study the (2+1-D Broer-Kaup equation. The composite periodic breather wave, the exact composite kink breather wave and the solitary wave solutions are obtained by using the coupled degradation technique and the consistent Riccati expansion method. These results may help us to investigate some complex dynamical behaviors and the interaction between composite non-linear waves in high dimensional models

Newton-Raphson based modified Laplace Adomian decomposition method for solving quadratic Riccati differential equations

Directory of Open Access Journals (Sweden)

Mishra Vinod

2016-01-01

Full Text Available Numerical Laplace transform method is applied to approximate the solution of nonlinear (quadratic Riccati differential equations mingled with Adomian decomposition method. A new technique is proposed in this work by reintroducing the unknown function in Adomian polynomial with that of well known Newton-Raphson formula. The solutions obtained by the iterative algorithm are exhibited in an infinite series. The simplicity and efficacy of method is manifested with some examples in which comparisons are made among the exact solutions, ADM (Adomian decomposition method, HPM (Homotopy perturbation method, Taylor series method and the proposed scheme.

On the properties of a variant of the Riccati system of equations

International Nuclear Information System (INIS)

Sarkar, Amartya; Guha, Partha; Bhattacharjee, J K; Mallik, A K; Ghose-Choudhury, Anindya; Leach, P G L

2012-01-01

A variant of the generalized Riccati system of equations is considered. It is shown that for α = 2n + 3 the system admits a bilagrangian description and the dynamics has a node at the origin, whereas for α much smaller than a critical value the dynamics is periodic, the origin being a centre. It is found that the solution changes from being periodic to aperiodic at a critical point, α c = 2√(2(n+1)), which is independent of the initial conditions. This behaviour is explained by finding a scaling argument via which the phase trajectories corresponding to different initial conditions collapse onto a single universal orbit. Numerical evidence for the transition is shown. Further, using a perturbative renormalization group argument, it is conjectured that the oscillator exhibits isochronous oscillations. The correctness of the conjecture is established numerically. (paper)

Riccati-parameter solutions of nonlinear second-order ODEs

International Nuclear Information System (INIS)

Reyes, M A; Rosu, H C

2008-01-01

It has been proven by Rosu and Cornejo-Perez (Rosu and Cornejo-Perez 2005 Phys. Rev. E 71 046607, Cornejo-Perez and Rosu 2005 Prog. Theor. Phys. 114 533) that for some nonlinear second-order ODEs it is a very simple task to find one particular solution once the nonlinear equation is factorized with the use of two first-order differential operators. Here, it is shown that an interesting class of parametric solutions is easy to obtain if the proposed factorization has a particular form, which happily turns out to be the case in many problems of physical interest. The method that we exemplify with a few explicitly solved cases consists in using the general solution of the Riccati equation, which contributes with one parameter to this class of parametric solutions. For these nonlinear cases, the Riccati parameter serves as a 'growth' parameter from the trivial null solution up to the particular solution found through the factorization procedure

On equations of motion on complex grassman manifold

International Nuclear Information System (INIS)

Berceanu, S.; Gheorghe, A.

1989-02-01

We investigate the equations of motion on the 'classical' phase space which corresponds to quantum state space in the case of the complex Grassmann manifold appearing in the Hartree-Fock problem. First and second degree polynomial Hamiltonians in bifermion operators are considered. The 'classical' motion corresponding to linear Hamiltonians is described by a Matrix Riccati equation.(authors)

A modified variable-coefficient projective Riccati equation method and its application to (2 + 1)-dimensional simplified generalized Broer-Kaup system

International Nuclear Information System (INIS)

Liu Qing; Zhu Jiamin; Hong Bihai

2008-01-01

A modified variable-coefficient projective Riccati equation method is proposed and applied to a (2 + 1)-dimensional simplified and generalized Broer-Kaup system. It is shown that the method presented by Huang and Zhang [Huang DJ, Zhang HQ. Chaos, Solitons and Fractals 2005; 23:601] is a special case of our method. The results obtained in the paper include many new formal solutions besides the all solutions found by Huang and Zhang

Strong convergence and convergence rates of approximating solutions for algebraic Riccati equations in Hilbert spaces

Science.gov (United States)

Ito, Kazufumi

1987-01-01

The linear quadratic optimal control problem on infinite time interval for linear time-invariant systems defined on Hilbert spaces is considered. The optimal control is given by a feedback form in terms of solution pi to the associated algebraic Riccati equation (ARE). A Ritz type approximation is used to obtain a sequence pi sup N of finite dimensional approximations of the solution to ARE. A sufficient condition that shows pi sup N converges strongly to pi is obtained. Under this condition, a formula is derived which can be used to obtain a rate of convergence of pi sup N to pi. The results of the Galerkin approximation is demonstrated and applied for parabolic systems and the averaging approximation for hereditary differential systems.

Parametric oscillators from factorizations employing a constant-shifted Riccati solution of the classical harmonic oscillator

Energy Technology Data Exchange (ETDEWEB)

Rosu, H.C., E-mail: hcr@ipicyt.edu.mx [IPICyT, Instituto Potosino de Investigacion Cientifica y Tecnologica, Apdo Postal 3-74 Tangamanga, 78231 San Luis Potosi, S.L.P. (Mexico); Khmelnytskaya, K.V. [Universidad Autonoma de Queretaro, Centro Universitario, Cerro de las Campanas s/n, C.P. 76010 Santiago de Queretaro, Qro. (Mexico)

2011-09-19

We determine the kind of parametric oscillators that are generated in the usual factorization procedure of second-order linear differential equations when one introduces a constant shift of the Riccati solution of the classical harmonic oscillator. The mathematical results show that some of these oscillators could be of physical nature. We give the solutions of the obtained second-order differential equations and the values of the shift parameter providing strictly periodic and antiperiodic solutions. We also notice that this simple problem presents parity-time (PT) symmetry. Possible applications are mentioned. -- Highlights: → A particular Riccati solution of the classical harmonic oscillator is shifted by a constant. → Such a solution is used in the factorization brackets to get different equations of motion. → The properties of the parametric oscillators obtained in this way are examined.

The constrained discrete-time state-dependent Riccati equation technique for uncertain nonlinear systems

Science.gov (United States)

Chang, Insu

The objective of the thesis is to introduce a relatively general nonlinear controller/estimator synthesis framework using a special type of the state-dependent Riccati equation technique. The continuous time state-dependent Riccati equation (SDRE) technique is extended to discrete-time under input and state constraints, yielding constrained (C) discrete-time (D) SDRE, referred to as CD-SDRE. For the latter, stability analysis and calculation of a region of attraction are carried out. The derivation of the D-SDRE under state-dependent weights is provided. Stability of the D-SDRE feedback system is established using Lyapunov stability approach. Receding horizon strategy is used to take into account the constraints on D-SDRE controller. Stability condition of the CD-SDRE controller is analyzed by using a switched system. The use of CD-SDRE scheme in the presence of constraints is then systematically demonstrated by applying this scheme to problems of spacecraft formation orbit reconfiguration under limited performance on thrusters. Simulation results demonstrate the efficacy and reliability of the proposed CD-SDRE. The CD-SDRE technique is further investigated in a case where there are uncertainties in nonlinear systems to be controlled. First, the system stability under each of the controllers in the robust CD-SDRE technique is separately established. The stability of the closed-loop system under the robust CD-SDRE controller is then proven based on the stability of each control system comprising switching configuration. A high fidelity dynamical model of spacecraft attitude motion in 3-dimensional space is derived with a partially filled fuel tank, assumed to have the first fuel slosh mode. The proposed robust CD-SDRE controller is then applied to the spacecraft attitude control system to stabilize its motion in the presence of uncertainties characterized by the first fuel slosh mode. The performance of the robust CD-SDRE technique is discussed. Subsequently

Hartree--Fock density matrix equation

International Nuclear Information System (INIS)

Cohen, L.; Frishberg, C.

1976-01-01

An equation for the Hartree--Fock density matrix is discussed and the possibility of solving this equation directly for the density matrix instead of solving the Hartree--Fock equation for orbitals is considered. Toward that end the density matrix is expanded in a finite basis to obtain the matrix representative equation. The closed shell case is considered. Two numerical schemes are developed and applied to a number of examples. One example is given where the standard orbital method does not converge while the method presented here does

Maximal imaginery eigenvalues in optimal systems

Directory of Open Access Journals (Sweden)

David Di Ruscio

1991-07-01

Full Text Available In this note we present equations that uniquely determine the maximum possible imaginary value of the closed loop eigenvalues in an LQ-optimal system, irrespective of how the state weight matrix is chosen, provided a real symmetric solution of the algebraic Riccati equation exists. In addition, the corresponding state weight matrix and the solution to the algebraic Riccati equation are derived for a class of linear systems. A fundamental lemma for the existence of a real symmetric solution to the algebraic Riccati equation is derived for this class of linear systems.

Minimal solution for inconsistent singular fuzzy matrix equations

Directory of Open Access Journals (Sweden)

M. Nikuie

2013-10-01

Full Text Available The fuzzy matrix equations $Ailde{X}=ilde{Y}$ is called a singular fuzzy matrix equations while the coefficients matrix of its equivalent crisp matrix equations be a singular matrix. The singular fuzzy matrix equations are divided into two parts: consistent singular matrix equations and inconsistent fuzzy matrix equations. In this paper, the inconsistent singular fuzzy matrix equations is studied and the effect of generalized inverses in finding minimal solution of an inconsistent singular fuzzy matrix equations are investigated.

On matrix fractional differential equations

Directory of Open Access Journals (Sweden)

Adem Kılıçman

2017-01-01

Full Text Available The aim of this article is to study the matrix fractional differential equations and to find the exact solution for system of matrix fractional differential equations in terms of Riemann–Liouville using Laplace transform method and convolution product to the Riemann–Liouville fractional of matrices. Also, we show the theorem of non-homogeneous matrix fractional partial differential equation with some illustrative examples to demonstrate the effectiveness of the new methodology. The main objective of this article is to discuss the Laplace transform method based on operational matrices of fractional derivatives for solving several kinds of linear fractional differential equations. Moreover, we present the operational matrices of fractional derivatives with Laplace transform in many applications of various engineering systems as control system. We present the analytical technique for solving fractional-order, multi-term fractional differential equation. In other words, we propose an efficient algorithm for solving fractional matrix equation.

On matrix fractional differential equations

OpenAIRE

Adem Kılıçman; Wasan Ajeel Ahmood

2017-01-01

The aim of this article is to study the matrix fractional differential equations and to find the exact solution for system of matrix fractional differential equations in terms of Riemann–Liouville using Laplace transform method and convolution product to the Riemann–Liouville fractional of matrices. Also, we show the theorem of non-homogeneous matrix fractional partial differential equation with some illustrative examples to demonstrate the effectiveness of the new methodology. The main objec...

Delay-differential equations and the Painlevé transcendents

Science.gov (United States)

Grammaticos, B.; Ramani, A.; Moreira, I. C.

1993-07-01

We apply the recently proposed integrability criterion for differential-difference systems (that blends the classical Painlevé analysis with singularity confinement for discrete systems) to a class of first-order differential-delay equations. Our analysis singles out the family of bi-Riccati equations, as integrability candidates. Among these equations that pass the test some are integrable in a straightforward way (usually by reduction to a standard Riccati equation for some transformed variable) while the remaining ones define new hysterodifferential forms of the Painlevé transcendental equations.

The modified extended Fan's sub-equation method and its application to (2 + 1)-dimensional dispersive long wave equation

International Nuclear Information System (INIS)

Yomba, Emmanuel

2005-01-01

By using a modified extended Fan's sub-equation method, we have obtained new and more general solutions including a series of non-travelling wave and coefficient function solutions namely: soliton-like solutions, triangular-like solutions, single and combined non-degenerative Jacobi elliptic wave function-like solutions for the (2 + 1)-dimensional dispersive long wave equation. The most important achievement of this method lies on the fact that, we have succeeded in one move to give all the solutions which can be previously obtained by application of at least four methods (method using Riccati equation, or first kind elliptic equation, or auxiliary ordinary equation, or generalized Riccati equation as mapping equation)

New analytic solutions of stochastic coupled KdV equations

International Nuclear Information System (INIS)

Dai Chaoqing; Chen Junlang

2009-01-01

In this paper, firstly, we use the exp-function method to seek new exact solutions of the Riccati equation. Then, with the help of Hermit transformation, we employ the Riccati equation and its new exact solutions to find new analytic solutions of the stochastic coupled KdV equation in the white noise environment. As some special examples, some analytic solutions can degenerate into these solutions reported in open literatures.

Traveling wave solutions for two nonlinear evolution equations with nonlinear terms of any order

International Nuclear Information System (INIS)

Feng Qing-Hua; Zhang Yao-Ming; Meng Fan-Wei

2011-01-01

In this paper, based on the known first integral method and the Riccati sub-ordinary differential equation (ODE) method, we try to seek the exact solutions of the general Gardner equation and the general Benjamin—Bona—Mahoney equation. As a result, some traveling wave solutions for the two nonlinear equations are established successfully. Also we make a comparison between the two methods. It turns out that the Riccati sub-ODE method is more effective than the first integral method in handling the proposed problems, and more general solutions are constructed by the Riccati sub-ODE method. (general)

The extended Fan's sub-equation method and its application to KdV-MKdV, BKK and variant Boussinesq equations

International Nuclear Information System (INIS)

Yomba, Emmanuel

2005-01-01

An extended Fan's sub-equation method is used for constructing exact travelling wave solutions of nonlinear partial differential equations (NLPDEs). The key idea of this method is to take full advantage of the general elliptic equation involving five parameters which has more new solutions and whose degeneracies can lead to special sub-equations involving three parameters. More new solutions are obtained for KdV-MKdV, Broer-Kaup-Kupershmidt (BKK) and variant Boussinesq equations. Then we present a technique which not only gives us a clear relation among this general elliptic equation and other sub-equations involving three parameters (Riccati equation, first kind elliptic equation, auxiliary ordinary equation, generalized Riccati equation and so on), but also provides an approach to construct new exact solutions to NLPDEs

New types of exact solutions for a breaking soliton equation

International Nuclear Information System (INIS)

Mei Jianqin; Zhang Hongqing

2004-01-01

In this paper based on a system of Riccati equations, we present a newly generally projective Riccati equation expansion method and its algorithm, which can be used to construct more new exact solutions of nonlinear differential equations in mathematical physics. A typical breaking soliton equation is chosen to illustrate our algorithm such that more families of new exact solutions are obtained, which contain soliton-like solutions and periodic solutions. This algorithm can also be applied to other nonlinear differential equations

New traveling wave solutions to AKNS and SKdV equations

International Nuclear Information System (INIS)

Ozer, Teoman

2009-01-01

We analyze the traveling wave solutions of Ablowitz-Kaup-Newell-Segur (AKNS) and Schwarz-Korteweg-de Vries (SKdV) equations. As the solution method for differential equations we consider the improved tanh approach. This approach provides to transform the partial differential equation into the ordinary differential equation and then obtain the new families of exact solutions based on the solutions of the Riccati equation. The different values of the coefficients of the Riccati equation allow us to obtain new type of traveling wave solutions to AKNS and SKdV equations.

On Solutions for Linear and Nonlinear Schrödinger Equations with Variable Coefficients: A Computational Approach

Directory of Open Access Journals (Sweden)

Gabriel Amador

2016-05-01

Full Text Available In this work, after reviewing two different ways to solve Riccati systems, we are able to present an extensive list of families of integrable nonlinear Schrödinger (NLS equations with variable coefficients. Using Riccati equations and similarity transformations, we are able to reduce them to the standard NLS models. Consequently, we can construct bright-, dark- and Peregrine-type soliton solutions for NLS with variable coefficients. As an important application of solutions for the Riccati equation with parameters, by means of computer algebra systems, it is shown that the parameters change the dynamics of the solutions. Finally, we test numerical approximations for the inhomogeneous paraxial wave equation by the Crank-Nicolson scheme with analytical solutions found using Riccati systems. These solutions include oscillating laser beams and Laguerre and Gaussian beams.

A generalization of the simplest equation method and its application to (3+1)-dimensional KP equation and generalized Fisher equation

International Nuclear Information System (INIS)

Zhao, Zhonglong; Zhang, Yufeng; Han, Zhong; Rui, Wenjuan

2014-01-01

In this paper, the simplest equation method is used to construct exact traveling solutions of the (3+1)-dimensional KP equation and generalized Fisher equation. We summarize the main steps of the simplest equation method. The Bernoulli and Riccati equation are used as simplest equations. This method is straightforward and concise, and it can be applied to other nonlinear partial differential equations

About the solvability of matrix polynomial equations

OpenAIRE

Netzer, Tim; Thom, Andreas

2016-01-01

We study self-adjoint matrix polynomial equations in a single variable and prove existence of self-adjoint solutions under some assumptions on the leading form. Our main result is that any self-adjoint matrix polynomial equation of odd degree with non-degenerate leading form can be solved in self-adjoint matrices. We also study equations of even degree and equations in many variables.

Some remarks on unilateral matrix equations

International Nuclear Information System (INIS)

Cerchiai, Bianca L.; Zumino, Bruno

2001-01-01

We briefly review the results of our paper LBNL-46775: We study certain solutions of left-unilateral matrix equations. These are algebraic equations where the coefficients and the unknown are square matrices of the same order, or, more abstractly, elements of an associative, but possibly noncommutative algebra, and all coefficients are on the left. Recently such equations have appeared in a discussion of generalized Born-Infeld theories. In particular, two equations, their perturbative solutions and the relation between them are studied, applying a unified approach based on the generalized Bezout theorem for matrix polynomials

Oscillation theory of linear differential equations

Czech Academy of Sciences Publication Activity Database

Došlý, Ondřej

2000-01-01

Roč. 36, č. 5 (2000), s. 329-343 ISSN 0044-8753 R&D Projects: GA ČR GA201/98/0677 Keywords : discrete oscillation theory %Sturm-Liouville equation%Riccati equation Subject RIV: BA - General Mathematics

OSCILLATION CRITERIA FOR FORCED SUPERLINEAR DIFFERENCE EQUATIONS

Institute of Scientific and Technical Information of China (English)

无

2008-01-01

Using Riccati transformation techniques,some oscillation criteria for the forced second-order superlinear difference equations are established.These criteria are dis- crete analogues of the criteria for differential equations proposed by Yan.

Simple equation method for nonlinear partial differential equations and its applications

Directory of Open Access Journals (Sweden)

Taher A. Nofal

2016-04-01

Full Text Available In this article, we focus on the exact solution of the some nonlinear partial differential equations (NLPDEs such as, Kodomtsev–Petviashvili (KP equation, the (2 + 1-dimensional breaking soliton equation and the modified generalized Vakhnenko equation by using the simple equation method. In the simple equation method the trial condition is the Bernoulli equation or the Riccati equation. It has been shown that the method provides a powerful mathematical tool for solving nonlinear wave equations in mathematical physics and engineering problems.

Approximate Solution of LR Fuzzy Sylvester Matrix Equations

Directory of Open Access Journals (Sweden)

Xiaobin Guo

2013-01-01

Full Text Available The fuzzy Sylvester matrix equation AX~+X~B=C~ in which A,B are m×m and n×n crisp matrices, respectively, and C~ is an m×n LR fuzzy numbers matrix is investigated. Based on the Kronecker product of matrices, we convert the fuzzy Sylvester matrix equation into an LR fuzzy linear system. Then we extend the fuzzy linear system into two systems of linear equations according to the arithmetic operations of LR fuzzy numbers. The fuzzy approximate solution of the original fuzzy matrix equation is obtained by solving the crisp linear systems. The existence condition of the LR fuzzy solution is also discussed. Some examples are given to illustrate the proposed method.

New multiple soliton solutions to the general Burgers-Fisher equation and the Kuramoto-Sivashinsky equation

International Nuclear Information System (INIS)

Chen Huaitang; Zhang Hongqing

2004-01-01

A generalized tanh function method is used for constructing exact travelling wave solutions of nonlinear partial differential equations in a unified way. The main idea of this method is to take full advantage of the Riccati equation which has more new solutions. More new multiple soliton solutions are obtained for the general Burgers-Fisher equation and the Kuramoto-Sivashinsky equation

Exact solutions of the Drinfel'd–Sokolov–Wilson equation using ...

Indian Academy of Sciences (India)

1Department of Engineering Mathematics and Physics, Higher Institute of Engineering, .... ordinary differential equation (NLODE) of the form. P .... 3.1 The Bäcklund transformation of Riccati equation method applied to DSW equation.

Exact traveling wave solutions of the Boussinesq equation

International Nuclear Information System (INIS)

Ding Shuangshuang; Zhao Xiqiang

2006-01-01

The repeated homogeneous balance method is used to construct exact traveling wave solutions of the Boussinesq equation, in which the homogeneous balance method is applied to solve the Riccati equation and the reduced nonlinear ordinary differential equation, respectively. Many new exact traveling wave solutions of the Boussinesq equation are successfully obtained

Weierstrass Elliptic Function Solutions to Nonlinear Evolution Equations

International Nuclear Information System (INIS)

Yu Jianping; Sun Yongli

2008-01-01

This paper is based on the relations between projection Riccati equations and Weierstrass elliptic equation, combined with the Groebner bases in the symbolic computation. Then the novel method for constructing the Weierstrass elliptic solutions to the nonlinear evolution equations is given by using the above relations

A new sub-equation method applied to obtain exact travelling wave solutions of some complex nonlinear equations

International Nuclear Information System (INIS)

Zhang Huiqun

2009-01-01

By using a new coupled Riccati equations, a direct algebraic method, which was applied to obtain exact travelling wave solutions of some complex nonlinear equations, is improved. And the exact travelling wave solutions of the complex KdV equation, Boussinesq equation and Klein-Gordon equation are investigated using the improved method. The method presented in this paper can also be applied to construct exact travelling wave solutions for other nonlinear complex equations.

Efficient Implementation of the Riccati Recursion for Solving Linear-Quadratic Control Problems

DEFF Research Database (Denmark)

Frison, Gianluca; Jørgensen, John Bagterp

2013-01-01

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 typically the main computational effort at each iteration....... In this paper, we compare a number of solvers for an extended formulation of the LQ control problem: a Riccati recursion based solver can be considered the best choice for the general problem with dense matrices. Furthermore, we present a novel version of the Riccati solver, that makes use of the Cholesky...... factorization of the Pn matrices to reduce the number of flops. When combined with regularization and mixed precision, this algorithm can solve large instances of the LQ control problem up to 3 times faster than the classical Riccati solver....

A note on Verhulst's logistic equation and related logistic maps

International Nuclear Information System (INIS)

Gutierrez, M Ranferi; Reyes, M A; Rosu, H C

2010-01-01

We consider the Verhulst logistic equation and a couple of forms of the corresponding logistic maps. For the case of the logistic equation we show that using the general Riccati solution only changes the initial conditions of the equation. Next, we consider two forms of corresponding logistic maps reporting the following results. For the map x n+1 = rx n (1 - x n ) we propose a new way to write the solution for r = -2 which allows better precision of the iterative terms, while for the map x n+1 - x n = rx n (1 - x n+1 ) we show that it behaves identically to the logistic equation from the standpoint of the general Riccati solution, which is also provided herein for any value of the parameter r.

Asymptotic properties for half-linear difference equations

Czech Academy of Sciences Publication Activity Database

Cecchi, M.; Došlá, Z.; Marini, M.; Vrkoč, Ivo

2006-01-01

Roč. 131, č. 4 (2006), s. 347-363 ISSN 0862-7959 R&D Projects: GA ČR(CZ) GA201/04/0580 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear second order difference equation * nonoscillatory solutions * Riccati difference equation Subject RIV: BA - General Mathematics

Differential equations extended to superspace

International Nuclear Information System (INIS)

Torres, J.; Rosu, H.C.

2003-01-01

We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)

Differential equations extended to superspace

Energy Technology Data Exchange (ETDEWEB)

Torres, J. [Instituto de Fisica, Universidad de Guanajuato, A.P. E-143, Leon, Guanajuato (Mexico); Rosu, H.C. [Instituto Potosino de Investigacion Cientifica y Tecnologica, A.P. 3-74, Tangamanga, San Luis Potosi (Mexico)

2003-07-01

We present a simple SUSY Ns = 2 superspace extension of the differential equations in which the sought solutions are considered to be real superfields but maintaining the common derivative operators and the coefficients of the differential equations unaltered. In this way, we get self consistent systems of coupled differential equations for the components of the superfield. This procedure is applied to the Riccati equation, for which we obtain in addition the system of coupled equations corresponding to the components of the general superfield solution. (Author)

A note on Verhulst's logistic equation and related logistic maps

Energy Technology Data Exchange (ETDEWEB)

Gutierrez, M Ranferi; Reyes, M A [Depto de Fisica, Universidad de Guanajuato, Apdo. Postal E143, 37150 Leon, Gto. (Mexico); Rosu, H C, E-mail: hcr@ipicyt.edu.m [IPICyT, Instituto Potosino de Investigacion Cientifica y Tecnologica, Apdo Postal 3-74 Tangamanga, 78231 San Luis PotosI (Mexico)

2010-05-21

We consider the Verhulst logistic equation and a couple of forms of the corresponding logistic maps. For the case of the logistic equation we show that using the general Riccati solution only changes the initial conditions of the equation. Next, we consider two forms of corresponding logistic maps reporting the following results. For the map x{sub n+1} = rx{sub n}(1 - x{sub n}) we propose a new way to write the solution for r = -2 which allows better precision of the iterative terms, while for the map x{sub n+1} - x{sub n} = rx{sub n}(1 - x{sub n+1}) we show that it behaves identically to the logistic equation from the standpoint of the general Riccati solution, which is also provided herein for any value of the parameter r.

A Riccati model for Denmark Strait overflow variability

Science.gov (United States)

Käse, R. H.

2006-10-01

A controlled volume box model of the western basins of the Nordic Seas for water denser than 1027.8 kg m-3 is constructed, where accumulation in volume ($\\frac{dV}{dt) is driven by net imbalances between prescribed net inflow from the northern, eastern and top boundaries (Qs) and hydraulically limited outflow through the Denmark Strait. The resulting Riccati equation is solved analytically for filling and flushing experiments with constant Qs and numerically for stochastic forcing Qs(t). For small perturbations to Qs with white noise spectrum, the overflow response is red noise with a time scale between 5 and 15 years depending on the mean interface height and area. For Qs proportional to the NAO index, the overflow is positively correlated with the NAO. A 140 years integration reveals variations in the overflow between 2.5 Sv in the 1970s and a maximum of 4 Sv in the 1990s. Hydraulic transport calculations from hydrographic data north of Iceland show good agreement with the model hindcast.

Galois action on solutions of a differential equation

NARCIS (Netherlands)

Hendriks, Peter A.; Put, Marius van der

Consider a second-order differential equation of the form y '' + ay' + by = O with a,b is an element of Q(x). Kovacic's algorithm tries to compute a solution of the associated Riccati equation that is algebraic and of minimal degree over (Q) over bar(x). The coefficients of the monic irreducible

An algorithm for solving an arbitrary triangular fully fuzzy Sylvester matrix equations

Science.gov (United States)

Daud, Wan Suhana Wan; Ahmad, Nazihah; Malkawi, Ghassan

2017-11-01

Sylvester matrix equations played a prominent role in various areas including control theory. Considering to any un-certainty problems that can be occurred at any time, the Sylvester matrix equation has to be adapted to the fuzzy environment. Therefore, in this study, an algorithm for solving an arbitrary triangular fully fuzzy Sylvester matrix equation is constructed. The construction of the algorithm is based on the max-min arithmetic multiplication operation. Besides that, an associated arbitrary matrix equation is modified in obtaining the final solution. Finally, some numerical examples are presented to illustrate the proposed algorithm.

Forced oscillation of hyperbolic equations with mixed nonlinearities

Directory of Open Access Journals (Sweden)

Yutaka Shoukaku

2012-04-01

Full Text Available In this paper, we consider the mixed nonlinear hyperbolic equations with forcing term via Riccati inequality. Some sufficient conditions for the oscillation are derived by using Young inequality and integral averaging method.

String beta function equations from c=1 matrix model

CERN Document Server

Dhar, A; Wadia, S R; Dhar, Avinash; Mandal, Gautam; Wadia, Spenta R

1995-01-01

We derive the \\sigma-model tachyon \\beta-function equation of 2-dimensional string theory, in the background of flat space and linear dilaton, working entirely within the c=1 matrix model. The tachyon \\beta-function equation is satisfied by a \\underbar{nonlocal} and \\underbar{nonlinear} combination of the (massless) scalar field of the matrix model. We discuss the possibility of describing the `discrete states' as well as other possible gravitational and higher tensor backgrounds of 2-dimensional string theory within the c=1 matrix model. We also comment on the realization of the W-infinity symmetry of the matrix model in the string theory. The present work reinforces the viewpoint that a nonlocal (and nonlinear) transform is required to extract the space-time physics of 2-dimensional string theory from the c=1 matrix model.

Minimal parameter solution of the orthogonal matrix differential equation

Science.gov (United States)

Bar-Itzhack, Itzhack Y.; Markley, F. Landis

1990-01-01

As demonstrated in this work, all orthogonal matrices solve a first order differential equation. The straightforward solution of this equation requires n sup 2 integrations to obtain the element of the nth order matrix. There are, however, only n(n-1)/2 independent parameters which determine an orthogonal matrix. The questions of choosing them, finding their differential equation and expressing the orthogonal matrix in terms of these parameters are considered. Several possibilities which are based on attitude determination in three dimensions are examined. It is shown that not all 3-D methods have useful extensions to higher dimensions. It is also shown why the rate of change of the matrix elements, which are the elements of the angular rate vector in 3-D, are the elements of a tensor of the second rank (dyadic) in spaces other than three dimensional. It is proven that the 3-D Gibbs vector (or Cayley Parameters) are extendable to other dimensions. An algorithm is developed emplying the resulting parameters, which are termed Extended Rodrigues Parameters, and numerical results are presented of the application of the algorithm to a fourth order matrix.

Exact solutions for the higher-order nonlinear Schoerdinger equation in nonlinear optical fibres

International Nuclear Information System (INIS)

Liu Chunping

2005-01-01

First, by using the generally projective Riccati equation method, many kinds of exact solutions for the higher-order nonlinear Schoerdinger equation in nonlinear optical fibres are obtained in a unified way. Then, some relations among these solutions are revealed

The generalised Sylvester matrix equations over the generalised bisymmetric and skew-symmetric matrices

Science.gov (United States)

Dehghan, Mehdi; Hajarian, Masoud

2012-08-01

A matrix P is called a symmetric orthogonal if P = P T = P -1. A matrix X is said to be a generalised bisymmetric with respect to P if X = X T = PXP. It is obvious that any symmetric matrix is also a generalised bisymmetric matrix with respect to I (identity matrix). By extending the idea of the Jacobi and the Gauss-Seidel iterations, this article proposes two new iterative methods, respectively, for computing the generalised bisymmetric (containing symmetric solution as a special case) and skew-symmetric solutions of the generalised Sylvester matrix equation ? (including Sylvester and Lyapunov matrix equations as special cases) which is encountered in many systems and control applications. When the generalised Sylvester matrix equation has a unique generalised bisymmetric (skew-symmetric) solution, the first (second) iterative method converges to the generalised bisymmetric (skew-symmetric) solution of this matrix equation for any initial generalised bisymmetric (skew-symmetric) matrix. Finally, some numerical results are given to illustrate the effect of the theoretical results.

A New Pseudoinverse Matrix Method For Balancing Chemical Equations And Their Stability

International Nuclear Information System (INIS)

Risteski, Ice B.

2008-01-01

In this work is given a new pseudoniverse matrix method for balancing chemical equations. Here offered method is founded on virtue of the solution of a Diophantine matrix equation by using of a Moore-Penrose pseudoinverse matrix. The method has been tested on several typical chemical equations and found to be very successful for the all equations in our extensive balancing research. This method, which works successfully without any limitations, also has the capability to determine the feasibility of a new chemical reaction, and if it is feasible, then it will balance the equation. Chemical equations treated here possess atoms with fractional oxidation numbers. Also, in the present work are introduced necessary and sufficient criteria for stability of chemical equations over stability of their extended matrices

Linear matrix differential equations of higher-order and applications

Directory of Open Access Journals (Sweden)

Mustapha Rachidi

2008-07-01

Full Text Available In this article, we study linear differential equations of higher-order whose coefficients are square matrices. The combinatorial method for computing the matrix powers and exponential is adopted. New formulas representing auxiliary results are obtained. This allows us to prove properties of a large class of linear matrix differential equations of higher-order, in particular results of Apostol and Kolodner are recovered. Also illustrative examples and applications are presented.

Online gaming for learning optimal team strategies in real time

Science.gov (United States)

Hudas, Gregory; Lewis, F. L.; Vamvoudakis, K. G.

2010-04-01

This paper first presents an overall view for dynamical decision-making in teams, both cooperative and competitive. Strategies for team decision problems, including optimal control, zero-sum 2-player games (H-infinity control) and so on are normally solved for off-line by solving associated matrix equations such as the Riccati equation. However, using that approach, players cannot change their objectives online in real time without calling for a completely new off-line solution for the new strategies. Therefore, in this paper we give a method for learning optimal team strategies online in real time as team dynamical play unfolds. In the linear quadratic regulator case, for instance, the method learns the Riccati equation solution online without ever solving the Riccati equation. This allows for truly dynamical team decisions where objective functions can change in real time and the system dynamics can be time-varying.

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.$ 

The matrix nonlinear Schrodinger equation in dimension 2

DEFF Research Database (Denmark)

Zuhan, L; Pedersen, Michael

2001-01-01

In this paper we study the existence of global solutions to the Cauchy problem for the matrix nonlinear Schrodinger equation (MNLS) in 2 space dimensions. A sharp condition for the global existence is obtained for this equation. This condition is in terms of an exact stationary solution...... of a semilinear elliptic equation. In the scalar case, the MNLS reduces to the well-known cubic nonlinear Schrodinger equation for which existence of solutions has been studied by many authors. (C) 2001 Academic Press....

Three Interpretations of the Matrix Equation Ax = b

Science.gov (United States)

Larson, Christine; Zandieh, Michelle

2013-01-01

Many of the central ideas in an introductory undergraduate linear algebra course are closely tied to a set of interpretations of the matrix equation Ax = b (A is a matrix, x and b are vectors): linear combination interpretations, systems interpretations, and transformation interpretations. We consider graphic and symbolic representations for each,…

Lax representations for matrix short pulse equations

Science.gov (United States)

Popowicz, Z.

2017-10-01

The Lax representation for different matrix generalizations of Short Pulse Equations (SPEs) is considered. The four-dimensional Lax representations of four-component Matsuno, Feng, and Dimakis-Müller-Hoissen-Matsuno equations are obtained. The four-component Feng system is defined by generalization of the two-dimensional Lax representation to the four-component case. This system reduces to the original Feng equation, to the two-component Matsuno equation, or to the Yao-Zang equation. The three-component version of the Feng equation is presented. The four-component version of the Matsuno equation with its Lax representation is given. This equation reduces the new two-component Feng system. The two-component Dimakis-Müller-Hoissen-Matsuno equations are generalized to the four-parameter family of the four-component SPE. The bi-Hamiltonian structure of this generalization, for special values of parameters, is defined. This four-component SPE in special cases reduces to the new two-component SPE.

Asymptotic behavior of second-order impulsive differential equations

Directory of Open Access Journals (Sweden)

Haifeng Liu

2011-02-01

Full Text Available In this article, we study the asymptotic behavior of all solutions of 2-th order nonlinear delay differential equation with impulses. Our main tools are impulsive differential inequalities and the Riccati transformation. We illustrate the results by an example.

Supersimetrías y Parasimetrías de la Ecuación de Riccati

OpenAIRE

José Socorro; Marco A. Reyes

2012-01-01

Se describe cómo la ecuación de Riccati nos permite definir de manera coloquial diferentes simetrías en ecuaciones de la física matemática. En particular utilizamos la definición de supersimetría en Mecánica Cuántica, obtenida al desarrollar la segunda solución de Riccati para las ecuaciones del modelo, para introducir la parasimetría de otras ecuaciones de la física matemática y de otras áreas de las ciencias.

Supersimetrías y Parasimetrías de la Ecuación de Riccati

Directory of Open Access Journals (Sweden)

José Socorro

2012-02-01

Full Text Available Se describe cómo la ecuación de Riccati nos permite definir de manera coloquial diferentes simetrías en ecuaciones de la física matemática. En particular utilizamos la definición de supersimetría en Mecánica Cuántica, obtenida al desarrollar la segunda solución de Riccati para las ecuaciones del modelo, para introducir la parasimetría de otras ecuaciones de la física matemática y de otras áreas de las ciencias.

The Extended Fractional Subequation Method for Nonlinear Fractional Differential Equations

OpenAIRE

Zhao, Jianping; Tang, Bo; Kumar, Sunil; Hou, Yanren

2012-01-01

An extended fractional subequation method is proposed for solving fractional differential equations by introducing a new general ansätz and Bäcklund transformation of the fractional Riccati equation with known solutions. Being concise and straightforward, this method is applied to the space-time fractional coupled Burgers’ equations and coupled MKdV equations. As a result, many exact solutions are obtained. It is shown that the considered method provides a very effective, convenient, and powe...

Symbolic computation of exact solutions expressible in rational formal hyperbolic and elliptic functions for nonlinear partial differential equations

International Nuclear Information System (INIS)

Wang Qi; Chen Yong

2007-01-01

With the aid of symbolic computation, some algorithms are presented for the rational expansion methods, which lead to closed-form solutions of nonlinear partial differential equations (PDEs). The new algorithms are given to find exact rational formal polynomial solutions of PDEs in terms of Jacobi elliptic functions, solutions of the Riccati equation and solutions of the generalized Riccati equation. They can be implemented in symbolic computation system Maple. As applications of the methods, we choose some nonlinear PDEs to illustrate the methods. As a result, we not only can successfully obtain the solutions found by most existing Jacobi elliptic function methods and Tanh-methods, but also find other new and more general solutions at the same time

An Innovative Approach to Balancing Chemical-Reaction Equations: A Simplified Matrix-Inversion Technique for Determining The Matrix Null Space

OpenAIRE

Thorne, Lawrence R.

2011-01-01

I propose a novel approach to balancing equations that is applicable to all chemical-reaction equations; it is readily accessible to students via scientific calculators and basic computer spreadsheets that have a matrix-inversion application. The new approach utilizes the familiar matrix-inversion operation in an unfamiliar and innovative way; its purpose is not to identify undetermined coefficients as usual, but, instead, to compute a matrix null space (or matrix kernel). The null space then...

Loop equations for multi-cut matrix models

International Nuclear Information System (INIS)

Akemann, G.

1995-03-01

The loop equation for the complex one-matrix model with a multi-cut structure is derived and solved in the planar limit. An iterative scheme for higher genus contributions to the free energy and the multi-loop correlators is presented for the two-cut model, where explicit results are given up to and including genus two. The double-scaling limit is analyzed and the relation to the one-cut solution of the hermitian and complex one-matrix model is discussed. (orig.)

Oscillation of a class of fractional differential equations with damping term.

Science.gov (United States)

Qin, Huizeng; Zheng, Bin

2013-01-01

We investigate the oscillation of a class of fractional differential equations with damping term. Based on a certain variable transformation, the fractional differential equations are converted into another differential equations of integer order with respect to the new variable. Then, using Riccati transformation, inequality, and integration average technique, some new oscillatory criteria for the equations are established. As for applications, oscillation for two certain fractional differential equations with damping term is investigated by the use of the presented results.

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

Science.gov (United States)

Hanks, Brantley R.; Skelton, Robert E.

1991-01-01

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

P A M Dirac meets M G Krein: matrix orthogonal polynomials and Dirac's equation

International Nuclear Information System (INIS)

Duran, Antonio J; Gruenbaum, F Alberto

2006-01-01

The solution of several instances of the Schroedinger equation (1926) is made possible by using the well-known orthogonal polynomials associated with the names of Hermite, Legendre and Laguerre. A relativistic alternative to this equation was proposed by Dirac (1928) involving differential operators with matrix coefficients. In 1949 Krein developed a theory of matrix-valued orthogonal polynomials without any reference to differential equations. In Duran A J (1997 Matrix inner product having a matrix symmetric second order differential operator Rocky Mt. J. Math. 27 585-600), one of us raised the question of determining instances of these matrix-valued polynomials going along with second order differential operators with matrix coefficients. In Duran A J and Gruenbaum F A (2004 Orthogonal matrix polynomials satisfying second order differential equations Int. Math. Res. Not. 10 461-84), we developed a method to produce such examples and observed that in certain cases there is a connection with the instance of Dirac's equation with a central potential. We observe that the case of the central Coulomb potential discussed in the physics literature in Darwin C G (1928 Proc. R. Soc. A 118 654), Nikiforov A F and Uvarov V B (1988 Special Functions of Mathematical Physics (Basle: Birkhauser) and Rose M E 1961 Relativistic Electron Theory (New York: Wiley)), and its solution, gives rise to a matrix weight function whose orthogonal polynomials solve a second order differential equation. To the best of our knowledge this is the first instance of a connection between the solution of the first order matrix equation of Dirac and the theory of matrix-valued orthogonal polynomials initiated by M G Krein

Novel loop-like solitons for the generalized Vakhnenko equation

International Nuclear Information System (INIS)

Zhang Min; Ma Yu-Lan; Li Bang-Qing

2013-01-01

A non-traveling wave solution of a generalized Vakhnenko equation arising from the high-frequent wave motion in a relaxing medium is derived via the extended Riccati mapping method. The solution includes an arbitrary function of an independent variable. Based on the solution, two hyperbolic functions are chosen to construct new solitons. Novel single-loop-like and double-loop-like solitons are found for the equation

Universal formats for nonlinear ordinary differential systems

International Nuclear Information System (INIS)

Kerner, E.H.

1981-01-01

It is shown that very general nonlinear ordinary differential systems (embracing all that arise in practice) may, first, be brought down to polynomial systems (where the nonlinearities occur only as polynomials in the dependent variables) by introducing suitable new variables into the original system; second, that polynomial systems are reducible to ''Riccati systems,'' where the nonlinearities are quadratic at most; third, that Riccati systems may be brought to elemental universal formats containing purely quadratic terms with simple arrays of coefficients that are all zero or unity. The elemental systems have representations as novel types of matrix Riccati equations. Different starting systems and their associated Riccati systems differ from one another, at the final elemental level, in order and in initial data, but not in format

Matrix form of Legendre polynomials for solving linear integro-differential equations of high order

Science.gov (United States)

Kammuji, M.; Eshkuvatov, Z. K.; Yunus, Arif A. M.

2017-04-01

This paper presents an effective approximate solution of high order of Fredholm-Volterra integro-differential equations (FVIDEs) with boundary condition. Legendre truncated series is used as a basis functions to estimate the unknown function. Matrix operation of Legendre polynomials is used to transform FVIDEs with boundary conditions into matrix equation of Fredholm-Volterra type. Gauss Legendre quadrature formula and collocation method are applied to transfer the matrix equation into system of linear algebraic equations. The latter equation is solved by Gauss elimination method. The accuracy and validity of this method are discussed by solving two numerical examples and comparisons with wavelet and methods.

OSCILLATION OF A SECOND-ORDER HALF-LINEAR NEUTRAL DAMPED DIFFERENTIAL EQUATION WITH TIME-DELAY

Institute of Scientific and Technical Information of China (English)

无

2012-01-01

In this paper,the oscillation for a class of second-order half-linear neutral damped differential equation with time-delay is studied.By means of Yang-inequality,the generalized Riccati transformation and a certain function,some new sufficient conditions for the oscillation are given for all solutions to the equation.

Solution of Moving Boundary Space-Time Fractional Burger’s Equation

Directory of Open Access Journals (Sweden)

E. A.-B. Abdel-Salam

2014-01-01

Full Text Available The fractional Riccati expansion method is used to solve fractional differential equations with variable coefficients. To illustrate the effectiveness of the method, the moving boundary space-time fractional Burger’s equation is studied. The obtained solutions include generalized trigonometric and hyperbolic function solutions. Among these solutions, some are found for the first time. The linear and periodic moving boundaries for the kink solution of the Burger’s equation are presented graphically and discussed.

Matrix Methods for Solving Hartree-Fock Equations in Atomic Structure Calculations and Line Broadening

Directory of Open Access Journals (Sweden)

Thomas Gomez

2018-04-01

Full Text Available Atomic structure of N-electron atoms is often determined by solving the Hartree-Fock equations, which are a set of integro-differential equations. The integral part of the Hartree-Fock equations treats electron exchange, but the Hartree-Fock equations are not often treated as an integro-differential equation. The exchange term is often approximated as an inhomogeneous or an effective potential so that the Hartree-Fock equations become a set of ordinary differential equations (which can be solved using the usual shooting methods. Because the Hartree-Fock equations are an iterative-refinement method, the inhomogeneous term relies on the previous guess of the wavefunction. In addition, there are numerical complications associated with solving inhomogeneous differential equations. This work uses matrix methods to solve the Hartree-Fock equations as an integro-differential equation. It is well known that a derivative operator can be expressed as a matrix made of finite-difference coefficients; energy eigenvalues and eigenvectors can be obtained by using linear-algebra packages. The integral (exchange part of the Hartree-Fock equation can be approximated as a sum and written as a matrix. The Hartree-Fock equations can be solved as a matrix that is the sum of the differential and integral matrices. We compare calculations using this method against experiment and standard atomic structure calculations. This matrix method can also be used to solve for free-electron wavefunctions, thus improving how the atoms and free electrons interact. This technique is important for spectral line broadening in two ways: it improves the atomic structure calculations, and it improves the motion of the plasma electrons that collide with the atom.

Covariant field equations, gauge fields and conservation laws from Yang-Mills matrix models

International Nuclear Information System (INIS)

Steinacker, Harold

2009-01-01

The effective geometry and the gravitational coupling of nonabelian gauge and scalar fields on generic NC branes in Yang-Mills matrix models is determined. Covariant field equations are derived from the basic matrix equations of motions, known as Yang-Mills algebra. Remarkably, the equations of motion for the Poisson structure and for the nonabelian gauge fields follow from a matrix Noether theorem, and are therefore protected from quantum corrections. This provides a transparent derivation and generalization of the effective action governing the SU(n) gauge fields obtained in [1], including the would-be topological term. In particular, the IKKT matrix model is capable of describing 4-dimensional NC space-times with a general effective metric. Metric deformations of flat Moyal-Weyl space are briefly discussed.

Dynamic behavior and chaos control in a complex Riccati-type map ...

African Journals Online (AJOL)

This paper is devoted to analyze the dynamic behavior of a Riccati- type map with complex variables and complex parameters. Fixed points and their asymptotic stability are studied. Lyapunov exponent is computed to indicate chaos. Bifurcation and chaos are discussed. Chaotic behavior of the map has been controlled by ...

Iterative Algorithm for Solving a Class of Quaternion Matrix Equation over the Generalized (P,Q-Reflexive Matrices

Directory of Open Access Journals (Sweden)

Ning Li

2013-01-01

Full Text Available The matrix equation ∑l=1uAlXBl+∑s=1vCsXTDs=F, which includes some frequently investigated matrix equations as its special cases, plays important roles in the system theory. In this paper, we propose an iterative algorithm for solving the quaternion matrix equation ∑l=1uAlXBl+∑s=1vCsXTDs=F over generalized (P,Q-reflexive matrices. The proposed iterative algorithm automatically determines the solvability of the quaternion matrix equation over generalized (P,Q-reflexive matrices. When the matrix equation is consistent over generalized (P,Q-reflexive matrices, the sequence {X(k} generated by the introduced algorithm converges to a generalized (P,Q-reflexive solution of the quaternion matrix equation. And the sequence {X(k} converges to the least Frobenius norm generalized (P,Q-reflexive solution of the quaternion matrix equation when an appropriate initial iterative matrix is chosen. Furthermore, the optimal approximate generalized (P,Q-reflexive solution for a given generalized (P,Q-reflexive matrix X0 can be derived. The numerical results indicate that the iterative algorithm is quite efficient.

PERTURBATION ESTIMATES FOR THE MAXIMAL SOLUTION OF A NONLINEAR MATRIX EQUATION

Directory of Open Access Journals (Sweden)

Vejdi I. Hasanov

2017-06-01

Full Text Available In this paper a nonlinear matrix equation is considered. Perturba- tion estimations for the maximal solution of the considered equation are obtained. The results are illustrated by the use of numerical ex- amples.

CPDES2: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in two dimensions

Science.gov (United States)

Anderson, D. V.; Koniges, A. E.; Shumaker, D. E.

1988-11-01

Many physical problems require the solution of coupled partial differential equations on two-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils and allows for general couplings between all of the component PDE's and it automatically generates the matrix structures needed to perform the algorithm. The resulting sparse matrix equation is solved by either the preconditioned conjugate gradient (CG) method or by the preconditioned biconjugate gradient (BCG) algorithm. An arbitrary number of component equations are permitted only limited by available memory. In the sub-band representation used, we generate an algorithm that is written compactly in terms of indirect indices which is vectorizable on some of the newer scientific computers.

CPDES3: A preconditioned conjugate gradient solver for linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions

Science.gov (United States)

Anderson, D. V.; Koniges, A. E.; Shumaker, D. E.

1988-11-01

Many physical problems require the solution of coupled partial differential equations on three-dimensional domains. When the time scales of interest dictate an implicit discretization of the equations a rather complicated global matrix system needs solution. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximations employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils and allows for general couplings between all of the component PDE's and it automatically generates the matrix structures needed to perform the algorithm. The resulting sparse matrix equation is solved by either the preconditioned conjugate gradient (CG) method or by the preconditioned biconjugate gradient (BCG) algorithm. An arbitrary number of component equations are permitted only limited by available memory. In the sub-band representation used, we generate an algorithm that is written compactly in terms of indirect induces which is vectorizable on some of the newer scientific computers.

Some optimal considerations in attitude control systems. [evaluation of value of relative weighting between time and fuel for relay control law

Science.gov (United States)

Boland, J. S., III

1973-01-01

The conventional six-engine reaction control jet relay attitude control law with deadband is shown to be a good linear approximation to a weighted time-fuel optimal control law. Techniques for evaluating the value of the relative weighting between time and fuel for a particular relay control law is studied along with techniques to interrelate other parameters for the two control laws. Vehicle attitude control laws employing control moment gyros are then investigated. Steering laws obtained from the expression for the reaction torque of the gyro configuration are compared to a total optimal attitude control law that is derived from optimal linear regulator theory. This total optimal attitude control law has computational disadvantages in the solving of the matrix Riccati equation. Several computational algorithms for solving the matrix Riccati equation are investigated with respect to accuracy, computational storage requirements, and computational speed.

Bifurcations of solitary wave solutions for (two and three)-dimensional nonlinear partial differential equation in quantum and magnetized plasma by using two different methods

Science.gov (United States)

Khater, Mostafa M. A.; Seadawy, Aly R.; Lu, Dianchen

2018-06-01

In this research, we study new two techniques that called the extended simple equation method and the novel (G‧/G) -expansion method. The extended simple equation method depend on the auxiliary equation (dϕ/dξ = α + λϕ + μϕ2) which has three ways for solving depends on the specific condition on the parameters as follow: When (λ = 0) this auxiliary equation reduces to Riccati equation, when (α = 0) this auxiliary equation reduces to Bernoulli equation and when (α ≠ 0, λ ≠ 0, μ ≠ 0) we the general solutions of this auxiliary equation while the novel (G‧/G) -expansion method depends also on similar auxiliary equation (G‧/G)‧ = μ + λ(G‧/G) + (v - 1)(G‧/G) 2 which depend also on the value of (λ2 - 4 (v - 1) μ) and the specific condition on the parameters as follow: When (λ = 0) this auxiliary equation reduces to Riccati equation, when (μ = 0) this auxiliary equation reduces to Bernoulli equation and when (λ2 ≠ 4 (v - 1) μ) we the general solutions of this auxiliary equation. This show how both of these auxiliary equation are special cases of Riccati equation. We apply these methods on two dimensional nonlinear Kadomtsev-Petviashvili Burgers equation in quantum plasma and three-dimensional nonlinear modified Zakharov-Kuznetsov equation of ion-acoustic waves in a magnetized plasma. We obtain the exact traveling wave solutions of these important models and under special condition on the parameters, we get solitary traveling wave solutions. All calculations in this study have been established and verified back with the aid of the Maple package program. The executed method is powerful, effective and straightforward for solving nonlinear partial differential equations to obtain more and new solutions.

Bessel equation as an operator identity's matrix element in quantum mechanics

International Nuclear Information System (INIS)

Fan Hongyi; Li Chao

2004-01-01

We study the well-known Bessel equation itself in the framework of quantum mechanics. We show that the Bessel equation is a spontaneous result of an operator identity's matrix element in some definite entangled state representations, which is a fresh look. Application of this operator formalism in the Hankel transform of Laplace equation is presented

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

Oscillation of second order neutral dynamic equations with distributed delay

Directory of Open Access Journals (Sweden)

Qiaoshun Yang

2016-06-01

Full Text Available In this paper, we establish new oscillation criteria for second order neutral dynamic equations with distributed delay by employing the generalized Riccati transformation. The obtained theorems essentially improve the oscillation results in the literature. And two examples are provided to illustrate to the versatility of our main results.

P A M Dirac meets M G Krein: matrix orthogonal polynomials and Dirac's equation

Energy Technology Data Exchange (ETDEWEB)

Duran, Antonio J [Departamento de Analisis Matematico, Universidad de Sevilla, Apdo (PO BOX) 1160, 41080 Sevilla (Spain); Gruenbaum, F Alberto [Department of Mathematics, University of California, Berkeley, CA 94720 (United States)

2006-04-07

The solution of several instances of the Schroedinger equation (1926) is made possible by using the well-known orthogonal polynomials associated with the names of Hermite, Legendre and Laguerre. A relativistic alternative to this equation was proposed by Dirac (1928) involving differential operators with matrix coefficients. In 1949 Krein developed a theory of matrix-valued orthogonal polynomials without any reference to differential equations. In Duran A J (1997 Matrix inner product having a matrix symmetric second order differential operator Rocky Mt. J. Math. 27 585-600), one of us raised the question of determining instances of these matrix-valued polynomials going along with second order differential operators with matrix coefficients. In Duran A J and Gruenbaum F A (2004 Orthogonal matrix polynomials satisfying second order differential equations Int. Math. Res. Not. 10 461-84), we developed a method to produce such examples and observed that in certain cases there is a connection with the instance of Dirac's equation with a central potential. We observe that the case of the central Coulomb potential discussed in the physics literature in Darwin C G (1928 Proc. R. Soc. A 118 654), Nikiforov A F and Uvarov V B (1988 Special Functions of Mathematical Physics (Basle: Birkhauser) and Rose M E 1961 Relativistic Electron Theory (New York: Wiley)), and its solution, gives rise to a matrix weight function whose orthogonal polynomials solve a second order differential equation. To the best of our knowledge this is the first instance of a connection between the solution of the first order matrix equation of Dirac and the theory of matrix-valued orthogonal polynomials initiated by M G Krein.

Exact solution of some linear matrix equations using algebraic methods

Science.gov (United States)

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.

Nonoscillation of half-linear dynamic equations

Czech Academy of Sciences Publication Activity Database

Matucci, S.; Řehák, Pavel

2010-01-01

Roč. 60, č. 5 (2010), s. 1421-1429 ISSN 0898-1221 R&D Projects: GA AV ČR KJB100190701 Grant - others:GA ČR(CZ) GA201/07/0145 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear dynamic equation * time scale * (non)oscillation * Riccati technique Subject RIV: BA - General Mathematics Impact factor: 1.472, year: 2010 http://www.sciencedirect.com/science/article/pii/S0898122110004384

Tetrahedron equations and the relativistic S-matrix of straight-strings in 2+1-dimensions

International Nuclear Information System (INIS)

Zamolodchikov, A.B.

1981-01-01

The quantum S-matrix theory of straight-strings (infinite one-dimensioanl objects like straight domain walls) in 2 + 1-dimensions is considered. The S-matrix is supposed to be purely elastic and factorized. The tetrahedron equations (which are the factorization conditions) are investigated for the special two-colour model. The relativistic three-string S-matrix, which apparently satisfies this tetrahedron equation, is proposed. (orig.)

Comparison of nonlinearities in oscillation theory of half-linear differential equations

Czech Academy of Sciences Publication Activity Database

Řehák, Pavel

2008-01-01

Roč. 121, č. 2 (2008), s. 93-105 ISSN 0236-5294 R&D Projects: GA AV ČR KJB100190701 Institutional research plan: CEZ:AV0Z10190503 Keywords : half-linear differential equation * comparison theorem * Riccati technique Subject RIV: BA - General Mathematics Impact factor: 0.317, year: 2008

The integrability of an extended fifth-order KdV equation with Riccati ...

Indian Academy of Sciences (India)

method was extended to investigate variable coefficient NLEEs, which included gener- alized KdV equation, generalized modified KdV equation and generalized Boussinesq equation [11,12]. It is well known that KdV equation models a variety of nonlinear phenomena, including ion-acoustic waves in plasmas and shallow ...

O(N)-matrix difference equations and a nested Bethe ansatz

International Nuclear Information System (INIS)

Babujian, Hrachya M; Foerster, Angela; Karowski, Michael

2012-01-01

A system of O(N)-matrix difference equations is solved by means of the off-shell version of the nested algebraic Bethe ansatz. In the nesting process, a new object, the Π-matrix, is introduced to overcome the complexities of the O(N)-group structure. The highest weight property of the solutions is proved and some explicit examples are discussed. (paper)

Riccati inequality, disconjugacy, and reciprocity principle for linear Hamiltonian dynamic systems

Czech Academy of Sciences Publication Activity Database

Hilscher, R.; Řehák, Pavel

2003-01-01

Roč. 12, č. 1 (2003), s. 171-189 ISSN 1056-2176 R&D Projects: GA ČR GA201/01/0079; GA ČR GP201/01/P041 Institutional research plan: CEZ:AV0Z1019905; CEZ:AV0Z1019905 Keywords : linear Hamiltonian dynamic systems * disconjugacy * Riccati inequality Subject RIV: BA - General Mathematics Impact factor: 0.256, year: 2002

Oscillation and asymptotic properties of a class of second-order Emden-Fowler neutral differential equations.

Science.gov (United States)

Wang, Rui; Li, Qiqiang

2016-01-01

We consider a class of second-order Emden-Fowler equations with positive and nonpositve neutral coefficients. By using the Riccati transformation and inequalities, several oscillation and asymptotic results are established. Some examples are given to illustrate the main results.

New Solutions of Three Nonlinear Space- and Time-Fractional Partial Differential Equations in Mathematical Physics

International Nuclear Information System (INIS)

Yao Ruo-Xia; Wang Wei; Chen Ting-Hua

2014-01-01

Motivated by the widely used ansätz method and starting from the modified Riemann—Liouville derivative together with a fractional complex transformation that can be utilized to transform nonlinear fractional partial differential equations to nonlinear ordinary differential equations, new types of exact traveling wave solutions to three important nonlinear space- and time-fractional partial differential equations are obtained simultaneously in terms of solutions of a Riccati equation. The results are new and first reported in this paper. (general)

Solution of the scattering T matrix equation in discrete complex momentum space

International Nuclear Information System (INIS)

Rawitscher, G.H.; Delic, G.

1984-01-01

The scattering solution to the Lippmann-Schwinger equation is expanded into a set of spherical Bessel functions of complex wave numbers, K/sub j/, with j = 1,2 , . . . , M. The value of each K/sub j/ is determined from the condition that the spherical Bessel function smoothly matches onto an asymptotically outgoing spherical Hankel (or Coulomb) function of the correct physical wave number at a matching point R. The spherical Bessel functions thus determined are Sturmian functions, and they form a complete set in the interval 0 to R. The coefficients of the expansion of the scattering function are determined by matrix inversion of a linear set of algebraic equations, which are equivalent to the solution of the T-matrix equation in complex momentum space. In view of the presence of a matching radius, no singularities are encountered for the Green's functions, and the inclusion of Coulomb potentials offers no computational difficulties. Three numerical examples are performed in order to illustrate the convergence of the elastic scattering matrix S with M. One of these consists of a set of coupled equations which describe the breakup of a deuteron as it scatters from the nucleus on 58 Ni. A value of M of 15 or less is found sufficient to reproduce the exact S matrix element to an accuracy of four figures after the decimal point

Revolving scheme for solving a cascade of Abel equations in dynamics of planar satellite rotation

Directory of Open Access Journals (Sweden)

Sergey V. Ershkov

2017-05-01

Full Text Available The main objective for this research was the analytical exploration of the dynamics of planar satellite rotation during the motion of an elliptical orbit around a planet. First, we revisit the results of J. Wisdom et al. (1984, in which, by the elegant change of variables (considering the true anomaly f as the independent variable, the governing equation of satellite rotation takes the form of an Abel ordinary differential equation (ODE of the second kind, a sort of generalization of the Riccati ODE. We note that due to the special character of solutions of a Riccati-type ODE, there exists the possibility of sudden jumping in the magnitude of the solution at some moment of time. In the physical sense, this jumping of the Riccati-type solutions of the governing ODE could be associated with the effect of sudden acceleration/deceleration in the satellite rotation around the chosen principle axis at a definite moment of parametric time. This means that there exists not only a chaotic satellite rotation regime (as per the results of J. Wisdom et al. (1984, but a kind of gradient catastrophe (Arnold, 1992 could occur during the satellite rotation process. We especially note that if a gradient catastrophe could occur, this does not mean that it must occur: such a possibility depends on the initial conditions. In addition, we obtained asymptotical solutions that manifest a quasi-periodic character even with the strong simplifying assumptions e→0, p=1, which reduce the governing equation of J. Wisdom et al. (1984 to a kind of Beletskii’s equation.

Preconditioned Krylov and Gauss-Seidel solutions of response matrix equations

International Nuclear Information System (INIS)

Lewis, E.E.; Smith, M.A.; Yang, W.S.; Wollaber, A.

2011-01-01

The use of preconditioned Krylov methods is examined as an alternative to the partitioned matrix acceleration applied to red-black Gauss Seidel (RBGS) iteration that is presently used in the variational nodal code, VARIANT. We employ the GMRES algorithm to treat non-symmetric response matrix equations. A pre conditioner is formulated for the within-group diffusion equation which is equivalent to partitioned matrix acceleration of RBGS iterations. We employ the pre conditioner, which closely parallels two-level p multigrid, to improve RBGS and GMRES algorithms. Of the accelerated algorithms, GMRES converges with less computational effort than RBGS and therefore is chosen for further development. The p multigrid pre conditioner requires response matrices with two or more degrees of freedom (DOF) per interface that are polynomials, which are both orthogonal and hierarchical. It is therefore not directly applicable to very fine mesh calculations that are both slow to converge and that are often modeled with response matrices with only one DOF per interface. Orthogonal matrix aggregation (OMA) is introduced to circumvent this difficulty by combining N×N fine mesh response matrices with one DOF per interface into a coarse mesh response matrix with N orthogonal DOF per interface. Numerical results show that OMA used alone or in combination with p multigrid preconditioning substantially accelerates GMRES solutions. (author)

Preconditioned Krylov and Gauss-Seidel solutions of response matrix equations

Energy Technology Data Exchange (ETDEWEB)

Lewis, E.E., E-mail: e-lewis@northwestern.edu [Department of Mechanical Engineering, Northwestern University, Evanston, IL (United States); Smith, M.A.; Yang, W.S.; Wollaber, A., E-mail: masmith@anl.gov, E-mail: wsyang@anl.gov, E-mail: wollaber@lanl.gov [Nuclear Engineering Division, Argonne National Laboratory, Argonne, IL (United States)

2011-07-01

The use of preconditioned Krylov methods is examined as an alternative to the partitioned matrix acceleration applied to red-black Gauss Seidel (RBGS) iteration that is presently used in the variational nodal code, VARIANT. We employ the GMRES algorithm to treat non-symmetric response matrix equations. A pre conditioner is formulated for the within-group diffusion equation which is equivalent to partitioned matrix acceleration of RBGS iterations. We employ the pre conditioner, which closely parallels two-level p multigrid, to improve RBGS and GMRES algorithms. Of the accelerated algorithms, GMRES converges with less computational effort than RBGS and therefore is chosen for further development. The p multigrid pre conditioner requires response matrices with two or more degrees of freedom (DOF) per interface that are polynomials, which are both orthogonal and hierarchical. It is therefore not directly applicable to very fine mesh calculations that are both slow to converge and that are often modeled with response matrices with only one DOF per interface. Orthogonal matrix aggregation (OMA) is introduced to circumvent this difficulty by combining N×N fine mesh response matrices with one DOF per interface into a coarse mesh response matrix with N orthogonal DOF per interface. Numerical results show that OMA used alone or in combination with p multigrid preconditioning substantially accelerates GMRES solutions. (author)

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.

Darboux transformations for (1+2)-dimensional Fokker-Planck equations with constant diffusion matrix

International Nuclear Information System (INIS)

Schulze-Halberg, Axel

2012-01-01

We construct a Darboux transformation for (1+2)-dimensional Fokker-Planck equations with constant diffusion matrix. Our transformation is based on the two-dimensional supersymmetry formalism for the Schrödinger equation. The transformed Fokker-Planck equation and its solutions are obtained in explicit form.

Low-temperature random matrix theory at the soft edge

International Nuclear Information System (INIS)

Edelman, Alan; Persson, Per-Olof; Sutton, Brian D.

2014-01-01

“Low temperature” random matrix theory is the study of random eigenvalues as energy is removed. In standard notation, β is identified with inverse temperature, and low temperatures are achieved through the limit β → ∞. In this paper, we derive statistics for low-temperature random matrices at the “soft edge,” which describes the extreme eigenvalues for many random matrix distributions. Specifically, new asymptotics are found for the expected value and standard deviation of the general-β Tracy-Widom distribution. The new techniques utilize beta ensembles, stochastic differential operators, and Riccati diffusions. The asymptotics fit known high-temperature statistics curiously well and contribute to the larger program of general-β random matrix theory

On the reduction of the multidimensional stationary Schroedinger equation to a first-order equation and its relation to the pseudoanalytic function theory

Energy Technology Data Exchange (ETDEWEB)

Kravchenko, Vladislav V [Departmento de Telecomunicaciones, SEPI, Escuela Superior de IngenierIa Mecanica y Electrica, Instituto Politecnico Nacional, CP 07738 Mexico DF (Mexico)

2005-01-28

Given a particular solution of a one-dimensional stationary Schroedinger equation this equation of second order can be reduced to a first-order linear ordinary differential equation. This is done with the aid of an auxiliary Riccati differential equation. In the present work we show that the same fact is true in a multidimensional situation also. For simplicity we consider the case of two or three independent variables. One particular solution of the stationary Schroedinger equation allows us to reduce this second-order equation to a linear first-order quaternionic differential equation. As in the one-dimensional case this is done with the aid of an auxiliary quaternionic Riccati equation. The resulting first-order quaternionic equation is equivalent to the static Maxwell system and is closely related to the Dirac equation. In the case of two independent variables it is the well-known Vekua equation from theory of pseudoanalytic (or generalized analytic) functions. Nevertheless, we show that even in this case it is very useful to consider not only complex valued functions, solutions of the Vekua equation, but complete quaternionic functions. In this way the first-order quaternionic equation represents two separate Vekua equations, one of which gives us solutions of the Schroedinger equation and the other one can be considered as an auxiliary equation of a simpler structure. Moreover for the auxiliary equation we always have the corresponding Bers generating pair (F, G), the base of the Bers theory of pseudoanalytic functions, and what is very important, the Bers derivatives of solutions of the auxiliary equation give us solutions of the main Vekua equation and as a consequence of the Schroedinger equation. Based on this fact we obtain an analogue of the Cauchy integral theorem for solutions of the stationary Schroedinger equation. Other results from theory of pseudoanalytic functions can be written for solutions of the Schroedinger equation. Moreover, for an ample

New exact solutions to MKDV-Burgers equation and (2 + 1)-dimensional dispersive long wave equation via extended Riccati equation method

International Nuclear Information System (INIS)

Kong Cuicui; Wang Dan; Song Lina; Zhang Hongqing

2009-01-01

In this paper, with the aid of symbolic computation and a general ansaetz, we presented a new extended rational expansion method to construct new rational formal exact solutions to nonlinear partial differential equations. In order to illustrate the effectiveness of this method, we apply it to the MKDV-Burgers equation and the (2 + 1)-dimensional dispersive long wave equation, then several new kinds of exact solutions are successfully obtained by using the new ansaetz. The method can also be applied to other nonlinear partial differential equations.

A Riccati-Based Interior Point Method for Efficient Model Predictive Control of SISO Systems

DEFF Research Database (Denmark)

Hagdrup, Morten; Johansson, Rolf; Bagterp Jørgensen, John

2017-01-01

model parts separate. The controller is designed based on the deterministic model, while the Kalman filter results from the stochastic part. The controller is implemented as a primal-dual interior point (IP) method using Riccati recursion and the computational savings possible for SISO systems...

Numerical solution of quadratic matrix equations for free vibration analysis of structures

Science.gov (United States)

Gupta, K. K.

1975-01-01

This paper is concerned with the efficient and accurate solution of the eigenvalue problem represented by quadratic matrix equations. Such matrix forms are obtained in connection with the free vibration analysis of structures, discretized by finite 'dynamic' elements, resulting in frequency-dependent stiffness and inertia matrices. The paper presents a new numerical solution procedure of the quadratic matrix equations, based on a combined Sturm sequence and inverse iteration technique enabling economical and accurate determination of a few required eigenvalues and associated vectors. An alternative procedure based on a simultaneous iteration procedure is also described when only the first few modes are the usual requirement. The employment of finite dynamic elements in conjunction with the presently developed eigenvalue routines results in a most significant economy in the dynamic analysis of structures.

One-parameter families of supersymmetric isospectral potentials from Riccati solutions in function composition form

Energy Technology Data Exchange (ETDEWEB)

Rosu, Haret C., E-mail: hcr@ipicyt.edu.mx [IPICYT, Instituto Potosino de Investigacion Cientifica y Tecnologica, Camino a la presa San José 2055, Col. Lomas 4a Sección, 78216 San Luis Potosí, S.L.P. (Mexico); Mancas, Stefan C., E-mail: mancass@erau.edu [Department of Mathematics, Embry–Riddle Aeronautical University, Daytona Beach, FL 32114-3900 (United States); Chen, Pisin, E-mail: pisinchen@phys.ntu.edu.tw [Leung Center for Cosmology and Particle Astrophysics (LeCosPA) and Department of Physics, National Taiwan University, Taipei 10617, Taiwan (China)

2014-04-15

In the context of supersymmetric quantum mechanics, we define a potential through a particular Riccati solution of the composition form (F∘f)(x)=F(f(x)) and obtain a generalized Mielnik construction of one-parameter isospectral potentials when we use the general Riccati solution. Some examples for special cases of F and f are given to illustrate the method. An interesting result is obtained in the case of a parametric double well potential generated by this method, for which it is shown that the parameter of the potential controls the heights of the localization probability in the two wells, and for certain values of the parameter the height of the localization probability can be higher in the smaller well. -- Highlights: •Function-composition generalization of parametric isospectral potentials is presented. •Mielnik one-parameter family of harmonic potentials is obtained as a particular case. •Graphical discussion of regular and singular regions in the parameter space is given.

Factorizations of rational matrix functions with application to discrete isomonodromic transformations and difference Painleve equations

International Nuclear Information System (INIS)

Dzhamay, Anton

2009-01-01

We study factorizations of rational matrix functions with simple poles on the Riemann sphere. For the quadratic case (two poles) we show, using multiplicative representations of such matrix functions, that a good coordinate system on this space is given by a mix of residue eigenvectors of the matrix and its inverse. Our approach is motivated by the theory of discrete isomonodromic transformations and their relationship with difference Painleve equations. In particular, in these coordinates, basic isomonodromic transformations take the form of the discrete Euler-Lagrange equations. Secondly we show that dPV equations, previously obtained in this context by D Arinkin and A Borodin, can be understood as simple relationships between the residues of such matrices and their inverses.

The modified alternative (G'/G)-expansion method to nonlinear evolution equation: application to the (1+1)-dimensional Drinfel'd-Sokolov-Wilson equation.

Science.gov (United States)

Akbar, M Ali; Mohd Ali, Norhashidah Hj; Mohyud-Din, Syed Tauseef

2013-01-01

Over the years, (G'/G)-expansion method is employed to generate traveling wave solutions to various wave equations in mathematical physics. In the present paper, the alternative (G'/G)-expansion method has been further modified by introducing the generalized Riccati equation to construct new exact solutions. In order to illustrate the novelty and advantages of this approach, the (1+1)-dimensional Drinfel'd-Sokolov-Wilson (DSW) equation is considered and abundant new exact traveling wave solutions are obtained in a uniform way. These solutions may be imperative and significant for the explanation of some practical physical phenomena. It is shown that the modified alternative (G'/G)-expansion method an efficient and advance mathematical tool for solving nonlinear partial differential equations in mathematical physics.

New Schemes for Positive Real Truncation

Directory of Open Access Journals (Sweden)

Kari Unneland

2007-07-01

Full Text Available Model reduction, based on balanced truncation, of stable and of positive real systems are considered. An overview over some of the already existing techniques are given: Lyapunov balancing and stochastic balancing, which includes Riccati balancing. A novel scheme for positive real balanced truncation is then proposed, which is a combination of the already existing Lyapunov balancing and Riccati balancing. Using Riccati balancing, the solution of two Riccati equations are needed to obtain positive real reduced order systems. For the suggested method, only one Lyapunov equation and one Riccati equation are solved in order to obtain positive real reduced order systems, which is less computationally demanding. Further it is shown, that in order to get positive real reduced order systems, only one Riccati equation needs to be solved. Finally, this is used to obtain positive real frequency weighted balanced truncation.

Lyapunov Functions and Solutions of the Lyapunov Matrix Equation for Marginally Stable Systems

DEFF Research Database (Denmark)

Kliem, Wolfhard; Pommer, Christian

2000-01-01

We consider linear systems of differential equations $I \\ddot{x}+B \\dot{x}+C{x}={0}$ where $I$ is the identity matrix and $B$ and $C$ are general complex $n$ x $n$ matrices. Our main interest is to determine conditions for complete marginalstability of these systems. To this end we find solutions...... of the Lyapunov matrix equation and characterize the set of matrices $(B, C)$ which guarantees marginal stability. The theory is applied to gyroscopic systems, to indefinite damped systems, and to circulatory systems, showing how to choose certain parameter matrices to get sufficient conditions for marginal...... stability.Comparison is made with some known results for equations with real system matrices.Moreover more general cases are investigated and several examples are given....

Matrix-type multiple reciprocity boundary element method for solving three-dimensional two-group neutron diffusion equations

International Nuclear Information System (INIS)

Itagaki, Masafumi; Sahashi, Naoki.

1997-01-01

The multiple reciprocity boundary element method has been applied to three-dimensional two-group neutron diffusion problems. A matrix-type boundary integral equation has been derived to solve the first and the second group neutron diffusion equations simultaneously. The matrix-type fundamental solutions used here satisfy the equation which has a point source term and is adjoint to the neutron diffusion equations. A multiple reciprocity method has been employed to transform the matrix-type domain integral related to the fission source into an equivalent boundary one. The higher order fundamental solutions required for this formulation are composed of a series of two types of analytic functions. The eigenvalue itself is also calculated using only boundary integrals. Three-dimensional test calculations indicate that the present method provides stable and accurate solutions for criticality problems. (author)

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

The solution space of the unitary matrix model string equation and the Sato Grassmannian

International Nuclear Information System (INIS)

Anagnostopoulos, K.N.; Bowick, M.J.; Schwarz, A.

1992-01-01

The space of all solutions to the string equation of the symmetric unitary one-matrix model is determined. It is shown that the string equations is equivalent to simple conditions on points V 1 and V 2 in the big cell Gr (0) of the Sato Grassmannian Gr. This is a consequence of a well-defined continuum limit in which the string equation has the simple form [P, 2 - ]=1, with P and 2 - 2x2 matrices of differential operators. These conditions on V 1 and V 2 yield a simple system of first order differential equations whose analysis determines the space of all solutions to the string equation. This geometric formulation leads directly to the Virasoro constraints L n (n≥0), where L n annihilate the two modified-KdV τ-functions whose product gives the partition function of the Unitary Matrix Model. (orig.)

Polynomial two-parameter eigenvalue problems and matrix pencil methods for stability of delay-differential equations

NARCIS (Netherlands)

Jarlebring, E.; Hochstenbach, M.E.

2009-01-01

Several recent methods used to analyze asymptotic stability of delay-differential equations (DDEs) involve determining the eigenvalues of a matrix, a matrix pencil or a matrix polynomial constructed by Kronecker products. Despite some similarities between the different types of these so-called

Numerical Solution of Multiterm Fractional Differential Equations Using the Matrix Mittag–Leffler Functions

Directory of Open Access Journals (Sweden)

Marina Popolizio

2018-01-01

Full Text Available Multiterm fractional differential equations (MTFDEs nowadays represent a widely used tool to model many important processes, particularly for multirate systems. Their numerical solution is then a compelling subject that deserves great attention, not least because of the difficulties to apply general purpose methods for fractional differential equations (FDEs to this case. In this paper, we first transform the MTFDEs into equivalent systems of FDEs, as done by Diethelm and Ford; in this way, the solution can be expressed in terms of Mittag–Leffler (ML functions evaluated at matrix arguments. We then propose to compute it by resorting to the matrix approach proposed by Garrappa and Popolizio. Several numerical tests are presented that clearly show that this matrix approach is very accurate and fast, also in comparison with other numerical methods.

Computation of rational solutions for a first-order nonlinear differential equation

Directory of Open Access Journals (Sweden)

Djilali Behloul

2011-09-01

Full Text Available In this article, we study differential equations of the form $y'=sum A_i(xy^i/sum B_i(xy^i$ which can be elliptic, hyperbolic, parabolic, Riccati, or quasi-linear. We show how rational solutions can be computed in a systematic manner. Such results are most likely to find applications in the theory of limit cycles as indicated by Gine et al [4].

Monodromy of the matrix Schroedinger equations and Darboux transformations

CERN Document Server

Goncharenko, V M

1998-01-01

A Schroedinger operator L=-d sup 2 /dz sup 2 +U(z) with a matrix-valued rational potential U(z) is said to have trivial monodromy if all the solutions of the corresponding Schroedinger equations L psi=lambda psi are single-valued in the complex plane z is an element of C for any lambda. A local criterion of this property in terms of the Laurent coefficients of the potential U near its singularities, which are assumed to be regular, is found. It is proved that any such operator with a potential vanishing at infinity can be obtained by a matrix analogue of the Darboux transformation from the Schroedinger operator L sub o =-d sup 2 /dz sup 2. This generalizes the well known Duistermaat-Gruenbaum result to the matrix case and gives the explicit description of the Schroedinger operators with trivial monodromy in this case. (author)

Alpha particle spectroscopy for CR-39 detector utilizing matrix of energy equations

Energy Technology Data Exchange (ETDEWEB)

Awad, E.M. [Department of General Sciences, Yanbu Industrial College, PO Box 30436, Madinat Yanbu Al-Sinaiya (Saudi Arabia); Physics Department, Faculty of Science, Menofia University, Shebin El-Koom (Egypt)], E-mail: ayawad@yahoo.com; Soliman, A.A. [Department of Mathematics, Faculty of Education (AL-Arish), Suez Canal University, AL-Arish 45111 (Egypt); Department of Mathematics, Teacher' s College (Bisha), King Khalid University, Bisha, PO Box 551 (Saudi Arabia)], E-mail: asoliman_99@yahoo.com; Rammah, Y.S. [Physics Department, Faculty of Science, Menofia University, Shebin El-Koom (Egypt)

2007-10-01

A method for determining alpha-particle energy using CR-39 detector by utilizing matrix of energy equation was described. The matrix was composed from two axes; the track minor axis (m) and diameter of etched out track end (d) axis of some selected elliptical tracks. The energy E in (m,d) coordinate was approximated by matrix of energy equations given by: E{sub k}={sigma}{sub i,j=0}{sup 2}a{sub ij}d{sub k}{sup i}m{sub k}{sup j}, which was identified using two different approaches. First, i and j were treated as power exponents for d and m. The adjusting parameters values a{sub ij} were obtained and the energy of a given track was deduced directly from it. Second, i and j were treated as indices of some chosen tracks that were fitted to obtain iso-energy curves that were superimposed on m-d scatter plot as calibration curves. The energy between any two successive iso-energy curves in this case was assumed varied linearly with d for a given m. The energy matrix in both cases was solved numerically. Results of the two approaches were compared.

Oscillation and non-oscillation criterion for Riemann–Weber type half-linear differential equations

Directory of Open Access Journals (Sweden)

Petr Hasil

2016-08-01

Full Text Available By the combination of the modified half-linear Prüfer method and the Riccati technique, we study oscillatory properties of half-linear differential equations. Taking into account the transformation theory of half-linear equations and using some known results, we show that the analysed equations in the Riemann–Weber form with perturbations in both terms are conditionally oscillatory. Within the process, we identify the critical oscillation values of their coefficients and, consequently, we decide when the considered equations are oscillatory and when they are non-oscillatory. As a direct corollary of our main result, we solve the so-called critical case for a certain type of half-linear non-perturbed equations.

An Operational Matrix Technique for Solving Variable Order Fractional Differential-Integral Equation Based on the Second Kind of Chebyshev Polynomials

Directory of Open Access Journals (Sweden)

Jianping Liu

2016-01-01

Full Text Available An operational matrix technique is proposed to solve variable order fractional differential-integral equation based on the second kind of Chebyshev polynomials in this paper. The differential operational matrix and integral operational matrix are derived based on the second kind of Chebyshev polynomials. Using two types of operational matrixes, the original equation is transformed into the arithmetic product of several dependent matrixes, which can be viewed as an algebraic system after adopting the collocation points. Further, numerical solution of original equation is obtained by solving the algebraic system. Finally, several examples show that the numerical algorithm is computationally efficient.

A Globally Convergent Matrix-Free Method for Constrained Equations and Its Linear Convergence Rate

Directory of Open Access Journals (Sweden)

Min Sun

2014-01-01

Full Text Available A matrix-free method for constrained equations is proposed, which is a combination of the well-known PRP (Polak-Ribière-Polyak conjugate gradient method and the famous hyperplane projection method. The new method is not only derivative-free, but also completely matrix-free, and consequently, it can be applied to solve large-scale constrained equations. We obtain global convergence of the new method without any differentiability requirement on the constrained equations. Compared with the existing gradient methods for solving such problem, the new method possesses linear convergence rate under standard conditions, and a relax factor γ is attached in the update step to accelerate convergence. Preliminary numerical results show that it is promising in practice.

A high-performance Riccati based solver for tree-structured quadratic programs

DEFF Research Database (Denmark)

Frison, Gianluca; Kouzoupis, Dimitris; Diehl, Moritz

2017-01-01

the online solution of such problems challenging and the development of tailored solvers crucial. In this paper, an interior point method is presented that can solve Quadratic Programs (QPs) arising in multi-stage MPC efficiently by means of a tree-structured Riccati recursion and a high-performance linear...... algebra library. A performance comparison with code-generated and general purpose sparse QP solvers shows that the computation times can be significantly reduced for all problem sizes that are practically relevant in embedded MPC applications. The presented implementation is freely available as part...

Nonlocal symmetries, solitary waves and cnoidal periodic waves of the (2+1)-dimensional breaking soliton equation

Science.gov (United States)

Zou, Li; Tian, Shou-Fu; Feng, Lian-Li

2017-12-01

In this paper, we consider the (2+1)-dimensional breaking soliton equation, which describes the interaction of a Riemann wave propagating along the y-axis with a long wave along the x-axis. By virtue of the truncated Painlevé expansion method, we obtain the nonlocal symmetry, Bäcklund transformation and Schwarzian form of the equation. Furthermore, by using the consistent Riccati expansion (CRE), we prove that the breaking soliton equation is solvable. Based on the consistent tan-function expansion, we explicitly derive the interaction solutions between solitary waves and cnoidal periodic waves.

A Riccati Based Homogeneous and Self-Dual Interior-Point Method for Linear Economic Model Predictive Control

DEFF Research Database (Denmark)

Sokoler, Leo Emil; Frison, Gianluca; Edlund, Kristian

2013-01-01

In this paper, we develop an efficient interior-point method (IPM) for the linear programs arising in economic model predictive control of linear systems. The novelty of our algorithm is that it combines a homogeneous and self-dual model, and a specialized Riccati iteration procedure. We test...

Massively parallel red-black algorithms for x-y-z response matrix equations

International Nuclear Information System (INIS)

Hanebutte, U.R.; Laurin-Kovitz, K.; Lewis, E.E.

1992-01-01

Recently, both discrete ordinates and spherical harmonic (S n and P n ) methods have been cast in the form of response matrices. In x-y geometry, massively parallel algorithms have been developed to solve the resulting response matrix equations on the Connection Machine family of parallel computers, the CM-2, CM-200, and CM-5. These algorithms utilize two-cycle iteration on a red-black checkerboard. In this work we examine the use of massively parallel red-black algorithms to solve response matric equations in three dimensions. This longer term objective is to utilize massively parallel algorithms to solve S n and/or P n response matrix problems. In this exploratory examination, however, we consider the simple 6 x 6 response matrices that are derivable from fine-mesh diffusion approximations in three dimensions

S-matrix approach to the equation of state of dilute nuclear matter

Indian Academy of Sciences (India)

2014-04-01

matrix framework, a method is presented to calculate the equation of state of dilute warm nuclear matter. The result is a model-independent virial series for the pressure and density that systematically includes contributions from ...

Solution of the Multigroup-Diffusion equation by the response matrix method

International Nuclear Information System (INIS)

Oliveira, C.R.E.

1980-10-01

A preliminary analysis of the response matrix method is made, considering its application to the solution of the multigroup diffusion equations. The one-dimensional formulation is presented and used to test some flux expansions, seeking the application of the method to the two-dimensional problem. This formulation also solves the equations that arise from the integro-differential synthesis algorithm. The slow convergence of the power method, used to solve the eigenvalue problem, and its acceleration by means of the Chebyshev polynomial method, are also studied. An algorithm for the estimation of the dominance ratio is presented, based on the residues of two successive iteration vectors. This ratio, which is not known a priori, is fundamental for the efficiency of the method. Some numerical problems are solved, testing the 1D formulation of the response matrix method, its application to the synthesis algorithm and also, at the same time, the algorithm to accelerate the source problem. (Author) [pt

The ESS and replicator equation in matrix games under time constraints.

Science.gov (United States)

Garay, József; Cressman, Ross; Móri, Tamás F; Varga, Tamás

2018-06-01

Recently, we introduced the class of matrix games under time constraints and characterized the concept of (monomorphic) evolutionarily stable strategy (ESS) in them. We are now interested in how the ESS is related to the existence and stability of equilibria for polymorphic populations. We point out that, although the ESS may no longer be a polymorphic equilibrium, there is a connection between them. Specifically, the polymorphic state at which the average strategy of the active individuals in the population is equal to the ESS is an equilibrium of the polymorphic model. Moreover, in the case when there are only two pure strategies, a polymorphic equilibrium is locally asymptotically stable under the replicator equation for the pure-strategy polymorphic model if and only if it corresponds to an ESS. Finally, we prove that a strict Nash equilibrium is a pure-strategy ESS that is a locally asymptotically stable equilibrium of the replicator equation in n-strategy time-constrained matrix games.

Boltzmann’s Six-Moment One-Dimensional Nonlinear System Equations with the Maxwell-Auzhan Boundary Conditions

Directory of Open Access Journals (Sweden)

A. Sakabekov

2016-01-01

Full Text Available We prove existence and uniqueness of the solution of the problem with initial and Maxwell-Auzhan boundary conditions for nonstationary nonlinear one-dimensional Boltzmann’s six-moment system equations in space of functions continuous in time and summable in square by a spatial variable. In order to obtain a priori estimation of the initial and boundary value problem for nonstationary nonlinear one-dimensional Boltzmann’s six-moment system equations we get the integral equality and then use the spherical representation of vector. Then we obtain the initial value problem for Riccati equation. We have managed to obtain a particular solution of this equation in an explicit form.

Solving Eigenvalue response matrix equations with Jacobian-Free Newton-Krylov methods

International Nuclear Information System (INIS)

Roberts, Jeremy A.; Forget, Benoit

2011-01-01

The response matrix method for reactor eigenvalue problems is motivated as a technique for solving coarse mesh transport equations, and the classical approach of power iteration (PI) for solution is described. The method is then reformulated as a nonlinear system of equations, and the associated Jacobian is derived. A Jacobian-Free Newton-Krylov (JFNK) method is employed to solve the system, using an approximate Jacobian coupled with incomplete factorization as a preconditioner. The unpreconditioned JFNK slightly outperforms PI, and preconditioned JFNK outperforms both PI and Steffensen-accelerated PI significantly. (author)

Parallel Implementation of Riccati Recursion for Solving Linear-Quadratic Control Problems

DEFF Research Database (Denmark)

Frison, Gianluca; Jørgensen, John Bagterp

2013-01-01

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...... an alternative version of the Riccati recursion solver for LQ control problems is presented. The performance of both the classical and the alternative version is analyzed from a theoretical as well as a numerical point of view, and the alternative version is found to be approximately 50% faster than...

Newton's method for solving a quadratic matrix equation with special coefficient matrices

International Nuclear Information System (INIS)

Seo, Sang-Hyup; Seo, Jong Hyun; Kim, Hyun-Min

2014-01-01

We consider the iterative method for solving a quadratic matrix equation with special coefficient matrices which arises in the quasi-birth-death problem. In this paper, we show that the elementwise minimal positive solvents to quadratic matrix equations can be obtained using Newton's method. We also prove that the convergence rate of the Newton iteration is quadratic if the Fréchet derivative at the elementwise minimal positive solvent is nonsingular. However, if the Fréchet derivative is singular, the convergence rate is at least linear. Numerical experiments of the convergence rate are given.(This is summarized a paper which is to appear in Honam Mathematical Journal.)

Soliton and periodic solutions for higher order wave equations of KdV type (I)

International Nuclear Information System (INIS)

Khuri, S.A.

2005-01-01

The aim of the paper is twofold. First, a new ansaetze is introduced for the construction of exact solutions for higher order wave equations of KdV type (I). We show the existence of a class of discontinuous soliton solutions with infinite spikes. Second, the projective Riccati technique is implemented as an alternate approach for obtaining new exact solutions, solitary solutions, and periodic wave solutions

A perturbative solution to metadynamics ordinary differential equation

Science.gov (United States)

Tiwary, Pratyush; Dama, James F.; Parrinello, Michele

2015-12-01

Metadynamics is a popular enhanced sampling scheme wherein by periodic application of a repulsive bias, one can surmount high free energy barriers and explore complex landscapes. Recently, metadynamics was shown to be mathematically well founded, in the sense that the biasing procedure is guaranteed to converge to the true free energy surface in the long time limit irrespective of the precise choice of biasing parameters. A differential equation governing the post-transient convergence behavior of metadynamics was also derived. In this short communication, we revisit this differential equation, expressing it in a convenient and elegant Riccati-like form. A perturbative solution scheme is then developed for solving this differential equation, which is valid for any generic biasing kernel. The solution clearly demonstrates the robustness of metadynamics to choice of biasing parameters and gives further confidence in the widely used method.

A perturbative solution to metadynamics ordinary differential equation.

Science.gov (United States)

Tiwary, Pratyush; Dama, James F; Parrinello, Michele

2015-12-21

Metadynamics is a popular enhanced sampling scheme wherein by periodic application of a repulsive bias, one can surmount high free energy barriers and explore complex landscapes. Recently, metadynamics was shown to be mathematically well founded, in the sense that the biasing procedure is guaranteed to converge to the true free energy surface in the long time limit irrespective of the precise choice of biasing parameters. A differential equation governing the post-transient convergence behavior of metadynamics was also derived. In this short communication, we revisit this differential equation, expressing it in a convenient and elegant Riccati-like form. A perturbative solution scheme is then developed for solving this differential equation, which is valid for any generic biasing kernel. The solution clearly demonstrates the robustness of metadynamics to choice of biasing parameters and gives further confidence in the widely used method.

Matrix Solution of Coupled Differential Equations and Looped Car Following Models

Science.gov (United States)

McCartney, Mark

2008-01-01

A simple mathematical model for the behaviour of how vehicles follow each other along a looped stretch of road is described. The resulting coupled first order differential equations are solved using appropriate matrix techniques and the physical significance of the model is discussed. A number possible classroom exercises are suggested to help…

On the solution of Stein's equation and Fisher information matrix of an ARMAX process

NARCIS (Netherlands)

Klein, A.; Spreij, P.

2004-01-01

The main goal of this paper consists in expressing the solution of a Stein equation in terms of the Fisher information matrix (FIM) of a scalar ARMAX process. A condition for expressing the FIM in terms of a solution to a Stein equation is also set forth. Such interconnections can be derived when a

Delay differential equations via the matrix Lambert W function and bifurcation analysis: application to machine tool chatter.

Science.gov (United States)

Yi, Sun; Nelson, Patrick W; Ulsoy, A Galip

2007-04-01

In a turning process modeled using delay differential equations (DDEs), we investigate the stability of the regenerative machine tool chatter problem. An approach using the matrix Lambert W function for the analytical solution to systems of delay differential equations is applied to this problem and compared with the result obtained using a bifurcation analysis. The Lambert W function, known to be useful for solving scalar first-order DDEs, has recently been extended to a matrix Lambert W function approach to solve systems of DDEs. The essential advantages of the matrix Lambert W approach are not only the similarity to the concept of the state transition matrix in lin ear ordinary differential equations, enabling its use for general classes of linear delay differential equations, but also the observation that we need only the principal branch among an infinite number of roots to determine the stability of a system of DDEs. The bifurcation method combined with Sturm sequences provides an algorithm for determining the stability of DDEs without restrictive geometric analysis. With this approach, one can obtain the critical values of delay, which determine the stability of a system and hence the preferred operating spindle speed without chatter. We apply both the matrix Lambert W function and the bifurcation analysis approach to the problem of chatter stability in turning, and compare the results obtained to existing methods. The two new approaches show excellent accuracy and certain other advantages, when compared to traditional graphical, computational and approximate methods.

Explicit solutions of the cubic matrix nonlinear Schrödinger equation

International Nuclear Information System (INIS)

Demontis, Francesco; Mee, Cornelis van der

2008-01-01

In this paper, we derive a class of explicit solutions, global in (x, t) is an element of R 2 , of the focusing matrix nonlinear Schrödinger equation using straightforward linear algebra. We obtain both the usual and multiple pole multisoliton solutions as well as a new class of solutions exponentially decaying as x → ±∞

Loop equations and topological recursion for the arbitrary-$\\beta$ two-matrix model

CERN Document Server

Bergère, Michel; Marchal, Olivier; Prats-Ferrer, Aleix

2012-01-01

We write the loop equations for the $\\beta$ two-matrix model, and we propose a topological recursion algorithm to solve them, order by order in a small parameter. We find that to leading order, the spectral curve is a "quantum" spectral curve, i.e. it is given by a differential operator (instead of an algebraic equation for the hermitian case). Here, we study the case where that quantum spectral curve is completely degenerate, it satisfies a Bethe ansatz, and the spectral curve is the Baxter TQ relation.

Direct Yaw Control of Vehicle using State Dependent Riccati Equation with Integral Terms

Directory of Open Access Journals (Sweden)

SANDHU, F.

2016-05-01

Full Text Available Direct yaw control of four-wheel vehicles using optimal controllers such as the linear quadratic regulator (LQR and the sliding mode controller (SMC either considers only certain parameters constant in the nonlinear equations of vehicle model or totally neglect their effects to obtain simplified models, resulting in loss of states for the system. In this paper, a modified state-dependent Ricatti equation method obtained by the simplification of the vehicle model is proposed. This method overcomes the problem of the lost states by including state integrals. The results of the proposed system are compared with the sliding mode slip controller and state-dependent Ricatti equation method using high fidelity vehicle model in the vehicle simulation software package, Carsim. Results show 38% reduction in the lateral velocity, 34% reduction in roll and 16% reduction in excessive yaw by only increasing the fuel consumption by 6.07%.

On the solution of a rational matrix equation arising in G-networks

NARCIS (Netherlands)

B. Meini (Beatrice); T. Nesti (Tommaso)

2017-01-01

textabstractWe consider the problem of solving a rational matrix equation arising in the solution of G-networks. We propose and analyze two numerical methods: a fixed point iteration and the Newton–Raphson method. The fixed point iteration is shown to be globally convergent with linear convergence

Uniformly constructing combinatorial solutions, combining a rational function with hyperbolic or trigonometric functions, for the (2+1) dimensional Broer-Kaup-Kupershmidt equation

International Nuclear Information System (INIS)

Liu Qing; Wang Zihua

2010-01-01

According to two dependent rational solutions to a generalized Riccati equation together with the equation itself, a rational-exponent solution to a nonlinear partial differential equation can be constructed. By selecting different parameter values in the rational-exponent solution, many families of combinatorial solutions combined with a rational function such as hyperbolic functions or trigonometric functions, are rapidly derived. This method is applied to the Whitham-Broer-Kaup equation and a series of combinatorial solutions are obtained, showing that this method is a more concise and efficient approach and can uniformly construct many types of combined solutions to nonlinear partial differential equations.

Some Matrix Iterations for Computing Generalized Inverses and Balancing Chemical Equations

Directory of Open Access Journals (Sweden)

Farahnaz Soleimani

2015-11-01

Full Text Available An application of iterative methods for computing the Moore–Penrose inverse in balancing chemical equations is considered. With the aim to illustrate proposed algorithms, an improved high order hyper-power matrix iterative method for computing generalized inverses is introduced and applied. The improvements of the hyper-power iterative scheme are based on its proper factorization, as well as on the possibility to accelerate the iterations in the initial phase of the convergence. Although the effectiveness of our approach is confirmed on the basis of the theoretical point of view, some numerical comparisons in balancing chemical equations, as well as on randomly-generated matrices are furnished.

New Jacobian Matrix and Equations of Motion for a 6 d.o.f Cable-Driven Robot

Directory of Open Access Journals (Sweden)

Ali Afshari

2007-03-01

Full Text Available In this paper, we introduce a new method and new motion variables to study kinematics and dynamics of a 6 d.o.f cable-driven robot. Using these new variables and Lagrange equations, we achieve new equations of motion which are different in appearance and several aspects from conventional equations usually used to study 6 d.o.f cable robots. Then, we introduce a new Jacobian matrix which expresses kinematical relations of the robot via a new approach and is basically different from the conventional Jacobian matrix. One of the important characteristics of the new method is computational efficiency in comparison with the conventional method. It is demonstrated that using the new method instead of the conventional one, significantly reduces the computation time required to determine workspace of the robot as well as the time required to solve the equations of motion.

Torsional Vibration of a Shafting System under Electrical Disturbances

Directory of Open Access Journals (Sweden)

Ling Xiang

2012-01-01

Full Text Available Torsional vibration responses of a nonlinear shafting system are studied by a modified Riccati torsional transfer matrix combining with the Newmark-β method. Firstly, the system is modeled as a chain consisting of an elastic spring with concentrated mass points, from which a multi-segment lumped mass model is established. Secondly, accumulated errors are eliminated from the eigenfrequencies and responses of the system's torsional vibration by this newly developed procedure. The incremental transfer matrix method, combining the modified Riccati torsional transfer matrix with Newmark-β method, is further applied to solve the dynamical equations for the torsional vibration of the nonlinear shafting system. Lastly, the shafting system of a turbine-generator is employed as an illustrating example, and simulation analysis has been performed on the transient responses of the shaft's torsional vibrations during typical power network disturbances, such as three-phase short circuit, two-phase short circuit and asynchronous juxtaposition. The results validate the present method and are instructive for the design of a turbo-generator shaft.

Multi-matrix loop equations: algebraic and differential structures and an approximation based on deformation quantization

International Nuclear Information System (INIS)

Krishnaswami, Govind S.

2006-01-01

Large-N multi-matrix loop equations are formulated as quadratic difference equations in concatenation of gluon correlations. Though non-linear, they involve highest rank correlations linearly. They are underdetermined in many cases. Additional linear equations for gluon correlations, associated to symmetries of action and measure are found. Loop equations aren't differential equations as they involve left annihilation, which doesn't satisfy the Leibnitz rule with concatenation. But left annihilation is a derivation of the commutative shuffle product. Moreover shuffle and concatenation combine to define a bialgebra. Motivated by deformation quantization, we expand concatenation around shuffle in powers of q, whose physical value is 1. At zeroth order the loop equations become quadratic PDEs in the shuffle algebra. If the variation of the action is linear in iterated commutators of left annihilations, these quadratic PDEs linearize by passage to shuffle reciprocal of correlations. Remarkably, this is true for regularized versions of the Yang-Mills, Chern-Simons and Gaussian actions. But the linear equations are underdetermined just as the loop equations were. For any particular solution, the shuffle reciprocal is explicitly inverted to get the zeroth order gluon correlations. To go beyond zeroth order, we find a Poisson bracket on the shuffle algebra and associative q-products interpolating between shuffle and concatenation. This method, and a complementary one of deforming annihilation rather than product are shown to give over and underestimates for correlations of a gaussian matrix model

Symbolic computation and solitons of the nonlinear Schroedinger equation in inhomogeneous optical fiber media

International Nuclear Information System (INIS)

Li Biao; Chen Yong

2007-01-01

In this paper, the inhomogeneous nonlinear Schroedinger equation with the loss/gain and the frequency chirping is investigated. With the help of symbolic computation, three families of exact analytical solutions are presented by employing the extended projective Riccati equation method. From our results, many previous known results of nonlinear Schroedinger equation obtained by some authors can be recovered by means of some suitable selections of the arbitrary functions and arbitrary constants. Of optical and physical interests, soliton propagation and soliton interaction are discussed and simulated by computer, which include snake-soliton propagation and snake-solitons interaction, boomerang-like soliton propagation and boomerang-like solitons interaction, dispersion managed (DM) bright (dark) soliton propagation and DM solitons interaction

Optimal control of transverse mode coupling instability based on the two particle model

International Nuclear Information System (INIS)

Ogata, Atsushi

1985-01-01

The optimal regulator design technique is applied to asymptotically stabilize the transverse mode coupling instability of a storage ring. The state equations are based on the two particle model. These are a pair of equation sets, one for the first and one for the second half of the synchrotron phase. Each set consists of first-order difference equations in vector-matrix form, with time step equal to the revolution time of the ring. Solution of the discrete Riccati equation gives the optimal gain matrix of the transverse feedback. Computer simulations are carried out to verify its effectiveness. Some modifications necessary to apply it to the real accelerator operation are made. The old methods, the classical output feedback and the reactive feedback, are interpreted from the viewpoint of the optimal control. (orig.)

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

Science.gov (United States)

Hanks, Brantley R.; Skelton, Robert E.

1991-01-01

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

Hierarchical matrices algorithms and analysis

CERN Document Server

Hackbusch, Wolfgang

2015-01-01

This self-contained monograph presents matrix algorithms and their analysis. The new technique enables not only the solution of linear systems but also the approximation of matrix functions, e.g., the matrix exponential. Other applications include the solution of matrix equations, e.g., the Lyapunov or Riccati equation. The required mathematical background can be found in the appendix. The numerical treatment of fully populated large-scale matrices is usually rather costly. However, the technique of hierarchical matrices makes it possible to store matrices and to perform matrix operations approximately with almost linear cost and a controllable degree of approximation error. For important classes of matrices, the computational cost increases only logarithmically with the approximation error. The operations provided include the matrix inversion and LU decomposition. Since large-scale linear algebra problems are standard in scientific computing, the subject of hierarchical matrices is of interest to scientists ...

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

Science.gov (United States)

El-Gebeily, M.; Yushau, B.

2008-01-01

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

Closed form solutions of two time fractional nonlinear wave equations

Directory of Open Access Journals (Sweden)

M. Ali Akbar

2018-06-01

Full Text Available In this article, we investigate the exact traveling wave solutions of two nonlinear time fractional wave equations. The fractional derivatives are described in the sense of conformable fractional derivatives. In addition, the traveling wave solutions are accomplished in the form of hyperbolic, trigonometric, and rational functions involving free parameters. To investigate such types of solutions, we implement the new generalized (G′/G-expansion method. The extracted solutions are reliable, useful and suitable to comprehend the optimal control problems, chaotic vibrations, global and local bifurcations and resonances, furthermore, fission and fusion phenomena occur in solitons, the relativistic energy-momentum relation, scalar electrodynamics, quantum relativistic one-particle theory, electromagnetic interactions etc. The results reveal that the method is very fruitful and convenient for exploring nonlinear differential equations of fractional order treated in theoretical physics. Keywords: Traveling wave solution, Soliton, Generalized (G′/G-expansion method, Time fractional Duffing equation, Time fractional Riccati equation

Some Matrix Iterations for Computing Generalized Inverses and Balancing Chemical Equations

OpenAIRE

Soleimani, Farahnaz; Stanimirovi´c, Predrag; Soleymani, Fazlollah

2015-01-01

An application of iterative methods for computing the Moore–Penrose inverse in balancing chemical equations is considered. With the aim to illustrate proposed algorithms, an improved high order hyper-power matrix iterative method for computing generalized inverses is introduced and applied. The improvements of the hyper-power iterative scheme are based on its proper factorization, as well as on the possibility to accelerate the iterations in the initial phase of the convergence. Although the ...

Time discretization of the point kinetic equations using matrix exponential method and First-Order Hold

International Nuclear Information System (INIS)

Park, Yujin; Kazantzis, Nikolaos; Parlos, Alexander G.; Chong, Kil To

2013-01-01

Highlights: • Numerical solution for stiff differential equations using matrix exponential method. • The approximation is based on First Order Hold assumption. • Various input examples applied to the point kinetics equations. • The method shows superior useful and effective activity. - Abstract: A system of nonlinear differential equations is derived to model the dynamics of neutron density and the delayed neutron precursors within a point kinetics equation modeling framework for a nuclear reactor. The point kinetic equations are mathematically characterized as stiff, occasionally nonlinear, ordinary differential equations, posing significant challenges when numerical solutions are sought and traditionally resulting in the need for smaller time step intervals within various computational schemes. In light of the above realization, the present paper proposes a new discretization method inspired by system-theoretic notions and technically based on a combination of the matrix exponential method (MEM) and the First-Order Hold (FOH) assumption. Under the proposed time discretization structure, the sampled-data representation of the nonlinear point kinetic system of equations is derived. The performance of the proposed time discretization procedure is evaluated using several case studies with sinusoidal reactivity profiles and multiple input examples (reactivity and neutron source function). It is shown, that by applying the proposed method under a First-Order Hold for the neutron density and the precursor concentrations at each time step interval, the stiffness problem associated with the point kinetic equations can be adequately addressed and resolved. Finally, as evidenced by the aforementioned detailed simulation studies, the proposed method retains its validity and accuracy for a wide range of reactor operating conditions, including large sampling periods dictated by physical and/or technical limitations associated with the current state of sensor and

The response matrix discrete ordinates solution to the 1D radiative transfer equation

International Nuclear Information System (INIS)

Ganapol, Barry D.

2015-01-01

The discrete ordinates method (DOM) of solution to the 1D radiative transfer equation has been an effective method of solution for nearly 70 years. During that time, the method has experienced numerous improvements as numerical and computational techniques have become more powerful and efficient. Here, we again consider the analytical solution to the discrete radiative transfer equation in a homogeneous medium by proposing a new, and consistent, form of solution that improves upon previous forms. Aided by a Wynn-epsilon convergence acceleration, its numerical evaluation can achieve extreme precision as demonstrated by comparison with published benchmarks. Finally, we readily extend the solution to a heterogeneous medium through the star product formulation producing a novel benchmark for closed form Henyey–Greenstein scattering as an example. - Highlights: • Presents a new solution to the RTE called the response matrix DOM (RM/DOM). • Solution representations avoid the instability common in exponential solutions. • Explicit form in terms of matrix hyperbolic functions. • Extreme accuracy through Wynn-epsilon acceleration checked by published benchmarks. • Provides a more transparent numerical evaluation than found previously

Green's matrix for a second-order self-adjoint matrix differential operator

International Nuclear Information System (INIS)

Sisman, Tahsin Cagri; Tekin, Bayram

2010-01-01

A systematic construction of the Green's matrix for a second-order self-adjoint matrix differential operator from the linearly independent solutions of the corresponding homogeneous differential equation set is carried out. We follow the general approach of extracting the Green's matrix from the Green's matrix of the corresponding first-order system. This construction is required in the cases where the differential equation set cannot be turned to an algebraic equation set via transform techniques.

Hamiltonians and variational principles for Alfvén simple waves

International Nuclear Information System (INIS)

Webb, G M; Hu, Q; Roux, J A le; Dasgupta, B; Zank, G P

2012-01-01

The evolution equations for the magnetic field induction B with the wave phase for Alfvén simple waves are expressed as variational principles and in the Hamiltonian form. The evolution of B with the phase (which is a function of the space and time variables) depends on the generalized Frenet–Serret equations, in which the wave normal n (which is a function of the phase) is taken to be tangent to a curve X, in a 3D Cartesian geometry vector space. The physical variables (the gas density, fluid velocity, gas pressure and magnetic field induction) in the wave depend only on the phase. Three approaches are developed. One approach exploits the fact that the Frenet equations may be written as a 3D Hamiltonian system, which can be described using the Nambu bracket. It is shown that B as a function of the phase satisfies a modified version of the Frenet equations, and hence the magnetic field evolution equations can be expressed in the Hamiltonian form. A second approach develops an Euler–Poincaré variational formulation. A third approach uses the Frenet frame formulation, in which the hodograph of B moves on a sphere of constant radius and uses a stereographic projection transformation due to Darboux. The equations for the projected field components reduce to a complex Riccati equation. By using a Cole–Hopf transformation, the Riccati equation reduces to a linear second order differential equation for the new variable. A Hamiltonian formulation of the second order differential equation then allows the system to be written in the Hamiltonian form. Alignment dynamics equations for Alfvén simple waves give rise to a complex Riccati equation or, equivalently, to a quaternionic Riccati equation, which can be mapped onto the Riccati equation obtained by stereographic projection. (paper)

Nonlocal Symmetries, Consistent Riccati Expansion, and Analytical Solutions of the Variant Boussinesq System

Science.gov (United States)

Feng, Lian-Li; Tian, Shou-Fu; Zhang, Tian-Tian; Zhou, Jun

2017-07-01

Under investigation in this paper is the variant Boussinesq system, which describes the propagation of surface long wave towards two directions in a certain deep trough. With the help of the truncated Painlevé expansion, we construct its nonlocal symmetry, Bäcklund transformation, and Schwarzian form, respectively. The nonlocal symmetries can be localised to provide the corresponding nonlocal group, and finite symmetry transformations and similarity reductions are computed. Furthermore, we verify that the variant Boussinesq system is solvable via the consistent Riccati expansion (CRE). By considering the consistent tan-function expansion (CTE), which is a special form of CRE, the interaction solutions between soliton and cnoidal periodic wave are explicitly studied.

Homotopy perturbation transform method for pricing under pure diffusion models with affine coefficients

Directory of Open Access Journals (Sweden)

Claude Rodrigue Bambe Moutsinga

2018-01-01

Full Text Available Most existing multivariate models in finance are based on diffusion models. These models typically lead to the need of solving systems of Riccati differential equations. In this paper, we introduce an efficient method for solving systems of stiff Riccati differential equations. In this technique, a combination of Laplace transform and homotopy perturbation methods is considered as an algorithm to the exact solution of the nonlinear Riccati equations. The resulting technique is applied to solving stiff diffusion model problems that include interest rates models as well as two and three-factor stochastic volatility models. We show that the present approach is relatively easy, efficient and highly accurate.

The Green's matrix and the boundary integral equations for analysis of time-harmonic dynamics of elastic helical springs.

Science.gov (United States)

Sorokin, Sergey V

2011-03-01

Helical springs serve as vibration isolators in virtually any suspension system. Various exact and approximate methods may be employed to determine the eigenfrequencies of vibrations of these structural elements and their dynamic transfer functions. The method of boundary integral equations is a meaningful alternative to obtain exact solutions of problems of the time-harmonic dynamics of elastic springs in the framework of Bernoulli-Euler beam theory. In this paper, the derivations of the Green's matrix, of the Somigliana's identities, and of the boundary integral equations are presented. The vibrational power transmission in an infinitely long spring is analyzed by means of the Green's matrix. The eigenfrequencies and the dynamic transfer functions are found by solving the boundary integral equations. In the course of analysis, the essential features and advantages of the method of boundary integral equations are highlighted. The reported analytical results may be used to study the time-harmonic motion in any wave guide governed by a system of linear differential equations in a single spatial coordinate along its axis. © 2011 Acoustical Society of America

Modifying a numerical algorithm for solving the matrix equation X + AX T B = C

Science.gov (United States)

Vorontsov, Yu. O.

2013-06-01

Certain modifications are proposed for a numerical algorithm solving the matrix equation X + AX T B = C. By keeping the intermediate results in storage and repeatedly using them, it is possible to reduce the total complexity of the algorithm from O( n 4) to O( n 3) arithmetic operations.

Oscillation results for certain fractional difference equations

Directory of Open Access Journals (Sweden)

Zhiyun WANG

2017-08-01

Full Text Available Fractional calculus is a theory that studies the properties and application of arbitrary order differentiation and integration. It can describe the physical properties of some systems more accurately, and better adapt to changes in the system, playing an important role in many fields. For example, it can describe the process of tumor growth (growth stimulation and growth inhibition in biomedical science. The oscillation of solutions of two kinds of fractional difference equations is studied, mainly using the proof by contradiction, that is, assuming the equation has a nonstationary solution. For the first kind of equation, the function symbol is firstly determined, and by constructing the Riccati function, the difference is calculated. Then the condition of the function is used to satisfy the contradiction, that is, the assumption is false, which verifies the oscillation of the solution. For the second kind of equation with initial condition, the equivalent fractional sum form of the fractional difference equation are firstly proved. With considering 0<α≤1 and α>1, respectively, by using the properties of Stirling formula and factorial function, the contradictory is got through enhanced processing, namely the assuming is not established, and the sufficient condition for the bounded solutions of the fractional difference equation is obtained. The above results will optimize the relevant conclusions and enrich the relevant results. The results are applied to the specific equations, and the oscillation of the solutions of equations is proved.

Optimal pole shifting controller for interconnected power system

International Nuclear Information System (INIS)

Yousef, Ali M.; Kassem, Ahmed M.

2011-01-01

Research highlights: → Mathematical model represents a power system which consists of synchronous machine connected to infinite bus through transmission line. → Power system stabilizer was designed based on optimal pole shifting controller. → The system performances was tested through load disturbances at different operating conditions. → The system performance with the proposed optimal pole shifting controller is compared with the conventional pole placement controller. → The digital simulation results indicated that the proposed controller has a superior performance. -- Abstract: Power system stabilizer based on optimal pole shifting is proposed. An approach for shifting the real parts of the open-loop poles to any desired positions while preserving the imaginary parts is presented. In each step of this approach, it is required to solve a first-order or a second-order linear matrix Lyapunov equation for shifting one real pole or two complex conjugate poles, respectively. This presented method yields a solution, which is optimal with respect to a quadratic performance index. The attractive feature of this method is that it enables solutions of the complex problem to be easily found without solving any non-linear algebraic Riccati equation. The present power system stabilizer is based on Riccati equation approach. The control law depends on finding the feedback gain matrix, and then the control signal is synthesized by multiplying the state variables of the power system with determined gain matrix. The gain matrix is calculated one time only, and it works over wide range of operating conditions. To validate the power of the proposed PSS, a linearized model of a simple power system consisted of a single synchronous machine connected to infinite bus bar through transmission line is simulated. The studied power system is subjected to various operating points and power system parameters changes.

Optimal pole shifting controller for interconnected power system

Energy Technology Data Exchange (ETDEWEB)

Yousef, Ali M., E-mail: drali_yousef@yahoo.co [Electrical Eng. Dept., Faculty of Engineering, Assiut University (Egypt); Kassem, Ahmed M., E-mail: kassem_ahmed53@hotmail.co [Control Technology Dep., Industrial Education College, Beni-Suef University (Egypt)

2011-05-15

Research highlights: {yields} Mathematical model represents a power system which consists of synchronous machine connected to infinite bus through transmission line. {yields} Power system stabilizer was designed based on optimal pole shifting controller. {yields} The system performances was tested through load disturbances at different operating conditions. {yields} The system performance with the proposed optimal pole shifting controller is compared with the conventional pole placement controller. {yields} The digital simulation results indicated that the proposed controller has a superior performance. -- Abstract: Power system stabilizer based on optimal pole shifting is proposed. An approach for shifting the real parts of the open-loop poles to any desired positions while preserving the imaginary parts is presented. In each step of this approach, it is required to solve a first-order or a second-order linear matrix Lyapunov equation for shifting one real pole or two complex conjugate poles, respectively. This presented method yields a solution, which is optimal with respect to a quadratic performance index. The attractive feature of this method is that it enables solutions of the complex problem to be easily found without solving any non-linear algebraic Riccati equation. The present power system stabilizer is based on Riccati equation approach. The control law depends on finding the feedback gain matrix, and then the control signal is synthesized by multiplying the state variables of the power system with determined gain matrix. The gain matrix is calculated one time only, and it works over wide range of operating conditions. To validate the power of the proposed PSS, a linearized model of a simple power system consisted of a single synchronous machine connected to infinite bus bar through transmission line is simulated. The studied power system is subjected to various operating points and power system parameters changes.

Numerical solution of multiband k.p model for tunnelling in type-II heterostructures

Directory of Open Access Journals (Sweden)

A.E. Botha

2010-01-01

Full Text Available A new and very general method was developed for calculating the charge and spin-resolved electron tunnelling in type-II heterojunctions. Starting from a multiband k.p description of the bulk energy-band structure, a multiband k.p Riccati equation was derived. The reflection and transmission coefficients were obtained for each channel by integrating the Riccati equation over the entire heterostructure. Numerical instability was reduced through this method, in which the multichannel log-derivative of the envelope function matrix, rather than the envelope function itself, was propagated. As an example, a six-band k.p Hamiltonian was used to calculate the current-voltage characteristics of a 10-nm wide InAs/ GaSb/InAs single quantum well device which exhibited negative differential resistance at room temperature. The calculated current as a function of applied (bias voltage was found to be in semiquantitative agreement with the experiment, a result which indicated that inelastic transport mechanisms do not contribute significantly to the valley currents measured in this particular device.

Second level semi-degenerate fields in W{sub 3} Toda theory: matrix element and differential equation

Energy Technology Data Exchange (ETDEWEB)

Belavin, Vladimir [I.E. Tamm Department of Theoretical Physics, P.N. Lebedev Physical Institute,Leninsky Avenue 53, 119991 Moscow (Russian Federation); Department of Quantum Physics, Institute for Information Transmission Problems,Bolshoy Karetny per. 19, 127994 Moscow (Russian Federation); Moscow Institute of Physics and Technology,Dolgoprudnyi, 141700 Moscow region (Russian Federation); Cao, Xiangyu [LPTMS, CNRS (UMR 8626), Université Paris-Saclay,15 rue Georges Clémenceau, 91405 Orsay (France); Estienne, Benoit [LPTHE, CNRS and Université Pierre et Marie Curie, Sorbonne Universités,4 Place Jussieu, 75252 Paris Cedex 05 (France); Santachiara, Raoul [LPTMS, CNRS (UMR 8626), Université Paris-Saclay,15 rue Georges Clémenceau, 91405 Orsay (France)

2017-03-02

In a recent study we considered W{sub 3} Toda 4-point functions that involve matrix elements of a primary field with the highest-weight in the adjoint representation of sl{sub 3}. We generalize this result by considering a semi-degenerate primary field, which has one null vector at level two. We obtain a sixth-order Fuchsian differential equation for the conformal blocks. We discuss the presence of multiplicities, the matrix elements and the fusion rules.

Splitting of the rate matrix as a definition of time reversal in master equation systems

International Nuclear Information System (INIS)

Liu Fei; Le, Hong

2012-01-01

Motivated by recent progress in nonequilibrium fluctuation relations, we present a generalized time reversal for stochastic master equation systems with discrete states, which is defined as a splitting of the rate matrix into irreversible and reversible parts. An immediate advantage of this definition is that a variety of fluctuation relations can be attributed to different matrix splittings. Additionally, we find that the accustomed total entropy production formula and conditions of the detailed balance must be modified appropriately to account for the reversible rate part, which was previously ignored. (paper)

The {P,Q,k+1}-Reflexive Solution to System of Matrix Equations AX=C, XB=D

Directory of Open Access Journals (Sweden)

Chang-Zhou Dong

2015-01-01

Full Text Available Let P∈Cm×m and Q∈Cn×n be Hermitian and {k+1}-potent matrices; that is, Pk+1=P=P⁎ and Qk+1=Q=Q⁎, where ·⁎ stands for the conjugate transpose of a matrix. A matrix X∈Cm×n is called {P,Q,k+1}-reflexive (antireflexive if PXQ=X (PXQ=-X. In this paper, the system of matrix equations AX=C and XB=D subject to {P,Q,k+1}-reflexive and antireflexive constraints is studied by converting into two simpler cases: k=1 and k=2. We give the solvability conditions and the general solution to this system; in addition, the least squares solution is derived; finally, the associated optimal approximation problem for a given matrix is considered.

Modified method of simplest equation: Powerful tool for obtaining exact and approximate traveling-wave solutions of nonlinear PDEs

Science.gov (United States)

Vitanov, Nikolay K.

2011-03-01

We discuss the class of equations ∑i,j=0mAij(u){∂iu}/{∂ti}∂+∑k,l=0nBkl(u){∂ku}/{∂xk}∂=C(u) where Aij( u), Bkl( u) and C( u) are functions of u( x, t) as follows: (i) Aij, Bkl and C are polynomials of u; or (ii) Aij, Bkl and C can be reduced to polynomials of u by means of Taylor series for small values of u. For these two cases the above-mentioned class of equations consists of nonlinear PDEs with polynomial nonlinearities. We show that the modified method of simplest equation is powerful tool for obtaining exact traveling-wave solution of this class of equations. The balance equations for the sub-class of traveling-wave solutions of the investigated class of equations are obtained. We illustrate the method by obtaining exact traveling-wave solutions (i) of the Swift-Hohenberg equation and (ii) of the generalized Rayleigh equation for the cases when the extended tanh-equation or the equations of Bernoulli and Riccati are used as simplest equations.

The Split Coefficient Matrix method for hyperbolic systems of gasdynamic equations

Science.gov (United States)

Chakravarthy, S. R.; Anderson, D. A.; Salas, M. D.

1980-01-01

The Split Coefficient Matrix (SCM) finite difference method for solving hyperbolic systems of equations is presented. This new method is based on the mathematical theory of characteristics. The development of the method from characteristic theory is presented. Boundary point calculation procedures consistent with the SCM method used at interior points are explained. The split coefficient matrices that define the method for steady supersonic and unsteady inviscid flows are given for several examples. The SCM method is used to compute several flow fields to demonstrate its accuracy and versatility. The similarities and differences between the SCM method and the lambda-scheme are discussed.

A self-consistent nodal method in response matrix formalism for the multigroup diffusion equations

International Nuclear Information System (INIS)

Malambu, E.M.; Mund, E.H.

1996-01-01

We develop a nodal method for the multigroup diffusion equations, based on the transverse integration procedure (TIP). The efficiency of the method rests upon the convergence properties of a high-order multidimensional nodal expansion and upon numerical implementation aspects. The discrete 1D equations are cast in response matrix formalism. The derivation of the transverse leakage moments is self-consistent i.e. does not require additional assumptions. An outstanding feature of the method lies in the linear spatial shape of the local transverse leakage for the first-order scheme. The method is described in the two-dimensional case. The method is validated on some classical benchmark problems. (author)

Parallel Computation of the Jacobian Matrix for Nonlinear Equation Solvers Using MATLAB

Science.gov (United States)

Rose, Geoffrey K.; Nguyen, Duc T.; Newman, Brett A.

2017-01-01

Demonstrating speedup for parallel code on a multicore shared memory PC can be challenging in MATLAB due to underlying parallel operations that are often opaque to the user. This can limit potential for improvement of serial code even for the so-called embarrassingly parallel applications. One such application is the computation of the Jacobian matrix inherent to most nonlinear equation solvers. Computation of this matrix represents the primary bottleneck in nonlinear solver speed such that commercial finite element (FE) and multi-body-dynamic (MBD) codes attempt to minimize computations. A timing study using MATLAB's Parallel Computing Toolbox was performed for numerical computation of the Jacobian. Several approaches for implementing parallel code were investigated while only the single program multiple data (spmd) method using composite objects provided positive results. Parallel code speedup is demonstrated but the goal of linear speedup through the addition of processors was not achieved due to PC architecture.

Optimal linear-quadratic control of coupled parabolic-hyperbolic PDEs

Science.gov (United States)

Aksikas, I.; Moghadam, A. Alizadeh; Forbes, J. F.

2017-10-01

This paper focuses on the optimal control design for a system of coupled parabolic-hypebolic partial differential equations by using the infinite-dimensional state-space description and the corresponding operator Riccati equation. Some dynamical properties of the coupled system of interest are analysed to guarantee the existence and uniqueness of the solution of the linear-quadratic (LQ)-optimal control problem. A state LQ-feedback operator is computed by solving the operator Riccati equation, which is converted into a set of algebraic and differential Riccati equations, thanks to the eigenvalues and the eigenvectors of the parabolic operator. The results are applied to a non-isothermal packed-bed catalytic reactor. The LQ-optimal controller designed in the early portion of the paper is implemented for the original nonlinear model. Numerical simulations are performed to show the controller performances.

A Discrete-Time Recurrent Neural Network for Solving Rank-Deficient Matrix Equations With an Application to Output Regulation of Linear Systems.

Science.gov (United States)

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.

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.

Graphical method for determination of critical dimensions and neutron flux distribution in multi zone nuclear reactors; Graficka metoda za odredjivanje kriticnih dimenzija i raspodele fluksa kod multiregionalnih nuklearnih reaktora

Energy Technology Data Exchange (ETDEWEB)

Radanovic, Lj; Bingulac, S; Lazarevic, B; Matausek, M [The Institute of Nuclear Sciences Boris Kidric, Vinca, Beograd (Yugoslavia)

1964-07-01

This paper describes the graphical method for calculating the neutron flux distribution by using normalized Riccati equations. It was shown that the solutions of adequately normalized Riccati equations could be used as standard curves for determining the critical dimensions and radial flux distribution in multi zone nuclear reactors. the methos is applicable irrelevant of the number and position of the region in the core.

The Basics of Anisotropy-Based Analysis of Discrete Time-Invariant Systems

Directory of Open Access Journals (Sweden)

I. R. Belov

2017-01-01

Full Text Available When investigating a behavior of dynamical systems, we should take into account the external noises, which have an effect on the system. The article introduces a concept of the anisotropy-based norm of the system as one of the ways to describe the measure of the effect of external disturbances on the system. The definition of the anisotropic norm includes some concepts from information theory, such as relative entropy and anisotropy. The theoretical section at the beginning of the article describes these definitions. The considered norm of the system can be evaluated in several ways. The article examines two methods - in the frequency domain and in the state space. To find the norm in the state space it is necessary to find the solution of the Riccati equation. This problem is rather laborious. So the algorithms to avoid the solution of Riccati equation are used in application of anisotropy-based norm’s evaluation methods. The principle of these algorithms is replacement of Riccati equation by the system of linear matrix inequalities. The software implementation of methods under consideration is designed using the MATLAB packages. The calculation results of the anisotropy-based norm for a given linear discrete system are obtained using this program. The article presents these results as graphs.This article enters into the Master's qualifying work "Basic quality criteria in the theory of linear systems". In this paper we consider various quality criteria, the solution of the optimal control problem for each of them, and compare the results obtained for different criteria. The anisotropy-based norm considered in the article is one of the quality criteria.

Solution of quadratic matrix equations for free vibration analysis of structures.

Science.gov (United States)

Gupta, K. K.

1973-01-01

An efficient digital computer procedure and the related numerical algorithm are presented herein for the solution of quadratic matrix equations associated with free vibration analysis of structures. Such a procedure enables accurate and economical analysis of natural frequencies and associated modes of discretized structures. The numerically stable algorithm is based on the Sturm sequence method, which fully exploits the banded form of associated stiffness and mass matrices. The related computer program written in FORTRAN V for the JPL UNIVAC 1108 computer proves to be substantially more accurate and economical than other existing procedures of such analysis. Numerical examples are presented for two structures - a cantilever beam and a semicircular arch.

General concepts of multichannel collision theory and their translation into the matrix formulation of few-body integral equations

International Nuclear Information System (INIS)

Sandhas, W.

1978-01-01

In the N-body problem mappings between channel states and scattering states are studied. It is shown in particular that the (2sup(N-1)-1) two-fragment MOELLER operators introduced on the whole Hilbert space are sufficient to provide all multi-fragment scattering states. Hence, each of these states is uniquely determined by (2sup(N-1)-1) Lippmann-Schwinger (LS) equations. Rewriting every set of LS equations as one matrix equation, current few-body approaches are incorporated in a rather natural way. The typical uniqueness questions of such coupled systems are discussed, and it is shown that Faddeev-type couplings lead to unique equations for arbitrary N. (author)

General concepts of multichannel collision theory and their translation into the matrix formulation of few-body integral equations

International Nuclear Information System (INIS)

Sandhas, W.

1978-04-01

In the N-body problem mappings between channel states and scattering states are studied. It is shown in particular that the (2sup(N-1)-1) two-fragment Moeller operators introduced on the whole Hilbert space are sufficient to provide all multifragment scattering states. Hence, each of these states is uniquly determined by (2sup(N-1)-1) Lippmann-Schwinger (LS) equations. Rewriting every set of LS equations as one matrix equation, current few-body approaches are incorporated in a rather natural way. The typical uniqueness questions of such coupled systems are discussed, and it si shown that Faddeev-type couplings lead to unique equations for arbitrary N. (orig.) [de

Nonlocal Symmetries, Conservation Laws and Interaction Solutions of the Generalised Dispersive Modified Benjamin-Bona-Mahony Equation

Science.gov (United States)

Yan, Xue-Wei; Tian, Shou-Fu; Dong, Min-Jie; Wang, Xiu-Bin; Zhang, Tian-Tian

2018-05-01

We consider the generalised dispersive modified Benjamin-Bona-Mahony equation, which describes an approximation status for long surface wave existed in the non-linear dispersive media. By employing the truncated Painlevé expansion method, we derive its non-local symmetry and Bäcklund transformation. The non-local symmetry is localised by a new variable, which provides the corresponding non-local symmetry group and similarity reductions. Moreover, a direct method can be provided to construct a kind of finite symmetry transformation via the classic Lie point symmetry of the normal prolonged system. Finally, we find that the equation is a consistent Riccati expansion solvable system. With the help of the Jacobi elliptic function, we get its interaction solutions between solitary waves and cnoidal periodic waves.

H{infinity} Filtering for Dynamic Compensation of Self-Powered Neutron Detectors - A Linear Matrix Inequality Based Method -

Energy Technology Data Exchange (ETDEWEB)

Park, M.G.; Kim, Y.H.; Cha, K.H.; Kim, M.K. [Korea Electric Power Research Institute, Taejon (Korea)

1999-07-01

A method is described to develop and H{infinity} filtering method for the dynamic compensation of self-powered neutron detectors normally used for fixed incore instruments. An H{infinity} norm of the filter transfer matrix is used as the optimization criteria in the worst-case estimation error sense. Filter modeling is performed for both continuous- and discrete-time models. The filter gains are optimized in the sense of noise attenuation level of H{infinity} setting. By introducing Bounded Real Lemma, the conventional algebraic Riccati inequalities are converted into Linear Matrix Inequalities (LMIs). Finally, the filter design problem is solved via the convex optimization framework using LMIs. The simulation results show that remarkable improvements are achieved in view of the filter response time and the filter design efficiency. (author). 15 refs., 4 figs., 3 tabs.

Leader-Following Consensus Stability of Discrete-Time Linear Multiagent Systems with Observer-Based Protocols

Directory of Open Access Journals (Sweden)

Bingbing Xu

2013-01-01

Full Text Available We consider the leader-following consensus problem of discrete-time multiagent systems on a directed communication topology. Two types of distributed observer-based consensus protocols are considered to solve such a problem. The observers involved in the proposed protocols include full-order observer and reduced-order observer, which are used to reconstruct the state variables. Two algorithms are provided to construct the consensus protocols, which are based on the modified discrete-time algebraic Riccati equation and Sylvester equation. In light of graph and matrix theory, some consensus conditions are established. Finally, a numerical example is provided to illustrate the obtained result.

Spatial charge motion on an uniform density matrix-general equations in opened and closed circuits

International Nuclear Information System (INIS)

Aguiar Monsanto, S. de.

1983-01-01

The motion of a space charge cloud embedded in a matrix of constant immobile charge density is studied in open as well as in closed circuit. In the first case, open circuit, the solution is almost trivial as compared as the other one in which, after some work, the problem is reduced to an ordinary differential equation. The method of solution is parallel to that employed in the study of monopolar free space charge motion. The voltage and the current produced by a system with no net charge but with unbalanced local charge density were calculated using the general equations derived in the first part of the work. (Author) [pt

Octonionic matrix representation and electromagnetism

Energy Technology Data Exchange (ETDEWEB)

Chanyal, B. C. [Kumaun University, S. S. J. Campus, Almora (India)

2014-12-15

Keeping in mind the important role of octonion algebra, we have obtained the electromagnetic field equations of dyons with an octonionic 8 x 8 matrix representation. In this paper, we consider the eight - dimensional octonionic space as a combination of two (external and internal) four-dimensional spaces for the existence of magnetic monopoles (dyons) in a higher-dimensional formalism. As such, we describe the octonion wave equations in terms of eight components from the 8 x 8 matrix representation. The octonion forms of the generalized potential, fields and current source of dyons in terms of 8 x 8 matrix are discussed in a consistent manner. Thus, we have obtained the generalized Dirac-Maxwell equations of dyons from an 8x8 matrix representation of the octonion wave equations in a compact and consistent manner. The generalized Dirac-Maxwell equations are fully symmetric Maxwell equations and allow for the possibility of magnetic charges and currents, analogous to electric charges and currents. Accordingly, we have obtained the octonionic Dirac wave equations in an external field from the matrix representation of the octonion-valued potentials of dyons.

Unbiased minimum variance estimator of a matrix exponential function. Application to Boltzmann/Bateman coupled equations solving

International Nuclear Information System (INIS)

Dumonteil, E.; Diop, C. M.

2009-01-01

This paper derives an unbiased minimum variance estimator (UMVE) of a matrix exponential function of a normal wean. The result is then used to propose a reference scheme to solve Boltzmann/Bateman coupled equations, thanks to Monte Carlo transport codes. The last section will present numerical results on a simple example. (authors)

Pseudospectral operational matrix for numerical solution of single and multiterm time fractional diffusion equation

OpenAIRE

GHOLAMI, SAEID; BABOLIAN, ESMAIL; JAVIDI, MOHAMMAD

2016-01-01

This paper presents a new numerical approach to solve single and multiterm time fractional diffusion equations. In this work, the space dimension is discretized to the Gauss$-$Lobatto points. We use the normalized Grunwald approximation for the time dimension and a pseudospectral successive integration matrix for the space dimension. This approach shows that with fewer numbers of points, we can approximate the solution with more accuracy. Some examples with numerical results in tables and fig...

The string difference equation of the D = 1 matrix model and W1+∞ symmetry of the KP hierarchy

International Nuclear Information System (INIS)

Awada, M.A.; Sin, S.J.

1992-01-01

In this paper, the authors give a connection between the D = 1 matrix model and the generalized KP hierarchy. First, the authors find a difference equation satisfied by F, the Legendre transformation of the free energy of the D = 1 matrix model on a circle of radius R. Then the authors show that it is a special case of the difference equation of the generalized KP hierarchy with its zero mode identified with the scaling variable of the D = 1 string theory. The authors propose that the massive D = 1 matrix model is described by the generalized KP hierarchy, which implies the manifest integrability of D = 1 string theory. The authors also show that the (generalized) KP hierarchy has an underlying W 1 + ∞ symmetry. By reduction, we prove that the generalized KdV hierarchy has a subalgebra of the above symmetry which again forms a W 1+ ∞ . The authors argue that there are no W constraints in D = 1 string theory, which is in contrast to D 1 + ∞ constraints

Entropy Measures as Geometrical Tools in the Study of Cosmology

Directory of Open Access Journals (Sweden)

Gilbert Weinstein

2017-12-01

Full Text Available Classical chaos is often characterized as exponential divergence of nearby trajectories. In many interesting cases these trajectories can be identified with geodesic curves. We define here the entropy by S = ln χ ( x with χ ( x being the distance between two nearby geodesics. We derive an equation for the entropy, which by transformation to a Riccati-type equation becomes similar to the Jacobi equation. We further show that the geodesic equation for a null geodesic in a double-warped spacetime leads to the same entropy equation. By applying a Robertson–Walker metric for a flat three-dimensional Euclidean space expanding as a function of time, we again reach the entropy equation stressing the connection between the chosen entropy measure and time. We finally turn to the Raychaudhuri equation for expansion, which also is a Riccati equation similar to the transformed entropy equation. Those Riccati-type equations have solutions of the same form as the Jacobi equation. The Raychaudhuri equation can be transformed to a harmonic oscillator equation, and it has been shown that the geodesic deviation equation of Jacobi is essentially equivalent to that of a harmonic oscillator. The Raychaudhuri equations are strong geometrical tools in the study of general relativity and cosmology. We suggest a refined entropy measure applicable in cosmology and defined by the average deviation of the geodesics in a congruence.

Interval Oscillation Criteria for Super-Half-Linear Impulsive Differential Equations with Delay

Directory of Open Access Journals (Sweden)

Zhonghai Guo

2012-01-01

Full Text Available We study the following second-order super-half-linear impulsive differential equations with delay [r(tφγ(x′(t]′+p(tφγ(x(t-σ+q(tf(x(t-σ=e(t, t≠τk, x(t+=akx(t, x′(t+=bkx′(t, t=τk, where t≥t0∈ℝ, φ*(u=|u|*-1u, σ is a nonnegative constant, {τk} denotes the impulsive moments sequence with τ1σ. By some classical inequalities, Riccati transformation, and two classes of functions, we give several interval oscillation criteria which generalize and improve some known results. Moreover, we also give two examples to illustrate the effectiveness and nonemptiness of our results.

A matrix formalism to solve interface condition equations in a reactor system

Energy Technology Data Exchange (ETDEWEB)

Matausek, M V [Boris Kidric Institute of Nuclear Sciences Vinca, Beograd (Yugoslavia)

1970-05-15

When a nuclear reactor or a reactor lattice cell is treated by an approximate procedure to solve the neutron transport equation, as the last computational step often appears a problem of solving systems of algebraic equations stating the interface and boundary conditions for the neutron flux moments. These systems have usually the coefficient matrices of the block-bi diagonal type, containing thus a large number of zero elements. In the present report it is shown how such a system can be solved efficiently accounting for all the zero elements both in the coefficient matrix and in the free term vector. The procedure is presented here for the case of multigroup P{sub 3} calculation of neutron flux distribution in a cylindrical reactor lattice cell. Compared with the standard gaussian elimination method, this procedure is more advantageous both in respect to the number of operations needed to solve a given problem and in respect to the computer memory storage requirements. A similar formalism can also be applied for other approximate methods, for instance for multigroup diffusion treatment of a multi zone reactor. (author)

Upper and Lower Bounds of Frequency Interval Gramians for a Class of Perturbed Linear Systems

DEFF Research Database (Denmark)

Shaker, Hamid Reza

2012-01-01

if the system is controllable or observable, but also it is required to know the degree of controllability or observability of the system. Gramian matrices were introduced to address this issue by providing a quantitative measure for controllability and observability. In many applications, the information...... of uncertain systems. In this paper, we derive upper and lower bounds of frequency interval gramians under perturbations of an A-matrix in the state-space form. These bounds are obtained by solving algebraic Riccati equations. The results are further used to obtain upper and lower bounds of the frequency...

Gravitational radiation from nearly Newtonian systems

International Nuclear Information System (INIS)

Kirk, E.M.

1989-09-01

A method of examining gravitational radiation from nearly Newtonian systems is presented. Using the Cartan formulation of Newtonian gravity, a one parameter family of space-times which have a strict Newtonian limit is constructed. An expression for the initial null data in terms of the Newtonian potential is obtained in the Newtonian limit. Using this, the problem is formulated as a series in the Newtonian parameter. The series expansions for the sources of the Bianchi identities are obtained to third order in both the vacuum and non-vacuum cases. A simple technique is presented for determining whether a particular source term gives rise to asymptotically flat null data. The far field quadrupole formula is derived in a leading approximation and a method for obtaining error bounds is discussed. Additionally, a method for solving Einstein's equations is shown. This involves expressing the Ricci identities as a matrix, Riccati equation and a system of linear matrix equations. A comparison of the formalisms of Bondi and Newman Penrose is presented and explicit correspondences between the supersurface constrain equations and the Ricci identities are shown. (author)

Special Cases in Optimization Problems for Stationary Linear Closed-Loop Systems

Science.gov (United States)

Aliev, F. A.; Larin, Vladimir Borisovich

2003-03-01

Consideration is given to problems of solving the algebraic Riccati equation (ARE)— J-factorization of matrix polynomials and J-factorization of rational matrices—to which traditional solution algorithms are not applicable. In this connection, solution algorithms for these problems are discussed where the eigenvalues of the Hamiltonian matrix corresponding to the ARE and the zeros of matrix polynomials are located on the imaginary axis. Moreover, a procedure is set forth for asymptotic expansion of a stabilizing solution of the ARE in the neighborhood of a point at which the ARE has no stabilizing solution. It is shown how this expansion can be used for constructing canonical J-factorization of matrix polynomials that is nearly a noncanonical J-factorization. It is pointed out that the algorithms described can be implemented with the help of MATLAB routines

A Novel Operational Matrix of Caputo Fractional Derivatives of Fibonacci Polynomials: Spectral Solutions of Fractional Differential Equations

Directory of Open Access Journals (Sweden)

Waleed M. Abd-Elhameed

2016-09-01

Full Text Available Herein, two numerical algorithms for solving some linear and nonlinear fractional-order differential equations are presented and analyzed. For this purpose, a novel operational matrix of fractional-order derivatives of Fibonacci polynomials was constructed and employed along with the application of the tau and collocation spectral methods. The convergence and error analysis of the suggested Fibonacci expansion were carefully investigated. Some numerical examples with comparisons are presented to ensure the efficiency, applicability and high accuracy of the proposed algorithms. Two accurate semi-analytic polynomial solutions for linear and nonlinear fractional differential equations are the result.

Chemical Equation Balancing.

Science.gov (United States)

Blakley, G. R.

1982-01-01

Reviews mathematical techniques for solving systems of homogeneous linear equations and demonstrates that the algebraic method of balancing chemical equations is a matter of solving a system of homogeneous linear equations. FORTRAN programs using this matrix method to chemical equation balancing are available from the author. (JN)

User's guide for SAMMY: a computer model for multilevel r-matrix fits to neutron data using Bayes' equations

International Nuclear Information System (INIS)

Larson, N.M.; Perey, F.G.

1980-11-01

A method is described for determining the parameters of a model from experimental data based upon the utilization of Bayes' theorem. This method has several advantages over the least-squares method as it is commonly used; one important advantage is that the assumptions under which the parameter values have been determined are more clearly evident than in many results based upon least squares. Bayes' method has been used to develop a computer code which can be utilized to analyze neutron cross-section data by means of the R-matrix theory. The required formulae from the R-matrix theory are presented, and the computer implementation of both Bayes' equations and R-matrix theory is described. Details about the computer code and compelte input/output information are given

A Matrix Method Based on the Fibonacci Polynomials to the Generalized Pantograph Equations with Functional Arguments

Directory of Open Access Journals (Sweden)

Ayşe Betül Koç

2014-01-01

Full Text Available A pseudospectral method based on the Fibonacci operational matrix is proposed to solve generalized pantograph equations with linear functional arguments. By using this method, approximate solutions of the problems are easily obtained in form of the truncated Fibonacci series. Some illustrative examples are given to verify the efficiency and effectiveness of the proposed method. Then, the numerical results are compared with other methods.

CRE Solvability, Nonlocal Symmetry and Exact Interaction Solutions of the Fifth-Order Modified Korteweg-de Vries Equation

Science.gov (United States)

Cheng, Wen-Guang; Qiu, De-Qin; Yu, Bo

2017-06-01

This paper is concerned with the fifth-order modified Korteweg-de Vries (fmKdV) equation. It is proved that the fmKdV equation is consistent Riccati expansion (CRE) solvable. Three special form of soliton-cnoidal wave interaction solutions are discussed analytically and shown graphically. Furthermore, based on the consistent tanh expansion (CTE) method, the nonlocal symmetry related to the consistent tanh expansion (CTE) is investigated, we also give the relationship between this kind of nonlocal symmetry and the residual symmetry which can be obtained with the truncated Painlevé method. We further study the spectral function symmetry and derive the Lax pair of the fmKdV equation. The residual symmetry can be localized to the Lie point symmetry of an enlarged system and the corresponding finite transformation group is computed. Supported by National Natural Science Foundation of China under Grant No. 11505090, and Research Award Foundation for Outstanding Young Scientists of Shandong Province under Grant No. BS2015SF009

Factorization and the synthesis of optimal feedback kernels for differential-delay systems

Science.gov (United States)

Milman, Mark M.; Scheid, Robert E.

1987-01-01

A combination of ideas from the theories of operator Riccati equations and Volterra factorizations leads to the derivation of a novel, relatively simple set of hyperbolic equations which characterize the optimal feedback kernel for the finite-time regulator problem for autonomous differential-delay systems. Analysis of these equations elucidates the underlying structure of the feedback kernel and leads to the development of fast and accurate numerical methods for its computation. Unlike traditional formulations based on the operator Riccati equation, the gain is characterized by means of classical solutions of the derived set of equations. This leads to the development of approximation schemes which are analogous to what has been accomplished for systems of ordinary differential equations with given initial conditions.

Efficient propagation of the hierarchical equations of motion using the matrix product state method

Science.gov (United States)

Shi, Qiang; Xu, Yang; Yan, Yaming; Xu, Meng

2018-05-01

We apply the matrix product state (MPS) method to propagate the hierarchical equations of motion (HEOM). It is shown that the MPS approximation works well in different type of problems, including boson and fermion baths. The MPS method based on the time-dependent variational principle is also found to be applicable to HEOM with over one thousand effective modes. Combining the flexibility of the HEOM in defining the effective modes and the efficiency of the MPS method thus may provide a promising tool in simulating quantum dynamics in condensed phases.

Optimal control of coupled parabolic-hyperbolic non-autonomous PDEs: infinite-dimensional state-space approach

Science.gov (United States)

Aksikas, I.; Moghadam, A. Alizadeh; Forbes, J. F.

2018-04-01

This paper deals with the design of an optimal state-feedback linear-quadratic (LQ) controller for a system of coupled parabolic-hypebolic non-autonomous partial differential equations (PDEs). The infinite-dimensional state space representation and the corresponding operator Riccati differential equation are used to solve the control problem. Dynamical properties of the coupled system of interest are analysed to guarantee the existence and uniqueness of the solution of the LQ-optimal control problem and also to guarantee the exponential stability of the closed-loop system. Thanks to the eigenvalues and eigenfunctions of the parabolic operator and also the fact that the hyperbolic-associated operator Riccati differential equation can be converted to a scalar Riccati PDE, an algorithm to solve the LQ control problem has been presented. The results are applied to a non-isothermal packed-bed catalytic reactor. The LQ optimal controller designed in the early portion of the paper is implemented for the original non-linear model. Numerical simulations are performed to show the controller performances.

A Fast Newton-Shamanskii Iteration for a Matrix Equation Arising from M/G/1-Type Markov Chains

Directory of Open Access Journals (Sweden)

Pei-Chang Guo

2017-01-01

Full Text Available For the nonlinear matrix equations arising in the analysis of M/G/1-type and GI/M/1-type Markov chains, the minimal nonnegative solution G or R can be found by Newton-like methods. We prove monotone convergence results for the Newton-Shamanskii iteration for this class of equations. Starting with zero initial guess or some other suitable initial guess, the Newton-Shamanskii iteration provides a monotonically increasing sequence of nonnegative matrices converging to the minimal nonnegative solution. A Schur decomposition method is used to accelerate the Newton-Shamanskii iteration. Numerical examples illustrate the effectiveness of the Newton-Shamanskii iteration.

Schwinger-Dyson loop equations as the w1+∞-like constraints for hermitian multi-matrix chain model at finite N

International Nuclear Information System (INIS)

Cheng, Yi-Xin

1992-01-01

The Schwinger-Dyson loop equations for the hermitian multi-matrix chain models at finite N, are derived from the Ward identities of the partition functional under the infinitesimal field transformations. The constraint operators W n (m) satisfy the w 1+∞ -like algebra up to a linear combination of the lower spin operators. We find that the all the higher spin constraints are reducible to the Virasoro-type constraints for all the matrix chain models. (author)

Straightening: existence, uniqueness and stability

Science.gov (United States)

Destrade, M.; Ogden, R. W.; Sgura, I.; Vergori, L.

2014-01-01

One of the least studied universal deformations of incompressible nonlinear elasticity, namely the straightening of a sector of a circular cylinder into a rectangular block, is revisited here and, in particular, issues of existence and stability are addressed. Particular attention is paid to the system of forces required to sustain the large static deformation, including by the application of end couples. The influence of geometric parameters and constitutive models on the appearance of wrinkles on the compressed face of the block is also studied. Different numerical methods for solving the incremental stability problem are compared and it is found that the impedance matrix method, based on the resolution of a matrix Riccati differential equation, is the more precise. PMID:24711723

A solver for General Unilateral Polynomial Matrix Equation with Second-Order Matrices Over Prime Finite Fields

Science.gov (United States)

Burtyka, Filipp

2018-03-01

The paper firstly considers the problem of finding solvents for arbitrary unilateral polynomial matrix equations with second-order matrices over prime finite fields from the practical point of view: we implement the solver for this problem. The solver’s algorithm has two step: the first is finding solvents, having Jordan Normal Form (JNF), the second is finding solvents among the rest matrices. The first step reduces to the finding roots of usual polynomials over finite fields, the second is essentially exhaustive search. The first step’s algorithms essentially use the polynomial matrices theory. We estimate the practical duration of computations using our software implementation (for example that one can’t construct unilateral matrix polynomial over finite field, having any predefined number of solvents) and answer some theoretically-valued questions.

Globally exponential stability condition of a class of neural networks with time-varying delays

International Nuclear Information System (INIS)

Liao, T.-L.; Yan, J.-J.; Cheng, C.-J.; Hwang, C.-C.

2005-01-01

In this Letter, the globally exponential stability for a class of neural networks including Hopfield neural networks and cellular neural networks with time-varying delays is investigated. Based on the Lyapunov stability method, a novel and less conservative exponential stability condition is derived. The condition is delay-dependent and easily applied only by checking the Hamiltonian matrix with no eigenvalues on the imaginary axis instead of directly solving an algebraic Riccati equation. Furthermore, the exponential stability degree is more easily assigned than those reported in the literature. Some examples are given to demonstrate validity and excellence of the presented stability condition herein

Integrable discretization s of derivative nonlinear Schroedinger equations

International Nuclear Information System (INIS)

Tsuchida, Takayuki

2002-01-01

We propose integrable discretizations of derivative nonlinear Schroedinger (DNLS) equations such as the Kaup-Newell equation, the Chen-Lee-Liu equation and the Gerdjikov-Ivanov equation by constructing Lax pairs. The discrete DNLS systems admit the reduction of complex conjugation between two dependent variables and possess bi-Hamiltonian structure. Through transformations of variables and reductions, we obtain novel integrable discretizations of the nonlinear Schroedinger (NLS), modified KdV (mKdV), mixed NLS, matrix NLS, matrix KdV, matrix mKdV, coupled NLS, coupled Hirota, coupled Sasa-Satsuma and Burgers equations. We also discuss integrable discretizations of the sine-Gordon equation, the massive Thirring model and their generalizations. (author)

Numerical study of Langevin equation in twisted Eguchi-Kawai model: distribution of eigenvalues of the plaquette matrix

International Nuclear Information System (INIS)

Migdal, A.A.; Polikarpov, M.I.; Veselov, A.I.; Yurov, V.P.

1983-01-01

The Langevin equation for the lattice theory with arbitrary gauge group is derived. The four-dimensional twisted Eguchi-Kawai model is investigated numerically. The results for the plaquette energy agree with those of the known Monte Carlo calculations. The new result is the distribution of eigenvalues of the plaquette matrix. In the strong coupling phase this distribution is smooth, whereas in the weak coupling phase a gap is clearly seen

Investigation of the existence and uniqueness of extremal and positive definite solutions of nonlinear matrix equations

Directory of Open Access Journals (Sweden)

Abdel-Shakoor M Sarhan

2016-05-01

Full Text Available Abstract We consider two nonlinear matrix equations X r ± ∑ i = 1 m A i ∗ X δ i A i = I $X^{r} \\pm \\sum_{i = 1}^{m} A_{i}^{*}X^{\\delta_{i}}A_{i} = I$ , where − 1 < δ i < 0 $- 1 < \\delta_{i} < 0$ , and r, m are positive integers. For the first equation (plus case, we prove the existence of positive definite solutions and extremal solutions. Two algorithms and proofs of their convergence to the extremal positive definite solutions are constructed. For the second equation (negative case, we prove the existence and the uniqueness of a positive definite solution. Moreover, the algorithm given in (Duan et al. in Linear Algebra Appl. 429:110-121, 2008 (actually, in (Shi et al. in Linear Multilinear Algebra 52:1-15, 2004 for r = 1 $r = 1$ is proved to be valid for any r. Numerical examples are given to illustrate the performance and effectiveness of all the constructed algorithms. In Appendix, we analyze the ordering on the positive cone P ( n ‾ $\\overline{P(n}$ .

Application of wavelets to singular integral scattering equations

International Nuclear Information System (INIS)

Kessler, B.M.; Payne, G.L.; Polyzou, W.N.

2004-01-01

The use of orthonormal wavelet basis functions for solving singular integral scattering equations is investigated. It is shown that these basis functions lead to sparse matrix equations which can be solved by iterative techniques. The scaling properties of wavelets are used to derive an efficient method for evaluating the singular integrals. The accuracy and efficiency of the wavelet transforms are demonstrated by solving the two-body T-matrix equation without partial wave projection. The resulting matrix equation which is characteristic of multiparticle integral scattering equations is found to provide an efficient method for obtaining accurate approximate solutions to the integral equation. These results indicate that wavelet transforms may provide a useful tool for studying few-body systems

A Numerical Scheme for Ordinary Differential Equations Having Time Varying and Nonlinear Coefficients Based on the State Transition Matrix

Science.gov (United States)

Bartels, Robert E.

2002-01-01

A variable order method of integrating initial value ordinary differential equations that is based on the state transition matrix has been developed. The method has been evaluated for linear time variant and nonlinear systems of equations. While it is more complex than most other methods, it produces exact solutions at arbitrary time step size when the time variation of the system can be modeled exactly by a polynomial. Solutions to several nonlinear problems exhibiting chaotic behavior have been computed. Accuracy of the method has been demonstrated by comparison with an exact solution and with solutions obtained by established methods.

Kraus operator solutions to a fermionic master equation describing a thermal bath and their matrix representation

Science.gov (United States)

Xiang-Guo, Meng; Ji-Suo, Wang; Hong-Yi, Fan; Cheng-Wei, Xia

2016-04-01

We solve the fermionic master equation for a thermal bath to obtain its explicit Kraus operator solutions via the fermionic state approach. The normalization condition of the Kraus operators is proved. The matrix representation for these solutions is obtained, which is incongruous with the result in the book completed by Nielsen and Chuang [Quantum Computation and Quantum Information, Cambridge University Press, 2000]. As especial cases, we also present the Kraus operator solutions to master equations for describing the amplitude-decay model and the diffusion process at finite temperature. Project supported by the National Natural Science Foundation of China (Grant No. 11347026), the Natural Science Foundation of Shandong Province, China (Grant Nos. ZR2013AM012 and ZR2012AM004), and the Research Fund for the Doctoral Program and Scientific Research Project of Liaocheng University, Shandong Province, China.

Nonlinear recurrent neural networks for finite-time solution of general time-varying linear matrix equations.

Science.gov (United States)

Xiao, Lin; Liao, Bolin; Li, Shuai; Chen, Ke

2018-02-01

In order to solve general time-varying linear matrix equations (LMEs) more efficiently, this paper proposes two nonlinear recurrent neural networks based on two nonlinear activation functions. According to Lyapunov theory, such two nonlinear recurrent neural networks are proved to be convergent within finite-time. Besides, by solving differential equation, the upper bounds of the finite convergence time are determined analytically. Compared with existing recurrent neural networks, the proposed two nonlinear recurrent neural networks have a better convergence property (i.e., the upper bound is lower), and thus the accurate solutions of general time-varying LMEs can be obtained with less time. At last, various different situations have been considered by setting different coefficient matrices of general time-varying LMEs and a great variety of computer simulations (including the application to robot manipulators) have been conducted to validate the better finite-time convergence of the proposed two nonlinear recurrent neural networks. Copyright © 2017 Elsevier Ltd. All rights reserved.

The Application Strategy of Iterative Solution Methodology to Matrix Equations in Hydraulic Solver Package, SPACE

International Nuclear Information System (INIS)

Na, Y. W.; Park, C. E.; Lee, S. Y.

2009-01-01

main object of this work is not to investigate the whole transient behavior of the models at hand but to focus on the behavior of numerical solutions part of the sparse asymmetric matrix equations in the transient of hydraulic system. It is outside of the scope of this work to improve the diagonal dominance or to pre-condition the matrix in the process of differencing and linearizing the governing equation, even though it is better to do it that way before applying the solver if there is any efficient way available

Exact perturbation theory of multiphoton processes at high intensities. [Schroedinger equation, perturbation theory, matrix

Energy Technology Data Exchange (ETDEWEB)

Faisal, F H.M. [Bielefeld Univ. (Germany, F.R.). Fakultaet fuer Physik

1976-06-11

In this work the perturbation theory for multiphoton processes at high intensities is investigated and it is described an analytical method of summing the perturbation series to extract the contribution from all terms that give rise to the absorption of N photons by an atomic system. The method is first applied to the solution of a simple model problem and the result is confirmed by direct integration of the model Schroedinger equation. The usual lowest (nonvanishing)-order perturbation-theoretical calculation is also carried out for this model to demonstrate explicitly that the full result correctly reproduces that of the lowest-order theory in the limit of low intensity. The method is then extended to the case of an atomic system with well-developed spectrum (e.g. H atom) and the N-photon T-matrix is derived in terms of a ''photon matrix'' asub(N), for which a three-term recurrence relation is established. Next, from the vantage point of the general result obtained here, A probe is made into the nature of several approximate nonperturbative solutions that have appeared in the literature in the past. It is shown here that their applicability is severely restricted by the requirement of the essential spectral degeneracy of the atomic system. Finally, appendix A outlines a prescription of computing the photon matrix asub(N), which (as in the usual lowest-order perturbation-theoretical calculation)requires a knowledge of the eigenfunctions and eigenvalues of the atomic Hamiltonian only.

Quantum mechanics in matrix form

CERN Document Server

Ludyk, Günter

2018-01-01

This book gives an introduction to quantum mechanics with the matrix method. Heisenberg's matrix mechanics is described in detail. The fundamental equations are derived by algebraic methods using matrix calculus. Only a brief description of Schrödinger's wave mechanics is given (in most books exclusively treated), to show their equivalence to Heisenberg's matrix  method. In the first part the historical development of Quantum theory by Planck, Bohr and Sommerfeld is sketched, followed by the ideas and methods of Heisenberg, Born and Jordan. Then Pauli's spin and exclusion principles are treated. Pauli's exclusion principle leads to the structure of atoms. Finally, Dirac´s relativistic quantum mechanics is shortly presented. Matrices and matrix equations are today easy to handle when implementing numerical algorithms using standard software as MAPLE and Mathematica.

Universal shocks in the Wishart random-matrix ensemble.

Science.gov (United States)

Blaizot, Jean-Paul; Nowak, Maciej A; Warchoł, Piotr

2013-05-01

We show that the derivative of the logarithm of the average characteristic polynomial of a diffusing Wishart matrix obeys an exact partial differential equation valid for an arbitrary value of N, the size of the matrix. In the large N limit, this equation generalizes the simple inviscid Burgers equation that has been obtained earlier for Hermitian or unitary matrices. The solution, through the method of characteristics, presents singularities that we relate to the precursors of shock formation in the Burgers equation. The finite N effects appear as a viscosity term in the Burgers equation. Using a scaling analysis of the complete equation for the characteristic polynomial, in the vicinity of the shocks, we recover in a simple way the universal Bessel oscillations (so-called hard-edge singularities) familiar in random-matrix theory.

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.

Delayed Stochastic Linear-Quadratic Control Problem and Related Applications

Directory of Open Access Journals (Sweden)

Li Chen

2012-01-01

stochastic differential equations (FBSDEs with Itô’s stochastic delay equations as forward equations and anticipated backward stochastic differential equations as backward equations. Especially, we present the optimal feedback regulator for the time delay system via a new type of Riccati equations and also apply to a population optimal control problem.

Matrix calculus

CERN Document Server

Bodewig, E

1959-01-01

Matrix Calculus, Second Revised and Enlarged Edition focuses on systematic calculation with the building blocks of a matrix and rows and columns, shunning the use of individual elements. The publication first offers information on vectors, matrices, further applications, measures of the magnitude of a matrix, and forms. The text then examines eigenvalues and exact solutions, including the characteristic equation, eigenrows, extremum properties of the eigenvalues, bounds for the eigenvalues, elementary divisors, and bounds for the determinant. The text ponders on approximate solutions, as well

Robust and scalable hierarchical matrix-based fast direct solver and preconditioner for the numerical solution of elliptic partial differential equations

KAUST Repository

Chavez, Gustavo Ivan

2017-07-10

This dissertation introduces a novel fast direct solver and preconditioner for the solution of block tridiagonal linear systems that arise from the discretization of elliptic partial differential equations on a Cartesian product mesh, such as the variable-coefficient Poisson equation, the convection-diffusion equation, and the wave Helmholtz equation in heterogeneous media. The algorithm extends the traditional cyclic reduction method with hierarchical matrix techniques. The resulting method exposes substantial concurrency, and its arithmetic operations and memory consumption grow only log-linearly with problem size, assuming bounded rank of off-diagonal matrix blocks, even for problems with arbitrary coefficient structure. The method can be used as a standalone direct solver with tunable accuracy, or as a black-box preconditioner in conjunction with Krylov methods. The challenges that distinguish this work from other thrusts in this active field are the hybrid distributed-shared parallelism that can demonstrate the algorithm at large-scale, full three-dimensionality, and the three stressors of the current state-of-the-art multigrid technology: high wavenumber Helmholtz (indefiniteness), high Reynolds convection (nonsymmetry), and high contrast diffusion (inhomogeneity). Numerical experiments corroborate the robustness, accuracy, and complexity claims and provide a baseline of the performance and memory footprint by comparisons with competing approaches such as the multigrid solver hypre, and the STRUMPACK implementation of the multifrontal factorization with hierarchically semi-separable matrices. The companion implementation can utilize many thousands of cores of Shaheen, KAUST\\'s Haswell-based Cray XC-40 supercomputer, and compares favorably with other implementations of hierarchical solvers in terms of time-to-solution and memory consumption.

Group-theoretical deduction of a dyadic Tamm-Dancoff equation by using a matrix-valued generator coordinate

International Nuclear Information System (INIS)

Nishiyama, Seiya; Morita, Hiroyuki; Ohnishi, Hiromasa

2004-01-01

The traditional Tamm-Dancoff (TD) method is one of the standard procedures for solving the Schroedinger equation of fermion many-body systems. However, it meets a serious difficulty when an instability occurs in the symmetry-adapted ground state of the independent particle approximation (IPA) and when the stable IPA ground state becomes of broken symmetry. If one uses the stable but broken symmetry IPA ground state as the starting approximation, TD wave functions also become of broken symmetry. On the contrary, if we start from a symmetry-adapted but unstable wave function, the convergence of the TD expansion becomes bad. Thus, the requirements of symmetry and rapid convergence are not in general compatible in the conventional TD expansion of the systems with strong collective correlations. Along the same line as Fukutome's, we give a group-theoretical deduction of a U(n) dyadic TD equation by using a matrix-valued generator coordinate

The harmonic oscillator and the position dependent mass Schroedinger equation: isospectral partners and factorization operators

International Nuclear Information System (INIS)

Morales, J.; Ovando, G.; Pena, J. J.

2010-01-01

One of the most important scientific contributions of Professor Marcos Moshinsky has been his study on the harmonic oscillator in quantum theory vis a vis the standard Schroedinger equation with constant mass [1]. However, a simple description of the motion of a particle interacting with an external environment such as happen in compositionally graded alloys consist of replacing the mass by the so-called effective mass that is in general variable and dependent on position. Therefore, honoring in memoriam Marcos Moshinsky, in this work we consider the position-dependent mass Schrodinger equations (PDMSE) for the harmonic oscillator potential model as former potential as well as with equi-spaced spectrum solutions, i.e. harmonic oscillator isospectral partners. To that purpose, the point canonical transformation method to convert a general second order differential equation (DE), of Sturm-Liouville type, into a Schroedinger-like standard equation is applied to the PDMSE. In that case, the former potential associated to the PDMSE and the potential involved in the Schroedinger-like standard equation are related through a Riccati-type relationship that includes the equivalent of the Witten superpotential to determine the exactly solvable positions-dependent mass distribution (PDMD)m(x). Even though the proposed approach is exemplified with the harmonic oscillator potential, the procedure is general and can be straightforwardly applied to other DEs.

Functional equations in matrix normed spaces

Indian Academy of Sciences (India)

The abstract characterization given for linear spaces of bounded Hilbert space operators in terms of ... effect on operator algebra theory (see [12]). .... of functional equations for the proof of new fixed point theorems with applications. By.

Numerical algebra, matrix theory, differential-algebraic equations and control theory festschrift in honor of Volker Mehrmann

CERN Document Server

Bollhöfer, Matthias; Kressner, Daniel; Mehl, Christian; Stykel, Tatjana

2015-01-01

This edited volume highlights the scientific contributions of Volker Mehrmann, a leading expert in the area of numerical (linear) algebra, matrix theory, differential-algebraic equations and control theory. These mathematical research areas are strongly related and often occur in the same real-world applications. The main areas where such applications emerge are computational engineering and sciences, but increasingly also social sciences and economics. This book also reflects some of Volker Mehrmann's major career stages. Starting out working in the areas of numerical linear algebra (his first full professorship at TU Chemnitz was in "Numerical Algebra," hence the title of the book) and matrix theory, Volker Mehrmann has made significant contributions to these areas ever since. The highlights of these are discussed in Parts I and II of the present book. Often the development of new algorithms in numerical linear algebra is motivated by problems in system and control theory. These and his later major work on ...

Real-time optical laboratory solution of parabolic differential equations

Science.gov (United States)

Casasent, David; Jackson, James

1988-01-01

An optical laboratory matrix-vector processor is used to solve parabolic differential equations (the transient diffusion equation with two space variables and time) by an explicit algorithm. This includes optical matrix-vector nonbase-2 encoded laboratory data, the combination of nonbase-2 and frequency-multiplexed data on such processors, a high-accuracy optical laboratory solution of a partial differential equation, new data partitioning techniques, and a discussion of a multiprocessor optical matrix-vector architecture.

Optimal reduction of flexible dynamic system

International Nuclear Information System (INIS)

Jankovic, J.

1994-01-01

Dynamic system reduction is basic procedure in various problems of active control synthesis of flexible structures. In this paper is presented direct method for system reduction by explicit extraction of modes included in reduced model form. Criterion for optimal system discrete approximation in synthesis reduced dynamic model is also presented. Subjected method of system decomposition is discussed in relation to the Schur method of solving matrix algebraic Riccati equation as condition for system reduction. By using exposed method procedure of flexible system reduction in addition with corresponding example is presented. Shown procedure is powerful in problems of active control synthesis of flexible system vibrations

Extended biorthogonal matrix polynomials

Directory of Open Access Journals (Sweden)

Ayman Shehata

2017-01-01

Full Text Available The pair of biorthogonal matrix polynomials for commutative matrices were first introduced by Varma and Tasdelen in [22]. The main aim of this paper is to extend the properties of the pair of biorthogonal matrix polynomials of Varma and Tasdelen and certain generating matrix functions, finite series, some matrix recurrence relations, several important properties of matrix differential recurrence relations, biorthogonality relations and matrix differential equation for the pair of biorthogonal matrix polynomials J(A,B n (x, k and K(A,B n (x, k are discussed. For the matrix polynomials J(A,B n (x, k, various families of bilinear and bilateral generating matrix functions are constructed in the sequel.

Unitarity or asymptotic completeness equations and analytic structure of the S matrix and Green functions

International Nuclear Information System (INIS)

Iagolnitzer, D.

1983-11-01

Recent axiomatic results on the (non holonomic) analytic structure of the multiparticle S matrix and Green functions are reviewed and related general conjectures are described: (i) formal expansions of Green functions in terms of (holonomic) Feynman-type integrals in which each vertex represents an irreducible kernel, and (ii) ''graph by graph unitarity'' and other discontinuity formulae of the latter. These conjectures are closely linked with unitarity or asymptotic completeness equations, which they yield in a formal sense. In constructive field theory, a direct proof of the first conjecture (together with an independent proof of the second) would thus imply, as a first step, asymptotic completeness in that sense

Parametric optimal control of uncertain systems under an optimistic value criterion

Science.gov (United States)

Li, Bo; Zhu, Yuanguo

2018-01-01

It is well known that the optimal control of a linear quadratic model is characterized by the solution of a Riccati differential equation. In many cases, the corresponding Riccati differential equation cannot be solved exactly such that the optimal feedback control may be a complex time-oriented function. In this article, a parametric optimal control problem of an uncertain linear quadratic model under an optimistic value criterion is considered for simplifying the expression of optimal control. Based on the equation of optimality for the uncertain optimal control problem, an approximation method is presented to solve it. As an application, a two-spool turbofan engine optimal control problem is given to show the utility of the proposed model and the efficiency of the presented approximation method.

Matrix Krylov subspace methods for image restoration

Directory of Open Access Journals (Sweden)

khalide jbilou

2015-09-01

Full Text Available In the present paper, we consider some matrix Krylov subspace methods for solving ill-posed linear matrix equations and in those problems coming from the restoration of blurred and noisy images. Applying the well known Tikhonov regularization procedure leads to a Sylvester matrix equation depending the Tikhonov regularized parameter. We apply the matrix versions of the well known Krylov subspace methods, namely the Least Squared (LSQR and the conjugate gradient (CG methods to get approximate solutions representing the restored images. Some numerical tests are presented to show the effectiveness of the proposed methods.

New and More General Rational Formal Solutions to (2+1)-Dimensional Toda System

International Nuclear Information System (INIS)

Bai Chenglin

2007-01-01

With the aid of computerized symbolic computation and Riccati equation rational expansion approach, some new and more general rational formal solutions to (2+1)-dimensional Toda system are obtained. The method used here can also be applied to solve other nonlinear differential-difference equation or equations.

Liouville equation of relativistic charged fermion

International Nuclear Information System (INIS)

Wang Renchuan; Zhu Dongpei; Huang Zhuoran; Ko Che-ming

1991-01-01

As a form of density martrix, the Wigner function is the distribution in quantum phase space. It is a 2 X 2 matrix function when one uses it to describe the non-relativistic fermion. While describing the relativistic fermion, it is usually represented by 4 x 4 matrix function. In this paper authors obtain a Wigner function for the relativistic fermion in the form of 2 x 2 matrix, and the Liouville equation satisfied by the Wigner function. this equivalent to the Dirac equation of changed fermion in QED. The equation is also equivalent to the Dirac equation in the Walecka model applied to the intermediate energy nuclear collision while the nucleon is coupled to the vector meson only (or taking mean field approximation for the scalar meson). Authors prove that the 2 x 2 Wigner function completely describes the quantum system just the same as the relativistic fermion wave function. All the information about the observables can be obtained with above Wigner function

Associative Yang-Baxter equation for quantum (semi-)dynamical R-matrices

International Nuclear Information System (INIS)

Sechin, Ivan; Zotov, Andrei

2016-01-01

In this paper we propose versions of the associative Yang-Baxter equation and higher order R-matrix identities which can be applied to quantum dynamical R-matrices. As is known quantum non-dynamical R-matrices of Baxter-Belavin type satisfy this equation. Together with unitarity condition and skew-symmetry it provides the quantum Yang-Baxter equation and a set of identities useful for different applications in integrable systems. The dynamical R-matrices satisfy the Gervais-Neveu-Felder (or dynamical Yang-Baxter) equation. Relation between the dynamical and non-dynamical cases is described by the IRF (interaction-round-a-face)-Vertex transformation. An alternative approach to quantum (semi-)dynamical R-matrices and related quantum algebras was suggested by Arutyunov, Chekhov, and Frolov (ACF) in their study of the quantum Ruijsenaars-Schneider model. The purpose of this paper is twofold. First, we prove that the ACF elliptic R-matrix satisfies the associative Yang-Baxter equation with shifted spectral parameters. Second, we directly prove a simple relation of the IRF-Vertex type between the Baxter-Belavin and the ACF elliptic R-matrices predicted previously by Avan and Rollet. It provides the higher order R-matrix identities and an explanation of the obtained equations through those for non-dynamical R-matrices. As a by-product we also get an interpretation of the intertwining transformation as matrix extension of scalar theta function likewise R-matrix is interpreted as matrix extension of the Kronecker function. Relations to the Gervais-Neveu-Felder equation and identities for the Felder’s elliptic R-matrix are also discussed.

Integrable systems of partial differential equations determined by structure equations and Lax pair

International Nuclear Information System (INIS)

Bracken, Paul

2010-01-01

It is shown how a system of evolution equations can be developed both from the structure equations of a submanifold embedded in three-space as well as from a matrix SO(6) Lax pair. The two systems obtained this way correspond exactly when a constraint equation is selected and imposed on the system of equations. This allows for the possibility of selecting the coefficients in the second fundamental form in a general way.

Parametric Approach to Trajectory Tracking Control of Robot Manipulators

Directory of Open Access Journals (Sweden)

Shijie Zhang

2013-01-01

Full Text Available The mathematic description of the trajectory of robot manipulators with the optimal trajectory tracking problem is formulated as an optimal control problem, and a parametric approach is proposed for the optimal trajectory tracking control problem. The optimal control problem is first solved as an open loop optimal control problem by using a time scaling transform and the control parameterization method. Then, by virtue of the relationship between the optimal open loop control and the optimal closed loop control along the optimal trajectory, a practical method is presented to calculate an approximate optimal feedback gain matrix, without having to solve an optimal control problem involving the complex Riccati-like matrix differential equation coupled with the original system dynamics. Simulation results of 2-link robot manipulator are presented to show the effectiveness of the proposed method.

Direct solution of the biharmonic equation on rectangular regions and the Poisson equation on irregular regions

International Nuclear Information System (INIS)

Buzbee, B.L.; Dorr, F.W.

1974-01-01

The discrete biharmonic equation on a rectangular region and the discrete Poisson equation on an irregular region can be treated as modifications to matrix problems with very special structure. It is shown how to use the direct method of matrix decomposition to formulate an effective numerical algorithm for these problems. For typical applications the operation count is O(N 3 ) for an N x N grid. Numerical comparisons with other techniques are included. (U.S.)

Massive Asynchronous Parallelization of Sparse Matrix Factorizations

Energy Technology Data Exchange (ETDEWEB)

Chow, Edmond [Georgia Inst. of Technology, Atlanta, GA (United States)

2018-01-08

Solving sparse problems is at the core of many DOE computational science applications. We focus on the challenge of developing sparse algorithms that can fully exploit the parallelism in extreme-scale computing systems, in particular systems with massive numbers of cores per node. Our approach is to express a sparse matrix factorization as a large number of bilinear constraint equations, and then solving these equations via an asynchronous iterative method. The unknowns in these equations are the matrix entries of the factorization that is desired.

An active vibration control method of bridge structures by the ...

African Journals Online (AJOL)

The 'Linearization of Feedback Gain Matrix' (LFGM) method proposed by the author to achieve this objective is obtained by activating some useful components of the Riccati matrix while others are strongly attenuated according to the expected output-results and requirements of control mechanism. The present algorithm is ...

A bias correction for covariance estimators to improve inference with generalized estimating equations that use an unstructured correlation matrix.

Science.gov (United States)

Westgate, Philip M

2013-07-20

Generalized estimating equations (GEEs) are routinely used for the marginal analysis of correlated data. The efficiency of GEE depends on how closely the working covariance structure resembles the true structure, and therefore accurate modeling of the working correlation of the data is important. A popular approach is the use of an unstructured working correlation matrix, as it is not as restrictive as simpler structures such as exchangeable and AR-1 and thus can theoretically improve efficiency. However, because of the potential for having to estimate a large number of correlation parameters, variances of regression parameter estimates can be larger than theoretically expected when utilizing the unstructured working correlation matrix. Therefore, standard error estimates can be negatively biased. To account for this additional finite-sample variability, we derive a bias correction that can be applied to typical estimators of the covariance matrix of parameter estimates. Via simulation and in application to a longitudinal study, we show that our proposed correction improves standard error estimation and statistical inference. Copyright © 2012 John Wiley & Sons, Ltd.

Darboux transformations and linear parabolic partial differential equations

International Nuclear Information System (INIS)

Arrigo, Daniel J.; Hickling, Fred

2002-01-01

Solutions for a class of linear parabolic partial differential equation are provided. These solutions are obtained by first solving a system of (n+1) nonlinear partial differential equations. This system arises as the coefficients of a Darboux transformation and is equivalent to a matrix Burgers' equation. This matrix equation is solved using a generalized Hopf-Cole transformation. The solutions for the original equation are given in terms of solutions of the heat equation. These results are applied to the (1+1)-dimensional Schroedinger equation where all bound state solutions are obtained for a 2n-parameter family of potentials. As a special case, the solutions for integral members of the regular and modified Poeschl-Teller potentials are recovered. (author). Letter-to-the-editor

On Triple Product and Rational Solutions of Yang—Baxter Equation

International Nuclear Information System (INIS)

Zhang Chun-Hong; Jia Xiao-Yu; Li Min-Li; Wu Ke; Zhao Wei-Zhong

2014-01-01

The Yang—Baxter equation is reinvestigated in the framework of triple system. By requiring the rational R matrix of the Yang—Baxter equation satisfying the generalized Filippov condition, we derive a relation with respect to the rational R matrix. Moreover the case of the super Yang—Baxter equation is also investigated. (general)

Discrete Ordinates Approximations to the First- and Second-Order Radiation Transport Equations

International Nuclear Information System (INIS)

FAN, WESLEY C.; DRUMM, CLIFTON R.; POWELL, JENNIFER L. email wcfan@sandia.gov

2002-01-01

The conventional discrete ordinates approximation to the Boltzmann transport equation can be described in a matrix form. Specifically, the within-group scattering integral can be represented by three components: a moment-to-discrete matrix, a scattering cross-section matrix and a discrete-to-moment matrix. Using and extending these entities, we derive and summarize the matrix representations of the second-order transport equations

Discrete Ordinates Approximations to the First- and Second-Order Radiation Transport Equations

CERN Document Server

Fan, W C; Powell, J L

2002-01-01

The conventional discrete ordinates approximation to the Boltzmann transport equation can be described in a matrix form. Specifically, the within-group scattering integral can be represented by three components: a moment-to-discrete matrix, a scattering cross-section matrix and a discrete-to-moment matrix. Using and extending these entities, we derive and summarize the matrix representations of the second-order transport equations.

Soliton solutions for ABS lattice equations: I. Cauchy matrix approach

Science.gov (United States)

Nijhoff, Frank; Atkinson, James; Hietarinta, Jarmo

2009-10-01

In recent years there have been new insights into the integrability of quadrilateral lattice equations, i.e. partial difference equations which are the natural discrete analogues of integrable partial differential equations in 1+1 dimensions. In the scalar (i.e. single-field) case, there now exist classification results by Adler, Bobenko and Suris (ABS) leading to some new examples in addition to the lattice equations 'of KdV type' that were known since the late 1970s and early 1980s. In this paper, we review the construction of soliton solutions for the KdV-type lattice equations and use those results to construct N-soliton solutions for all lattice equations in the ABS list except for the elliptic case of Q4, which is left to a separate treatment.

Finite Difference Schemes as Algebraic Correspondences between Layers

Science.gov (United States)

Malykh, Mikhail; Sevastianov, Leonid

2018-02-01

For some differential equations, especially for Riccati equation, new finite difference schemes are suggested. These schemes define protective correspondences between the layers. Calculation using these schemes can be extended to the area beyond movable singularities of exact solution without any error accumulation.

Channel Estimation and Information Symbol Detection for DS-UWB Communication Systems

Directory of Open Access Journals (Sweden)

Wei Wang

2014-01-01

estimation, the one-step predictor of information symbol is used and the estimation error is also considered as a multiplicative noise. The solutions to the above two problems are obtained by solving a couple of Riccati equations together with two Lyapunov equations.

Statistically and Computationally Efficient Estimating Equations for Large Spatial Datasets

KAUST Repository

Sun, Ying; Stein, Michael L.

2014-01-01

For Gaussian process models, likelihood based methods are often difficult to use with large irregularly spaced spatial datasets, because exact calculations of the likelihood for n observations require O(n3) operations and O(n2) memory. Various approximation methods have been developed to address the computational difficulties. In this paper, we propose new unbiased estimating equations based on score equation approximations that are both computationally and statistically efficient. We replace the inverse covariance matrix that appears in the score equations by a sparse matrix to approximate the quadratic forms, then set the resulting quadratic forms equal to their expected values to obtain unbiased estimating equations. The sparse matrix is constructed by a sparse inverse Cholesky approach to approximate the inverse covariance matrix. The statistical efficiency of the resulting unbiased estimating equations are evaluated both in theory and by numerical studies. Our methods are applied to nearly 90,000 satellite-based measurements of water vapor levels over a region in the Southeast Pacific Ocean.

Statistically and Computationally Efficient Estimating Equations for Large Spatial Datasets

KAUST Repository

Sun, Ying

2014-11-07

For Gaussian process models, likelihood based methods are often difficult to use with large irregularly spaced spatial datasets, because exact calculations of the likelihood for n observations require O(n3) operations and O(n2) memory. Various approximation methods have been developed to address the computational difficulties. In this paper, we propose new unbiased estimating equations based on score equation approximations that are both computationally and statistically efficient. We replace the inverse covariance matrix that appears in the score equations by a sparse matrix to approximate the quadratic forms, then set the resulting quadratic forms equal to their expected values to obtain unbiased estimating equations. The sparse matrix is constructed by a sparse inverse Cholesky approach to approximate the inverse covariance matrix. The statistical efficiency of the resulting unbiased estimating equations are evaluated both in theory and by numerical studies. Our methods are applied to nearly 90,000 satellite-based measurements of water vapor levels over a region in the Southeast Pacific Ocean.

Analytic solutions of the multigroup space-time reactor kinetics equations

International Nuclear Information System (INIS)

Lee, C.E.; Rottler, S.

1986-01-01

The development of analytical and numerical solutions to the reactor kinetics equations is reviewed. Analytic solutions of the multigroup space-time reactor kinetics equations are developed for bare and reflected slabs and spherical reactors for zero flux, zero current and extrapolated endpoint boundary conditions. The material properties of the reactors are assumed constant in space and time, but spatially-dependent source terms and initial conditions are investigated. The system of partial differential equations is reduced to a set of linear ordinary differential equations by the Laplace transform method. These equations are solved by matrix Green's functions yielding a general matrix solution for the neutron flux and precursor concentration in the Laplace transform space. The detailed pole structure of the Laplace transform matrix solutions is investigated. The temporally- and spatially-dependent solutions are determined from the inverse Laplace transform using the Cauchy residue theorem, the theorem of Frobenius, a knowledge of the detailed pole structure and matrix operators. (author)

Bandwidth Optimization of Normal Equation Matrix in Bundle Block Adjustment in Multi-baseline Rotational Photography

Directory of Open Access Journals (Sweden)

WANG Xiang

2016-02-01

Full Text Available A new bandwidth optimization method of normal equation matrix in bundle block adjustment in multi-baseline rotational close range photography by image index re-sorting is proposed. The equivalent exposure station of each image is calculated by its object space coverage and the relationship with other adjacent images. Then, according to the coordinate relations between equivalent exposure stations, new logical indices of all images are computed, based on which, the optimized bandwidth value can be obtained. Experimental results show that the bandwidth determined by our proposed method is significantly better than its original value, thus the operational efficiency, as well as the memory consumption of multi-baseline rotational close range photography in real-data applications, is optimized to a certain extent.

Highly efficient parallel direct solver for solving dense complex matrix equations from method of moments

Directory of Open Access Journals (Sweden)

Yan Chen

2017-03-01

Full Text Available Based on the vectorised and cache optimised kernel, a parallel lower upper decomposition with a novel communication avoiding pivoting scheme is developed to solve dense complex matrix equations generated by the method of moments. The fine-grain data rearrangement and assembler instructions are adopted to reduce memory accessing times and improve CPU cache utilisation, which also facilitate vectorisation of the code. Through grouping processes in a binary tree, a parallel pivoting scheme is designed to optimise the communication pattern and thus reduces the solving time of the proposed solver. Two large electromagnetic radiation problems are solved on two supercomputers, respectively, and the numerical results demonstrate that the proposed method outperforms those in open source and commercial libraries.

Flavored quantum Boltzmann equations

International Nuclear Information System (INIS)

Cirigliano, Vincenzo; Lee, Christopher; Ramsey-Musolf, Michael J.; Tulin, Sean

2010-01-01

We derive from first principles, using nonequilibrium field theory, the quantum Boltzmann equations that describe the dynamics of flavor oscillations, collisions, and a time-dependent mass matrix in the early universe. Working to leading nontrivial order in ratios of relevant time scales, we study in detail a toy model for weak-scale baryogenesis: two scalar species that mix through a slowly varying time-dependent and CP-violating mass matrix, and interact with a thermal bath. This model clearly illustrates how the CP asymmetry arises through coherent flavor oscillations in a nontrivial background. We solve the Boltzmann equations numerically for the density matrices, investigating the impact of collisions in various regimes.

A numerical algorithm to find all feedback Nash equilibria in scalar affine quadratic differential games

NARCIS (Netherlands)

Engwerda, Jacob

2015-01-01

This note deals with solving scalar coupled algebraic Riccati equations. These equations arise in finding linear feedback Nash equilibria of the scalar N-player affine quadratic differential game. A numerical procedure is provided to compute all the stabilizing solutions. The main idea is to

Inverse scattering transform for the time dependent Schroedinger equation with applications to the KPI equation

Energy Technology Data Exchange (ETDEWEB)

Xin, Zhou [Wisconsin Univ., Madison (USA). Dept. of Mathematics

1990-03-01

For the direct-inverse scattering transform of the time dependent Schroedinger equation, rigorous results are obtained based on an operator-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution. (orig.).

Inverse scattering transform for the time dependent Schroedinger equation with applications to the KPI equation

International Nuclear Information System (INIS)

Zhou Xin

1990-01-01

For the direct-inverse scattering transform of the time dependent Schroedinger equation, rigorous results are obtained based on an operator-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution. (orig.)

A Numerical Approximation Framework for the Stochastic Linear Quadratic Regulator on Hilbert Spaces

Energy Technology Data Exchange (ETDEWEB)

Levajković, Tijana, E-mail: tijana.levajkovic@uibk.ac.at, E-mail: t.levajkovic@sf.bg.ac.rs; Mena, Hermann, E-mail: hermann.mena@uibk.ac.at [University of Innsbruck, Department of Mathematics (Austria); Tuffaha, Amjad, E-mail: atufaha@aus.edu [American University of Sharjah, Department of Mathematics (United Arab Emirates)

2017-06-15

We present an approximation framework for computing the solution of the stochastic linear quadratic control problem on Hilbert spaces. We focus on the finite horizon case and the related differential Riccati equations (DREs). Our approximation framework is concerned with the so-called “singular estimate control systems” (Lasiecka in Optimal control problems and Riccati equations for systems with unbounded controls and partially analytic generators: applications to boundary and point control problems, 2004) which model certain coupled systems of parabolic/hyperbolic mixed partial differential equations with boundary or point control. We prove that the solutions of the approximate finite-dimensional DREs converge to the solution of the infinite-dimensional DRE. In addition, we prove that the optimal state and control of the approximate finite-dimensional problem converge to the optimal state and control of the corresponding infinite-dimensional problem.

An efficient technique for the point reactor kinetics equations with Newtonian temperature feedback effects

International Nuclear Information System (INIS)

Nahla, Abdallah A.

2011-01-01

Highlights: → An efficient technique for the nonlinear reactor kinetics equations is presented. → This method is based on Backward Euler or Crank Nicholson and fundamental matrix. → Stability of efficient technique is defined and discussed. → This method is applied to point kinetics equations of six-groups of delayed neutrons. → Step, ramp, sinusoidal and temperature feedback reactivities are discussed. - Abstract: The point reactor kinetics equations of multi-group of delayed neutrons in the presence Newtonian temperature feedback effects are a system of stiff nonlinear ordinary differential equations which have not any exact analytical solution. The efficient technique for this nonlinear system is based on changing this nonlinear system to a linear system by the predicted value of reactivity and solving this linear system using the fundamental matrix of the homogenous linear differential equations. The nonlinear point reactor kinetics equations are rewritten in the matrix form. The solution of this matrix form is introduced. This solution contains the exponential function of a variable coefficient matrix. This coefficient matrix contains the unknown variable, reactivity. The predicted values of reactivity in the explicit form are determined replacing the exponential function of the coefficient matrix by two kinds, Backward Euler and Crank Nicholson, of the rational approximations. The nonlinear point kinetics equations changed to a linear system of the homogenous differential equations. The fundamental matrix of this linear system is calculated using the eigenvalues and the corresponding eigenvectors of the coefficient matrix. Stability of the efficient technique is defined and discussed. The efficient technique is applied to the point kinetics equations of six-groups of delayed neutrons with step, ramp, sinusoidal and the temperature feedback reactivities. The results of these efficient techniques are compared with the traditional methods.

The development of the deterministic nonlinear PDEs in particle physics to stochastic case

Science.gov (United States)

Abdelrahman, Mahmoud A. E.; Sohaly, M. A.

2018-06-01

In the present work, accuracy method called, Riccati-Bernoulli Sub-ODE technique is used for solving the deterministic and stochastic case of the Phi-4 equation and the nonlinear Foam Drainage equation. Also, the control on the randomness input is studied for stability stochastic process solution.

New Families of Rational Form Solitary Wave Solutions to (2+1)-Dimensional Broer-Kaup-Kupershmidt System

International Nuclear Information System (INIS)

Wang Qi; Li Biao; Zhang Hongqing; Chen Yong

2005-01-01

Taking the (2+1)-dimensional Broer-Kaup-Kupershmidt system as a simple example, some families of rational form solitary wave solutions, triangular periodic wave solutions, and rational wave solutions are constructed by using the Riccati equation rational expansion method presented by us. The method can also be applied to solve more nonlinear partial differential equation or equations.

Stability theory for dynamic equations on time scales

CERN Document Server

Martynyuk, Anatoly A

2016-01-01

This monograph is a first in the world to present three approaches for stability analysis of solutions of dynamic equations. The first approach is based on the application of dynamic integral inequalities and the fundamental matrix of solutions of linear approximation of dynamic equations. The second is based on the generalization of the direct Lyapunovs method for equations on time scales, using scalar, vector and matrix-valued auxiliary functions. The third approach is the application of auxiliary functions (scalar, vector, or matrix-valued ones) in combination with differential dynamic inequalities. This is an alternative comparison method, developed for time continuous and time discrete systems. In recent decades, automatic control theory in the study of air- and spacecraft dynamics and in other areas of modern applied mathematics has encountered problems in the analysis of the behavior of solutions of time continuous-discrete linear and/or nonlinear equations of perturbed motion. In the book “Men of Ma...

Cluster-enriched Yang-Baxter equation from SUSY gauge theories

Science.gov (United States)

Yamazaki, Masahito

2018-04-01

We propose a new generalization of the Yang-Baxter equation, where the R-matrix depends on cluster y-variables in addition to the spectral parameters. We point out that we can construct solutions to this new equation from the recently found correspondence between Yang-Baxter equations and supersymmetric gauge theories. The S^2 partition function of a certain 2d N=(2,2) quiver gauge theory gives an R-matrix, whereas its FI parameters can be identified with the cluster y-variables.

Generating Nice Linear Systems for Matrix Gaussian Elimination

Science.gov (United States)

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…

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

Science.gov (United States)

Frisch, H.

1994-01-01

computations: the eigenvalues and eigenvectors of real matrices; the relative stability of a given matrix; matrix factorization; the solution of linear constant coefficient vector-matrix algebraic equations; the controllability properties of a linear time-invariant system; the steady-state covariance matrix of an open-loop stable system forced by white noise; and the transient response of continuous linear time-invariant systems. The control law design routines of ORACLS implement some of the more common techniques of time-invariant LQG methodology. For the finite-duration optimal linear regulator problem with noise-free measurements, continuous dynamics, and integral performance index, a routine is provided which implements the negative exponential method for finding both the transient and steady-state solutions to the matrix Riccati equation. For the discrete version of this problem, the method of backwards differencing is applied to find the solutions to the discrete Riccati equation. A routine is also included to solve the steady-state Riccati equation by the Newton algorithms described by Klein, for continuous problems, and by Hewer, for discrete problems. Another routine calculates the prefilter gain to eliminate control state cross-product terms in the quadratic performance index and the weighting matrices for the sampled data optimal linear regulator problem. For cases with measurement noise, duality theory and optimal regulator algorithms are used to calculate solutions to the continuous and discrete Kalman-Bucy filter problems. Finally, routines are included to implement the continuous and discrete forms of the explicit (model-in-the-system) and implicit (model-in-the-performance-index) model following theory. These routines generate linear control laws which cause the output of a dynamic time-invariant system to track the output of a prescribed model. In order to apply ORACLS, the user must write an executive (driver) program which inputs the problem coefficients

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

Science.gov (United States)

Armstrong, E. S.

1994-01-01

computations: the eigenvalues and eigenvectors of real matrices; the relative stability of a given matrix; matrix factorization; the solution of linear constant coefficient vector-matrix algebraic equations; the controllability properties of a linear time-invariant system; the steady-state covariance matrix of an open-loop stable system forced by white noise; and the transient response of continuous linear time-invariant systems. The control law design routines of ORACLS implement some of the more common techniques of time-invariant LQG methodology. For the finite-duration optimal linear regulator problem with noise-free measurements, continuous dynamics, and integral performance index, a routine is provided which implements the negative exponential method for finding both the transient and steady-state solutions to the matrix Riccati equation. For the discrete version of this problem, the method of backwards differencing is applied to find the solutions to the discrete Riccati equation. A routine is also included to solve the steady-state Riccati equation by the Newton algorithms described by Klein, for continuous problems, and by Hewer, for discrete problems. Another routine calculates the prefilter gain to eliminate control state cross-product terms in the quadratic performance index and the weighting matrices for the sampled data optimal linear regulator problem. For cases with measurement noise, duality theory and optimal regulator algorithms are used to calculate solutions to the continuous and discrete Kalman-Bucy filter problems. Finally, routines are included to implement the continuous and discrete forms of the explicit (model-in-the-system) and implicit (model-in-the-performance-index) model following theory. These routines generate linear control laws which cause the output of a dynamic time-invariant system to track the output of a prescribed model. In order to apply ORACLS, the user must write an executive (driver) program which inputs the problem coefficients

Rational parametrisation of normalised Stiefel manifolds, and explicit non-'t Hooft solutions of the Atiyah-Drinfeld-Hitchin-Manin instanton matrix equations for Sp(n)

International Nuclear Information System (INIS)

McCarthy, P.J.

1981-01-01

It is proved that normalised Stiefel manifolds admit a rational parametrisation which generalises Cayley's parametrisation of the unitary groups. Applying (the quaternionic case of) this parametrisation to the Atiyah-Drinfeld-Hitchin-Manin (ADHM) instanton matrix equations, large families of new explicit rational solutions emerge. In particular, new explicit non-'t Hooft solutions are presented. (orig.)

How Should Equation Balancing Be Taught?

Science.gov (United States)

Porter, Spencer K.

1985-01-01

Matrix methods and oxidation-number methods are currently advocated and used for balancing equations. This article shows how balancing equations can be introduced by a third method which is related to a fundamental principle, is easy to learn, and is powerful in its application. (JN)

Modeling and Simulation of Matrix Converter

DEFF Research Database (Denmark)

Liu, Fu-rong; Klumpner, Christian; Blaabjerg, Frede

2005-01-01

This paper discusses the modeling and simulation of matrix converter. Two models of matrix converter are presented: one is based on indirect space vector modulation and the other is based on power balance equation. The basis of these two models is• given and the process on modeling is introduced...

Hierarchical matrix techniques for the solution of elliptic equations

KAUST Repository

Chávez, Gustavo

2014-05-04

Hierarchical matrix approximations are a promising tool for approximating low-rank matrices given the compactness of their representation and the economy of the operations between them. Integral and differential operators have been the major applications of this technology, but they can be applied into other areas where low-rank properties exist. Such is the case of the Block Cyclic Reduction algorithm, which is used as a direct solver for the constant-coefficient Poisson quation. We explore the variable-coefficient case, also using Block Cyclic reduction, with the addition of Hierarchical Matrices to represent matrix blocks, hence improving the otherwise O(N2) algorithm, into an efficient O(N) algorithm.

Inverse scattering transform for the time dependent Schrödinger equation with applications to the KPI equation

Science.gov (United States)

Zhou, Xin

1990-03-01

For the direct-inverse scattering transform of the time dependent Schrödinger equation, rigorous results are obtained based on an opertor-triangular-factorization approach. By viewing the equation as a first order operator equation, similar results as for the first order n x n matrix system are obtained. The nonlocal Riemann-Hilbert problem for inverse scattering is shown to have solution.

The dressing factor and crossing equations

International Nuclear Information System (INIS)

Arutyunov, Gleb; Frolov, Sergey

2009-01-01

We utilize the DHM integral representation for the BES dressing factor of the world-sheet S-matrix of the AdS 5 x S 5 light-cone string theory, and the crossing equations to fix the principal branch of the dressing factor on the rapidity torus. The results obtained are further used, in conjunction with the fusion procedure, to determine the bound state dressing factor of the mirror theory. We convincingly demonstrate that the mirror bound state S-matrix found in this way does not depend on the internal structure of a bound state solution employed in the fusion procedure. This welcome feature is in perfect parallel to string theory, where the corresponding bound state S-matrix has no bearing on bound state constituent particles as well. The mirror bound state S-matrix we found provides the final missing piece in setting up the TBA equations for the AdS 5 x S 5 mirror theory.

From polymers to quantum gravity: Triple-scaling in rectangular random matrix models

International Nuclear Information System (INIS)

Myers, R.C.; Periwal, V.

1993-01-01

Rectangular NxM matrix models can be solved in several qualitatively distinct large-N limits, since two independent parameters govern the size of the matrix. Regarded as models of random surfaces, these matrix models interpolate between branched polymer behaviour and two-dimensional quantum gravity. We solve such models in a 'triple-scaling' regime in this paper, with N and M becoming large independently. A correspondence between phase transitions and singularities of mappings from R 2 to R 2 is indicated. At different critical points, the scaling behaviour is determined by (i) two decoupled ordinary differential equations; (ii) an ordinary differential equation and a finite-difference equation; or (iii) two coupled partial differential equations. The Painleve II equation arises (in conjunction with a difference equation) at a point associated with branched polymers. For critical points described by partial differential equations, there are dual weak-coupling/strong-coupling expansions. It is conjectured that the new physics is related to microscopic topology fluctuations. (orig.)

Microlocal study of S-matrix singularity structure

International Nuclear Information System (INIS)

Kawai, Takahiro; Kyoto Univ.; Stapp, H.P.

1975-01-01

Support is adduced for two related conjectures of simplicity of the analytic structure of the S-matrix and related function; namely, Sato's conjecture that the S-matrix is a solution of a maximally over-determined system of pseudo-differential equations, and our conjecture that the singularity spectrum of any bubble diagram function has the conormal structure with respect to a canonical decomposition of the solutions of the relevant Landau equations. This latter conjecture eliminates the open sets of allowed singularities that existing procedures permit. (orig.) [de

Gravitational lensing by eigenvalue distributions of random matrix models

Science.gov (United States)

Martínez Alonso, Luis; Medina, Elena

2018-05-01

We propose to use eigenvalue densities of unitary random matrix ensembles as mass distributions in gravitational lensing. The corresponding lens equations reduce to algebraic equations in the complex plane which can be treated analytically. We prove that these models can be applied to describe lensing by systems of edge-on galaxies. We illustrate our analysis with the Gaussian and the quartic unitary matrix ensembles.

Z4-symmetric factorized S-matrix in two space-time dimensions

International Nuclear Information System (INIS)

Zamolodchikov, A.B.

1979-01-01

The factorized S-matrix with internal symmetry Z 4 is constructed in two space-time dimensions. The two-particle amplitudes are obtained by means of solving the factorization, unitarity and analyticity equations. The solution of factorization equations can be expressed in terms of elliptic functions. The S-matrix cotains the resonance poles naturally. The simple formal relation between the general factorized S-matrices and the Baxter-type lattice transfer matrices is found. In the sense of this relation the Z 4 -symmetric S-matrix corresponds to the Baxter transfer matrix itself. (orig.)

SIMPLE FLUID IN AN ATTRACTIVE, DISORDERED POLYDISPERSE MATRIX

Directory of Open Access Journals (Sweden)

T.Patsahan

2004-01-01

Full Text Available The extension of the replica Ornstein-Zernike (ROZ equations is applied to the study of the structural properties of a Lennard-Jones fluid confined in an attractive polydisperse disordered matrix. The ROZ equations in combination with the orthogonal polynomial expansions for the correlation functions are used. The radial distribution functions are calculated for the adsorbed fluid at different temperatures. The effect of matrix polydispersity on the excess internal energy is considered in our study as well.

Explicit treatment of N-body correlations within a density-matrix formalism

International Nuclear Information System (INIS)

Shun-Jin, W.; Cassing, W.

1985-01-01

The nuclear many-body problem is reformulated in the density-matrix approach such that n-body correlations are separated out from the reduced density matrix rho/sub n/. A set of equations for the time evolution of the n-body correlations c/sub n/ is derived which allows for physically transparent truncations with respect to the order of correlations. In the stationary limit (c/sub n/ = 0) a restriction to two-body correlations yields a generalized Bethe-Goldstone equation a restriction to body correlations yields generalized Faddeev equations in the density-matrix formulation. Furthermore it can be shown that any truncation of the set of equations (c/sub n/ = 0, n>m) is compatible with conservation laws, a quality which in general is not fulfilled if higher order correlations are treated perturbatively

On a novel matrix method for three-dimensional photoelasticity

International Nuclear Information System (INIS)

Theocaris, P.S.; Gdoutos, E.E.

1978-01-01

A non-destructive method for the photoelastic determination of three-dimensional stress distributions, based on the Mueller and Jones calculi, is developed. The differential equations satisfied by the Stokes and Jones vectors, when a polarized light beam passes through a photoelastic model, presenting rotation of the secondary principal stress directions, are established in matrix form. The Peano-Baker method is used for the solution of these differential equations in a matrix series form, establishing the elements of the Mueller and Jones matrices of the photoelastic model. These matrices are experimentally determined by using different wavelengths in conjunction with Jones' 'equivalence theorem'. The Neumann equations are immediately deduced from the above-mentioned differential equations. (orig.) [de

Elementary matrix algebra

CERN Document Server

Hohn, Franz E

2012-01-01

This complete and coherent exposition, complemented by numerous illustrative examples, offers readers a text that can teach by itself. Fully rigorous in its treatment, it offers a mathematically sound sequencing of topics. The work starts with the most basic laws of matrix algebra and progresses to the sweep-out process for obtaining the complete solution of any given system of linear equations - homogeneous or nonhomogeneous - and the role of matrix algebra in the presentation of useful geometric ideas, techniques, and terminology.Other subjects include the complete treatment of the structur

CPDS3, Coupled 3-D Partial Differential Equation Solution

International Nuclear Information System (INIS)

Anderson, D.V.; Koniges, A.E.; Shumaker, D.E.

1992-01-01

1 - Description of program or function: CPDES3 solves the linear asymmetric matrix equations arising from coupled partial differential equations in three dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES3 allows each spatial operator to have 7, 15, 19, or 27 point stencils, permits general couplings between all of the component PDE's, and automatically generates the matrix structures needed to perform the algorithm. 2 - Method of solution: The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned bi-conjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used. 3 - Restrictions on the complexity of the problem: The discretization of the coupled three-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective

Maxwell's equations, quantum physics and the quantum graviton

International Nuclear Information System (INIS)

Gersten, Alexander; Moalem, Amnon

2011-01-01

Quantum wave equations for massless particles and arbitrary spin are derived by factorizing the d'Alembertian operator. The procedure is extensively applied to the spin one photon equation which is related to Maxwell's equations via the proportionality of the photon wavefunction Ψ to the sum E + iB of the electric and magnetic fields. Thus Maxwell's equations can be considered as the first quantized one-photon equation. The photon wave equation is written in two forms, one with additional explicit subsidiary conditions and second with the subsidiary conditions implicitly included in the main equation. The second equation was obtained by factorizing the d'Alembertian with 4×4 matrix representation of 'relativistic quaternions'. Furthermore, scalar Lagrangian formalism, consistent with quantization requirements is developed using derived conserved current of probability and normalization condition for the wavefunction. Lessons learned from the derivation of the photon equation are used in the derivation of the spin two quantum equation, which we call the quantum graviton. Quantum wave equation with implicit subsidiary conditions, which factorizes the d'Alembertian with 8×8 matrix representation of relativistic quaternions, is derived. Scalar Lagrangian is formulated and conserved probability current and wavefunction normalization are found, both consistent with the definitions of quantum operators and their expectation values. We are showing that the derived equations are the first quantized equations of the photon and the graviton.

Notes on Mayer expansions and matrix models

International Nuclear Information System (INIS)

Bourgine, Jean-Emile

2014-01-01

Mayer cluster expansion is an important tool in statistical physics to evaluate grand canonical partition functions. It has recently been applied to the Nekrasov instanton partition function of N=2 4d gauge theories. The associated canonical model involves coupled integrations that take the form of a generalized matrix model. It can be studied with the standard techniques of matrix models, in particular collective field theory and loop equations. In the first part of these notes, we explain how the results of collective field theory can be derived from the cluster expansion. The equalities between free energies at first orders is explained by the discrete Laplace transform relating canonical and grand canonical models. In a second part, we study the canonical loop equations and associate them with similar relations on the grand canonical side. It leads to relate the multi-point densities, fundamental objects of the matrix model, to the generating functions of multi-rooted clusters. Finally, a method is proposed to derive loop equations directly on the grand canonical model

Hierarchical matrix techniques for the solution of elliptic equations

KAUST Repository

Chá vez, Gustavo; Turkiyyah, George; Yokota, Rio; Keyes, David E.

2014-01-01

Hierarchical matrix approximations are a promising tool for approximating low-rank matrices given the compactness of their representation and the economy of the operations between them. Integral and differential operators have been the major

Multi-cut solutions in Chern-Simons matrix models

Science.gov (United States)

Morita, Takeshi; Sugiyama, Kento

2018-04-01

We elaborate the Chern-Simons (CS) matrix models at large N. The saddle point equations of these matrix models have a curious structure which cannot be seen in the ordinary one matrix models. Thanks to this structure, an infinite number of multi-cut solutions exist in the CS matrix models. Particularly we exactly derive the two-cut solutions at finite 't Hooft coupling in the pure CS matrix model. In the ABJM matrix model, we argue that some of multi-cut solutions might be interpreted as a condensation of the D2-brane instantons.

The relativistic electron wave equation

International Nuclear Information System (INIS)

Dirac, P.A.M.

1977-08-01

The paper was presented at the European Conference on Particle Physics held in Budapest between the 4th and 9th July of 1977. A short review is given on the birth of the relativistic electron wave equation. After Schroedinger has shown the equivalence of his wave mechanics and the matrix mechanics of Heisenberg, a general transformation theory was developed by the author. This theory required a relativistic wave equation linear in delta/delta t. As the Klein--Gordon equation available at this time did not satisfy this condition the development of a new equation became necessary. The equation which was found gave the value of the electron spin and magnetic moment automatically. (D.P.)

Nonstatic, self-consistent πN t matrix in nuclear matter

International Nuclear Information System (INIS)

Van Orden, J.W.

1984-01-01

In a recent paper, a calculation of the self-consistent πN t matrix in nuclear matter was presented. In this calculation the driving term of the self-consistent equation was chosen to be a static approximation to the free πN t matrix. In the present work, the earlier calculation is extended by using a nonstatic, fully-off-shell free πN t matrix as a starting point. Right-hand pole and cut contributions to the P-wave πN amplitudes are derived using a Low expansion and include effects due to recoil of the interacting πN system as well as the transformation from the πN c.m. frame to the nuclear rest frame. The self-consistent t-matrix equation is rewritten as two integral equations which modify the pole and cut contributions to the t matrix separately. The self-consistent πN t matrix is calculated in nuclear matter and a nonlocal optical potential is constructed from it. The resonant contribution to the optical potential is found to be broadened by 20% to 50% depending on pion momentum and is shifted upward in energy by approximately 10 MeV in comparison to the first-order optical potential. Modifications to the nucleon pole contribution are found to be negligible

Tensor-GMRES method for large sparse systems of nonlinear equations

Science.gov (United States)

Feng, Dan; Pulliam, Thomas H.

1994-01-01

This paper introduces a tensor-Krylov method, the tensor-GMRES method, for large sparse systems of nonlinear equations. This method is a coupling of tensor model formation and solution techniques for nonlinear equations with Krylov subspace projection techniques for unsymmetric systems of linear equations. Traditional tensor methods for nonlinear equations are based on a quadratic model of the nonlinear function, a standard linear model augmented by a simple second order term. These methods are shown to be significantly more efficient than standard methods both on nonsingular problems and on problems where the Jacobian matrix at the solution is singular. A major disadvantage of the traditional tensor methods is that the solution of the tensor model requires the factorization of the Jacobian matrix, which may not be suitable for problems where the Jacobian matrix is large and has a 'bad' sparsity structure for an efficient factorization. We overcome this difficulty by forming and solving the tensor model using an extension of a Newton-GMRES scheme. Like traditional tensor methods, we show that the new tensor method has significant computational advantages over the analogous Newton counterpart. Consistent with Krylov subspace based methods, the new tensor method does not depend on the factorization of the Jacobian matrix. As a matter of fact, the Jacobian matrix is never needed explicitly.

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

Explicit solutions to the generalized Sylvester matrix equation AX- XF = BY%广义Sylvester矩阵方程AX-XF=BY的显式解

Institute of Scientific and Technical Information of China (English)

周彬; 段广仁

2006-01-01

A complete, general and explicit solution to the generalized Sylvester matrix equation AX-XF = BY, with F being an arbitrary square matrix, is investigated. The proposed solution is in an extremely neat form represented by a controllability matrix of the matrix pair (A,B), a symmetric operator and an observability matrix of the matrix pair (Z,F), where Z is an arbitrary matrix used to denote the degree of freedom in the solution. Furthermore, based on the Faddeev - Leverrier algorithm, an equivalent form of the proposed solution is established. At the same time, an equivalent form of the solutions proposed in [ 13 ] is also induced. These results provide great convenience to the analysis and design problems in control systems. The results proposed in this note is a further discussion of the results proposed in [ 13 ].%给出了广义Sylvester矩阵方程AX-XF=BY当F为任意矩阵时的一种完全的解析通解.该通解由矩阵对(A,B)构成的能控性矩阵,一个对称算子矩阵和矩阵对(Z,F)构成的能观性矩阵组成,这里Z是一个任意的参数矩阵,用来表征该方程的解的自由度.利用著名的Levverrier算法,该解析解的一个等价形式被给出.给出的结果是参考文献[13]的推广,在[13]中F被假设为友矩阵.

Quasi-exact evaluation of time domain MFIE MOT matrix elements

KAUST Repository

Shi, Yifei; Bagci, Hakan; Shanker, Balasubramaniam; Lu, Mingyu; Michielssen, Eric

2013-01-01

A previously proposed quasi-exact scheme for evaluating matrix elements resulting from the marching-on-in-time (MOT) discretization of the time domain electric field integral equation (EFIE) is extended to matrix entries resulting from the discretization of its magnetic field integral equation (MFIE) counterpart. Numerical results demonstrate the accuracy of the scheme as well as the late-time stability of the resulting MOT-MFIE solver. © 2013 IEEE.

Quasi-exact evaluation of time domain MFIE MOT matrix elements

KAUST Repository

Shi, Yifei

2013-07-01

A previously proposed quasi-exact scheme for evaluating matrix elements resulting from the marching-on-in-time (MOT) discretization of the time domain electric field integral equation (EFIE) is extended to matrix entries resulting from the discretization of its magnetic field integral equation (MFIE) counterpart. Numerical results demonstrate the accuracy of the scheme as well as the late-time stability of the resulting MOT-MFIE solver. © 2013 IEEE.

The Optimization on Ranks and Inertias of a Quadratic Hermitian Matrix Function and Its Applications

Directory of Open Access Journals (Sweden)

Yirong Yao

2013-01-01

Full Text Available We solve optimization problems on the ranks and inertias of the quadratic Hermitian matrix function subject to a consistent system of matrix equations and . As applications, we derive necessary and sufficient conditions for the solvability to the systems of matrix equations and matrix inequalities , and in the Löwner partial ordering to be feasible, respectively. The findings of this paper widely extend the known results in the literature.

BKP and CKP revisited: the odd KP system

International Nuclear Information System (INIS)

Dimakis, Aristophanes; Müller-Hoissen, Folkert

2009-01-01

By restricting a linear system for the KP hierarchy to those independent variables t n with odd n, its compatibility (Zakharov-Shabat conditions) leads to the 'odd KP hierarchy'. The latter consists of pairs of equations for two dependent variables, taking values in an (typically noncommutative) associative algebra. If the algebra is commutative, the odd KP hierarchy is known to admit reductions to the BKP and the CKP hierarchy. We approach the odd KP hierarchy and its relation to BKP and CKP in different ways, and address the question of whether noncommutative versions of the BKP and the CKP equation (and some of their reductions) exist. In particular, we derive a functional representation of a linear system for the odd KP hierarchy, which in the commutative case produces functional representations of the BKP and CKP hierarchies in terms of a tau function. Furthermore, we consider a functional representation of the KP hierarchy that involves a second (auxiliary) dependent variable and features the odd KP hierarchy directly as a subhierarchy. A method to generate large classes of exact solutions to the KP hierarchy from solutions to a linear matrix ODE system, via a hierarchy of matrix Riccati equations, then also applies to the odd KP hierarchy, and this in turn can be exploited, in particular, to obtain solutions to the BKP and CKP hierarchies

MATLAB Simulation of Gradient-Based Neural Network for Online Matrix Inversion

Science.gov (United States)

Zhang, Yunong; Chen, Ke; Ma, Weimu; Li, Xiao-Dong

This paper investigates the simulation of a gradient-based recurrent neural network for online solution of the matrix-inverse problem. Several important techniques are employed as follows to simulate such a neural system. 1) Kronecker product of matrices is introduced to transform a matrix-differential-equation (MDE) to a vector-differential-equation (VDE); i.e., finally, a standard ordinary-differential-equation (ODE) is obtained. 2) MATLAB routine "ode45" is introduced to solve the transformed initial-value ODE problem. 3) In addition to various implementation errors, different kinds of activation functions are simulated to show the characteristics of such a neural network. Simulation results substantiate the theoretical analysis and efficacy of the gradient-based neural network for online constant matrix inversion.

Numerical Simulation of Antennas with Improved Integral Equation Method

International Nuclear Information System (INIS)

Ma Ji; Fang Guang-You; Lu Wei

2015-01-01

Simulating antennas around a conducting object is a challenge task in computational electromagnetism, which is concerned with the behaviour of electromagnetic fields. To analyze this model efficiently, an improved integral equation-fast Fourier transform (IE-FFT) algorithm is presented in this paper. The proposed scheme employs two Cartesian grids with different size and location to enclose the antenna and the other object, respectively. On the one hand, IE-FFT technique is used to store matrix in a sparse form and accelerate the matrix-vector multiplication for each sub-domain independently. On the other hand, the mutual interaction between sub-domains is taken as the additional exciting voltage in each matrix equation. By updating integral equations several times, the whole electromagnetic system can achieve a stable status. Finally, the validity of the presented method is verified through the analysis of typical antennas in the presence of a conducting object. (paper)

The S-matrix of superstring field theory

International Nuclear Information System (INIS)

Konopka, Sebastian

2015-01-01

We show that the classical S-matrix calculated from the recently proposed superstring field theories give the correct perturbative S-matrix. In the proof we exploit the fact that the vertices are obtained by a field redefinition in the large Hilbert space. The result extends to include the NS-NS subsector of type II superstring field theory and the recently found equations of motions for the Ramond fields. In addition, our proof implies that the S-matrix obtained from Berkovits’ WZW-like string field theory then agrees with the perturbative S-matrix to all orders.

One-dimensional integral equations for a system of three identical particles in the boundary condition models and the possibility of changing the off-shell behaviour of the two-particle t-matrix

International Nuclear Information System (INIS)

Efimov, V.N.; Schulz, H.

1976-01-01

It is shown that in the framework of the boundary condition models (BCM) for the two-particle interaction the Schroedinger equation for the system of three identical bosons can be reduced to the one-dimensional integral equation in an exact way. The method used for obtaining such an equation is based on a special consideration of the two-particle off-shell wave functions. The binding energy of the simple three-particle system is calculated. It is indicated that by means of the equation obtained it is possible to change the off-shell behaviour of the two-particle t-matrix and therefore to simulate three particle effects. (Auth.)

Quadratic Functionals with General Boundary Conditions

International Nuclear Information System (INIS)

Dosla, Z.; Dosly, O.

1997-01-01

The purpose of this paper is to give the Reid 'Roundabout Theorem' for quadratic functionals with general boundary conditions. In particular, we describe the so-called coupled point and regularity condition introduced in terms of Riccati equation solutions

Shifted Legendre method with residual error estimation for delay linear Fredholm integro-differential equations

Directory of Open Access Journals (Sweden)

Şuayip Yüzbaşı

2017-03-01

Full Text Available In this paper, we suggest a matrix method for obtaining the approximate solutions of the delay linear Fredholm integro-differential equations with constant coefficients using the shifted Legendre polynomials. The problem is considered with mixed conditions. Using the required matrix operations, the delay linear Fredholm integro-differential equation is transformed into a matrix equation. Additionally, error analysis for the method is presented using the residual function. Illustrative examples are given to demonstrate the efficiency of the method. The results obtained in this study are compared with the known results.

Analytical representation for solution of the neutron point kinetics equation with time-dependent reactivity and free of the stiffness character

International Nuclear Information System (INIS)

Silva, Milena Wollmann da

2013-01-01

In this work, we report a genuine analytical representation for the solution of the neutron point kinetics equation free of the stiffness character, assuming that the reactivity is a continuous and sectionally continuous function of time. To this end, we initially cast the point kinetics equation in a first order linear differential equation. Next, we split the corresponding matrix as a sum of a diagonal matrix with a matrix, whose components contain the off-diagonal elements. Next, expanding the neutron density and the delayed neutron precursors concentrations in a truncated series, and replacing these expansions in the matrix equation, we come out with an equation, which allows to construct a recursive system, a first order matrix differential equation with source. The fundamental characteristic of this system relies on the fact that the corresponding matrix is diagonal, meanwhile the source term is written in terms of the matrix with the off-diagonal components. Further, the first equation of the recursive system has no source and satisfies the initial conditions. On the other hand, the remaining equations satisfy the null initial condition. Due to the diagonal feature of the matrix, we attain analytical solutions for these recursive equations. We also mention that we evaluate the results for any time value, without the analytical continuity because the purposed solution is free on the stiffness character. Finally, we present numerical simulations and comparisons against literature results, considering specific the applications for the following reactivity functions: constant, step, ramp, and sine. (author)

Planck constant as spectral parameter in integrable systems and KZB equations

Science.gov (United States)

Levin, A.; Olshanetsky, M.; Zotov, A.

2014-10-01

We construct special rational gl N Knizhnik-Zamolodchikov-Bernard (KZB) equations with Ñ punctures by deformation of the corresponding quantum gl N rational R-matrix. They have two parameters. The limit of the first one brings the model to the ordinary rational KZ equation. Another one is τ. At the level of classical mechanics the deformation parameter τ allows to extend the previously obtained modified Gaudin models to the modified Schlesinger systems. Next, we notice that the identities underlying generic (elliptic) KZB equations follow from some additional relations for the properly normalized R-matrices. The relations are noncommutative analogues of identities for (scalar) elliptic functions. The simplest one is the unitarity condition. The quadratic (in R matrices) relations are generated by noncommutative Fay identities. In particular, one can derive the quantum Yang-Baxter equations from the Fay identities. The cubic relations provide identities for the KZB equations as well as quadratic relations for the classical r-matrices which can be treated as halves of the classical Yang-Baxter equation. At last we discuss the R-matrix valued linear problems which provide gl Ñ CM models and Painlevé equations via the above mentioned identities. The role of the spectral parameter plays the Planck constant of the quantum R-matrix. When the quantum gl N R-matrix is scalar ( N = 1) the linear problem reproduces the Krichever's ansatz for the Lax matrices with spectral parameter for the gl Ñ CM models. The linear problems for the quantum CM models generalize the KZ equations in the same way as the Lax pairs with spectral parameter generalize those without it.

The eigenvalues of the SN transport matrix

International Nuclear Information System (INIS)

Ourique, L.E.; Vilhena, M.T. de

2005-01-01

In a recent paper, we analyze the dependence of the eigenvalues of the S N matrix transport, associated with the system of linear differential equations that corresponds to the S N approximations of the transport equation [1]. By considering a control parameter, we have shown that there exist some bifurcation points. This means that the solutions of S N approximations change from oscillatory to non-oscillatory behavior, a different approach of the study by [2]. Nowadays, the one-dimensional transport equation and related problems have been a source of new techniques for solving particular cases as well the development of analytical methods that search aspects of existence and uniqueness of the solutions [3], [4]. In this work, we generalize the results shown in [1], searching for a model of the distribution of the bifurcation points of the S N matrix transport, studying the one-dimensional case in a slab, with anisotropic differential cross section of order 3. The result indicates that the bifurcation points obey a certain rule of distribution. Beside that, the condition number of the matrix transport increases too much in the neighborhood of these points, as we have seen in [1]. (author)

Suboptimal Regulation of a Class of Bilinear Interconnected Systems with Finite-Time Sliding Planning Horizons

Directory of Open Access Journals (Sweden)

M. de la Sen

2008-01-01

Full Text Available This paper focuses on the suboptimization of a class of multivariable discrete-time bilinear systems consisting of interconnected bilinear subsystems with respect to a linear quadratic optimal regulation criterion which involves the use of state weighting terms only. Conditions which ensure the controllability of the overall system are given as a previous requirement for optimization. Three transformations of variables are made on the system equations in order to implement the scheme on an equivalent linear system. This leads to an equivalent representation of the used quadratic performance index that involves the appearance of quadratic weighting terms related to both transformed input and state variables. In this way, a Riccati-matrix sequence, allowing the synthesis of a standard feedback control law, is obtained. Finally, the proposed control scheme is tested on realistic examples.

Matrix-exponential description of radiative transfer

International Nuclear Information System (INIS)

Waterman, P.C.

1981-01-01

By appling the matrix-exponential operator technique to the radiative-transfer equation in discrete form, new analytical solutions are obtained for the transmission and reflection matrices in the limiting cases x >1, where x is the optical depth of the layer. Orthongonality of the eigenvectors of the matrix exponential apparently yields new conditions for determining. Chandrasekhar's characteristic roots. The exact law of reflection for the discrete eigenfunctions is also obtained. Finally, when used in conjuction with the doubling method, the matrix exponential should result in reduction in both computation time and loss of precision

Nonadiabatic quantum Vlasov equation for Schwinger pair production

International Nuclear Information System (INIS)

Kim, Sang Pyo; Schubert, Christian

2011-01-01

Using Lewis-Riesenfeld theory, we derive an exact nonadiabatic master equation describing the time evolution of the QED Schwinger pair-production rate for a general time-varying electric field. This equation can be written equivalently as a first-order matrix equation, as a Vlasov-type integral equation, or as a third-order differential equation. In the last version it relates to the Korteweg-de Vries equation, which allows us to construct an exact solution using the well-known one-soliton solution to that equation. The case of timelike delta function pulse fields is also briefly considered.

CPDES2, Coupled 2-D Partial Differential Equation Solution

International Nuclear Information System (INIS)

Anderson, D.V.; Koniges, A.E.; Shumaker, D.E.

1992-01-01

1 - Description of program or function: CPDES2 solves the linear asymmetric equations arising from coupled partial differential equations in two dimensions. The exact form of the matrix depends on the choice of spatial grids and on the finite element or finite difference approximation employed. CPDES2 allows each spatial operator to have 5 or 9 point stencils, permits general coupling between all of the component PDE's, and automatically generates the matrix structures needed to perform the algorithm. 2 - Method of solution: The resulting sparse matrix equation with a complicated sub-band structure and generally asymmetric is solved by either the preconditioned conjugate gradient (CG) method or the preconditioned bi-conjugate gradient (BCG) algorithm. BCG enjoys faster convergence in most cases but in rare instances diverges. Then, CG iterations must be used. 3 - Restrictions on the complexity of the problem: The discretization of the coupled two-dimensional PDE's and their boundary conditions must result in an operator stencil which fits in the Cray2 memory. In addition, the matrix must possess a reasonable amount of diagonal dominance for the preconditioning technique to be effective

Some functional solutions of the Yang-Baxter equation

International Nuclear Information System (INIS)

Stoyanov, D.Ts.

1994-09-01

A general functional definition of the infinite dimensional quantum R-matrix satisfying the Yang-Baxter equation is given. A procedure for extracting a finite dimensional R-matrix from the general definition is demonstrated in a particular case when the group SU(2) takes place. (author). 6 refs

Efficient sparse matrix-matrix multiplication for computing periodic responses by shooting method on Intel Xeon Phi

Science.gov (United States)

Stoykov, S.; Atanassov, E.; Margenov, S.

2016-10-01

Many of the scientific applications involve sparse or dense matrix operations, such as solving linear systems, matrix-matrix products, eigensolvers, etc. In what concerns structural nonlinear dynamics, the computations of periodic responses and the determination of stability of the solution are of primary interest. Shooting method iswidely used for obtaining periodic responses of nonlinear systems. The method involves simultaneously operations with sparse and dense matrices. One of the computationally expensive operations in the method is multiplication of sparse by dense matrices. In the current work, a new algorithm for sparse matrix by dense matrix products is presented. The algorithm takes into account the structure of the sparse matrix, which is obtained by space discretization of the nonlinear Mindlin's plate equation of motion by the finite element method. The algorithm is developed to use the vector engine of Intel Xeon Phi coprocessors. It is compared with the standard sparse matrix by dense matrix algorithm and the one developed by Intel MKL and it is shown that by considering the properties of the sparse matrix better algorithms can be developed.

On a new solution of Chew-Low type equations

International Nuclear Information System (INIS)

Rerikh, K.V.

1985-01-01

The system is investigated of Chew-Low type equations, defined by the crossing symmetry matrix A(1, 1) for S-matrix elements as dfunctions of the uniformmizing variable w, in terms of which this system is a system of nonlinear difference equations. The quadratic Cremona transformation for unknown functions reducing the initial equations to a vey simple form is found. New particular solutions are obtained as functions of variable exp(cw). Existence of the new first integral γ(w) that is an even periodical function of w is estalished. The structure of the general solution depending on γ(w) and the relation of the found particular solutions with the first integral are discussed

Matrix integral solutions to the discrete KP hierarchy and its Pfaffianized version

International Nuclear Information System (INIS)

Lafortune, Stéphane; Li, Chun-Xia

2016-01-01

Matrix integrals used in random matrix theory for the study of eigenvalues of Hermitian ensembles have been shown to provide τ -functions for several hierarchies of integrable equations. In this article, we extend this relation by showing that such integrals can also provide τ -functions for the discrete KP hierarchy and a coupled version of the same hierarchy obtained through the process of Pfaffianization. To do so, we consider the first equation of the discrete KP hierarchy, the Hirota–Miwa equation. We write the Wronskian determinant solutions to the Hirota–Miwa equation and consider a particular form of matrix integrals, which we show is an example of those Wronskian solutions. The argument is then generalized to the whole hierarchy. A similar strategy is used for the Pfaffianized version of the hierarchy except that in that case, the solutions are written in terms of Pfaffians rather than determinants. (paper)

Reformulation of the Hermitean 1-matrix model as an effective field theory

Energy Technology Data Exchange (ETDEWEB)

Klitz, Alexander

2009-07-15

The formal Hermitean 1-matrix model is shown to be equivalent to an effective field theory. The correlation functions and the free energy of the matrix model correspond directly to the correlation functions and the free energy of the effective field theory. The loop equation of the field theory coupling constants is stated. Despite its length, this loop equation is simpler than the loop equations in the matrix model formalism itself since it does not contain operator inversions in any sense, but consists instead only of derivative operators and simple projection operators. Therefore the solution of the loop equation could be given for an arbitrary number of cuts up to the fifth order in the topological expansion explicitly. Two different methods of obtaining the contributions to the free energy of the higher orders are given, one depending on an operator H and one not depending on it. (orig.)

Positivity of Fundamental Matrix and Exponential Stability of Delay Differential System

Directory of Open Access Journals (Sweden)

Alexander Domoshnitsky

2014-01-01

Full Text Available The classical Wazewski theorem established that nonpositivity of all nondiagonal elements pij  (i≠j,  i,j=1,…,n is necessary and sufficient for nonnegativity of the fundamental (Cauchy matrix and consequently for applicability of the Chaplygin approach of approximate integration for system of linear ordinary differential equations xi′t+∑j=1n‍pijtxjt=fit,   i=1,…,n. Results on nonnegativity of the Cauchy matrix for system of delay differential equations xi′t+∑j=1n‍pijtxjhijt=fit,   i=1,…,n, which were based on nonpositivity of all diagonal elements, were presented in the previous works. Then examples, which demonstrated that nonpositivity of nondiagonal coefficients pij is not necessary for systems of delay equations, were found. In this paper first sufficient results about nonnegativity of the Cauchy matrix of the delay system without this assumption are proven. A necessary condition of nonnegativity of the Cauchy matrix is proposed. On the basis of these results on nonnegativity of the Cauchy matrix, necessary and sufficient conditions of the exponential stability of the delay system are obtained.

R-matrix arising from affine Hecke algebras and its application to Macdonald's difference operators

International Nuclear Information System (INIS)

Kato, Shinichi

1994-01-01

We shall give a certain trigonometric R-matrix associated with each root system by using affine Hecke algebras. From this R-matrix, we derive a quantum Knizhnik-Zamolodchikov equation after Cherednik, and show that the solutions of this KZ equation yield eigenfunctions of Macdonald's difference operators. (orig.)

Integrable boundary conditions and modified Lax equations

International Nuclear Information System (INIS)

Avan, Jean; Doikou, Anastasia

2008-01-01

We consider integrable boundary conditions for both discrete and continuum classical integrable models. Local integrals of motion generated by the corresponding 'transfer' matrices give rise to time evolution equations for the initial Lax operator. We systematically identify the modified Lax pairs for both discrete and continuum boundary integrable models, depending on the classical r-matrix and the boundary matrix

The R-matrix theory

International Nuclear Information System (INIS)

Descouvemont, P; Baye, D

2010-01-01

The different facets of the R-matrix method are presented pedagogically in a general framework. Two variants have been developed over the years: (i) The 'calculable' R-matrix method is a calculational tool to derive scattering properties from the Schroedinger equation in a large variety of physical problems. It was developed rather independently in atomic and nuclear physics with too little mutual influence. (ii) The 'phenomenological' R-matrix method is a technique to parametrize various types of cross sections. It was mainly (or uniquely) used in nuclear physics. Both directions are explained by starting from the simple problem of scattering by a potential. They are illustrated by simple examples in nuclear and atomic physics. In addition to elastic scattering, the R-matrix formalism is applied to inelastic and radiative-capture reactions. We also present more recent and more ambitious applications of the theory in nuclear physics.

Quantum Gross-Pitaevskii Equation

Directory of Open Access Journals (Sweden)

Jutho Haegeman, Damian Draxler, Vid Stojevic, J. Ignacio Cirac, Tobias J. Osborne, Frank Verstraete

2017-07-01

Full Text Available We introduce a non-commutative generalization of the Gross-Pitaevskii equation for one-dimensional quantum gasses and quantum liquids. This generalization is obtained by applying the time-dependent variational principle to the variational manifold of continuous matrix product states. This allows for a full quantum description of many body system ---including entanglement and correlations--- and thus extends significantly beyond the usual mean-field description of the Gross-Pitaevskii equation, which is known to fail for (quasi one-dimensional systems. By linearizing around a stationary solution, we furthermore derive an associated generalization of the Bogoliubov -- de Gennes equations. This framework is applied to compute the steady state response amplitude to a periodic perturbation of the potential.

On separable Pauli equations

International Nuclear Information System (INIS)

Zhalij, Alexander

2002-01-01

We classify (1+3)-dimensional Pauli equations for a spin-(1/2) particle interacting with the electro-magnetic field, that are solvable by the method of separation of variables. As a result, we obtain the 11 classes of vector-potentials of the electro-magnetic field A(t,x(vector sign))=(A 0 (t,x(vector sign)), A(vector sign)(t,x(vector sign))) providing separability of the corresponding Pauli equations. It is established, in particular, that the necessary condition for the Pauli equation to be separable into second-order matrix ordinary differential equations is its equivalence to the system of two uncoupled Schroedinger equations. In addition, the magnetic field has to be independent of spatial variables. We prove that coordinate systems and the vector-potentials of the electro-magnetic field providing the separability of the corresponding Pauli equations coincide with those for the Schroedinger equations. Furthermore, an efficient algorithm for constructing all coordinate systems providing the separability of Pauli equation with a fixed vector-potential of the electro-magnetic field is developed. Finally, we describe all vector-potentials A(t,x(vector sign)) that (a) provide the separability of Pauli equation, (b) satisfy vacuum Maxwell equations without currents, and (c) describe non-zero magnetic field

Interior Gradient Estimates for Nonuniformly Parabolic Equations II

Directory of Open Access Journals (Sweden)

Lieberman Gary M

2007-01-01

Full Text Available We prove interior gradient estimates for a large class of parabolic equations in divergence form. Using some simple ideas, we prove these estimates for several types of equations that are not amenable to previous methods. In particular, we have no restrictions on the maximum eigenvalue of the coefficient matrix and we obtain interior gradient estimates for so-called false mean curvature equation.

Equation for disentangling time-ordered exponentials with arbitrary quadratic generators

International Nuclear Information System (INIS)

Budanov, V.G.

1987-01-01

In many quantum-mechanical constructions, it is necessary to disentangle an operator-valued time-ordered exponential with time-dependent generators quadratic in the creation and annihilation operators. By disentangling, one understands the finding of the matrix elements of the time-ordered exponential or, in a more general formulation. The solution of the problem can also be reduced to calculation of a matrix time-ordered exponential that solves the corresponding classical problem. However, in either case the evolution equations in their usual form do not enable one to take into account explicitly the symmetry of the system. In this paper the methods of Weyl analysis are used to find an ordinary differential equation on a matrix Lie algebra that is invariant with respect to the adjoint action of the dynamical symmetry group of a quadratic Hamiltonian and replaces the operator evolution equation for the Green's function

Matrix analysis of electrical machinery

CERN Document Server

Hancock, N N

2013-01-01

Matrix Analysis of Electrical Machinery, Second Edition is a 14-chapter edition that covers the systematic analysis of electrical machinery performance. This edition discusses the principles of various mathematical operations and their application to electrical machinery performance calculations. The introductory chapters deal with the matrix representation of algebraic equations and their application to static electrical networks. The following chapters describe the fundamentals of different transformers and rotating machines and present torque analysis in terms of the currents based on the p

Numerical solutions of multi-order fractional differential equations by Boubaker polynomials

Directory of Open Access Journals (Sweden)

Bolandtalat A.

2016-01-01

Full Text Available In this paper, we have applied a numerical method based on Boubaker polynomials to obtain approximate numerical solutions of multi-order fractional differential equations. We obtain an operational matrix of fractional integration based on Boubaker polynomials. Using this operational matrix, the given problem is converted into a set of algebraic equations. Illustrative examples are are given to demonstrate the efficiency and simplicity of this technique.

Improved WKB radial wave functions in several bases

International Nuclear Information System (INIS)

Durand, B.; Durand, L.; Department of Physics, University of Wisconsin, Madison, Wisconsin 53706)

1986-01-01

We develop approximate WKB-like solutions to the radial Schroedinger equation for problems with an angular momentum barrier using Riccati-Bessel, Coulomb, and harmonic-oscillator functions as basis functions. The solutions treat the angular momentum singularity near the origin more accurately in leading approximation than the standard WKB solutions based on sine waves. The solutions based on Riccati-Bessel and free Coulomb wave functions continue smoothly through the inner turning point and are appropriate for scattering problems. The solutions based on oscillator and bound Coulomb wave functions incorporate both turning points smoothly and are particularly appropriate for bound-state problems; no matching of piecewise solutions using Airy functions is necessary

Multi-component bi-Hamiltonian Dirac integrable equations

Energy Technology Data Exchange (ETDEWEB)

Ma Wenxiu [Department of Mathematics and Statistics, University of South Florida, Tampa, FL 33620-5700 (United States)], E-mail: mawx@math.usf.edu

2009-01-15

A specific matrix iso-spectral problem of arbitrary order is introduced and an associated hierarchy of multi-component Dirac integrable equations is constructed within the framework of zero curvature equations. The bi-Hamiltonian structure of the obtained Dirac hierarchy is presented be means of the variational trace identity. Two examples in the cases of lower order are computed.

Yang–Mills equations on conformally connected torsion-free 4-manifolds with different signatures

Directory of Open Access Journals (Sweden)

Vyacheslav A. Luk'yanov

2017-12-01

Full Text Available In this paper we study spaces of conformal torsion-free connection of dimension 4 whose connection matrix satisfies the Yang–Mills equations. Here we generalize and strengthen the results obtained by us in previous articles, where the angular metric of these spaces had Minkowski signature. The generalization is that here we investigate the spaces of all possible metric signatures, and the enhancement is due to the fact that additional attention is paid to calculating the curvature matrix and establishing the properties of its components. It is shown that the Yang–Mills equations on 4-manifolds of conformal torsion-free connection for an arbitrary signature of the angular metric are reduced to Einstein's equations, Maxwell's equations and the equality of the Bach tensor of the angular metric and the energy-momentum tensor of the skew-symmetric charge tensor. It is proved that if the Weyl tensor is zero, the Yang–Mills equations have only self-dual or anti-self-dual solutions, i.e the curvature matrix of a conformal connection consists of self-dual or anti-self-dual external 2-forms. With the Minkowski signature (antiself-dual external 2-forms can only be zero. The components of the curvature matrix are calculated in the case when the angular metric of an arbitrary signature is Einstein, and the connection satisfies the Yang–Mills equations. In the Euclidean and pseudo-Euclidean 4-spaces we give some particular self-dual and anti-self-dual solutions of the Maxwell equations, to which all the Yang–Mills equations are reduced in this case.

A hierarchy of Liouville integrable discrete Hamiltonian equations

Energy Technology Data Exchange (ETDEWEB)

Xu Xixiang [College of Science, Shandong University of Science and Technology, Qingdao 266510 (China)], E-mail: xixiang_xu@yahoo.com.cn

2008-05-12

Based on a discrete four-by-four matrix spectral problem, a hierarchy of Lax integrable lattice equations with two potentials is derived. Two Hamiltonian forms are constructed for each lattice equation in the resulting hierarchy by means of the discrete variational identity. A strong symmetry operator of the resulting hierarchy is given. Finally, it is shown that the resulting lattice equations are all Liouville integrable discrete Hamiltonian systems.

Variational principles and Heisenberg matrix mechanics

International Nuclear Information System (INIS)

Klein, A.; Li, C.-T.

1979-01-01

If in Heisenberg's equations of motion for a problem in quantum mechanics (or quantum field theory) one studies matrix elements in the energy representation and by use of completeness conditions expresses the equations solely in terms of matrix elements of the canonical variables, and if one does likewise with the associated kinematical constraints (commutation relations), one arrives at a formulation - largely unexplored hitherto - which can be exploited for both practical and theoretical development. In this contribution, the above theme is developed within the framework of one-dimensional problems. It is shown how this formulation, both dynamics and kinematics, can be derived from a new variational principle, indeed from an entire class of such principles. A powerful method of diagonalizing the Hamiltonians by means of computations utilizing these equations is described. The variational method is shown to be particularly useful for the study of the regime of large quantum numbers. The usual WKB approximation is seen to be contained as well as a basis for the study of systematic corrections to it. Further applications in progress are mentioned. (Auth.)

Solution of fractional-order differential equations based on the operational matrices of new fractional Bernstein functions

Directory of Open Access Journals (Sweden)

M.H.T. Alshbool

2017-01-01

Full Text Available An algorithm for approximating solutions to fractional differential equations (FDEs in a modified new Bernstein polynomial basis is introduced. Writing x→xα(0<α<1 in the operational matrices of Bernstein polynomials, the fractional Bernstein polynomials are obtained and then transformed into matrix form. Furthermore, using Caputo fractional derivative, the matrix form of the fractional derivative is constructed for the fractional Bernstein matrices. We convert each term of the problem to the matrix form by means of fractional Bernstein matrices. A basic matrix equation which corresponds to a system of fractional equations is utilized, and a new system of nonlinear algebraic equations is obtained. The method is given with some priori error estimate. By using the residual correction procedure, the absolute error can be estimated. Illustrative examples are included to demonstrate the validity and applicability of the presented technique.

Inverse Interval Matrix: A Survey

Czech Academy of Sciences Publication Activity Database

Rohn, Jiří; Farhadsefat, R.

2011-01-01

Roč. 22, - (2011), s. 704-719 E-ISSN 1081-3810 R&D Projects: GA ČR GA201/09/1957; GA ČR GC201/08/J020 Institutional research plan: CEZ:AV0Z10300504 Keywords : interval matrix * inverse interval matrix * NP-hardness * enclosure * unit midpoint * inverse sign stability * nonnegative invertibility * absolute value equation * algorithm Subject RIV: BA - General Mathematics Impact factor: 0.808, year: 2010 http://www.math.technion.ac.il/iic/ ela / ela -articles/articles/vol22_pp704-719.pdf

Hecke algebraic properties of dynamical R-matrices. Application to related quantum matrix algebras

International Nuclear Information System (INIS)

Khadzhiivanov, L.K.; Todorov, I.T.; Isaev, A.P.; Pyatov, P.N.; Ogievetskij, O.V.

1998-01-01

The quantum dynamical Yang-Baxter (or Gervais-Neveu-Felder) equation defines an R-matrix R cap (p), where p stands for a set of mutually commuting variables. A family of SL (n)-type solutions of this equation provides a new realization of the Hecke algebra. We define quantum antisymmetrizers, introduce the notion of quantum determinant and compute the inverse quantum matrix for matrix algebras of the type R cap (p) a 1 a 2 = a 1 a 2 R cap. It is pointed out that such a quantum matrix algebra arises in the operator realization of the chiral zero modes of the WZNW model

Interconverting the matrix and principal-meridional representations of dioptric power and reduced vergence.

Science.gov (United States)

Harris, W F

2000-11-01

Converting the traditional representation of power as sphere, cylinder and axis to the dioptric power matrix F is usually performed by means of Long's equations and the reverse process by means of Keating's equations. It is sometimes useful to be able to convert directly between the matrix and power expressed in terms of principal powers F1 and F2 along corresponding principal meridians at angles a1 and a2. The equations for interconverting F and the principal-meridional representation expressed as F1(a1)F2 are presented here. Equivalent equations allow direct interconversion of the reduced vergence matrix L and the principal-meridional representation of vergence L1(a1)L2. Vergence becomes infinite at line and point focuses. Similarly effective power and back- and front-vertex power are infinite for some systems. Nevertheless it is possible unambiguously to represent infinite vergence and vertex power in principal-meridional form. However, information is usually lost in these infinite cases when the principal-meridional representation is converted to the matrix representation, and the former is not recoverable from the latter. As a consequence the matrix representation is usually unsatisfactory for vergences and vertex powers that are infinite. On the other hand, the principal-meridional representation of vergence and power is always satisfactory. If one adopts the position that effective powers and vertex powers are really vergences rather than powers then one concludes that the matrix provides a satisfactory representation for powers of thin systems in general but not for vergences. Implied by a vergence at a point is an interval of Sturm. The equations for characterizing the interval from the reduced vergence are presented.

Linear orbit parameters for the exact equations of motion

International Nuclear Information System (INIS)

Parzen, G.

1995-01-01

This paper defines the beta function and other linear orbit parameters using the exact equations of motion. The β, α and ψ functions are redefined using the exact equations. Expressions are found for the transfer matrix and the emittance. The differential equations for η = x/β 1/2 is found. New relationships between α, β, ψ and ν are derived

Isotropic stars in general relativity

International Nuclear Information System (INIS)

Mak, M.K.; Harko, T.

2013-01-01

We present a general solution of the Einstein gravitational field equations for the static spherically symmetric gravitational interior space-time of an isotropic fluid sphere. The solution is obtained by transforming the pressure isotropy condition, a second order ordinary differential equation, into a Riccati type first order differential equation, and using a general integrability condition for the Riccati equation. This allows us to obtain an exact non-singular solution of the interior field equations for a fluid sphere, expressed in the form of infinite power series. The physical features of the solution are studied in detail numerically by cutting the infinite series expansions, and restricting our numerical analysis by taking into account only n=21 terms in the power series representations of the relevant astrophysical parameters. In the present model all physical quantities (density, pressure, speed of sound etc.) are finite at the center of the sphere. The physical behavior of the solution essentially depends on the equation of state of the dense matter at the center of the star. The stability properties of the model are also analyzed in detail for a number of central equations of state, and it is shown that it is stable with respect to the radial adiabatic perturbations. The astrophysical analysis indicates that this solution can be used as a realistic model for static general relativistic high density objects, like neutron stars. (orig.)

Series: Utilization of Differential Equations and Methods for Solving Them in Medical Physics (3).

Science.gov (United States)

Murase, Kenya

2016-01-01

In this issue, simultaneous differential equations were introduced. These differential equations are often used in the field of medical physics. The methods for solving them were also introduced, which include Laplace transform and matrix methods. Some examples were also introduced, in which Laplace transform and matrix methods were applied to solving simultaneous differential equations derived from a three-compartment kinetic model for analyzing the glucose metabolism in tissues and Bloch equations for describing the behavior of the macroscopic magnetization in magnetic resonance imaging.In the next (final) issue, partial differential equations and various methods for solving them will be introduced together with some examples in medical physics.

Inverse scattering scheme for the Dirac equation at fixed energy

International Nuclear Information System (INIS)

Leeb, H.; Lehninger, H.; Schilder, C.

2001-01-01

Full text: Based on the concept of generalized transformation operators a new hierarchy of Dirac equations with spherical symmetric scalar and fourth component vector potentials is presented. Within this hierarchy closed form expressions for the solutions, the potentials and the S-matrix can be given in terms of solutions of the original Dirac equation. Using these transformations an inverse scattering scheme has been constructed for the Dirac equation which is the analog to the rational scheme in the non-relativistic case. The given method provides for the first time an inversion scheme with closed form expressions for the S-matrix for non-relativistic scattering problems with central and spin-orbit potentials. (author)

Mean-variance portfolio selection for defined-contribution pension funds with stochastic salary.

Science.gov (United States)

Zhang, Chubing

2014-01-01

This paper focuses on a continuous-time dynamic mean-variance portfolio selection problem of defined-contribution pension funds with stochastic salary, whose risk comes from both financial market and nonfinancial market. By constructing a special Riccati equation as a continuous (actually a viscosity) solution to the HJB equation, we obtain an explicit closed form solution for the optimal investment portfolio as well as the efficient frontier.

Mean-Variance Portfolio Selection for Defined-Contribution Pension Funds with Stochastic Salary

OpenAIRE

Chubing Zhang

2014-01-01

This paper focuses on a continuous-time dynamic mean-variance portfolio selection problem of defined-contribution pension funds with stochastic salary, whose risk comes from both financial market and nonfinancial market. By constructing a special Riccati equation as a continuous (actually a viscosity) solution to the HJB equation, we obtain an explicit closed form solution for the optimal investment portfolio as well as the efficient frontier.

Continuous affine processes

DEFF Research Database (Denmark)

Buchardt, Kristian

2016-01-01

Affine processes possess the property that expectations of exponential affine transformations are given by a set of Riccati differential equations, which is the main feature of this popular class of processes. In this paper we generalise these results for expectations of more general transformati...

Block circulant and block Toeplitz approximants of a class of spatially distributed systems-An LQR perspective

NARCIS (Netherlands)

Iftime, Orest V.

2012-01-01

In this paper block circulant and block Toeplitz long strings of MIMO systems with finite length are compared with their corresponding infinite-dimensional spatially invariant systems. The focus is on the convergence of the sequence of solutions to the control Riccati equations and the convergence

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

NARCIS (Netherlands)

Opmeer, MR; Curtain, RF

2004-01-01

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

Two derivations of the master equation of quantum Brownian motion

International Nuclear Information System (INIS)

Halliwell, J J

2007-01-01

Central to many discussion of decoherence is a master equation for the reduced density matrix of a massive particle experiencing scattering from its surrounding environment, such as that of Joos and Zeh. Such master equations enjoy a close relationship with spontaneous localization models, like the GRW model. The aim of this paper is to present two derivations of the master equation. The first derivation is a pedagogical model designed to illustrate the origins of the master equation as simply as possible, focusing on physical principles and without the complications of S-matrix theory. This derivation may serve as a useful tutorial example for students attempting to learn this subject area. The second is the opposite: a very general derivation using non-relativistic many-body field theory. It reduces to the equation of the type given by Joos and Zeh in the one-particle sector, but correcting certain numerical factors which have recently become significant in connection with experimental tests of decoherence. This master equation also emphasizes the role of local number density as the 'preferred basis' for decoherence in this model

Two derivations of the master equation of quantum Brownian motion

Energy Technology Data Exchange (ETDEWEB)

Halliwell, J J [Blackett Laboratory, Imperial College, London SW7 2BZ (United Kingdom)

2007-03-23

Central to many discussion of decoherence is a master equation for the reduced density matrix of a massive particle experiencing scattering from its surrounding environment, such as that of Joos and Zeh. Such master equations enjoy a close relationship with spontaneous localization models, like the GRW model. The aim of this paper is to present two derivations of the master equation. The first derivation is a pedagogical model designed to illustrate the origins of the master equation as simply as possible, focusing on physical principles and without the complications of S-matrix theory. This derivation may serve as a useful tutorial example for students attempting to learn this subject area. The second is the opposite: a very general derivation using non-relativistic many-body field theory. It reduces to the equation of the type given by Joos and Zeh in the one-particle sector, but correcting certain numerical factors which have recently become significant in connection with experimental tests of decoherence. This master equation also emphasizes the role of local number density as the 'preferred basis' for decoherence in this model.

The universal R-matrix and its associated quantum algebra as functionals of the classical r-matrix: the sl2 case

International Nuclear Information System (INIS)

Freidel, L.; Maillet, J.M.

1992-09-01

Using a geometrical approach to the quantum Yang-Baxter equation, the quantum algebra U h (sl 2 ) and its universal quantum R-matrix are explicitly constructed as functionals of the associated classical r-matrix. In this framework, the quantum algebra U h (sl 2 ) is naturally imbedded in the universal enveloping algebra of the sl 2 current algebra. (author) 13 refs

Frozen and broken color: a matrix Schroedinger equation in the semiclassical limit

International Nuclear Information System (INIS)

Orbach, H.S.

1981-01-01

We consider the case of frozen color, i.e, where global color symmetry remains exact, but where colored states have a mass large compared to color-singlet mesons. Using semiclassical WKB formalism, we construct the spectrum of bound states. In order to determine the charge of the constituents, we then consider deep-inelastic scattering of an external probe (e.g., lepton) from our one-dimensional meson. We calculate explicitly the structure function, W, in the WKB limit and show how Lipkin's mechanism is manifested, as well as how scaling behavior comes. We derive the WKB formalism as a special case of a method of obtaining WKB type solutions for generalized Schroedinger equations for which the Hamiltonian is an arbitrary matrix function of any number of pairs of canonical operators. We generalize these considerations to the case of broken color symmetry - but where the breaking is not so strong as to allow low-lying states to have a large amount of mixing with the colored states. In this case, the degeneracy of excited colored states can be broken. We find that local excitation of color guarantees global excitation of color; i.e., if at a given energy colored semiclassical states can be constructed with size comparable to that of the ground state wave function, colored states of that energy will also exist in the spectrum of the full theory and will be observed. However, global existence of color does not guarantee the excitation of colored states via deep-inelastic processes

Calculating Relativistic Transition Matrix Elements for Hydrogenic Atoms Using Monte Carlo Methods

Science.gov (United States)

Alexander, Steven; Coldwell, R. L.

2015-03-01

The nonrelativistic transition matrix elements for hydrogen atoms can be computed exactly and these expressions are given in a number of classic textbooks. The relativistic counterparts of these equations can also be computed exactly but these expressions have been described in only a few places in the literature. In part, this is because the relativistic equations lack the elegant simplicity of the nonrelativistic equations. In this poster I will describe how variational Monte Carlo methods can be used to calculate the energy and properties of relativistic hydrogen atoms and how the wavefunctions for these systems can be used to calculate transition matrix elements.

Direct integration of the S-matrix applied to rigorous diffraction

International Nuclear Information System (INIS)

Iff, W; Lindlein, N; Tishchenko, A V

2014-01-01

A novel Fourier method for rigorous diffraction computation at periodic structures is presented. The procedure is based on a differential equation for the S-matrix, which allows direct integration of the S-matrix blocks. This results in a new method in Fourier space, which can be considered as a numerically stable and well-parallelizable alternative to the conventional differential method based on T-matrix integration and subsequent conversions from the T-matrices to S-matrix blocks. Integration of the novel differential equation in implicit manner is expounded. The applicability of the new method is shown on the basis of 1D periodic structures. It is clear however, that the new technique can also be applied to arbitrary 2D periodic or periodized structures. The complexity of the new method is O(N 3 ) similar to the conventional differential method with N being the number of diffraction orders. (fast track communication)

Structure-preserving algorithms for the Duffing equation

International Nuclear Information System (INIS)

Gang Tieqiang; Mei Fengxiang; Xie Jiafang

2008-01-01

In this paper, the dissipative and the forced terms of the Duffing equation are considered as the perturbations of nonlinear Hamiltonian equations and the perturbational effect is indicated by parameter ε. Firstly, based on the gradient-Hamiltonian decomposition theory of vector fields, by using splitting methods, this paper constructs structure-preserving algorithms (SPAs) for the Duffing equation. Then, according to the Liouville formula, it proves that the Jacobian matrix determinants of the SPAs are equal to that of the exact flow of the Duffing equation. However, considering the explicit Runge–Kutta methods, this paper finds that there is an error term of order p+1 for the Jacobian matrix determinants. The volume evolution law of a given region in phase space is discussed for different algorithms, respectively. As a result, the sum of Lyapunov exponents is exactly invariable for the SPAs proposed in this paper. Finally, through numerical experiments, relative norm errors and absolute energy errors of phase trajectories of the SPAs and the Heun method (a second-order Runge–Kutta method) are compared. Computational results illustrate that the SPAs are evidently better than the Heun method when ε is small or equal to zero. (general)

Optimal Sliding Mode Controllers for Attitude Stabilization of Flexible Spacecraft

Directory of Open Access Journals (Sweden)

Chutiphon Pukdeboon

2011-01-01

Full Text Available The robust optimal attitude control problem for a flexible spacecraft is considered. Two optimal sliding mode control laws that ensure the exponential convergence of the attitude control system are developed. Integral sliding mode control (ISMC is applied to combine the first-order sliding mode with optimal control and is used to control quaternion-based spacecraft attitude manoeuvres with external disturbances and an uncertainty inertia matrix. For the optimal control part the state-dependent Riccati equation (SDRE and optimal Lyapunov techniques are employed to solve the infinite-time nonlinear optimal control problem. The second method of Lyapunov is used to guarantee the stability of the attitude control system under the action of the proposed control laws. An example of multiaxial attitude manoeuvres is presented and simulation results are included to verify the usefulness of the developed controllers.

Computing with linear equations and matrices

International Nuclear Information System (INIS)

Churchhouse, R.F.

1983-01-01

Systems of linear equations and matrices arise in many disciplines. The equations may accurately represent conditions satisfied by a system or, more likely, provide an approximation to a more complex system of non-linear or differential equations. The system may involve a few or many thousand unknowns and each individual equation may involve few or many of them. Over the past 50 years a vast literature on methods for solving systems of linear equations and the associated problems of finding the inverse or eigenvalues of a matrix has been produced. These lectures cover those methods which have been found to be most useful for dealing with such types of problem. References are given where appropriate and attention is drawn to the possibility of improved methods for use on vector and parallel processors. (orig.)

Integrable peakon equations with cubic nonlinearity

International Nuclear Information System (INIS)

Hone, Andrew N W; Wang, J P

2008-01-01

We present a new integrable partial differential equation found by Vladimir Novikov. Like the Camassa-Holm and Degasperis-Procesi equations, this new equation admits peaked soliton (peakon) solutions, but it has nonlinear terms that are cubic, rather than quadratic. We give a matrix Lax pair for V Novikov's equation, and show how it is related by a reciprocal transformation to a negative flow in the Sawada-Kotera hierarchy. Infinitely many conserved quantities are found, as well as a bi-Hamiltonian structure. The latter is used to obtain the Hamiltonian form of the finite-dimensional system for the interaction of N peakons, and the two-body dynamics (N = 2) is explicitly integrated. Finally, all of this is compared with some analogous results for another cubic peakon equation derived by Zhijun Qiao. (fast track communication)

Discrete-ordinate method with matrix exponential for a pseudo-spherical atmosphere: Vector case

International Nuclear Information System (INIS)

Doicu, A.; Trautmann, T.

2009-01-01

The paper is devoted to the extension of the matrix-exponential formalism for the scalar radiative transfer to the vector case. Using basic results of the theory of matrix-exponential functions we provide a compact and versatile formulation of the vector radiative transfer. As in the scalar case, we operate with the concept of the layer equation incorporating the level values of the Stokes vector. The matrix exponentials which enter in the expression of the layer equation are computed by using the matrix eigenvalue method and the Pade approximation. A discussion of the computational efficiency of the proposed method for both an aerosol-loaded atmosphere as well as a cloudy atmosphere is also provided

On renormalization group flow in matrix model

International Nuclear Information System (INIS)

Gao, H.B.

1992-10-01

The renormalization group flow recently found by Brezin and Zinn-Justin by integrating out redundant entries of the (N+1)x(N+1) Hermitian random matrix is studied. By introducing explicitly the RG flow parameter, and adding suitable counter terms to the matrix potential of the one matrix model, we deduce some interesting properties of the RG trajectories. In particular, the string equation for the general massive model interpolating between the UV and IR fixed points turns out to be a consequence of RG flow. An ambiguity in the UV region of the RG trajectory is remarked to be related to the large order behaviour of the one matrix model. (author). 7 refs

Matrix-oriented implementation for the numerical solution of the partial differential equations governing flows and transport in porous media

KAUST Repository

Sun, Shuyu

2012-09-01

In this paper we introduce a new technique for the numerical solution of the various partial differential equations governing flow and transport phenomena in porous media. This method is proposed to be used in high level programming languages like MATLAB, Python, etc., which show to be more efficient for certain mathematical operations than for others. The proposed technique utilizes those operations in which these programming languages are efficient the most and keeps away as much as possible from those inefficient, time-consuming operations. In particular, this technique is based on the minimization of using multiple indices looping operations by reshaping the unknown variables into one-dimensional column vectors and performing the numerical operations using shifting matrices. The cell-centered information as well as the face-centered information are shifted to the adjacent face-center and cell-center, respectively. This enables the difference equations to be done for all the cells at once using matrix operations rather than within loops. Furthermore, for results post-processing, the face-center information can further be mapped to the physical grid nodes for contour plotting and stream lines constructions. In this work we apply this technique to flow and transport phenomena in porous media. © 2012 Elsevier Ltd.

The coordinate system transformation of a serial kinematic structures and use in the derivation of systems motion equations

Directory of Open Access Journals (Sweden)

Zátopek Jiří

2016-01-01

Full Text Available This text discusses the use of transformation matrices to determine the motion equations of the complex mechanical structure. Use of the transformation matrix does not apply only to motion equations but has the general use in relative positions determine of objects in the 3D space. Analysed model is divided into seven physical objects, the transformation matrix and the corresponding inertia/pseudo-inertia matrix is included in each of them. This matrices are strictly necessary to the system dynamic description using the matrix form of Lagrange Equations of the Second Type. Another possibility to use the transformation matrix is shown in the camera system measurement. Model was designed in 3D CAD system SolidWorks, MATLAB was used for the mathematical calculations.

Computer programs for the solution of systems of linear algebraic equations

Science.gov (United States)

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.

A Numerical Algorithm to find All Scalar Feedback Nash Equilibria

NARCIS (Netherlands)

Engwerda, J.C.

2013-01-01

Abstract: In this note we generalize a numerical algorithm presented in [9] to calculate all solutions of the scalar algebraic Riccati equations that play an important role in finding feedback Nash equilibria of the scalar N-player linear affine-quadratic differential game. The algorithm is based on

Response matrix properties and convergence implications for an interface-current nodal formulation

International Nuclear Information System (INIS)

Yang, W.S.

1995-01-01

An analytic study was performed of the properties and the associated convergence implications of the response matrix equations derived via the widely used nodal expansion method. By using the DIF3D nodal formulation in hexagonal-z geometry as a concrete example, an analytic expression for the response matrix is first derived by using the hexagonal prism symmetry transformations. The spectral radius of the local response matrix is shown to be always 2 -norm of the response matrix is shown to be ∞ -norm is not always 2 - and l ∞ -norms of the response matrix are found to increase as the removal cross section decreases. On the other hand, for a given removal cross section, each of these matrix norms takes its minimum at a certain diffusion coefficient and increases as the diffusion coefficient deviates from this value. Based on these matrix norms, sufficient conditions for the convergence of the iteration schemes for solving the response matrix equations are discussed. The range of node-height-to-hexagon-pitch ratios that guarantees a positive solution is derived as a function of the diffusion coefficient and the removal cross section

Matrix inversion tomosynthesis improvements in longitudinal x-ray slice imaging

International Nuclear Information System (INIS)

Dobbines, J.T. III.

1990-01-01

This patent describes a tomosynthesis apparatus. It comprises: an x-ray tomography machine for producing a plurality of x-ray projection images of a subject including an x-ray source, and detection means; and processing means, connected to receive the plurality of projection images, for: shifting and reconstructing the projection x-ray images to obtain a tomosynthesis matrix of images T; acquiring a blurring matrix F having components which represent out-of-focus and in-focus components of the matrix T; obtaining a matrix P representing only in-focus components of the imaged subject by solving a matrix equation including the matrix T and the matrix F; correcting the matrix P for low spatial frequency components; and displaying images indicative of contents of the matrix P

Solution of the Lyapunov matrix equation for a system with a time-dependent stiffness matrix

DEFF Research Database (Denmark)

Pommer, Christian; Kliem, Wolfhard

2004-01-01

The stability of the linearized model of a rotor system with non-symmetric strain and axial loads is investigated. Since we are using a fixed reference system, the differential equations have the advantage to be free of Coriolis and centrifugal forces. A disadvantage is nevertheless the occurrence...

Nonlinear integral equations for the sausage model

Science.gov (United States)

Ahn, Changrim; Balog, Janos; Ravanini, Francesco

2017-08-01

The sausage model, first proposed by Fateev, Onofri, and Zamolodchikov, is a deformation of the O(3) sigma model preserving integrability. The target space is deformed from the sphere to ‘sausage’ shape by a deformation parameter ν. This model is defined by a factorizable S-matrix which is obtained by deforming that of the O(3) sigma model by a parameter λ. Clues for the deformed sigma model are provided by various UV and IR information through the thermodynamic Bethe ansatz (TBA) analysis based on the S-matrix. Application of TBA to the sausage model is, however, limited to the case of 1/λ integer where the coupled integral equations can be truncated to a finite number. In this paper, we propose a finite set of nonlinear integral equations (NLIEs), which are applicable to generic value of λ. Our derivation is based on T-Q relations extracted from the truncated TBA equations. For a consistency check, we compute next-leading order corrections of the vacuum energy and extract the S-matrix information in the IR limit. We also solved the NLIE both analytically and numerically in the UV limit to get the effective central charge and compared with that of the zero-mode dynamics to obtain exact relation between ν and λ. Dedicated to the memory of Petr Petrovich Kulish.

Improved Asymmetric Cipher Based on Matrix Power Function with Provable Security

Directory of Open Access Journals (Sweden)

Eligijus Sakalauskas

2017-01-01

Full Text Available The improved version of the author’s previously declared asymmetric cipher protocol based on matrix power function (MPF is presented. Proposed modification avoids discrete logarithm attack (DLA which could be applied to the previously declared protocol. This attack allows us to transform the initial system of MPF equations to so-called matrix multivariate quadratic (MMQ system of equations, which is a system representing a subclass of multivariate quadratic (MQ systems of equations. We are making a conjecture that avoidance of DLA in protocol, presented here, should increase its security, since an attempt to solve the initial system of MPF equations would appear to be no less complex than solving the system of MMQ equations. No algorithms are known to solve such a system of equations. Security parameters and their secure values are defined. Security analysis against chosen plaintext attack (CPA and chosen ciphertext attack (CCA is presented. Measures taken to prevent DLA attack increase the security of this protocol with respect to the previously declated protocol.

Raney Distributions and Random Matrix Theory

Science.gov (United States)

Forrester, Peter J.; Liu, Dang-Zheng

2015-03-01

Recent works have shown that the family of probability distributions with moments given by the Fuss-Catalan numbers permit a simple parameterized form for their density. We extend this result to the Raney distribution which by definition has its moments given by a generalization of the Fuss-Catalan numbers. Such computations begin with an algebraic equation satisfied by the Stieltjes transform, which we show can be derived from the linear differential equation satisfied by the characteristic polynomial of random matrix realizations of the Raney distribution. For the Fuss-Catalan distribution, an equilibrium problem characterizing the density is identified. The Stieltjes transform for the limiting spectral density of the singular values squared of the matrix product formed from inverse standard Gaussian matrices, and standard Gaussian matrices, is shown to satisfy a variant of the algebraic equation relating to the Raney distribution. Supported on , we show that it too permits a simple functional form upon the introduction of an appropriate choice of parameterization. As an application, the leading asymptotic form of the density as the endpoints of the support are approached is computed, and is shown to have some universal features.

Reduction of Under-Determined Linear Systems by Sparce Block Matrix Technique

DEFF Research Database (Denmark)

Tarp-Johansen, Niels Jacob; Poulsen, Peter Noe; Damkilde, Lars

1996-01-01

numerical stability of the aforementioned reduction. Moreover the coefficient matrix for the equilibrium equations is typically very sparse. The objective is to deal efficiently with the full pivoting reduction of sparse rectangular matrices using a dynamic storage scheme based on the block matrix concept.......Under-determined linear equation systems occur in different engineering applications. In structural engineering they typically appear when applying the force method. As an example one could mention limit load analysis based on The Lower Bound Theorem. In this application there is a set of under......-determined equilibrium equation restrictions in an LP-problem. A significant reduction of computer time spent on solving the LP-problem is achieved if the equilib rium equations are reduced before going into the optimization procedure. Experience has shown that for some structures one must apply full pivoting to ensure...

Mean-Variance Portfolio Selection for Defined-Contribution Pension Funds with Stochastic Salary

Directory of Open Access Journals (Sweden)

Chubing Zhang

2014-01-01

Full Text Available This paper focuses on a continuous-time dynamic mean-variance portfolio selection problem of defined-contribution pension funds with stochastic salary, whose risk comes from both financial market and nonfinancial market. By constructing a special Riccati equation as a continuous (actually a viscosity solution to the HJB equation, we obtain an explicit closed form solution for the optimal investment portfolio as well as the efficient frontier.

Iterative Methods for Solving Nonlinear Parabolic Problem in Pension Saving Management

Science.gov (United States)

Koleva, M. N.

2011-11-01

In this work we consider a nonlinear parabolic equation, obtained from Riccati like transformation of the Hamilton-Jacobi-Bellman equation, arising in pension saving management. We discuss two numerical iterative methods for solving the model problem—fully implicit Picard method and mixed Picard-Newton method, which preserves the parabolic characteristics of the differential problem. Numerical experiments for comparison the accuracy and effectiveness of the algorithms are discussed. Finally, observations are given.

Mean-Variance Portfolio Selection for Defined-Contribution Pension Funds with Stochastic Salary

Science.gov (United States)

Zhang, Chubing

2014-01-01

This paper focuses on a continuous-time dynamic mean-variance portfolio selection problem of defined-contribution pension funds with stochastic salary, whose risk comes from both financial market and nonfinancial market. By constructing a special Riccati equation as a continuous (actually a viscosity) solution to the HJB equation, we obtain an explicit closed form solution for the optimal investment portfolio as well as the efficient frontier. PMID:24782667

Simultaneous Balancing and Model Reduction of Switched Linear Systems

OpenAIRE

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

2011-01-01

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

Hecke symmetries and characteristic relations on reflection equation algebras

International Nuclear Information System (INIS)

Gurevich, D.I.; Pyatov, P.N.

1996-01-01

We discuss how properties of Hecke symmetry (i.e., Hecke type R-matrix) influence the algebraic structure of the corresponding Reflection Equation (RE) algebra. Analogues of the Newton relations and Cayley-Hamilton theorem for the matrix of generators of the RE algebra related to a finite rank even Hecke symmetry are derived. 10 refs

Comparison of the iterated equation of motion approach and the density matrix formalism for the quantum Rabi model

Science.gov (United States)

Kalthoff, Mona; Keim, Frederik; Krull, Holger; Uhrig, Götz S.

2017-05-01

The density matrix formalism and the equation of motion approach are two semi-analytical methods that can be used to compute the non-equilibrium dynamics of correlated systems. While for a bilinear Hamiltonian both formalisms yield the exact result, for any non-bilinear Hamiltonian a truncation is necessary. Due to the fact that the commonly used truncation schemes differ for these two methods, the accuracy of the obtained results depends significantly on the chosen approach. In this paper, both formalisms are applied to the quantum Rabi model. This allows us to compare the approximate results and the exact dynamics of the system and enables us to discuss the accuracy of the approximations as well as the advantages and the disadvantages of both methods. It is shown to which extent the results fulfill physical requirements for the observables and which properties of the methods lead to unphysical results.

Dynamical TAP equations for non-equilibrium Ising spin glasses

DEFF Research Database (Denmark)

Roudi, Yasser; Hertz, John

2011-01-01

We derive and study dynamical TAP equations for Ising spin glasses obeying both synchronous and asynchronous dynamics using a generating functional approach. The system can have an asymmetric coupling matrix, and the external fields can be time-dependent. In the synchronously updated model, the TAP...... equations take the form of self consistent equations for magnetizations at time t+1, given the magnetizations at time t. In the asynchronously updated model, the TAP equations determine the time derivatives of the magnetizations at each time, again via self consistent equations, given the current values...... of the magnetizations. Numerical simulations suggest that the TAP equations become exact for large systems....

On the universal R-matrix for the Izergin-Korepin model

International Nuclear Information System (INIS)

Boos, Herman; Goehmann, Frank; Kluemper, Andreas; Nirov, Khazret S; Razumov, Alexander V

2011-01-01

We continue our exercises with the universal R-matrix based on the Khoroshkin and Tolstoy formula. Here we present our results for the case of the twisted affine Kac-Moody Lie algebra of type A (2) 2 . Our interest in this case is inspired by the fact that the Tzitzeica equation is associated with A (2) 2 in a similar way as the sine-Gordon equation is related to A (1) 1 . The fundamental spin-chain Hamiltonian is constructed systematically as the logarithmic derivative of the transfer matrix. L-operators of two types are obtained by using q-deformed oscillators. (paper)

An iterative method to invert the LTSn matrix

Energy Technology Data Exchange (ETDEWEB)

Cardona, A.V.; Vilhena, M.T. de [UFRGS, Porto Alegre (Brazil)

1996-12-31

Recently Vilhena and Barichello proposed the LTSn method to solve, analytically, the Discrete Ordinates Problem (Sn problem) in transport theory. The main feature of this method consist in the application of the Laplace transform to the set of Sn equations and solve the resulting algebraic system for the transport flux. Barichello solve the linear system containing the parameter s applying the definition of matrix invertion exploiting the structure of the LTSn matrix. In this work, it is proposed a new scheme to invert the LTSn matrix, decomposing it in blocks and recursively inverting this blocks.

Self-similar solutions of the modified nonlinear schrodinger equation

International Nuclear Information System (INIS)

Kitaev, A.V.

1986-01-01

This paper considers a 2 x 2 matrix linear ordinary differential equation with large parameter t and irregular singular point of fourth order at infinity. The leading order of the monodromy data of this equation is calculated in terms of its coefficients. Isomonodromic deformations of the equation are self-similar solutions of the modified nonlinear Schrodinger equation, and therefore inversion of the expressions obtained for the monodromy data gives the leading term in the time-asymptotic behavior of the self-similar solution. The application of these results to the type IV Painleve equation is considered in detail

Chern-Simons matrix models and unoriented strings

International Nuclear Information System (INIS)

Halmagyi, Nick; Yasnov, Vadim

2004-01-01

For matrix models with measure on the Lie algebra of SO/Sp, the sub-leading free energy is given by F 1 (S) ±{1/4}({δF 0 (S)}/{δS}). Motivated by the fact that this relationship does not hold for Chern-Simons theory on S 3 , we calculate the sub-leading free energy in the matrix model for this theory, which is a Gaussian matrix model with Haar measure on the group SO/Sp. We derive a quantum loop equation for this matrix model and then find that F 1 is an integral of the leading order resolvent over the spectral curve. We explicitly calculate this integral for quadratic potential and find agreement with previous studies of SO/Sp Chern-Simons theory. (author)

On the Curvature and Heat Flow on Hamiltonian Systems

Directory of Open Access Journals (Sweden)

Ohta Shin-ichi

2014-01-01

Full Text Available We develop the differential geometric and geometric analytic studies of Hamiltonian systems. Key ingredients are the curvature operator, the weighted Laplacian, and the associated Riccati equation.We prove appropriate generalizations of the Bochner-Weitzenböck formula and Laplacian comparison theorem, and study the heat flow.

Perturbation theory for continuous stochastic equations

International Nuclear Information System (INIS)

Chechetkin, V.R.; Lutovinov, V.S.

1987-01-01

The various general perturbational schemes for continuous stochastic equations are considered. These schemes have many analogous features with the iterational solution of Schwinger equation for S-matrix. The following problems are discussed: continuous stochastic evolution equations for probability distribution functionals, evolution equations for equal time correlators, perturbation theory for Gaussian and Poissonian additive noise, perturbation theory for birth and death processes, stochastic properties of systems with multiplicative noise. The general results are illustrated by diffusion-controlled reactions, fluctuations in closed systems with chemical processes, propagation of waves in random media in parabolic equation approximation, and non-equilibrium phase transitions in systems with Poissonian breeding centers. The rate of irreversible reaction X + X → A (Smoluchowski process) is calculated with the use of general theory based on continuous stochastic equations for birth and death processes. The threshold criterion and range of fluctuational region for synergetic phase transition in system with Poissonian breeding centers are also considered. (author)

On the exact solution for the multi-group kinetic neutron diffusion equation in a rectangle

International Nuclear Information System (INIS)

Petersen, C.Z.; Vilhena, M.T.M.B. de; Bodmann, B.E.J.

2011-01-01

In this work we consider the two-group bi-dimensional kinetic neutron diffusion equation. The solution procedure formalism is general with respect to the number of energy groups, neutron precursor families and regions with different chemical compositions. The fast and thermal flux and the delayed neutron precursor yields are expanded in a truncated double series in terms of eigenfunctions that, upon insertion into the kinetic equation and upon taking moments, results in a first order linear differential matrix equation with source terms. We split the matrix appearing in the transformed problem into a sum of a diagonal matrix plus the matrix containing the remaining terms and recast the transformed problem into a form that can be solved in the spirit of Adomian's recursive decomposition formalism. Convergence of the solution is guaranteed by the Cardinal Interpolation Theorem. We give numerical simulations and comparisons with available results in the literature. (author)

The General Analytic Solution of a Functional Equation of Addition Type

OpenAIRE

Braden, H. W.; Buchstaber, V. M.

1995-01-01

The general analytic solution to the functional equation $$ \\phi_1(x+y)= { { \\biggl|\\matrix{\\phi_2(x)&\\phi_2(y)\\cr\\phi_3(x)&\\phi_3(y)\\cr}\\biggr|} \\over { \\biggl|\\matrix{\\phi_4(x)&\\phi_4(y)\\cr\\phi_5(x)&\\phi_5(y)\\cr}\\biggr|} } $$ is characterised. Up to the action of the symmetry group, this is described in terms of Weierstrass elliptic functions. We illustrate our theory by applying it to the classical addition theorems of the Jacobi elliptic functions and the functional equations $$ \\phi_1(x+...

Acoustic 3D modeling by the method of integral equations

Science.gov (United States)

Malovichko, M.; Khokhlov, N.; Yavich, N.; Zhdanov, M.

2018-02-01

This paper presents a parallel algorithm for frequency-domain acoustic modeling by the method of integral equations (IE). The algorithm is applied to seismic simulation. The IE method reduces the size of the problem but leads to a dense system matrix. A tolerable memory consumption and numerical complexity were achieved by applying an iterative solver, accompanied by an effective matrix-vector multiplication operation, based on the fast Fourier transform (FFT). We demonstrate that, the IE system matrix is better conditioned than that of the finite-difference (FD) method, and discuss its relation to a specially preconditioned FD matrix. We considered several methods of matrix-vector multiplication for the free-space and layered host models. The developed algorithm and computer code were benchmarked against the FD time-domain solution. It was demonstrated that, the method could accurately calculate the seismic field for the models with sharp material boundaries and a point source and receiver located close to the free surface. We used OpenMP to speed up the matrix-vector multiplication, while MPI was used to speed up the solution of the system equations, and also for parallelizing across multiple sources. The practical examples and efficiency tests are presented as well.

Solution of Large Systems of Linear Equations with Quadratic or Non-Quadratic Matrices and Deconvoiution of Spectra

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.

A New Algorithm for System of Integral Equations

Directory of Open Access Journals (Sweden)

Abdujabar Rasulov

2014-01-01

Full Text Available We develop a new algorithm to solve the system of integral equations. In this new method no need to use matrix weights. Beacause of it, we reduce computational complexity considerable. Using the new algorithm it is also possible to solve an initial boundary value problem for system of parabolic equations. To verify the efficiency, the results of computational experiments are given.

Parallel algorithms for computation of the manipulator inertia matrix

Science.gov (United States)

Amin-Javaheri, Masoud; Orin, David E.

1989-01-01

The development of an O(log2N) parallel algorithm for the manipulator inertia matrix is presented. It is based on the most efficient serial algorithm which uses the composite rigid body method. Recursive doubling is used to reformulate the linear recurrence equations which are required to compute the diagonal elements of the matrix. It results in O(log2N) levels of computation. Computation of the off-diagonal elements involves N linear recurrences of varying-size and a new method, which avoids redundant computation of position and orientation transforms for the manipulator, is developed. The O(log2N) algorithm is presented in both equation and graphic forms which clearly show the parallelism inherent in the algorithm.

Non-markovian boltzmann equation

International Nuclear Information System (INIS)

Kremp, D.; Bonitz, M.; Kraeft, W.D.; Schlanges, M.

1997-01-01

A quantum kinetic equation for strongly interacting particles (generalized binary collision approximation, ladder or T-matrix approximation) is derived in the framework of the density operator technique. In contrast to conventional kinetic theory, which is valid on large time scales as compared to the collision (correlation) time only, our approach retains the full time dependencies, especially also on short time scales. This means retardation and memory effects resulting from the dynamics of binary correlations and initial correlations are included. Furthermore, the resulting kinetic equation conserves total energy (the sum of kinetic and potential energy). The second aspect of generalization is the inclusion of many-body effects, such as self-energy, i.e., renormalization of single-particle energies and damping. To this end we introduce an improved closure relation to the Bogolyubov endash Born endash Green endash Kirkwood endash Yvon hierarchy. Furthermore, in order to express the collision integrals in terms of familiar scattering quantities (Mo/ller operator, T-matrix), we generalize the methods of quantum scattering theory by the inclusion of medium effects. To illustrate the effects of memory and damping, the results of numerical simulations are presented. copyright 1997 Academic Press, Inc

Equations of motion for a spectrum-generating algebra: Lipkin-Meshkov-Glick model

International Nuclear Information System (INIS)

Rosensteel, G; Rowe, D J; Ho, S Y

2008-01-01

For a spectrum-generating Lie algebra, a generalized equations-of-motion scheme determines numerical values of excitation energies and algebra matrix elements. In the approach to the infinite particle number limit or, more generally, whenever the dimension of the quantum state space is very large, the equations-of-motion method may achieve results that are impractical to obtain by diagonalization of the Hamiltonian matrix. To test the method's effectiveness, we apply it to the well-known Lipkin-Meshkov-Glick (LMG) model to find its low-energy spectrum and associated generator matrix elements in the eigenenergy basis. When the dimension of the LMG representation space is 10 6 , computation time on a notebook computer is a few minutes. For a large particle number in the LMG model, the low-energy spectrum makes a quantum phase transition from a nondegenerate harmonic vibrator to a twofold degenerate harmonic oscillator. The equations-of-motion method computes critical exponents at the transition point

Symmetries of the 2D magnetic particle imaging system matrix

International Nuclear Information System (INIS)

Weber, A; Knopp, T

2015-01-01

In magnetic particle imaging (MPI), the relation between the particle distribution and the measurement signal can be described by a linear system of equations. For 1D imaging, it can be shown that the system matrix can be expressed as a product of a convolution matrix and a Chebyshev transformation matrix. For multidimensional imaging, the structure of the MPI system matrix is not yet fully explored as the sampling trajectory complicates the physical model. It has been experimentally found that the MPI system matrix rows have symmetries and look similar to the tensor products of Chebyshev polynomials. In this work we will mathematically prove that the 2D MPI system matrix has symmetries that can be used for matrix compression. (paper)

An inherently parallel method for solving discretized diffusion equations

International Nuclear Information System (INIS)

Eccleston, B.R.; Palmer, T.S.

1999-01-01

A Monte Carlo approach to solving linear systems of equations is being investigated in the context of the solution of discretized diffusion equations. While the technique was originally devised decades ago, changes in computer architectures (namely, massively parallel machines) have driven the authors to revisit this technique. There are a number of potential advantages to this approach: (1) Analog Monte Carlo techniques are inherently parallel; this is not necessarily true to today's more advanced linear equation solvers (multigrid, conjugate gradient, etc.); (2) Some forms of this technique are adaptive in that they allow the user to specify locations in the problem where resolution is of particular importance and to concentrate the work at those locations; and (3) These techniques permit the solution of very large systems of equations in that matrix elements need not be stored. The user could trade calculational speed for storage if elements of the matrix are calculated on the fly. The goal of this study is to compare the parallel performance of Monte Carlo linear solvers to that of a more traditional parallelized linear solver. The authors observe the linear speedup that they expect from the Monte Carlo algorithm, given that there is no domain decomposition to cause significant communication overhead. Overall, PETSc outperforms the Monte Carlo solver for the test problem. The PETSc parallel performance improves with larger numbers of unknowns for a given number of processors. Parallel performance of the Monte Carlo technique is independent of the size of the matrix and the number of processes. They are investigating modifications to the scheme to accommodate matrix problems with positive off-diagonal elements. They are also currently coding an on-the-fly version of the algorithm to investigate the solution of very large linear systems

Complete factorisation and analytic solutions of generalized Lotka-Volterra equations

Science.gov (United States)

Brenig, L.

1988-11-01

It is shown that many systems of nonlinear differential equations of interest in various fields are naturally imbedded in a new family of differential equations. This family is invariant under nonlinear transformations based on the concept of matrix power of a vector. Each equation belonging to that family can be brought into a factorized canonical form for which integrable cases can be easily identified and solutions can be found by quadratures.

TEMPEST-2, Thermalization Program for Neutron Spectra and Multigroup Cross-Sections

International Nuclear Information System (INIS)

Gowins, G.

1984-01-01

Description of problem or function: TEMPEST2 is a neutron thermalization program based upon the Wigner-Wilkins approximation for light moderators and the Wilkins approximation for heavy moderators. A Maxwellian distribution may also be used. The model used may be selected as a function of energy. The second-order differential equations are integrated directly rather than transformed to the Riccati equation. The program provides microscopic and macroscopic cross-section averages over the thermal neutron spectrum

Stochastic Linear Quadratic Optimal Control Problems

International Nuclear Information System (INIS)

Chen, S.; Yong, J.

2001-01-01

This paper is concerned with the stochastic linear quadratic optimal control problem (LQ problem, for short) for which the coefficients are allowed to be random and the cost functional is allowed to have a negative weight on the square of the control variable. Some intrinsic relations among the LQ problem, the stochastic maximum principle, and the (linear) forward-backward stochastic differential equations are established. Some results involving Riccati equation are discussed as well

Matrix biorthogonal polynomials on the unit circle and non-Abelian Ablowitz-Ladik hierarchy

International Nuclear Information System (INIS)

Cafasso, Mattia

2009-01-01

Adler and van Moerbeke (2001 Commun. Pure Appl. Math. 54 153-205) described a reduction of the 2D-Toda hierarchy called the Toeplitz lattice. This hierarchy turns out to be equivalent to the one originally described by Ablowitz and Ladik (1975 J. Math. Phys. 16 598-603) using semidiscrete zero- curvature equations. In this paper, we obtain the original semidiscrete zero-curvature equations starting directly from the Toeplitz lattice and we generalize these computations to the matrix case. This generalization leads us to the semidiscrete zero-curvature equations for the non-Abelian (or multicomponent) version of the Ablowitz-Ladik equations (Gerdzhikov and Ivanov 1982 Theor. Math. Phys. 52 676-85). In this way, we extend the link between biorthogonal polynomials on the unit circle and the Ablowitz-Ladik hierarchy to the matrix case.

Exact many-body dynamics with stochastic one-body density matrix evolution

International Nuclear Information System (INIS)

Lacroix, D.

2004-05-01

In this article, we discuss some properties of the exact treatment of the many-body problem with stochastic Schroedinger equation (SSE). Starting from the SSE theory, an equivalent reformulation is proposed in terms of quantum jumps in the density matrix space. The technical details of the derivation a stochastic version of the Liouville von Neumann equation are given. It is shown that the exact Many-Body problem could be replaced by an ensemble of one-body density evolution, where each density matrix evolves according to its own mean-field augmented by a one-body noise. (author)

Theory of reflection reflection and transmission of electromagnetic, particle and acoustic waves

CERN Document Server

Lekner, John

2016-01-01

This book deals with the reflection of electromagnetic and particle waves by interfaces. The interfaces can be sharp or diffuse. The topics of the book contain absorption, inverse problems, anisotropy, pulses and finite beams, rough surfaces, matrix methods, numerical methods,  reflection of particle waves and neutron reflection. Exact general results are presented, followed by long wave reflection, variational theory, reflection amplitude equations of the Riccati type, and reflection of short waves. The Second Edition of the Theory of Reflection is an updated and much enlarged revision of the 1987 monograph. There are new chapters on periodically stratified media, ellipsometry, chiral media, neutron reflection and reflection of acoustic waves. The chapter on anisotropy is much extended, with a complete treatment of the reflection and transmission properties of arbitrarily oriented uniaxial crystals. The book gives a systematic and unified treatment reflection and transmission of electromagnetic and particle...

Distributed Secure Coordinated Control for Multiagent Systems Under Strategic Attacks.

Science.gov (United States)

Feng, Zhi; Wen, Guanghui; Hu, Guoqiang

2017-05-01

This paper studies a distributed secure consensus tracking control problem for multiagent systems subject to strategic cyber attacks modeled by a random Markov process. A hybrid stochastic secure control framework is established for designing a distributed secure control law such that mean-square exponential consensus tracking is achieved. A connectivity restoration mechanism is considered and the properties on attack frequency and attack length rate are investigated, respectively. Based on the solutions of an algebraic Riccati equation and an algebraic Riccati inequality, a procedure to select the control gains is provided and stability analysis is studied by using Lyapunov's method.. The effect of strategic attacks on discrete-time systems is also investigated. Finally, numerical examples are provided to illustrate the effectiveness of theoretical analysis.

On the quantum-mechanical Fokker-Planck and Kramers-Chandrasekhar equation

International Nuclear Information System (INIS)

Balazs, N.L.

1978-01-01

In the classical theory of Brownian motion the Langevin equation can be considered as an infinitesimal transformation between the coordinates and momenta of a Brownian particle, given probabilistically, since the impulse appearing is characterized by a Gaussian random process. This probabilistic infinitesimal transformation generates a streaming on the distribution function, expressed by the classical Fokker-Planck and Kramers-Chandrasekhar equations. If the laws obeyed by the Brownian particle are quantum mechanical, the Langevin equation can be reinterpreted as an operator relation expressing an infinitesimal transformation of these operators. Since the impulses are independent of the coordinates and momenta one can think of them as c numbers described by a Gaussian random process. The so resulting infinitesimal operator transformation induces a streaming on the density matrix. One may associate, according to Weyl, functions with operators. The function associated with the density matrix is the Wigner function. Expressing, then, these operator relations in terms of these functions the streaming can be expressed as a continuity equation of the Wigner function. It is found that in this parametrization the extra terms which appear are the same as in the classical theory, augmenting the usual Wigner equation. (Auth.)

Differential-difference equations associated with the fractional Lax operators

Energy Technology Data Exchange (ETDEWEB)

Adler, V E [LD Landau Institute for Theoretical Physics, 1A Ak. Semenov, Chernogolovka 142432 (Russian Federation); Postnikov, V V, E-mail: adler@itp.ac.ru, E-mail: postnikofvv@mail.ru [Sochi Branch of Peoples' Friendship University of Russia, 32 Kuibyshev str., 354000 Sochi (Russian Federation)

2011-10-14

We study integrable hierarchies associated with spectral problems of the form P{psi} = {lambda}Q{psi}, where P and Q are difference operators. The corresponding nonlinear differential-difference equations can be viewed as inhomogeneous generalizations of the Bogoyavlensky-type lattices. While the latter turn into the Korteweg-de Vries equation under the continuous limit, the lattices under consideration provide discrete analogs of the Sawada-Kotera and Kaup-Kupershmidt equations. The r-matrix formulation and several of the simplest explicit solutions are presented. (paper)

Effect of fillers on parameters of dry and swollen polymer matrix networks

Directory of Open Access Journals (Sweden)

Stojčeva-Radovanović Blaga

2002-01-01

Full Text Available The effect of nano- and micro- particle size of SiO2 on dry and swollen parameter network of the polymer matrix blends of acrylontrile-butadiene (NBR and chlorosulphonated polyethylene (CSM such as: volume and mass degree of swelling Rv and Rw; volume fraction of NBR-CSM polymer matrix in swollen gel V2 elasticity modulus G; interaction parameter between NBR-CSM polymer matrix and solvent λ and crossiinking density ν, was tested. The influence of nano-and micro- particle size of SiO2 on physical and mechanical properties, as well as effectiveness volume ratio of filiers in NBR-CSM polymer matrix at 300% elongation was tested using Einstein-Quth-Gold equation. The Kraus equation for swelling test of NBR-CSM polymer matrix containing nano- and micro- particle size of SiO2. Test results have shown that a greater interaction of nano-particie size of SiO2 with NBR-CSM polymer matrix, and possible chemical bonding, than the one of micro-silica was a consequence of a greater contact area. This results in better physical and mechanical properties.

Random-Effects Models for Meta-Analytic Structural Equation Modeling: Review, Issues, and Illustrations

Science.gov (United States)

Cheung, Mike W.-L.; Cheung, Shu Fai

2016-01-01

Meta-analytic structural equation modeling (MASEM) combines the techniques of meta-analysis and structural equation modeling for the purpose of synthesizing correlation or covariance matrices and fitting structural equation models on the pooled correlation or covariance matrix. Both fixed-effects and random-effects models can be defined in MASEM.…

On the fusion matrix of the N=1 Neveu-Schwarz blocks

OpenAIRE

Hadasz, Leszek

2007-01-01

We propose an exact form of the fusion matrix of the Neveu-Schwarz blocks that appear in a correlation function of four super-primary fields. Orthogonality relation satisfied by this matrix is equivalent to the bootstrap equation for the four-point super-primary correlator in the N=1 supersymmetric Liouville theory.

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

Science.gov (United States)

Szyszka, Barbara

2013-10-01

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

Chew-Low equations as Cremoma transformations

International Nuclear Information System (INIS)

Rerikh, K.V.

1982-01-01

The Chew-Low equations for the p-wave pion-nucleon scattering with the crossing-symmetry matrix (3x3) are investigated in their well-known formulation as a system of nonlinear difference equations. These equations interpreted as geometrical transformations are shown to be a special case of the Cremona transformaions. Using the properties of the Cremona transformations we obtain the general 3-parametric functional equation on invariant algebraic and nonalgebraic curves in the space solutions of the Chew- Low equations. It is proved that there exists only one invariant algebraic curve, the parabola corresponding to the well-known solution. Analysis of the general functional equation on invariant nonalgebraic curves makes it possible to select in addition to this parabola 3 invariant forms defining implicitly 3 nonalgebraic curves and to concretize for them the general equation by means of fixing the parameters. From the transformational properties of the invariant forms with respect to the Cremona transformations, there follows an important result that the ration of these forms in proper powers is the general integral of the nonlinear system of the Chew-Low equations, which is an even antiperiodic function. The structure of the second general integral is given and the functional equations which determinne this integral are presented [ru

A random matrix model of relaxation

International Nuclear Information System (INIS)

Lebowitz, J L; Pastur, L

2004-01-01

We consider a two-level system, S 2 , coupled to a general n level system, S n , via a random matrix. We derive an integral representation for the mean reduced density matrix ρ(t) of S 2 in the limit n → ∞, and we identify a model of S n which possesses some of the properties expected for macroscopic thermal reservoirs. In particular, it yields the Gibbs form for ρ(∞). We also consider an analog of the van Hove limit and obtain a master equation (Markov dynamics) for the evolution of ρ(t) on an appropriate time scale

On the fusion matrix of the N = 1 Neveu-Schwarz blocks

International Nuclear Information System (INIS)

Hadasz, Leszek

2007-01-01

We propose an exact form of the fusion matrix of the Neveu-Schwarz blocks that appear in a correlation function of four super-primary fields. Orthogonality relation satisfied by this matrix is equivalent to the bootstrap equation for the four-point super-primary correlator in the N = 1 supersymmetric Liouville field theory

On the fusion matrix of the N = 1 Neveu-Schwarz blocks

Science.gov (United States)

Hadasz, Leszek

2007-12-01

We propose an exact form of the fusion matrix of the Neveu-Schwarz blocks that appear in a correlation function of four super-primary fields. Orthogonality relation satisfied by this matrix is equivalent to the bootstrap equation for the four-point super-primary correlator in the N = 1 supersymmetric Liouville field theory.

The finite element response matrix method

International Nuclear Information System (INIS)

Nakata, H.; Martin, W.R.

1983-02-01

A new technique is developed with an alternative formulation of the response matrix method implemented with the finite element scheme. Two types of response matrices are generated from the Galerkin solution to the weak form of the diffusion equation subject to an arbitrary current and source. The piecewise polynomials are defined in two levels, the first for the local (assembly) calculations and the second for the global (core) response matrix calculations. This finite element response matrix technique was tested in two 2-dimensional test problems, 2D-IAEA benchmark problem and Biblis benchmark problem, with satisfatory results. The computational time, whereas the current code is not extensively optimized, is of the same order of the well estabilished coarse mesh codes. Furthermore, the application of the finite element technique in an alternative formulation of response matrix method permits the method to easily incorporate additional capabilities such as treatment of spatially dependent cross-sections, arbitrary geometrical configurations, and high heterogeneous assemblies. (Author) [pt

Computer Applications in Balancing Chemical Equations.

Science.gov (United States)

Kumar, David D.

2001-01-01

Discusses computer-based approaches to balancing chemical equations. Surveys 13 methods, 6 based on matrix, 2 interactive programs, 1 stand-alone system, 1 developed in algorithm in Basic, 1 based on design engineering, 1 written in HyperCard, and 1 prepared for the World Wide Web. (Contains 17 references.) (Author/YDS)

Bright, dark and singular optical solitons in a cascaded system

International Nuclear Information System (INIS)

Zhou, Qin; Zhu, Qiuping; Yu, Hua; Liu, Yaxian; Wei, Chun; Yao, Ping; Bhrawy, Ali H; Biswas, Anjan

2015-01-01

This work studies nonlinear dynamics of optical solitons in a cascaded system with Kerr law nonlinearity and spatio-temporal dispersion. The mathematical model that describes the propagation of optical solitons through a cascaded system is given by the vector-coupled nonlinear Schrödinger equation. It is investigated analytically using three integration algorithms. The Jacobian elliptic equation expansion method, Bernoulli equation expansion approach and Riccati equation expansion scheme are the integration tools of this model that are recruited to extract singular, bright and dark solitons. The restrictions that need to hold for the existence of these solitons are derived. (paper)

[Evaluation of the Peusner's coefficients matrix for polymeric membrane and ternary non-electrolyte solutions].

Science.gov (United States)

Jasik-Slęzak, Jolanta; Slęzak-Prochazka, Izabella; Slęzak, Andrzej

2014-01-01

A system of network forms of Kedem-Katchalsky (K-K) equations for ternary non-electrolyte solutions is made of eight matrix equations containing Peusner's coefficients R(ij), L(ij), H(ij), W(ij), K(ij), N(ij), S(ij) or P(ij) (i, j ∈ {1, 2, 3}). The equations are the result of symmetric or hybrid transformation of the classic form of K-K equations by the use of methods of Peusner's network thermodynamics (PNT). Calculating concentration dependences of the determinant of Peusner's coefficients matrixes R(ij), L(ij), H(ij), W(ij), S(ij), N(ij), K(ij) and P(ij) (i, j ∈ {1, 2, 3}). The material used in the experiment was a hemodialysis Nephrophan membrane with specified transport properties (L(p), σ, Ω) in aqueous glucose and ethanol solution. The method involved equations for determinants of the matrixes coefficients R(ij), L(ij), H(ij), W(ij), S(ij), N(ij), K(ij) or P(ij) (i, j ∈ {1, 2, 3}). The objective of calculations were dependences of determinants of Peusner's coeffcients matrixes R(ij), L(ij), H(ij), W(ij), S(ij), N(ij), K(ij) or P(ij) (i, j ∈ {1, 2, 3}) within the conditions of solution homogeneity upon an average concentration of one component of solution in the membrane (C1) with a determined value of the second component (C2). The method of calculating the determinants of Peusner's coeffcients matrixes R(ij), L(ij), H(ij), W(ij), S(ij), N(ij), K(ij) or P(ij) (i, j ∈ {1, 2, 3}) is a new tool that may be applicable in studies on membrane transport. Calculations showed that the coefficients are sensitive to concentration and composition of solutions separated by a polymeric membrane.

The optimal solution of a non-convex state-dependent LQR problem and its applications.

Directory of Open Access Journals (Sweden)

Xudan Xu

Full Text Available This paper studies a Non-convex State-dependent Linear Quadratic Regulator (NSLQR problem, in which the control penalty weighting matrix [Formula: see text] in the performance index is state-dependent. A necessary and sufficient condition for the optimal solution is established with a rigorous proof by Euler-Lagrange Equation. It is found that the optimal solution of the NSLQR problem can be obtained by solving a Pseudo-Differential-Riccati-Equation (PDRE simultaneously with the closed-loop system equation. A Comparison Theorem for the PDRE is given to facilitate solution methods for the PDRE. A linear time-variant system is employed as an example in simulation to verify the proposed optimal solution. As a non-trivial application, a goal pursuit process in psychology is modeled as a NSLQR problem and two typical goal pursuit behaviors found in human and animals are reproduced using different control weighting [Formula: see text]. It is found that these two behaviors save control energy and cause less stress over Conventional Control Behavior typified by the LQR control with a constant control weighting [Formula: see text], in situations where only the goal discrepancy at the terminal time is of concern, such as in Marathon races and target hitting missions.

Dynamic Modeling and Analysis of the Large-Scale Rotary Machine with Multi-Supporting

Directory of Open Access Journals (Sweden)

Xuejun Li

2011-01-01

Full Text Available The large-scale rotary machine with multi-supporting, such as rotary kiln and rope laying machine, is the key equipment in the architectural, chemistry, and agriculture industries. The body, rollers, wheels, and bearings constitute a chain multibody system. Axis line deflection is a vital parameter to determine mechanics state of rotary machine, thus body axial vibration needs to be studied for dynamic monitoring and adjusting of rotary machine. By using the Riccati transfer matrix method, the body system of rotary machine is divided into many subsystems composed of three elements, namely, rigid disk, elastic shaft, and linear spring. Multiple wheel-bearing structures are simplified as springs. The transfer matrices of the body system and overall transfer equation are developed, as well as the response overall motion equation. Taken a rotary kiln as an instance, natural frequencies, modal shape, and response vibration with certain exciting axis line deflection are obtained by numerical computing. The body vibration modal curves illustrate the cause of dynamical errors in the common axis line measurement methods. The displacement response can be used for further measurement dynamical error analysis and compensation. The response overall motion equation could be applied to predict the body motion under abnormal mechanics condition, and provide theory guidance for machine failure diagnosis.

The Gibbs-Thomson equation for a spherical coherent precipitate with applications to nucleation

International Nuclear Information System (INIS)

Rottman, C.; Voorhees, P.W.; Johnson, W.C.

1988-01-01

The conditions for interfacial thermodynamic equilibrium form the basis for the derivation of a number of basic equations in materials science, including the various forms of the Gibbs-Thomson equation. The equilibrium conditions pertaining to a curved interface in a two-phase fluid system are well-known. In contrast, the conditions for thermodynamic equilibrium at a curved interface in nonhydrostatically stressed solids have only recently been examined. These conditions can be much different from those at a fluid interface and, as a result, the Gibbs-Thomson equation appropriate to coherent solids is likely to be considerably different from that for fluids. In this paper, the authors first derive the conditions necessary for thermodynamic equilibrium at the precipitate-matrix interface of a coherent spherical precipitate. The authors' derivation of these equilibrium conditions includes a correction to the equilibrium conditions of Johnson and Alexander for a spherical precipitate in an isotropic matrix. They then use these conditions to derive the dependence of the interfacial precipitate and matrix concentrations on precipitate radius (Gibbs-Thomson equation) for a such a precipitate. In addition, these relationships are then used to calculate the critical radius for the nucleation of a coherent misfitting precipitate

The QCD spin chain S matrix

International Nuclear Information System (INIS)

Ahn, Changrim; Nepomechie, Rafael I.; Suzuki, Junji

2008-01-01

Beisert et al. have identified an integrable SU(2,2) quantum spin chain which gives the one-loop anomalous dimensions of certain operators in large N c QCD. We derive a set of nonlinear integral equations (NLIEs) for this model, and compute the scattering matrix of the various (in particular, magnon) excitations

Combinatorics of Generalized Bethe Equations

Science.gov (United States)

Kozlowski, Karol K.; Sklyanin, Evgeny K.

2013-10-01

A generalization of the Bethe ansatz equations is studied, where a scalar two-particle S-matrix has several zeroes and poles in the complex plane, as opposed to the ordinary single pole/zero case. For the repulsive case (no complex roots), the main result is the enumeration of all distinct solutions to the Bethe equations in terms of the Fuss-Catalan numbers. Two new combinatorial interpretations of the Fuss-Catalan and related numbers are obtained. On the one hand, they count regular orbits of the permutation group in certain factor modules over {{Z}^M}, and on the other hand, they count integer points in certain M-dimensional polytopes.