Pseudo-random signal construction method and system based on arbitrary square wave controllable phase modulation superposition

Abstract

The present invention belongs to the field of geophysical electromagnetic exploration technology, and provides a pseudo-random signal construction method and system based on arbitrary square wave controllable phase modulation superposition. The present invention can obtain a pseudo-random signal containing any specified frequency components and the amplitudes of these specified frequency components are basically consistent by phase adjustment and superposition of square waves containing different frequency components, and hard limiting of the superimposed waveform. In the implementation process of the present invention, the number of calculations required for a single cycle is relatively small. Compared with the signal construction method that needs to solve complex equations or perform large-scale optimization, the total time consumption of the present invention is likely to be shorter.

Description

Technical Field

[0001] The present invention belongs to the technical field of geophysical electromagnetic exploration, and in particular relates to a pseudo-random signal construction method and system based on arbitrary square wave controllable phase modulation superposition. Background Art

[0002] The statements in this section merely provide background information related to the present invention and do not necessarily constitute prior art.

[0003] Frequency domain electromagnetic method is a type of geophysical exploration method widely used in the fields of mineral exploration, oil and gas, groundwater exploration and geological structure exploration. It is a geophysical exploration method based on the differences in conductivity, magnetic permeability and dielectricity of rocks, minerals (stones) and surrounding media in the earth's crust. It uses the principle of electromagnetic induction to observe and study the frequency distribution law of artificial or naturally formed electromagnetic fields, and then solve various geological problems. It is also called frequency domain electromagnetic induction method. Frequency domain electromagnetic method uses artificial or natural alternating electromagnetic fields as the source. When the alternating electromagnetic field propagates in the underground medium in the form of waves, due to electromagnetic induction, the observed value of the ground electromagnetic field will contain information on the resistivity distribution of the underground medium. At the same time, due to the skin effect of the electromagnetic field, electromagnetic field signals of different periods have different penetration depths. Therefore, by studying the frequency response of the earth to the electromagnetic field, the distribution law of rock resistivity at different depths underground can be obtained. According to the source of the field, there are two types: passive source frequency domain electromagnetic method and active source frequency domain electromagnetic method. The former includes magnetotelluric method, very low frequency electromagnetic method, etc., and the latter includes controlled source audio magnetotelluric method, wide area electromagnetic method, etc.

[0004] In the active source frequency domain electromagnetic method, an artificial alternating electromagnetic field is required as the field source, and a loop is formed by a generator, a long wire, a grounding electrode and the earth between the grounding electrodes. Using an alternating current with specific spectral characteristics in the loop to excite the alternating electromagnetic field is a method often used in actual work. In order to ensure the effect of exploration and improve the efficiency of exploration, the current waveform (or the transmitted signal waveform) emitted in the field source loop is very important. At present, the research on the transmitted waveform is mostly focused on how to increase the available frequency components in the waveform, or how to generate a waveform with specified spectral characteristics more quickly, but there is still a lack of research on the arbitrary customization of waveform spectral characteristics. Summary of the invention

[0005] In order to solve the above problems, the present invention proposes a pseudo-random signal construction method and system based on controllable phase-modulated superposition of arbitrary square waves. The present invention can obtain a pseudo-random signal containing any specified frequency components and whose amplitudes are basically consistent by phase adjustment and superposition of square waves containing different frequency components and hard limiting of the superimposed waveform.

[0006] According to some embodiments, the present invention adopts the following technical solutions:

[0007] A method for constructing a pseudo-random signal based on arbitrary square wave controllable phase modulation superposition includes the following steps:

[0008] (1) Setting the parameters of the waveform to be generated, including the frequency group, waveform length, and sampling rate;

[0009] (2) Construct a square wave with the lowest frequency in the frequency group using the initial phase as the initial waveform;

[0010] (3) Perform a fast Fourier transform on the initial waveform to obtain its spectrum and then obtain the coefficient of the specified frequency;

