Multi-order Arnoldi-based model order reduction of second-order time-delay systems

In this paper, we discuss the Krylov subspace-based model order reduction methods of second-order systems with time delays, and present two structure-preserving methods for model order reduction of these second-order systems, which avoid to convert the second-order systems into first-order ones. One method is based on a Krylov subspace by using the Taylor series expansion, the other method is based on the Laguerre series expansion. These two methods are used in the multi-order Arnoldi algorithm to construct the projection matrices. The resulting reduced models can not only preserve the structure of the original systems, but also can match a certain number of approximate moments or Laguerre expansion coefficients. The effectiveness of the proposed methods is demonstrated by two numerical examples.     

Compact thermoelectric coupled models of micromachined thermal sensors using trajectory piecewise-linear model order reduction

Numerical simulation has become an effective way to design and optimize micromachined thermal sensors. To improve the speed of development process, fast simulation is indispensable. This paper investigates compact models, which are accurate low order representations of high order finite element models, for micromachined thermal sensors. Because the thermal field and the electric field act on and influence each other, the compact models should be established from the thermoelectric coupled full models. Thermoelectric coupling and temperature dependent resistivity make the problems strong nonlinearity. Therefore a powerful nonlinear model order reduction method, named trajectory piecewise-linear (TPWL) method, is employed. Its core idea is approximating the nonlinear model piecewise-linearly along a training trajectory. The performance of the TPWL method and fast TPWL method are compared. And the effects of linearization point number and local reduced basis order on accuracy, efficiency and size of the TPWL compact models are studied. Moreover, the expandability of the TPWL compact models is also discussed. Results show that the TPWL compact models are suitable for the design and optimization of micromachined thermal sensors.     

基于Krylov-Schur重启技术的Arnoldi模型降阶方法

Krylov子空间模型降阶方法是模型降阶中的典型方法之一,Arnoldi模型降阶方法是这类方法中的一类基本方法。运用重正交化的Arnoldi算法得到[r]步Arnoldi分解;执行Krylov-Schur重启过程,导出基于Krylov-Schur重启技术的Arnoldi模型降阶方法。运用此方法对大规模线性时不变系统进行降阶,得到具有较高近似精度的稳定的降阶系统,从而改善了Krylov子空间降阶方法不能保持降阶系统稳定性的不足。数值算例验证了此方法是行之有效的。     

Model order reduction of MIMO bilinear systems by multi-order Arnoldi method

In this paper, we present a time domain model order reduction method for multi-input multi-output (MIMO) bilinear systems by general orthogonal polynomials. The proposed method is based on a multi-order Arnoldi algorithm applied to construct the projection matrix. The resulting reduced model can match a desired number of expansion coefficient terms of the original system. The approximate error estimate of the reduced model is given. And we also briefly discuss the stability preservation of the reduced model in some cases. Additionally, in combination with Krylov subspace methods, we propose a two-sided projection method to generate reduced models which capture properties of the original system in the time and frequency domain simultaneously. The effectiveness of the proposed methods is demonstrated by two numerical examples.     

互连线高效时域梯形差分模型降阶算法

为了进一步提高现有互连电路模型降阶方法的精度和效率,提出一种基于时域梯形法差分的互连线模型降阶方法.首先将互连电路的时域方程用梯形法差分离散后获得一种关于状态变量的递推关系,形成了一个非齐次Krylov子空间;然后利用非齐次Arnoldi算法求得非齐次Krylov子空间的正交基,再通过正交基对原始系统进行投影得到降阶系统.该算法可以保证时域差分后降阶系统和原始系统的状态变量在离散时间点的匹配,保证时域降阶精度,同时也保证了降阶过程的数值稳定性及降阶系统的无源性.与现有的时域模型降阶方法相比,文中算法可降低计算复杂度;与频域降阶方法相比,由于避免了时频域转换误差,其在时域具有更高的精度.     

Identification of multi-model LPV models with two scheduling variables

In order to model complex industrial processes, this work studies the identification of linear parameter varying (LPV) models with two scheduling variables. The LPV model is parameterized as blended linear models, which is also called multi-model structure. Several weighting functions, linear, polynomial and Gaussian functions, are used and compared. The usefulness of the method is tested using a high purity distillation column model in a case study. The case study shows that a good fit of identification data is not enough to verify model quality and can even be misleading in nonlinear process identification; other measures related to process knowledge should be used in model validation. The case study also shows that commonly used LPV model based on parameter interpolation can fail for the high purity distillation column. Finally, several pitfalls in nonlinear process identification are pointed out.     

Model order reduction of second-order systems with nonlinear stiffness using Krylov subspace methods and their symmetric transfer functions

