Phase estimation method, signal processing channel and storage medium

Abstract

A phase estimation method, a signal processing channel and a storage medium are applied to the technical field of signal processing, and the phase estimation method comprises the following steps: obtaining an original input signal and a to-be-measured phase table of a coarse interval; performing a traversal test based on the to-be-tested phase table of the coarse interval to obtain M1 amplitude estimation values corresponding to the M1 phase values; obtaining a target phase corresponding to an amplitude estimation target value in the M1 amplitude estimation values, and generating a fine-interval to-be-measured phase table according to the target phase; and performing traversal test based on the fine-interval phase table to be tested to obtain M2 amplitude estimation values corresponding to the M2 phase values, and obtaining the phase corresponding to the minimum amplitude estimation value based on the M2 amplitude estimation values and the corresponding phase values. The phase estimation is performed by analyzing the amplitude estimation value of the signal to be tested, so that the phase estimation process is simplified, the complexity of phase estimation is reduced based on the combination of the coarse-interval phase meter to be tested and the fine-interval phase meter to be tested, and the estimated phase can be obtained more accurately.

Description

Technical Field

[0001] The present application relates to the field of signal processing technology, and in particular to a phase estimation method, a signal processing channel and a storage medium. Background Art

[0002] As an important device for converting digital signals into analog signals, the Digital-to-Analog Converter (DAC) is widely used in all walks of life. However, due to the influence of nonlinear factors such as mismatch of internal components of the DAC, the output of the DAC always contains other signals besides the main signal of the input signal, which in turn affects the spurious free dynamic range (SFDR) of the DAC. Among them, the harmonics and spurs corresponding to the main signal in the output signal of the DAC have a greater impact on SFDR.

[0003] In the current technical solution, other analytical instruments are mainly used to analyze the analog signal output by DAC to directly analyze the harmonic phase of harmonics or the spurious phase of spurious signals. However, this requires high accuracy of the analytical instrument, and the cost of high-precision analytical instruments is also relatively high. In addition, the direct analysis process is complicated and not suitable for large-scale applications. Therefore, a new technical solution is needed to obtain the phase of harmonics or spurious signals generated by DAC. Summary of the invention

[0004] The main technical problem solved by the present application is how to simplify the acquisition of the phase of harmonics or spurious generated by DAC.

[0005] According to the first aspect, an embodiment provides a phase estimation method, including:

[0006] Obtaining an original input signal and a coarse-interval phase table to be measured; wherein the coarse-interval phase table to be measured includes M1 phase values, and M1 is a positive integer greater than 1;

[0007] Perform traversal test on the phase table to be tested based on the coarse interval: generate a cancellation signal based on the harmonic frequency and harmonic amplitude value of the harmonics of the original input signal or based on the spurious frequency and spurious amplitude value of the spurious of the original input signal, and any phase value in the corresponding phase table to be tested; superimpose the cancellation signal and the original input signal to obtain a digital signal to be tested; input the digital signal to be tested into a DAC for digital-to-analog conversion to obtain an analog signal to be tested, perform spectrum analysis on the analog signal to be tested to obtain an amplitude estimation value of the harmonic frequency or spurious frequency corresponding to any phase value; traverse the cancellation signal corresponding to each phase value in the corresponding phase table to be tested to obtain each amplitude estimation value corresponding to each phase value;

[0008] Acquire M1 amplitude estimation values ​​corresponding to the M1 phase values, and determine a target phase corresponding to an amplitude estimation target value among the M1 amplitude estimation values;

[0009] Generate a fine-interval phase table to be measured according to the target phase; wherein the fine-interval phase table to be measured includes M2 phase values, and M2 is a positive integer greater than 1;

[0010] The traversal test is performed based on the fine-interval phase table to be tested to obtain M2 amplitude estimation values ​​corresponding to the M2 phase values, and the phase corresponding to the minimum amplitude estimation value is obtained based on the M2 amplitude estimation values ​​and the corresponding phase values.

[0011] In some embodiments, obtaining a target phase corresponding to an amplitude estimation target value in the M1 amplitude estimation values ​​includes:

[0012] Obtaining an average value of the M1 amplitude estimation values;

[0013] Determining whether the average value is greater than the harmonic amplitude value or the spurious amplitude value of the original input signal;

[0014] If it is greater than, then the phase corresponding to the maximum value among the M1 amplitude estimation values ​​is obtained and used as the target phase corresponding to the amplitude estimation target value;

[0015] If it is less than, the phase corresponding to the minimum value among the M1 amplitude estimation values ​​is obtained and used as the target phase corresponding to the amplitude estimation target value.

[0016] In some embodiments, obtaining the phase corresponding to the minimum amplitude estimation value based on the M2 amplitude estimation values ​​and the corresponding phase values ​​includes:

[0017] If the phase corresponding to the minimum value among the M1 amplitude estimation values ​​is obtained, then the phase corresponding to the minimum value among the M2 amplitude estimation values ​​is determined to be the phase corresponding to the minimum amplitude estimation value;

[0018] If the phase corresponding to the maximum value among the M1 amplitude estimation values ​​is obtained, the phase corresponding to the maximum value among the M2 amplitude estimation values ​​is obtained, and the phase corresponding to the minimum amplitude estimation value is determined based on the phase corresponding to the maximum value.

[0019] In some embodiments, determining the phase corresponding to the minimum amplitude estimate value according to the phase corresponding to the maximum value includes:

[0020] Shifting the phase corresponding to the maximum value among the M2 amplitude estimation values ​​by a preset phase so that the shifted phase corresponds to the amplitude estimation minimum value;

[0021] The shifted phase is output, where the shifted phase is the phase corresponding to the minimum value of the amplitude estimation.

[0022] In some embodiments, if the phase value corresponding to the maximum value among the M2 amplitude estimation values ​​is between [0, π), the preset phase of the translation is π;

[0023] If the phase value corresponding to the maximum value among the M2 amplitude estimation values ​​is between [π, 2π), the preset phase of the translation is -π.

[0024] In some embodiments, the coarse-interval phase table to be measured is generated by:

[0025] Dividing M1 phase values ​​from the phase range [a, a+π] according to a preset first phase interval to generate the coarse interval phase table to be measured; wherein a is a phase greater than or equal to 0 and less than or equal to π;

[0026] or,

[0027] An estimated phase empirical value is obtained, and a phase table to be measured of the coarse interval is generated according to the estimated phase empirical value.

