Circuit stable state simulation calculation method and device, medium and program product

Abstract

The invention provides a circuit stable state simulation calculation method and device, a medium and a program product. The method comprises the steps that a circuit equation and a steady-state equation based on target variables are constructed according to a circuit structure and the Kirchhoff's law; performing pre-simulation through time domain transient simulation to obtain a first initial value of the target variable during iteration; performing iteration processing on the circuit equation through a Newton iteration method in a time domain to obtain a first iteration equation, and obtaining a calculation value of the target variable at each sampling point according to the first initial value; reducing the value of the target variable at each sampling point in the period into the selection of a first initial value through a targeting method; performing iteration processing on the steady-state equation through a Newton iteration method to obtain a target initial value of a target variable in a convergence state; and substituting the target initial value into the first iterative equation to obtain a stable periodic solution of the circuit. According to the scheme, the stable state of the circuit can be quickly calculated without waiting for a system to reach a stable operation state, and the compatibility of circuit nonlinearity can be improved.

Description

Technical Field

[0001] The present invention relates to the technical field of circuit stable state calculation, and in particular to a circuit stable state simulation calculation method, equipment, medium and program product. Background Art

[0002] Modern society has increasingly higher requirements for the stability of signal reception in communication technology devices, including smartphones, computers, and other wireless communication devices. Therefore, it is particularly important to test the stability of integrated circuits faster and more accurately to determine the state of these devices in stable mode. Existing integrated circuit stability testing methods usually use standard continuous time domain transient simulation until the state of the device reaches stability. However, this method requires a lot of time to wait for the system to reach a stable operating state from startup when applied to communication equipment, and there is no guarantee that results will be obtained. On the other hand, although there are methods to calculate the circuit stability state from the frequency domain, this method is extremely poorly applicable to nonlinear circuits, and today's circuits are becoming more and more complex and more nonlinear. Therefore, a direct method to quickly calculate the circuit stability state in the time domain is necessary. Summary of the invention

[0003] The purpose of the present invention is to provide a circuit stable state simulation calculation method, device, medium and program product, which can quickly calculate the circuit stable state without spending a lot of time waiting for the system to reach a stable operating state from startup, and directly perform calculations in the time domain, thereby increasing compatibility with nonlinear circuits.

[0004] The technical solution provided by the present invention is as follows: In a first aspect, the present application provides a circuit stable state simulation calculation method for stability testing of an integrated circuit, comprising the steps of: Constructing a circuit equation based on a target variable and a steady-state equation satisfied by the target variable when the circuit to be calculated reaches a steady state according to the circuit structure of the circuit to be calculated and Kirchhoff's law; Pre-simulating the circuit to be calculated by time-domain transient simulation to obtain a first initial value of the target variable during iteration; Iteratively processing the circuit equation in the time domain by Newton iteration method to obtain a first iterative equation within a cycle, and obtaining a calculated value of the target variable at each sampling point within the cycle according to the first initial value; The value of the target variable at each sampling point in the period is reduced to the selection of the first initial value by a shooting method, so as to obtain a derivative relationship between the value of the target variable at each sampling point and the first initial value; Iteratively processing the steady-state equation by Newton iteration method, and obtaining the target initial value of the target variable when in a convergent state according to the calculated value of the target variable at each sampling point and the derivative relationship between the value of the target variable at each sampling point and the first initial value; Substituting the target initial value into the first iterative equation, a stable periodic solution of the circuit to be calculated is obtained.

[0005] In some implementations, the pre-simulating the circuit to be calculated by time-domain transient simulation to obtain the first initial value of the target variable during iteration specifically includes: Performing transient simulation on the first preset number of cycles of the circuit to be calculated, and determining whether the circuit to be calculated reaches a steady state; If it is determined that the circuit to be calculated reaches a steady state, the simulation is stopped, and a stable periodic solution of the circuit to be calculated is calculated; If it is determined that the circuit to be calculated has not reached a steady state, further determining whether the simulation waveform has a periodic convergence trend; If it is determined that the simulation waveform does not have a periodic convergence trend, the simulation is stopped, and the simulation value of the target variable at the last moment is selected as the first initial value; If it is determined that the simulation waveform has a periodic convergence trend, the preset simulation period is extended. If the circuit to be calculated reaches a steady state after the preset simulation period, the simulation is stopped. If the simulation waveform still does not converge after the preset simulation period, the simulation value of the target variable at the last moment is selected as the first initial value.