In this investigation, Model Order Reduction (MOR) of second-order systems having cubic nonlinearity in stiffness is developed for the first time using Krylov subspace methods and the associated symmetric transfer functions. In doing so, new second-order Krylov subspaces will be defined for MOR procedure which avoids the need to transform the second-order system to its state space form and thus the main characteristics of the second-order system such as symmetry and positive definiteness of mass and stiffness matrices will be preserved. To show the efficacy of the presented method, three examples will be considered as practical case studies. The first example is a nonlinear shear-beam building model subjected to a seismic disturbance. The second and third examples are nonlinear longitudinal vibration of a rod and vibration of a cantilever beam resting on a nonlinear elastic foundation, respectively. Simulation results in all cases show good accuracy of the vibrational response of the reduced order models when compared with the original ones while reducing the computational load.     

Model-order reductions for MIMO systems using global Krylov subspace methods

This paper presents theoretical foundations of global Krylov subspace methods for model order reductions. This method is an extension of the standard Krylov subspace method for multiple-inputs multiple-outputs (MIMO) systems. By employing the congruence transformation with global Krylov subspaces, both one-sided Arnoldi and two-sided Lanczos oblique projection methods are explored for both single expansion point and multiple expansion points. In order to further reduce the computational complexity for multiple expansion points, adaptive-order multiple points moment matching algorithms, or the so-called rational Krylov space method, are also studied. Two algorithms, including the adaptive-order rational global Arnoldi (AORGA) algorithm and the adaptive-order global Lanczos (AOGL) algorithm, are developed in detail. Simulations of practical dynamical systems will be conducted to illustrate the feasibility and the efficiency of proposed methods.     

An integrated state space partition and optimal control method of multi-model for nonlinear systems based on hybrid systems

Multilinear model approach turns out to be an ideal candidate for dealing with nonlinear systems control problem. However, how to identify the optimal active state subspace of each linear subsystem is an open problem due to that the closed-loop performance of nonlinear systems interacts with these subspaces ranges. In this paper, a new systematic method of integrated state space partition and optimal control of multi-model for nonlinear systems based on hybrid systems is initially proposed, which can deal with the state space partition and associated optimal control simultaneously and guarantee an overall performance of nonlinear systems consequently. The proposed method is based on the framework of hybrid systems which synthesizes the multilinear model, produced by nonlinear systems, in a unified criterion and poses a two-level structure. At the upper level, the active state subspace of each linear subsystem is determined under the optimal control index of a hybrid system over infinite horizon, which is executed off-line. At the low level, the optimal control is implemented online via solving the optimal control of hybrid system over finite horizon. The finite horizon optimal control problem is numerically computed by simultaneous method for speeding up computation. Meanwhile, the model mismatch produced by simultaneous method is avoided by using the strategy of receding-horizon. Simulations on CSTR (Continuous Stirred Tank Reactor) confirm that a superior performance can be obtained by using the presented method.     

基于二次仿射传播聚类的非线性系统多模型LSSVM建模

针对复杂非线性系统单模型建模存在计算量大和精度差的问题,提出一种采用仿射传播聚类的多模型LSSVM建模方法,通过仿射传播聚类对样本数据按输入集和输出集二次聚类划分,并分别建立LSSVM子模型,非线性系统模型通过子模型加权合成.将该方法应用于两电机变频调速系统的张力和速度模型辨识,仿真结果表明,该建模方法具有较高的精度,能准确拟合系统的非线性特性.     

A parameterised model order reduction method for parametric systems based on Laguerre polynomials

This paper presents a Laguerre polynomials-based parametrised model order reduction method for the parametric system in time domain. The method allows that the parametric dependence in system matrices is nonaffine. The method is presented via reducing an approximate polynomial parametric system based on Taylor expansion and Laguerre polynomials, resulting in a parametric reduced system that can accurately approximate the time response of the original parametric system over a wide range of parameter. The reduced parametric system obtained by proposed method can be implemented by two algorithms. Algorithm 1 is a direct way that is suitable for single-input multi-output parametric systems. Algorithm 2 is presented based on a connection to the Krylov subspace, which is efficient and suitable for multi-input multi-output parametric systems. The effectiveness of the proposed method is illustrated with two benchmarks in practical applications.     

A Krylov‐subspace technique for the simulation of integrated RF/microwave subsystems driven by digitally modulated carriers