[0028] In some embodiments, obtaining the estimated phase empirical value includes:

[0029] Obtaining a phase estimation value corresponding to a harmonic frequency or a spurious frequency of a known input signal after the known input signal is input into the DAC;

[0030] Based on the frequency relationship between the known input signal and the original input signal, the estimated phase empirical value corresponding to the harmonic of the original input signal is obtained from the phase estimation value; or,

[0031] Based on the frequency relationship between the known input signal and the original input signal, and the frequency relationship between the spurious frequency of the known input signal and the spurious frequency of the original input signal, the estimated phase empirical value corresponding to the spur of the original input signal is obtained from the phase estimation value.

[0032] In some embodiments, generating the coarse-interval phase table to be measured according to the estimated phase empirical value includes:

[0033] According to the preset second phase interval, M2 phase values ​​are divided from the nearby range [bc, b+c] of the estimated phase empirical value to generate the coarse interval phase table to be tested; wherein b is the estimated phase empirical value and c is the preset phase value.

[0034] In some embodiments, the phase interval between the M1 phase values ​​is greater than the phase interval between the M2 phase values; and / or, the phase range included in the M1 phase values ​​is greater than the phase range included in the M2 phase values, and the value of M1 is smaller than the value of M2.

[0035] In some embodiments, the performing spectrum analysis on the simulated signal to be tested includes:

[0036] Inputting the simulated signal to be tested into a spectrum analyzer for spectrum analysis;

[0037] or,

[0038] The simulated signal to be tested is input into an analog-to-digital converter for analog-to-digital conversion and then subjected to Fourier transform processing for spectrum analysis.

[0039] In some embodiments, for harmonics, the corresponding frequency of the cancellation signal is the harmonic frequency of the harmonics of the original input signal; for spurious signals, the corresponding frequency of the cancellation signal is the spurious frequency of the spurious signals of the original input signal.

[0040] According to a second aspect, an embodiment provides a signal processing channel for implementing the phase estimation method as described in the first aspect, wherein the signal processing channel comprises at least two DDS phase accumulators, and the signal processing circuits where each of the DDS phase accumulators is located are connected in parallel, and the number of parallel connections is defined as an integer H greater than 1;

[0041] One of the signal processing circuits including a DDS phase accumulator generates the original input signal after phase-to-amplitude conversion processing and attenuation processing, and the remaining signal processing circuits including DDS phase accumulators generate cancellation signals corresponding to harmonics or cancellation signals corresponding to spurious signals based on the corresponding phase values ​​in the phase table to be measured through phase-to-amplitude conversion processing.

[0042] In some embodiments, the signal processing channel is used to process multiple harmonics, or multiple spurious signals, or at least one harmonic and at least one spurious signal. If the value of H is ≥3, the remaining other signal processing circuits including the DDS phase accumulator are used to generate cancellation signals corresponding to multiple harmonics, or cancellation signals corresponding to multiple spurious signals, or cancellation signals corresponding to at least one harmonic and at least one spurious signal respectively.

[0043] The maximum number of times the remaining signal processing circuits including the DDS phase accumulator are called is equal to twice the sum of the number of harmonics and the number of spurious signals.

[0044] In some embodiments, the signal processing channel also includes a phase counter and a memory, wherein the phase counter is used to accumulate and count the M1 phase values ​​or M2 phase values ​​corresponding to the cancellation signal during traversal, and the memory is used to store the corresponding phase values ​​in the phase table to be tested and the corresponding amplitude estimation values.

[0045] In some embodiments, obtaining a total number L of harmonics and spurious signals to be processed;

[0046] Divide twice the total number L by H-1 and round up to get the number of times the signal processing channel is called;

[0047] According to the number of times the signal processing channel is called, the phase counter and the memory are reset to clear the amplitude estimation value and the corresponding phase value stored in history.

[0048] According to the third aspect, an embodiment provides a computer-readable storage medium, on which a program is stored, and the program can be executed by a processor to implement the method as described in the first aspect.

[0049] According to the phase estimation method, signal processing channel and storage medium of the above-mentioned embodiment, after the original input signal and the cancellation signal are superimposed, if the phase of the harmonic frequency or spurious frequency corresponding to the original input signal is in phase with the phase of the harmonic frequency or spurious frequency corresponding to the cancellation signal, the amplitude estimation value of the harmonic frequency or spurious frequency corresponding to the phase reaches the maximum value, and vice versa, if the phase is inverted, the amplitude estimation value of the harmonic frequency or spurious frequency corresponding to the phase reaches the minimum value. Among them, the coarse-interval phase table to be tested is first traversed to obtain the target phase of the approximate phase position, and then the fine-interval phase table to be tested is generated based on the target phase, and then the fine-interval phase table to be tested is traversed again to obtain the phase corresponding to the minimum amplitude estimation value. Since the phase corresponding to the minimum amplitude estimation value can be obtained mainly by superimposing the cancellation signal and the original input signal, and analyzing the amplitude estimation value corresponding to the harmonic frequency or spurious frequency in the simulated signal to be tested, the phase estimation process is simplified. At the same time, a coarse estimation of the phase is performed based on the coarse interval phase table to be measured, and a fine estimation of the phase is performed based on the fine interval phase table to be measured. The combination of the two reduces the complexity of the phase estimation and can obtain the estimated phase more accurately. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] Figure 1 A schematic diagram of the structure of a signal processing channel according to an embodiment;

[0051] Figure 2 A schematic diagram of the structure of a signal processing channel in another embodiment;

[0052] Figure 3A phase estimation flow chart of a signal processing channel of an embodiment;

[0053] Figure 4 A recording diagram of harmonic amplitudes of an embodiment;

[0054] Figure 5 A recording diagram of harmonic amplitudes of another embodiment;

[0055] Figure 6 A flow chart of a phase estimation method according to an embodiment. DETAILED DESCRIPTION

[0056] The present application is further described in detail below by specific embodiments in conjunction with the accompanying drawings. Wherein similar elements in different embodiments use associated similar element numbers. In the following embodiments, many detailed descriptions are intended to enable the present application to be better understood. However, those skilled in the art can easily recognize that some of the features can be omitted in different situations, or can be replaced by other elements, materials, and methods. In some cases, some operations related to the present application are not shown or described in the specification, in order to avoid the core part of the present application being overwhelmed by too much description, and for those skilled in the art, it is not necessary to describe these related operations in detail, and the related operations can be fully understood according to the description in the specification and the general technical knowledge in the art.

