一种线性系统中基于误差控制的矢量拟合模型降阶方法.pdf

(19)中华人民共和国国家知识产权局(12)发明专利申请(10)申请公布号CN106372348A(43)申请公布日2017.02.01(21)申请号201610808296.9(22)申请日2016.09.07(71)申请人北京航空航天大学地址100191北京市海淀区学院路37号(72)发明人谢树果张卫东苏东林阎照文(74)专利代理机构北京科迪生专利代理有限责任公司11251代理人成金玉卢纪(51)Int.Cl.G06F17/50(2006.01)权利要求书2页说明书6页附图4页(54)发明名称一种线性系统中基于误差控制的矢量拟合模型降阶方法(57)摘要本发明公开了一种线性系统中基于误差控制的矢量拟合模型降阶方法，针对传统矢量拟合在宽频带下，尤其是具有谐振点的高动态线性系统，所建极点-留数模型阶数往往很高。本发明引入误差控制函数通过不断迭代的方式选取出包含线性系统重要特性的数据，通过忽略某些代表细节处的微小变化数据，从而在牺牲一定建模精度情况下有效地降低所建模型的阶数。在得到满足建模精度要求的情况下，获得低阶极点-留数宏模型。然后将极点-留数宏模型进一步转化为状态空间模型，然后再进行无源性处理，最后输出所建极点-留数宏模型。本发明实现了模型精度和复杂度的平衡，提高电磁兼容电磁发射的仿真效率。CN106372348ACN106372348A权利要求书1/2页1.一种线性系统中基于误差控制的矢量拟合模型降阶方法，其特征在于包括如下步骤：步骤一：通过仿真获取设定带宽下线性系统端口网络参数；步骤二：在整个设定宽带内从步骤一种获取的网络参数中均匀地选取几对数据作为初始建模数组；步骤三：采用矢量拟合算法对建模数组进行建模；步骤四：计算采用矢量拟合算法所建模型与建模数组之间的均方根误差σ，若均方根误差σ≤σ0，则输出矢量拟合模型，σ0为预先设置的均方根误差；否则，增加阶数，重复步骤三，直到满足σ≤σ0；步骤五：在得到矢量拟合模型后，计算误差函数δ(si)；步骤六：若max(δ(si))≤δ0则输出宽频带的极点-留数模型，进入步骤六中的无源性判定及增强，δ0为预先设定的值；若max(δ(si))＞δ0，则增加新的数据到建模数组，重新进入步骤三迭代循环；经过迭代循环，预先设定的值σ0和δ0都被满足，则输出的极点-留数宏模型为：其中，s＝jω表示复频率变量；ω为角频率，am为极点留数宏模型HFit(s)第m个极点值，rm为极点留数宏模型HFit(s)第m个留数值，d为极点留数宏模型HFit(s)的常数项值，N表示极点留数宏模型HFit(s)的阶数；步骤六：无源性判定及增强，经过多次迭代修正后最终保证模型满足无源性要求；最后输出高精度低阶极点-留数宏模型HFit：其中，Δrm为极点留数宏模型HFit(s)第m个留数rm的微扰动值，Δd为极点留数宏模型HFit(s)常数项d的微扰动值。2.根据权利要求1所述的线性系统中基于误差控制的矢量拟合模型降阶方法，其特征在于：所述步骤一中，采用全波仿真软件HFSS、CST或者MicrowaveOffice对进行仿真。3.根据权利要求1所述的线性系统中基于误差控制的矢量拟合模型降阶方法，其特征在于：所述步骤五中，误差函数δ(si)的计算公式如下：其中HRaw(si)代表仿真得到的线性系统端口网络参数数据，HFit(si)代表矢量拟合所建宏模型,i＝1,2,...K。4.根据权利要求1所述的线性系统中基于误差控制的矢量拟合模型降阶方法，其特征在于：所述步骤五中，所增加的新数据选取在δ(si)最大处。5.根据权利要求1所述的线性系统中基于误差控制的矢量拟合模型降阶方法，其特征在于：所述步骤六具体实现为：将得到极点-留数宏模型转化为状态空间模型：2CN106372348A权利要求书2/2页其中，x为状态变量，u为输入向量，y为输出向量，表示对x求导，A、B、C、和D为状态方程系数矩阵；由状态方程系数矩阵构建无源性判定矩阵J＝(A-B(D-I)-1C)(A-B(D+I)-1C)；计算矩阵J的特征值λi，若λi为负实数，则状态空间模型不满足无源性；反之，状态空间模型满足无源性；当状态空间模型不满足无源性时，通过对留数rm和常数项d进行扰动处理增强其无源性；使得无源性判定矩阵J的特征值λi不存在负实数；经过多次迭代修正后最终保证模型满足无源性要求；最后输出高精度低阶极点-留数宏模型HFit。3CN106372348A说明书1/6页一种线性系统中基于误差控制的矢量拟合模型降阶方法技术领域[0001]本发明涉及一种线性系统中基于误差控制的矢量拟合模型降阶方法，属于电磁兼容技术领域，为实现电磁兼容量化设计仿真奠定基础。背景技术[0002]随着电子系统向着高频率、高速化、小型化、多功能化和集成化的方向发展，电磁兼容性问题在电子系统