[0011] (4) constructing the phase of the next square wave according to the angle and amplitude of the coefficient and the ratio of the lowest frequency amplitude in the frequency group, and constructing a square wave with a base frequency of the specified frequency with the phase;

[0012] (5) Add the constructed square wave to the current waveform to obtain a new waveform;

[0013] (6) Repeat steps (3) to (5) until all frequencies in the frequency group are added to the new waveform;

[0014] (7) After hard limiting the final new waveform, a fast Fourier transform is performed to obtain its spectrum, and the coefficients of all frequencies of the corresponding frequency group are taken out and the amplitudes are calculated to obtain an array. Based on the array, the coefficient of variation and the theoretical energy ratio are calculated;

[0015] (8) If the coefficient of variation and the theoretical energy ratio meet the requirements, the final new waveform is output; otherwise, starting from the second frequency in the frequency group, steps (3) to (7) are repeated until the coefficient of variation and the theoretical energy ratio meet the requirements.

[0016] As an optional implementation, in step (1), the frequencies in the frequency group are set arbitrarily, the total number is an odd number, and the frequencies are arranged from low to high.

[0017] As an optional implementation, in step (1), the sampling rate is at least twice the highest frequency in the frequency group, and the waveform length is the ratio of the sampling rate to the greatest common divisor of all frequencies in the frequency group.

[0018] As an optional implementation, in step (2), the initial phase may be any phase.

[0019] As an optional implementation, in step (4), the process of constructing the phase of the next square wave according to the ratio of the angle and amplitude of the coefficient to the lowest frequency amplitude in the frequency group includes: ;

[0020] in, is the angle of the coefficient, is the ratio of the amplitude of the coefficient to the amplitude of the lowest frequency in the frequency group, is the phase of the next square wave.

[0021] As an optional implementation, in step (5), the process of constructing a square wave includes: creating Array of rows and 2 columns , used to save the position index of the square wave with a value of -1, Determined by the following formula: ;

[0022] right No. The first row is filled with , the second column is filled with elements , Represents rounding to the nearest integer;

[0023] Will All elements plus ;

[0024] Create a one-dimensional array with all values ​​1 , whose length is The maximum value in the greater of

[0025] right Each row of , starting with the first column element, ending with the second column element, and using 1 as the interval to construct an index sequence. The value of the corresponding index position is set to -1;

[0026] From a one-dimensional array Delete the beginning of the length that exceeds The required square wave is obtained;

[0027] in, is the waveform length, is the sampling rate, is the calculated phase, is the fundamental frequency, is the total number of frequencies in the frequency group.

[0028] As an optional implementation, in step (7), the process of hard limiting the final new waveform includes: setting the part of the final new waveform that is greater than 1 to 1, and setting the part that is less than -1 to -1.

[0029] As an optional implementation, in step (7), the waveform Perform fast Fourier transform, obtain its spectrum, take out the coefficients of all frequencies in the corresponding frequency group and calculate the amplitude to get an array , is the total number of frequencies in the frequency group, and its coefficient of variation is calculated by the following formula and theoretical energy ratio :

[0030] ;

[0031] ;

[0032] in For array The coefficient of variation reflects the unevenness of these amplitudes, and the theoretical energy ratio reflects the ratio of the average value of these amplitudes to the theoretical maximum average amplitude.

[0033] A pseudo-random signal construction system based on arbitrary square wave controllable phase modulation superposition, comprising:

[0034] A parameter configuration module, configured to set the parameters of the waveform to be generated, the parameters including frequency group, waveform length, and sampling rate;

[0035] An initial construction module is configured to construct a square wave with a fundamental frequency as the lowest frequency in the frequency group with an initial phase as an initial waveform;

[0036] A transform module is configured to perform a fast Fourier transform on the initial waveform to obtain its spectrum and then obtain coefficients of a specified frequency;

