CN114004191A - Method, system, device and medium for extracting delay circuit macro model - Google Patents
Method, system, device and medium for extracting delay circuit macro model Download PDFInfo
- Publication number
- CN114004191A CN114004191A CN202111161438.4A CN202111161438A CN114004191A CN 114004191 A CN114004191 A CN 114004191A CN 202111161438 A CN202111161438 A CN 202111161438A CN 114004191 A CN114004191 A CN 114004191A
- Authority
- CN
- China
- Prior art keywords
- matrix
- delay circuit
- parameters
- transfer function
- macro model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 239000011159 matrix material Substances 0.000 claims abstract description 150
- 238000012546 transfer Methods 0.000 claims abstract description 37
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 230000006870 function Effects 0.000 claims description 35
- 238000004590 computer program Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 7
- 238000004422 calculation algorithm Methods 0.000 claims description 6
- 238000000605 extraction Methods 0.000 claims description 5
- 230000003071 parasitic effect Effects 0.000 claims description 5
- 230000008054 signal transmission Effects 0.000 claims description 4
- 238000009795 derivation Methods 0.000 claims description 3
- 238000000638 solvent extraction Methods 0.000 claims description 3
- 238000004088 simulation Methods 0.000 abstract description 15
- 238000004458 analytical method Methods 0.000 abstract description 12
- 230000009467 reduction Effects 0.000 description 8
- 230000008569 process Effects 0.000 description 6
- 238000010586 diagram Methods 0.000 description 4
- 238000013178 mathematical model Methods 0.000 description 4
- 230000014509 gene expression Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 150000001875 compounds Chemical class 0.000 description 2
- 230000003111 delayed effect Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000012938 design process Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/39—Circuit design at the physical level
- G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/32—Circuit design at the digital level
- G06F30/33—Design verification, e.g. functional simulation or model checking
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Design And Manufacture Of Integrated Circuits (AREA)
Abstract
The invention discloses a method for extracting a macro model of a delay circuit, which comprises the following steps: obtaining a delay circuit of a macro model to be extracted; embedding the delay circuit into preset parameters to obtain a delay circuit with parameters; calculating a transfer function of the delay circuit with the parameters and calculating all moments of the transfer function; constructing a projection matrix according to all moments of the transfer function; projecting the delay circuit with parameters in the projection matrix to obtain a macro model. The invention also discloses a system, a computer device and a readable storage medium. The scheme provided by the invention solves the problem of simulation analysis complexity in a high-speed integrated circuit by reducing the order of the model, and can reduce the calculation amount under the condition of keeping a certain calculation precision.
Description
Technical Field
The invention relates to the field of integration, in particular to a method, a system, equipment and a storage medium for extracting a macro model of a delay circuit.
Background
In the design process of the integrated circuit, one of the essential links is to perform simulation verification on the integrated circuit. With the increasing scale of integrated circuits, simulation verification is an important factor restricting the development of integrated circuits. The EDA software also provides various methods for simulation analysis of large-scale integrated circuits to quickly and accurately perform simulation verification on the designed circuits. However, as the operating frequency of the chip is higher and the number of components is higher, simulation analysis becomes a great challenge in the design process of the integrated circuit. For simulation analysis of large-scale analog integrated circuits, accurate mathematical modeling of original circuits is required, and various effects in the circuit operation process, such as transmission delay of signals, parasitic parameters of devices, and the like, need to be taken into consideration. The mathematical models of integrated circuits that take into account a number of factors are large and complex. It is time consuming to directly solve the exception, sometimes even impossible, which presents a significant challenge to EDA software.
Establishing a macro model of an original large-scale circuit becomes an important means for solving the problems of circuit simulation and analysis. Model reduction is an effective way to find circuit macro models. The main purpose of model order reduction is to provide redundant information in the original large-scale system and find an approximate smaller order reduction system, which can well approximate the input-output relationship of the original system and can retain the main properties of the original system, such as passivity, stability and the like. Therefore, the scale of the original system can be reduced to a great extent, and the simulation analysis difficulty and complexity of the original system are reduced. The main mechanism is that the main information of the high-dimensional space can be often represented through the low-dimensional space, and the original system is approximately represented in the low-dimensional space by constructing a proper low-dimensional space.
High speed interconnect circuits need to take the delay of the signal into account, so their mathematical model is typically a system of delayed differential equations. Directly solving a large-scale delay differential equation set is time-consuming, and the simulation analysis timeliness of the original delay circuit is seriously influenced. At present, a macro model establishing method for a large-scale delay circuit system mainly comprises a Taylor expansion method and a Laguerre expansion method, but the methods are not high enough in precision, and the obtained macro model cannot well approximate the input and output characteristics of an original system and the important properties of stability passivity and the like, so that the effective macro model establishing method for the delay circuit system is a problem which needs to be solved urgently in the field of circuit analysis and simulation at present.
Integrated circuit interconnect systems can generally be mathematically modeled using Differential Equations (differentials) but as the operating frequencies of integrated circuits increase, the propagation time of signals on the interconnect lines needs to be taken into account to meet timing requirements in high speed circuits, so that the entire integrated circuit interconnect system is built up as a system of delayed Differential Equations. The simulation analysis of the circuit needs to solve the delay differential equation set, and in the case of a large-scale integrated circuit interconnection system, the solution of the delay differential equation set is usually time-consuming and seriously affects the development cycle of the integrated circuit design. Particularly, in a differential algebra delay system, the solving process is more complex, and the working efficiency of the EDA software is seriously influenced.
Disclosure of Invention
In view of the above, in order to overcome at least one aspect of the above problems, an embodiment of the present invention provides a method for extracting a macro model of a delay circuit, including:
obtaining a delay circuit of a macro model to be extracted;
embedding the delay circuit into preset parameters to obtain a delay circuit with parameters;
calculating a transfer function of the delay circuit with the parameters and calculating all moments of the transfer function;
constructing a projection matrix according to all moments of the transfer function;
projecting the delay circuit with parameters in the projection matrix to obtain a macro model.
In some embodiments, in obtaining the delay circuit of the macro model to be extracted, the delay circuit is:
wherein the matrix E ∈ Rn×nIs a singular matrix formed by inductance L and capacitance C, the matrix A belongs to Rn×nIs composed of resistance, conductance values and correlation matrix between devices, matrix Ad∈Rn×nFormed by values of parasitic parameters in the interconnection line, B ∈ Rn×1Is an input matrix, u (t) is an input, C ∈ R1×nFor the output matrix, the unknown variable x (t) e Rn×1N is the order of the circuit system, namely the unknown quantity in the whole delay circuit, and tau is the signal transmission delay.
In some embodiments, the delay circuit is embedded with preset parameters to obtain a delay circuit with parameters, the delay circuit with parameters is:
In some embodiments, calculating a transfer function of the delay circuit with the parameters and calculating all moments of the transfer function further comprises:
calculating a laplace transform of the delay circuit with the parameter, wherein the laplace transform of the delay circuit with the parameter is H (s, epsilon) ═ C (s (E + epsilon E)0)-A-Ade-sτ)-1B;
Computing a Taylor expansion of the Laplace transformed delay circuit with the parameters, wherein the Taylor expansion is
Using (s, epsilon) as parameter pair Performing Taylor expansion to obtainWherein, thereinIs the (i, j) th moment of the transfer function H (s, epsilon) and the matrix C is constant, thenCan be expressed as
in some embodiments, constructing a projection matrix from all moments of the transfer function further comprises:
Constructing a second-order Krylov subspace about the R matrix:
wherein i is 1,2, …, p0+1,j=1,2,…,q0+1;
When j is 1, an orthogonal basis matrix V of a first row matrix in the matrix R is calculated1:
when j is 2, …, q0When +1, calculate the orthogonal base matrix L formed by the jth row matrix in the matrix RjThe j-th row matrix satisfies the following second-order Krylov subspace relationship:
Orthogonalizing the jth row matrix by using a second-order Arnoldi algorithm to obtain an orthogonal base matrixFor matrixPartitioning:wherein
Calculating an orthogonal base matrix L of a jth row matrix in a matrix Rj:
colspan{Q2}=colspan{Lj}
Using the first j-1 row orthogonal matrix Vj-1And LjConstruction of a New matrix [ V ]j-1 Lj]To matrix [ V ]j-1 Lj]Orthogonalizing to obtain VjAccording to j 2, …, q0+1 calculation matrix R front q0Orthogonal basis matrix of +1 row matrix, Vj=orth([Vj-1Lj])。
In some embodiments, projecting the delay circuit with parameters in the projection matrix to obtain a macro model, further comprises:
using projection matrices VjReducing the order of the delay circuit with parameters to obtain a macro model:
Based on the same inventive concept, according to another aspect of the present invention, an embodiment of the present invention further provides a delay circuit macro model extraction system, including:
the acquisition module is configured to acquire a delay circuit of the macro model to be extracted;
the parameter embedding module is configured to embed the delay circuit into preset parameters to obtain the delay circuit with parameters;
a calculation module configured to calculate a transfer function of the delay circuit with the parameters and to calculate all moments of the transfer function;
a construction module configured to construct a projection matrix from all moments of the transfer function;
a projection module configured to project the delay circuit with parameters in the projection matrix to obtain a macro model.
Based on the same inventive concept, according to another aspect of the present invention, an embodiment of the present invention further provides a computer apparatus, including:
at least one processor; and
a memory storing a computer program operable on the processor, wherein the processor executes the program to perform any of the steps of the delay circuit macro model extraction method described above.
Based on the same inventive concept, according to another aspect of the present invention, an embodiment of the present invention further provides a computer-readable storage medium storing a computer program, which when executed by a processor performs the steps of any one of the methods for extracting a delay circuit macro model as described above.
The invention has one of the following beneficial technical effects: the scheme provided by the invention solves the problem of simulation analysis complexity in a high-speed integrated circuit by reducing the order of the model, and can reduce the calculation amount under the condition of keeping a certain calculation precision.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other embodiments can be obtained by using the drawings without creative efforts.
Fig. 1 is a schematic flowchart of a method for extracting a macro model of a delay circuit according to an embodiment of the present invention;
FIG. 2 is a diagram illustrating relative errors between a reduced order system and an original system according to an embodiment of the present invention;
fig. 3 is a schematic structural diagram of a delay circuit macro model extraction system according to an embodiment of the present invention;
FIG. 4 is a schematic structural diagram of a computer device provided in an embodiment of the present invention;
fig. 5 is a schematic structural diagram of a computer-readable storage medium according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the following embodiments of the present invention are described in further detail with reference to the accompanying drawings.
It should be noted that all expressions using "first" and "second" in the embodiments of the present invention are used for distinguishing two entities with the same name but different names or different parameters, and it should be noted that "first" and "second" are merely for convenience of description and should not be construed as limitations of the embodiments of the present invention, and they are not described in any more detail in the following embodiments.
According to an aspect of the present invention, an embodiment of the present invention provides a method for extracting a macro model of a delay circuit, as shown in fig. 1, which may include the steps of:
s1, obtaining a delay circuit of the macro model to be extracted;
s2, embedding the delay circuit into preset parameters to obtain a delay circuit with parameters;
s3, calculating the transfer function of the delay circuit with parameters and calculating all moments of the transfer function;
s4, constructing a projection matrix according to all moments of the transfer function;
s5, projecting the delay circuit with parameters in the projection matrix to obtain a macro model.
In some embodiments, in obtaining the delay circuit of the macro model to be extracted, the delay circuit is:
wherein the matrix E ∈ Rn×nIs a singular matrix composed of inductance L and capacitance C, and n belongs to Rn×nIs composed of resistance, conductance values and correlation matrix between devices, matrix Ad∈Rn×nFormed by values of parasitic parameters in the interconnection line, B ∈ Rn×1Is an input matrix, u (t) is an input, C ∈ R1×nFor the output matrix, the unknown variable x (t) e Rn×1As nodes in a circuitVoltage or branch current, n is the order of the circuit system, that is, the unknown quantity in the whole delay circuit, and τ is the signal transmission delay.
In some embodiments, the delay circuit is embedded with preset parameters to obtain a delay circuit with parameters, the delay circuit with parameters is:
In some embodiments, calculating a transfer function of the delay circuit with the parameters and calculating all moments of the transfer function further comprises:
calculating a laplace transform of the delay circuit with the parameter, wherein the laplace transform of the delay circuit with the parameter is H (s, epsilon) ═ C (s (E + epsilon E)0)-A-Ade-sτ)-1B;
Computing a Taylor expansion of the Laplace transformed delay circuit with the parameters, wherein the Taylor expansion is
Using (s, epsilon) as parameter pair Performing Taylor expansion to obtainWherein, thereinIs the (i, j) th moment of the transfer function H (s, epsilon) and the matrix C is constant, thenCan be expressed as
in some embodiments, constructing a projection matrix from all moments of the transfer function further comprises:
Constructing a second-order Krylov subspace about the R matrix:
wherein i is 1,2, …, p0+1,j=1,2,…,q0+1;
When j is 1, an orthogonal basis matrix V of a first row matrix in the matrix R is calculated1:
when j is 2, …, q0When +1, calculate the orthogonal base matrix L formed by the jth row matrix in the matrix RjThe j-th row matrix satisfies the following second-order Krylov subspace relationship:
Orthogonalizing the jth row matrix by using a second-order Arnoldi algorithm to obtain an orthogonal base matrixFor matrixPartitioning:wherein
Calculating an orthogonal base matrix L of a jth row matrix in a matrix Rj:
colspan{Q2}=colspan{Lj}
Using the first j-1 row orthogonal matrix Vj-1And LjConstruction of a New matrix [ V ]j-1 Lj]To matrix [ V ]j-1 Lj]Orthogonalizing to obtain VjAccording to j 2, …, q0+1 calculation matrix R front q0Orthogonal basis matrix of +1 row matrix, Vj=orth([Vj-1Lj])。
In some embodiments, projecting the delay circuit with parameters in the projection matrix to obtain a macro model, further comprises:
using projection matrices VjReducing the order of the delay circuit with parameters to obtain a macro model:
In the embodiment of the invention, in order to solve the problem of complexity of simulation analysis in a high-speed integrated circuit, a novel model order reduction method for reducing the calculated amount under the condition of keeping a certain calculation precision is provided.
Firstly, establishing a mathematical model of an original large-scale delay circuit, wherein the mathematical model is a delay differential algebraic system:
in the formula (1), the matrix E is formed by Rn×nIs a singular matrix formed by inductance L and capacitance C, the matrix A belongs to Rn×nIs composed of resistance, conductance values and correlation matrix between devices, matrix Ad∈Rn×nIs formed by the values of the parasitic parameters in the interconnect lines. B is belonged to Rn×1Is an input matrix, u (t) is an input, C ∈ R1×nIs an output matrix; unknown variable x (t) e Rn×1Is the node voltage or branch current in the circuit. n is the order of the circuit system, namely the unknown quantity number of the whole system equation (1), and tau is the signal transmission delay.
Because E is a singular matrix, the solution of the differential algebraic system, especially the differential algebraic system with high index, is time-consuming, so that the conversion of the differential algebraic system into the ordinary differential system is an effective method for solving the problem. The invention uses epsilon embedding technology to convert the differential algebraic delay system into a differential delay system with parameters, and uses a model order reduction method to reduce the order of the differential algebraic delay system, and the reduced system can be well approximate to an original large-scale differential algebraic system.
Considering a differential algebraic delay system as in equation (1), where the coefficient matrix E is a singular matrix, which is converted to a nonsingular matrix using an epsilon embedding technique, a delay system of the form (2) with parameters epsilon can be obtained, where the values of the parameters epsilon satisfy epsilon >0 and epsilon < 1.
Here matrix E0Is composed ofn1+n2=n,Is of size n2The unit matrix of (2). Laplace transformation of system (2)
H(s,ε)=C(s(E+εE0)-A-Ade-sτ)-1B (3)
And (3) referring to a transfer function of the system (2), wherein the transfer function is an expression form of the system (2) in a frequency domain, next, extracting important information of the system (2) from the formula (3), and performing model reduction on the parameter system (2) so as to extract a macro model of the original circuit system.
Taylor expansion is carried out on the formula (3), a Taylor expansion coefficient is called as a moment of the system (3), the moment contains important information of the original system (2), and as long as the macro model contains most moment information of the original system, input and output information and basic properties of the original large-scale circuit system can be well reflected. The taylor expansion of the exponential term in equation (3) is:
taking the first three items as e-sτApproximation of
Bringing formula (4) into (3) yields:
taylor expansion is performed on equation (5) using (s, ε) as a parameter:
whereinIs the (i, j) th moment of the transfer function H (s, epsilon) and the matrix C is constant, thenCan be expressed as
When i >0, j >0
After all moments in the transfer function (6) are obtained, the original system (2) is reduced by a moment matching model reduction technique. Moment matching technique toConstructing a projection matrix V based on the moment information, and further constructing a projection matrix V based on the moment informationThe original system (2) is projected in a V-spanned space, whereby a reduced-order system of the original system (2), i.e. a macro-model of the original circuitry, is obtained. The procedure for solving the projection matrix V is given below, anda matrix written in the form:
constructing a second-order Krylov subspace about the R matrix:
wherein i is 1,2, …, p0+1,j=1,2,…,q0+1
Obtaining an orthogonal basis matrix V1:
when j is 2,3, … q0When +1, the space spanned by the j-th row is marked as:
In the formula (10), the compound represented by the formula (10),
(10) the formula can be constructed as a second order Krylov subspace
I.e. the orthogonal basis matrix L formed by the jth row matrixjThe j-th row matrix satisfies the following second-order Krylov subspace relationship:
(11) in the formula
The present invention uses a second order Arnoldi methodBy orthogonalizing equation (11), an orthogonal basis matrix can be obtainedFor orthogonal matrixIs divided into blocks to obtain
In formula (12):
according to the formulas (11) and (8), the compounds are obtained
colspan{Q2}=colspan{Lj}
So far, an orthogonal matrix L can be obtainedjAccording to the obtained orthogonal matrix Vj-1And LjConstruction of a New matrix [ V ]j-1Lj]To matrix [ V ]j-1 Lj]Orthogonalizing to obtain Vj
Vj=orth([Vj-1 Lj])
VjIs a space base matrix formed by the matrix of the first j rows in the formula (9). Sequentially calculating in a recursion way to obtain an orthogonal matrix V formed by all the matrixes in the formula (9), and obtaining a reduced-order matrix by taking V as a system (2)
In the formula (13)
The system (13) is a reduced order system of the system (2), wherein the value of r is according to p in the formula (9)0And q is0And (4) determining.
The system (13) can well approximate the original delay differential algebraic system (1), and the system (13) is called a macro model of the original large-scale delay circuit, and the macro model can keep important properties of the input-output relation, passivity and the like of the original circuit.
In the RLC circuit, a newly proposed algorithm and a reference algorithm reduce the original system from 607 th order to 36 th order, and a macro model is an efficient approximate system of the original large-scale system, so that an error exists in the approximate process. The invention compares the output of the obtained macro model with the output of the original large-scale system, and calculates the relative error of the output of the macro model, and the relative error calculation mode is as follows:
as shown in FIG. 2, y is the output of the original system, yrIs the output of the macro model. It can be observed that the relative error between the output of the macro model obtained by the present invention and the original system output is controlled at 10-5And the higher accuracy is achieved.
The scheme provided by the invention solves the problem of simulation analysis complexity in a high-speed integrated circuit by reducing the order of the model, and can reduce the calculation amount under the condition of keeping a certain calculation precision. Therefore, model reduction is carried out on the large-scale differential algebraic delay system, and a differential algebraic system with a smaller scale is searched to approximate the original large-scale system, so that the main properties of the original large-scale system are reserved.
Based on the same inventive concept, according to another aspect of the present invention, an embodiment of the present invention further provides a delay circuit macro model extraction system 400, as shown in fig. 3, including:
an obtaining module 401 configured to obtain a delay circuit of a macro model to be extracted;
a parameter embedding module 402 configured to embed the delay circuit into a preset parameter to obtain a delay circuit with a parameter;
a calculation module 403 configured to calculate a transfer function of the delay circuit with the parameters and to calculate all moments of the transfer function;
a construction module 404 configured to construct a projection matrix from all moments of the transfer function;
a projection module 405 configured to project the delay circuit with parameters in the projection matrix to obtain a macro model.
Based on the same inventive concept, according to another aspect of the present invention, as shown in fig. 4, an embodiment of the present invention further provides a computer apparatus 501, including:
at least one processor 520; and
the memory 510, the memory 510 stores a computer program 511 that is executable on the processor, and the processor 520 executes the program to perform any of the steps of the method for extracting a macro model of a delay circuit as described above.
Based on the same inventive concept, according to another aspect of the present invention, as shown in fig. 5, an embodiment of the present invention further provides a computer-readable storage medium 601, where the computer-readable storage medium 601 stores computer program instructions 610, and the computer program instructions 610, when executed by a processor, perform the steps of any one of the methods for extracting a delay circuit macro model as above.
Finally, it should be noted that, as will be understood by those skilled in the art, all or part of the processes of the methods of the above embodiments may be implemented by a computer program, which may be stored in a computer-readable storage medium, and when executed, may include the processes of the embodiments of the methods described above.
Further, it should be appreciated that the computer-readable storage media (e.g., memory) herein can be either volatile memory or nonvolatile memory, or can include both volatile and nonvolatile memory.
Those of skill would further appreciate that the various illustrative logical blocks, modules, circuits, and algorithm steps described in connection with the disclosure herein may be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate this interchangeability of hardware and software, various illustrative components, blocks, modules, circuits, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as software or hardware depends upon the particular application and design constraints imposed on the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the disclosed embodiments of the present invention.
The foregoing is an exemplary embodiment of the present disclosure, but it should be noted that various changes and modifications could be made herein without departing from the scope of the present disclosure as defined by the appended claims. The functions, steps and/or actions of the method claims in accordance with the disclosed embodiments described herein need not be performed in any particular order. Furthermore, although elements of the disclosed embodiments of the invention may be described or claimed in the singular, the plural is contemplated unless limitation to the singular is explicitly stated.
It should be understood that, as used herein, the singular forms "a", "an" and "the" are intended to include the plural forms as well, unless the context clearly supports the exception. It should also be understood that "and/or" as used herein is meant to include any and all possible combinations of one or more of the associated listed items.
The numbers of the embodiments disclosed in the embodiments of the present invention are merely for description, and do not represent the merits of the embodiments.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, and the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.
Those of ordinary skill in the art will understand that: the discussion of any embodiment above is meant to be exemplary only, and is not intended to intimate that the scope of the disclosure, including the claims, of embodiments of the invention is limited to these examples; within the idea of an embodiment of the invention, also technical features in the above embodiment or in different embodiments may be combined and there are many other variations of the different aspects of the embodiments of the invention as described above, which are not provided in detail for the sake of brevity. Therefore, any omissions, modifications, substitutions, improvements, and the like that may be made without departing from the spirit and principles of the embodiments of the present invention are intended to be included within the scope of the embodiments of the present invention.
Claims (10)
1. A method for extracting a macro model of a delay circuit is characterized by comprising the following steps:
obtaining a delay circuit of a macro model to be extracted;
embedding the delay circuit into preset parameters to obtain a delay circuit with parameters;
calculating a transfer function of the delay circuit with the parameters and calculating all moments of the transfer function;
constructing a projection matrix according to all moments of the transfer function;
projecting the delay circuit with parameters in the projection matrix to obtain a macro model.
2. The method of claim 1, wherein in obtaining the delay circuit of the macro model to be extracted, the delay circuit is:
wherein the matrix E ∈ Rn×nIs a singular matrix formed by inductance L and capacitance C, the matrix A belongs to Rn×nIs composed of resistance, conductance values and correlation matrix between devices, matrix Ad∈Rn×nFormed by values of parasitic parameters in the interconnection lines, B∈Rn×1Is an input matrix, u (t) is an input, C ∈ R1×nFor the output matrix, the unknown variable x (t) e Rn×1N is the order of the circuit system, namely the unknown quantity in the whole delay circuit, and tau is the signal transmission delay.
4. The method of claim 3, wherein calculating a transfer function of the delay circuit with the parameters and calculating all moments of the transfer function, further comprises:
calculating a laplace transform of the delay circuit with the parameter, wherein the laplace transform of the delay circuit with the parameter is H (s, epsilon) ═ C (s (E + epsilon E)0)-A-Ade-sτ)-1B;
Computing a Taylor expansion of the Laplace transformed delay circuit with the parameters, wherein the Taylor expansion is
Using (s, epsilon) as parameter pairPerforming Taylor expansion to obtainWherein, thereinIs the (i, j) th moment of the transfer function H (s, epsilon) and the matrix C is constant, thenCan be expressed as
6. the method of claim 5, wherein constructing a projection matrix from all moments of the transfer function further comprises:
Constructing a second-order Krylov subspace about the R matrix:
wherein i is 1,2, …, p0+1,j=1,2,…,q0+1;
When j is 1, an orthogonal basis matrix V of a first row matrix in the matrix R is calculated1:
when j is 2, …, q0When +1, calculate the orthogonal base matrix L formed by the jth row matrix in the matrix RjThe j-th row matrix satisfies the following second-order Krylov subspace relationship:
wherein i is 1,2, …, p0+1;
Orthogonalizing the jth row matrix by using a second-order Amoldi algorithm to obtain an orthogonal base matrixFor matrixPartitioning:whereinCalculating an orthogonal base matrix L of a jth row matrix in a matrix Rj:
colspan{Q2}=colspan{Lj}
Using the first j-1 row orthogonal matrix Vj-1And LjConstruction of a New matrix [ V ]j-1 Lj]To matrix [ V ]j-1 Lj]Orthogonalizing to obtain VjAccording to j 2, …, q0+1 calculation matrix R front q0Orthogonal basis matrix of +1 row matrix, Vj=orth([Vj-1 Lj])。
7. The method of claim 6, wherein projecting the delay circuit with parameters in the projection matrix to obtain a macro model, further comprises:
using projection matrices VjReducing the order of the delay circuit with parameters to obtain a macro model:
8. A delay circuit macro model extraction system, comprising:
the acquisition module is configured to acquire a delay circuit of the macro model to be extracted;
the parameter embedding module is configured to embed the delay circuit into preset parameters to obtain the delay circuit with parameters;
a calculation module configured to calculate a transfer function of the delay circuit with the parameters and to calculate all moments of the transfer function;
a construction module configured to construct a projection matrix from all moments of the transfer function;
a projection module configured to project the delay circuit with parameters in the projection matrix to obtain a macro model.
9. A computer device, comprising:
at least one processor; and
memory storing a computer program operable on the processor, wherein the processor executes the program to perform the steps of the method according to any of claims 1-7.
10. A computer-readable storage medium, in which a computer program is stored which, when being executed by a processor, is adapted to carry out the steps of the method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111161438.4A CN114004191A (en) | 2021-09-30 | 2021-09-30 | Method, system, device and medium for extracting delay circuit macro model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111161438.4A CN114004191A (en) | 2021-09-30 | 2021-09-30 | Method, system, device and medium for extracting delay circuit macro model |
Publications (1)
Publication Number | Publication Date |
---|---|
CN114004191A true CN114004191A (en) | 2022-02-01 |
Family
ID=79922179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111161438.4A Pending CN114004191A (en) | 2021-09-30 | 2021-09-30 | Method, system, device and medium for extracting delay circuit macro model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114004191A (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115496019A (en) * | 2022-11-17 | 2022-12-20 | 浪潮电子信息产业股份有限公司 | Simulation analysis method, system, equipment and storage medium of circuit |
WO2023179379A1 (en) * | 2022-03-24 | 2023-09-28 | 苏州浪潮智能科技有限公司 | Simulation method and system for nonlinear delay circuit system, and medium |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070255538A1 (en) * | 2006-04-27 | 2007-11-01 | Chang Gung University | Method of developing an analogical VLSI macro model in a global Arnoldi algorithm |
CN112861457A (en) * | 2021-02-10 | 2021-05-28 | 山东英信计算机技术有限公司 | Model order reduction method, device and medium for delay circuit system |
-
2021
- 2021-09-30 CN CN202111161438.4A patent/CN114004191A/en active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20070255538A1 (en) * | 2006-04-27 | 2007-11-01 | Chang Gung University | Method of developing an analogical VLSI macro model in a global Arnoldi algorithm |
CN112861457A (en) * | 2021-02-10 | 2021-05-28 | 山东英信计算机技术有限公司 | Model order reduction method, device and medium for delay circuit system |
Non-Patent Citations (1)
Title |
---|
苑伟政等: "《MEMS集成设计技术及应用》", 30 November 2014, pages: 5 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2023179379A1 (en) * | 2022-03-24 | 2023-09-28 | 苏州浪潮智能科技有限公司 | Simulation method and system for nonlinear delay circuit system, and medium |
CN115496019A (en) * | 2022-11-17 | 2022-12-20 | 浪潮电子信息产业股份有限公司 | Simulation analysis method, system, equipment and storage medium of circuit |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Vrudhula et al. | Hermite polynomial based interconnect analysis in the presence of process variations | |
Ratzlaff et al. | RICE: Rapid interconnect circuit evaluator | |
Raghavan et al. | AWE-inspired | |
CN114004191A (en) | Method, system, device and medium for extracting delay circuit macro model | |
Shook et al. | MLParest: Machine learning based parasitic estimation for custom circuit design | |
Lee et al. | HiPRIME: Hierarchical and passivity preserved interconnect macromodeling engine for RLKC power delivery | |
US20050160387A1 (en) | Methods to generate state space models by closed forms and transfer functions by recursive algorithms for RLC interconnect and transmission line and their model reduction and simulations | |
Shi | Graph-pair decision diagram construction for topological symbolic circuit analysis | |
Heydari et al. | Model reduction of variable-geometry interconnects using variational spectrally-weighted balanced truncation | |
TWI575394B (en) | Characterizing cell using input waveform geneartion considering different circuit topoloiges | |
Li et al. | Model order reduction of parameterized interconnect networks via a two-directional Arnoldi process | |
CN110456188A (en) | The stability of power system detection system and method for sparse polynomial chaos expansion | |
CN116151324A (en) | RC interconnection delay prediction method based on graph neural network | |
CN102915385B (en) | A kind of Interconnect model order reducing method based on time domain trapezoidal method difference | |
Aridhi et al. | Towards improving simulation of analog circuits using model order reduction | |
CN114169113A (en) | Transmission line time domain equivalent macro model generation method based on delay extraction | |
CN103678738B (en) | Interconnect model order reducing method based on time domain Multi-step Integration | |
Qi et al. | Wideband passive multiport model order reduction and realization of RLCM circuits | |
WO2022120078A1 (en) | Modelling timing behavior using augmented sensitivity data for physical parameters | |
Voss et al. | Model order reduction for nonlinear differential algebraic equations in circuit simulation | |
Guerra et al. | Approximate symbolic analysis of hierarchically decomposed analog circuits | |
US9946824B2 (en) | Efficient Ceff model for gate output slew computation in early synthesis | |
US9009636B2 (en) | Analog circuit simulator and analog circuit verification method | |
US11568117B2 (en) | Generating simulation-friendly compact physical models for passive structures | |
Hao et al. | Symbolic moment computation for statistical analysis of large interconnect networks |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |