A new energy unit impedance correction method and system based on parameter fitting

Abstract

The application discloses a new energy unit impedance correction method and system based on parameter fitting, and belongs to the technical field of impedance measurement error correction of new energy units. The method comprises the following steps: measuring the impedance of a new energy unit, obtaining impedance amplitude data arrays and impedance phase data arrays under each frequency; calculating the impedance amplitude data arrays and the impedance phase data arrays to obtain impedance real part data arrays and impedance imaginary part data arrays; fitting the equation coefficients of the impedance real part data arrays and the equation coefficients of the impedance imaginary part data arrays to obtain fitted impedance real parts and impedance imaginary parts under each frequency; and calculating corrected impedance amplitudes and impedance phases under each frequency based on the fitted impedance real parts and the fitted impedance imaginary parts. The application can effectively suppress random errors caused by factors such as sampling random errors and harmonic interference on the impedance measurement results.

Description

Technical Field

[0001] This invention relates to the field of impedance measurement error correction technology for new energy power units, and more specifically, to a method and system for impedance correction of new energy power units based on parameter fitting. Background Technology

[0002] The large-scale penetration of new energy power generation systems, such as wind and solar power, into the power system has brought significant challenges to the stable operation of the power system and is prone to causing harmonic oscillations. Furthermore, the control strategies for new energy generator units are complex, with each manufacturer employing its own control methods, and these methods involve multiple control parameters, all of which can affect the stable operation of the new energy generator units under different grid operating conditions. Therefore, how to assess the stability of new energy generator units is a problem that urgently needs to be solved.

[0003] Impedance stability theory can analyze and evaluate the stable operation capability of renewable energy generating units based on their externally exhibited impedance characteristics. Even if the internal control strategies and parameters of the renewable energy generating unit are unknown, its impedance can be measured by injecting external disturbances. However, the impedance measurement results of renewable energy generating units are affected by random factors such as sampling errors and harmonics, leading to uncertain errors in the amplitude and phase of the directly measured impedance. Reasonably correcting these factors for impedance measurement errors in renewable energy generating units is of significant value for impedance measurement and stable operation capability analysis. Summary of the Invention

[0004] To address the above problems, this invention proposes a method for impedance correction of new energy generating units based on parameter fitting, comprising:

[0005] The impedance of the new energy unit is measured to obtain impedance amplitude data arrays and impedance phase data arrays at various frequencies;

[0006] The impedance magnitude data array and the impedance phase data array are calculated to obtain the impedance real part data array and the impedance imaginary part data array;

[0007] The coefficients of the impedance real part data array equation and the impedance imaginary part data array equation are fitted to obtain the fitted real part and imaginary part of the impedance at each frequency.

[0008] Based on the fitted real and imaginary parts of the impedance, the corrected impedance amplitude and impedance phase at each frequency are calculated.

[0009] Optionally, calculations performed on the impedance magnitude data array and the impedance phase data array include: sine calculations and cosine calculations.

[0010] Optionally, the method also includes:

[0011] The impedance angular frequency is standardized, as are the real and imaginary impedance data arrays.

[0012] Optionally, fitting is performed on the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation, including:

[0013] An impedance transfer function is established by substituting the standardized impedance angular frequency, real part data array, and imaginary part data array into the transfer function, and extracting the equations for the real part data array and the imaginary part data array respectively. The coefficients of the impedance real part data array and the impedance imaginary part data array are fitted using the elements in the impedance real part data array and the impedance imaginary part data array as dependent variables.

[0014] Furthermore, this invention also proposes an impedance correction system for new energy generating units based on parameter fitting, comprising:

[0015] The measurement unit is used to measure the impedance of the new energy unit and obtain impedance amplitude data arrays and impedance phase data arrays at various frequencies.

[0016] The calculation unit is used to calculate the impedance magnitude data array and the impedance phase data array to obtain the impedance real part data array and the impedance imaginary part data array.

