[go: up one dir, main page]
More Web Proxy on the site http://driver.im/

CN114251214A - Fractional order power system chaotic state judgment method and device - Google Patents

Fractional order power system chaotic state judgment method and device Download PDF

Info

Publication number
CN114251214A
CN114251214A CN202111500825.6A CN202111500825A CN114251214A CN 114251214 A CN114251214 A CN 114251214A CN 202111500825 A CN202111500825 A CN 202111500825A CN 114251214 A CN114251214 A CN 114251214A
Authority
CN
China
Prior art keywords
power system
differential equation
fractional order
water turbine
order differential
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
CN202111500825.6A
Other languages
Chinese (zh)
Other versions
CN114251214B (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.)
719th Research Institute of CSIC
Original Assignee
719th Research Institute of CSIC
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 719th Research Institute of CSIC filed Critical 719th Research Institute of CSIC
Priority to CN202111500825.6A priority Critical patent/CN114251214B/en
Publication of CN114251214A publication Critical patent/CN114251214A/en
Application granted granted Critical
Publication of CN114251214B publication Critical patent/CN114251214B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03B—MACHINES OR ENGINES FOR LIQUIDS
    • F03B11/00—Parts or details not provided for in, or of interest apart from, the preceding groups, e.g. wear-protection couplings, between turbine and generator
    • F03B11/008—Measuring or testing arrangements
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03B—MACHINES OR ENGINES FOR LIQUIDS
    • F03B11/00—Parts or details not provided for in, or of interest apart from, the preceding groups, e.g. wear-protection couplings, between turbine and generator
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03B—MACHINES OR ENGINES FOR LIQUIDS
    • F03B13/00—Adaptations of machines or engines for special use; Combinations of machines or engines with driving or driven apparatus; Power stations or aggregates
    • F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F03—MACHINES OR ENGINES FOR LIQUIDS; WIND, SPRING, OR WEIGHT MOTORS; PRODUCING MECHANICAL POWER OR A REACTIVE PROPULSIVE THRUST, NOT OTHERWISE PROVIDED FOR
    • F03B—MACHINES OR ENGINES FOR LIQUIDS
    • F03B3/00—Machines or engines of reaction type; Parts or details peculiar thereto
    • F03B3/12—Blades; Blade-carrying rotors
    • F03B3/121—Blades, their form or construction
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00—Computing arrangements based on specific mathematical models
    • G06N7/08—Computing arrangements based on specific mathematical models using chaos models or non-linear system models
    • Y—GENERAL 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/20—Hydro energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Data Mining & Analysis (AREA)
  • Chemical & Material Sciences (AREA)
  • Mechanical Engineering (AREA)
  • Combustion & Propulsion (AREA)
  • Algebra (AREA)
  • Software Systems (AREA)
  • Computing Systems (AREA)
  • Databases & Information Systems (AREA)
  • Operations Research (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Nonlinear Science (AREA)
  • Other Liquid Machine Or Engine Such As Wave Power Use (AREA)
  • Control Of Water Turbines (AREA)

Abstract

The invention discloses a fractional order power system chaotic state judgment method and a device, which relate to the technical field of water turbine control and comprise the following steps: step S10, expanding the fractional order differential equation of the hydraulic turbine power system to an integer order differential equation according to the integer order Jacobian matrix of the fractional order differential equation of the hydraulic turbine power system; step S20, solving an expanded integer order differential equation according to the initial rotor angle, the rotation speed, the driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain a Lyapunov exponent; and step S30, judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent. According to the invention, after the fractional order differential equation of the water turbine power system is expanded into the executable integer order differential equation, the calculation feasibility is good, the calculated amount is small, the oscillation is stable, the convergence speed is high, the computed Lyapunov exponent can directly judge the chaotic state of the fractional order water turbine power system, and the judgment result is accurate and efficient.

Description

Fractional order power system chaotic state judgment method and device
Technical Field
The invention relates to the technical field of water turbine control, in particular to a method and a device for judging a chaotic state of a fractional order power system.
Background
A non-linear system refers to a system in which the state and output variables of the system cannot be described in a linear relationship under the influence of external conditions. In a deterministic system, there is seemingly random irregular motion whose behavior appears as uncertainty, unrepeatable, unpredictable, a chaotic phenomenon. The chaos is the inherent characteristic of the nonlinear power system and is a ubiquitous phenomenon of the nonlinear system.
The Lyapunov exponent is the most direct and reliable quantitative index for judging whether a nonlinear power system is in a chaotic state or not, but the current Lyapunov exponents are all obtained through an integer-order power system, and a related acquisition method is lacked for the Lyapunov exponent of a fractional-order power system, so that the chaotic state of the system cannot be directly judged from the perspective of the Lyapunov exponent, and only the indexes such as initial condition sensitivity, chaotic attractor and the like are indirectly proved.
Disclosure of Invention
The embodiment of the invention provides a method and a device for judging a chaotic state of a fractional order power system, which are used for solving the technical problem that the chaotic state of the fractional order power system is difficult to judge in the related technology.
In a first aspect, a method for judging a chaotic state of a fractional order power system is provided, where the method includes:
expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system;
solving an expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain a Lyapunov exponent;
and judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
In some embodiments, the step of expanding the fractional order differential equation of the turbine power system to an integer order differential equation according to an integer order jacobian matrix of the fractional order differential equation of the turbine power system comprises:
establishing a fractional order differential equation of a water turbine power system:
Figure BDA0003402517490000021
wherein mu is (0-1)]Of any order in between, DμAs a differential operator, δ is the rotor angle of the turbine, ω is the rotational speed of the turbine, mtIs the driving moment of the water turbine, y is the servomotor stroke of the water turbine, KdIs a differential control coefficient.
In some embodiments, the step of expanding the fractional order differential equation of the water turbine power system to an integer order differential equation based on an integer order jacobian of the fractional order differential equation of the water turbine power system comprises:
and expanding the fractional order differential equation of the water turbine power system to an integer order differential equation by using a prediction-correction algorithm.
In some embodiments, the step of expanding the fractional order differential equation to the integer order differential equation of the turbine power system using a pre-estimation-correction algorithm comprises:
the fractional order differential equation to the integer order differential equation of the expansion turbine power system is as follows:
Figure BDA0003402517490000022
Figure BDA0003402517490000031
Figure BDA0003402517490000032
wherein Γ is a gamma function, and h is KdAlpha and beta are intermediate process parameters.
In some embodiments, the step of solving an expanded integer order differential equation according to the rotor angle, the rotational speed, the driving torque, the servomotor stroke, and the differential control coefficient of the initial water turbine to obtain the lyapunov exponent includes:
substituting the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine into an expanded integer order differential equation;
with h as KdThe step length of calculation is subjected to iterative calculation to obtain the following KdA varied lyapunov index.
In some embodiments, the term h is KdThe step length of calculation is subjected to iterative calculation to obtain the following KdA step of varying the Lyapunov exponent comprising:
k is 0.01 ═ hdThe step length of calculation is subjected to iterative calculation to obtain the following KdA varied lyapunov index.
In a second aspect, a device for determining a chaotic state of a fractional order power system is provided, where the device includes:
the expansion unit is used for expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system;
the solving unit is used for solving the expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain the Lyapunov exponent;
and the judging unit is used for judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
In some embodiments, the determining means further comprises:
the establishing unit is used for establishing a fractional order differential equation of the water turbine power system:
Figure BDA0003402517490000041
wherein mu is (0-1)]Of any order in between, DμAs a differential operator, δ is the rotor angle of the turbine, ω is the rotational speed of the turbine, mtIs the driving moment of the water turbine, y is the servomotor stroke of the water turbine, KdIs a differential control coefficient.
In a third aspect, a computer device is provided, comprising: the device comprises a memory and a processor, wherein at least one instruction is stored in the memory, and the at least one instruction is loaded and executed by the processor so as to realize the judgment method of the fractional order power system chaotic state.
In a fourth aspect, a computer-readable storage medium is provided, which stores computer instructions, and when the computer instructions are executed by a computer, the computer is enabled to execute the aforementioned method for determining the chaotic state of a fractional order power system.
The technical scheme provided by the invention has the beneficial effects that:
the method and the device for judging the chaotic state of the fractional order power system in the embodiment of the invention comprise the steps of firstly expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system, then solving the expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain a Lyapunov exponent, and finally judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent. According to the invention, after the fractional order differential equation of the water turbine power system is expanded into the executable integer order differential equation, the calculation feasibility is good, the calculated amount is small, the oscillation is stable, the convergence speed is high, the computed Lyapunov exponent can directly judge the chaotic state of the fractional order water turbine power system, and the judgment result is accurate and efficient.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for determining a chaotic state of a fractional order power system according to an embodiment of the present invention;
fig. 2 is a schematic flowchart of the step S20 implemented in fig. 1 according to an embodiment of the present invention;
FIG. 3 is a Lyapunov exponent diagram with order μ taken to be 1, according to an embodiment of the present invention;
FIG. 4 is a Lyapunov exponent diagram of a conventional integer order turbine power system;
FIG. 5 is a Lyapunov exponent plot with order μ taken to be 0.9, provided by an embodiment of the present invention;
FIG. 6 is a Lyapunov exponent plot with order μ taken to be 0.8, provided by an embodiment of the present invention;
fig. 7 is a schematic structural diagram of a fractional order power system chaotic state judgment device according to an embodiment of the present invention;
fig. 8 is a schematic structural diagram of a computer device according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be obtained by a person skilled in the art without any inventive step based on the embodiments of the present invention, are within the scope of the present invention.
The embodiment of the invention provides a method for judging a chaotic state of a fractional order power system, which can solve the technical problem that the chaotic state of the existing fractional order power system is difficult to judge.
Referring to fig. 1, an embodiment of the present invention provides a method for determining a chaotic state of a fractional order power system, where the method includes:
and step S10, expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to the integer order Jacobian matrix of the fractional order differential equation of the water turbine power system.
Specifically, before the step of expanding the fractional order differential equation of the hydraulic turbine power system to the integer order differential equation according to the integer order jacobian matrix of the fractional order differential equation of the hydraulic turbine power system, the method includes:
establishing a fractional order differential equation of a water turbine power system:
Figure BDA0003402517490000061
wherein mu is (0-1)]Of any order in between, DμAs a differential operator, δ is the rotor angle of the turbine, ω is the rotational speed of the turbine, mtIs the driving moment of the water turbine, y is the servomotor stroke of the water turbine, KdIs a differential control coefficient.
Further, the step of expanding the fractional order differential equation of the hydraulic turbine power system to the integer order differential equation according to the integer order jacobian matrix of the fractional order differential equation of the hydraulic turbine power system includes:
and expanding the fractional order differential equation of the water turbine power system to an integer order differential equation by using a prediction-correction algorithm.
Further, the step of expanding the fractional order differential equation of the water turbine power system to an integer order differential equation using a pre-estimation-correction algorithm includes:
the fractional order differential equation to the integer order differential equation of the expansion turbine power system is as follows:
Figure BDA0003402517490000071
Figure BDA0003402517490000072
Figure BDA0003402517490000073
wherein Γ is a gamma function, and h is KdAlpha and beta are intermediate process parameters.
And expanding the fractional order differential equation to the required integer order differential equation by utilizing a pre-estimation-correction algorithm, ensuring that the structure of the integer order differential equation of the Jacobian matrix part is not changed, and facilitating subsequent solution and calculation to obtain the Lyapunov exponent.
And step S20, solving the expanded integer order differential equation according to the initial rotor angle, the rotation speed, the driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain the Lyapunov exponent.
The step of solving an expanded integer order differential equation according to the rotor angle, the rotating speed, the driving moment, the servomotor stroke and the differential control coefficient of the initial water turbine to obtain the Lyapunov exponent comprises the following steps:
and substituting the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine into the expanded integer order differential equation.
With h as KdThe step length of calculation is subjected to iterative calculation to obtain the following KdA varied lyapunov index. Preferably, the calculation step h takes 0.01.
Specifically, equations (3) and (4) are expanded and substituted into equation (2), and the reconstructed integer order differential equation of the water turbine power system is expressed as follows:
Figure BDA0003402517490000081
in equation (5):
Figure BDA0003402517490000082
Figure BDA0003402517490000091
Figure BDA0003402517490000101
Figure BDA0003402517490000102
wherein i1Taking 0 and 1; i.e. i2Taking 0, 1 and 2; i.e. i3Taking 0, 1, 2 and 3; i.e. i4Taking 0, 1, 2, 3 and 4; i.e. i5Take 0, 1, 2, 3, 4, 5.
Referring to fig. 2, the rotor angle, the rotational speed, the driving moment, the servomotor stroke and the differential control coefficient of the initial water turbine are substituted into the expanded integer order differential equation, numerical iteration solution is carried out on the reconstructed equation (5), and the Lyapunov exponent is calculated by referring to the Lyapunov calculation method (wolf calculation method) of the existing integer order differential equation.
And step S30, judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
When mu is 1.0, the fractional order water turbine power system is equivalent to an integer order water turbine power system in the embodiment of the invention, and K isdReferring to fig. 3 and 4, the lyapunov exponent calculated by the fractional order power system chaotic state judgment method is basically consistent with the lyapunov exponent calculated by the existing integer order water turbine power system chaotic state judgment method in the range of 7-10, which shows that the computed lyapunov exponent of the fractional order power system chaotic state judgment method is high in accuracy and can be used for directly judging the fractional order water turbine power system chaotic state. Among them, as can be seen from FIGS. 3 and 4, at KdNear 7.8, the turbine power system will alternate continuously between chaotic motion and stable periodic motion, and when K isd>8.0 later, LE1 is more than 0, the water turbine power system will have complex chaos behavior until KdWhen the LE1 is larger than zero, the water turbine power system has chaotic motion, and is in a destabilization state, namely 10.0.
When mu is 0.9 and 0.8, the water turbine power system is in fractional order, and the Lyapunov exponent obtained by adopting the fractional order power system chaotic state judgment method is shown in figures 5 and 6. When mu is 0.9, when K isd>After 9.5, LE1 is larger than 0, and the water turbine power system has complex chaotic behavior. When mu is 0.9, when K isd>After 8.8, LE1 is larger than 0, and the water turbine power system has complex chaotic behavior.
The method for judging the chaotic state of the fractional order power system in the embodiment of the invention comprises the steps of firstly expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system, then solving the expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the main moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain a Lyapunov exponent, and finally judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent. According to the invention, after the fractional order differential equation of the water turbine power system is expanded into the executable integer order differential equation, the calculation feasibility is good, the calculated amount is small, the oscillation is stable, the convergence speed is high, the computed Lyapunov exponent can directly judge the chaotic state of the fractional order water turbine power system, and the judgment result is accurate and efficient.
Referring to fig. 7, an embodiment of the present invention further provides a device for determining a chaotic state of a fractional order power system, where the device includes:
the expansion unit is used for expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system;
the solving unit is used for solving the expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain the Lyapunov exponent;
and the judging unit is used for judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
The device for determining the chaotic state of the fractional order power system provided by the above embodiment may be implemented in the form of a computer program, and the computer program may be run on a computer device as shown in fig. 8.
An embodiment of the present invention further provides a computer device, including: the method comprises the steps of storing at least one instruction in a memory, a processor and a network interface which are connected through a system bus, wherein the at least one instruction is loaded and executed by the processor so as to realize all or part of the steps of the fractional order power system chaotic state judgment method.
The network interface is used for performing network communication, such as sending distributed tasks. Those skilled in the art will appreciate that the architecture shown in fig. 8 is merely a block diagram of some of the structures associated with the inventive arrangements and is not intended to limit the computing devices to which the inventive arrangements may be applied, as a particular computing device may include more or less components than those shown, or may combine certain components, or have a different arrangement of components.
The Processor may be a CPU, other general purpose Processor, a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf Programmable Gate Array (FPGA) or other Programmable logic device, discrete Gate or transistor logic device, discrete hardware component, or the like. The general purpose processor may be a microprocessor or the processor may be any conventional processor or the like, the processor being the control center of the computer device and the various interfaces and lines connecting the various parts of the overall computer device.
The memory may be used to store computer programs and/or modules, and the processor may implement various functions of the computer device by running or executing the computer programs and/or modules stored in the memory, as well as by invoking data stored in the memory. The memory may mainly include a storage program area and a storage data area, wherein the storage program area may store an operating system, an application program required by at least one function (such as a video playing function, an image playing function, etc.), and the like; the storage data area may store data (such as video data, image data, etc.) created according to the use of the cellular phone, etc. In addition, the memory may include high speed random access memory, and may also include non-volatile memory, such as a hard disk, a memory, a plug-in hard disk, a Smart Media Card (SMC), a Secure Digital (SD) Card, a Flash memory Card (Flash Card), at least one magnetic disk storage device, a Flash memory device, or other volatile solid state storage device.
Wherein, in one embodiment, the processor is configured to execute a computer program stored in the memory to implement the steps of:
and step S10, expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to the integer order Jacobian matrix of the fractional order differential equation of the water turbine power system.
And step S20, solving the expanded integer order differential equation according to the initial rotor angle, the rotation speed, the driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain the Lyapunov exponent.
And step S30, judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
The embodiment of the invention also provides a computer readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, all steps or part of steps of the method for judging the chaotic state of the fractional order power system are realized.
The embodiment of the present invention may implement all or part of the foregoing processes, and may also be implemented by instructing related hardware by a computer program, where the computer program may be stored in a computer-readable storage medium, and when the computer program is executed by a processor, the computer program may implement the steps of the foregoing methods. Wherein the computer program comprises computer program code, which may be in the form of source code, object code, an executable file or some intermediate form, etc. The computer readable medium may include: any entity or device capable of carrying computer program code, recording medium, U-disk, removable hard disk, magnetic disk, optical disk, computer memory, Read-Only memory (ROM), Random Access Memory (RAM), electrical carrier wave signals, telecommunications signals, software distribution media, and the like. It should be noted that the computer readable medium may contain other components which may be suitably increased or decreased as required by legislation and patent practice in jurisdictions, for example, in some jurisdictions, in accordance with legislation and patent practice, the computer readable medium does not include electrical carrier signals and telecommunications signals.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, server, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, optical storage, and the like) having computer-usable program code embodied therein.
It should be noted that, in this document, the terms "comprises," "comprising," or any other variation thereof, are intended to cover a non-exclusive inclusion, such that a process, method, article, or system that comprises a list of elements does not include only those elements but may include other elements not expressly listed or inherent to such process, method, article, or system. Without further limitation, an element defined by the phrase "comprising an … …" does not exclude the presence of other like elements in a process, method, article, or system that comprises the element.
The above-mentioned serial numbers in the embodiments of the present invention are merely for description and do not represent the merits of the embodiments.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing are merely exemplary embodiments of the present invention, which enable those skilled in the art to understand or practice the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.

