[发明专利]一种抑制架空线频域传播矩阵矢量拟合过程中留数极点比的方法及系统

权利要求书

1.一种抑制架空线频域传播矩阵矢量拟合过程中留数极点比的方法，其特征在于，所述方法包括：

步骤1，利用矢量拟合方法获取频域传播矩阵中每个元素对应的至少一个初始留数极点比值；其中，每个元素对应的极点个数相同；

步骤2，对于每个元素，当存在对应的初始留数极点比值大于第一预设阈值时，对矢量拟合过程中该元素对应的系数矩阵进行奇异值分解，以获取该元素对应的第一奇异值矩阵；

步骤3，对于每个元素，计算该元素对应的第一奇异值矩阵中最大的奇异值与预设精度阈值的乘积，并将该元素对应的第一奇异值矩阵中的小于所述乘积的奇异值置零，以获取该元素对应的第二奇异值矩阵；

步骤4，分别利用每个元素对应的第二奇异值矩阵确定该元素对应的更新的系数矩阵；

步骤5，分别根据每个元素对应的更新的系数矩阵，利用最小二乘法确定该元素对应的留数向量；

步骤6，分别去除每个元素对应的留数值为0的留数对应的极点，并分别计算剩余的极点中每个极点对应的留数与该极点的比值；

步骤7，判断所述剩余的极点中每个极点对应的留数与该极点的比值是否均小于第一预设阈值；若是，则表示抑制成功，能够利用所述剩余的极点和对应的留数确定以有理式方式表示每个元素的频域传播矩阵，以保证时域仿真的稳定性；反之，则表示抑制失败，并返回步骤2。

2.根据权利要求1所述的方法，其特征在于，所述利用矢量拟合方法获取频域传播矩阵中每个元素对应的至少一个初始留数极点比值，包括：

通过矢量拟合方法确定频域传播矩阵中用于分别表示每个元素的有理分式，并分别根据每个元素对应的有理分式确定该元素对应的至少一个初始留数极点比值。

3.根据权利要求1所述的方法，其特征在于，所述第一预设阈值为100；所述预设精度阈值为10^-5。

4.根据权利要求1所述的方法，其特征在于，利用如下方式获取第一奇异值矩阵：

A＝USV^T，

其中，S为第一奇异值矩阵；A为系数矩阵；A矩阵为m×n大小，U矩阵为m×n大小，S矩阵为n×n大小，V矩阵为n×n大小。

5.根据权利要求1所述的方法，其特征在于，所述分别根据每个元素对应的更新的系数矩阵，利用最小二乘法确定该元素对应的留数向量，包括：

对于每个元素，利用该元素对应的更新的系数矩阵构建方程组A′x＝b，并利用最小二乘法求解，以确定该元素对应的留数向量；

其中，A′为某个元素对应的更新的系数矩阵，x为该元素对应的留数向量，x中的数据即为留数；b为常数项，b中的数据即为该元素在某一频率采样点的响应。

6.一种抑制架空线频域传播矩阵矢量拟合过程中留数极点比的系统，其特征在于，所述系统包括：

初始留数极点比值获取单元，用于利用矢量拟合方法获取频域传播矩阵中每个元素对应的至少一个初始留数极点比值；其中，每个元素对应的极点个数相同；

第一奇异值矩阵获取单元，用于对于每个元素，当存在对应的初始留数极点比值大于第一预设阈值时，对矢量拟合过程中该元素对应的系数矩阵进行奇异值分解，以获取该元素对应的第一奇异值矩阵；

第二奇异值矩阵获取单元，用于对于每个元素，计算该元素对应的第一奇异值矩阵中最大的奇异值与预设精度阈值的乘积，并将该元素对应的第一奇异值矩阵中的小于所述乘积的奇异值置零，以获取该元素对应的第二奇异值矩阵；

更新的系数矩阵确定单元，用于分别利用每个元素对应的第二奇异值矩阵确定该元素对应的更新的系数矩阵；

留数向量确定单元，用于分别根据每个元素对应的更新的系数矩阵，利用最小二乘法确定该元素对应的留数向量；

留数极点比值计算单元，用于分别去除每个元素对应的留数值为0的留数对应的极点，并分别计算剩余的极点中每个极点对应的留数与该极点的比值；

判断单元，用于判断所述剩余的极点中每个极点对应的留数与该极点的比值是否均小于第一预设阈值；若是，则表示抑制成功，能够利用所述剩余的极点和对应的留数确定以有理式方式表示每个元素的频域传播矩阵，以保证时域仿真的稳定性；反之，则表示抑制失败，并返回第一奇异值矩阵获取单元。