CN110148946A - A kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor - Google Patents

A kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor Download PDF

Info

Publication number
CN110148946A
CN110148946A CN201910441049.3A CN201910441049A CN110148946A CN 110148946 A CN110148946 A CN 110148946A CN 201910441049 A CN201910441049 A CN 201910441049A CN 110148946 A CN110148946 A CN 110148946A
Authority
CN
China
Prior art keywords
equation
confactor
formula
node
capacitance sensor
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201910441049.3A
Other languages
Chinese (zh)
Other versions
CN110148946B (en
Inventor
林立亨
董树锋
朱承治
唐坤杰
宋永华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Zhejiang University ZJU
State Grid Zhejiang Electric Power Co Ltd
Original Assignee
Zhejiang University ZJU
State Grid Zhejiang Electric Power Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Zhejiang University ZJU, State Grid Zhejiang Electric Power Co Ltd filed Critical Zhejiang University ZJU
Priority to CN201910441049.3A priority Critical patent/CN110148946B/en
Publication of CN110148946A publication Critical patent/CN110148946A/en
Application granted granted Critical
Publication of CN110148946B publication Critical patent/CN110148946B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/04Circuit arrangements for ac mains or ac distribution networks for connecting networks of the same frequency but supplied from different sources
    • H02J3/06Controlling transfer of power between connected networks; Controlling sharing of load between connected networks
    • H02J3/382
    • H02J3/383
    • H02J3/386
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/388Islanding, i.e. disconnection of local power supply from the network
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/50Photovoltaic [PV] energy
    • Y02E10/56Power conversion systems, e.g. maximum power point trackers
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00Energy generation through renewable energy sources
    • Y02E10/70Wind energy
    • Y02E10/76Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a kind of isolated island micro-capacitance sensor tidal current computing methods solved based on two step of confactor, this method comprises the following steps: rewriting to former power flow equation, all constant terms are moved on on the right side of equation, and retain PV node voltage magnitude, voltage term is replaced with system frequency ω by logarithmic form;Introduce confactor y, u, by Load flow calculation it is equations turned be relation function between one group of underdetermined equation, one group of overdetermined equation and one group of confactor, set trend initial value, the equation after above-mentioned conversion is solved by two-step method, wherein the first step by it is transformed it is equations turned find for least square problem closer to the linearisation point that really solves, improve algorithm to the robustness of initial value, second step solves following iteration step variable directly to reduce calculation amount.Method of the invention can effectively improve convergence and robustness.

Description

A kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor
Technical field
The present invention relates to technical field of power systems more particularly to isolated island micro-capacitance sensor Load flow calculation fields, propose one kind The isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor.
Background technique
Micro-capacitance sensor Load flow calculation is an important field of research as micro-capacitance sensor stability analysis, the basis distributed rationally. In micro-grid connection operation, Load flow calculation is similar to distribution power system load flow calculation.And the micro-capacitance sensor of isolated operation is controlled in equity Under system, balance nodes are not present in system, and there are the DG of sagging control, need to solve to system frequency, therefore traditional tide Flow calculation methodologies are no longer applicable in, and need to study the algorithm for being more suitable for isolated island micro-capacitance sensor Load flow calculation.
From the point of view of current isolated island micro-capacitance sensor Load flow calculation result of study, Part Methods are using the thought of optimization to trend side Cheng Jinhang is solved, such as BFGS Trust Region Algorithm, Levenberg-Marquardt (LM) algorithm, but such algorithm has ginseng Number is excessive, the problem for adjusting ginseng complicated, and there are tail portion effects for LM algorithm, it is difficult to adapt to the calculating of high-precision requirement.Another kind is thought Road is former Power Flow Problem to be decomposed into traditional Load flow calculation and sagging node updates two sub-problems, but convergence rate is slower.For This, it is necessary to propose a kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor, it is auxiliary by introducing two Help the factor, by former non-linear power flow equation be decomposed into one group owe constant linear equation, one group of determined linear equation and one group of auxiliary to Relation function between amount is iterated solution to transformed equation in two steps.Through example contrast verification, which has convergence The feature that speed is fast, strong robustness and calculating time are short.
Summary of the invention
The technical problem to be solved by the present invention is to which existing technical solution is improved and improved, one kind is provided and is based on The isolated island micro-capacitance sensor tidal current computing method that two step of confactor solves, to improve the convergence rate and robustness that calculate.
In order to solve the above technical problems, the present invention adopts the following technical scheme:
Include the following steps:
1, Load flow calculation equation is rewritten
Two step derivation algorithm of confactor there are certain requirements the form of power flow equation, need to carry out former power flow equation It rewrites.
Isolated island micro-capacitance sensor power flow equation can abstract representation are as follows:
The x ∈ of F (x)=0 Rn
If n node is shared in system, wherein there is nPQA PQ node, nPVA PV node and nDA sagging node, then formula (1) Contain n+n for onePQ+nDThe n+n of a location variablePQ+nDTie up equation;
By moving on to constant terms all in above formula on the right side of equation, and retain PV node voltage magnitude, be rewritten as Lower form:
F1(x)=p
In formula:
p0For in power flow equation constant term form vector, and be subsequent convenience of calculation, the electricity of PV node is remained in p Pressure amplitude value indicates by its quadratic form,
For vector x, will wherein voltage term be replaced with system frequency ω by logarithmic form:
X=[δ2,…,δn,a1,…,an,w]T
In formula: δiFor node i phase angle, ai=lnVi=2lnUi, w=ln ω.
2, confactor is introduced
Confactor y, u are introduced,
In formula: item U second from the bottomD,iFor sagging node voltage amplitude, n is sharedD?;α, β are active and reactive load factor, For the node i at each branch both ends, j enables Kij=UiUj cosδij、Lij=UiUj sinδij, wherein δijij
In formula: aij=ai+aj, δijij
By Load flow calculation it is equations turned be relationship letter between one group of underdetermined equation, one group of overdetermined equation and one group of confactor Number, it may be assumed that
In formula: nonlinear transformation f () reversible between confactor, inverse transformation y=f-1(u) concrete form is as follows:
3, trend initial value design
Set flow state initial guess x0, the number of iterations k=0 is taken, convergence precision ε, maximum number of iterations k are setmax, System frequency ω0, voltage magnitude U0With phase angle initial value δ0, take y0=f-1(Cx0)。
4, two-step method solves equation
The equation after above-mentioned conversion is solved by two-step method, wherein the first step will be transformed equations turned for most Small two multiply problem to find closer to the linearisation point really solved, improve algorithm to the robustness of initial value, second step directly solves Following iteration walks variable to reduce calculation amount.
Step 1:
1) it solves following formula and obtains vector λ:
(EET) λ=p-Eyk
2) it calculates
3)And calculate Jacobian matrix
Step 2:
1) it solves following formula and obtains xk+1:
2) y is updatedk+1=f-1(Cxk+1);
If 3) | | p-Eyk+1||< ε then restrains, otherwise k=k+1, and return step 1 continues iteration.
Beneficial achievement of the invention is: solving equation by the introducing of confactor and using two-step method, works as initial value x0When deviateing farther out with the true solution of equation, classical Newton method, which is easy to appear, does not restrain or restrains slower situation.The present invention By constructing a least square problem, finds and meet under constraint as close possible to the linearisation point really solved, can effectively mention Better utility and robustness.
Detailed description of the invention
Fig. 1 is a kind of specific implementation example diagram of the invention: containing 115 node isolated island micro-capacitance sensor example figures.
Fig. 2 is the flow chart of the isolated island micro-capacitance sensor tidal current computing method solved the present invention is based on two step of confactor.
Fig. 3 is the present invention and Newton method (N-R), the adaptive LM algorithm (A-LM) of single step, three step LM algorithms (MTLM) not With convergence number comparison diagram under precision.
Specific embodiment
In the following with reference to the drawings and specific embodiments, the present invention will be described in detail method and specific steps, it should be understood that these implementations Example is only illustrative of the invention and is not intended to limit the scope of the invention, after the present invention has been read, those skilled in the art The application range as defined in the appended claims is fallen within to the modification of various equivalent forms of the invention.
The present embodiment illustrates specific method of the present invention by taking large-scale 115 node example of micro-capacitance sensor as an example.The example has 115 Node, 118 branches contain 3 Wind turbines, 2 photovoltaic cells and 8 gas turbines, as shown in Fig. 1.
1, Load flow calculation equation is rewritten
Isolated island micro-capacitance sensor power flow equation can abstract representation are as follows:
The x ∈ of F (x)=0 Rn (12)
If n node is shared in system, wherein there is nPQA PQ node, nPVA PV node and nDA sagging node, then formula (1) Contain n+n for onePQ+nDThe n+n of a location variablePQ+nDTie up equation.
By moving on to all constant terms in formula (1) on the right side of equation, and retain PV node voltage magnitude, formula (1) can be changed It is written as following form:
F1(x)=p (13)
In formula:
P=[p0,V1,…,VnPV]T (14)
p0For in power flow equation constant term form vector, and be subsequent convenience of calculation, the electricity of PV node is remained in p Pressure amplitude value indicates by its quadratic form,
For vector x, will wherein voltage term and system frequency be replaced by logarithmic form:
X=[δ2,…,δn,a1,…,an,w]T (15)
In formula: δiFor node i phase angle, ai=ln Vi=2 ln Ui, w=ln ω.
2, confactor is introduced
Introduce confactor y:
In formula: item U second from the bottomD,iFor sagging node voltage amplitude, n is sharedD?;α, β are active and reactive load factor, For the node i at each branch both ends, j enables Kij=UiUj cosδij、Lij=UiUj sinδij, wherein δijij
Introduce confactor u:
In formula: item second from the bottom is identical as formula (2), represents sagging node;aij=ai+aj, δijij
Using confactor, former power flow equation can be become following form:
In formula: nonlinear transformation f () reversible between confactor, inverse transformation y=f-1(u) concrete form is as follows:
3, trend initial value design
Set flow state initial guess x0, the number of iterations k=0 is taken, convergence precision ε, maximum number of iterations k are setmax, System frequency ω0, voltage magnitude U0With phase angle initial value δ0, take y0=f-1(Cx0)。
4, two-step method solves equation
After completing to the transformation of power flow equation, two-step method can be used, transformed equation is solved, specific steps are such as Under:
Step 1:
1) it solves following formula and obtains vector λ:
(EET) λ=p-Eyk (20)
2) it calculates
3)And the Jacobian matrix of calculating formula (8)
Step 2:
1) it solves following formula and obtains xk+1:
2) y is updatedk+1=f-1(Cxk+1)。
If 3) | | p-Eyk+1||< ε then restrains, otherwise k=k+1, and return step 1 continues iteration.
It will be calculated in 115 node large size micro-capacitance sensors using after the method for the present invention with Newton method (N-R), the adaptive LM of single step Method (A-LM), three step LM methods (MTLM) compare, recorded under initial value variation algorithms of different not restrain number as shown in table 1 (number of run is 100 times):
1 initial value of table changes lower 4 kinds of algorithmic statement situations comparison
It is as shown in table 2 to expend the time in 1000 Monte Carlo simulations for 4 kinds of algorithms of comparison simultaneously:
2 Monte Carlo simulation of table expends the time