[0017] The fitting unit is used to fit the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation to obtain the fitted real part and imaginary part of the impedance at each frequency.

[0018] The correction unit is used to calculate the corrected impedance amplitude and impedance phase at each frequency based on the fitted real and imaginary parts of the impedance.

[0019] Optionally, calculations performed on the impedance magnitude data array and the impedance phase data array include: sine calculations and cosine calculations.

[0020] Optionally, the fitting unit is also used for:

[0021] The impedance angular frequency is standardized, as are the real and imaginary impedance data arrays.

[0022] Optionally, fitting is performed on the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation, including:

[0023] An impedance transfer function is established by substituting the standardized impedance angular frequency, real part data array, and imaginary part data array into the transfer function, and extracting the equations for the real part data array and the imaginary part data array respectively. The coefficients of the impedance real part data array and the impedance imaginary part data array are fitted using the elements in the impedance real part data array and the impedance imaginary part data array as dependent variables.

[0024] In another aspect, the present invention also provides a computing device, comprising: one or more processors;

[0025] A processor is used to execute one or more programs;

[0026] When the one or more programs are executed by the one or more processors, the method described above is implemented.

[0027] In another aspect, the present invention also provides a computer-readable storage medium having a computer program stored thereon, which, when executed, implements the method described above.

[0028] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0029] This invention provides an impedance correction method for new energy power generation units based on parameter fitting, comprising: measuring the impedance of the new energy power generation unit to obtain impedance amplitude data arrays and impedance phase data arrays at various frequencies; calculating the impedance amplitude data arrays and impedance phase data arrays to obtain impedance real part data arrays and impedance imaginary part data arrays; fitting the equation coefficients of the impedance real part data arrays and impedance imaginary part data arrays to obtain the fitted impedance real part and impedance imaginary part at various frequencies; and calculating the corrected impedance amplitude and impedance phase at various frequencies based on the fitted impedance real part and impedance imaginary part. This invention can effectively suppress random errors in impedance measurement results caused by factors such as sampling random errors and harmonic interference. Attached Figure Description

[0030] Figure 1 This is a flowchart of the method of the present invention;

[0031] Figure 2 This is a flowchart of an embodiment of the method of the present invention;

[0032] Figure 3 The present invention provides the results of d-axis impedance correction for a photovoltaic power generation system inverter in an embodiment of the method.

[0033] Figure 4 This is the q-axis impedance correction result of the photovoltaic power generation system inverter in the embodiment of the method of the present invention;

[0034] Figure 5This is the positive sequence impedance correction result of the inverter in the direct-drive wind turbine power generation system in the embodiment of the method of the present invention;

[0035] Figure 6 This is a structural diagram of the system of the present invention. Detailed Implementation

[0036] Exemplary embodiments of the invention will now be described with reference to the accompanying drawings. However, the invention may be embodied in many different forms and is not limited to the embodiments described herein. These embodiments are provided to fully and completely disclose the invention and to fully convey its scope to those skilled in the art. The terminology used in the exemplary embodiments illustrated in the drawings is not intended to limit the invention. In the drawings, the same units / elements are referred to by the same reference numerals.

[0037] Unless otherwise stated, the terms used herein (including technical terms) have their common meaning as understood by one of ordinary skill in the art. Furthermore, it is understood that terms defined in commonly used dictionaries should be understood to have a meaning consistent with the context of their relevant field, and not to be interpreted as having an idealized or overly formal meaning.

[0038] Example 1:

[0039] This invention proposes an impedance correction method for new energy generating units based on parameter fitting, such as... Figure 1 As shown, it includes:

[0040] Step 1: Measure the impedance of the new energy unit and obtain the impedance amplitude data array and impedance phase data array at each frequency;

[0041] Step 2: Calculate the impedance amplitude data array and impedance phase data array to obtain the impedance real part data array and impedance imaginary part data array;

[0042] Step 3: Fit the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation to obtain the fitted real part and imaginary part of the impedance at each frequency;

