A Circuit Steady-State Simulation Calculation Method, Device, Medium, and Program Product

Abstract

The present invention provides a method, device, medium and program product for simulating and calculating the stable state of a circuit. The method includes: constructing a circuit equation and a steady-state equation based on target variables according to the circuit structure and Kirchhoff's law; performing pre-simulation through time-domain transient simulation to obtain a first initial value of the target variable during iteration; performing iterative processing on the circuit equation by the Newton iteration method in the time domain to obtain a first iterative equation, and obtaining the calculated values of the target variable at each sampling point according to the first initial value; reducing the values of the target variable at each sampling point within a period to the selection of the first initial value by the shooting method; performing iterative processing on the steady-state equation by the Newton iteration method to obtain a target initial value of the target variable when in a converged state; substituting the target initial value into the first iterative equation to obtain a stable periodic solution of the circuit. This solution can quickly calculate the stable state of the circuit without waiting for the system to reach a stable operating state, and can increase the compatibility with circuit nonlinearity.

Description

Technical Field

[0001] The present invention relates to the technical field of circuit steady state calculation, and particularly to a circuit steady state simulation calculation method, device, medium and program product. Background Art

[0002] In modern society, there are increasingly high requirements for the stability of signal reception of communication technology devices, including smart phones, 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 the stable mode. The existing integrated circuit stability test methods usually use standard continuous time-domain transient simulation until the state of the device reaches stability. However, this method applied to communication devices requires a lot of time to wait for the system to reach the stable operating state from startup, and it cannot ensure that results will be obtained. On the other hand, although there are already methods in the frequency domain to calculate the circuit steady state, this method has extremely poor applicability to non-linear circuits, and today's circuits are becoming more and more complex and tend to be non-linear. Therefore, a direct method for quickly calculating the circuit steady state in the time domain is necessary. Summary of the Invention

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

[0004] The technical solution provided by the present invention is as follows:

[0005] In a first aspect, the present application provides a circuit steady state simulation calculation method for the stability test of integrated circuits, including the steps of:

[0006] Construct 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 the steady state according to the circuit structure of the circuit to be calculated and Kirchhoff's law;

[0007] Perform pre-simulation on the circuit to be calculated through time-domain transient simulation to obtain the first initial value of the target variable during iteration;

[0008] Perform iterative processing on the circuit equation by the Newton iteration method in the time domain to obtain the first iterative equation within a period, and obtain the calculated values of the target variable at each sampling point within the period according to the first initial value;

[0009] Reduce the values of the target variable at each sampling point within the period to the selection of the first initial value through the shooting method, and obtain the derivative relationship between the values of the target variable at each sampling point and the first initial value;

[0010] Iteratively process the steady-state equation by Newton's iterative method, and obtain the target initial value of the target variable in the converged state according to the calculated values of the target variable at each sampling point and the derivative relationship between the values of the target variable at each sampling point and the first initial value;

[0011] Substitute the target initial value into the first iterative equation to obtain the stable periodic solution of the circuit to be calculated.

[0012] In some embodiments, 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:

[0013] Perform transient simulation on the first preset number of cycles of the circuit to be calculated, and determine whether the circuit to be calculated reaches a steady state;

[0014] If it is determined that the circuit to be calculated reaches a steady state, stop the simulation and calculate the stable periodic solution of the circuit to be calculated;

[0015] If it is determined that the circuit to be calculated does not reach a steady state, further determine whether there is a periodic convergence trend in the simulation waveform;

[0016] If it is determined that there is no periodic convergence trend in the simulation waveform, stop the simulation and select the simulation value of the target variable at the last moment as the first initial value;

[0017] If it is determined that there is a periodic convergence trend in the simulation waveform, extend the preset simulation period. If the circuit to be calculated reaches a steady state after the preset simulation period, stop the simulation. If the simulation waveform still does not converge after the preset simulation period, select the simulation value of the target variable at the last moment as the first initial value.

[0018] In some embodiments, iteratively processing the circuit equation by Newton's iterative method in the time domain to obtain the first iterative equation within a period specifically includes:

