基于时域多步积分的互连线模型降阶方法

基于时域多步积分的互连线模型降阶方法
【专利摘要】本发明属于集成电路领域，涉及一种基于时域多步积分的互连线模型降阶方法；该方法包括步骤：首先读取互连线电路的特性数据并利用改进节点电压法建立对应的时域方程；然后对原始电路的时域方程进行多步积分得到关于状态变量的二阶递推关系；进而通过二次Arnoldi方法得到投影矩阵，再通过投影矩阵对原始时域方程进行投影得到降阶系统，最后用多步积分法对降阶系统进行离散求得时域输出。本方法可保证时域积分后降阶系统和原始系统的状态变量在离散时间点的匹配以及时域降阶精度，同时保证降阶过程的数值稳定性及降阶系统的无源性。本发明的方法复杂度低、精度高。
【专利说明】基于时域多步积分的互连线模型降阶方法
【技术领域】
[0001]本发明属于集成电路领域，涉及一种可对互连线电路进行快速仿真的模型降阶方法，具体涉及一种基于时域多步积分的互连线模型降阶方法。
【背景技术】
[0002]随着集成电路设计和制造技术的不断进步，互连线已经成为影响集成电路的工作性能和可靠性的关键因素；尤其是日益增加的信号频率和电路规模对互连电路建模和分析带来了很大挑战。当前，模型降阶技术已经成为大规模互连电路分析的主流方法。
[0003]传统的模型降阶方法可分为频域模型降阶和时域模型降阶两种。在频域模型降阶方法中，最经典的模型降阶方法有AWE方法，Kry I ο V子空间类方法，例如PVL方法、PRIMA方法。然而，频域模型降阶方法无法保证降阶后的电路在时域的精度；频域逼近的误差转换到时域会放大，频域很小的误差在时域可能会产生很大的误差。为了解决上述问题，近年来直接在时域进行降阶的时域模型降阶方法不断被提出。有研究提出采用基于Chebyshev多项式的时域降阶方法(文献4)，以及提出采用基于小波配置的时域降阶方法(文献5)，但时域模型降阶方法在求解正交多项式展开系统复杂度过高，难以用于大规模系统的模型降阶。
[0004]文献6提出了一种基于时域梯形法差分模型降阶算法，该算法结合了现有的频域模型降阶方法和时域模型降阶方法的优势，与现有的时域模型降阶方法相比，该算法计算复杂度极大降低，与现有的频域模型降阶方法相比，该算法在时域有更高的精度；但是，该算法只对单一的输入有效。
[0005]文献7提出了一种基于时域单步积分的模型降阶方法，该算法克服了基于时域梯形法差分模型降阶方法只对单一输入有效的缺点，对各种输入都有效，同时具有比现有的时域模型降阶方法计算复杂度低和比现有的频域模型降阶方法精度高的优点。
[0006]目前，需要一种基于时域单步积分的模型降阶方法，该方法仅采用单步法来对积分进行离散，其降阶精度和效率可进一步提高。
[0007]与本发明有关的参考文献有:
[0008][I]L.T.Pillage and R.A.Rohrer, “Asymptotic Waveform Evaluation forTiming Analysis”，IEEE Trans.Computer-Aided Design, vol.9，pp.352 - 366，Apr.1990.[0009][2]P.Feldmann and R.W.Freund, “Efficient Linear Circuit Analysis byPadevia Lanczos process，，，IEEE Trans.Computer-Aided Design, vol.14，pp.639 -649，Mayl995.[0010][3]Odabasiogluj M.Celik and L.Pileggij “PRIMA:Passive Reduced-OrderInterconnect Macromodeling Algorithm”，IEEE Trans.0n CAD of Integrated Circuitsand Systems, vol.17，n0.8，pp.645 - 654，Aug.1998.[0011][4] Janet Meiling Wang, Chia-Chi Chu，Qingjian Yu and Ernest S.Kuh，”0nProjection-based Algorithms for Mode 1-order-reduction of Interconnects，，，IEEEtrans.Circuits and Systems, vol.49，n0.11，pp.1563—1585，Nov.2002.[0012][5] Xuan Zeng, Lihong Feng, Yangfeng Su, Wei Cai, Dian Zhou andCharles Chiang, ，，Time Domain Model Order Reduction by Wavelet collocationmethod，，，pp.1-6，March6, IEEE/ACM Design Automation and Test in Europe, 2006.[0013][6] Hou Limin, Yang fan, Zeng Xuan.“An Efficient T ime-domainTrapezoidal Difference Based Model Order Reduction Method forInterconnect Circuits，，，Journal of Computer-Aided Design & ComputerGraphics, vol.24, pp.683-689，May 2012
[0014][7]侯丽敏，”互连线高效时域模型降阶算法研究”，硕士论文，复旦大学，2012.。
【发明内容】

