回归分析方法在油田抽油机井管理中的应用

本文在抽油机井大量现场测试数据的基础上，建立了抽油机井系统效率、排量系数这两个重要指标与诸影响因素之间的回归方程。以此为基础，作者提出了一种抽油机井系统效率预测及参数优化设计的新方法，该方法在油田中应用已表现出良好的经济效益。     

Intelligent control of chaos using linear feedback controller and neural network identifier

A method for controlling chaos when the mathematical model of the system is unknown is presented in this paper. The controller is designed by the pole placement algorithm which provides a linear feedback control method. For calculating the feedback gain, a neural network is used for identification of the system from which the Jacobian of the system in its fixed point can be approximated. The weights of the neural network are adjusted online by the gradient descent algorithm in which the difference between the system output and the network output is considered as the error to be decreased. The method is applied on both discrete-time and continuous-time systems. For continuous-time systems, equivalent discrete-time systems are constructed by using the Poincare map concept. Two discrete-time systems and one continuous-time system are tested as examples for simulation and the results show good functionality of the proposed method. It can be concluded that the chaos in systems with unknown dynamics may be eliminated by the presented intelligent control system based on pole placement and neural network.     

非线性系统运动稳定性的一种判据

对于自治的非线性系统来说,只要其线性部分系数矩阵的特征值不属于临界情形,其无扰运动在其足够小的邻域内的稳定性完全可以由其线性部分的特征值确定.关于线性系统的稳定性,已有不少简单易行的判别方法,而关于非线性系统的稳定性,很多数学家和力学家作了大量的研究工作;但大都是针对特殊类型的非线性系统解决了一些问题,直到现在为止,还没有普遍适用于任何的非线性系统的简单易行的判别方法.本文所给的是判别非线性系统稳定性的充要条件,常用的克拉索夫斯基方法只是这一方法的一个特例[1],[2].     

Partitioning of invariant imbedding systems

The choice of partitioning the system matrix for a system of N linear ordinary differential equations may determine the ease or difficulty of obtaining a solution by the invariant imbedding method of Scott. This paper shows how the configuration and partitioning of the system matrix is reflected in the fundamental matrix. The partitioning of the fundamental matrix is the key to the ease or difficulty of obtaining a solution. If the fundamental matrix is known for a given system matrix configuration and partitioning, then the fundamental matrix associated with a new system matrix configuration may be derived by the same row and column interchanges that transformed the old system matrix into the new system matrix. The fundamental matrix for the new system matrix does not have to be recalculated anew from the Kronecker delta initial conditions.     

An averaging theorem for two-point boundary value problems with applications to optimal control

An asymptotic result is obtained for a two-point boundary value problem for a vector system of nonlinear ordinary differential equations involving “fast” and “slow” inputs. The asymptotically limiting system is obtained by an averaging procedure. Using this result, an approximate analysis of the original system may be carried out by considering two lower-order systems each involving only one time scale. It is shown that some optimal control problems for systems with multiple time scales may be analyzed by this method.     

A Method for the Numerical Solution of Boundary Value Problems Governed by Second-order Elliptic Systems

A method for the solution of boundary value problems governedby a system of second-order elliptic partial differential equationsis derived. The method is obtained by reducing the boundaryvalue problem to integral equations which may be solved numerically.     

Stability analysis of controlled multiple-link robotic manipulator systems with time delays

Controlled robotic manipulator systems, normally operating in an acceptable manner, may become unstable, when time delays exist in the feedback loop. This instability may cause oscillations and thus deteriorate the system performance. In this paper we introduce a method to analyze the performance of robotic manipulators with small time delays in the feedback loop. Specifically, we may predict the limits of a stable system Operation, i.e., the robustness and expected system error bounds, as a function of the time delay and system parameters. The analysis method uses nonlinear control theory concepts, such as Lyapunov's second method, and takes advantage of norms to establish operational bounds. Previous research of time-delay systems analysis is also reviewed, in general, and pertaining to robotic controls, in particular. The results of the analysis can be displayed graphically and may be used to adjust the controller gains and trajectories to ensure satisfactory operation. Predictably, the system response is slowed as the time delay increases. The validity of the analysis is demonstrated with a numerical example of a two-link manipulator using a computed-torque controller. The time responses of the example, run on a computer-generated model, supports the analysis results and predicted performance.     

Mathematical modelling of a quasi-linear feedback system

This paper develops a method for system identification of a two input-one output quasi-linear system modelling the dynamics of a non-linear element of a feedback system. The identification method, which may be performed in either the time domain or the frequency domain, is designed to use normal operating records without the need to introduce special test inputs.     

Determination of the Chain-Like Non-Linear Multi-Degree-of-Freedom Systems Constant Parameters under Dynamical Complex Loads