Claims (10)

1. A fractional order power system chaotic state judgment method is characterized by comprising the following steps:
expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system;
solving an expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain a Lyapunov exponent;
and judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
2. The method for determining the chaotic state of the fractional order power system according to claim 1, wherein the step of expanding the fractional order differential equation of the hydraulic turbine power system to an integer order differential equation according to the integer order jacobian matrix of the fractional order differential equation of the hydraulic turbine power system comprises:
establishing a fractional order differential equation of a water turbine power system:
Figure FDA0003402517480000011
wherein mu is (0-1)]Of any order in between, DμAs a differential operator, δ is the rotor angle of the turbine, ω is the rotational speed of the turbine, mtIs the driving moment of the water turbine, y is the servomotor stroke of the water turbine, KdIs a differential control coefficient.
3. The method for determining the chaotic state of the fractional order power system according to claim 2, wherein the step of expanding the fractional order differential equation of the hydraulic turbine power system to an integer order differential equation according to an integer order jacobian matrix of the fractional order differential equation of the hydraulic turbine power system comprises:
and expanding the fractional order differential equation of the water turbine power system to an integer order differential equation by using a prediction-correction algorithm.
4. The method for determining the chaotic state of the fractional order power system according to claim 3, wherein the step of expanding the fractional order differential equation of the water turbine power system to an integer order differential equation using a pre-estimation-correction algorithm comprises:
the fractional order differential equation to the integer order differential equation of the expansion turbine power system is as follows:
Figure FDA0003402517480000021
Figure FDA0003402517480000022
Figure FDA0003402517480000023
wherein Γ is a gamma function, and h is KdAlpha and beta are intermediate process parameters.
5. The method for judging the chaotic state of the fractional order power system according to claim 4, wherein the step of solving an expanded integer order differential equation according to the rotor angle, the rotating speed, the driving moment, the servomotor stroke and the differential control coefficient of the initial water turbine to obtain the Lyapunov exponent comprises:
substituting the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine into an expanded integer order differential equation;
with h as KdThe step length of calculation is subjected to iterative calculation to obtain the following KdA varied lyapunov index.
6. The method for determining the chaotic state of a fractional order power system as claimed in claim 5, wherein h is KdThe step size of the calculation of (a) is iteratively calculated,to obtain as KdA step of varying the Lyapunov exponent comprising:
k is 0.01 ═ hdThe step length of calculation is subjected to iterative calculation to obtain the following KdA varied lyapunov index.
7. A device for judging the chaotic state of a fractional order power system is characterized by comprising:
the expansion unit is used for expanding the fractional order differential equation of the water turbine power system to an integer order differential equation according to an integer order Jacobian matrix of the fractional order differential equation of the water turbine power system;
the solving unit is used for solving the expanded integer order differential equation according to the initial rotor angle, the initial rotating speed, the initial driving moment, the servomotor stroke and the differential control coefficient of the water turbine to obtain the Lyapunov exponent;
and the judging unit is used for judging the chaotic state of the fractional order water turbine power system according to the Lyapunov exponent.
8. The apparatus for determining a chaotic status of a fractional order power system as claimed in claim 7, further comprising:
the establishing unit is used for establishing a fractional order differential equation of the water turbine power system:
Figure FDA0003402517480000031
wherein mu is (0-1)]Of any order in between, DμAs a differential operator, δ is the rotor angle of the turbine, ω is the rotational speed of the turbine, mtIs the driving moment of the water turbine, y is the servomotor stroke of the water turbine, KdIs a differential control coefficient.
9. A computer device, comprising: the device comprises a memory and a processor, wherein at least one instruction is stored in the memory, and is loaded and executed by the processor to realize the fractional order power system chaotic state judgment method in any one of claims 1 to 6.
10. A computer-readable storage medium characterized by: the computer storage medium stores computer instructions that, when executed by a computer, cause the computer to perform the fractional order power system chaotic state determination method of any one of claims 1 to 6.
CN202111500825.6A 2021-12-09 2021-12-09 Fractional order power system chaotic state judgment method and device Active CN114251214B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111500825.6A CN114251214B (en) 2021-12-09 2021-12-09 Fractional order power system chaotic state judgment method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111500825.6A CN114251214B (en) 2021-12-09 2021-12-09 Fractional order power system chaotic state judgment method and device