[0006] In some implementations, the iterative processing of the circuit equation in the time domain by Newton's iteration method to obtain a first iterative equation within a period specifically includes: Divide the simulation cycle into a number of evenly distributed sampling points; Differentiating the circuit equation using an implicit Euler method to obtain a differential equation; The differential equation is iteratively processed by Newton's iteration method to obtain the first iterative equation.

[0007] In some embodiments, it further comprises: The first iterative equation in the period is solved by LU decomposition to obtain the calculated value of the target variable at each sampling point and several intermediate calculated quantities, where the intermediate calculated quantities are the derivative values ​​of the target variable and the related parameters having a functional relationship with the target variable.

[0008] In some embodiments, the method of reducing the value of the target variable at each sampling point within a period to the selection of the first initial value by the targeting method specifically includes: Deriving the first initial value on both sides of the differential equation and expanding it by the chain rule to obtain a derivative equation related to the first initial value; Substitute the intermediate calculated value into the derivative equation to obtain the derivative relationship between the value of the target variable at each sampling point and the first initial value.

[0009] In some embodiments, the iterative processing of the steady-state equation by Newton's iteration method comprises: Iterate the steady-state equation by Newton iteration method to obtain a second iterative equation based on the Jacobian matrix of the target variable; The Jacobian matrix is ​​solved by the derivative relationship between the value of the target variable at each sampling point and the first initial value, and then the second iterative equation is solved to obtain the iterative value of the first initial value.

[0010] In some embodiments, it further comprises: Substitute the iterative value of the first initial value as the updated first initial value into the second iterative equation for iteration until the first initial value converges to obtain the target initial value.

[0011] In a second aspect, the present application provides a computer device, comprising a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of a circuit stable state simulation calculation method described in the first aspect.

[0012] In a third aspect, the present application provides a computer storage medium having a computer program or instructions stored thereon, which, when executed by a processor, implements the steps of a circuit stable state simulation calculation method described in the first aspect.

[0013] In a fourth aspect, the present application provides a computer program product, including a computer program or instructions, which, when executed by a processor, implements the steps of a circuit stable state simulation calculation method described in the first aspect.

[0014] The circuit stable state simulation calculation method, device, medium and program product provided by the present invention can quickly calculate the circuit stable state without spending a lot of time waiting for the system to reach a stable operating state from startup, and the calculation is performed directly in the time domain, which increases the compatibility with nonlinear circuits. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The preferred implementation scheme will be described below in a clear and understandable manner with reference to the accompanying drawings to further illustrate the above-mentioned characteristics, technical features, advantages and implementation methods of the present solution.

[0016] Figure 1It is a schematic diagram of the overall process of an embodiment of the present invention; Figure 2 is a flow chart of an embodiment of the present invention; Figure 3 It is a flow chart of an embodiment of the present invention. DETAILED DESCRIPTION

[0017] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the specific implementation methods of the present invention will be described below with reference to the accompanying drawings. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings and other implementation methods can be obtained based on these drawings without creative work.

[0018] In order to simplify the drawings, only the parts related to the present invention are schematically shown in each figure, and they do not represent the actual structure of the product. In addition, in order to simplify the drawings and facilitate understanding, in some figures, only one of the parts with the same structure or function is schematically drawn or marked. In this article, "one" not only means "only one", but also means "more than one".

[0019] Modern society has increasingly higher requirements for the stability of signal reception in communication technology devices, including smartphones, computers, and other wireless communication devices. Therefore, it is particularly important to test the stability of integrated circuits faster and more accurately to determine the state of these devices in stable mode. Existing integrated circuit stability testing methods usually use standard continuous time domain transient simulation until the device state reaches stability.

[0020] However, this method needs to spend a lot of time to wait for the system to reach a stable operating state from startup when applied to communication equipment, and there is no guarantee that the result will be obtained. On the other hand, although there are methods to calculate the circuit stability state from the frequency domain, this method is extremely poorly applicable to nonlinear circuits, and today's circuits are becoming more and more complex and more nonlinear. Therefore, a method that can directly and quickly calculate the circuit stability state in the time domain is necessary.

