66results about How to "Solving Convergence Problems" patented technology

The invention provides an underwater autonomous robot fixed depth control method based on reinforcement learning and belongs to the underwater robot control field. The method comprises steps that firstly, a Markoff decision process model for underwater autonomous robot fixed depth control is constructed, and expressions of a state variable, a control variable, a transfer model and a one-step lossfunction for underwater autonomous robot fixed depth control are respectively acquired; the decision network and the evaluation network are respectively established; through reinforcement learning, the decision network and the evaluation network can be continuously updated whenever an underwater autonomous robot moves forwards for each step during fixed depth control training till convergence; thefinal decision network for fixed depth control is acquired. The method is advantaged in that fixed depth control on the underwater autonomous robot under the condition of a completely-unknown underwater autonomous robot dynamic model is realized, and the practical value is high.

The invention discloses a hybrid-domain full wave form inversion method of central processing unit (CPU) / graphics processing unit (GPU) synergetic parallel computing. Compared with a traditional method, the method can be adopted to conduct the CPU / GPU synergetic parallel computing so as to obviously improve computing efficiency. The method enables a forward part of full wave form inversion to be placed to a time domain, namely forward is conducted in the time domain, the forward is converted to be conducted as the inversion in a frequency domain by the discrete fourier transform (DFT) being utilized, namely the DFT is adopted to pick wave field components corresponding to inversion frequency, and the inversion is conducted in the frequency domain from low frequency to high frequency. The method effectively resolves the problem of astringency of a standard time domain method, and avoids the problem that internal memory of wave form inversion of a standard three-dimensional frequency domain and a Laplace domain cannot be satisfied. The method is few in step, reduces computing expenses, and due to the fact the method can be adopted to conduct the CPU / GPU synergetic parallel computing, largely improves speed-up ratio of the computing method.

The invention discloses a virtualization network boundary data flow gathering method and apparatus and relates to the technical field of information safety. The apparatus disclosed by the invention comprises: a safety strategy module for maintaining and managing network boundary safety strategies of different tenants, a virtualization network data flow capturing module for monitoring a network data flow on a virtual switch and capturing the data packets of network interfaces specified by the safety strategies, and a network boundary data flow filter screen module for establishing, according to the network boundary safety strategies of the different tenants, filter screens based on network data flow filtering for the tenants respectively, filtering the captured network data packets through the filter screens of the corresponding tenants, encapsulating the network data flows which are obtained after filtering and belong to the safety domain boundaries of the tenants, and sending corresponding network safety products to the tenants. The invention also discloses the virtualization network boundary data flow gathering method. The technical scheme of the application effectively solves problems of gathering virtualization network safety domain boundary network data flows in a multi-tenant environment.

The invention discloses a rectangular coordinate Newton method load flow calculation method with a changeable Jacobian matrix. The method includes the following steps that original data input and voltage initialization are conducted; a node admittance matrix is formed; power and voltage deviations are calculated, and the maximum amount of unbalance delta Wmax is obtained; the Jacobian matrix J is formed; a correction equation is solved and a real part e and an imaginary part f of voltage are corrected; node and circuit data are output. According to the rectangular coordinate Newton method load flow calculation method, a Jacobian matrix calculation method different from that used in the following iteration processes is adopted in the initial iteration process, and the convergence problem of rectangular coordinate Newton method load flow calculation in analyzing a system with a small-impedance branch circuit is solved. When misconvergence happens with conventional rectangular coordinate Newton method load flow calculation, the rectangular coordinate Newton method load flow calculation method with the changeable Jacobian matrix can achieve reliable convergence, and the number of iterations is fewer compared with the prior art. The rectangular coordinate Newton method load flow calculation method with the changeable Jacobian matrix can effectively solve the convergence problem of the conventional rectangular coordinate Newton method load flow calculation in analyzing the system with the small-impedance circuit branch, and meanwhile load flow calculation can be performed on normal systems, and adverse effects are avoided.