[0057] In addition, the features, operations or characteristics described in the specification can be combined in any appropriate manner to form various implementations. At the same time, the steps or actions in the method description can also be interchanged or adjusted in a manner that is obvious to those skilled in the art. Therefore, the various sequences in the specification and the drawings are only for the purpose of clearly describing a certain embodiment and are not meant to be a required sequence, unless otherwise specified that a certain sequence must be followed.

[0058] The serial numbers of the components in this document, such as "first", "second", etc., are only used to distinguish the objects described and do not have any order or technical meaning. The "connection" and "coupling" mentioned in this application, unless otherwise specified, include direct and indirect connections (couplings).

[0059] In some embodiments of the present application, when two signals of the same frequency are superimposed, if the two signals are in phase, the superimposed amplitude will increase, and if the two signals are in anti-phase, the superimposed amplitude will cancel out. The original input signal and the cancellation signal are obtained, and then the cancellation signal and the original input signal are superimposed, and the DAC performs digital-to-analog conversion, and finally the analog signal output by the digital-to-analog converter is analyzed to obtain the amplitude estimation value corresponding to the harmonic frequency or the spurious frequency. If the phase of the harmonic frequency or the spurious frequency corresponding to the original input signal is in phase with the phase of the harmonic frequency or the spurious frequency corresponding to the cancellation signal, the amplitude estimation value of the harmonic frequency or the spurious frequency corresponding to the phase reaches the maximum value, otherwise if the phase is in anti-phase, the amplitude estimation value of the harmonic frequency or the spurious frequency corresponding to the phase reaches the minimum value. Among them, the coarse-interval phase table to be measured is first traversed to obtain the target phase of the approximate phase position, and then the fine-interval phase table to be measured is generated based on the target phase, and the fine-interval phase table to be measured is traversed again to obtain the phase corresponding to the minimum amplitude estimation value. Since the phase corresponding to the minimum amplitude estimation value can be obtained by superimposing the cancellation signal and the original input signal and analyzing the amplitude estimation value corresponding to the harmonic frequency or spurious frequency in the simulated test signal, the phase estimation process is simplified. At the same time, the coarse estimation of the phase based on the coarse interval test phase table and the fine estimation of the phase based on the fine interval test phase table reduce the complexity of the phase estimation and can obtain the estimated phase more accurately.

[0060] Some embodiments provide a signal processing channel, please refer to Figure 1 The signal processing channel includes at least two DDS phase accumulators 10, and the signal processing circuits where each DDS phase accumulator 10 is located are connected in parallel, and the number of parallel connections is defined as an integer H greater than 1.

[0061] The DDS phase accumulator 10 is used to complete the accumulation of phase values ​​according to the received frequency control word to obtain a phase accumulation value. Figure 2 In some embodiments, the components within the signal processing circuit where each DDS phase accumulator 10 is located are connected in series. For example, the DDS phase accumulator 10 includes a phase adder and a phase register connected in series. The phase register is used to output the current accumulated phase to the phase adder. The phase adder adds the frequency control word and the accumulated phase output by the phase register, and sends the result of the addition to the phase register for the next phase accumulation. In some embodiments, the frequency control word can be generated based on a processor, wherein the processor can also be used to adjust the initial phase of the phase register. In some embodiments, the processor can be implemented based on a device with data processing capabilities such as a CPU, an FPGA, and a single-chip microcomputer. Among them, the frequency control word is determined by the frequency of the input signal, the clock frequency, and the total number of bits of the DDS phase accumulator 10.

[0062] In some embodiments, the signal processing circuit where each DDS phase accumulator 10 is located can be configured with a phase amplitude mapping module 20, and the phase amplitude mapping module 20 is used to obtain the corresponding signal amplitude according to the phase accumulation value output by the corresponding DDS phase accumulator 10, and output the corresponding digital waveform signal. In some embodiments, the phase amplitude mapping module 20 can obtain the signal amplitude corresponding to the phase accumulation value based on a phase lookup table, for example, the phase lookup table can be a sine lookup table. In some embodiments, the phase amplitude mapping module 20 can also calculate the signal amplitude corresponding to the phase accumulation value in real time based on algorithms such as coordinate rotation.

[0063] In some embodiments, the signal processing circuits where the DDS phase accumulators 10 are located may also share a phase-amplitude mapping module 20 , for example, the phase-amplitude mapping module 20 is multiplexed between the DDS phase accumulators 10 .

[0064] In some embodiments, the digital waveform signal output by the phase amplitude mapping module 20 can be converted from digital to analog by the digital-to-analog conversion unit 30 to output an analog waveform signal. In some embodiments, the digital-to-analog conversion unit 30 can realize the digital-to-analog conversion of the signal based on an analog-to-digital converter. In some embodiments, the digital-to-analog conversion unit 30 can also include a signal adder, which is used to superimpose multiple digital waveform signals and then output them to the analog-to-digital converter to realize the digital-to-analog conversion of the superimposed signal.

[0065] In the above embodiment, each signal processing circuit including a DDS phase accumulator 10 is used to generate a digital waveform signal. Therefore, each time the signal processing circuit is adjusted, a digital waveform signal can be generated to correspond to a harmonic or spurious cancellation signal, and each time the signal processing channel is adjusted, a plurality of digital waveform signals corresponding to the number of DDS phase accumulators 10 can be generated to correspond to multiple harmonics, or multiple spurious signals, or spurious and harmonic cancellation signals.

[0066] The above is some description of the signal processing channel, please refer to Figure 3 ,The following describes the process of phase estimation based on the signal processing channel in detail.

[0067] The original input signal is obtained, wherein one of the signal processing channels includes a signal processing circuit of a DDS phase accumulator 10, and generates the original input signal after phase-amplitude conversion processing and attenuation processing.

[0068] In some embodiments, the signal processing circuit may further include an attenuator, so that attenuation processing of the signal is performed based on the attenuator.

[0069] In some embodiments, the harmonics of the original input signal may include one harmonic or multiple harmonics, for example, the second harmonic or the third harmonic, or multiple harmonics such as the second harmonic and the third harmonic, so the cancellation signal may be a signal corresponding to the second harmonic or the third harmonic, or a signal corresponding to multiple harmonics such as the second harmonic and the third harmonic. In some embodiments, the spurious of the original input signal may include one spurious or multiple spurious, so the cancellation signal may be a signal corresponding to one spurious, or a signal corresponding to multiple spurious. In this embodiment, the harmonics and spurious of the original input signal refer to the harmonics and spurious generated after the original input signal is input into the DAC.

