Abstract:
A method for efficiently estimating crosstalk noise of nanometer VLSI interconnects is provided. In the invention, nanometer VLSI interconnects are modeled as nonuniform distributed RLC coupled trees. The efficiency and the accuracy of moment computation of distributed lines can be shown that outperform those of lumped ones. The inductive crosstalk noise waveform can be accurately estimated in an efficient manner using the linear time moment computation technique in conjunction with the projection-based order reduction method. Recursive formulas of moment computations for coupled RC trees are derived with considering both self inductances and mutual inductances. Also, analytical formulas of voltage moments at each node will be derived explicitly. These formulas can be efficiently implemented for crosstalk estimations.

Description:
FIELD OF THE INVENTION  
       [0001]     The present invention relates generally to a technique that use moment computation to estimate crosstalk noise of nanometer VLSI interconnects, and more particularly to the moment computation of nonuniform distributed RLC coupled trees and the projection-based order reduction method.  
       BACKGROUND OF THE INVENTION  
       [0002]     With the improvement of semiconductor producing technology, interconnect models get more and more important in the design flow of VLSI. While due to quality consideration, increasing timing frequency, decreasing ascending time, increasing online density and use of low-resistance materials, etc. in circuit design, designers should take into consideration of inductance effect at time of constructing interconnects. Also, due to the creation of nanotechnology in recent years, the importance of mutual inductance increases gradually. In the condition that RC is merely taken into consideration, the estimation of crosstalk noise is not accurate (refer to Interconnect Analysis and Synthesis written by C. K. Chang, J. Lillis, S. Lin and N. H. Chang, published by John Wiley and Sons Inc. in 2000). More and longer parallel nets will multiply capacitance coupling and increase current change on victim nets. Inductance added onto aggressor nets will cause more overshooting voltage, and further exacerbate noise on victim nets. The above two phenomena will lead to error in traditional estimation of crosstalk noise. According to the actual condition, inductance is necessary to be introduced to interconnect models to construct RLC coupled trees.  
         [0003]     In the prior technology of estimating crosstalk noise, simulation is generally carried out in circuit. However, though the result of simulation carried in VLSI interconnects is considerably accurate, the computation complexity increases correspondingly. In order to solve this problem, the so-called Model-Order Reduction is gradually adopted in prior technology. (Provided by L. T. Pillage and R. A. Rohrer in “Asymptotic Waveform Evaluation for Timing Analysis,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 9, No. 4, pp. 352-366 in 1990, P. Feldmann and R. W. Freund in “Efficient Linear Circuit Analysis by Pade Approximation Via the Lanczos Process,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 14, No. 5, pp. 639-649 in 1995, and A. Odabasioglu, M. Celik and L. T. Pileggi in “PRIMA: Passive Reduced-order Interconnect Macromodeling Algorithm,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 17, No. 8, pp. 645-653 in 1998). This invention adopts order reduction method to estimate crosstalk noise. However, though the order reduction method can decrease the complexity of noise estimating computation, the computation is still too complicated in the process of noise optimization.  
         [0004]     Among different prior technologies of model order reduction, most of them adopt moment matching method in crosstalk noise estimation in interconnects. In consideration of computation efficiency, many traditional methods to estimate crosstalk noise are developed in RLC coupled trees. The traditional technology includes One-Pole Model, 1P (A. Vittal, L. H. Cheng, M. Marek-Sadowska, K. P. Wang and S. Yang, “Crosstalk in VLSI Interconnects,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 18, pp. 1817-1824 (1999), and A. Vittal and M. Marek-Sadowska, “Crosstalk Recuction for VLSI,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 16, pp. 290-298 (1997), Modified One-Pole Model, M1P (Q. Yu and E. S. Kuh, “Moment Computation of Lumped and Distributed Coupled RC Trees with Application to Delay and Corsstalk Estimation,” Proceedings of the IEEE, Vol. 89, No. 5, pp. 772-788 (2001); Two-Pole Model, 2P (M. Kuhlmann and S. S. Sapatnekar, “Exact and Efficient Corsstalk Estimation,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 20, No. 7, pp. 858-866 (2001), and Q. Yu and E. S. Kuh, “Moment Computation of Lumped and Distributed Coupled RC Trees with Application to Delay and Crosstalk Estimation,” Proceedings of the IEEE, Vol. 89, No. 5, pp. 772-788 (2001), Stable Two-Pole Model, S2P (E. Acar, A. Odabasioglu, M. Celik and L. T. Pileggi, “S2P: A Stable 2-Pole RC Delay and Coupling Noise Metric,” Proceeding 9 th  Great Lakes VLSI Symposium, March 1999, pp. 60-63), and Stable Three-Pole Model, S3P (Q. Yu and E. S. Kuh, “Moment Computation of Lumped and Distributed Coupled RC Trees with Application to Delay and Crosstalk Estimation,” Proceedings of the IEEE, Vol. 89, No. 5, pp. 772-788 (2001). The difference from the general model order reduction is that the above prior methods simply estimate peak value and its time of crosstalk noise, instead of waveform of crosstalk noise. Others existing American patents U.S. Pat. No. 5,481,695, U.S. Pat. No. 5,568,395, U.S. Pat. No. 5,596,506, U.S. Pat. No. 6,253,355, U.S. Pat. No. 6,253,359, U.S. Pat. No. 6,507,935, U.S. Pat. No. 6,536,022 and U.S. Pat. No. 6,662,149, etc. provide the application of crosstalk noise estimation. However, crosstalk noise in interconnects is probably non-monotonic waveform. None of the above estimating methods of capacitance crosstalk noise is suitable to estimate inductance crosstalk noise.  
         [0005]     Some existing traditional technologies provide delay and noise formula with considering self inductance and mutual inductance. However, this model only applies to double parallel net (Y. Cao, X. Huang, D. Sylvester, N. Chang and C. Hu, “A New Analytical Delay and Delay and Noise Formulas with Considering Self Inductances and Mutual Inductances,” Proceedings of IEDM 2000, 2000, pp. 823-826); Other traditional technologies provide analytic formula of RLC transmission line computation delay and overshooting voltage, but without research of the influence of inductance on crosstalk noise analysis (M. H. Chowdhury, Y, I. Ismail, C. V. Kashyap and B. L. Krauter, “Performance Analysis of Deep Sub Micro VLSI Circuits in the Presence of Self and Mutual Inductance,” Proceedings of ISCAS 2002, 2002, pp. 197-200); In the existing technologies, recursive algorithm is provided to compute RLC tress moment in linear time (C. L. Ratzlaff and L. T. Pillage, “RICE: Rapid Interconnect Circuit Evaluation Using AWE,” IEEE Transaction on Computer-Aided Design of Integrated Circuits and Systems, Vol. 13, No. 6, pp. 763-776 (1994), Q. Yu and E. S. Kuh, “Exact Moment Matching Model of Transmission Lines and Application to Interconnect Delay Estimation,” IEEE Transaction on VLSI Symposium, Vol. 3, No. 2, pp. 311-322 (1995). However, this technology does not provide moment formula of coupled circuit. The inventor&#39;s previous application, “Method of VLSI to estimate crosstalk noise in lumped RIC coupled interconnects” provided an algorithm to estimate crosstalk noise in circuit by means of lumped RLC model with an aim at RLC coupled trees in VLSI. However, section number of lumped circuit should be increased to make the result of simulation more accurate, and this, on the contrary, increases the program EMS memory loads and the whole computation time. In addition, the inventor&#39;s previous application, “Designing method and proof of nanometer VLSI to estimate crosstalk noise in distributed RIC coupled interconnects” once provided an algorithm to estimate crosstalk noise in circuit by means of unitary and uniform RLC model with an aim at RLC coupled trees in nanometer VLSI. However, in the existing circuit designing flow, the design of part of circuit adopts nonuniform distributed interconnects to optimize circuit operation. Therefore, unitary and uniform disturbed interconnects fail to analyze this special design.  
         [0006]     This invention aims at nonuniform distributed coupled RLC trees and carries out crosstalk noise estimation. The traditional technologies merely concentrate on distributed circuit that is necessary to circuit simulation. For example, R. Achar and M. S. Nakhla “Simulation on High-Speed Interconnects,” Proceeding IEEE, Vol. 89, No. 5, pp. 693-728, in 2001, A. C. Cangellaris, S. Pasha, J. L. Prince and M. Celik, “A New Discrete Transmission Line Model for Passive Model Order Reduction and Macromodeling of High-Speed Interconnections,” IEEE Transaction on Advanced Packing, Vol. 22, No. 3, pp. 356-364, in 1999, M. Celik and A. C. Cangellaris, “Simulation of Dispersive Multiconductor Transmission Lines by Pade Approximation Via the Lanczos Process,” IEEE Trans. Microwave. Theory Tech., Vol. 44, No. 12, pp. 2525-2533, in 1996, M. Celik and L. T. Pileggi, “Simulation of Lossy Multiconductor Transmission Lines Using Backward Euler Integration,” IEEE Trans. Circuits Syst. I-Fundam. Theor. AppI., Vol. 45, No. 3, pp. 238-243, in 1998, P. K. Gunupudi, R. Khazaka, M. S. Nakhla, T. Smy, and D. Celo, “Passive Parameterized Time-Domain Macromodels for High-Speed Transmission-Line Networks,” IEEE Trans. Microwave Theory Tech., Vol. 51, No. 12, pp. 2347-2354, in 2003, J. M. Wang, C. C. Chu, Q. Yu, and E. S. Kuh, “On Projection-Based Algorithms for Model-Order Reduction of Interconnects,” IEEE Trans. Circuits Syst. I-Fundam. Theor. Appl., Vol. 49, No. 11, pp. 1563-1585, in 2002 and Q. Xu and P. Mazumder, “Accurate Modeling of Lossy Nonuniform Transmission Lines by Using Differential Quadrature Methods,” IEEE Trans. Microwave Theory Tech., Vol. 50, No. 10, pp. 22233-2246 in 2002. However, up to now, there is still no efficient nonuniform distributed circuit model to estimate noise.  
       SUMMARY OF THE INVENTION  
       [0007]     This invention provides nonuniform distributed RLC coupled trees interconnects in nanometer VLSI, adopts moment matching method to efficiently estimate crosstalk noise while abandoning lumped circuit model with the traditional subsection method so as to shorten the computation time in the process of circuit simulation, and adopts nonuniform distributed model so as to be more approximate to the actual circuit design. Though in substance distributed circuit is a limitless series system, it approximately seems to be a limited series system in the form of multinomial moment model. Voltage and current moment model in nonuniform distributed circuit are approximate to a coordinate function multinomial, while the circuit parameter can be computed by means of data interpolation method. The moment model of every nonuniform distributed RLC coupled circuit includes one resistance, one independent power source and two independent voltage sources, which can reflect the information of resistance, capacitance, coupled capacitance, inductance, mutual inductance and moment. This invention additionally provides an order reduction moment computation formula. Each coefficient on multinomial moment model can be computed by order reduction computation and crosstalk estimating method of stable nonuniform distributed RLC coupled trees can be constructed by using projection-based order reduction method to compute the crosstalk noise of this order reduction series model, which can be regarded as the estimated value of crosstalk noise in the original circuit.  
       DETAILED DESCRIPTION OF THE INVENTION  
       [0008]      FIG. 1  is the input and output block diagram of computation device to implement this invention. This computation method should provide three input files, including circuit input signal  10 , series q  12  of order reduction mothod and circuit parameter  14  of nonuniform ditributed circuit in RLC coupled trees. Firstly, in step  16 , it is to compute the system moment {X 1 ,X 2 , . . . ,X q } of series q as the computation foundation of moment matching method, and then in step  18 , it is to construct {circumflex over (M)} and {circumflex over (N)} matrix in series q order reduction method. In step  20 , it is to compute the coefficient {b 1 ,b 2 , . . . ,b q } of equation |{circumflex over (N)}+s{circumflex over (M)}|=1+b 1 s+b 2 s 2 + . . . +b q s q . In step  22 , it is to update the moment values according to the input signals. In step  24 , as to the series reduction transformation function  
           V   ^     ⁡     (   s   )       =         a   0     +       a   1     ⁢   s     +   ⋯   +       a     q   -   2       ⁢     s     q   -   1             1   +       b   1     ⁢   s     +   ⋯   +       b     q   -   1       ⁢     s     q   -   1         +       b   q     ⁢     s   q                 
 of q poles, it is to compute the coefficient {a 0 ,a 1 , . . . ,a q-2 } based on the moment values in step  24  with moment matching method. In step  26 , it is to show {circumflex over (V)}(s) in the form of pole-residue,  
             V   ^     ⁡     (   s   )       =         k   1       s   -     p   1         +       k   2       s   -     p   2         +   ⋯   +       k   q       s   -     p   q             ,         
 and then to transform it into {circumflex over (v)}(t)=k 1 e p     1     1 +k 2 e p     2     1 + . . . +k q e p     q     1  by Inverse Laplace. Finally in step  28 , it is use {circumflex over (v)}(t) to estimate crosstalk voltage peak value. In step  30 , the computation is finished. Each step is expatiated as follows: 
 
         [0009]     This invention intends to use projection-base order reduction method to solve the above stability problem so as to create a stable order reduction method to solve the problem of crosstalk noise. This technique uses congruence transformation to project vector of the original n dimension to vector of order deduction q dimension, and q&lt;&lt;n, among which, q is determined in step  12 .  
         [0010]     In the nonuniform distributed circuit model provided by this invention, the crosstalk waveform can be expressed as {circumflex over (v)}(t)=k 1 e p     1     1 +k 2 e p     2     1 + . . . +k q e p     q     1 , among which, k i  and p 1 (1≦i≦q) are the pole and residue of q pole order reduction model {circumflex over (V)}(s). In order to ensure pole stability, it can be computed by using the root of equation |s{circumflex over (M)}+{circumflex over (N)}|=0, among which, matrix {circumflex over (M)} and matrix {circumflex over (N)} are obtained by integrated-congruence transforming matrix M and matrix N in the matrix of modified nodal analysis (MNA) by the prior technique.  
         [0011]     The prior technology provides that limitless series distributed circuit is simulated in limited series model by using integrated-congruence transform. By means of this technology, MNA formula can be expressed as the following formula:  
                   (       s   ⁢       [             M   ^     d         0       0           0       C       0           0       0       L         ]       ︸   M         +       [             N   ^     d           -     A   d           0             A   d   T         G         A   l             0         -     A   l   T           R         ]       ︸   N         )     ⁢       [               X   ^     d     ⁡     (   s   )                   V   n     ⁡     (   s   )                   I   L     ⁡     (   s   )             ]       ︸     X   ⁡     (   s   )             =         [               0             A   s                   0         ]       ︸   b       ⁢       V   s     ⁡     (   s   )           ,           (   1   )             
 
         [0012]     According to the circuit parameter provided in step  14 , matrix M includes matrix {circumflex over (M)} d , lumped capacitance matrix C and lumped inductance matrix L in nonuniform distributed order reduction method. Matrix N includes matrix {circumflex over (N)} d , lumped resistance matrix R, lumped conductance matrix G and A d , A l  incident matrix in nonuniform distributed order reduction model to balance Kirchhoff&#39;s Current Law (KCL) equation. Matrix X(s) is the transformation function of system condition variable, including system condition variable {circumflex over (X)} d (s) in nonuniform distributed order reduction model, voltage vector V n (s) in node and current vector I L (s) in resistance-inductance branch; Matrix b includes incident matrix A s  showing the connecting method of input signal V n (s) and circuit model. In formula (1), (s{circumflex over (M)} d +{circumflex over (N)} d ){circumflex over (X)} d (s)=A d V n (s) represents the circuit condition formula in nonuniform distributed order reduction model.  
         [0013]     X(s) expands in Taylor Series when frequency s=0 and series k system moment vector is X k [{circumflex over (X)} d,k  V u,k  I l,k ], among which, {circumflex over (X)} d,k , V n,k  and I L,k  represent the system moment of {circumflex over (X)} d (s), V n (s) and I L (s) respectively in series k. While the former q system moment can all be computed in step  16 .  FIG. 2  is the flow chart to compute transformation function moment of system condition variables, which will be expatiated in step  16 . The computation technique is expatiated as follows:  
         [0000]     Moment Model in Nonuniform Distributed RLC Coupled Circuit  
         [0014]     One group of RLC coupled trees includes several independent RLC decoupled trees, coupled capacitance and mutual inductance. Each RLC decoupled tree includes floating resistance and self inductance, as well as capacitance that connects tree node and the ground. If the root of one independent RLC tree connects with one input voltage source, this tree is called aggressor tree. On the contrary, if the root of this RLC tree directly connects with the ground, this tree is called victim tree. If self inductance and mutual inductance are deleted from the circuit, it turns into the regular RC tree circuit model in the traditional estimating technique of crosstalk noise. In this invention, coupled interconnects are transformed into RLC coupled trees to analyze crosstalk noise.  
         [0015]     The symbols are now detailed to demonstrate the complete RLC coupled trees. In consideration of N nonuniform distributed coupled transmission lines in  FIG. 3  and one small section of RLC tree T i  in a group of RLC coupled trees in  FIG. 4 , n j   i  is the j node in tree T i , F(n j   i ) is the father node of node f j   i . Line j   i  presents the nonuniform distributed RLC model between node n i   j  and F(n j   i ), among which, x=0 and x=d present the immediate end and remote end respectively. While r j   i (x), l j   i (x) and c j   i (x) represent resistance, inductance and capacitance of unit length respectively on Line j   i , and the conductance g j   i (x) of unit length is supposed to be neglected; i j   i (0,s) and v j   i (0,s) represent current and voltage multinomial of Line j   i  in the immediate end respectively, while i j   i (d,s) and v j   i (d,s) represent current and voltage multinomial of Line j   i  in the remote end respectively, c j,j     1     i,i     1   (x) and m j,j     1     i,i     1   (x) represent the coupled capacitance and mutual inductance of unit length between Line j   i  and Line j     1     i     1   ; cc j   i (x) and mm j   i (x) represent the aggregation of coupled capacitance and mutual inductance of Line j   i  respectively. P jk   i  represents the routes P j   i  of root from n j   i  to T i , and the common route of root route P k   i  from n k   i  to T i . The aggregation of n j   i  ancestor nodes is defined as A(n j   i ), including route P F(j)   i , i.e. all nodes in root route from F(n j   i ) to T i  so as to make Â(n j   i )={{A(n j   i )−n 0   i }∪n j   i }. On the contrary, D(n j   i )={n x   i |n j   i εA(n x   i )} represents the aggregation of n j   i  descendant nodes, which is other defined as {circumflex over (D)}(n j   i )={n j   i ∪D(n j   i )}. Generally speaking, coupling affect is not limited to take function on two nearest lines, especially inductance coupling affect. Therefore, this invention provides a model covering the general coupling condition. Each aggregation of cc j   i (x) and mm j   i (x) may probably include many coupled capacitance and mutual inductance, all of which can be solved efficiently.  
         [0016]     The voltage transformation function on node n j   i  is defined as V j   i (s), and the transformation function of current passing by n j   i  is defined as I j   i (s). V 0   i (s)=V s   i  represents the voltage of root n 0   i  in circuit trees, among which, V s   i  represents the voltage source connecting between root of tree T i  (i.e. n 0   i ) and the ground. In case V s   i =1, tree T i  is regarded as an aggressor tree. On the contrary, tree T i  can be regarded as a victim tree. V j   i (s) and I j   i (s) expands in Taylor Series in case s=0, then  
           V   j   i     ⁡     (   s   )       =       ∑     k   =   0     ∞     ⁢       V     j   ,   k     i     ⁢     s   k     ⁢           ⁢   and           
             I   j   i     ⁡     (   s   )       =       ∑     k   =   0     ∞     ⁢       I     j   ,   k     i     ⁢     s   k           ,       
 
 among which, V j,k   i  is called the voltage moment in series k of V j   i (s), and I j,k   i  is called the current moment in series k of I j   i (s). The voltage moment −V j,l   i  in the first series on node n j   i  is the common Elmore delay model. This invention will compute the moment V j,k   i  and I j,k   i  in series k according to each node n j   i  in tree structure. 
 
 Moment Computation in Nonuniform Distributed RLC Coupled Tree Interconnects 
 
         [0017]     This invention intends to transform the lumped circuit between RLC coupled trees n j   i  and its father node F(n j   i ) in the prior technology (the previous application “Method of VLSI to estimate crosstalk noise in lumped RIC coupled interconnects” by the inventor) into nonuniform distributed circuit Line j   i . Make v j   i (x,s), i j   i (x,s) and i e     j     i (x,s) represent the transformation functions to input signals of the voltage, current and capacitance current of any coordinate point x in Line j   i  respectively, among which, x=0 represents the immediate end of the line, x=d represents the remote end of the line. Making use of Laplace Transformation, the Telegrapher&#39;s Equation of the relation of voltage and current in Line j   i  can be shown as follows:  
                   ∂       v   j   i     ⁡     (     x   ,   s     )           ∂   x       =         -     (         r   j   i     ⁡     (   x   )       +       sl   j   i     ⁡     (   x   )         )       ⁢       i   j   i     ⁡     (     x   ,   s     )         -       ∑     mm   j   i       ⁢     (         sm     j   ,     j   1         i   ,     i   1         ⁡     (   x   )       ⁢       i     j   1       i   1       ⁡     (     x   ,   s     )         )           ,     
     ⁢               ∂       i   j   i     ⁡     (     x   ,   s     )           ∂   x       =         -       sc   jT   i     ⁡     (   x   )         ⁢       v   j   i     ⁡     (     x   ,   s     )         +       ∑     cc   j   i       ⁢     (         sc     j   ,     j   1         i   ,     i   1         ⁡     (   x   )       ⁢       v     j   1       i   1       ⁡     (     x   ,   s     )         )                     =       i     c   j     i     ⁡     (     x   ,   s     )                       (   2   )             
 
 Among which,  
           c   jT   i     ⁡     (   x   )       =         c   j   i     ⁡     (   x   )       +       ∑     cc   j   i       ⁢       c     j   ,     j   1         i   ,     i   1         ⁡     (   x   )               
 
 (x) represents all capacitance values on Line j   i , including self grounding capacitance and coupled capacitance aggregation. In consideration of  FIG. 2  in step  102 , v j,k   i (x), i j,k   i (x) and i c     j,k     i (x) represent the moments in series k in case v j   i (x,s), i j   i (x,s) and i c     j     i (x,s) expands in Taylor&#39;s Series respectively. In case k=0, the capacitance in circuit can take equivalent effect as open circuit. Therefore, the current moment in series 0 i c     i,0     i (x)=i j,0   i (d)=0, while the voltage moment in series 0 v j,0   i (x)=V j,0   i (0)=V s   i ; In case k&gt;0, the moment in series k is as follows:  
                   i     c     j   ,   k       i     ⁡     (   z   )       =           c   jT   i     ⁡     (   z   )       ⁢       v     j   ,     k   -   1       i     ⁡     (   z   )         -       ∑     cc   j   i       ⁢     (         c     j   ,     j   1         i   ,     i   1         ⁡     (   z   )       ⁢       v       j   1     ,     k   -   1         i   1       ⁡     (   z   )         )           ,           (   3   )                     i     j   ,   k     i     ⁡     (   x   )       =         i     j   ,   k     i     ⁡     (   d   )       +       ∫   x   d     ⁢         i     c     j   ,   k       i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z             ,           (   4   )                           v     j   ,   k     i     ⁡     (   x   )       =       ⁢         v     j   ,   k     i     ⁡     (   0   )       -       ∫   0   x     ⁢         r   j   i     ⁡     (   z   )       ⁢       i     c     j   ,   k       i     ⁡     (   z   )           -         R   j   i     ⁡     (   x   )       ⁢       i     j   ,   k     i     ⁡     (   x   )         -                     ⁢         ∫   0   x     ⁢         l   j   i     ⁡     (   z   )       ⁢       i     c     j   ,     k   -   1         i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z         -         L   j   i     ⁡     (   x   )       ⁢       i     j   ,     k   -   1       i     ⁡     (   x   )         -                     ⁢       ∑     mm   j   i       ⁢     (         ∫   0   x     ⁢         m     j   ,     j   1         i   ,     i   1         ⁡     (   z   )       ⁢       i       j   1     ,     k   -   1         i   1       ⁡     (   z   )       ⁢           ⁢     ⅆ   z         +         M     j   ,     j   1         i   ,     i   1         ⁡     (   x   )       ⁢       i       j   1     ,     k   -   1         i   1       ⁡     (   x   )           )               ,           (   5   )             
 
 Among which,  
             R   j   i     ⁡     (   x   )       =       ∫   0   x     ⁢         r   j   i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z           ,     
     ⁢         L   j   i     ⁡     (   x   )       =       ∫   0   x     ⁢         l   j   i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z           ,   and       
           M     j   ,     j   1         i   ,     i   1         ⁡     (   x   )       =       ∫   0   x     ⁢         m     j   ,     j   1         i   ,     i   1         ⁡     (   z   )       ⁢           ⁢     ⅆ   z             
 
 represent the progressive resistance, inductance and mutual inductance in the position of x on Line j   i  respectively. Formula (3) and (4) can be deduced by Kirchhoff&#39;s Current Law (KCL), and formula (5) can be deducted by Kirchhoff&#39;s Voltage Law (KVL). 
 
         [0018]     In order to simplify formula (4) and (5), the circuit current moment i c     j,k     i (x) and voltage moment v j,k   i (x) to multinomial in step  104  are approximated:  
                   i     c     j   ,   k       i     ⁡     (   x   )       =       ∑     n   =   0       m   k       ⁢       α     j   ,   k   ,   n     i     ⁢     x   n           ,     
     ⁢         v     j   ,   k     i     ⁡     (   x   )       =       ∑     n   =   0       p   k       ⁢       β     j   ,   k   ,   n     i     ⁢       x   n     .                   (   6   )             
 
 In addition, all circuit parameters, such as r j   i (x), l j   i (x), c j   i (x), cc j   i (x) and mm j   i (x) all approximate to q term multinomial, among which, each coefficient can be computed by Interpolation Technique. Therefore, step  106  is to compute the following multinomial from the multinomial multiplication integral in formula (5) with analytic method:  
             ∫   0   x     ⁢         r   j   i     ⁡     (   z   )       ⁢       I     c     j   ,   k       i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z         =       ∑     n   =   0       q   +     m   k     +   2       ⁢       a     j   ,   k   ,   n     i     ⁢     x   n           ,     
     ⁢           R   j   i     ⁡     (   x   )       ⁢       I     j   ,   k     i     ⁡     (   x   )         -=       ∑     n   =   0       q   +     m   k     +   2       ⁢       b     j   ,   k   ,   n     i     ⁢     x   n           ,     
     ⁢         ∫   0   x     ⁢         l   j   i     ⁡     (   z   )       ⁢       I     c     j   ,     k   -   1         i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z         =       ∑     n   =   0       q   +     m     k   -   1       +   2       ⁢       c     j   ,     k   -   1     ,   n     i     ⁢     x   n           ,     
     ⁢           L   j   i     ⁡     (   x   )       ⁢       I     j   ,     k   -   1       i     ⁡     (   x   )         =       ∑     n   =   0       q   +     m     k   -   1       +   2       ⁢       d     j   ,     k   -   1     ,   n     i     ⁢     x   n           ,     
     ⁢         ∫   0   x     ⁢         m     j   ,     j   1         i   ,     i   1         ⁡     (   z   )       ⁢       I     c       j   1     ,     k   -   1           i   1       ⁡     (   z   )       ⁢           ⁢     ⅆ   z         =       ∑     n   =   0       q   +     m     k   -   1       +   2       ⁢       e     j   ,     k   -   1     ,   n     i     ⁢     x   n           ,     
     ⁢           M     j   ,     j   1         i   ,     i   1         ⁡     (   x   )       ⁢       I       j   1     ,     k   -   1         i   1       ⁡     (   x   )         =       ∑     n   =   0       q   +     m     k   -   1       +   2       ⁢       f     j   ,     k   -   1     ,   n     i     ⁢       x   n     .               
 
         [0019]     It should be noted that all coefficients can be computed by means of recursive moment computation in the prior technology (the previous application “Method of VLSI to estimate crosstalk noise in lumped RIC coupled interconnects” by the inventor); In formula (6), it can be seen that the multinomial in series 0 i c     j,0     i (x) and  j,0   i (x)=V s   i  in case k=0, that is α j,0,0   i =0 and β j,0,0   i =V s   i  so m 0 =p 0 =0; m k =q+(k−1)(2q+2) and p k =k(2q+2) in case k&gt;0.  
         [0020]     In step  108 , it is to estimate whether a and β in coefficients can be computed, otherwise, return to step  106  by recursive computation. Step  110  is to finish the computation of current transformation function moment i c     j,k     i (x) and voltage transformation function moment v j,k   i (x) in system variable.  
         [0000]     Establishment of Matrix MNA of Simplified and Stable Pole Model  
         [0021]     In consideration of  FIG. 1  in step  18 , make the congruence transformation matrix Q=└X 0  X 1  . . . X q-1 ┘, then matrix MNA of the order reduction model can be computed by using the formula {circumflex over (M)}=Q T MQ and {circumflex over (N)}=Q T NQ. Make the elements in row k and row l in matrix {circumflex over (M)} and {circumflex over (N)} be {circumflex over (m)} k,l =X k-1   T MX l-1  and {circumflex over (n)} k,l =X k-1   T NX l-1  to observe the relation of different elements in {circumflex over (M)} and {circumflex over (N)}. We can discover the following two phenomena:  
         [0022]     1. {circumflex over (m)} ij =−X l-1   T NX j =−{circumflex over (n)} i,j+1 ;  
         [0023]     2. {circumflex over (m)} ij =X j−1   T MX l-1 =−X j−1   T NX i =−{circumflex over (n)} j,i+1 .  
         [0024]     We can see from the prior technology (the previous application “Method of VLSI to estimate crosstalk noise in lumped RIC coupled interconnects” by the inventor), the steps to compute different elements in matrix {circumflex over (N)} can be further simplified. By observing the elements in the first line and the first row in matrix {circumflex over (N)}, we can discover the following relation: 
        1. {circumflex over (n)} i1 =0;     2. {circumflex over (n)} i1 (i&gt;1), can be shown as I i,j−1   a , which is the same as the current moment in series (i−1) of node n 1   a  coming to the aggressor tree T a ;     3. {circumflex over (n)} 1i =−{circumflex over (n)} 1i .        
 
         [0028]     It can be computed by inserting the data into X k =[{circumflex over (X)} d,k  V n,k  I L,k ] T : 
 
 {circumflex over (m)}   k,l =−( {circumflex over (X)}   d,k-1   T   {circumflex over (M)}   d   {circumflex over (X)}   d,l-1   +V   n,k-1   T   CV   n,l-1   +I   L,k-1   T LI L,k-1 )   (7) 
 
         [0029]     We can discover the following relation by observing formula (1): 
 
NX 0 =b 
 
 NX   i+1   =−MX   i , for i=0,1,  . . . ,q′ 
 
 It can be computed by inserting the data into {circumflex over (n)} k,l =X k-1   T NX l-1   
                         n   ^       k   ,   l       =       -     X     k   -   1     T       ⁢     MX     l   -   2                     =     -     (           X   ^       d   ,     k   -   1       T     ⁢       M   ^     d     ⁢       X   ^       d   ,     l   -   2           +       V     n   ,     k   -   1       T     ⁢     CV     n   ,     l   -   2           +       I     L   ,     k   -   1       T     ⁢     LI     L   ,     l   -   2                         ,           (   8   )             
 
         [0030]      FIG. 5  is the detailed flow chart of Step  18 . From the above discussion, we can get the relation of different elements in matrix {circumflex over (M)} and matrix {circumflex over (N)}. Step  152  is the results of V n,k-1   T CV n,l-1 , V n,k-1   T CV m,l-1 , I L,k-1   T LI L,l-1  and I L,k-1   T LI L,l-2  in formula (7) and (8) related to the lumped circuit by means of the prior technology (the previous application “Method of VLSI to estimate crosstalk noise in lumped RIC coupled interconnects” by the inventor), as shown in step  156 ; While in step  158 , {circumflex over (X)} d,k-1   T {circumflex over (M)} d {circumflex over (X)} d,l-1  and {circumflex over (X)} d,k-1   T {circumflex over (M)} d {circumflex over (X)} d,l-2  are related to nonuniform distributed circuit. After computing different coefficients of system variable transformation function in step  16 , insert them into the formula to get the result. It should be noted that {circumflex over (X)} d,k-1   T {circumflex over (M)} d {circumflex over (X)} d,l-2  can be computed by the influence caused by different nonuniform distributed circuits, while the result of Line j   i  is as follows:  
             ∫   0   d     ⁢         v     j   ,     k   -   1       i     ⁡     (   z   )       ⁢       i     C     j   ,     l   -   1         i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z         =       ∑     n   =   0         p     k   -   1       +     m     l   -   1           ⁢         g     j   ,   n     i       n   +   1       ⁢     d     n   +   1             ,     
     ⁢         ∫   0   d     ⁢         i     j   ,     k   -   1       i     ⁡     (   z   )       ⁢     (           L   j   i     ⁡     (   z   )       ⁢       i     j   ,     l   -   3       i     ⁡     (   z   )         +       ∑     mm   j   i       ⁢     (         M     j   ,     j   1         i   ,     i   1         ⁡     (   z   )       ⁢       i         j   1     ⁢   l     -   3       i   1       ⁡     (   z   )         )         )     ⁢           ⁢     ⅆ   z         =       ∑     n   =   0       q   +     m     k   -   1       +     m     l   -   3       +   3       ⁢         h     j   ,   n     i       n   +   1       ⁢       d     n   +   1       .                 
         [0031]     Alike, the coefficients g j,n   i  and h j,n   i  in the above formula can be computed from multinomial by multiplication integral with analytic method.  
         [0032]     In consideration of the circuit of two grounding capacitances and one coupled capacitance in  FIG. 6 ( a ), in case k&gt;0, the output currents of n j   i  and n j     1     i     1    are as follows: 
 
 i   c     j,k     i =( c   j,0   i   +c   j,j     1     i,i     1   ) v   j,k-1   i   −c   j,j     1     i,i     1     v   j     1     ,k-1   i     1   , 
 
 i   c     j,k   =( c   j     1     ,0   i     1     +c   j,j     1     i,i     1   ) v   j     1     ,k-1   i     1     −c   j,j     1     i,i     1     v   j     1     ,k-1   i     1   , 
 
         [0033]     Therefore, coupled capacitance can be regarded as two current sources. When there are many coupled capacitances in circuit, the model of each decoupled current moment is as follows:  
         i     c     j   ,   k       i     =         c   jT   i     ⁢     v     j   ,     k   -   1       i       -       ∑       c     j   ,     j   1         i   ,     i   1         ∈     cc   j   i         ⁢       c     j   ,     j   1         i   ,     i   1         ⁢       v       j   1     ,     k   -   1         i   1       .               
 
         [0034]     The current moment i j,k   i  in series k is the aggregation of capacitance current source in series k corresponding to each node after node  
             n     j   ,   k     i     ·     i     j   ,   k     i       =       ∑       n   x   i     ∈       D   ^     ⁡     (     n   j   i     )           ⁢     i     C     x   ,   k       i         ,       
 
 then the equivalent circuit of coupled capacitance is as shown in  FIG. 6 ( b ). 
 
         [0035]     Finally, in step  160 , the moment model in RLC coupled circuit can be established, as shown in  FIG. 7 . In case x=0, we can get from formula (4) that 
 
 i   j,k   i (0)= i   j,k   i ( d )+ J   j,k   i    (9) 
 
 Among which,  
         J     j   ,   k     i     =         ∫   d   0     ⁢         I     c     j   ,   k       i     ⁡     (   z   )       ⁢           ⁢     ⅆ   z         =       ∑     n   =   0       m   k       ⁢       α     j   ,   k   ,   n     i       n   +   1               
 
 represents the aggregation of all capacitance currents on Line j   i . Alike, formula (5) can be also expressed as  
                   v     j   ,   k     i     ⁡     (   d   )       =         V     j   ,   k     i     ⁡     (   0   )       -       ∑     n   =   0       q   +     m   k     +   2       ⁢       b     j   ,   k   ,   n     i     ⁢     d   n         -     E     LM     i   ,   j       i     -     E     j   ,   k     i         ⁢     
     ⁢       Among   ⁢           ⁢   which     ,     
     ⁢       E     LM     j   ,   k       i     =       ∑     n   =   0       q   +     m     k   -   1       +   2       ⁢       (       d     j   ,     k   -   1     ,   n     i     +       ∑     mm   j   i       ⁢     f     j   ,     k   -   1     ,   n     i         )     ⁢     d   n           ,     
     ⁢       E     j   ,   k     i     =         ∑     n   =   0       q   +     m   k     +   2       ⁢       a     j   ,   k   ,   n     i     ⁢     d   n         +       ∑     n   =   0       q   +     m     k   -   1       +   2       ⁢       (       c     j   ,     k   -   1     ,   n     i     +       ∑     mm   j   i       ⁢     e     j   ,     k   -   1     ,   n     i         )     ⁢       d   n     .                       (   10   )             
 
 E j,k   i  and E LM     j,k     i  represent the sudden decrease of voltage in voltage moment v j,k   i (0) when the capacitance current moment in series k and series k-1 passes by the resistance and induction on Line j   i . When the circuit is RC coupled circuit, E LM     j,k   =0. 
 
 Combination of Nonuniform Distributed Coupled Circuit with RLC Coupled Trees 
 
         [0036]     In this invention, line(n j   i ) is used to represent interconnects between n j   i  and F(n j   i ). In case line(n j   i )=1, it is to represent that there is one net between n j   i  and F(n j   i ), otherwise, line(n j   i )=0. R j   i  and L j   i  are the resistance and inductance on line(n j   i ). C j,0   i  is the grounding capacitance of n j   i ; C j,j     1     i,i     1    is the coupled capacitance between n j   i  and n j     1     i     i   ; M j,j     1     i,j     1    is the mutual inductance between L j   i  and L j     1     i     1   ; CC j   i  is the aggregation that take effect of coupled capacitance with n j   i ; MM j   i  is the aggregation that take effect of mutual inductance with L j   i ; S(n j   i ) is the aggregation of descendant nodes after n j   i .  
         [0037]     In the computation of circuit model moment, in order to process nonuniform distributed coupled circuit at the same time, the current moment I j,k   i  in series k in lumped circuit in the prior technology (the previous application “Method of VLSI to estimate crosstalk noise in lumped RIC coupled interconnects” by the inventor) is applied and modified as:  
                 I     j   ,   k     i     =       I     C     j   ,   k       i     +       ∑       n   y   i     ∈     S   ⁡     (     n   j   i     )           ⁢     I     y   ,   k     i           ⁢     
     ⁢       Among   ⁢           ⁢   which     ,     
     ⁢       I     y   ,   k     i     =       I     y   ,   k     i     +       line   ⁡     (     n   y   i     )       ·       J     y   ,   k     i     .                     (   11   )             
 
         [0038]     Each current moment can be computed in the direction from leaf node in T i  to root node. The relation between voltage moments V j,k   i  and V F(i),k   i  is as follows:  
               V     j   ,   k     i     =       V       F   ⁡     (   j   )       ,   k     i     -       R   j   i     ⁢     I     j   ,     k   -   1       i       -       L   j   i     ⁢     I     j   ,     k   -   1       i       -       ∑     MM   j   i       ⁢       M     j   ,     j   1         i   ,     i   1         ⁢     I       j   1     ,     k   -   1         i   1           -     line   ⁢           ⁢       (     n   j   i     )     ·       (       E     j   ,   k     i     +     E     LM     i   ,   j       i       )     .                   (   12   )             
 
         [0039]     The complexity of recursive moment computations in nonuniform distributed circuit provided by this invention is O(nk 2 ), among which, n is the number of nodes in tree model. On the other hand, the computation complexity applied in lumped circuit model is o(mk), among which, m is the number of nodes in the lumped circuit model. Generally speaking, in order to make the result of simulation of lumped circuit more accurate, it is to make m&gt;nk, so the complexity of model moment computation of nonuniform distributed circuit provided in this invention is less than that of the computation in lumped circuit model.  
         [0000]     Update Moment Values According to the Input Signals  
         [0040]     In the previous moment computation, input waveform is supposed to expand to the frequency domain under the step function. However, the input signals in step  10  are probably random signals, which make the transformation function in step  22  after moment update as follows:  
         v   ⁡     (   s   )       =         m   1   ′     ⁢   s     +       m   2   ′     ⁢     s   2       +       m   3   ′     ⁢     s   3       +       m   4   ′     ⁢     s   4       +       m   5   ′     ⁢     s   5     ⁢     ⋯   ⁢           .             
 
 For example, if the input signal in step  10  is ramp function, it can be expressed as:  
           v   ⁡     (   t   )       =         t   τ     ⁢     u   ⁡     (   t   )         -       t   τ     ⁢     u   ⁡     (     t   -   τ     )         +     u   ⁡     (     t   -   τ     )           ,       
 
 Among which, u(t) represents series function and 1/τ is the ramp rate of ramp function. After x(t) processes Laplace Transform, it can conclude:  
         V   ⁡     (   s   )       =       1   s     ⁢       (     1   -       τ   2     ⁢   s     +         τ   2     6     ⁢     s   2       -         τ   3     24     ⁢     s   3       +         τ   4     120     ⁢     s   4       +   ⋯     ⁢           )     .           
 
 After coefficient matching, it can conclude:  
                 m   1   ′     =     m   1                   m   2   ′     =       m   2     -       τ   2     ⁢     m   1                       m   3   ′     =       m   3     -       τ   2     ⁢     m   2       +         τ   2     6     ⁢     m   1                       m   4   ′     =       m   4     -       τ   2     ⁢     m   3       +         τ   2     6     ⁢     m   2       -       τ   24     ⁢     m   1                       m   5   ′     =       m   5     -       τ   2     ⁢     m   4       +         τ   2     6     ⁢     m   3       -         τ   3     24     ⁢     m   2       +         τ   4     120     ⁢     m   1                   ⋮         .       
 
 After moment update computation, it can conclude the voltage moment of each node in interconnects under random waveform input. 
 
 Crosstalk Noise Estimation in Nonuniform Distributed RLC Coupled Trees 
 
         [0041]     In step  20 , it is to apply matrix {circumflex over (N)} and {circumflex over (M)} in step  18  to compute the coefficient {b 1 b 2 , . . . ,b q } of formula |{circumflex over (N)}+s{circumflex over (M)}|=1+b 1 s+b 2 s 2 + . . . +b q s q . Later in step  24 , make the order reduction formula {circumflex over (V)}(s) of q pole as follows:  
                 V   ^     ⁡     (   s   )       =           a   0     +       a   1     ⁢   s     +   ⋯   +       a     q   -   2       ⁢     s     q   -   1             1   +       b   1     ⁢   s     +   ⋯   +       b     q   -   1       ⁢     s     q   -   1         +       b   q     ⁢     s   q           .             (   13   )             
 
         [0042]     Therefore, when time t approximates to 0 or 8, its approximate crosstalk noise {circumflex over (v)}(t)=0. It is to make use of 2q-1 moments {V 1 ,V 2 , . . . ,V 2q-1 } of the original model to compute the unknown coefficient a i (0≦i≦q-2).  
         [0043]     In step  26 , formula (13) is shown in the pole-residue form:  
                   V   ^     ⁡     (   s   )       =         k   1       s   -     p   1         +       k   2       s   -     p   2         +   ⋯   +       k   q       s   -     p   q             ,           (   14   )             
 
 Among which, p i ,i=1,2, . . . , q is the pole of {circumflex over (V)}(s), k i  is the residue corresponding to each pole p i . It can be concluded by Inverse Laplace Transformation: 
 
 {circumflex over (v)} ( t )= k   i   e   p     1     i   +k   2   e   p     2     i   + . . . +k   q   e   p     q     i . 
 
 If crosstalk {circumflex over (V)}(s) reaches to the peak value in case t=t m , then {circumflex over (v)}′(t n )=0 and {circumflex over (v)}″(t m )&lt;0. v(t m ) is the estimated value of required crosstalk noise. 
 
 Simple Implementing Case 
 
         [0044]     In order to prove the correctness of computation provided in this invention,  FIG. 8  provides a model with three coupled circuits to research the crosstalk estimation technique by using this model. However, the application of this computation is not limited to this structure. In the figure, rectangle □ represents root nodes in the tree model, and round ∘ represents leaf node in the tree model.  
         [0045]     2003 International Technology Roadmap of Semiconductors (ITRS) is introduced to the circuit parameter in the circuit model, among which, under the 90 nanometer semiconductor producing technology, the coefficient of line resistance is 22 mΩ-μm and the coefficient of dielectric value is 3.1. In the implementing case of this invention, it is to suppose that with the same width 10.88 μm, same thickness 0.58 μm and same height from the substrate 0.58 μm of all unitary and uniform metal lines, the resistance in unit length of metal line is 3.5 mΩ/μm and the grounding capacitance in unit length is 0.516 fF/μm after computation. Now it is to suppose that in nonuniform distributed circuit: 
 
Line resistance is 3.50−8.53·10 −3 ×+1.05·10 −4 x 2  mΩ/˜m, 
 
Grounding capacitance is 0.55+3.31·10 −3 x−1.32·10 −5 x 2  fF/μm. 
 
 In addition, this implementing case adopts the inductance parameter in unit length of unitary and uniform circuit and introduces the data from the prior technology (Provided by A. Deutsch et al., “When are Transmission-Line Effects Important for On-Chip Interconnections?,” IEEE Trans. Microwave Theory Tech., Vol. 45, No. 10, pp. 1836-1846, in 1997). The inductance is 0.347 pH/μm. Now it is to suppose that in ununiform circuits, 
 
The inductance is 0.27−6.60·10 −4 x+8.09·10 −6 x 2  pH/μm. 
 
         [0046]     In this implementing case, coupled capacitance is supposed to be 0.47+6.61·10 −3 x−2.63·10 −5 x 2  fF/μm and the inductance is 0.12+6.60·10 −4 x−8.09·10 −6 x 2  pH/μm to prove the correctness of estimation device of crosstalk noise in this invention. Finally, the loading capacitance is supposed to be 50 fF. Noise peak values and their occurring time in different circuits should be taken into consideration due to the difference in structures, including length, coupled position, effective driving resistance and ascending time, etc. As shown in  FIG. 8 , there are five lengths of coupled line in net  1  L 1 ={2,3,4,5}(mm), while there are also five lengths of coupled line in net  2 , among which, the later is shorter than the former. In addition, other branches in FIGS.  8 ( b ) and ( c ) are all 1 mm. In the testing case, the topology of net  1  remains fixed, while the coupled positions of net  2  will change: moving from the immediate end of net  1  to the remote end of line  1  with space of 1 mm. In each testing case, line  1  and line  2  are activated independently. In addition, the immediate ends of two lines connect with four effective driving resistances respectively: 3O-3O, 3O-30O, 30O-3O and 30O-30O. Moreover, the voltage source connecting on the aggressor tree includes two ascending times: the ramp functions 0.02 ns and 0.2 ns, and its unit intensity is normalized. Therefore, the computation provided in this invention may totally apply in 1640 testing cases.  
         [0047]     In this invention, it is to make comparison with the traditional one-pole (1P) model and two-pole (2P) model, as well as three-pole (3P) model, four-pole (4P) model, five-pole (5P) and six-pole (6P) model in this invention. Table I lists the absolute error and comparative error by comparing the simulation result of crosstalk peak value and commercial software HSPICE, among which, resistance, capacitance and inductance are set to be put into sections per 20 μm by HSPICE. Among 1640 testing cases, there are 40 cases with unstable poles in 1P model; there are 15 cases with unstable poles in 2P model. In order to compare the efficiency and correctness of moment computation complexity in distributed model and lumped model, Table II lists the simulation results and their comparative errors of moment computation time in S6P lumped RLC trees, among which, the testing case is to cut the length of 1 mm into different sections.  
         [0048]     In this invention, the phenomena observed will conclude in the following items: 
        1. The model by applying the method provided by this invention outperforms the traditional 1P and 2P models. Therefore, these traditional models are not suitable to apply in RLC coupled trees. Moreover, increase of the series of order reduction model will make the result more accurate.     2. In the absolute error list of S3P model in Table I, we can find that the average error is less than 10%, which seems to be quite suitable to crosstalk noise estimation. However, the comparative error of S3P model is not as accurate as expected. By observing the simulation result, we can find the computation efficiency and estimation accuracy of S6P model can get better balance point.  FIG. 9  shows the crosstalk waveform of Spice, S3P, S4P and S6P in case the coupled circuit L 1 =L 2 =1 mm. We can clearly find that the waveform of S6P model is more accurate than the waveform of S3P and Spice.     3. In Table I, the stimulated computation time of S6P distributed circuit is 29.56 seconds with comparative error of 6.38%. While Table II shows that under the same comparative error, the simulation of lumped RLC circuit takes 902.13 seconds. Obviously, the efficiency and correctness of distributed circuit model are better than the lumped one.        
 
         [0052]     In short, this invention provides a method for efficiently estimating crosstalk noise of nanometer VLSI interconnects, which can quickly estimate crosstalk noise in circuit nodes by cooperating with the present VLSI design flow. In this invention, VLSI interconnects are regarded to be RLC coupled trees including nonuniform distributed circuits and lumped ones, and projection-based recursive formulas of moment computations is provided to estimate the crosstalk noise waveform of circuit inductance.  
                                                                           TABLE I                           Comparison list of absolute errors and comparative       errors of crosstalk noise waveform                1P   2P   S3P   S4P   S5P   S6P                        Absolute                               errors (%)       Maximum   104.21   82.30   20.16   25.53   19.77   12.81       Average   16.27   12.66   4.57   3.16   2.23   1.26       Minimum   0.04   0.18   0.02   0   0   0       Comparative       errors       (%)       Maximum   312.44   182.24   57.85   73.48   63.57   49.89       Average   63.52   43.08   20.33   15.93   11.34   6.38       Minimum   0.12   1.49   0.03   0.02   0   0                  
 
         [0053]    
       
         
               
             
               
               
             
               
               
               
               
               
               
             
               
               
               
               
               
               
             
           
               
                 TABLE II 
               
             
             
               
                   
               
               
                   
               
               
                 S6P model is expressed in lumped RLC coupled 
               
               
                 trees. 1 mm is cut into different sections to make comparison of 
               
               
                 computation time and their comparative errors. 
               
             
          
           
               
                   
                 Number of sections 
               
             
          
           
               
                   
                 2 
                 3 
                 4 
                 5 
                 6 
               
               
                   
                   
               
             
          
           
               
                 Time 
                 32.34 
                 56.35 
                 63.46 
                 81.03 
                 95.22 
               
               
                 (second) 
               
               
                 Comparative 
                 8.80 
                 7.68 
                 7.18 
                 6.97 
                 6.80 
               
               
                 errors 
               
               
                 (%) 
               
               
                 Number 
                 10 
                 20 
                 30 
                 40 
                 50 
               
               
                 of 
               
               
                 sections 
               
               
                 Time 
                 167.06 
                 325.23 
                 478.96 
                 657.48 
                 902.13 
               
               
                 (second) 
               
               
                 Comparative 
                 6.63 
                 6.46 
                 6.42 
                 6.40 
                 6.39 
               
               
                 errors 
               
               
                 (%) 
               
               
                   
               
             
          
         
       
     
     
    
     BRIEF DESCRIPTION OF THE DRAWINGS  
       [0054]      FIG. 1  is the input and output block diagram of computation device of implementing this invention.  
         [0055]      FIG. 2  is the computation flow chart to compute system moment.  
         [0056]      FIG. 3  is the nonuniform distributed transmission lines in coupled trees.  
         [0057]      FIG. 4  is the classic expression of RLC coupled trees T i .  
         [0058]      FIG. 5  is the flow chart of accomplishing nonuniform distributed circuit.  
         [0059]      FIG. 6  is the conditions of coupled capacitance between two nodes: (a) is the original circuit; (b) is the equivalent moment model.  
         [0060]      FIG. 7  is the moment model of nonuniform distributed RLC coupled circuit.  
         [0061]      FIG. 8  is the condition of three RLC coupled trees with two parallel nets: (a) two parallel nets, (b) tree  1 , (c) tree  2 , among which, there are five lengths of coupled line of net  1  L 1 ={2,3,4,5}(mm), while there are also five lengths of coupled line of net  2  L 2 ={2,3,4,5}(mm), among which, the later is shorter than the former.  
         [0062]      FIG. 9  is the simulation of crosstalk noise peak value of Spice, S3P, S4P and S6P.