Claims (5)

1. a kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor, which is characterized in that including walking as follows It is rapid:
1) Load flow calculation equation is rewritten;
2) introduce confactor, by Load flow calculation it is equations turned be one group of underdetermined equation, one group of overdetermined equation and one group of auxiliary because Relation function between son;
3) trend initial value is set;
4) two-step method solves equation: solving to the equation after above-mentioned conversion, the first step first equations turned is by transformed Least square problem, second step directly solve following iteration step variable.
2. the isolated island micro-capacitance sensor tidal current computing method according to claim 1 solved based on two step of confactor, feature It is, rewriting Load flow calculation equation described in step 1), specifically: isolated island micro-capacitance sensor power flow equation can abstract representation are as follows:
The x ∈ of F (x)=0 Rn (1)
If n node is shared in system, wherein there is nPQA PQ node, nPVA PV node and nDA sagging node, then formula (1) is one It is a to contain n+nPQ+nDThe n+n of a location variablePQ+nDTie up equation;
By moving on to all constant terms in formula (1) on the right side of equation, and retain PV node voltage magnitude, formula (1) can be rewritten as Following form:
F1(x)=p (2)
In formula:
p0For in power flow equation constant term form vector, and be subsequent convenience of calculation, the voltage amplitude of PV node is remained in p Value, is indicated by its quadratic form,
For vector x, will wherein voltage term be replaced with system frequency ω by logarithmic form:
X=[δ2,…,δn,a1,…,an,w]T (4)
In formula: δiFor node i phase angle, ai=lnVi=2lnUi, w=ln ω.
3. the isolated island micro-capacitance sensor tidal current computing method according to claim 2 solved based on two step of confactor, feature It is, step 2) specifically:
Introduce confactor y:
In formula: item U second from the bottomD,iFor sagging node voltage amplitude, n is sharedD?;α, β are active and reactive load factor, for The node i at each branch both ends, j enable Kij=UiUjcosδij、Lij=UiUjsinδij, wherein δijij
Introduce confactor u:
In formula: aij=ai+aj, δijij
Using confactor, former power flow equation can be become following form:
In formula: nonlinear transformation f () reversible between confactor, inverse transformation y=f-1(u) concrete form is as follows:
4. the isolated island micro-capacitance sensor tidal current computing method according to claim 3 solved based on two step of confactor, feature It is, sets flow state initial guess x0, the number of iterations k=0 is taken, convergence precision ε, maximum number of iterations k are setmax, system Frequencies omega0, voltage magnitude U0With phase angle initial value δ0, take y0=f-1(Cx0)。
5. the isolated island micro-capacitance sensor tidal current computing method according to claim 4 solved based on two step of confactor, feature It is, after completing to the transformation of power flow equation, transformed equation is solved using two-step method, the specific steps are as follows:
Step 1:
1) it solves following formula and obtains vector λ:
(EET) λ=p-Eyk (9)
2) it calculates
3)And the Jacobian matrix of calculating formula (8)
Step 2:
1) it solves following formula and obtains xk+1:
2) y is updatedk+1=f-1(Cxk+1);
If 3) | | p-Eyk+1||< ε then restrains, otherwise k=k+1, and return step 1 continues iteration.
CN201910441049.3A 2019-05-24 2019-05-24 Island microgrid load flow calculation method based on two-step solution of auxiliary factors Active CN110148946B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910441049.3A CN110148946B (en) 2019-05-24 2019-05-24 Island microgrid load flow calculation method based on two-step solution of auxiliary factors

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910441049.3A CN110148946B (en) 2019-05-24 2019-05-24 Island microgrid load flow calculation method based on two-step solution of auxiliary factors