[0015]本发明的目的是克服现有技术的缺陷和不足，提供一种基于时域多步积分的互连线模型降阶方法。本方法可保证时域积分后降阶系统和原始系统的状态变量在离散时间点的匹配，保证时域降阶精度，同时也保证了降阶过程的数值稳定性及降阶系统的无源性。
[0016]具体而言，本发明的基于时域多步积分的互连线模型降阶方法，其特征在于，其步骤如下(如图1所示):
[0017]步骤一:读取互连线电路的特性数据及输入激励；
[0018]步骤二:利用改进节点分析方法(MNA)建立互连线电路的时域方程；
[0019]步骤三:用多步积分方法对互连线电路的时域方程进行离散，得到二次Arnoldi递推关系，该递推关系形成了一个二次Krylov子空间；
[0020]步骤四:利用二次Arnoldi算法构造步骤三产生的递推关系的正交投影矩阵Vq e RNXn, n=N,其中η为降阶系统的阶数，N为原始系统的阶数；
[0021]步骤五:利用正交投影矩阵\ e RNXn，对互连线电路的时域方程进行合同变换获得η阶的降阶系统；
[0022]步骤六:利用多步积分法数值求解降阶系统的时域输出。
[0023]本发明中，采用基于多步积分的方法对时域方程进行离散(如图2所示)，基于单步积分的方法仅利用当前时刻信息逼近积分(如图3所示)，基于多步积分的方法利用多个数据点构造出一条曲线，用该条曲线张成的面积来逼近原函数的积分，相比于基于单步积分的方法更能精确地表示被积函数的积分值，其精度高于现有技术(如图2所示)中的矩形逼近。
[0024]本发明中，采用二次Arnoldi算法产生正交投影矩阵，再通过投影矩阵对原始时域矩阵进行投影得到降阶系统；能保证时域积分后降阶系统和原始系统的状态变量在离散时间点的匹配，及时域降阶精度。
[0025]本发明中，直接在时域上进行降阶，可消除时频转换引入的误差。
[0026]本发明中，采用二次Arnoldi算法构造二次Krylov子空间的正交基,二次Arnoldi方法其数值稳定性良好。
[0027]本发明中，基于合同变换获得降阶系统，可保证降阶系统的无源性。
[0028]本发明中，直接利用二次Arnoldi算法求解投影矩阵，算法复杂度极大降低。
[0029]本发明所述的方法与现有技术相比，具有如下优点:
[0030](I)高的降阶精度[0031]首先，本发明采用基于多步积分的方法对时域方程进行离散(如图2所示)，基于单步积分的方法仅利用当前时刻信息逼近积分(如图3所示)，基于多步积分的方法利用多个数据点构造出一条曲线，用该条曲线张成的面积来逼近原函数的积分，相比于基于单步积分的方法更能精确地表示被积函数的积分值，其精度高于图2中的矩形逼近，因此，相比于单步积分的方法，本发明采用的基于多步积分模型降阶方法要有更高的精度；
[0032]其次，本发明采用二次Arnoldi算法产生正交投影矩阵，再通过投影矩阵对原始时域矩阵进行投影得到降阶系统；该方法可保证时域积分后降阶系统和原始系统的状态变量在离散时间点的匹配，保证时域降阶精度，因此，本发明具有较高的降阶精度；
[0033]此外，本发明由于直接在时域上进行降阶，可消除时频转换引入的误差；相比于频域降阶方法，本发明在时域有更高的精度；
[0034](2)良好的数值稳定性
[0035]本发明采用二次Arnoldi算法来构造二次Krylov子空间的正交基；二次Arnoldi方法是数值稳定的，因此，本基于时域多步积分的模型降阶方法具有良好的数值稳定性。
[0036](3)保证无源性
[0037]本发明提出的基于时域多步积分的模型降阶方法，基于合同变换来获得降阶系统；经过合同变换得到降阶系统可保证降阶系统的无源性，因此，本发明的基于时域多步积分的模型降阶方法得到的降阶系统可以保持原系统的无源性；
[0038](4)低的计算复杂度
[0039]现有时域模型降阶方法在求解正交多项式展开系数时复杂度过高，本发明直接利用二次Arnoldi算法求解投影矩阵，算法复杂度极大降低。
[0040]本发明方法能保证时域积分后降阶系统和原始系统的状态变量在离散时间点的匹配以及时域降阶精度，和降阶过程的数值稳定性及降阶系统的无源性。本发明比现有的时域模型降阶方法复杂度低和比现有的频域模型降阶方法精度高，尤其与时域单步积分的模型降阶方法相比，可在保证与之计算复杂度相当的基础上，达到更高的精度。
【专利附图】