[0070] The coarse-interval phase table to be measured is obtained, wherein the coarse-interval phase table to be measured includes M1 phase values, where M1 is a positive integer greater than 1.

[0071] In some embodiments, M1 phase values ​​uniformly distributed in the phase range [a, a+π] are divided according to the first phase interval to generate a coarse interval phase table to be measured. Wherein, a is a phase greater than or equal to 0 and less than or equal to π, and in some embodiments, a can take a value of 0. In this embodiment, after the original input signal is superimposed with the cancellation signal, there is a specific relationship between the amplitude estimation value corresponding to the harmonic or spurious and the phase value in the phase table to be measured, for example, when the phase value is the same as the phase of the corresponding harmonic or spurious, the amplitude estimation value reaches the maximum value, and conversely, when the phase is reversed, the amplitude estimation value reaches the minimum value, and the change of the amplitude estimation value is symmetrical based on the change of the phase value in the phase table to be measured, that is, the change of the amplitude estimation value corresponding to the phase value in any π range is symmetrical with the change of the amplitude estimation value corresponding to the phase value in another π range. Therefore, based on any π range, a coarse interval phase table to be measured is generated, and phase estimation can be performed based on symmetry, thereby reducing the amount of calculation to simplify the complexity of the process. In some embodiments, the selected phase may be formed into a corresponding phase table to be measured and stored, so that subsequent adjustments may be made directly based on the phase in the phase table to be measured.

[0072] In some embodiments, an estimated phase empirical value is obtained, and a coarse-interval phase table to be measured is generated based on the estimated phase empirical value. In some embodiments, M1 phase values ​​with uniform distribution can be selected based on the front and rear ranges of the estimated phase empirical value to generate a coarse-interval phase table to be measured, for example, the estimated phase empirical value is taken as the center position, and M1 phase values ​​with uniform distribution are selected in a preset range greater than the estimated phase empirical value and a preset range less than the estimated phase empirical value. In some embodiments, for each original input signal without an estimated phase empirical value, a coarse-interval phase table to be measured generated based on the phase range [a, a+π] can be shared, and for each original input signal with an estimated phase empirical value, a coarse-interval phase table to be measured generated based on the estimated phase empirical value is used.

[0073] In some embodiments, after the known input signal is input into the DAC, the phase estimation value corresponding to the harmonic frequency or spurious frequency of the known input signal is obtained, wherein, since there is a certain relationship between the frequency of the signal and the phase of the signal, the phase corresponding to the harmonic frequency or spurious frequency of the original input signal can be estimated based on the frequency of the known input signal and the phase estimation value corresponding to the harmonic frequency or spurious frequency of the known input signal. For example, if the frequency between the known input signal and the original input signal is different, based on the frequency relationship between the known input signal and the original input signal, the estimated phase empirical value corresponding to the harmonic of the original input signal is obtained from the phase estimation value corresponding to the harmonic, or based on the frequency relationship between the known input signal and the original input signal, and the frequency relationship between the spurious frequency of the spurious of the known input signal and the spurious frequency of the spurious of the original input signal, the estimated phase empirical value corresponding to the spurious of the original input signal is obtained from the phase estimation value corresponding to the spurious. For example, if the frequency between the known input signal and the original input signal is the same, the phase estimation value corresponding to the harmonic is used as the estimated phase empirical value corresponding to the harmonic of the original input signal. For example, if the spurious frequency of the spurious input signal is known to be the same as the spurious frequency of the original input signal, the phase estimation value corresponding to the spurious frequency is used as the estimated phase empirical value corresponding to the spurious of the original input signal.

[0074] In the above embodiment, when signals of the same frequency are input to the same digital-to-analog converter, the phases corresponding to the generated harmonics are substantially consistent. Therefore, for a known input signal for which phase estimation has been completed, the estimated phase position obtained can be used for subsequent processing without the need for repeated acquisition. In addition, for an original input signal for which phase estimation has not been completed, the phase of the harmonic corresponding to the original input signal can be predicted based on the phase of the harmonic of the known input signal, and then the phase of the original input signal can be estimated based on the predicted phase, thereby further simplifying the acquisition complexity and improving the accuracy of the phase estimation.

[0075] In some embodiments, when a coarse-interval phase table to be measured is generated based on the estimated phase empirical value, M2 phase values ​​are divided from the nearby range [bc, b+c] of the estimated phase empirical value according to the preset second phase interval to generate a coarse-interval phase table to be measured, where b is the estimated phase empirical value and c is the preset phase value. In this embodiment, since the estimated phase empirical value has been obtained, a coarse-interval phase table to be measured can be generated within a smaller phase range. In some embodiments, the preset second phase interval can be the same as or different from the preset first phase interval.

[0076] The traversal test is performed based on the coarse-interval phase table to be tested. That is, the phase values ​​in the coarse-interval phase table to be tested are selected in turn, and the corresponding cancellation signal is generated based on the phase value, and is superimposed with the original input signal and input into the DAC, and then the spectrum analysis of the analog signal output by the DAC is performed to obtain the amplitude estimation value of the harmonic frequency or spurious frequency corresponding to the phase value.

[0077] In some embodiments, when obtaining the cancellation signal, the remaining signal processing circuits including the DDS phase accumulator 10 in the signal processing channel are subjected to phase-amplitude conversion processing based on the phase value in the corresponding phase table to be measured to generate a cancellation signal corresponding to the harmonic or a cancellation signal corresponding to the spurious.

[0078] In some embodiments, before obtaining the cancellation signal, the original input signal may be input into a digital-to-analog converter for digital-to-analog conversion to obtain an analog input signal corresponding to the output of the digital-to-analog converter, and then the analog input signal may be subjected to spectrum analysis to obtain an estimated amplitude value of the harmonics of the original input signal. Alternatively, the spurious frequency and spurious amplitude value of the original input signal may be obtained, for example, based on the spurious frequency and spurious amplitude value of the original input signal known by other measuring equipment.

[0079] In some embodiments, the analog input signal may be converted into a digital signal based on an analog-to-digital converter to obtain a corresponding digital signal, and then the digital signal may be transformed from the time domain to the frequency domain based on a processor or an external data processing device, such as performing a fast Fourier transform (FFT) process, and then analyzing the transformed frequency domain signal to obtain an amplitude estimate of the harmonic of the original input signal. Since the frequency of the FFT may have a deviation, the obtained amplitude estimate may be selected as the amplitude corresponding to the harmonic frequency of the harmonic, or may be selected as the maximum amplitude within a frequency range around the harmonic frequency.