[0021] Based on Newton's shooting method, this application designs a new solution idea to solve the iterative equation, and reduces the calculation problem of the value of the target variable at each sampling point in the iteration process to the selection problem of the initial value of the target variable, so as to calculate the working condition of the circuit in a stable state. On the one hand, this method avoids the large amount of time required for the traditional continuous time-domain transient simulation for the device to reach a stable state from the startup state, greatly reducing the waiting time; on the other hand, this method is a calculation directly performed in the time domain, which increases the compatibility with nonlinear circuits. At the same time, this application adds the optimization of the adaptive selection of initial values ​​during initial value iteration, as well as the optimization of the new matrix solution method in the initial value iteration calculation process, which reduces the number of calculations and iterations, and can further accelerate convergence and calculation speed. Below, this scheme will be described in detail with reference to the accompanying drawings: In one embodiment, the reference specification Figure 1 The present application provides a circuit stable state simulation calculation method for stability testing of integrated circuits, comprising the steps of: S100, constructing a circuit equation based on the target variable and a steady-state equation satisfied by the target variable when the circuit to be calculated reaches a steady state according to the circuit structure of the circuit to be calculated and Kirchhoff's law; S200, pre-simulating the circuit to be calculated by time-domain transient simulation to obtain a first initial value of the target variable during iteration; S300, iteratively processing the circuit equation in the time domain by Newton iteration method to obtain a first iterative equation within a cycle, and obtaining a calculated value of the target variable at each sampling point within the cycle according to the first initial value; S400, reducing the value of the target variable at each sampling point in the period to the selection of the first initial value by a shooting method, and obtaining a derivative relationship between the value of the target variable at each sampling point and the first initial value; S500, iteratively processing the steady-state equation by Newton iteration method, and obtaining the target initial value of the target variable when in a convergent state according to the calculated value of the target variable at each sampling point and the derivative relationship between the value of the target variable at each sampling point and the first initial value; S600 , substituting the target initial value into the first iterative equation to obtain a stable periodic solution of the circuit to be calculated.

[0022] This application is used for stability testing of integrated circuits or communication equipment. According to the circuit structure of the integrated circuit or communication equipment, the circuit equation and steady-state equation of the target variable can be constructed. For example, taking the node voltage of the integrated circuit as the target variable, according to the circuit structure of the integrated circuit and Kirchhoff's theorem, the circuit equation can be constructed: (1) in, For external input source, Node voltage, is the voltage of the relevant node, is the change in charge of the relevant node. , and There is a corresponding functional relationship. According to the circuit structure, the corresponding functional equation can be listed, so it can usually be , Recorded as , to indicate the relationship between them.

[0023] In a specific implementation, in order to highlight the nonlinearity of the circuit, the two can be taken as the simplest nonlinear function ,Right now , .Pick All of the above variables change with time. To find the stable state of the circuit is to find the stable solution of this equation with respect to period T that satisfies: (2) That is, the target variable is in a convergent state.

[0024] Before performing simulation calculation, this solution first pre-simulates the circuit to be calculated through time domain transient simulation, and uses the calculated value of the target variable obtained by pre-simulation as the first initial value in subsequent iterations. After obtaining the circuit equation of the target variable, the circuit equation is iteratively processed in the time domain through Newton iteration method to obtain the first iterative equation within the cycle, and calculation is performed based on the first initial value to obtain the calculated value of the target variable at each sampling point within the cycle.