The invention discloses a load flow calculation method of a rectangular coordinate newton method. The load flow calculation method comprises the following steps: inputting original data and initializing a voltage; forming a node admittance matrix; calculating power and a voltage deviation; forming a jacobian matrix J; solving a modified equation and modifying a real part e and an imaginary part f of the voltage; and outputting nodes and branch data. The calculation formula of partial elements (i=j) of the jacobian matrix obtained in a manner of load flow calculation by utilizing the rectangular coordinate newton method is modified and improved, so that the problem of astringency of the load flow calculation of the rectangular coordinate newton method during the analysis of branch systems containing small impedance is solved. As compared with the fact that the astringency of the small impedance branch systems cannot be realized by adopting the load flow calculation of the conventional rectangular coordinate newton method, the reliable astringency of the small impedance branch systems can be realized by utilizing the load flow calculation disclosed by the invention; and the load flow calculation to normal systems can be carried out at the same time, so that harmful effects are avoided.

The invention discloses a crowd calculation quality control method and device. The crowd calculation quality control method includes partitioning crowdsourcing task set to be processed into a plurality of subtask sets; selecting test tasks from the subtask sets and determining difficulty coefficient and accurate answer of the test tasks; determining the first subtask set and the second subtask set; determining the difficulty coefficient of subtasks in the first subtask set according to the difficulty coefficient of the test tasks; assessing ability coefficient of workers according to the test tasks, determining subtask completion accuracy of the workers according to the ability coefficient of the workers and the subtasks in the first subtask set and acquiring completion results of the subtasks in the second subtask set. By the crowd calculation quality control method, completion accuracy in crowdsourcing tasks is improved and cost is controlled to some extent.

The invention discloses a rapid decomposing method load flow calculating method suitable for a branch system comprising small impedance. The method comprises the step of inputting original data and initializing voltage, the step of forming a node admittance matrix, the step of forming a coefficient matrix B' and a coefficient matrix B'' of a modified equation and decomposing a factor table, the step of setting iteration count and a convergence mark, the step of calculating the inactive power unbalance amount delta Q, the step of judging whether the absolute value delta Qmax is smaller than the convergence accuracy epsilon or not, the step of calculating the active power unbalance amount delta P, the step of judging whether the absolute value delta Pmax is smaller than the convergence accuracy epsilon or not, the step of calculating the balance node power and the inactive power of a PV node, and the step of calculating the branch power. Part of elements in the coefficient matrix B' of the modified equation of a rapid decomposing method are modified, and the convergence problem produced during analyzing the branch system comprising the small impedance with the resistance not equal to zero in rapid decomposing load flow calculation is solved. When a convention method is adopted to calculate non-convergence, the rapid decomposing method load flow calculating method can carry out convergence reliably.

The invention discloses a rectangular coordinate Newton method load flow calculation method for a small-impedance-branch-containing power grid based on a compensation method. Before performing load flow calculation, a small-impedance branch is processed firstly; the small-impedance branch is changed into two series branches, wherein the impedance of one branch is z1 which is equal to r<ij>+j(x<ij>+x<c>) while the impedance of the other branch is z2 which is equal to -jx<c>, wherein x<c> is a normal electric reactance average value x<av> of an electric power system. One small-impedance branch is changed into two branches with relatively high impedance through series compensation, so that the convergence of load flow calculation can be improved. By adoption of the load flow calculation method, the small-impedance branch of any impedance can be calculated, or even a branch with zero impedance also can be calculated; the small-impedance branch data is processed before the load flow calculation is performed; and the data is preprocessed, so that a load flow calculation program is not changed, so that the rectangular coordinate Newton method load flow calculation method is particularly suitable for transforming an old load flow calculation program.

The invention discloses a rectangular coordinate newton method power flow calculation method based on an improved jacobian matrix; a jacobian matrix calculation method which is different between the first time of iteration and the rest time of iteration process is adopted, so that the problem of convergence of the rectangular coordinate newton method power flow calculation in a small-impedance branch power system is solved. During the first iteration, a PQ node adopts a current vector real part and virtual part obtained by an injected power given value to calculate jacobian elements; a PV nodeadopts a current vector real part and virtual part obtained by an active power given value and a reactive power given value to calculate jacobian elements, without using a power supply reactive powergiven value, so that the adverse effect to the power flow calculation convergence caused by the PV node power supply reactive given value can be avoided. When the conventional rectangular coordinatenewton method is adopted for calculation, the flow calculation is not converged; and compared with the prior art, the method disclosed by the invention has the advantages that the convergence is reliable, and the number of iterations is less than that of the prior art. The method can also carry out power flow calculation on a normal electric power system without adverse effects.

