基于二重Arnoldi过程的分段二次化算法在非线性电路系统模型降阶中的应用研究.pdfVIP

主星皇王兰垒皇堕兰．墨竺堂垒苎±垄旦堡垒丝兰叁 I d 基于二重Arno i过程的分段二次化算法在非线性电路 系统模型降阶中的应用 饶建翔 蒋耀林 71 (西安交通大学理学院信息与系统研究所 西安 0049) 摘要： 本文提出了一种新的模型降阶算法．这种算法利用分段二次化的方法对非线性电路系统进行r近 似，同时构造T-重的Krylov子空间并利用Arnoldi过程来对近似系统进行降阶．我们还给出了相关定理以 及算法的性能分析．数值算例表明，经过降阶后的较小的电路模型能准确反映原系统的相关瞬态情况，相比 于类似算法，其精确性得到了提高 关键词： 模型降阶二重Arnoldi过程Krylov子空间Taylor级数 一、引言 模型降阶这一思想，最初出现在E．J．Davison的论文中，其发展已经有30多年历史．到目 前为止，众多的模型降阶算法基本上是关于线性系统模型降阶的，更确切的说，大多数是适 用于LTI的模型降阶．一般来说，模型降阶算法可以大致分为两类：一类就是基于Krylov子 空间技术的方法；另外一类是基于奇异值分解的平衡性方法，前者瞬态分析时近似效果较好， 运算相对简单：而后者的全局精确性较好，误差容易估计，但是一个缺点就是运算量较大． 非线性系统由于其结构复杂，基于Krylov子空间的算法都不能直接应用于系统．在很长 一段时间里对于非线性系统的模型降阶算法发展缓慢． 本文算法是基于非线性函数多点Taylor展开的分段二次化方法的，利用二重Amoldi过程 构造出正交矩阵对近似模型降阶，最终我们通过一种加权的办法得到原电路系统的降阶模型． 全文余下的结构如下：第二部分是对算法的详细的介绍，在论文的第三部分是算法相关定理 以及性能分析。第四部分的数值算例对算法进行了验证． 二、算法介绍 我们考虑这样一个常见的大型非线性电路系统： l C掣+，(x(f))=bu(O， (1) ‘ dr 【，；。r，(f)． 这里的C∈R“”，为常数矩阵，为了方便研究问题，本文假设系统为SISO的，而且c是非 奇异的．x∽∈R”为电路中结点电压参数，f【《0)是电路中非线性元件作用F的函数，函数 u(O是由电压源或电流源所控制的激励源． 为了方便研究问题起见，本文假设f(x)是任意阶可导的，对于满足以上条件的非线性函数 “x1，我们记给定的初始状态为Xo．一般情况下，在这些点对f(x)的Taylor展开为： (2) ，(z)=f(xo)+以(x—Xo)+÷旷(扛一Xo)o扛一矗))+一 VI一401 ——————————兰兰坐蔓王兰生垒塑兰堡竺兰全笙±!!旦堡叁堕苎壅 @盘 阵．从而我们得到了近似的二次化系统(见文[1】)： (3) IY=c’x(f)． 将系统中的微分方程分解为两个线性的子微分方程 c型兰：：』堕+4。。。。。=6“(，)+Aoxo-，(z。) (4) z-xo)。(x。：，一x。)’ (5) f!!i：；；!堕+4。x。r：=一W((x。E 中微分方程的解是足够近似的，具体证明在下节给出 对于。‘一一，,r—l 4)，我们构造这样一个Krylov子空间