[0037] A parameter update calculation module is configured to construct a phase of the next square wave according to the angle and amplitude of the coefficient and the ratio of the lowest frequency amplitude in the frequency group, and to construct a square wave with a base frequency of the specified frequency with the phase;

[0038] A waveform superposition module is configured to add the constructed square wave to the current waveform to obtain a new waveform;

[0039] An iteration module is configured to repeatedly call the transformation module, the parameter update calculation module and the waveform superposition module until all frequencies in the frequency group are added to the new waveform;

[0040] The evaluation index calculation module is configured to perform hard limiting on the final new waveform, then perform fast Fourier transform, obtain its spectrum, extract the coefficients of all frequencies of the corresponding frequency group and calculate the amplitude to obtain an array, and calculate the coefficient of variation and theoretical energy proportion according to the array;

[0041] The judgment module is configured to output the final new waveform if the coefficient of variation and the theoretical energy ratio meet the requirements, otherwise, starting from the second frequency in the frequency group, repeatedly call the transformation module to the evaluation index calculation module until the coefficient of variation and the theoretical energy ratio meet the requirements.

[0042] A computer-readable storage medium is used to store computer instructions, and when the computer instructions are executed by a processor, the steps in the above method are completed.

[0043] An electronic device comprises a memory and a processor, and computer instructions stored in the memory and executed on the processor. When the computer instructions are executed by the processor, the steps in the above method are completed.

[0044] Compared with the prior art, the present invention has the following beneficial effects:

[0045] (1) The frequencies in the given frequency group when constructing the signal in the present invention can be arbitrary and do not need to satisfy the specified multiple or distribution relationship. The only restriction on the frequency group in the present invention is that the total number of frequencies must be an odd number.

[0046] (2) In the complete implementation process of the present invention, the number of calculations required for a single cycle is relatively small. Compared with the signal construction method that requires solving complex equations or performing large-scale optimization, the total time consumption of the present invention is likely to be shorter.

[0047] In order to make the above-mentioned objects, features and advantages of the present invention more obvious and easy to understand, preferred embodiments are given below and described in detail with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0048] The accompanying drawings in the specification, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0049] Figure 1 A flow chart of a method for constructing a pseudo-random signal based on controllable phase modulation superposition of arbitrary square waves in one or more embodiments of the present invention;

[0050] Figure 2 The initial waveform and its spectrum diagram constructed in one or more embodiments of the present invention, wherein (a) is a waveform diagram and (b) is a spectrum diagram;

[0051] Figure 3 The waveform and its spectrum diagram obtained by superimposing a second square wave on the initial waveform in one or more embodiments of the present invention, wherein (a) is a waveform diagram and (b) is a spectrum diagram;

[0052] Figure 4 For Figure 3The waveform and its spectrum obtained by continuing to superimpose all the remaining frequencies in the frequency group on the basis of the waveform shown, where (a) is the waveform diagram and (b) is the spectrum diagram;

[0053] Figure 5 For Figure 4 The waveform and its spectrum obtained after hard limiting of the waveform shown, where (a) is the waveform diagram and (b) is the spectrum diagram;

[0054] Figure 6 For Figure 5 The waveform shown is a waveform that meets the assumed quality requirements after repeated cyclic superposition processing and its spectrum diagram, where (a) is the waveform diagram and (b) is the spectrum diagram. DETAILED DESCRIPTION

[0055] The present invention will be further described below in conjunction with the accompanying drawings and embodiments.

[0056] It should be noted that the following detailed descriptions are all illustrative and intended to provide further explanation of the present invention. Unless otherwise specified, all technical and scientific terms used herein have the same meanings as those commonly understood by those skilled in the art to which the present invention belongs.

[0057] It should be noted that the terms used herein are only for describing specific embodiments and are not intended to limit exemplary embodiments according to the present invention. As used herein, unless the context clearly indicates otherwise, the singular form is also intended to include the plural form. In addition, it should be understood that when the terms "comprising" and / or "including" are used in this specification, it indicates the presence of features, steps, operations, devices, components and / or combinations thereof.

