Harmonic interaction characteristic analysis method, device and equipment based on chaotic characteristic, storage medium and program product
By constructing an equivalent circuit model of the harmonic source and calculating the chaotic interaction factor, the nonlinear interaction problem between harmonic sources in the multi-harmonic source system is solved, and the system stability and power quality are improved.
Patent Information
- Application Number
- CN202510853647.7
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2025-06-24
- Publication Date
- 2025-10-03
AI Technical Summary
Traditional harmonic modeling and analysis methods cannot effectively consider the nonlinear interactions between harmonic sources in multi-harmonic source systems, resulting in low system stability. In particular, when the load changes and the control strategy is frequently adjusted, the harmonic emission level is difficult to characterize.
By obtaining the harmonic impedances of multiple harmonic sources, an equivalent circuit model of the transmission line is constructed, the series impedance and parallel admittance are determined, and the chaotic interaction factor is calculated in combination with the Lyapunov exponent to analyze the interaction characteristics between the harmonic sources.
It realizes the quantitative evaluation of harmonic interaction between harmonic sources, improves the stability and power quality of multi-harmonic source systems, and optimizes the design of power systems.
Smart Images

Figure CN120749698A_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the field of new energy grid-connected technology, and in particular to a method, device, equipment, storage medium and program product for analyzing harmonic interaction characteristics based on chaos characteristics. Background Art
[0002] Driven by the dual goals of "dual carbon" and the strategy of "building a new power system," the dual-high characteristics of the power system are becoming increasingly prominent. Distributed renewable energy, electric vehicles, flexible loads, and other types of sources connected to the power system through power electronic conversion devices such as inverters and converters inject significant amounts of harmonics into the grid during operation due to their nonlinear characteristics and pulse width modulation (PWM) strategies. The harmonics emitted by these multiple harmonic sources interact and couple with each other, exacerbating uncertainty in power system analysis and making it difficult to characterize harmonic emission levels.
[0003] Traditional harmonic modeling and analysis methods are often based on the assumption of a single harmonic source, ignoring the impact of multi-source coupling and failing to fully consider the nonlinear interactions between harmonic sources in multi-harmonic source systems (such as inverters). This is especially true when the multi-harmonic source system is subject to load changes, grid-connected operations, and frequent adjustments to control strategies. The harmonic interactions between harmonic sources often exhibit complex dynamic characteristics and nonlinear behaviors, resulting in low stability in the multi-harmonic source system.
[0004] Therefore, how to improve the stability of the multi-harmonic source system has become an urgent problem to be solved. Summary of the Invention
[0005] The embodiments of the present application provide a method, apparatus, device, storage medium, and program product for analyzing harmonic interaction characteristics based on chaotic characteristics, which can be beneficial for improving the stability of a multi-harmonic source system.
[0006] In a first aspect, an embodiment of the present application provides a method for analyzing harmonic interaction characteristics based on chaotic characteristics, the method comprising:
[0007] Obtaining harmonic impedances corresponding to a plurality of harmonic sources respectively; the plurality of harmonic sources include a first harmonic source and a second harmonic source;
[0008] Based on a pre-constructed equivalent circuit model of a transmission line between a plurality of harmonic sources, the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source are determined, and based on the series impedance and the parallel admittance, the equivalent connection impedance between the first harmonic source and the second harmonic source is determined;
[0009] The first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov index are input into the preset chaos interaction factor calculation model to obtain the chaos interaction factor between the first harmonic source and the second harmonic source. The chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect;
[0010] Based on the chaotic interaction factor, the interaction characteristic analysis results between the first harmonic source and the second harmonic source are determined.
[0011] In one embodiment, the harmonic impedances corresponding to the multiple harmonic sources are determined in the following manner: for each harmonic source among the multiple harmonic sources, a target harmonic signal is injected into the input end of the harmonic source through a signal generator; the target harmonic signal is determined based on the operating frequency of the harmonic source; voltage response data and current response data corresponding to the output end of the harmonic source under the target harmonic signal are obtained; the amplitude of the harmonic voltage is determined based on the voltage response data, and the amplitude of the harmonic current is determined based on the current response data; and the harmonic impedance corresponding to the harmonic source is determined based on the amplitude of the harmonic voltage and the amplitude of the harmonic current.
[0012] In one embodiment, determining the harmonic impedance corresponding to the harmonic source based on the amplitude of the harmonic voltage and the amplitude of the harmonic current includes: determining a quotient of the amplitude of the harmonic voltage and the amplitude of the harmonic current; determining a difference between a first phase angle of the harmonic voltage and a second phase angle of the harmonic current; and determining the harmonic impedance corresponding to the harmonic source based on the quotient and the difference.
[0013] In one embodiment, determining an equivalent connection impedance between a first harmonic source and a second harmonic source based on the series impedance and the parallel admittance includes: determining an input impedance of a transmission line based on the series impedance and the parallel admittance; obtaining capacitive reactance and inductive reactance of a passive filter between the first harmonic source and the second harmonic source, and determining an equivalent impedance of the passive filter based on the capacitive reactance and the inductive reactance; and determining an equivalent connection impedance between the first harmonic source and the second harmonic source based on the input impedance and the equivalent impedance.
[0014] In one embodiment, the input impedance of the transmission line is determined based on the series impedance and the shunt admittance, including: determining the line characteristic impedance of the transmission line based on the series impedance and the shunt admittance; and determining the input impedance of the transmission line based on the line characteristic impedance and the line propagation coefficient.
[0015] In one embodiment, the preset chaotic interaction factor calculation model is as follows:
[0016]
[0017] Among them, H ijIt represents the chaotic interaction factor between the i-th harmonic source and the j-th harmonic source; Z ij (h) represents the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source under the h-th harmonic; Z i (h) represents the harmonic impedance of the i-th harmonic source under the h-th harmonic; Z j (h) represents the harmonic impedance of the jth harmonic source under the hth harmonic; It represents the Lyapunov exponent between the i-th harmonic source and the j-th harmonic source under the h-th harmonic.
[0018] In a second aspect, the present application provides a device for analyzing harmonic interaction characteristics based on chaotic characteristics, the device comprising:
[0019] An acquisition module, configured to acquire harmonic impedances corresponding to a plurality of harmonic sources, including a first harmonic source and a second harmonic source;
[0020] a determination module, configured to determine the series impedance and parallel admittance of the transmission line between the first harmonic source and the second harmonic source based on a pre-constructed equivalent circuit model of the transmission line between the plurality of harmonic sources, and to determine the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance and the parallel admittance;
[0021] a processing module, configured to input the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov exponent into a preset chaos interaction factor calculation model to obtain a chaos interaction factor between the first harmonic source and the second harmonic source. The chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect;
[0022] The determination module is further used to determine the analysis result of the interaction characteristics between the first harmonic source and the second harmonic source based on the chaotic interaction factor.
[0023] In a third aspect, the present application provides a computer device comprising a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are performed:
[0024] Obtaining harmonic impedances corresponding to a plurality of harmonic sources respectively; the plurality of harmonic sources include a first harmonic source and a second harmonic source;
[0025] Based on a pre-constructed equivalent circuit model of a transmission line between a plurality of harmonic sources, the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source are determined, and based on the series impedance and the parallel admittance, the equivalent connection impedance between the first harmonic source and the second harmonic source is determined;
[0026] The first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov index are input into the preset chaos interaction factor calculation model to obtain the chaos interaction factor between the first harmonic source and the second harmonic source. The chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect;
[0027] Based on the chaotic interaction factor, the interaction characteristic analysis results between the first harmonic source and the second harmonic source are determined.
[0028] In a fourth aspect, the present application further provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements the following steps:
[0029] Obtaining harmonic impedances corresponding to a plurality of harmonic sources respectively; the plurality of harmonic sources include a first harmonic source and a second harmonic source;
[0030] Based on a pre-constructed equivalent circuit model of a transmission line between a plurality of harmonic sources, the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source are determined, and based on the series impedance and the parallel admittance, the equivalent connection impedance between the first harmonic source and the second harmonic source is determined;
[0031] The first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov index are input into the preset chaos interaction factor calculation model to obtain the chaos interaction factor between the first harmonic source and the second harmonic source. The chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect;
[0032] Based on the chaotic interaction factor, the interaction characteristic analysis results between the first harmonic source and the second harmonic source are determined.
[0033] In a fifth aspect, the present application further provides a computer program product, comprising a computer program, which, when executed by a processor, implements the following steps:
[0034] Obtaining harmonic impedances corresponding to a plurality of harmonic sources respectively; the plurality of harmonic sources include a first harmonic source and a second harmonic source;
[0035] Based on a pre-constructed equivalent circuit model of a transmission line between a plurality of harmonic sources, the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source are determined, and based on the series impedance and the parallel admittance, the equivalent connection impedance between the first harmonic source and the second harmonic source is determined;
[0036] The first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov index are input into the preset chaos interaction factor calculation model to obtain the chaos interaction factor between the first harmonic source and the second harmonic source. The chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect;
[0037] Based on the chaotic interaction factor, the interaction characteristic analysis results between the first harmonic source and the second harmonic source are determined.
[0038] The above-mentioned harmonic interaction characteristic analysis method, device, equipment and storage medium based on chaotic characteristics, the computer equipment can obtain the harmonic impedances corresponding to multiple harmonic sources respectively; the first harmonic source and the second harmonic source included in the multiple harmonic sources; based on the pre-constructed equivalent circuit model of the transmission line between the multiple harmonic sources, the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source are determined, and based on the series impedance and the parallel admittance, the equivalent connection impedance between the first harmonic source and the second harmonic source is determined; the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance and the Lyapunov index are input into the preset chaos interaction factor calculation model to obtain the chaos interaction factor between the first harmonic source and the second harmonic source, and the chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect; based on the chaos interaction factor, the interaction characteristic analysis result between the first harmonic source and the second harmonic source is determined. By using this method, by introducing the Lyapunov index, a chaos indicator, the complexity of harmonic interaction between harmonic sources (i.e., the chaotic interaction factor) can be quantitatively evaluated, and the analysis results of the interaction characteristics between two harmonic sources can be determined based on the chaotic interaction factor. This can help optimize the power system based on the analysis results, solve the harmonic uncertainty problem caused by nonlinear and dynamic factors in multi-harmonic sources, and thus help improve the stability of the multi-harmonic source system. BRIEF DESCRIPTION OF THE DRAWINGS
[0039] In order to more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the following briefly introduces the drawings required for use in the embodiments of the present application or related technical descriptions. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying any creative work.
[0040] Figure 1 This is a schematic diagram of an application scenario of a harmonic interaction characteristic analysis method based on chaotic characteristics provided in an embodiment of the present application;
[0041] Figure 2This is a flow chart of a harmonic interaction characteristic analysis method based on chaotic characteristics provided by an embodiment of the present application;
[0042] Figure 3 is a schematic diagram of an equivalent circuit model of a transmission line provided in this application;
[0043] Figure 4 is a schematic diagram of an equivalent connection impedance model provided in an embodiment of the present application;
[0044] Figure 5 1 is a flow chart of another harmonic interaction characteristic analysis method based on chaotic characteristics provided in an embodiment of the present application;
[0045] Figure 6 Schematic diagram of a harmonic interaction characteristic analysis device based on chaotic characteristics provided in an embodiment of the present application;
[0046] Figure 7 It is a structural diagram of a computer device provided in an embodiment of the present application. DETAILED DESCRIPTION
[0047] In order to make the purpose, technical solutions and advantages of this application more clear, the following further describes this application in detail with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain this application and are not intended to limit this application.
[0048] The following introduces the application scenarios of the harmonic interaction characteristic analysis method based on chaotic characteristics provided in the embodiments of the present application.
[0049] See Figure 1 , Figure 1 Schematic diagram of an application scenario of a harmonic interaction characteristic analysis method based on chaotic characteristics provided by an embodiment of the present application. Figure 1 As shown, the system includes a computer device 101 and a database server 102, which communicate with each other via a network. The database server 102 stores harmonic impedances corresponding to multiple harmonic sources. The harmonic impedance corresponding to each harmonic source can be determined by the computer device 101 and sent to the database server 102.
[0050] The computer device 101 can obtain the harmonic impedances corresponding to multiple harmonic sources respectively from the database server 102; the first harmonic source and the second harmonic source included in the multiple harmonic sources; based on the pre-constructed equivalent circuit model of the transmission line between the multiple harmonic sources, determine the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source, and based on the series impedance and the parallel admittance, determine the equivalent connection impedance between the first harmonic source and the second harmonic source; the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance and the Lyapunov index are input into the preset chaos interaction factor calculation model to obtain the chaos interaction factor between the first harmonic source and the second harmonic source, and the chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaos effect; based on the chaos interaction factor, determine the interaction characteristic analysis result between the first harmonic source and the second harmonic source. By using this method, by introducing the Lyapunov index, a chaos indicator, the complexity of harmonic interaction between harmonic sources (i.e., the chaotic interaction factor) can be quantitatively evaluated, and the analysis results of the interaction characteristics between two harmonic sources can be determined based on the chaotic interaction factor. This can help optimize the power system based on the analysis results, solve the harmonic uncertainty problem caused by nonlinear and dynamic factors in multi-harmonic sources, and thus help improve the stability of the multi-harmonic source system.
[0051] Optionally, computer device 101 can be a terminal device or a server. The terminal devices mentioned here may include, but are not limited to, various personal computers, laptops, smartphones, tablets, IoT devices, and portable wearable devices. IoT devices may include smart speakers, smart TVs, smart air conditioners, smart car devices, projectors, etc. Portable wearable devices may include smart watches, smart bracelets, head-mounted devices, etc. Head-mounted devices may include virtual reality (VR) devices, augmented reality (AR) devices, smart glasses, etc. The server mentioned here may be an independent physical server, a server cluster or distributed system composed of multiple physical servers, or a cloud server providing cloud computing services, etc., without limitation here.
[0052] See Figure 2 , Figure 2 1 is a flow chart of a method for analyzing harmonic interaction characteristics based on chaotic characteristics provided by an embodiment of the present application. The method can be executed by a computer device (for example, the computer device 100 described above). Figure 2 As shown, the harmonic interaction characteristic analysis method based on chaotic characteristics may include but is not limited to the following steps:
[0053] S201: Acquire harmonic impedances corresponding to a plurality of harmonic sources, including a first harmonic source and a second harmonic source.
[0054] A harmonic source refers to any device or load in the power system that generates harmonic currents or voltages. Harmonic sources include, but are not limited to, inverters, rectifiers, and transformers.
[0055] Among them, harmonic impedance refers to the impedance characteristics of the power system or electrical equipment under the action of harmonic voltage of a specific frequency. It is usually expressed as a complex number, including the combined effects of resistance (R), inductance (L) and capacitance (C).
[0056] S202. Based on a pre-constructed equivalent circuit model of a transmission line between multiple harmonic sources, determine the series impedance and parallel admittance of the transmission line between the first harmonic source and the second harmonic source, and based on the series impedance and parallel admittance, determine the equivalent connection impedance between the first harmonic source and the second harmonic source.
[0057] The transmission line equivalent circuit model is a mathematical model used to simplify the analysis of the electrical characteristics of transmission lines. It converts actual distributed parameters (resistance, inductance, capacitance, and conductance) into lumped parameter circuits. The construction of a transmission line equivalent circuit model depends primarily on the line length and the frequency characteristics of the line unit parameters. Depending on the line length and frequency, transmission line equivalent circuit models can be divided into short-line models, medium-length line models (π-type / T-type), and long-line distributed parameter models.
[0058] In some embodiments, the pre-built equivalent circuit model of the transmission line between multiple harmonic sources can be as follows: Figure 3 As shown, Figure 3 This is a schematic diagram of a transmission line equivalent circuit model provided by this application. Computer equipment is based on Figure 3 , determining the series impedance and parallel admittance of the transmission line between the first harmonic source and the second harmonic source can be done by first determining the line propagation coefficient, and then determining the series impedance and parallel admittance of the transmission line between the first harmonic source and the second harmonic source based on the line propagation coefficient.
[0059] Optionally, the computer device may use the following formula (1) to determine the line propagation coefficient.
[0060] (1)
[0061] In formula (1), h represents the harmonic order; It represents the line propagation coefficient under the hth harmonic; R represents the resistance of the transmission line per unit length; X represents the reactance of the transmission line per unit length; B represents the susceptance of the transmission line per unit length; L represents the line length of the transmission line; j represents the imaginary part of the complex number.
[0062] Optionally, when determining the series impedance of the computer device based on the line propagation coefficient, the following formula (2) may be used.
[0063] (2)
[0064] In formula (2), It represents the series impedance under the hth harmonic; It represents the line propagation coefficient under the hth harmonic, which can be determined by the above formula (1). The physical meanings of other parameters can be found in the explanation of the physical meanings of each parameter in the above formula (1), which will not be repeated here.
[0065] Optionally, when the computer device determines the parallel admittance based on the line propagation coefficient, the following formula (3) may be used.
[0066] (3)
[0067] In formula (3), It represents the parallel admittance under the hth harmonic; It represents the line propagation coefficient under the hth harmonic, which can be determined by the above formula (1). The physical meanings of other parameters can be found in the explanation of the physical meanings of each parameter in the above formula (1), which will not be repeated here.
[0068] In some embodiments, the computer device determines the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance and the parallel admittance, and can determine the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance, the parallel admittance, the capacitive reactance and the inductive reactance of the passive filter between the first harmonic source and the second harmonic source.
[0069] S203. Input the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov index into a preset chaos interaction factor calculation model to obtain a chaos interaction factor between the first harmonic source and the second harmonic source.
[0070] Among them, the chaos interaction factor is used to characterize the influence of the first harmonic source on the second harmonic source under the chaos effect.
[0071] S204. Determine an analysis result of interaction characteristics between the first harmonic source and the second harmonic source based on the chaotic interaction factor.
[0072] In an optional embodiment, the computer device determines the interaction characteristics between the first harmonic source and the second harmonic source based on the chaos interaction factor, which can be: when the chaos interaction factor is greater than 1, the interaction characteristic between the first harmonic source and the second harmonic source is determined to be that the degree of interaction influence of the first harmonic source on the second harmonic source is strong; when the chaos interaction factor is less than 1, the interaction characteristic is determined to be that the degree of interaction influence of the first harmonic source on the second harmonic source is weak; when the chaos interaction factor is equal to 1, the interaction characteristic is determined to be that the first harmonic source has no influence on the second harmonic source.
[0073] Among them, when the interactive influence of the first harmonic source on the second harmonic source is strong, the harmonics of the first harmonic source propagate and amplify in the multi-harmonic source system, thereby affecting other harmonic sources; when the interactive influence of the first harmonic source on the second harmonic source is weak, the harmonics of the first harmonic source gradually attenuate in the multi-harmonic source system and are eventually effectively suppressed; when the first harmonic source has no influence on the second harmonic source, the transmission and attenuation of the harmonics of the first harmonic source are balanced, and the harmonics are in a stable state. At this time, harmonic interaction analysis can be further performed in combination with impedance.
[0074] In an embodiment of the present application, a computer device can obtain harmonic impedances corresponding to multiple harmonic sources respectively; a first harmonic source and a second harmonic source included in the multiple harmonic sources; based on a pre-constructed equivalent circuit model of a transmission line between the multiple harmonic sources, determine the series impedance and parallel admittance of the transmission line between the first harmonic source and the second harmonic source, and based on the series impedance and the parallel admittance, determine the equivalent connection impedance between the first harmonic source and the second harmonic source; the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance and the Lyapunov index are input into a preset chaos interaction factor calculation model to obtain a chaos interaction factor between the first harmonic source and the second harmonic source, and the chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaos effect; based on the chaos interaction factor, determine the interaction characteristic analysis results between the first harmonic source and the second harmonic source. This method, by introducing the Lyapunov index as a chaos indicator, can quantitatively assess the complexity of harmonic interactions between harmonic sources (i.e., the chaotic interaction factor). Based on this chaotic interaction factor, the interaction characteristics analysis results between two harmonic sources can be determined. This analysis can be used to optimize the power system, addressing harmonic uncertainty issues caused by nonlinear and dynamic factors in multi-harmonic sources, thereby improving the stability of multi-harmonic source systems. Furthermore, it can also contribute to improved power quality.
[0075] In an optional embodiment, Figure 2In the harmonic interaction characteristic analysis method based on chaotic characteristics shown, the harmonic impedances corresponding to the multiple harmonic sources can be determined by a computer device through the following steps: for each harmonic source among the multiple harmonic sources, a target harmonic signal is injected into the input end of the harmonic source through a signal generator; the target harmonic signal is determined based on the operating frequency of the harmonic source; the voltage response data and current response data corresponding to the output end of the harmonic source under the target harmonic signal are obtained; the amplitude of the harmonic voltage is determined based on the voltage response data, and the amplitude of the harmonic current is determined based on the current response data; and the harmonic impedance corresponding to the harmonic source is determined based on the amplitude of the harmonic voltage and the amplitude of the harmonic current.
[0076] For example, the operating frequency of the harmonic source may be the 2nd or 3rd order.
[0077] In some embodiments, a computer device determines the harmonic impedance corresponding to a harmonic source based on the amplitude of the harmonic voltage and the amplitude of the harmonic current, which may include: determining a quotient of the amplitude of the harmonic voltage and the amplitude of the harmonic current; determining a difference between a first phase angle of the harmonic voltage and a second phase angle of the harmonic current; and determining the harmonic impedance corresponding to the harmonic source based on the quotient and the difference. The computer device may determine the harmonic impedance corresponding to the harmonic source as shown in the following formula (4).
[0078] (4)
[0079] In formula (4), f represents the injected target harmonic signal; Z(f) represents the harmonic impedance after the target harmonic signal f is injected; V p It represents the amplitude of harmonic voltage; I p It represents the amplitude of harmonic current; It represents the phase angle of harmonic voltage (first phase angle); It represents the phase angle of harmonic current (second phase angle).
[0080] By adopting this embodiment, the computer device can accurately determine the harmonic impedance corresponding to the harmonic source, thereby providing a data basis for subsequent harmonic interaction analysis.
[0081] In an optional embodiment, Figure 2 In the harmonic interaction characteristic analysis method based on chaotic characteristics shown, the computer device determines the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance and the parallel admittance, which may include: determining the input impedance of the transmission line based on the series impedance and the parallel admittance; obtaining the capacitive reactance and the inductive reactance of the passive filter between the first harmonic source and the second harmonic source, and determining the equivalent impedance of the passive filter based on the capacitive reactance and the inductive reactance; determining the equivalent connection impedance between the first harmonic source and the second harmonic source based on the input impedance and the equivalent impedance.
[0082] In some embodiments, a computer device determines the input impedance of a transmission line based on series impedance and shunt admittance, which may include: determining the line characteristic impedance of the transmission line based on the series impedance and shunt admittance; and determining the input impedance of the transmission line based on the line characteristic impedance and the line propagation coefficient.
[0083] Optionally, when the computer device determines the input impedance of the transmission line based on the series impedance and the shunt admittance, the following formula (5) may be used.
[0084] (5)
[0085] In formula (5), It represents the line characteristic impedance of the transmission line under the hth harmonic; It represents the series impedance under the hth harmonic, which can be determined by the above formula (2); It represents the parallel admittance under the hth harmonic and can be determined by the above formula (3).
[0086] Optionally, when the computer device determines the input impedance of the transmission line based on the line characteristic impedance and the line propagation coefficient, the following formula (6) may be used.
[0087] (6)
[0088] In formula (6), It represents the input impedance of the transmission line under the hth harmonic; It represents the line characteristic impedance of the transmission line under the hth harmonic, which can be determined by the above formula (5); It represents the line propagation coefficient under the hth harmonic, which can be determined by the above formula (1); L represents the line length of the transmission line.
[0089] In some embodiments, the capacitive reactance of the passive filter may be determined by a computer device using the following formula (7).
[0090] (7)
[0091] In formula (7), h represents the harmonic order; C F represents the capacitance of the passive filter; ω represents the angular frequency; X FCh It represents the capacitive reactance of the passive filter under the hth harmonic.
[0092] In some embodiments, the inductive reactance of the passive filter may be determined by a computer device using the following formula (8).
[0093] (8)
[0094] In formula (8), h represents the harmonic order; L F represents the inductance of the passive filter; ω represents the angular frequency; X FLh It represents the capacitive reactance of the passive filter under the hth harmonic.
[0095] For example, assuming that the passive filter is a parallel inductor-capacitor type filter, the equivalent impedance of the passive filter under the hth harmonic can be determined by the computer device using the following formula (9).
[0096] (9)
[0097] In formula (9), It represents the equivalent impedance of the passive filter under the hth harmonic; X FCh It represents the capacitive reactance of the passive filter under the hth harmonic, which can be determined by the above formula (7); X FLh represents the capacitive reactance of the passive filter under the hth harmonic, which can be determined by the above formula (8); j represents the imaginary part of the complex number.
[0098] In some embodiments, when the computer device determines the equivalent connection impedance between the first harmonic source and the second harmonic source based on the input impedance and the equivalent impedance, the following formula (10) may be used.
[0099] (6)
[0100] In formula (10), It represents the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source under the h-th harmonic; It represents the equivalent impedance of the passive filter under the hth harmonic, which can be determined by the above formula (9); It represents the equivalent impedance of the passive filter under the hth harmonic, which can be determined by the above formula (6).
[0101] Taking the interaction between the i-th harmonic source and the j-th harmonic source as an example, the computer device can also construct an equivalent connection impedance model between the i-th harmonic source and the j-th harmonic source as follows: Figure 4 shown. Figure 4 In, Z i (h) It represents the harmonic impedance of the i-th harmonic source under the h-th harmonic; Z j (h) It represents the harmonic impedance of the jth harmonic source under the hth harmonic; It represents the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source under the h-th harmonic.
[0102] By adopting this embodiment, the computer device can accurately determine the equivalent connection impedance between the first harmonic source and the second harmonic source, thereby providing data support for the subsequent determination of the chaotic interaction factor between the first harmonic source and the second harmonic source.
[0103] In an optional embodiment, Figure 2 In the harmonic interaction characteristic analysis method based on chaotic characteristics shown in FIG, the preset chaotic interaction factor calculation model can be shown as the following formula (11).
[0104] (11)
[0105] Formula (11), H ij It represents the chaotic interaction factor between the i-th harmonic source and the j-th harmonic source; Z ij (h) It represents the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source under the h-th harmonic; Z i (h) It represents the harmonic impedance of the i-th harmonic source under the h-th harmonic; Z j (h) It represents the harmonic impedance of the jth harmonic source under the hth harmonic; It represents the Lyapunov exponent between the i-th harmonic source and the j-th harmonic source under the h-th harmonic.
[0106] The Lyapunov exponent can be used to represent the sensitivity of the interaction process to disturbances, reflecting the chaotic characteristics of the interaction process. In some embodiments, the Lyapunov exponent can be calculated by a computer device through a numerical simulation method.
[0107] This implementation method introduces the Lyapunov index to quantitatively evaluate the complexity of harmonic interactions between harmonic sources (chaotic interaction factor), and determines the interaction characteristic analysis results between two harmonic sources based on the chaotic interaction factor. This helps predict whether harmonics are amplified or attenuated between multi-harmonic source devices, thereby optimizing power system design and improving power quality and system stability.
[0108] The following combination Figure 5 , the harmonic interaction characteristic analysis method based on chaotic characteristics provided by the embodiment of the present application is specifically described. For the convenience of description, Figure 5 The interaction between two harmonic sources is used as an example. Figure 5 , Figure 5 FIG. 1 is a flow chart of another method for analyzing harmonic interaction characteristics based on chaotic characteristics provided by an embodiment of the present application. Figure 5 As shown, the harmonic interaction characteristic analysis method based on chaotic characteristics may include but is not limited to the following steps:
[0109] S501: Determine a first harmonic impedance corresponding to an i-th harmonic source and a second harmonic impedance corresponding to a j-th harmonic source.
[0110] In some embodiments, the first harmonic impedance and the second harmonic impedance may both be determined by a computer device using the aforementioned formula (4).
[0111] S502: Construct an equivalent circuit model of the transmission line between the i-th harmonic source and the j-th harmonic source, and an equivalent model of the filtering link between the i-th harmonic source and the j-th harmonic source.
[0112] In some embodiments, the equivalent circuit model of the transmission line between the i-th harmonic source and the j-th harmonic source can be found in Figure 3 shown.
[0113] S503. Determine the series impedance and parallel admittance of the transmission line between the i-th harmonic source and the j-th harmonic source based on the equivalent circuit model of the transmission line, and determine the capacitive reactance and inductive reactance of the passive filter between the i-th harmonic source and the j-th harmonic source based on the equivalent model of the filtering link.
[0114] In some embodiments, the series impedance of the transmission line between the i-th harmonic source and the j-th harmonic source may be determined by a computer device using the aforementioned formula (2); and the parallel admittance may be determined by a computer device using the aforementioned formula (3).
[0115] In some embodiments, the capacitive reactance of the passive filter between the i-th harmonic source and the j-th harmonic source may be determined by a computer device using the aforementioned formula (7), and the inductive reactance may be determined by a computer device using the aforementioned formula (8).
[0116] S504 : Determine the input impedance of the transmission line based on the series impedance and the shunt admittance, and determine the equivalent impedance of the passive filter based on the capacitive reactance and the inductive reactance.
[0117] In some embodiments, when the computer device determines the input impedance of the transmission line based on the series impedance and the shunt admittance, the aforementioned formulas (5) and (6) may be used; when the computer device determines the equivalent impedance of the passive filter based on the capacitive reactance and the inductive reactance, the aforementioned formula (9) may be used.
[0118] S505 : Determine the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source based on the input impedance of the transmission line and the equivalent impedance of the passive filter.
[0119] In some embodiments, the computer device may use the aforementioned formula (10) when determining the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source based on the input impedance of the transmission line and the equivalent impedance of the passive filter.
[0120] S506 : Determine the chaotic interaction factor between the i-th harmonic source and the j-th harmonic source based on the first harmonic impedance, the second harmonic impedance, the equivalent connection impedance, and the Lyapunov exponent.
[0121] In some embodiments, when the computer device determines the chaotic interaction factor between the i-th harmonic source and the j-th harmonic source based on the first harmonic impedance, the second harmonic impedance, the equivalent connection impedance, and the Lyapunov index, the aforementioned formula (11) can be used.
[0122] S507 : Determine the analysis result of the interaction characteristics between the i-th harmonic source and the j-th harmonic source based on the chaotic interaction factor.
[0123] In some embodiments, the computer device is based on the chaotic interaction factor (denoted as H ij ), determine the interaction characteristic analysis results between the i-th harmonic source and the j-th harmonic source, which can be shown as follows:
[0124] (1) In H ij When the value is greater than 1, the analysis results show that harmonic source i has a strong interactive influence on harmonic source j, indicating that the harmonics of harmonic source i are propagating and amplifying in the multi-harmonic source system, further affecting other harmonic sources. In this case, users can focus on this interactive process and take timely measures to address it.
[0125] (2) In H ij When <1, the analysis results show that the interaction effect of harmonic source i on harmonic source j is weak, indicating that the harmonics of harmonic source i gradually decay in the multi-harmonic source system and are eventually effectively suppressed.
[0126] (3) In H ij = 1, the analysis result shows that the transmission and attenuation of the harmonics of harmonic source i are balanced, and the harmonics are in a stable state. In this case, further analysis is required in combination with impedance.
[0127] In an embodiment of the present application, the computer device analyzes the interaction characteristics between harmonic sources by introducing the Lyapunov index, a chaos indicator, and can quantitatively evaluate the complexity of harmonic interactions between harmonic sources, thereby solving the harmonic uncertainty problem caused by nonlinear and dynamic factors in multi-harmonic sources, which can be beneficial to improving the stability of the multi-harmonic source system.
[0128] It should be understood that, although the various steps in the flowcharts involved in the various embodiments described above are displayed in sequence according to the instructions of the arrows, these steps are not necessarily executed in sequence in the order indicated by the arrows. Unless otherwise specified herein, there is no strict order restriction on the execution of these steps, and these steps can be executed in other orders. Moreover, at least a portion of the steps in the flowcharts involved in the various embodiments described above can include multiple steps or multiple stages, and these steps or stages are not necessarily executed and completed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a portion of steps or stages in other steps.
[0129] Based on the same inventive concept, an embodiment of the present application further provides a chaotic-characteristics-based harmonic interaction characteristics analysis device for implementing the aforementioned chaotic-characteristics-based harmonic interaction characteristics analysis method. The solution provided by this device is similar to the solution described in the aforementioned method. Therefore, the specific limitations of one or more chaotic-characteristics-based harmonic interaction characteristics analysis device embodiments provided below can be found in the limitations of the chaotic-characteristics-based harmonic interaction characteristics analysis method described above and will not be repeated here.
[0130] See Figure 6 , Figure 6 Schematic diagram of a harmonic interaction characteristic analysis device based on chaotic characteristics provided by an embodiment of the present application. Figure 6 As shown, the harmonic interaction characteristic analysis device based on chaotic characteristics may include but is not limited to:
[0131] An acquisition module 601 is configured to acquire harmonic impedances corresponding to a plurality of harmonic sources, including a first harmonic source and a second harmonic source;
[0132] a determination module 602 for determining the series impedance and parallel admittance of the transmission line between the first harmonic source and the second harmonic source based on a pre-constructed equivalent circuit model of the transmission line between the plurality of harmonic sources, and determining the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance and the parallel admittance;
[0133] Processing module 603 is used to input the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov exponent into a preset chaos interaction factor calculation model to obtain a chaos interaction factor between the first harmonic source and the second harmonic source. The chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect;
[0134] The determination module 602 is further configured to determine an analysis result of an interaction characteristic between the first harmonic source and the second harmonic source based on the chaotic interaction factor.
[0135] In one embodiment, the determination module 602 is also used to inject a target harmonic signal into the input end of the harmonic source through a signal generator for each harmonic source among multiple harmonic sources; the target harmonic signal is determined based on the operating frequency of the harmonic source; the voltage response data and current response data corresponding to the output end of the harmonic source under the target harmonic signal are obtained; the amplitude of the harmonic voltage is determined based on the voltage response data, and the amplitude of the harmonic current is determined based on the current response data; and the harmonic impedance corresponding to the harmonic source is determined based on the amplitude of the harmonic voltage and the amplitude of the harmonic current.
[0136] In one embodiment, when the determination module 602 is used to determine the harmonic impedance corresponding to the harmonic source based on the amplitude of the harmonic voltage and the amplitude of the harmonic current, it is specifically used to: determine the quotient of the amplitude of the harmonic voltage and the amplitude of the harmonic current; determine the difference between the first phase angle of the harmonic voltage and the second phase angle of the harmonic current; and determine the harmonic impedance corresponding to the harmonic source based on the quotient and the difference.
[0137] In one embodiment, when the determination module 602 is used to determine the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance and the parallel admittance, it is specifically used to: determine the input impedance of the transmission line based on the series impedance and the parallel admittance; obtain the capacitive reactance and inductive reactance of the passive filter between the first harmonic source and the second harmonic source, and determine the equivalent impedance of the passive filter based on the capacitive reactance and the inductive reactance; determine the equivalent connection impedance between the first harmonic source and the second harmonic source based on the input impedance and the equivalent impedance.
[0138] In one embodiment, when the determination module 602 is used to determine the input impedance of the transmission line based on the series impedance and the shunt admittance, it is specifically used to: determine the line characteristic impedance of the transmission line based on the series impedance and the shunt admittance; and determine the input impedance of the transmission line based on the line characteristic impedance and the line propagation coefficient.
[0139] In one embodiment, the preset chaotic interaction factor calculation model is as follows:
[0140]
[0141] Among them, H ij It represents the chaotic interaction factor between the i-th harmonic source and the j-th harmonic source; Z ij (h) represents the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source under the h-th harmonic; Z i (h) represents the harmonic impedance of the i-th harmonic source under the h-th harmonic; Z j(h) represents the harmonic impedance of the jth harmonic source under the hth harmonic; It represents the Lyapunov exponent between the i-th harmonic source and the j-th harmonic source under the h-th harmonic.
[0142] Each module in the aforementioned chaotic-based harmonic interaction characteristics analysis device can be implemented in whole or in part through software, hardware, or a combination thereof. Each module can be embedded in or independent of a processor in a terminal device in hardware form, or can be stored in a memory in the terminal device in software form, so that the processor can call and execute the corresponding operations of each module.
[0143] In an exemplary embodiment, a computer device is provided. The computer device may be a terminal, and its internal structure diagram may be as shown in FIG. Figure 7 As shown. The computer device includes a processor, memory, an input / output interface, a communication interface, a display unit, and an input device. The processor, memory, and input / output interface are connected via a system bus, and the communication interface, display unit, and input device are connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and internal memory. The non-volatile storage medium stores an operating system and a computer program. The internal memory provides an environment for the operation of the operating system and computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with external terminals via wired or wireless means, and the wireless means can be implemented via Wi-Fi, a mobile cellular network, near-field communication (NFC), or other technologies. When executed by the processor, the computer program implements a method for analyzing harmonic interaction characteristics based on chaotic characteristics. The display unit of the computer device is used to form a visually visible image, and can be a display screen, a projection device, or a virtual reality imaging device. The display screen can be a liquid crystal display screen or an electronic ink display screen, and the input device of the computer device can be a touch layer covering the display screen, or a button, trackball or touchpad set on the computer device casing, or an external keyboard, touchpad or mouse.
[0144] Those skilled in the art will understand that Figure 7 The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different component arrangement.
[0145] In an exemplary embodiment, the present application provides a computer device including a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the steps in the above-mentioned method for analyzing harmonic interaction characteristics based on chaotic characteristics are implemented.
[0146] In an exemplary embodiment, the present application provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps in the above-mentioned method for analyzing harmonic interaction characteristics based on chaotic characteristics are implemented.
[0147] In an exemplary embodiment, the present application provides a computer program product, including a computer program, which implements the steps in the above-mentioned chaotic characteristic-based harmonic interaction characteristic analysis method when executed by a processor.
[0148] It should be noted that the data involved in this application (including but not limited to acquired data, data used for analysis, etc.) are all information and data authorized by the user or fully authorized by all parties, and the collection, use and processing of relevant data must comply with relevant regulations.
[0149] Those skilled in the art will understand that all or part of the processes in the above-mentioned embodiments can be implemented by instructing the relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. In particular, any reference to memory, database, or other media used in the embodiments provided in this application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The databases involved in the various embodiments provided herein may include at least one of a relational database and a non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the various embodiments provided herein may be, but are not limited to, general-purpose processors, central processing units (CPUs), graphics processing units (GPUs), digital signal processors (DSPs), programmable logic devices (PLDs), quantum computing-based data processing logic devices, artificial intelligence (AI) processors, and the like.
[0150] The technical features of the above embodiments can be combined arbitrarily. In order to make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.
[0151] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present application. It should be noted that a person of ordinary skill in the art may make various modifications and improvements without departing from the spirit of the present application, and these modifications and improvements fall within the scope of protection of the present application. Therefore, the scope of protection of the present application shall be determined by the appended claims.
Claims
1. A harmonic interaction characteristics analysis method based on chaotic characteristics, characterized in that: The method comprises: Acquire harmonic impedances corresponding to a plurality of harmonic sources respectively; the plurality of harmonic sources include a first harmonic source and a second harmonic source; Determining the series impedance and the parallel admittance of the transmission line between the first harmonic source and the second harmonic source based on a pre-constructed equivalent circuit model of the transmission line between the plurality of harmonic sources, and determining the equivalent connection impedance between the first harmonic source and the second harmonic source based on the series impedance and the parallel admittance; Inputting the first harmonic impedance corresponding to the first harmonic source, the second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and the Lyapunov index into a preset chaos interaction factor calculation model to obtain a chaos interaction factor between the first harmonic source and the second harmonic source, wherein the chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under the chaotic effect; Based on the chaotic interaction factor, an analysis result of interaction characteristics between the first harmonic source and the second harmonic source is determined.
2. The method according to claim 1, characterized in that The harmonic impedances corresponding to the multiple harmonic sources are determined in the following manner: For each of the plurality of harmonic sources, injecting a target harmonic signal into an input terminal of the harmonic source through a signal generator; the target harmonic signal is determined based on an operating frequency of the harmonic source; Acquiring voltage response data and current response data corresponding to the output end of the harmonic source under the target harmonic signal; determining a magnitude of a harmonic voltage based on the voltage response data, and determining a magnitude of a harmonic current based on the current response data; A harmonic impedance corresponding to the harmonic source is determined based on the amplitude of the harmonic voltage and the amplitude of the harmonic current.
3. The method according to claim 2, characterized in that The determining the harmonic impedance corresponding to the harmonic source based on the amplitude of the harmonic voltage and the amplitude of the harmonic current includes: determining a quotient of the amplitude of the harmonic voltage and the amplitude of the harmonic current; determining a difference between a first phase angle of the harmonic voltage and a second phase angle of the harmonic current; Based on the quotient and the difference, a harmonic impedance corresponding to the harmonic source is determined.
4. The method according to claim 1, wherein The determining, based on the series impedance and the parallel admittance, an equivalent connection impedance between the first harmonic source and the second harmonic source, includes: determining an input impedance of the transmission line based on the series impedance and the shunt admittance; Obtaining capacitive reactance and inductive reactance of a passive filter between the first harmonic source and the second harmonic source, and determining an equivalent impedance of the passive filter based on the capacitive reactance and the inductive reactance; An equivalent connection impedance between the first harmonic source and the second harmonic source is determined based on the input impedance and the equivalent impedance.
5. The method according to claim 4, characterized in that The determining the input impedance of the transmission line based on the series impedance and the shunt admittance includes: determining a line characteristic impedance of the transmission line based on the series impedance and the shunt admittance; The input impedance of the power transmission line is determined based on the line characteristic impedance and the line propagation coefficient.
6. The method according to claim 1, characterized in that The preset chaotic interaction factor calculation model is as follows: Among them, H ij It represents the chaotic interaction factor between the i-th harmonic source and the j-th harmonic source; Z ij (h) represents the equivalent connection impedance between the i-th harmonic source and the j-th harmonic source under the h-th harmonic; Z i (h) represents the harmonic impedance of the i-th harmonic source under the h-th harmonic; Z j (h) represents the harmonic impedance of the jth harmonic source under the hth harmonic; It represents the Lyapunov exponent between the i-th harmonic source and the j-th harmonic source under the h-th harmonic.
7. A harmonic interaction characteristic analysis device based on chaotic characteristics, characterized in that: The device comprises: An acquisition module, configured to acquire harmonic impedances corresponding to a plurality of harmonic sources, wherein the plurality of harmonic sources include a first harmonic source and a second harmonic source; a determination module, configured to determine, based on a pre-constructed equivalent circuit model of a transmission line between the plurality of harmonic sources, a series impedance and a parallel admittance of the transmission line between the first harmonic source and the second harmonic source, and determine, based on the series impedance and the parallel admittance, an equivalent connection impedance between the first harmonic source and the second harmonic source; a processing module, configured to input a first harmonic impedance corresponding to the first harmonic source, a second harmonic impedance corresponding to the second harmonic source, the equivalent connection impedance, and a Lyapunov exponent into a preset chaos interaction factor calculation model to obtain a chaos interaction factor between the first harmonic source and the second harmonic source, wherein the chaos interaction factor is used to characterize the degree of influence of the first harmonic source on the second harmonic source under a chaotic effect; The determination module is further configured to determine an analysis result of interaction characteristics between the first harmonic source and the second harmonic source based on the chaotic interaction factor.
8. A computer device, characterized in that: The method comprises a memory and a processor, wherein the memory stores a computer program, and the processor implements the steps of the method according to any one of claims 1 to 6 when executing the computer program.
9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.
10. A computer program product comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.