[0080] Then, a cancellation signal is generated based on the harmonic frequency and harmonic amplitude value of the harmonics of the original input signal or based on the spurious frequency and spurious amplitude value of the spurs of the original input signal and any phase value in the corresponding phase table to be measured.

[0081] In some embodiments, when generating a cancellation signal, for harmonics, the frequency of the corresponding cancellation signal is the harmonic frequency of the harmonic of the original input signal, and for spurious, the frequency of the corresponding cancellation signal is the spurious frequency of the spurious of the original input signal.

[0082] In some embodiments, when the original input signal is obtained by attenuation, the amplitude of the digital signal to be measured after the original input signal is superimposed with the cancellation signal can be made to meet the maximum input range of the digital-to-analog converter or be less than the maximum input range. In some embodiments, when the amplitude of the digital signal to be measured meets the maximum input range of the digital-to-analog converter, the digital-to-analog converter can be fully utilized to improve the subsequent data processing accuracy. In some embodiments, the amplitude of the original input signal attenuation can be controlled by a processor or by an external data processing device.

[0083] In some embodiments, the remaining other signal processing circuits including the DDS phase accumulator 10 may generate a cancellation signal corresponding to a harmonic or a cancellation signal corresponding to a spurious signal. In some embodiments, the remaining other signal processing circuits including the DDS phase accumulator 10 may generate cancellation signals corresponding to multiple harmonics, or cancellation signals corresponding to multiple spurious signals, or cancellation signals corresponding to at least one harmonic and at least one spurious signal. The cancellation signals may be generated separately based on multiple signal processing circuits including the DDS phase accumulator 10, or may be generated by multiplexing based on one signal processing circuit including the DDS phase accumulator 10.

[0084] For example, if the value of H is ≥3, the remaining other paths including the signal processing circuit of the DDS phase accumulator 10 are used to generate cancellation signals corresponding to multiple harmonics, or cancellation signals corresponding to multiple spurious signals, or cancellation signals corresponding to at least one harmonic and at least one spurious signal. The maximum number of times the remaining other paths including the signal processing circuit of the DDS phase accumulator 10 are called is equal to the sum of the number of harmonics and the number of spurious signals, wherein each time the signal processing circuit is called, a cancellation signal corresponding to a harmonic or a cancellation signal corresponding to a spurious signal is generated.

[0085] Then, the cancellation signal and the original input signal are superimposed to obtain a digital test signal, the digital test signal is input into the DAC for digital-to-analog conversion to obtain an analog test signal, and the analog test signal is subjected to spectrum analysis to obtain an amplitude estimation value of the harmonic frequency or spurious frequency corresponding to any phase value.

[0086] In some embodiments, the analog signal to be tested can also be converted into digital form based on an analog-to-digital converter to obtain a corresponding digital signal, and then the processor or external data processing device transforms the digital signal from the time domain to the frequency domain, such as performing a fast Fourier transform (FFT) process, and then analyzing the transformed frequency domain signal to obtain an amplitude estimation value corresponding to the harmonic frequency or spurious frequency in the analog signal to be tested. In some embodiments, the analog signal to be tested can also be directly input into a spectrum analyzer for spectrum analysis.

[0087] Finally, the cancellation signals corresponding to the phase values ​​in the corresponding phase table to be measured are traversed to obtain the amplitude estimation values ​​corresponding to the phase values.

[0088] In some embodiments, when performing a traversal test, the signal processing channel may also include a phase counter and a memory, the phase counter is used to accumulate and count the M phase values ​​corresponding to the cancellation signal during traversal, and the memory is used to store the corresponding phase values ​​in the phase table to be tested and their corresponding amplitude estimation values.

[0089] Among them, when performing traversal testing on the phase table to be tested based on the coarse interval, M1 amplitude estimation values ​​corresponding to M1 phase values ​​can be obtained, and the target phase corresponding to the amplitude estimation target value among the M1 amplitude estimation values ​​can be determined.

[0090] In some embodiments, when determining the target phase corresponding to the amplitude estimation target value among the M1 amplitude estimation values, the average value of the M1 amplitude estimation values ​​is obtained, and it is determined whether the average value is greater than the harmonic amplitude value of the original input signal. If greater, the phase corresponding to the maximum value among the M1 amplitude estimation values ​​is obtained and used as the target phase corresponding to the amplitude estimation target value. If less, the phase corresponding to the minimum value among the M1 amplitude estimation values ​​is obtained and used as the target phase corresponding to the amplitude estimation target value.

[0091] In this embodiment, since the amplitude estimation value changes with the phase value in a certain regular pattern, such as a sine wave-like change, first increasing and then decreasing, or first decreasing and then increasing, the maximum value in the amplitude estimation value is usually also the maximum value, and the minimum value in the amplitude estimation value is usually also the minimum value. Therefore, when the average value of the M1 amplitude estimation values ​​is greater than the harmonic amplitude value of the original input signal, it means that the M1 amplitude estimation values ​​are closer to the maximum value, and vice versa, when the average value of the M1 amplitude estimation values ​​is less than the harmonic amplitude value of the original input signal, it means that the M1 amplitude estimation values ​​are closer to the minimum value. Thus, the phase range where the phase to be estimated is located is determined according to the comparison result, and the phase search range is reduced, thereby reducing the complexity of the phase estimation process.

[0092] A fine-interval phase table to be measured is generated according to the target phase; wherein the fine-interval phase table to be measured includes M2 phase values, and M2 is a positive integer greater than 1.

[0093] In some embodiments, M2 uniformly distributed phase values ​​may be selected based on the front and rear range of the target phase to generate a fine-interval phase table to be measured, for example, with the target phase as the center position, and M2 uniformly distributed phase values ​​are selected in a preset range greater than the target phase and a preset range less than the target phase, respectively. In some embodiments, the phase interval between the M1 phase values ​​is greater than the phase interval between the M2 phase values, that is, the intervals between the phase values ​​in the fine-interval phase table to be measured are smaller, so that a more accurate phase estimate can be determined. In some embodiments, the phase range included in the M1 phase values ​​is greater than the phase range included in the M2 phase values, and the value of M1 is smaller than the value of M2, that is, the number of phase values ​​in the fine-interval phase table to be measured is greater and the phase range included is smaller, so that a more accurate phase estimate can be determined.