【附图说明】
[0041]图1是本发明基于时域多步积分的互连线模型降阶方法的流程图。
[0042]图2是基于单步积分方法的示意图。
[0043]图3是基于多步积分方法的不意图。
[0044]图4是阶数为3298，输入IGHz的脉冲信号的总线电路降阶到15阶时，本发明基于时域多步积分的模型降阶方法、基于小波配置的时域模型降阶方法、频域降阶方法PRIMA和基于时域单步积分的模型降阶方法的误差比较图。
【具体实施方式】
[0045]为使本发明的上述目的、特征和优点能够更加明显易理解，下面通过具体的实例进一步说明本发明。
[0046]实施例1
[0047]本发明的基于时域多步积分的互连线模型降阶方法的包括步骤(如图1所示):
[0048]步骤一:读取互连线电路的特性数据及输入激励；互连线电路的特性数据包括经过互连线寄生参数提取得到的电阻、电容和电感寄生网表；
[0049]步骤二:利用改进节点电压法建立互连线电路的时域方程(I):
[0050]
【权利要求】
1.一种高效的基于时域多步积分的互连线模型降阶方法，其特征在于:步骤如下: 步骤一:读取互连线电路的特性数据及输入激励； 步骤二:利用改进节点分析方法建立互连线电路的时域方程； 步骤三:用多步积分方法对互连线电路的时域方程进行离散，得到二次Arnoldi递推关系，该递推关系形成了一个二次Krylov子空间； 步骤四:利用二次Arnoldi算法构造步骤三产生的递推关系的正交投影矩阵V, e RNXn，n=N,其中η为降阶系统的阶数，N为原始系统的阶数； 步骤五:利用正交投影矩阵V, e RNXn，对互连线电路的时域方程进行合同变换获得η阶的降阶系统； 步骤六:利用多步积分法数值求解降阶系统的时域输出。
2.如权利要求1所述的方法，其特征在于:所述步骤二中，时域方程为:

3.如权利要求1所述的方法，其特征在于:所述步骤三中，用多步积分方法对互连线的时域方程进行离散，得到非齐次递推关系如下:
4.如权利要求1所述的方法，其特征在于:所述步骤三中，形成的递推关系形成二次Krylov子空间,利用二次Arnoldi算法求得二次Krylov子空间的正交基,然后通过正交基对原始系统进行投影得到降阶系统。
5.如权利要求1所述的方法，其特征在于:所述步骤五中利用正交投影矩阵V,得到的降阶系统为:
0m + Cx(f) = Bu(t) (8)
y(i) = llx(i) 其中= € = QlCQtt M = QlB I = QjL ?
6.如权利要求1所述的方法，其特征在于:所述步骤六中，多步积分法离散(8)，得:
【文档编号】G06F17/50GK103678738SQ201210332712
【公开日】2014年3月26日 申请日期:2012年9月9日 优先权日:2012年9月9日
【发明者】曾璇, 郭倞, 杨帆 申请人:复旦大学