[0019] Divide the simulation period into several uniformly distributed sampling points;

[0020] Perform differential processing on the circuit equation by the implicit Euler method to obtain a differential equation;

[0021] Perform iterative processing on the differential equation by Newton's iterative method to obtain the first iterative equation.

[0022] In some embodiments, it further includes:

[0023] Solve the first iterative equation within the period through LU decomposition to obtain the calculated values of the target variable at each sampling point and several intermediate calculated quantities, where the intermediate calculated quantities are the relevant parameters having a functional relationship with the target variable and the derivative values of the target variable.

[0024] In some embodiments, the selection of the first initial value by reducing the values of the target variable at each sampling point within the period through the shooting method specifically includes:

[0025] Take the derivative of both sides of the differential equation with respect to the first initial value and expand it through the chain rule of differentiation to obtain a derivative equation related to the first initial value;

[0026] Substitute the intermediate calculated quantities into the derivative equation to obtain the derivative relationship between the values of the target variable at each sampling point and the first initial value.

[0027] In some embodiments, the iterative processing of the steady-state equation through the Newton iteration method includes:

[0028] Iterate the steady-state equation through the Newton iteration method to obtain a second iterative equation based on the Jacobian matrix of the target variable;

[0029] Solve the Jacobian matrix through the derivative relationship between the values of the target variable at each sampling point and the first initial value, and then solve the second iterative equation to obtain the iterative value of the first initial value.

[0030] In some embodiments, it further includes:

[0031] Use the iterative value of the first initial value as the updated first initial value and substitute it into the second iterative equation for iteration until the first initial value converges to obtain the target initial value.

[0032] In a second aspect, the present application provides a computer device, including a memory, a processor, and a computer program stored on the memory, and the processor executes the computer program to implement the steps of a circuit steady-state simulation calculation method described in the first aspect.

[0033] In a third aspect, the present application provides a computer storage medium, on which a computer program or instruction is stored, and when the computer program or instruction is executed by a processor, it implements the steps of a circuit steady-state simulation calculation method described in the first aspect.

[0034] In a fourth aspect, the present application provides a computer program product, including a computer program or instruction, and when the computer program or instruction is executed by a processor, it implements the steps of a circuit steady-state simulation calculation method described in the first aspect.

[0035] A method, device, medium and program product for simulating and calculating the stable state of a circuit provided by the present invention can quickly calculate the stable state of the circuit without spending a large amount of time waiting for the system to reach a stable operating state from startup, and directly calculates in the time domain, increasing the compatibility with non-linear circuits. BRIEF DESCRIPTION OF THE DRAWINGS

[0036] The following will further illustrate the above-mentioned characteristics, technical features, advantages and their implementation manners of the present solution in a clear and understandable manner in combination with the drawings of the preferred embodiments.

[0037] Figure 1 is a schematic diagram of the overall process of an embodiment of the present invention;

[0038] Figure 2 is a schematic diagram of the process of an embodiment of the present invention;

[0039] Figure 3 is a schematic diagram of the process of an embodiment of the present invention. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0040] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the specific embodiments of the present invention will be described below with reference to the accompanying drawings. Obviously, the accompanying drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, without creative efforts, other drawings can be obtained according to these drawings and other embodiments can be obtained.

[0041] For the sake of simplicity of 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, for the sake of simplicity and easy understanding of the drawings, in some figures, components with the same structure or function are only schematically shown for one of them, or only one of them is marked. In this article, "one" not only means "only this one", but also means "more than one" situation.

[0042] In modern society, there are increasingly high requirements for the signal reception stability of communication technology devices, including smart phones, 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 the stable mode. The existing integrated circuit stability test methods usually use standard continuous time-domain transient simulation until the state of the device reaches stability.

[0043] However, applying this method to communication devices requires a large amount of time to wait for the system to reach a stable operating state from startup, and it cannot ensure obtaining results. On the other hand, although there are already methods in the frequency domain to calculate the stable state of circuits, this method has extremely poor applicability to non-linear circuits, and current circuits are becoming more and more complex and tend to be non-linear. Therefore, a method that can directly and quickly calculate the stable state of circuits in the time domain is necessary.

