[发明专利]一种用于集成电路的快速汉克尔变换方法及装置

权利要求书

1.一种用于集成电路的快速汉克尔变换方法，其特征在于，包括如下步骤：

基于计算集成电路任意层的覆铜区域上任意点电流源在其他层产生的场的并矢格林函数，获取所述并矢格林函数包含的所有贝塞尔积分，由此确定并矢格林函数中所用的贝塞尔函数的阶数；所述贝塞尔积分为包含贝塞尔函数的被积函数的无穷积分；

根据积分核函数与贝塞尔函数乘积形成的被积函数，对贝塞尔积分进行汉克尔变换；

对所述汉克尔变换进行离散化，得到无限长度滤波器的表达式；

根据所述积分核函数与贝塞尔函数乘积形成的被积函数的衰减速度，在预设的精度控制范围内对无限长度滤波器进行截断，得到有限长度滤波器的表达式，根据所述有限长度滤波器的表达式，得到有限长度滤波器的方程组；

根据所述并矢格林函数中所要计算积分的贝塞尔函数的类别和阶数，找到能给出积分的解析表达式的贝塞尔积分，构造出指定的贝塞尔积分的类别和阶数的贝塞尔积分-汉克尔变换对；

基于构造的贝塞尔积分-汉克尔变换对，代入所述有限长度滤波器的表达式，获得滤波器方程组对应的系数矩阵方程，并对所述滤波器系数矩阵方程进行改进；

对改进后的滤波器系数矩阵方程进行求解，得到滤波器系数；

根据所述滤波器系数、滤波器的采样点以及所述积分核函数，计算贝塞尔积分；

所述根据所述滤波器系数、滤波器的采样点以及所述积分核函数，计算贝塞尔积分，包括如下步骤：

根据采样点公式，输出有限长度滤波器的采样点；

将所述滤波器系数和采样点，进行映射变换，将原区间映射到端点值为1和有限长度的最大值的区间；

对所述积分核函数也进行相同的映射变换，将积分核函数的变量由原区间映射到端点值为1和有限长度的最大值的区间；

根据映射后求得的滤波器系数和采样点，计算所述积分核函数和贝塞尔函数形成的贝塞尔积分。

2.根据权利要求1所述的用于集成电路的快速汉克尔变换方法，其特征在于，所述积分核函数是由集成电路各层材料的电磁参数、各层厚度、集成电路的工作频率、以及源点到场点之间的距离确定的，所述源点到场点之间的距离为集成电路覆铜区域上连续电流离散的点电流源形成的源点到作用于其他层的场点之间的距离。

3.根据权利要求1所述的用于集成电路的快速汉克尔变换方法，其特征在于，所述对所述滤波器系数矩阵方程进行改进包括如下步骤：

根据所述有限长度滤波器的方程组，建立滤波器系数矩阵方程；

将所述滤波器系数矩阵方程两边同加一个列向量δh，h为线性滤波器系数的向量，δ为预先设定的阈值；

δ取一个足够小的值，使得列向量的模远小于输出向量的模；

忽略所述滤波器系数矩阵方程中输出端的列向量δh，得到改进后的滤波器系数矩阵方程。

4.根据权利要求3所述的用于集成电路的快速汉克尔变换方法，其特征在于，所述足够小的值为当前使用的数据类型下的机器精度。