[0058] In the absence of conflict, the embodiments in this application and the features in the embodiments may be combined with each other.

[0059] Embodiment 1

[0060] A method for constructing a pseudo-random signal based on arbitrary square wave controllable phase modulation superposition includes the following steps:

[0061] ① The frequency group that needs to be included in the given waveform (hereinafter referred to as "frequency group", the frequency in the frequency group can be given arbitrarily, but the total number must be an odd number, and the frequencies are arranged from low to high), the given waveform sampling rate (It must be at least twice the highest frequency in the frequency group. It is recommended to choose 360 ​​times the highest frequency in the frequency group). Calculate the length (Sampling rate Divide by the greatest common factor of all frequencies in the frequency group);

[0062] ② Construct a square wave with the lowest frequency in the frequency group at an arbitrary phase as the initial waveform, which is called waveform ;

[0063] ③ Waveform Perform fast Fourier transform to obtain its spectrum and then get the specified frequency The coefficient of ( is a complex number), where "specified frequency" refers to the next frequency in the frequency group sorted from low to high, and the amplitude of the lowest frequency in the frequency group;

[0064] ④According to the coefficient Angle , the ratio of the amplitude to the lowest frequency amplitude in the frequency group , get the phase of the next square wave :

[0065] ;

[0066] ⑤ Phase The base frequency is constructed as Square wave ;

[0067] ⑥ Square wave With waveform Add together to get a new waveform ;

[0068] ⑦ Repeat the above steps ③ to ⑥ until all frequencies in the frequency group are added to the waveform middle;

[0069] ⑧ Waveform Perform hard limiting to make the waveform The part greater than 1 is set to 1, and the part less than -1 is set to -1;

[0070] ⑨For waveform Perform fast Fourier transform, obtain its spectrum, take out the coefficients of all frequencies in the corresponding frequency group and calculate the amplitude to get an array , is the total number of frequencies in the frequency group, and its coefficient of variation is calculated by the following formula and theoretical energy ratio :

[0071] ;

[0072] ;

[0073] in For array The coefficient of variation reflects the unevenness of these amplitudes, and the theoretical energy ratio reflects the ratio of the average value of these amplitudes to the theoretical maximum amplitude.

[0074] ⑩ If the coefficient of variation and theoretical energy ratio If the quality requirements are met, the waveform is output Otherwise, starting from the second frequency in the frequency group, repeat steps ③ to ⑨ until the coefficient of variation and theoretical energy ratio Meet quality requirements.

[0075] The parameters required to construct a square wave are: length , sampling rate , Phase , fundamental frequency , the square wave is constructed as follows:

[0076] ①Create a Array of rows and 2 columns , used to save the position index of the square wave with a value of -1, Determined by the following formula:

[0077] ;

[0078] ②Yes No. The first row is filled with , the second column is filled with elements , Stands for rounding to the nearest integer.

[0079] ③Will All elements plus .