[0044] Based on the Newton shooting method, this application designs a new solution idea to solve the iterative equation, reduces the calculation problem of the values of the target variable at each sampling point during the iterative process to the selection problem of the initial value of the target variable, and thus calculates the operating conditions of the circuit in the stable state. On the one hand, this method avoids the large amount of time required for traditional continuous time-domain transient simulation for the device to reach the stable state from the startup state, greatly reducing the waiting time; on the other hand, this method performs calculations directly in the time domain, increasing the compatibility with non-linear circuits. At the same time, this application adds the optimization of the adaptive selection of the initial value during the initial value iteration, as well as the optimization of the new matrix solution method during the initial value iteration calculation process, reducing the number of calculations and the number of iterations, and being able to further accelerate the convergence and calculation speed. Next, this solution will be described in detail in conjunction with the accompanying drawings:

[0045] In one embodiment, referring to the attached Figure 1 of the specification, this application provides a method for simulating and calculating the stable state of a circuit for the stability test of an integrated circuit, including the steps:

[0046] S100. Construct 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;

[0047] S200. Perform pre-simulation on the circuit to be calculated through time-domain transient simulation to obtain the first initial value of the target variable during iteration;

[0048] S300. Iteratively process the circuit equation in the time domain through the Newton iteration method to obtain the first iterative equation within a period, and obtain the calculated values of the target variable at each sampling point within the period according to the first initial value;

[0049] S400. Reduce the values of the target variable at each sampling point within the period to the selection of the first initial value through the shooting method, and obtain the derivative relationship between the values of the target variable at each sampling point and the first initial value;

[0050] S500. Iteratively process the steady-state equation through the Newton iteration method, and obtain the target initial value of the target variable in the converged state according to the calculated values of the target variable at each sampling point and the derivative relationship between the values of the target variable at each sampling point and the first initial value;

[0051] S600. Substitute the target initial value into the first iteration equation to obtain the stable periodic solution of the circuit to be calculated.

[0052] This application is used for the stability test of integrated circuits or communication devices. According to the circuit structure of the integrated circuit or communication device, the circuit equation and the 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 as follows:

[0053] (1)

[0054] where is the external input source, is the node voltage, is the related node voltage, is the change in the charge of the related node. Since , and have a corresponding functional relationship, according to the circuit structure, the corresponding functional equation can be listed. Therefore, usually , are denoted as , to represent the relationship between them.

[0055] In a specific implementation, in order to highlight the nonlinearity of the circuit, the two can be taken as the simplest nonlinear function , that is , . Take . All the above variables change with time. Solving the stable state of the circuit is to solve the stable solution of this equation with respect to the period T that satisfies:

[0056] (2)

[0057] That is, the target variable is in a convergent state.

[0058] Before performing the simulation calculation in this solution, first perform a pre-simulation on the circuit to be calculated through a time-domain transient simulation, and use the calculated value of the target variable obtained from the pre-simulation as the first initial value for subsequent iteration. After obtaining the circuit equation of the target variable, perform iterative processing on the circuit equation through the Newton iteration method in the time domain to obtain the first iteration equation within the period, and calculate according to the first initial value to obtain the calculated values of the target variable at each sampling point within the period.

[0059] Considering that the calculated values of the target variable at each sampling point are directly affected by the first iteration equation and the first initial value, this solution reduces the values of the target variable at each sampling point within a period to the problem of selecting the first initial value through the shooting method, obtains the derivative relationship between the values of the target variable at each sampling point and the first initial value, and iteratively processes the steady-state equation through the Newton iteration method to obtain the second iteration equation, which can obtain the target initial value of the target variable in the convergent state based on the calculated values of the target variable at each sampling point and the derivative relationship between the values of the target variable at each sampling point and the first initial value. Substituting this target initial value into the first iteration equation for calculation can obtain the stable periodic solution of the circuit to be calculated. This calculation method is faster, can quickly calculate the stable state of the circuit without spending a lot of time waiting for the system to reach the stable operation state from the startup state, and this calculation method is directly calculated in the time domain, which can increase the compatibility with nonlinear circuits.