[0094] The above traversal test is performed based on the fine-interval phase table to be tested to obtain M2 amplitude estimation values ​​corresponding to M2 phase values, and the phase corresponding to the minimum amplitude estimation value is obtained based on the M2 amplitude estimation values ​​and the corresponding phase values.

[0095] In some embodiments, when obtaining the phase corresponding to the minimum amplitude estimation value based on M2 amplitude estimation values ​​and the corresponding phase values, if the phase corresponding to the minimum value among the M1 amplitude estimation values ​​is obtained, the phase corresponding to the minimum value among the M2 amplitude estimation values ​​is determined to be the phase corresponding to the minimum amplitude estimation value. If the phase corresponding to the maximum value among the M1 amplitude estimation values ​​is obtained, the phase corresponding to the maximum value among the M2 amplitude estimation values ​​is obtained, and the phase corresponding to the minimum amplitude estimation value is determined based on the phase corresponding to the maximum value.

[0096] In some embodiments, when determining the phase corresponding to the minimum value of the amplitude estimate based on the phase corresponding to the maximum value, the phase corresponding to the maximum value among the M2 amplitude estimates is shifted by a preset phase so that the shifted phase corresponds to the minimum value of the amplitude estimate, and the shifted phase is output, where the shifted phase is the phase corresponding to the minimum value of the amplitude estimate. In this embodiment, since the phase corresponding to the maximum value or the local maximum value is in antiphase with the phase corresponding to the local minimum value, the phase corresponding to the minimum value of the amplitude estimate can be determined based on the maximum value or the local maximum value.

[0097] In some embodiments, if the phase value corresponding to the maximum value among the M2 amplitude estimation values ​​is between [0, π), the preset phase of the translation is π, and if the phase value corresponding to the maximum value among the M2 amplitude estimation values ​​is between [π, 2π), the preset phase of the translation is -π. The phase after translation is the sum of the phase value corresponding to the maximum value and the preset phase of the translation. In this embodiment, based on the phase range in which the phase value corresponding to the maximum value is located, different addition and subtraction operations are performed on the phase value corresponding to the maximum value, so that the phase after translation is always a positive value, which is convenient for subsequent use and calculation.

[0098] In the above embodiments, although the phase of the final estimated output is called the estimated phase of the harmonic or the estimated phase of the spurious, it should be clear that it is the antiphase of the phase of the real harmonic, and there is a phase difference of π between it and the phase of the real harmonic. Therefore, the phase of the output can be directly used when generating the cancellation signal corresponding to the harmonic or spurious.

[0099] In the above embodiment, the phase estimation process first searches for the phase table to be measured with a coarse interval, and then searches for the phase table to be measured with a fine interval based on the search results of the phase table to be measured with a coarse interval. This search idea can reduce the overall number of searches (corresponding to the number of phase values ​​in the phase table) while also obtaining a more accurate phase estimation result.

[0100] In the above embodiments, the signal processing channel can process at least one harmonic or spurious wave, or can process harmonics and spurious waves. In some embodiments, the total number L of harmonics and spurious waves to be processed is obtained, and twice the total number L is divided by H-1 and rounded up to obtain the number of times the signal processing channel is called. According to the number of times the signal processing channel is called, the phase counter and the memory are reset to clear the amplitude estimation value and its corresponding phase value stored in the history. Among them, H-1 harmonics or spurious waves can be processed every two times the signal processing channel is called.

[0101] In some embodiments, a reset operation is required when calling a signal processing channel to avoid mutual interference. For example, if it is detected that the signal processing channel has been cleared before being called, no reset is required; if it is detected that the signal processing channel has not been cleared before being called, a reset and clearing operation is required. Therefore, when resetting the phase counter and the memory according to the number of times the signal processing channel is called, the number of times the signal processing channel is called can be the same as the number of resets. For example, a reset and clearing operation is performed before each call to the signal processing channel. When it is detected that the signal processing channel has been cleared, the number of times the signal processing channel is called can be less than the number of resets.

[0102] The following is an illustration through specific embodiments.

[0103] Assume that the input frequency f_in of the original input signal is 30MHz, then calculate the frequency control word based on the input frequency f_in, and then calculate the phase sequence output by the DDS phase accumulator 10, address the sine lookup table according to the output phase sequence, obtain the full-scale signal, and then multiply the full-scale signal by the attenuation coefficient α=0.95 and input it to the DAC. Assuming that the harmonics to be estimated are the second harmonic, the third harmonic and the fourth harmonic, their frequencies are f_2=2*f_in=60MHz, f_3=3*f_in=90MHz, and f_4=4*f_in=120MHz respectively.

[0104] Assume that the number of phase points selected from the phase range [0,π] is M1 = 5, then the phase value of the mth phase point in the coarse interval phase table to be measured is set to In some embodiments, when obtaining the estimated phase empirical value, a phase value is selected near the estimated phase empirical value. For example, if the estimated phase empirical value is π / 4, then

[0105] First, the original input signal is input into the DAC, and the ADC is used to sample the output signal of the DAC, and then FFT is performed to record the original amplitude estimation values ​​of the second harmonic f_2, third harmonic f_3, and fourth harmonic f_4 respectively. Since the frequency of FFT may have deviations, the recorded amplitude can be selected as the maximum value within a given range around f_2, f_3, and f_4.

[0106] A traversal test is performed based on the coarse-interval phase table to be tested, that is, a phase value θ_m is selected in sequence from the coarse-interval phase table to be tested, and the corresponding harmonic cancellation signal is generated according to the second harmonic frequency f_2, the third harmonic frequency f_3, the fourth harmonic frequency f_4 and the corresponding original harmonic amplitude value. The generation process of the cancellation signal is similar to that of the original input signal. Its input frequency is the frequency of the harmonic, and then the frequency control word and the phase sequence are calculated in sequence, and then the full-scale signal, that is, the original amplitude estimation value, is obtained. The cancellation signal is superimposed on the original input signal and input into the DAC to obtain the simulated phase test signal output by the DAC.