The invention discloses a fast decoupled power flow calculation method for serially compensating a power grid containing small impedance branches. Before power flow calculation is carried out, the small impedance branches of which resistance values are greater than reactance absolute values are first processed, and each small impedance branch is changed into a form of serially connecting the two branches, wherein impedance of one branch is that z1 equals to the sum of r, and j (x+xav), the impedance of the other branch is that z2 equals to -jxav, and the xav is a normal reactance mean of a power system. One small impedance branch is changed into two branches with greater impedance by serial compensation, so that the astringency of power flow calculation can be improved; the small impedance branches with any impedance can be calculated. According to the fast decoupled power flow calculation method disclosed by the invention, the small impedance branch with r being greater than an absolute value of x is only compensated, and the number of newly-added nodes and branches is effectively reduced; compared with an existing patent method, the fast decoupled power flow calculation method has the advantages that memory occupancy is less and the calculation speed is higher; data are preprocessed and power flow calculation procedures do need to be changed.

The invention discloses a Newton method load flow calculation method based on a node type and a modified Jacobian matrix, and solves a problem of convergence when rectangular coordinate Newton methodload flow calculation is used for solving a branch power system comprising a small impedance through employing a Jacobian matrix calculation method difference from each iteration process in the firstiteration. In the first iteration process, a PQ node employs real and virtual parts of a current phasor calculated through a node injection power given value for calculating Jacobian elements, and a PV node does not use the real and virtual parts of the current phasor for calculating the Jacobian elements, thereby avoiding the adverse impact on the load flow calculation convergence from a reactivepower given value of the PV node. When the rectangular coordinate Newton method load flow calculation is not convergent, the method can achieve the reliable convergence. Moreover, compared with the prior art, the number of iteration times is reduced. The method can perform the load flow calculation of a normal power system, and caused no adverse impact.

The invention discloses a wind resistance design method of a horizontal axis wind turbine blade based on strong wind action. The wind resistance design method comprises the following steps: selecting a blade model by utilizing Profili software; importing the data into Excel, and processing the data; importing the three-dimensional data into SW software; drawing an overall model of the wind turbine; carrying out grid division on the overall computational domain; carrying out grid independence verification; specifically setting numerical simulation calculation; and setting specific boundary conditions of simulation. Based on the strong wind load effect, the wind pressure, the speed, the vorticity and the speed distribution condition of a wake area of a wind turbine blade are involved, the page partitioning problem of a rotating domain and a static domain involved in the three-dimensional modeling and numerical simulation process of a wind turbine is solved, locally-encrypted unstructured grid division is adopted, when the calculation efficiency is improved, the influence on the calculation precision is small, the working time is shortened, and the calculation cost is saved; and a reference effect is provided for understanding the actual condition of the wind turbine under the action of the strong wind load and reasonably arranging the wind power plant.

The invention discloses a method for calculating cooperative power flow of an active power distribution network and a power transmission network based on a system of systems (SoS). The method comprises the steps of defining an electric power system to be the SoS, and building an SoS model; modeling the power transmission network and the active power distribution network to be two independent systems in the whole system; considering the relation between a distributed power supply and a load in the active power distribution network, and determining an operating state of the active power distribution network; respectively considering power flow calculation of the active power distribution system and the power transmission network system according to different operating states of the active power distribution network; taking boundary calculation results of the two systems as initial values of each iteration thereof until the boundary results converge, and acquiring a power flow calculation result of the whole electric power system. According to the invention, the SoS model for power flow calculation of the electric power system is built, the power transmission network and the active power distribution network are modeled to be two independent systems, and voltage and power of only boundary nodes are required to be interacted, so that complex power flow calculation of the electric power system is enabled to be simplified. In addition, distributed management is easy to be performed, and a certain basis is laid for optimization of a distributed electric power system.

The invention relates to a thermoelectric interface interaction method for integrated energy network dynamic hybrid simulation, and belongs to the technical field of integrated energy network digitalsimulation. The characteristics of different sub-network width time scales of the comprehensive energy network are fully considered, and the accuracy, effectiveness and convergence of comprehensive energy network hybrid simulation interface interaction of different calculation step lengths are achieved by utilizing a post-prediction correction method. Hybrid simulation electromechanical transientand thermal dynamic two-side thermoelectric interface interaction with different calculation step lengths can be effectively processed, effective interaction of thermoelectricity in hybrid simulationis ensured, and the wide time scale characteristic under system faults and disturbance is embodied. The method is easy to implement, is a practical and effective method for thermoelectric interface interaction of thermoelectric coupling equipment of an actual multi-energy system, is small in calculation amount, and is suitable for real-time and super-real-time dynamic simulation of an integrated energy network system.

