[发明专利]一种基于最优模型降阶阈值的架空线时域仿真方法及系统

权利要求书

1.一种基于最优模型降阶阈值的架空线时域仿真方法，其特征在于，所述架空线时域仿真方法包括：

利用矢量拟合方法，获取架空线频域传播函数矩阵中每个元素对应的初始留数极点比值集合；所述初始留数极点比值集合中包括多个初始留数极点比值；

获取留数极点比值阈值；

选取所述初始留数极点比值集合中的第一个大于所述留数极点比值阈值的初始留数极点比值对应的元素作为目标元素；

在矢量拟合过程中获取所述目标元素对应的系数矩阵，并对所述系数矩阵进行奇异值分解，得到所述目标元素对应的第一奇异值矩阵；

将所述第一奇异值矩阵中的奇异值按照数值大小排序，并对排序后的奇异值根据顺序进行编号，获得第二奇异值矩阵；其中，最大奇异值的编号为1，最小奇异值的编号为N_σ；N_σ表示所述第一奇异值矩阵中奇异值的个数；

根据所述第二奇异值矩阵采用黄金分割法确定第一阈值编号，具体包括：

根据所述第二奇异值矩阵，采用公式n_threshold＝I_L+α(I_R+I_L)+1确定第一阈值编号；其中，I_L＝1，I_R＝N_σ，α＝0.618，n_threshold表示所述第一阈值编号，n_threshold向下取整；

根据所述第二奇异值矩阵采用黄金分割法确定第二阈值编号，具体包括：

根据所述第二奇异值矩阵，采用公式n′_threshold＝I_L+α(I_R+I_L)确定第二阈值编号；其中，I_L＝1，I_R＝N_σ，α＝0.618，n′_threshold表示所述第二阈值编号，n′_threshold向下取整；

判断所述第一阈值编号是否等于所述第二阈值编号，获得第一判断结果；

若所述第一判断结果为所述第一阈值编号等于所述第二阈值编号，确定所述第一阈值编号对应的奇异值为最优模型降阶阈值；

根据所述最优模型降阶阈值对应的降阶模型进行架空线时域仿真；

若所述第一判断结果为所述第一阈值编号不等于所述第二阈值编号，将所述第二奇异值矩阵中编号大于所述第一阈值编号的奇异值置零，得到第三奇异值矩阵；

根据所述第三奇异值矩阵确定所述目标元素对应的修改系数矩阵；

将所述第二阈值编号的值赋值给所述第一阈值编号；

根据所述修改系数矩阵计算所述架空线频域传播函数矩阵中每个元素对应的修改留数极点比值集合；所述修改留数极点比值集合中包括多个修改留数极点比值；

根据所述修改留数极点比值集合更新所述第二阈值编号，返回所述判断所述第一阈值编号是否等于所述第二阈值编号的步骤。

2.根据权利要求1所述的架空线时域仿真方法，其特征在于，所述在矢量拟合过程中获取所述目标元素对应的系数矩阵，并对所述系数矩阵进行奇异值分解，得到所述目标元素对应的第一奇异值矩阵，具体包括：

在矢量拟合过程中，利用所述目标元素对应的极点确定所述目标元素对应的系数矩阵；

根据所述系数矩阵，采用公式A＝USV^T进行奇异值分解，确定所述目标元素对应的第一奇异值矩阵；其中，A表示所述系数矩阵，S表示所述第一奇异值矩阵，U表示由奇异值所对应的特征向量组成的左奇异向量矩阵，V表示由奇异值所对应的特征向量组成的右奇异向量矩阵，V^T表示矩阵V的转置。