[0043] Step 4: Based on the fitted real and imaginary parts of the impedance, calculate the corrected impedance amplitude and impedance phase at each frequency.

[0044] The calculations performed on the impedance amplitude data array and the impedance phase data array include: sine calculation and cosine calculation.

[0045] The methods also include:

[0046] The impedance angular frequency is standardized, as are the real and imaginary impedance data arrays.

[0047] The fitting of the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation includes:

[0048] An impedance transfer function is established by substituting the standardized impedance angular frequency, real part data array, and imaginary part data array into the transfer function, and extracting the equations for the real part data array and the imaginary part data array respectively. The coefficients of the impedance real part data array and the impedance imaginary part data array are fitted using the elements in the impedance real part data array and the impedance imaginary part data array as dependent variables.

[0049] The invention will be further illustrated below with specific examples:

[0050] The theoretical basis of this invention is as follows:

[0051] The impedance of a new energy unit can be solved in a rotating coordinate system, and the impedance measurement results are generally as follows:

[0052]

[0053] In the above formula, Z dd (s) represents the d-axis impedance, Z qq (s) represents the q-axis impedance, Z dq (s) and Z qd (s) represents the cross impedance form of the d-axis and q-axis.

[0054] The impedance of new energy generating units can also be solved in a positive / negative sequence coordinate system, and the impedance measurement results are generally as follows:

[0055]

[0056] In the above formula, Z pp (s) represents the positive sequence impedance, Z nn (s) represents the negative sequence impedance, Z pn (s) and Z np (s) represents the positive and negative sequence cross impedances.