[0107] When the original input signal is a signal related to θ_m, the corresponding simulated phase test signal is input into the DAC, and the ADC is used to sample the output signal of the DAC, and then FFT is performed, and the amplitude A_(m_2) corresponding to the second harmonic, the amplitude A_(m_3) corresponding to the third harmonic, and the amplitude A_(m_4) corresponding to the fourth harmonic are recorded, where m = 1, 2, ..., 5. For example, please refer to Figure 4 , where the phase numbers are divided into the corresponding phase values ​​in the phase table to be measured with a coarse interval of 1-5, and the estimated amplitude value is the recorded amplitude of the corresponding phase number.

[0108] After completing the traversal test of all phase points in the phase table to be tested with coarse intervals, the average values ​​of the amplitudes of all recorded second harmonics, third harmonics, and fourth harmonics A_(m_2) / A_(m_3) / A_(m_4), m=1,2,…,5 are respectively solved. If the average value is less than the corresponding original amplitude estimation value, the phase value corresponding to the minimum value is output, and the flag flag is set to 1. If the average value is greater than the corresponding original amplitude estimation value, the phase value corresponding to the maximum value is output, and the flag flag is set to 0.

[0109] Based on the phase value output above, the number of phase points M2=10 is selected, and the selection method is uniformly spaced sampling to form a fine-interval phase table to be tested. Due to the periodicity of the phase, if the phase value output above is a boundary point of [0,π], the phase can be extended and then uniformly spaced sampling can be performed. The above traversal test is performed based on the fine-interval phase table to be tested. For example, please refer to Figure 5 , where the phase numbers are divided into the corresponding phase values ​​in the phase table to be measured with fine intervals of 1-10, and the estimated amplitude value is the recorded amplitude of the corresponding phase number.

[0110] Determine whether flag is 1. If so, select the fine-interval phase table to be tested to perform the minimum value recorded during the phase traversal process, and select the corresponding phase value for output. If not, select the fine-interval phase table to be tested to perform the maximum value recorded during the phase traversal, and select the corresponding phase value to increase or subtract π for output.

[0111] The above is the description of the signal processing channel.

[0112] Some embodiments provide a phase estimation method, which can be applied to the above-mentioned signal processing channel. Figure 6 , the harmonic phase acquisition method includes the following steps:

[0113] Step 100: Obtain the original input signal and the coarse-interval phase table to be measured, wherein the coarse-interval phase table to be measured includes M1 phase values, where M1 is a positive integer greater than 1.

[0114] Step 110: Perform traversal test based on the coarse-interval phase table to be tested. Generate a cancellation signal based on the harmonic frequency and harmonic amplitude value of the harmonics of the original input signal or based on the spurious frequency and spurious amplitude value of the spurious of the original input signal, and any phase value in the corresponding phase table to be tested; superimpose the cancellation signal and the original input signal to obtain a digital signal to be tested; input the digital signal to be tested into a DAC for digital-to-analog conversion to obtain an analog signal to be tested, perform spectrum analysis on the analog signal to be tested, and obtain an amplitude estimation value of the harmonic frequency or spurious frequency corresponding to any phase value; traverse the cancellation signal corresponding to each phase value in the corresponding phase table to be tested to obtain each amplitude estimation value corresponding to each phase value.

[0115] Step 120: Determine the target phase corresponding to the amplitude estimation target value among the M1 amplitude estimation values, and generate a fine-interval phase table to be tested according to the target phase. In which, a traversal test is performed based on the coarse-interval phase table to be tested to obtain M1 amplitude estimation values ​​corresponding to the M1 phase values.

[0116] Step 130: Perform a traversal test based on the fine-interval phase table to be tested to obtain M2 amplitude estimation values ​​corresponding to M2 phase values, and obtain the phase corresponding to the minimum amplitude estimation value based on the M2 amplitude estimation values ​​and the corresponding phase values.

[0117] Some embodiments provide a computer-readable storage medium having a program stored thereon, which can be executed by a processor to implement the above-mentioned phase estimation method.

[0118] Those skilled in the art will appreciate that all or part of the functions of the various methods in the above-mentioned embodiments can be implemented by hardware or by computer programs. When all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can be stored in a computer-readable storage medium, and the storage medium can include: a read-only memory, a random access memory, a disk, an optical disk, a hard disk, etc., and the program is executed by a computer to implement the above-mentioned functions. For example, the program is stored in the memory of the device, and when the program in the memory is executed by the processor 30, all or part of the above-mentioned functions can be implemented. In addition, when all or part of the functions in the above-mentioned embodiments are implemented by computer programs, the program can also be stored in a storage medium such as a server, another computer, a disk, an optical disk, a flash disk or a mobile hard disk, and can be downloaded or copied and saved in the memory of the local device, or the system of the local device is updated, and when the program in the memory is executed by the processor 30, all or part of the functions in the above-mentioned embodiments can be implemented.

[0119] The above specific examples are used to illustrate the present application, which is only used to help understand the present application and is not intended to limit the present application. For technicians in the technical field to which the present application belongs, they can also make some simple deductions, deformations or substitutions based on the ideas of the present application.

Claims

1. A phase estimation method, characterized in that: include: Obtaining an original input signal and a coarse-interval phase table to be measured; wherein the coarse-interval phase table to be measured includes M1 phase values, and M1 is a positive integer greater than 1; Perform traversal test on the phase table to be tested based on the coarse interval: generate a cancellation signal based on the harmonic frequency and harmonic amplitude value of the harmonics of the original input signal or based on the spurious frequency and spurious amplitude value of the spurious of the original input signal, and any phase value in the corresponding phase table to be tested; superimpose the cancellation signal and the original input signal to obtain a digital signal to be tested; input the digital signal to be tested into a DAC for digital-to-analog conversion to obtain an analog signal to be tested, perform spectrum analysis on the analog signal to be tested to obtain an amplitude estimation value of the harmonic frequency or spurious frequency corresponding to any phase value; traverse the cancellation signal corresponding to each phase value in the corresponding phase table to be tested to obtain each amplitude estimation value corresponding to each phase value; Acquire M1 amplitude estimation values ​​corresponding to the M1 phase values, and determine a target phase corresponding to an amplitude estimation target value among the M1 amplitude estimation values; Generate a fine-interval phase table to be measured according to the target phase; wherein the fine-interval phase table to be measured includes M2 phase values, and M2 is a positive integer greater than 1; The traversal test is performed based on the fine-interval phase table to be tested to obtain M2 amplitude estimation values ​​corresponding to the M2 phase values, and the phase corresponding to the minimum amplitude estimation value is obtained based on the M2 amplitude estimation values ​​and the corresponding phase values.