Publications (2)

Publication Number Publication Date
CN114251214A true CN114251214A (en) 2022-03-29
CN114251214B CN114251214B (en) 2023-01-24

Family

ID=80794434

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111500825.6A Active CN114251214B (en) 2021-12-09 2021-12-09 Fractional order power system chaotic state judgment method and device

Country Status (1)

Country Link
CN (1) CN114251214B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20020015547A (en) * 2000-08-22 2002-02-28 심규석 recording medium for sensitivity treatment
US20070166168A1 (en) * 2006-01-16 2007-07-19 Volvo Construction Equipment Ab Control system for a work machine and method for controlling a hydraulic cylinder in a work machine
CN110361974A (en) * 2019-07-16 2019-10-22 华北水利水电大学 Turbine Governor System optimization method based on BP-FOA hybrid algorithm
CN111666628A (en) * 2020-05-13 2020-09-15 哈尔滨工业大学 Adaptive fault-tolerant control method for supercavitation navigation body
CN112989556A (en) * 2021-01-13 2021-06-18 华北电力大学 Transient stability discrimination method based on maximum Lyapunov index of switching system

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20020015547A (en) * 2000-08-22 2002-02-28 심규석 recording medium for sensitivity treatment
US20070166168A1 (en) * 2006-01-16 2007-07-19 Volvo Construction Equipment Ab Control system for a work machine and method for controlling a hydraulic cylinder in a work machine
CN110361974A (en) * 2019-07-16 2019-10-22 华北水利水电大学 Turbine Governor System optimization method based on BP-FOA hybrid algorithm
CN111666628A (en) * 2020-05-13 2020-09-15 哈尔滨工业大学 Adaptive fault-tolerant control method for supercavitation navigation body
CN112989556A (en) * 2021-01-13 2021-06-18 华北电力大学 Transient stability discrimination method based on maximum Lyapunov index of switching system

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李俊益: "分数阶水轮机调节系统分岔与混沌特性研究", 《中国博士学位论文全文数据库 基础科学辑》 *