The invention relates to a method for carrying out optimized speed classification for a digital integrated circuit with a transparent latch under the effect of process deviation, which belongs to the field of integrated circuits and comprises the following steps; computing the minimum clock cycle accumulative density distribution function of the digital integrated circuit with the transparent latch and determining an optimal clock cycle classification demarcation point and an optimal test sequence of the demarcation point according to the minimum clock cycle accumulative density distribution function to maximize the total income. The method can obtain the work clock cycle distribution of the transparent latch circuit through a random configuration method based on very low computation complexity and very high solution accuracy to avoid the convergence problem in a random arrival time solution; the optimal test sequence of the cycle classification demarcation point can be determined through an optimization method of which computation complexity only is 0 (n log n) so as to minimize the test cost; and when the position of the cycle classification demarcation point is determined by a greedy algorithm to maximize the circuit design income, the optimality of each iterative computation can be ensured in theory.

The invention discloses a power flow calculation method for a small impedance branch-containing power network employing a compensation method and a fast decomposition method. Before power flow calculation, a small impedance branch is firstly processed and the small impedance branch is changed into a form of series connection of two branches, wherein the impedance of one branch is z1=r+j(x+xc), the impedance of the other branch is z2=-jxc, and xc is a mean value of reactance absolute values of the branches of which the reactance absolute values are smaller than a reactance threshold in a power system, and is called a normal reactance mean value of the power system. One small impedance branch is changed into two branches with relatively high impedance through series compensation, so that the convergence of power flow calculation can be improved. The small impedance branch with any impedance or even the branch with zero impedance can be calculated. According to the power flow calculation method, the problem of convergence due to the fact that a small impedance branch-containing branch system is analyzed by a conventional power flow calculation method employing the fast decomposition method can be effectively solved, and meanwhile, power flow calculation can be carried out on a normal system. Therefore, adverse effects are avoided.

The invention discloses a small-impedance branch PQ endpoint variable Jacobian matrix newton-method power flow calculation method, which comprises the following steps: determining the type T of a branch connected by two-end nodes according to branch resistance and reactance; for a small-impedance branch PQ endpoint and other nodes, adopting different methods to calculate a Jacobian matrix; and solving a correction equation. In the iteration process, the small-impedance branch PQ endpoint adopts a Jacobian matrix calculating method different from those of other nodes, thereby solving the problem of convergence of rectangular coordinate newton power flow calculation when analyzing an electrical power system having a small-impedance branch. If non-convergence occurs when a conventional rectangular coordinate newton power flow calculation method is adopted, the method in the invention can realize reliable convergence and is smaller in iteration times than existing patent technology; and since the method not only can solve the problem of convergence of conventional rectangular coordinate newton power flow calculation when analyzing the electrical power system having the small-impedance branch, but also can carry out power flow calculation on a normal power system, no bad influence is caused.

The invention discloses a method for applying standards in an electromagnetic field edge element method. The method includes the steps that a three-dimensional model of an object to be analyzed is established, and the three-dimensional model is subjected to mesh generation; element matrixes corresponding to a magnetic vector potential double-rotation item,a coulomb standard item and an electric current density item of edge elements obtained through mesh generation are calculated, and integration is performed so that an overall finite element equation can be formed; corresponding boundary conditions are applied to the overall finite element equation; the overall finite element equation to which the corresponding boundary conditions are applied is solved, a unique magnetic vector potential basic solution is obtained, post-processing is performed on the obtained magnetic vector potential basic solution, and a corresponding electromagnetic field quantity solution is obtained. According to the method, on the premise that accuracy is guaranteed, the magnetic vector potential uniqueness problem in the edge element method is effectively solved; the establishment problem of a complex generating tree in a traditional method is avoided; meanwhile, the efficiency of the method is higher than that of a coulomb standard applying method adopted in the edge elements through a Lagrange multiplier method.