This article introduces a new approach to the analysis of nonlinear RF/microwave systems or subsystems described at the circuit level and excited by sinusoidal carriers modulated by arbitrary baseband signals. The circuit is simulated by a sequence of harmonic‐balance analyses based on a Krylov‐subspace method driven by an inexact Newton loop, and suitably modified to account for coupling with a finite number of preceding time instants. The Jacobian matrix of the nonlinear solving system is computed by an exact algorithm, and its structure is shown to be very well suited for the application of Krylov‐subspace techniques such as the GMRES method. In this way, problems with many millions of nodal unknowns may be efficiently tackled at the workstation level. ©1999 John Wiley & Sons, Inc. Int J RF and Microwave CAE 9: 490–505, 1999.     

Control-relevant nonlinearity measure and integrated multi-model control

A control-relevant nonlinearity measure (CRNM) method is proposed based on the gap metric and the gap metric stability margin to measure the nonlinear degree of a system once a linear control strategy is selected. Supported by the CRNM method, an integrated multi-model control framework is developed, in which the multi-model decomposition and local controller design are closely integrated, model redundancy is avoided, computational load is reduced, and dependency on a prior knowledge is reduced. Besides, a 1/δ gap-based weighting method is put forward to combine the local controllers. On one hand, the 1/δ gap-based weighting method has merely one tuning parameter and can be computed off-line; on the other hand, it is sensitive to the tuning parameter, flexible and easy to tune. Two continuous stirred tank reactor (CSTR) systems are investigated. Closed-loop simulations validate the effectiveness and benefits of the proposed integrated multi-model control approach based on CRNM.     

Model-order reduction of large-scale second-order MIMO dynamical systems via a block second-order Arnoldi method

In this paper, we present a structure-preserving model-order reduction method for solving large-scale second-order MIMO dynamical systems. It is a projection method based on a block second-order Krylov subspace. We use the block second-order Arnoldi (BSOAR) method to generate an orthonormal basis of the projection subspace. The reduced system preserves the second-order structure of the original system. Some theoretical results are given. Numerical experiments report the effectiveness of this method.     

A robust multi-model predictive controller for distributed parameter systems

In this work a robust nonlinear model predictive controller for nonlinear convection-diffusion-reaction systems is presented. The controller makes use of a collection of reduced order approximations of the plant (models) reconstructed on-line by projection methods on proper orthogonal decomposition (POD) basis functions. The model selection and model update step is based on a sufficient condition that determines the maximum allowable process-model mismatch to guarantee stable control performance despite process uncertainty and disturbances. Proofs on the existence of a sequence of feasible approximations and control stability are given.Since plant approximations are built on-line based on actual measurements, the proposed controller can be interpreted as a multi-model nonlinear predictive control (MMPC). The performance of the MMPC strategy is illustrated by simulation experiments on a problem that involves reactant concentration control of a tubular reactor with recycle.     

非线性系统预测控制的多模型方法

本文在非线性系统的线性化多模型基础上,引入多模型参考轨迹逼近期望轨迹,提出了一 种非线性系统预控制的多模型方法.仿真结果表明,这种方法是有效的.     

A sequential LQG approach to nonlinear tracking problem

A new approach, named the sequential LQG solution for nonlinear and nonstationary systems tracking problem, has been proposed in the article. Proposed method for tracking of a desired reference trajectory uses a state space model of the multivariable nonlinear/nonstationary system. The method involves three steps. The first step is the design of nominal trajectory using the predictive control technique. The second step is the sequential linearization of the nonlinear system around the fixed operating points, chosen in accordance to the plant dynamics changes. The third step involves the choice of the weighting matrices in the classical LQG controller design for the sequence of linearized models. The feasibility of the proposed method has been demonstrated through its application to the control of aircraft, represented by the six-degree-of-freedom model around the prespecified nominal trajectory.     

中心刚体-柔性梁刚柔耦合动力学模型降阶研究

有限单元法被广泛的采用来描述柔性体的弹性变形,然而有限元节点坐标数目庞大,将会给动力学方程求解带来巨大的计算负担.如何降低柔性体的自由度,是当前柔性多体系统动力学研究的一个重要命题.本文以中心刚体-柔性梁系统为例,采用Krylov方法和模态方法进行降价.然后分别采用有限元全模型、Krylov降阶模型和模态降阶模型,对中心刚体-柔性梁进行刚-柔耦合动力学仿真.仿真结果表明,与采用模态降阶方法相比,采用Krylov模型降阶方法只需要较低的自由度,就可以得到与采用有限元方法完全一致的结果.说明Krylov模型降阶方法能够有效的用于柔性多体系统的模型降价研究.     

Optimal control of a particular class of singularly perturbednonlinear discrete-time systems

This note studies a class of discrete-time nonlinear systems which depend on a small parameter. Using the singular perturbation theory in a systematic way, we give a trajectory approximation result based on the decomposition of the model into reduced and boundary layer models. This decomposition is used to analyze optimal control via maximum principle of such systems. As a result, significant order reduction of optimal control problems is achieved