[0057] The first step is to construct an array Z from the impedance amplitude data of the new energy generator units at various frequencies obtained through measurement. mag (ω k ) and phase data constitute an array A pha (ω k A second calculation is performed to obtain an array Z consisting of the real part of the impedance. imp_real (ω k The array Z consisting of the imaginary part and the data of the imaginary part. imp_imag (ω k The frequency information of the measured impedance is represented by ω. kThis represents the angular frequency array of the measured impedance data, which consists of a series of angular frequencies ω. k0 ~ω ka The composition, where 'a' represents the number of frequencies contained in the frequency array. Here, frequency f can also be used. k This indicates the frequency information of the measured impedance, and there is a corresponding relationship between the angular frequency and the frequency ω. k =2πf k The aforementioned secondary calculation method is as follows:

[0058] Z imp_real (ω k ) = Z mag cos[A pha (ω k )]

[0059] Z imp_imag (ω k ) = Z mag sin[A pha (ω k )]

[0060] In the above formula, cos and sin are used to calculate the cosine and sine, respectively.

[0061] The second step involves using an array Z composed of the real part of the impedance data obtained from the second calculation. imp_real (ω k The array Z consisting of the imaginary part and the data of the imaginary part. imp_imag (ω k When fitting the equation consisting of polynomials, the data was standardized.

[0062] The data standardization process involves standardizing each angular frequency. The specific standardization method is as follows:

[0063]

[0064] In the above formula, ω sd This represents the standardized angular frequency array, Mean represents the average value of the corresponding array, and std represents the standard deviation of the corresponding array.

[0065] Data standardization processing method for array Z imp_real (ω k ) and array Z imp_imag (ω k The standardization method is as follows:

[0066]

[0067] In the above formula, Z sd_real (ω k Z is an array consisting of the standardized real parts of the data.sd_imag (ω k ) is an array consisting of the imaginary part of the data after standardization. Max represents the maximum value of the response array, and Min represents the minimum value of the response array.

[0068] Thirdly, considering that impedance is a transfer function, the impedance Z in the rotating coordinate system... dd (s), Z qq (s), Z dq (s) and Z qd (s), and impedance Z in positive / negative sequence coordinates. pp (s), Z nn (s), Z pn (s) and Z np (s) can all be uniformly expressed in the complex frequency domain as:

[0069]

[0070] In the above formula, Z imp (s) is the transfer function of the impedance of the new energy unit, b0~b n c0~c n-1 Let be the correlation coefficient, s be the Laplace operator, and s = jω, where j is the imaginary operator and ω is the angular frequency variable.

[0071] Substitute s = jω into Z imp (s), and by extracting the real and imaginary parts of the equation, we can obtain the angular frequency array, real part data array, and imaginary part data array used for fitting. The equation is composed of high-order polynomials, and the equation is as follows:

[0072]

[0073] In the above formula, Z ce_real (ω) and Z ce_imag (ω) represents the real and imaginary parts of the impedance at the m-th angular frequency position obtained through polynomial calculation. km For the standardized array ω sd The m-th element in p r0 ~p rn q r0 ~q rn p i0 ~p in q i0 ~q in-1 , respectively, are the coefficients of the polynomial in the equation, and n is the order of the polynomial.

[0074] Using the angular frequency array ω sd The elements in the array Z are the independent variables. sd_real (ω k ) and array Zsd_imag (ω k The elements in ) are the dependent variables, and the coefficients p in the polynomial equation are... r0 ~p rn q r0 ~q rn p i0 ~p in q i0 ~q in-1 Perform a fitting operation and obtain the fitting results for these coefficients.

[0075] The fourth step involves using the obtained equation composed of higher-order polynomials and the fitted coefficients p. r0 ~p rn q r0 ~q rn p i0 ~p in q i0 ~q in-1 Calculate the real part Z of the impedance at each frequency. ce_real (ω) and the imaginary part Z ce_imag (ω), and reconstruct the real and imaginary parts for each frequency.

[0076] Based on the real and imaginary parts of the impedance obtained in the previous step, the amplitude and phase of the impedance at each frequency are calculated one by one. The method for calculating the impedance amplitude and phase is as follows:

[0077]

[0078]

[0079] In the above formula, Zmagc(ωkm) and Amagc(ωkm) are the impedance amplitude and phase information at each frequency after error correction.

[0080] The results of impedance correction of the d-axis and q-axis impedances of a photovoltaic inverter in a rotating coordinate system using the method of this invention are as follows: Figure 3 and Figure 4 As shown, the result of correcting the positive sequence impedance of the direct-drive wind turbine inverter using the method of the present invention is as follows: Figure 5 As shown.

[0081] The impedance of the new energy unit measured by this invention, according to the method provided by this invention, can effectively suppress random errors caused by factors such as sampling random errors and harmonic interference on the impedance measurement results, and can obtain specific equations characterizing the real and imaginary parts of the impedance under specific operating conditions of the new energy unit. It provides a convenient interpolation method for impedance measurement under actual discontinuous frequencies, and can directly calculate the impedance of the new energy unit at unmeasured frequencies, thereby helping to build a database of new energy units.

[0082] Example 2:

[0083] This invention also proposes a new energy unit impedance correction system 200 based on parameter fitting, such as... Figure 6 As shown, it includes:

[0084] The measurement unit 201 is used to measure the impedance of the new energy unit and obtain impedance amplitude data array and impedance phase data array at various frequencies;

[0085] The calculation unit 202 is used to calculate the impedance amplitude data array and the impedance phase data array to obtain the impedance real part data array and the impedance imaginary part data array.

[0086] Fitting unit 203 is used to fit the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation to obtain the fitted impedance real part and impedance imaginary part at each frequency.

[0087] The correction unit 204 is used to calculate the corrected impedance amplitude and impedance phase at each frequency based on the fitted real part and imaginary part of the impedance.

[0088] The calculations performed on the impedance amplitude data array and the impedance phase data array include: sine calculation and cosine calculation.

[0089] Among them, the fitting unit 203 is also used for:

[0090] The impedance angular frequency is standardized, as are the real and imaginary impedance data arrays.

[0091] The fitting of the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation includes:

[0092] An impedance transfer function is established by substituting the standardized impedance angular frequency, real part data array, and imaginary part data array into the transfer function, and extracting the equations for the real part data array and the imaginary part data array respectively. The coefficients of the impedance real part data array and the impedance imaginary part data array are fitted using the elements in the impedance real part data array and the impedance imaginary part data array as dependent variables.

[0093] This invention can effectively suppress random errors in impedance measurement results caused by factors such as sampling random errors and harmonic interference.

[0094] Example 3:

[0095] Based on the same inventive concept, this invention also provides a computer device, which includes a processor and a memory. The memory stores a computer program, which includes program instructions. The processor executes the program instructions stored in the computer storage medium. The processor may be a Central Processing Unit (CPU), or other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. It is the computing and control core of the terminal, suitable for implementing one or more instructions, specifically suitable for loading and executing one or more instructions in the computer storage medium to implement corresponding method flows or corresponding functions, thereby implementing the steps of the methods in the above embodiments.

[0096] Example 4:

[0097] Based on the same inventive concept, this invention also provides a storage medium, specifically a computer-readable storage medium (Memory), which is a memory device in a computer device used to store programs and data. It is understood that the computer-readable storage medium here can include both the built-in storage medium in the computer device and extended storage media supported by the computer device. The computer-readable storage medium provides storage space that stores the terminal's operating system. Furthermore, this storage space also stores one or more instructions suitable for loading and execution by a processor. These instructions can be one or more computer programs (including program code). It should be noted that the computer-readable storage medium here can be high-speed RAM or non-volatile memory, such as at least one disk storage device. The processor can load and execute one or more instructions stored in the computer-readable storage medium to implement the steps of the method in the above embodiments.

[0098] Those skilled in the art will understand that embodiments of the present invention can be provided as methods, systems, or computer program products. Therefore, the present invention can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention can take the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code. The solutions in the embodiments of the present invention can be implemented using various computer languages, such as the object-oriented programming language Java and the interpreted scripting language JavaScript.

[0099] This 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 block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can 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, generate instructions for implementing the flowchart illustrations and / or block diagrams. Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0100] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0101] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0102] Although preferred embodiments of the invention have been described, those skilled in the art, upon learning the basic inventive concept, can make other changes and modifications to these embodiments. Therefore, the appended claims are intended to be interpreted as including both the preferred embodiments and all changes and modifications falling within the scope of the invention.

[0103] Obviously, those skilled in the art can make various modifications and variations to this invention without departing from its spirit and scope. Therefore, if these modifications and variations fall within the scope of the claims of this invention and their equivalents, this invention also intends to include these modifications and variations.

Claims

1. A method for impedance correction of new energy generating units based on parameter fitting, characterized in that, The method includes: The impedance of the new energy unit is measured to obtain impedance amplitude and phase data arrays at various frequencies. The impedance of the new energy unit is solved in a positive / negative sequence coordinate system, and the impedance measurement results are as follows: Z pp (s) represents the positive sequence impedance, Z nn (s) represents the negative sequence impedance, Z pn (s) and Z np (s) represents the positive and negative sequence cross impedances; The impedance amplitude data array and impedance phase data array are calculated to obtain the impedance real part data array and impedance imaginary part data array: the array Z is composed of the impedance amplitude data of the new energy unit at each frequency obtained by measurement. mag (ω k ) and phase data constitute an array A pha (ω k A second calculation is performed to obtain an array Z consisting of the real part of the impedance. imp_real (ω k The array Z consisting of the imaginary part and the data of the imaginary part. imp_imag (ω k The frequency information of the measured impedance is represented by ω. k This indicates that the angular frequency array of the measured impedance data consists of a series of angular frequencies ω k0 ~ω ka Composition, where 'a' represents the number of frequencies contained in the frequency array, and the aforementioned secondary calculation method is as follows: , cos and sin are used for calculating cosine and sine, respectively. Fitting the equation coefficients of the real and imaginary impedance data arrays yields the fitted real and imaginary impedance parts at various frequencies. The array Z, composed of the real impedance data obtained through quadratic calculations, is then used. imp_real (ω k The array Z consisting of the imaginary part and the data of the imaginary part. imp_imag (ω k When fitting the polynomial equation, the data was standardized. The data standardization method standardized each angular frequency. The standardization method is as follows: , ω sd This represents the standardized angular frequency array, Mean represents the mean of the corresponding array, and std represents the standard deviation of the corresponding array. For array Z... imp_real (ω k ) and array Z imp_imag (ω k The standardization method is as follows: Z sd_real (ω k Z is an array consisting of the standardized real parts of the data. sd_imag (ω k The array consists of the standardized imaginary parts of the response array; Max represents the maximum value of the response array, and Min represents the minimum value of the response array; the impedance transfer function and the impedance Z in positive / negative order coordinate systems. pp (s), Z nn (s), Z pn (s) and Z np (s) can be uniformly expressed in the complex frequency domain as: Z imp (s) is the transfer function of the impedance of the new energy unit, b0~b n c0~c n-1 Let s be the correlation coefficient, s be the Laplace operator, and s = jω, where j is the imaginary operator and ω is the angular frequency variable; Substituting s = jω into Z imp (s), and by extracting the real and imaginary parts of the equation, we can obtain the equation composed of higher-order polynomials used to fit the angular frequency array, real part data, and imaginary part data array, as follows: Z ce_real (ω) and Z ce_imag (ω) represents the real and imaginary parts of the impedance at the m-th angular frequency obtained through polynomial calculation. km For the standardized array ω sd The m-th element in p r0 ~p rn q r0 ~q rn p i0 ~p in q i0 ~q in-1 These are the coefficients of the polynomial in the equation, where n is the order of the polynomial, and ω is the angular frequency array. sd The elements in the array are independent variables, respectively, represented by array Z. sd_real (ω k ) and array Z sd_imag (ω k The elements in ) are the dependent variables, and the coefficients p in the polynomial equation are... r0 ~p rn q r0 ~q rn p i0 ~p in q i0 ~q in-1 Perform a fit and obtain the fitting results for these coefficients; Based on the fitted real and imaginary parts of the impedance, the corrected impedance amplitude and phase at each frequency are calculated: using the obtained equation composed of higher-order polynomials and the fitted coefficients p... r0 ~p rn q r0 ~q rn p i0 ~p in q i0 ~q in-1 Calculate the real part Z of the impedance at each frequency. ce_real (ω) and the imaginary part Z ce_imag (ω), and reconstruct the real and imaginary parts of the impedance for each frequency. Based on the real and imaginary parts of the impedance obtained in the previous step, calculate the amplitude and phase of the impedance at each frequency one by one. The method for calculating the impedance amplitude and phase is as follows: , In the above formula, Zmagc(ωkm) and Amagc(ωkm) are the impedance amplitude and phase information at each frequency after error correction.

2. A parameter fitting-based impedance correction system for new energy generating units, used to implement the impedance correction method in claim 1, characterized in that, The system includes: The measurement unit is used to measure the impedance of the new energy unit and obtain impedance amplitude data arrays and impedance phase data arrays at various frequencies. The calculation unit is used to calculate the impedance magnitude data array and the impedance phase data array to obtain the impedance real part data array and the impedance imaginary part data array. The fitting unit is used to fit the coefficients of the impedance real part data array equation and the impedance imaginary part data array equation to obtain the fitted real part and imaginary part of the impedance at each frequency. The correction unit is used to calculate the corrected impedance amplitude and impedance phase at each frequency based on the fitted real and imaginary parts of the impedance.

3. The impedance correction system for new energy generating units according to claim 2, characterized in that, The fitting unit is also used for: The impedance angular frequency is standardized, as are the real and imaginary impedance data arrays.

4. A computer device, characterized in that, include: One or more processors; A processor is used to execute one or more programs; When the one or more programs are executed by the one or more processors, the method of claim 1 is implemented.

5. A computer-readable storage medium, characterized in that, It contains a computer program, which, when executed, implements the method as described in claim 1.

Images