2. The phase estimation method according to claim 1, characterized in that: The obtaining a target phase corresponding to the amplitude estimation target value in the M1 amplitude estimation values ​​includes: Obtaining an average value of the M1 amplitude estimation values; Determining whether the average value is greater than the harmonic amplitude value or the spurious amplitude value of the original input signal; If it is greater than, then the phase corresponding to the maximum value among the M1 amplitude estimation values ​​is obtained and used as the target phase corresponding to the amplitude estimation target value; If it is less than, the phase corresponding to the minimum value among the M1 amplitude estimation values ​​is obtained and used as the target phase corresponding to the amplitude estimation target value.

3. The phase estimation method according to claim 2, characterized in that: The obtaining, based on the M2 amplitude estimation values ​​and the corresponding phase values, a phase corresponding to the minimum amplitude estimation value, comprises: If the phase corresponding to the minimum value among the M1 amplitude estimation values ​​is obtained, then the phase corresponding to the minimum value among the M2 amplitude estimation values ​​is determined to be the phase corresponding to the minimum amplitude estimation value; If the phase corresponding to the maximum value among the M1 amplitude estimation values ​​is obtained, the phase corresponding to the maximum value among the M2 amplitude estimation values ​​is obtained, and the phase corresponding to the minimum amplitude estimation value is determined based on the phase corresponding to the maximum value.

4. The phase estimation method according to claim 3, characterized in that: The determining, according to the phase corresponding to the maximum value, the phase corresponding to the amplitude estimation minimum value comprises: Shifting the phase corresponding to the maximum value among the M2 amplitude estimation values ​​by a preset phase so that the shifted phase corresponds to the amplitude estimation minimum value; The shifted phase is output, where the shifted phase is the phase corresponding to the minimum value of the amplitude estimation.

5. The phase estimation method according to claim 4, characterized in that: If the phase value corresponding to the maximum value among the M2 amplitude estimation values ​​is between [0, π), the preset phase of the translation is π; If the phase value corresponding to the maximum value among the M2 amplitude estimation values ​​is between [π, 2π), the preset phase of the translation is -π.

6. The phase estimation method according to claim 1, wherein: The coarse interval phase table to be measured is generated in the following manner: Dividing M1 phase values ​​from the phase range [a, a+π] according to a preset first phase interval to generate the coarse interval phase table to be measured; wherein a is a phase greater than or equal to 0 and less than or equal to π; or, An estimated phase empirical value is obtained, and a phase table to be measured of the coarse interval is generated according to the estimated phase empirical value.

7. The phase estimation method according to claim 6, characterized in that: The obtaining of the estimated phase empirical value comprises: Obtaining a phase estimation value corresponding to a harmonic frequency or a spurious frequency of a known input signal after the known input signal is input into the DAC; Based on the frequency relationship between the known input signal and the original input signal, the estimated phase empirical value corresponding to the harmonic of the original input signal is obtained from the phase estimation value; or, Based on the frequency relationship between the known input signal and the original input signal, and the frequency relationship between the spurious frequency of the known input signal and the spurious frequency of the original input signal, the estimated phase empirical value corresponding to the spur of the original input signal is obtained from the phase estimation value.

8. The phase estimation method according to claim 6, characterized in that: The step of generating the coarse-interval phase table to be measured according to the estimated phase empirical value comprises: According to the preset second phase interval, M2 phase values ​​are divided from the nearby range [bc, b+c] of the estimated phase empirical value to generate the coarse interval phase table to be tested; wherein b is the estimated phase empirical value and c is the preset phase value.

9. The phase estimation method according to claim 1, characterized in that: The phase interval between the M1 phase values ​​is greater than the phase interval between the M2 phase values; and / or the phase range included in the M1 phase values ​​is greater than the phase range included in the M2 phase values, and the value of M1 is smaller than the value of M2.

10. The phase estimation method according to claim 1, characterized in that: The performing spectrum analysis on the simulated signal to be tested comprises: Inputting the simulated signal to be tested into a spectrum analyzer for spectrum analysis; or, The simulated signal to be tested is input into an analog-to-digital converter for analog-to-digital conversion and then subjected to Fourier transform processing for spectrum analysis.

11. The phase estimation method according to claim 1, characterized in that: For harmonics, the corresponding frequency of the cancellation signal is the harmonic frequency of the harmonics of the original input signal; for spurious signals, the corresponding frequency of the cancellation signal is the spurious frequency of the spurious signals of the original input signal.

12. A signal processing channel, characterized in that: Used to implement the phase estimation method according to any one of claims 1 to 11, the signal processing channel comprises at least two DDS phase accumulators, and the signal processing circuits where each of the DDS phase accumulators is located are connected in parallel, and the number of parallel connections is defined as an integer H greater than 1; One of the signal processing circuits including a DDS phase accumulator generates the original input signal after phase-to-amplitude conversion processing and attenuation processing, and the remaining signal processing circuits including DDS phase accumulators generate cancellation signals corresponding to harmonics or cancellation signals corresponding to spurious signals based on the corresponding phase values ​​in the phase table to be measured through phase-to-amplitude conversion processing.

13. The signal processing channel according to claim 12, characterized in that: The signal processing channel is used to process multiple harmonics, or multiple spurious signals, or at least one harmonic and at least one spurious signal. If the value of H is ≥3, the remaining signal processing circuits including the DDS phase accumulator are used to generate cancellation signals corresponding to multiple harmonics, or cancellation signals corresponding to multiple spurious signals, or cancellation signals corresponding to at least one harmonic and at least one spurious signal respectively. The maximum number of times the remaining signal processing circuits including the DDS phase accumulator are called is equal to twice the sum of the number of harmonics and the number of spurious signals.

14. The signal processing channel according to claim 13, characterized in that The signal processing channel also includes a phase counter and a memory. The phase counter is used to accumulate and count the M1 phase values ​​or M2 phase values ​​corresponding to the cancellation signal during traversal, and the memory is used to store the corresponding phase values ​​in the phase table to be tested and the corresponding amplitude estimation values.

15. The signal processing channel according to claim 14, characterized in that Obtain the total number L of harmonics and spurious signals to be processed; Divide twice the total number L by H-1 and round up to get the number of times the signal processing channel is called; According to the number of times the signal processing channel is called, the phase counter and the memory are reset to clear the amplitude estimation value and the corresponding phase value stored in history.

16. A computer-readable storage medium, characterized in that: The medium stores a program, which can be executed by a processor to implement the method according to any one of claims 1 to 11.

Images