[0025] Considering that the calculated value of the target variable at each sampling point is directly affected by the first iterative equation and the first initial value, this solution reduces the value of the target variable at each sampling point within the period to the selection problem of the first initial value through the shooting method, obtains the derivative relationship between the value of the target variable at each sampling point and the first initial value, and iterates the steady-state equation through the Newton iteration method to obtain the second iterative equation. The target initial value of the target variable when in the convergence state can be obtained according to the calculated value of the target variable at each sampling point and the derivative relationship between the value of the target variable at each sampling point and the first initial value. Substituting the target initial value into the first iterative equation for calculation, the stable periodic solution of the circuit to be calculated can be obtained. This calculation method has a faster calculation speed and can quickly calculate the stable state of the circuit without spending a lot of time waiting for the system to reach a stable operating state from the startup state. In addition, this calculation method is directly calculated in the time domain, which can increase compatibility with nonlinear circuits. In one embodiment, the reference specification Figure 2Based on the above-mentioned embodiment, the circuit to be calculated is pre-simulated by time domain transient simulation to obtain the first initial value of the target variable during iteration, which specifically includes: S210, performing transient simulation on the first preset number of cycles of the circuit to be calculated, and determining whether the circuit to be calculated reaches a steady state; S220, if it is determined that the circuit to be calculated reaches a steady state, the simulation is stopped, and a stable periodic solution of the circuit to be calculated is calculated; S230, if it is determined that the circuit to be calculated has not reached a steady state, further determining whether the simulation waveform has a periodic convergence trend; S240, if it is determined that the simulation waveform does not have a periodic convergence trend, the simulation is stopped, and the simulation value of the target variable at the last moment is selected as the first initial value; S250. If it is determined that the simulation waveform has a periodic convergence trend, the preset simulation period is extended. If the calculation circuit reaches a steady state after the preset simulation period, the simulation is stopped. If the simulation waveform still does not converge after the preset simulation period, the simulation value of the target variable at the last moment is selected as the first initial value.

[0026] When performing pre-simulation, transient simulation is performed for the first preset number of cycles of the circuit to be calculated. The preset number can be selected according to actual needs and scenarios, and this application does not limit it. For example, transient simulation is performed for the first three cycles of the circuit to be calculated. During pre-simulation, if the circuit to be calculated directly reaches a steady state, the simulation is stopped directly, and the stable periodic solution of the circuit to be calculated is directly calculated. If it is judged that the circuit to be calculated has not reached a steady state, in order to avoid complicating the calculation, it can be further judged whether the simulation waveform has a periodic convergence trend. If there is no periodic convergence trend, it means that the circuit will not reach a steady state in a short time, stop the simulation, and select the simulation value of the target variable at the last moment as the first initial value; if it is judged that the simulation waveform has a periodic convergence trend, it means that the surface may reach a steady state in a short time in the future, and the preset simulation cycle can be extended. If the circuit to be calculated reaches a steady state after the preset simulation cycle, the simulation is stopped. If the simulation waveform still does not converge after the preset simulation cycle, the simulation value of the target variable at the last moment is selected as the first initial value.

[0027] In one embodiment, the reference specification Figure 3 Based on the above-mentioned embodiment, the circuit equation is iteratively processed in the time domain by Newton iteration method to obtain the first iterative equation within the cycle, which specifically includes: S310, dividing the simulation period into a number of evenly distributed sampling points; S320, performing differentiation processing on the circuit equation by using an implicit Euler method to obtain a differential equation; S330, iteratively processing the differential equation by Newton's iteration method to obtain a first iterative equation; S340, solving the first iterative equation in the cycle through LU decomposition to obtain the calculated value of the target variable at each sampling point and several intermediate calculation quantities, where the intermediate calculation quantities are the derivative values ​​of the target variable and the related parameters having a functional relationship with the target variable.

[0028] The implicit Euler method, also known as the backward Euler method, is a method of numerical solution based on an implicit formula. The implicit formula cannot be solved directly. Generally, the initial value needs to be obtained using the explicit Euler formula, and then the implicit Euler formula is used for iterative solution. Therefore, the implicit formula is more complex than the explicit formula, but has better stability.

[0029] In a specific implementation, the period T is divided into M uniform sampling points, denoted as . , At the same time, remember for The calculated value in the iteration, for The estimated value in the iteration. is the time difference between two sampling points. Substituting the above parameters into formula (1), Using implicit Euler method for numerical differentiation, we can get the differential equation: (3) Use Newton's iteration method for the above differential equation, that is, , At the same time, for the convenience of notation, , , so the above formula (3) becomes: (4) In this case, due to , ,but , Therefore, formula (4) can be expressed as: .

[0030] By iterating formula (3), we can get Calculated The value of The values ​​of each sampling point in the entire cycle. At the same time, for each m, the intermediate calculation amount can also be obtained and These two intermediate quantities will be used in subsequent calculations, so it is necessary to store C and G obtained in each step. In this solution, in order to facilitate storage, when solving equation (4), you can choose to use the LU decomposition method to solve it, so that you can get The LU decomposition form is easy to store and saves storage space. At the same time, since the dimension of the matrix does not change with the number of iterations, the information about the matrix structure only needs to be stored once.