The invention discloses a hybrid-domain full wave form inversion method of central processing unit (CPU) / graphics processing unit (GPU) synergetic parallel computing. Compared with a traditional method, the method can be adopted to conduct the CPU / GPU synergetic parallel computing so as to obviously improve computing efficiency. The method enables a forward part of full wave form inversion to be placed to a time domain, namely forward is conducted in the time domain, the forward is converted to be conducted as the inversion in a frequency domain by the discrete fourier transform (DFT) being utilized, namely the DFT is adopted to pick wave field components corresponding to inversion frequency, and the inversion is conducted in the frequency domain from low frequency to high frequency. The method effectively resolves the problem of astringency of a standard time domain method, and avoids the problem that internal memory of wave form inversion of a standard three-dimensional frequency domain and a Laplace domain cannot be satisfied. The method is few in step, reduces computing expenses, and due to the fact the method can be adopted to conduct the CPU / GPU synergetic parallel computing, largely improves speed-up ratio of the computing method.

The invention discloses a radio frequency amplifier circuit design method and a radio frequency amplifier circuit, and the method comprises the steps: obtaining a bias condition when a radio frequency amplifier works in a specified class, and controlling the radio frequency amplifier to be biased to the specified class; obtaining a radio frequency small signal S parameter of the radio frequency amplifier under a specified class through testing; obtaining an input equivalent circuit based on the radio frequency small signal S parameter; the input equivalent circuit comprises a first inductor, a first resistor and a first capacitor which are connected in series; the load-pull impedance of the radio frequency amplifier under different frequencies is obtained through testing; based on the load-pull impedance, an output equivalent circuit is obtained; the output equivalent circuit comprises a second inductor, a second resistor and a second capacitor; the second resistor and the second capacitor are connected in parallel and then are connected in series with the second inductor. According to the design method, based on the application characteristics of the amplifier, the circuit design process is optimized, the convergence problem is solved, and the amplifier design efficiency is improved.

The invention discloses a method for pre-evaluating the line loss of an electric power system. The method comprises the following steps of: collecting original data; preparing data needed by direct current power flows; carrying out direct current power flow calculation to obtain the power flow of each branch; calculating a line loss power of each branch according to a line loss calculation model provided by the method; and calculating the line loss of the system. The method provided by the invention is used for calculating on the basis of a direct current power flow model, the method is simple, the equation is completely linearized and the solving speed is rapid. The method is sufficiently exact for pre-evaluation of the line loss of a new project of a power grid, especially exact for pre-evaluation of the line loss of a new project of a main network with high voltage grade and high ratio of line reactance to resistance (X / R).

The invention discloses a load flow calculation method with a small-impedance power grid serial-connection compensation right-angle coordinate Newton method. Before load flow calculation, small-impedance branches except those of which nodes at two ends are both PQ nodes and a branch transformation ratio k is 1 are firstly processed, then the small-impedance branches can be transformed into a mode that two branches are serially connected, the impedance of one branch is z1=r+j(x+xav), the impedance of the other branch is z2=-jxav, and xav is a normal reactance mean value of a power system. Through serial-connection compensation, small impedance can be transformed into two branches with relatively large impedance, so that the astringency of load flow calculation can be improved. Meanwhile, as serial-connection compensation on small-impedance branches of which nodes at two ends are both PQ nodes and the branch transformation ratio k is 1 are not implemented, the number of new nodes and branches can be effectively reduced, furthermore the memory space can be reduced, and the calculation speed can be increased. The load flow calculation method is applicable to transformation of conventional load flow calculation procedures.

The invention discloses a multicast routing optimization method based on multiple-upstream-interface IGMP / MLD (Internet Group Management Protocol / Multicast Listener Discovery Protocol) proxy. The method comprises the following steps: the multiple-upstream-interface IGMP / MLD proxy has a plurality of upstream interfaces; for a multicast group (S, G) sent by each multicast source, any source multicast group (*, G) subscribed by each receiver, or a specific source multicast group (S, G), the quantity of upstream interface in using is one and only one, and the upstream interface in using is called a practical upstream interface; an accessible upstream interface selected according to a routing table and policy routing is called a default upstream interface; and * indicates any multicast source including S. The method effectively solves the problem of tunnel convergence, improves the mobile multicast routing efficiency, and realizes the performance optimization.