[0080] ④Create a one-dimensional array with all values ​​1 , whose length is ( The maximum value of The larger one between.

[0081] ⑤From Delete the beginning of the length that exceeds part, and the required square wave is obtained.

[0082] Embodiment 2

[0083] This embodiment discloses a pseudo-random signal construction method based on arbitrary square wave controllable phase modulation superposition, by phase adjustment and superposition of square waves containing different frequency components, and hard limiting of the superposed waveform, to obtain a pseudo-random signal waveform containing specified frequency components, and the amplitudes of these specified frequency components are basically consistent. Figure 1 As shown, the method comprises the following steps:

[0084] Step 1: Given the parameters of the waveform to be generated, including the frequency group, waveform length, and sampling rate. The frequency group can be given arbitrarily according to actual needs. Here, the frequency group shown in Table 1 is taken as an example (for the convenience of observation, the table is represented by three groups, which are actually one-dimensional arrays arranged from small to large). The sampling rate is selected to be 360 ​​times the highest frequency ( ), the waveform length is (1 is the greatest common factor of all frequencies in the frequency group).

[0085] Table 1 Example frequency groups (unit: Hz)

[0086]

[0087] Step 2: Construct a square wave with a base frequency of 1 Hz (the lowest frequency in the frequency group) with 0 phase as the initial phase and use it as the waveform , its waveform and spectrum are as follows Figure 2 shown.

[0088] Step 3: From the frequency group, we know that the next frequency is 2Hz. The coefficient of this frequency is obtained from the spectrum. , calculate its angle , the ratio of the amplitude of the lowest frequency in the frequency group , according to the formula, the phase of the next square wave is obtained .

[0089] Step 4: Build Phase , a square wave with a base frequency of 2Hz, and compare it with the original waveform Superposition to obtain a new waveform , its waveform and spectrum are as follows Figure 3 As shown, it can be seen that the amplitudes of 1Hz and 2Hz in the current spectrum are equal.

[0090] Step 5: Repeat steps 3 and 4 until all frequencies in the frequency group are added to the waveform. The waveform and spectrum obtained at the end of this step are as follows: Figure 4 As shown, it can be seen that the amplitudes of all frequency groups in the current spectrum are equal.

[0091] Step 6: Hard limit the waveform obtained in step 5, with the upper limit of the absolute value of the amplitude being 1, that is, the part of the waveform greater than 1 is set to 1, and the part less than -1 is set to -1. The waveform and spectrum obtained at the end of this step are as follows: Figure 5 As shown, it can be seen that the hard limiting operation converts the waveform into a binary waveform and affects the energy distribution in the spectrum.

[0092] Step 7: Calculate the coefficient of variation according to the formula and theoretical energy ratio If the quality requirement is , , then the generated waveform meets the quality requirements and the generated waveform is output; if the quality requirement is , , then the waveform generated at this time does not meet the quality requirements, proceed to step 8.

[0093] Step 8: For the waveform obtained in step 7, repeat steps 3 to 7 starting from the second smallest frequency in the frequency group until the result meets the quality requirements. , In the case of, after repeating the above steps for 1 round, the result that meets the requirements is obtained ( , ), the results are as follows Figure 6 shown.

[0094] Thus, the construction of a binary signal waveform that meets the specified quality requirements is achieved given a frequency group.

[0095] Embodiment 3

[0096] A pseudo-random signal construction system based on arbitrary square wave controllable phase modulation superposition, comprising:

[0097] A parameter configuration module, configured to set the parameters of the waveform to be generated, the parameters including frequency group, waveform length, and sampling rate;

[0098] An initial construction module is configured to construct a square wave with a fundamental frequency as the lowest frequency in the frequency group with an initial phase as an initial waveform;

[0099] A transform module is configured to perform a fast Fourier transform on the initial waveform to obtain its spectrum and then obtain coefficients of a specified frequency;

[0100] A parameter update calculation module is configured to construct a phase of the next square wave according to the angle and amplitude of the coefficient and the ratio of the lowest frequency amplitude in the frequency group, and to construct a square wave with a base frequency of the specified frequency with the phase;

[0101] A waveform superposition module is configured to add the constructed square wave to the current waveform to obtain a new waveform;

[0102] An iteration module is configured to repeatedly call the transformation module, the parameter update calculation module and the waveform superposition module until all frequencies in the frequency group are added to the new waveform;

[0103] The evaluation index calculation module is configured to perform hard limiting on the final new waveform, then perform fast Fourier transform, obtain its spectrum, extract the coefficients of all frequencies of the corresponding frequency group and calculate the amplitude to obtain an array, and calculate the coefficient of variation and theoretical energy proportion according to the array;

[0104] The judgment module is configured to output the final new waveform if the coefficient of variation and the theoretical energy ratio meet the requirements, otherwise, starting from the second frequency in the frequency group, repeatedly call the transformation module to the evaluation index calculation module until the coefficient of variation and the theoretical energy ratio meet the requirements.

[0105] It should be understood by those skilled in the art that embodiments of the present invention may be provided as methods, systems, or computer program products. Therefore, the present invention may take the form of an all-hardware embodiment, an all-software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of one or more computer-usable storage media (including but not limited to disk storage, CD - ROM , optical storage, etc.).

[0106] The present invention is described with reference to flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to embodiments of the present invention. It should be understood that each process and / or block in the flowchart and / or block diagram, as well as the combination of processes and / or blocks in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to produce a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0107] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0108] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 A step that specifies a function in one or more boxes.