[0031] In one embodiment, based on the above embodiment, the value of the target variable at each sampling point in the period is reduced to the selection of the first initial value by a shooting method, which specifically includes: Derivative the first initial value on both sides of the differential equation and expand it by the chain rule to obtain the derivative equation related to the first initial value; Substitute the intermediate calculated value into the derivative equation to obtain the derivative relationship between the value of the target variable at each sampling point and the first initial value.

[0032] As can be seen from the previous text, The value over the entire cycle can be attributed to the initial value , so we can Considered as function, so that Find the derivative. Taking the derivative on both sides of equation (3) we get:

[0033] Using the chain rule, expand the above formula and replace the previous notation , Substituting the backshift term, we can get: (5) Through formula (5), we can get and The calculation relationship between them can be obtained by continuous calculation. , and the initial In addition, the intermediate quantity in formula (5) and It has been obtained in the previous article and can be directly quoted and used. At the same time, since the structure stored in this method is In the form of LU decomposition, this structure can be directly used to solve (5), saving the time of performing matrix decomposition again. In this example, formula (5) is: .

[0034] In one embodiment, based on the above-mentioned embodiment, the steady-state equation is iteratively processed by Newton iteration method, including: Iterate the steady-state equation by Newton iteration method to obtain the second iterative equation based on the Jacobian matrix of the target variable; The Jacobian matrix is ​​solved by the derivative relationship between the value of the target variable at each sampling point and the first initial value, and then the second iterative equation is solved to obtain the iterative value of the first initial value.

[0035] Specifically, as explained above Can be considered as The function of . Then the periodicity condition (2) can be expressed as: Using Newton's iteration method, we can get: (6) in, is the identity matrix, yes The Jacobian matrix is: , which is what we mentioned earlier , substituting into equation (6) can complete one iteration. In the general Newton shooting method, the Gaussian elimination method is required to solve equation (6). This method is extremely complex in multi-node circuits, and the speed is very slow and the convergence is also very poor. In this scheme, according to the definition of derivative, the design uses directional derivative to replace In the calculation When using This avoids the use of Gaussian elimination to solve the matrix equation.

[0036] According to the above method, one iteration of the initial value can be completed, that is, from calculate . use Repeat the above operation to perform another iteration until a certain , - If it is small enough, it can be considered that the initial value has converged, that is, the iterative value of the first initial value is substituted into the second iterative equation as the updated first initial value and iterated until the first initial value converges, and the target initial value can be obtained.

[0037] In addition, in order to prevent the occurrence of The value causes the matrix to be singular. ,in, is the inverse of the norm of the general residual, that is .

[0038] The target initial value is obtained by the above method , and then substitute the target initial value into the first iterative equation within the cycle to obtain the stable periodic solution of the circuit to be calculated This solution can quickly calculate the circuit stable state without spending a lot of time waiting for the system to reach a stable operating state from startup, and it is calculated directly in the time domain, which increases the compatibility with nonlinear circuits.

[0039] In one embodiment, the present application provides a computer device, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of a circuit stable state simulation calculation method of the aforementioned embodiment.

[0040] In one embodiment, the present application provides a computer storage medium having a computer program or instruction stored thereon, and when the computer program or instruction is executed by a processor, the steps of a circuit stable state simulation calculation method of the aforementioned embodiment are implemented.

[0041] In one embodiment, the present application provides a computer program product, including a computer program or instructions, which implements the steps of a circuit stable state simulation calculation method of the aforementioned embodiment when the computer program or instructions are executed by a processor.

[0042] It should be noted that the above embodiments can be freely combined as needed. The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the principle of the present invention, and these improvements and modifications should also be considered as the protection scope of the present invention.

Claims