The invention discloses an improved constrained multi-objective optimization problem solving method. The method comprises the following steps: 1, processing equality constraint conditions and inequality constraint conditions of a constrained multi-objective optimization problem by adopting a constraint processing method based on a constraint violation index; 2, using an improved multi-target bacterial population taxis algorithm based on Pareto dominance to solve and obtain a Pareto optimal solution set and a corresponding Pareto front of a multi-target optimization problem; 3, selecting an optimal compromise solution for constraining the multi-objective optimization problem from the Pareto optimal solution set obtained by solving by adopting a multi-objective decision method based on the relationship between the objective satisfaction and the objective weight, According to the invention, the problems that constraint processing is complex, algorithm convergence, solution set diversity and time efficiency are poor, and preference of decision makers cannot be met are solved, and the requirement for solving the constraint multi-objective optimization problem is met.

The invention relates to a multi-objective optimization method and system for regenerative type electric heating to participate in wind power accommodation. A multi-objective optimization model for regenerative type electric heating to participate in wind power accommodation is constructed, wherein the multi-objective optimization model comprises a maximum wind power accommodation objective function, a minimum regenerative type electric heating operation cost objective function and a minimum system carbon emission objective function, and the problems that existing regenerative type electric heating participates in wind power accommodation research and mainly focuses on benefit evaluation and economy, and the objective functions are single are solved; and an NSGA-III algorithm is used for solving the multi-objective model for regenerative type electric heating to participate in wind power accommodation, and in the aspect of the algorithm, when a high-dimensional multi-objective optimization problem is processed, the problems that a traditional algorithm is insufficient in convergence and solving efficiency and the like are solved. According to the method and the system, the optimization efficiency and the wind power accommodation capability are improved, the carbon emission is effectively reduced, and the equipment operation cost is reduced.

The invention discloses a crowd calculation quality control method and device. The crowd calculation quality control method includes partitioning crowdsourcing task set to be processed into a plurality of subtask sets; selecting test tasks from the subtask sets and determining difficulty coefficient and accurate answer of the test tasks; determining the first subtask set and the second subtask set; determining the difficulty coefficient of subtasks in the first subtask set according to the difficulty coefficient of the test tasks; assessing ability coefficient of workers according to the test tasks, determining subtask completion accuracy of the workers according to the ability coefficient of the workers and the subtasks in the first subtask set and acquiring completion results of the subtasks in the second subtask set. By the crowd calculation quality control method, completion accuracy in crowdsourcing tasks is improved and cost is controlled to some extent.

The invention discloses a power flow calculation method for a small impedance branch-containing power network employing a compensation method and a fast decomposition method. Before power flow calculation, a small impedance branch is firstly processed and the small impedance branch is changed into a form of series connection of two branches, wherein the impedance of one branch is z1=r+j(x+xc), the impedance of the other branch is z2=-jxc, and xc is a mean value of reactance absolute values of the branches of which the reactance absolute values are smaller than a reactance threshold in a power system, and is called a normal reactance mean value of the power system. One small impedance branch is changed into two branches with relatively high impedance through series compensation, so that the convergence of power flow calculation can be improved. The small impedance branch with any impedance or even the branch with zero impedance can be calculated. According to the power flow calculation method, the problem of convergence due to the fact that a small impedance branch-containing branch system is analyzed by a conventional power flow calculation method employing the fast decomposition method can be effectively solved, and meanwhile, power flow calculation can be carried out on a normal system. Therefore, adverse effects are avoided.

The invention discloses a rapid decomposing method load flow calculating method suitable for a branch system comprising small impedance. The method comprises the step of inputting original data and initializing voltage, the step of forming a node admittance matrix, the step of forming a coefficient matrix B' and a coefficient matrix B'' of a modified equation and decomposing a factor table, the step of setting iteration count and a convergence mark, the step of calculating the inactive power unbalance amount delta Q, the step of judging whether the absolute value delta Qmax is smaller than the convergence accuracy epsilon or not, the step of calculating the active power unbalance amount delta P, the step of judging whether the absolute value delta Pmax is smaller than the convergence accuracy epsilon or not, the step of calculating the balance node power and the inactive power of a PV node, and the step of calculating the branch power. Part of elements in the coefficient matrix B' of the modified equation of a rapid decomposing method are modified, and the convergence problem produced during analyzing the branch system comprising the small impedance with the resistance not equal to zero in rapid decomposing load flow calculation is solved. When a convention method is adopted to calculate non-convergence, the rapid decomposing method load flow calculating method can carry out convergence reliably.