[0109] The above description is only a preferred embodiment of the present invention and is not intended to limit the present invention. For those skilled in the art, the present invention may have various modifications and variations. Any modification, equivalent replacement, improvement, etc. made by those skilled in the art within the spirit and principle of the present invention without creative labor shall be included in the protection scope of the present invention.

Claims

1. A method for constructing a pseudo-random signal based on arbitrary square wave controllable phase modulation superposition, characterized in that: The following steps are involved: (1) Setting the parameters of the waveform to be generated, including the frequency group, waveform length, and sampling rate; (2) Construct a square wave with the lowest frequency in the frequency group using the initial phase as the initial waveform; (3) Perform a fast Fourier transform on the initial waveform to obtain its spectrum and then obtain the coefficient of the specified frequency; (4) constructing the phase of the next square wave according to the angle and amplitude of the coefficient and the ratio of the lowest frequency amplitude in the frequency group, and constructing a square wave with a base frequency of the specified frequency with the phase; (5) Add the constructed square wave to the current waveform to obtain a new waveform; (6) Repeat steps (3) to (5) until all frequencies in the frequency group are added to the new waveform; (7) After hard limiting the final new waveform, a fast Fourier transform is performed to obtain its spectrum, and the coefficients of all frequencies of the corresponding frequency group are taken out and the amplitudes are calculated to obtain an array. Based on the array, the coefficient of variation and the theoretical energy ratio are calculated; (8) If the coefficient of variation and the theoretical energy ratio meet the requirements, the final new waveform is output; otherwise, starting from the second frequency in the frequency group, steps (3) to (7) are repeated until the coefficient of variation and the theoretical energy ratio meet the requirements; In the step (1), the frequencies in the frequency group are set arbitrarily, the total number is an odd number, and the frequencies are arranged from low to high; The sampling rate is at least twice the highest frequency in the frequency group, and the waveform length is the ratio of the sampling rate to the greatest common factor of all frequencies in the frequency group; In step (4), the process of constructing the phase of the next square wave according to the ratio of the angle and amplitude of the coefficient to the lowest frequency amplitude in the frequency group includes: ; in, is the angle of the coefficient, is the ratio of the amplitude of the coefficient to the amplitude of the lowest frequency in the frequency group, is the phase of the next square wave; In step (5), the process of constructing a square wave includes: creating Array of rows and 2 columns , used to save the position index of the square wave with a value of -1, Determined by the following formula: ; right No. The first row is filled with , the second column is filled with elements , Represents rounding to the nearest integer; Will All elements plus ; Create a one-dimensional array with all values ​​1 , whose length is The maximum value in The greater of right Each row of , starting with the first column element, ending with the second column element, and using 1 as the interval to construct an index sequence. The value of the corresponding index position is set to -1; From a one-dimensional array Delete the beginning of the length that exceeds The required square wave is obtained; in, is the waveform length, is the sampling rate, is the calculated phase, is the fundamental frequency, is the total number of frequencies in the frequency group.

2. A method for constructing a pseudo-random signal based on arbitrary square wave controllable phase modulation superposition as claimed in claim 1, characterized in that: In the step (2), the initial phase is an arbitrary phase.

3. A method for constructing a pseudo-random signal based on arbitrary square wave controllable phase modulation superposition as claimed in claim 1, characterized in that: In step (7), the process of hard limiting the final new waveform includes: setting the part of the final new waveform that is greater than 1 to 1, and setting the part that is less than -1 to -1.