[0060] In one embodiment, referring to the attached drawings of the specification Figure 2 , on the basis of the foregoing embodiment, pre-simulate the circuit to be calculated through time-domain transient simulation to obtain the first initial value of the target variable during iteration, specifically including:

[0061] S210. Perform transient simulation on the first preset number of cycles of the circuit to be calculated, and judge whether the circuit to be calculated reaches a steady state;

[0062] S220. If it is judged that the circuit to be calculated reaches a steady state, stop the simulation and calculate the stable periodic solution of the circuit to be calculated;

[0063] S230. If it is judged that the circuit to be calculated does not reach a steady state, further judge whether there is a periodic convergence trend in the simulation waveform;

[0064] S240. If it is judged that there is no periodic convergence trend in the simulation waveform, stop the simulation and select the simulation value of the target variable at the last moment as the first initial value;

[0065] S250. If it is judged that there is a periodic convergence trend in the simulation waveform, extend the preset simulation period. If the circuit to be calculated reaches a steady state after the preset simulation period, stop the simulation. If the simulation waveform still does not converge after the preset simulation period, select the simulation value of the target variable at the last moment as the first initial value.

[0066] When performing pre-simulation, transient simulation is carried out on the first preset number of cycles of the circuit to be calculated. The preset number can be selected according to actual requirements and scenarios, and this application does not limit it. For example, transient simulation is carried out on the first three cycles of the circuit to be calculated. During pre-simulation, if the circuit to be calculated directly reaches the steady state, the simulation is directly stopped, and the steady-state periodic solution of the circuit to be calculated is directly calculated. If it is determined that the circuit to be calculated does not reach the steady state, in order to avoid complicating the calculation, it can be further determined whether there is a periodic convergence trend in the simulation waveform. If there is no periodic convergence trend, it means that the circuit will not reach the steady state in a short time, 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, it means that the steady state may be reached in a short time in the future, the preset simulation period can be extended. If the circuit to be calculated reaches the 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.

[0067] In one embodiment, referring to the attached specification Figure 3 , on the basis of the foregoing embodiment, the circuit equation is iteratively processed by the Newton iteration method in the time domain to obtain the first iterative equation within the period, specifically including:

[0068] S310. Divide the simulation period into a number of uniformly distributed sampling points;

[0069] S320. Differentiate the circuit equation by the implicit Euler method to obtain a differential equation;

[0070] S330. Iteratively process the differential equation by the Newton iteration method to obtain the first iterative equation;

[0071] S340. Solve the first iterative equation within the period by LU decomposition to obtain the calculated values of the target variable at each sampling point, and a number of intermediate calculated quantities. The intermediate calculated quantities are the relevant parameters with a functional relationship with the target variable and the derivative values of the target variable.

[0072] The implicit Euler method, also known as the backward Euler method, is a method for numerically solving according to an implicit formula. The implicit formula cannot be directly solved. Generally, the Euler explicit formula is used to obtain the initial value, and then the Euler implicit formula is used for iterative solution. Therefore, the implicit formula is more complex to calculate than the explicit formula, but has good stability.

[0073] In a specific implementation, the period T is divided into M uniformly sampled points, denoted as . , . At the same time, denote as the calculated value in the iteration, is the estimated value in the iteration. is the time difference between two sampling points. Substitute the above parameters into formula (1), using the implicit Euler method for numerical differentiation, the differential equation can be obtained:

[0074] (3)

[0075] Using the Newton iteration method for the above differential equation, that is ,

[0076] Meanwhile, for the convenience of notation, denote , , so the above formula (3) becomes:

[0077] (4)

[0078] In this example, since , , then , . Therefore, equation (4) can be expressed as:

[0079] .

