计算微波谐振电路时域响应的外推方法
【专利说明】计算微波谐振电路时域响应的外推方法
[0001] 方法领域
[0002] 本发明属于电磁学仿真计算与信号处理的交叉领域，具体是一种电磁特性的快速 仿真与外推计算方法。
[0003] 背景方法
[0004] 在信号处理领域，信号的外推方法可以用于实现信号参数的特征估计、雷达目标 探测识别、核磁共振快速成像等方法；在数值计算仿真领域，外推方法可以极大的加快计 算，有利于时间与内存资源的合理利用。在过去的几十年里，计算电磁学数值分析方法被广 泛运用于各种电磁问题的分析当中，对通讯、集成电路设计、射频器件与天线设计、电磁兼 容、电磁散射与目标识别等领域的发展起了极大地推动作用，不仅有利于降低实验研究的 成本与风险，缩短设计与开发周期，还对理论研究有很大帮助。
[0005] 国内外对于电磁响应信号参数估计和外推方法的研究已经发展多年，从1971年 C.E.Baum提出的奇点展开法，到后来发展起来的Prony算法、KT算法和Matrix-Pencil算 法及其它们的各种改进算法，已经可以从信噪比较低的环境中中提取信号特征，这类方法 的数据主要来自于实际测量，优点是只需要知道信号的一系列采样数据，即可推测出信号 特征，有利于未知源或目标的探测识别，缺点是受采样值影响较大，且外推范围有限；近年 来，也有研究者将基于数值仿真的全波分析方法与信号外推方法结合来加快计算过程，采 样数据来源于数值仿真计算，主要用于已知结构的响应特性快速分析，例如FDTD与改良的 MatrixPencil方法结合，来分析电路散射特性等，缺点主要是耦合度低，由于所采用的外 推方法是一种纯粹的信号处理过程，脱离目标结构本身特性，大多只适用于分析衰减与非 衰减的周期信号与指数型叠加信号，适用范围有一定的局限性，且受外推算法的影响，外推 范围仍然有限。
【发明内容】
[0006] 本发明的目的在于提供一种基于数值分析方法计算微波谐振电路响应的外推方 法，该方法提高了数值分析方法与外推方法的耦合度，结合目标本身结构特性外推计算，使 得外推范围大幅度提高，计算更加快速，外推结果更加可靠。
[0007] 实现本发明目的的方法解决方案为：一种计算微波谐振电路时域响应的外推方 法，步骤如下：
[0008] 第一步，采用全波分析的时域谱元法对使用曲六面体剖分的目标模型进行显式格 式的迭代求解。
[0009] 第二步，用一定间隔时间步(采样间隔需要满足奈奎斯特采样定律）采样计算的电 场构成一个正交矩阵V，用该矩阵对原问题进行降阶，求降阶系统的特征值问题，每一个采 样歩进行上述过程，直到每个采样步计算得到的重要特征值趋于稳定，预处理结束。
[0010] 第三步，用重要特征值对应的特征向量构成一个矩阵，该矩阵对预处理结束时 的降阶系统进行再降阶，即特征提取，可以只保留下原系统的重要特征，构造出用于外推计 算的小规模系统，通过迭代求解该系统，可以还原出原系统的解，从而实现快速而准确的外 推计算。
[0011] 本发明与现有方法相比，其显著优点：
[0012] (1)所采用的数值方法为时域谱元法，具有谱精度，由于产生的质量矩阵为块对角 矩阵，可以快速求逆，而且可以采用显式格式的差分迭代求解，求解速度较常用的时域电磁 学分析方法快，质量矩阵与刚度矩阵都满足对称性与正定半正定特性，使得广义特征值的 求解收敛速度更快。
[0013] (2)采用改进型的模态叠加法进行参数估计，使得提取特征参数时只需要求解小 规模的特征值问题，速度快，且降阶后的矩阵特性变好，参数提取更加准确。
[0014] (3)提取出来的特征参数用于重构一个小规模线性系统用于外推计算，这种方式 提高了外推方法与结构本身特性的耦合度，使得外推的范围大大增加，结果更加可信。
【附图说明】
[0015] 图1实例1卫星舱体模型与曲六面体剖分效果示意图。
[0016] 图2实例1使用外推方法与纯数值方法的时域电场结果对比。
[0017] 图3实例2含不连续结构的矩形波导模型示意图。
[0018] 图4实例2观察点1使用外推方法与纯数值方法的时域电场结果对比。
[0019] 图5实例2观察点2使用外推方法与纯数值方法的时域电场结果对比。
[0020] 图6实例2使用外推方法与纯数值方法的S参数结果对比。
【具体实施方式】
[0021] 下面结合附图对本发明作进一步详细描述。
[0022] 结合图1，本发明为一种计算微波谐振电路时域响应的外推方法，具体实施步骤如 下：
[0023] 第一步，时域谱元法微分方程的建立。
[0024] 对目标模型使用ANSYS软件建模，并使用20个节点的曲六面体剖分，得到模型剖 分后的单元与节点的信息文件，将每个20点表示的曲六面体映射成一个三轴坐标均从-1 到1的标准立方体，这里使用20个节点曲六面体的形函数为：
【主权项】
1. 一种计算微波谐振电路时域响应的外推方法，其特征在于，步骤如下： 第一步，构造一个待求解的时域谱元法线性系统方程； 对物理坐标系下的三维微波谐振类型的结构进行建模并采用曲六面体剖分，得到模型 剖分后的单元与节点信息，将物理坐标系下得到的曲六面体映射为三边坐标均为-1到1的 标准立方体，映射后的坐标系称为标准参量坐标系； 整个模型为求解域Ω，求解域Ω内的电场f满足电场时域波动方程：
(1) 其中ε为求解域内的介电常数，μ为磁导率，σ为电导率，J为电流密度； 对于方程（1)，采用Gauss-Lobatto-Legendre插值基函数展开物理坐标系单元内电场 E ；
(2) 其中，表示物理坐标系中的展开基函数，ei(t)为展开系数，N表示单元内展开基函数 的个数； 然后，对方程（1)使用Galerkin测试，得到：
其中i代表展开基函数索引，j代表测试基函数索引，f表示单元e内，$为单元e上 的残数加权积分，Galerkin测试要求在整个求解域Ω上的残数加权积分为0,即必须满足
最后，整理（3)式得到待求解的时域谱元法线性系统方程：
其中，向量为待求解向量，代表标准参量坐标系下的不同位置和方向的电场在某 时刻的值，J⑴为电流密度，这里作为激励源； (4)式中其他矩阵或向量形式为：
采用中心差分格式对方程（4)迭代求解：
求出标准参量坐标系下的电场^1后，通过雅克比矩阵[J]映射回物理坐标系，得到物 理坐标系下电场：
至此，完成了（4)式所示的时域谱元法线性系统方程的构造； 第二步，进行预处理得到重要特征构成的特征向量矩阵 首先，创建一个空的V矩阵,行为ge+1向量的元素个数,初始列数为O ; 对第一步构造的系统方程进行（5)式所示的时间步差分迭代求解，每一定时间间隔对 进行采样，将采样的g"+l向量追加为V矩阵的一列向量，对V矩阵正交化，V矩阵的列数 为k ; 将气〇用V矩阵展开，= 将其代入（4)式，方程两边同时左乘[V]T,得 到：
其中[VTTV]、[VTSV]、[VtRV]矩阵是规模为kXk的对称矩阵，k << N，N为向量 中元素的个数；方程(6)的阶数被降低到了 k阶，该降阶方程对应的特征多项式方程为：
求解（7)式所示的特征值方程，得到k对特征值λ与特征向量_，根据重要特征值的判 定准则来判断是否停止预处理，预处理停止后，取出能够表征重要特征的r个特征向量多， 构成^矩阵； 第三步，外推计算小规模系统； 用第二步最终构造的矩阵，展开方程（6)的解
(8) 其中k为（6)式中V矩阵的列数，r为重要特征值的个数，Xi为每个特征向量^的展开 系数； 将⑶式代入（6)式，方程两边同时左乘Φ/，
得到： 该（9)式所示的系统规模仅为rXr，r<k<<N; 迭代求解（9)式所示的小规模系统,得到，再根据
还原出系统（4)的解i(l)，从而实现快速准确的外推计算。
2. 根据权利要求1所述的计算微波谐振电路时域响应的外推方法，其特征在于：所述 步骤一中,物理坐标系下得到的曲六面体映射为三边坐标均为-1到1的标准立方体的步骤 如下： 使用20个节点的曲六面体形函数Pi表示一个曲六面体：
其中（ξ，η, ζ)表示标准参量坐标系下标准立方体内任意一点的坐标，（ξ i, Hi, ζ i) 是20个曲六面体节点映射到标准参量坐标系下的坐标； 标准参量坐标系下任意一点的坐标（ξ，n，O对应的物理坐标系下的坐标为 (X，y, ζ)，映射关系如下：
采用雅克比矩阵表示点（ξ，Π ，4)与（x，y，z)之间的映射关系：
3. 根据权利要求1所述的计算微波谐振电路时域响应的外推方法，其特征在于：所述 步骤一种，采用Gauss-Lobatto-Legendre插值基函数的步骤如下： 设为标准参量坐标系下的展开基函数，
通过映射关系，物理坐标系下基函数2)与标准参量坐标系下的基函数 tW. )之间的关系表不为：
其中φ^、Φ1、分别代表标准参量坐标系下三个方向ξ、η、ξ的基函数，r、s、t 则代表标准参量坐标系下某点的坐标，IJl为雅克比矩阵；Φ^ξ)、Φ3(Π )、是具 有正交性质的Gauss-Lobatto-Legendre基函数，表达式为：
其中~为Legendre多项式，4: (Ir)为Legendre多项式的一阶导数，Νξ代表基 函数的插值阶数，的微分形式如下：
4.根据权利要求1所述的计算微波谐振电路时域响应的外推方法，其特征在于：所述 步骤二中的判定准则如下： 对于权利要求1中第二步的预处理过程，每次采样构成正交化的V矩阵，求解（7)式得 到k个特征值及每个特征值对应的特征向量； 设置一个条件：
其中λ i代表第i个特征值，i = 1，2. . . k，η代表前一次采样，n+1代表当前采样，ε i 是设定的一个阈值，通常取〇. 001或者〇. 01，该条件称之为稳定性条件，满足该条件的特征 值为稳定的重要特征值； 再设定一个重要特征值个数M，只要满足稳定性条件的特征值个数r大于等于M，则停 止预处理过程，此时r个特征值对应的特征向量构成矩阵Φρ
【专利摘要】本发明公开了一种计算微波谐振电路时域响应的外推方法。该方法的实现分为两个过程。第一个过程为预处理过程，对三维微波谐振电路结构建模并使用曲六面体剖分，使用时域谱元法数值分析方法计算出前一段时间的响应，并通过改进的模态叠加方法提取出能够准确描述线性系统的一系列重要特征值，这些特征值对应的特征向量用以构成特征向量矩阵；第二个过程即外推计算过程，利用预处理过程最终得到的特征向量矩阵，对时域谱元法描述的线性系统进行特征提取降阶，重新形成一个只包含重要特征的小规模线性系统，该系统的时域外推求解不仅快速而且准确，能够适用于各种具有谐振、损耗与反射结构电磁响应的外推分析。
【IPC分类】G06F17-50
【公开号】CN104679919
【申请号】CN201310617627
【发明人】丁大志, 陈如山, 樊振宏, 牛荣鑫, 曾强, 盛亦军
【申请人】南京理工大学
【公开日】2015年6月3日
【申请日】2013年11月27日