Also Published As

Publication number Publication date
CN114251214B (en) 2023-01-24

Similar Documents

Publication Publication Date Title
US10929746B2 (en) Low-power hardware acceleration method and system for convolution neural network computation
US20190079975A1 (en) Scheduling method and system based on hybrid variable neighborhood search and gravitational search algorithm
CN110390075B (en) Matrix preprocessing method, device, terminal and readable storage medium
CN111882038A (en) Model conversion method and device
CN113867685A (en) Multiplier conversion method, device and equipment and readable storage medium
CN118520975A (en) Large language model training method and device, electronic equipment and storage medium
CN114251214B (en) Fractional order power system chaotic state judgment method and device
CN114819159A (en) Inference method, device, equipment and storage medium of deep learning model
US11886832B2 (en) Operation device and operation method
CN112861239A (en) Method, system, device and storage medium for initially balancing ground stress of numerical model
CN112507323A (en) Model training method and device based on unidirectional network and computing equipment
CN116912653A (en) Model training method and device and electronic equipment
CN115170917B (en) Image processing method, electronic device and storage medium
CN116842532A (en) Data processing method, device, computer equipment and computer readable storage medium
CN114254453B (en) Chaos control method and device for six-dimensional fractional order power system
CN111383156B (en) Image processing method and device, intelligent driving system and vehicle-mounted operation platform
CN113971453A (en) Data normalization processing method, storage medium and computer equipment
CN117792800B (en) Information verification method and system based on Internet of things security evaluation system
CN113434571B (en) Knowledge mining method, platform, system and storage medium
CN112445525B (en) Data processing method, related equipment and computer readable medium
CN118589933A (en) Method and device for calibrating MTPA current instruction of asynchronous motor
Mohamed A novel fast hybrid GCD computation algorithm
CN109284285B (en) Data processing method, device, computer equipment and computer readable storage medium
CN114676827A (en) Model acceleration method, device, integrated circuit and storage medium
CN114301088A (en) Inertia control method, device, equipment and medium of wind turbine generator and wind turbine generator

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