[0080] By performing iterative calculations on equation (3), can be calculated from to obtain the value of, thus obtaining the values of each sampling point over the entire period. Meanwhile, for each m, the intermediate calculation quantities and can also be obtained. These two intermediate quantities will be used in subsequent calculations, so it is necessary to store C and G obtained at each step. In this solution, for the convenience of storage, when solving equation (4), the LU decomposition method can be selected for solution, so that the LU decomposition form can be obtained, which is convenient for storage 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.

[0081] In one embodiment, based on the foregoing embodiment, reducing the values of the target variable at each sampling point within the period to the selection of the first initial value by the shooting method specifically includes:

[0082] Taking the derivative of the first initial value on both sides of the differential equation and expanding it through the chain rule of differentiation to obtain the derivative equation related to the first initial value;

[0083] Substitute the intermediate calculation results into the derivative equation to obtain the derivative relationship between the values of the target variable at each sampling point and the first initial value.

[0084] As can be seen from the foregoing, the values over the entire period can be attributed to the selection of the initial value , so can be regarded as a function of, and thus the derivative of can be obtained. Differentiating both sides of Equation (3) with respect to gives:

[0085]

[0086] Using the chain rule of differentiation, expand the above equation, and at the same time substitute the notations , and then transpose the terms to obtain:

[0087] (5)

[0088] From Equation (5), the calculation relationship between and can be obtained, so that the final can be obtained through continuous calculations, and the initial . In addition, the intermediate quantities and in Equation (5) have been obtained previously and can be directly cited. At the same time, since the stored structure in this method is in the form of the LU decomposition of , this structure can be directly used to solve Equation (5), saving the time for matrix decomposition again. In this example, Equation (5) is: .

[0089] In one embodiment, based on the foregoing embodiment, the steady-state equation is iteratively processed by the Newton iteration method, including:

[0090] Iterate the steady-state equation by the Newton iteration method to obtain a second iteration equation based on the Jacobian matrix of the target variable;

[0091] Solve the Jacobian matrix through the derivative relationship between the values of the target variable at each sampling point and the first initial value, and then solve the second iteration equation to obtain the iterative value of the first initial value.

[0092] Specifically, as described above can be regarded as a function of, denoted as . Thus, the periodic condition (2) can be expressed as: . Using the Newton iteration method for this equation, we can obtain:

[0093] (6)

[0094] Wherein, is the identity matrix, is the Jacobian matrix of, that is:

[0095] , which is also the in the previous text. Substituting it into equation (6) can complete one iteration. In the general Newton shooting method, solving equation (6) requires the use of Gaussian elimination method, which has extremely high complexity, slow speed, and poor convergence in multi-node circuits. In this solution, according to the definition of the derivative, the directional derivative is designed to replace . When calculating , is used to replace it, thus avoiding the use of Gaussian elimination method to solve the matrix equation.

[0096] According to the above method, one iteration of the initial value can be completed, that is, calculating from . Using to repeat the above operation can perform another iteration until for a certain , - is small enough, it can be considered that the initial value has converged, that is, the iterative value of the first initial value is used as the updated first initial value and substituted into the second iterative equation for iteration until the first initial value converges, and the target initial value can be obtained.

[0097] In addition, in order to prevent the situation of matrix singularity caused by the value after iteration, take , wherein, is the reciprocal of the norm of the general residual, that is .

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

[0099] In one embodiment, the present application provides a computer device, including a memory, a processor, and a computer program stored on the memory. The processor executes the computer program to implement the steps of a method for simulating and calculating the stable state of a circuit in the foregoing embodiment.

[0100] In one embodiment, the present application provides a computer storage medium, on which a computer program or instruction is stored, and when the computer program or instruction is executed by a processor, the steps of a circuit steady-state simulation calculation method in the foregoing embodiment are implemented.

[0101] In one embodiment, the present application provides a computer program product, including a computer program or instruction, and when the computer program or instruction is executed by a processor, the steps of a circuit steady-state simulation calculation method in the foregoing embodiment are implemented.

[0102] It should be noted that the above embodiments can be freely combined as needed. The above are only the preferred embodiments of the present invention. It should be pointed out that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements should also be regarded 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