A determination of mechanical properties of the real material systems under dynamical excitations is a very complex problem. The real system can be described by continuous model as well as a discrete model with a finite number of degree-of-freedom. Some mechanical systems like beams, rods, strings, tow-lines etc. have so-called chain-like configuration and may be treated as the dynamical one-dimensional cascade systems. In this paper, an identification method of internal interaction force, which can be arbitrary selected in discrete chain-like non-linear system, is presented. The method was tested on example of a two-mass system in which the interaction forces were modeled by spring-damper elements in a non-parallel arrangement. In experiment some non-sinusoidal complex loads were applied. (© 2009 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Eigenvalue estimates for a preconditioned Galerkin matrix arising from mixed finite element discretizations of viscous incompressible flows

The article deals with Galerkin matrices arising with finite element discretizations of the Navier–Stokes system. Usually these matrices are indefinite and nonsymmetric. They have to be preconditioned if a related linear system is to be solved efficiently by an iterative method. We consider preconditioning by a pressure mass matrix. It is shown how upper and lower bounds of the eigenvalues of a preconditioned Galerkin matrix may be found by variational arguments.     

一类带小时滞的非线性快慢系统的渐近解

研究了一类带小时滞的非线性快慢系统的初始值问题,在一定假设条件下,利用奇异摄动理论和校正函数法构造了该问题的形式渐近解,并利用微分不等式理论证明了渐近解的一致有效性.最后进行了算例分析,结果显示时滞能对快慢系统产生重要影响,并表明所述摄动方法是一个行之有效的近似解析方法.从而,可以利用得到的渐近解对系统的动力学行为进行更深层次地分析与研究.     

New exact solutions for a generalization of the Korteweg-de Vries equation (KdV6)

The generalized tanh-coth method is used to construct periodic and soliton solutions for a new integrable system, which has been derived from an integrable sixth-order nonlinear wave equation (KdV6). The system is formed by two equations. One of the equations may be considered as a Korteweg-de Vries equation with a source and the second equation is a third-order linear differential equation.     

Componentwise fast convergence in the solution of full-rank systems of nonlinear equations

The asymptotic convergence of parameterized variants of Newton’s method for the solution of nonlinear systems of equations is considered. The original system is perturbed by a term involving the variables and a scalar parameter which is driven to zero as the iteration proceeds. The exact local solutions to the perturbed systems then form a differentiable path leading to a solution of the original system, the scalar parameter determining the progress along the path. A path-following algorithm, which involves an inner iteration in which the perturbed systems are approximately solved, is outlined. It is shown that asymptotically, a single linear system is solved per update of the scalar parameter. It turns out that a componentwise Q-superlinear rate may be attained, both in the direct error and in the residuals, under standard assumptions, and that this rate may be made arbitrarily close to quadratic. Numerical experiments illustrate the results and we discuss the relationships that this method shares with interior methods in constrained optimization. Received: September 8, 2000 / Accepted: September 17, 2001?Published online February 14, 2002     

Stability analysis of impulsive control systems

This paper studies impulsive control systems. Several stability criteria are established by employing the method of Lyapunov functions. These criteria may be used for impulsive feedback control design. As an application, impulsive control of the Lorenz chaotic system is discussed. Numerical experiments are carried out for the control of the Lorenz system. It is shown that small and frequent impulses need to be used in order to stabilize the Lorenz system.     

基于Udwadia-Kalaba理论的Hamel嵌入法研究

Hamel嵌入法直接将约束嵌入到非约束运动的动能中去,从而避免使用Lagrange(拉格朗日)乘子.但这个简单、直观的方法却并不总是正确.Hamel认为这种方法可能导致错误的结果,然而他并没有给出Hamel嵌入法正确性的适用条件.在利用Udwadia Kalaba理论的基础上,提出了Hamel嵌入法成立的充要条件;指出了Rosenberg在Hamel嵌入法正确性研究中的不足,通过给出的具体算例可以看出,在完整约束下Hamel嵌入法可能不正确,而在非完整约束下也可能得出正确的结果;理论和实例分析表明,Hamel嵌入法是否成立除了与约束有关以外还与系统模型相关.     

A dynamic Stackelberg duopoly model with different strategies

In this paper, a duopoly Stackelberg model of competition on output is formulated. The firms announce plan products sequentially in planning phase and act simultaneously in production phase. For the duopoly Stackelberg model, a nonlinear dynamical system which describes the time evolution with different strategies is analyzed. We present results on existence, stability and local bifurcations of the equilibrium points. Numerical simulations demonstrate that the system with varying model parameters may drive to chaos and the loss of stability may be caused by period doubling bifurcations. It is also shown that the state variables feedback and parameter variation method can be used to keep the system from instability and chaos.     

Life Expectancy of Transient Microchaotic Behaviour

The application of digital control may lead to so-called transient chaotic behaviour. In the present paper, we analyse a simple model of a digitally controlled mechanical system, which may create such vibrations. As a consequence of the digital effects, i.e., the sampling and the round-off error, the behaviour of this system can be described by a one-dimensional piecewise linear map. The lifetime of chaotic transients is usually characterized by the so-called escape rate. In the literature, the reciprocal of the escape rate is considered to be the expected duration of the transient chaotic phenomenon. We claim that this approach is not always fruitful, and present a different way of calculating the mean lifetime in the case of one-dimensional piecewise linear maps. Our method might also be used to solve diffusion problems in one-dimensional models of periodic arrays.     

A new computational method for solving fully fuzzy linear systems of triangular fuzzy numbers

In this paper, a new computational method is proposed to solve fully fuzzy linear systems (FFLS) of triangular fuzzy numbers based on the computation of row reduced echelon form for solving the crisp linear system of equations. The method is illustrated by solving three numerical examples. As compared to the existing methods, the proposed method is easy to understand and to apply for solving FFLS occurring in real life situations and scientific applications. The primary advantage of the proposed method is that, by using it, the consistency of the FFLS can be checked easily and nature of the solutions of an FFLS can also be obtained which may be unique or infinite.     

Eigenvalue bounds for the Schur complement with a pressure convection–diffusion preconditioner in incompressible flow computations

If the stationary Navier–Stokes system or an implicit time discretization of the evolutionary Navier–Stokes system is linearized by a Picard iteration and discretized in space by a mixed finite element method, there arises a saddle point system which may be solved by a Krylov subspace method or an Uzawa type approach. For each of these resolution methods, it is necessary to precondition the Schur complement associated to the saddle point problem in question. In the work at hand, we give upper and lower bounds of the eigenvalues of this Schur complement under the assumption that it is preconditioned by a pressure convection–diffusion matrix.