1. A circuit stable state simulation calculation method, characterized in that: Used for stability testing of integrated circuits, including the following steps: Constructing a circuit equation based on a target variable and a steady-state equation satisfied by the target variable when the circuit to be calculated reaches a steady state according to the circuit structure of the circuit to be calculated and Kirchhoff's law; Pre-simulating the circuit to be calculated by time-domain transient simulation to obtain a first initial value of the target variable during iteration; Iteratively processing the circuit equation in the time domain by Newton iteration method to obtain a first iterative equation within a cycle, and obtaining a calculated value of the target variable at each sampling point within the cycle according to the first initial value; The value of the target variable at each sampling point in the period is reduced to the selection of the first initial value by a shooting method, so as to obtain a derivative relationship between the value of the target variable at each sampling point and the first initial value; Iteratively processing the steady-state equation by Newton iteration method, and obtaining the target initial value of the target variable when in a convergent state according to the calculated value of the target variable at each sampling point and the derivative relationship between the value of the target variable at each sampling point and the first initial value; Substituting the target initial value into the first iterative equation, a stable periodic solution of the circuit to be calculated is obtained.

2. A circuit stable state simulation calculation method according to claim 1, characterized in that: The pre-simulation of the circuit to be calculated by time-domain transient simulation to obtain the first initial value of the target variable during iteration specifically includes: Performing transient simulation on the first preset number of cycles of the circuit to be calculated, and determining whether the circuit to be calculated reaches a steady state; If it is determined that the circuit to be calculated reaches a steady state, the simulation is stopped, and a stable periodic solution of the circuit to be calculated is calculated; If it is determined that the circuit to be calculated has not reached a steady state, further determining whether the simulation waveform has a periodic convergence trend; If it is determined that the simulation waveform does not have a periodic convergence trend, the simulation is stopped, and the simulation value of the target variable at the last moment is selected as the first initial value; If it is determined that the simulation waveform has a periodic convergence trend, the preset simulation period is extended. If the circuit to be calculated reaches a steady state after the preset simulation period, the simulation is stopped. If the simulation waveform still does not converge after the preset simulation period, the simulation value of the target variable at the last moment is selected as the first initial value.

3. A circuit stable state simulation calculation method according to claim 1, characterized in that: The iterative processing of the circuit equation in the time domain by Newton's iteration method to obtain the first iterative equation within a period specifically includes: Divide the simulation cycle into a number of evenly distributed sampling points; Differentiating the circuit equation using an implicit Euler method to obtain a differential equation; The differential equation is iteratively processed by Newton's iteration method to obtain the first iterative equation.

4. A circuit stable state simulation calculation method according to claim 3, characterized in that: Also includes: The first iterative equation in the cycle is solved by LU decomposition to obtain the calculated value of the target variable at each sampling point and several intermediate calculation quantities, where the intermediate calculation quantities are the derivative values ​​of the target variable and the related parameters having a functional relationship with the target variable.

5. A circuit stable state simulation calculation method according to claim 4, characterized in that: The method of reducing the value of the target variable at each sampling point within the period to the first initial value by the shooting method specifically includes: Deriving the first initial value on both sides of the differential equation and expanding it by the chain rule to obtain a derivative equation related to the first initial value; Substitute the intermediate calculated value into the derivative equation to obtain the derivative relationship between the value of the target variable at each sampling point and the first initial value.

6. A circuit stable state simulation calculation method according to claim 1, characterized in that: The iterative processing of the steady-state equation by Newton's iteration method comprises: Iterate the steady-state equation by Newton iteration method to obtain a second iterative equation based on the Jacobian matrix of the target variable; The Jacobian matrix is ​​solved by the derivative relationship between the value of the target variable at each sampling point and the first initial value, and then the second iterative equation is solved to obtain the iterative value of the first initial value.

7. A circuit stable state simulation calculation method according to claim 6, characterized in that: Also includes: Substitute the iterative value of the first initial value as the updated first initial value into the second iterative equation for iteration until the first initial value converges to obtain the target initial value.

8. A computer device comprising a memory, a processor and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of a circuit stable state simulation calculation method according to any one of claims 1 to 7.

9. A computer storage medium having a computer program or instruction stored thereon, characterized in that: When the computer program or instruction is executed by a processor, the steps of a circuit stable state simulation calculation method according to any one of claims 1 to 7 are implemented.

10. A computer program product comprising a computer program or instructions, characterized in that When the computer program or instruction is executed by a processor, the steps of a circuit stable state simulation calculation method according to any one of claims 1 to 7 are implemented.

Images