4. A method for constructing a pseudo-random signal based on arbitrary square wave controllable phase modulation superposition as claimed in claim 1, characterized in that: In the step (7), the waveform Perform fast Fourier transform, obtain its spectrum, take out the coefficients of all frequencies in the corresponding frequency group and calculate the amplitude to get an array , is the total number of frequencies in the frequency group, and its coefficient of variation is calculated by the following formula and theoretical energy ratio : ； ； in For array The coefficient of variation reflects the unevenness of these amplitudes, and the theoretical energy ratio reflects the ratio of the average value of these amplitudes to the theoretical maximum amplitude.

5. A pseudo-random signal construction system based on arbitrary square wave controllable phase modulation superposition, characterized in that: include: A parameter configuration module is configured to set the parameters of the waveform to be generated, the parameters including the frequency group, the waveform length, and the sampling rate; the frequencies in the frequency group are arbitrarily set, the total number is an odd number, and the frequencies are arranged from low to high; the sampling rate is at least twice the highest frequency in the frequency group, and the waveform length is the ratio of the sampling rate to the greatest common factor of all the frequencies in the frequency group; An initial construction module is configured to construct a square wave with a fundamental frequency as the lowest frequency in the frequency group with an initial phase as an initial waveform; A transform module is configured to perform a fast Fourier transform on the initial waveform to obtain its spectrum and then obtain coefficients of a specified frequency; A parameter update calculation module is configured to construct a phase of the next square wave according to the angle and amplitude of the coefficient and the ratio of the lowest frequency amplitude in the frequency group, and to construct a square wave with a base frequency of the specified frequency with the phase; The process of constructing the phase of the next square wave according to the angle and amplitude of the coefficient and the ratio of the amplitude of the lowest frequency in the frequency group includes: ;in, is the angle of the coefficient, is the ratio of the amplitude of the coefficient to the amplitude of the lowest frequency in the frequency group, is the phase of the next square wave; The waveform superposition module is configured to add the constructed square wave to the current waveform to obtain a new waveform; the process of constructing the square wave includes: creating Array of rows and 2 columns , used to save the position index of the square wave with a value of -1, Determined by the following formula: ; right No. The first row is filled with , the second column is filled with elements , Represents rounding to the nearest integer; Will All elements plus ; Create a one-dimensional array with all values ​​1 , whose length is The maximum value in The greater of right Each row of , starting with the first column element, ending with the second column element, and using 1 as the interval to construct an index sequence. The value of the corresponding index position is set to -1; From a one-dimensional array Delete the beginning of the length that exceeds The required square wave is obtained; in, is the waveform length, is the sampling rate, is the calculated phase, is the fundamental frequency, is the total number of frequencies in the frequency group; An iteration module is configured to repeatedly call the transformation module, the parameter update calculation module and the waveform superposition module until all frequencies in the frequency group are added to the new waveform; The evaluation index calculation module is configured to perform hard limiting on the final new waveform, then perform fast Fourier transform, obtain its spectrum, extract the coefficients of all frequencies of the corresponding frequency group and calculate the amplitude to obtain an array, and calculate the coefficient of variation and theoretical energy proportion according to the array; The judgment module is configured to output the final new waveform if the coefficient of variation and the theoretical energy ratio meet the requirements, otherwise, starting from the second frequency in the frequency group, repeatedly call the transformation module to the evaluation index calculation module until the coefficient of variation and the theoretical energy ratio meet the requirements.

6. A computer-readable storage medium, characterized in that: Used to store computer instructions, which, when executed by a processor, complete the steps of the method according to any one of claims 1 to 4.

7. An electronic device, characterized in that: The method comprises a memory and a processor, and computer instructions stored in the memory and executed on the processor, wherein when the computer instructions are executed by the processor, the steps in the method according to any one of claims 1 to 4 are completed.

Images