Publications (2)

Publication Number Publication Date
CN110148946A true CN110148946A (en) 2019-08-20
CN110148946B CN110148946B (en) 2022-03-25

Family

ID=67591989

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910441049.3A Active CN110148946B (en) 2019-05-24 2019-05-24 Island microgrid load flow calculation method based on two-step solution of auxiliary factors

Country Status (1)

Country Link
CN (1) CN110148946B (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110635485A (en) * 2019-10-09 2019-12-31 湖南大学 Load flow model calculation method for multi-terminal direct current system self-adaptive droop control
CN112398134A (en) * 2020-11-04 2021-02-23 浙江大学 Transmission and distribution cooperative rapid robust load flow calculation method

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003084041A1 (en) * 2002-03-28 2003-10-09 Curtin University Of Technology Power conversion system and method of converting power
CN102684201A (en) * 2012-05-30 2012-09-19 华南理工大学 Voltage threshold probability-based reactive power optimizing method for grid containing wind power plant
CN104578157A (en) * 2015-01-04 2015-04-29 云南电网有限责任公司电力科学研究院 Load flow calculation method of distributed power supply connection power grid
CN104767198A (en) * 2015-03-27 2015-07-08 国家电网公司 Computing method of power grid flow feasible solution
CN105490269A (en) * 2015-12-30 2016-04-13 中国南方电网有限责任公司 WAMS measurement-based multi-region power system state estimation method and system
CN106877336A (en) * 2017-03-14 2017-06-20 长沙理工大学 A kind of AC and DC power system Continuation power flow for considering wind-powered electricity generation randomness
CN109274116A (en) * 2018-09-12 2019-01-25 湘潭大学 A kind of alternating current-direct current mixing micro-capacitance sensor continuous tide calculation method
US10594136B1 (en) * 2015-07-02 2020-03-17 Antonio Trias System and method to simulate controls in the holomorphic embedding loadflow calculation

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2003084041A1 (en) * 2002-03-28 2003-10-09 Curtin University Of Technology Power conversion system and method of converting power
CN102684201A (en) * 2012-05-30 2012-09-19 华南理工大学 Voltage threshold probability-based reactive power optimizing method for grid containing wind power plant
CN104578157A (en) * 2015-01-04 2015-04-29 云南电网有限责任公司电力科学研究院 Load flow calculation method of distributed power supply connection power grid
CN104767198A (en) * 2015-03-27 2015-07-08 国家电网公司 Computing method of power grid flow feasible solution
US10594136B1 (en) * 2015-07-02 2020-03-17 Antonio Trias System and method to simulate controls in the holomorphic embedding loadflow calculation
CN105490269A (en) * 2015-12-30 2016-04-13 中国南方电网有限责任公司 WAMS measurement-based multi-region power system state estimation method and system
CN106877336A (en) * 2017-03-14 2017-06-20 长沙理工大学 A kind of AC and DC power system Continuation power flow for considering wind-powered electricity generation randomness
CN109274116A (en) * 2018-09-12 2019-01-25 湘潭大学 A kind of alternating current-direct current mixing micro-capacitance sensor continuous tide calculation method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
CATALINA GÓMEZ-QUILES 等: ""Computation of Maximum Loading Points"", 《IEEE TRANSACTIONS ON POWER SYSTEMS》 *
严正 等: ""自适应Levenberg-Marquardt 方法提高潮流计算收敛性"", 《中 国 电 机 工 程 学 报》 *
彭寒梅: ""微电网潮流计算及可靠性评估研究"", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110635485A (en) * 2019-10-09 2019-12-31 湖南大学 Load flow model calculation method for multi-terminal direct current system self-adaptive droop control
CN112398134A (en) * 2020-11-04 2021-02-23 浙江大学 Transmission and distribution cooperative rapid robust load flow calculation method

Also Published As

Publication number Publication date
CN110148946B (en) 2022-03-25

Similar Documents

Publication Publication Date Title
CN110110413B (en) Structural topology optimization method based on material field reduction progression expansion
Jordehi Enhanced leader particle swarm optimisation (ELPSO): An efficient algorithm for parameter estimation of photovoltaic (PV) cells and modules
Ghatak et al. A fast and efficient load flow technique for unbalanced distribution system
Gözel et al. An analytical method for the sizing and siting of distributed generators in radial systems
Lagacé et al. Improving power flow convergence by Newton Raphson with a Levenberg-Marquardt method
Kumar et al. Global maximum power point tracking using variable sampling time and pv curve region shifting technique along with incremental conductance for partially shaded photovoltaic systems
CN110148946A (en) A kind of isolated island micro-capacitance sensor tidal current computing method solved based on two step of confactor
CN106684857A (en) Linear optimal power flow model containing unified power flow controller
CN108683191A (en) A kind of Three-phase Power Flow analysis method of droop control type isolated island micro-capacitance sensor
CN105718634A (en) Airfoil robust optimization design method based on non-probability interval analysis model
de Araujo et al. Sensitivity analysis of convergence characteristics in power flow methods for distribution systems
CN106339562A (en) Energy interband non-local quantum tunneling simulation method with current conservation characteristic
CN103915840A (en) Method for estimating state of large power grid based on Givens orthogonal increment line transformation
CN104993525A (en) Active power distribution network coordination optimization control method considering ZIP loads
Zeng et al. Parameter identification of PV cell via adaptive compass search algorithm
Selvam Computation of pressures on Texas tech building
Yousri et al. An efficient photovoltaic modeling using an Adaptive Fractional-order Archimedes Optimization Algorithm: Validation with partial shading conditions
CN102780220A (en) Power flow calculation method for electric distribution network comprising PV constant distributed generation
Ali et al. Calculating multiple loadability points in the power flow solution space
Fahim et al. Overview of maximum power point tracking techniques for PV system
CN104113079B (en) MPPT control method and system
CN110417022A (en) The power distribution network Three Phase Power Flow of matrix multiplication operation extraction Jacobi&#39;s element
CN107958129B (en) Algorithm for simulating microcosmic current distribution of zinc oxide piezoresistor
Idema et al. Fast Newton load flow
CN107086603A (en) A kind of Random-fuzzy Continuation power flow of power system containing DFIG

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant