CN117433622A - Signal processing method, device, equipment and storage medium - Google Patents
Signal processing method, device, equipment and storage medium Download PDFInfo
- Publication number
- CN117433622A CN117433622A CN202311380114.9A CN202311380114A CN117433622A CN 117433622 A CN117433622 A CN 117433622A CN 202311380114 A CN202311380114 A CN 202311380114A CN 117433622 A CN117433622 A CN 117433622A
- Authority
- CN
- China
- Prior art keywords
- signal
- initial
- determining
- amplitude ratio
- interference
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000003672 processing method Methods 0.000 title claims abstract description 38
- 238000004804 winding Methods 0.000 claims abstract description 20
- 238000000034 method Methods 0.000 claims description 29
- 238000004422 calculation algorithm Methods 0.000 claims description 19
- 230000010354 integration Effects 0.000 claims description 16
- 238000005070 sampling Methods 0.000 claims description 11
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 230000009466 transformation Effects 0.000 claims 2
- 230000001427 coherent effect Effects 0.000 abstract description 10
- 238000001514 detection method Methods 0.000 abstract description 4
- 238000004590 computer program Methods 0.000 description 17
- 238000010586 diagram Methods 0.000 description 15
- 230000006870 function Effects 0.000 description 13
- 238000005259 measurement Methods 0.000 description 13
- 238000012545 processing Methods 0.000 description 12
- 238000001228 spectrum Methods 0.000 description 12
- 230000010287 polarization Effects 0.000 description 11
- 230000015654 memory Effects 0.000 description 10
- 238000004364 calculation method Methods 0.000 description 7
- 230000008859 change Effects 0.000 description 6
- 230000010355 oscillation Effects 0.000 description 6
- 238000004891 communication Methods 0.000 description 5
- 230000004069 differentiation Effects 0.000 description 4
- 238000012937 correction Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- 230000008569 process Effects 0.000 description 3
- 238000009825 accumulation Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000006073 displacement reaction Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000014509 gene expression Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 230000007704 transition Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000036541 health Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000000758 substrate Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01H—MEASUREMENT OF MECHANICAL VIBRATIONS OR ULTRASONIC, SONIC OR INFRASONIC WAVES
- G01H9/00—Measuring mechanical vibrations or ultrasonic, sonic or infrasonic waves by using radiation-sensitive means, e.g. optical means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
- G06F18/15—Statistical pre-processing, e.g. techniques for normalisation or restoring missing data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2218/00—Aspects of pattern recognition specially adapted for signal processing
- G06F2218/02—Preprocessing
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Biology (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Probability & Statistics with Applications (AREA)
- Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)
Abstract
The disclosure provides a signal processing method, a device, equipment and a storage medium, which can be applied to the technical field of laser coherent detection. The signal processing method comprises the following steps: performing differential operation on a winding phase signal of an interference signal based on the current iteration times to obtain a differential signal of the winding phase signal, wherein the interference signal is obtained by heterodyne interference between a reference signal and a signal to be identified; based on the current iteration times, performing integral operation on the compensation signal of the differential signal to obtain an integral signal of the compensation signal; determining an initial unwrapping signal according to the integrated signal and a fitting signal of the integrated signal; determining a current amplitude ratio based on the initial unwind signal; and under the condition that the current amplitude ratio represents that the initial unwrapping signal meets the requirement, determining the signal to be identified according to the initial unwrapping signal corresponding to the previous iteration number.
Description
Technical Field
The present disclosure relates to the field of laser coherent detection technology, and more particularly, to a signal processing method, apparatus, device, and storage medium.
Background
The laser coherent vibration measuring system is a non-contact vibration measuring device, has the characteristics of high precision, high spatial resolution, non-contact measurement and the like, and is widely applied to the fields of automobile industry detection, structural health monitoring, precise device manufacturing and the like. In the laser coherent vibration measurement system, phase demodulation is one of key technical links of the laser coherent vibration measurement system. In the phase demodulation method, an arctangent function operation is often used to extract the phase, and the demodulated output signal thereof is often wrapped due to the characteristic of the arctangent function.
In the prior art, the problem that the recovery range is small and the recovery result precision is low exists when the output signal which is coiled is uncoiled.
Disclosure of Invention
In view of the foregoing, the present disclosure provides a signal processing method, apparatus, device, medium, and program product.
According to a first aspect of the present disclosure, there is provided a signal processing method comprising: performing differential operation on a winding phase signal of an interference signal based on the current iteration times to obtain a differential signal of the winding phase signal, wherein the interference signal is obtained by heterodyne interference between a reference signal and a signal to be identified; based on the current iteration times, performing integral operation on the compensation signal of the differential signal to obtain an integral signal of the compensation signal; determining an initial unwrapping signal according to the integrated signal and a fitting signal of the integrated signal; determining a current amplitude ratio based on the initial unwind signal; and under the condition that the current amplitude ratio represents that the initial unwrapping signal meets the requirement, determining the signal to be identified according to the initial unwrapping signal corresponding to the previous iteration number.
According to an embodiment of the present disclosure, determining a current amplitude ratio from an initial unwind signal includes: performing Fourier transform on the initial unwrapped signal to obtain a transformed signal; calculating the sum of the amplitudes of the components of the frequency of the transformed signal in a first predetermined frequency range, resulting in a first frequency sum, wherein the first predetermined frequency range may be determined from the frequency of the interference signal and the nyquist frequency of the initial unwrapped signal; calculating the sum of the amplitudes of the components of the transformed signal whose frequencies lie in a first predetermined frequency range, obtaining a second frequency sum, wherein the second predetermined frequency range can be based on the nyquist frequency of the initial unwrapped signal; and calculating the ratio of the first frequency sum to the second frequency sum to obtain the current amplitude ratio.
According to an embodiment of the present disclosure, determining an initial unwind signal from an integrated signal and a fitted signal of the integrated signal comprises: performing polynomial fitting on the integrated signal based on the current iteration times to obtain a fitted signal; and calculating the difference between the integrated signal and the fitting signal to obtain an initial unwrapping signal.
According to an embodiment of the present disclosure, the signal processing method further includes: according to the current iteration times, determining a historical amplitude ratio corresponding to the last iteration; and under the condition that the current amplitude ratio is larger than the historical amplitude ratio, determining that the initial unwinding signal meets the requirement.
According to an embodiment of the present disclosure, the signal processing method further includes: under the condition that the current amplitude ratio is smaller than the historical amplitude ratio, determining that the initial unwinding signal does not meet the requirement; and accumulating the current iteration times by 1 under the condition that the current amplitude ratio represents that the initial unwrapping signal does not meet the requirement, wherein the initial value of the current iteration times is 0.
According to an embodiment of the present disclosure, the signal processing method further includes: and determining a compensation signal according to the differential signal and a count value, wherein the count value is determined according to the historical count value and the differential signal corresponding to the last iteration.
According to an embodiment of the present disclosure, the signal processing method further includes: determining the quadrature component of the interference signal using a demodulation algorithm; the wrapped-around phase signal of the interference signal is determined using an arctangent function based on the ratio of the quadrature components.
According to an embodiment of the present disclosure, determining a compensation signal from a differential signal and a count value includes: analyzing the difference value of adjacent sampling points of the differential signal to obtain the jump direction and amplitude; determining the current count value according to the direction and the amplitude of the jump; and summing the historical count value and the current count value to obtain the count value.
According to an embodiment of the present disclosure, the signal processing method further includes: and carrying out orthogonal decomposition on the interference signal to obtain a first orthogonal component and a second orthogonal component, wherein the first orthogonal component is determined according to the interference signal and the sine of the local oscillator, and the second orthogonal component is determined according to the interference signal and the cosine of the local oscillator.
A second aspect of the present disclosure provides a signal processing apparatus comprising:
the signal difference module is used for carrying out difference operation on the winding phase signals of the interference signals based on the current iteration times to obtain difference signals of the winding phase signals, wherein the interference signals are obtained by heterodyne interference of the reference signals and the signals to be identified;
the signal integration module is used for performing integration operation on the compensation signal of the differential signal based on the current iteration times to obtain an integrated signal of the compensation signal;
the unwrapping signal determining module is used for determining an initial unwrapping signal according to the integrated signal and the fitting signal of the integrated signal;
the ratio determining module is used for determining a current amplitude ratio based on the initial unwinding signal;
the signal to be identified determining module is used for determining the signal to be identified according to the initial unwrapping signal corresponding to the previous iteration number under the condition that the current amplitude ratio represents that the initial unwrapping signal meets the requirement.
A third aspect of the present disclosure provides an electronic device, comprising: one or more processors; and a memory for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to perform the signal processing method described above.
A fourth aspect of the present disclosure also provides a computer-readable storage medium having stored thereon executable instructions that, when executed by a processor, cause the processor to perform the above-described signal processing method.
A fifth aspect of the present disclosure also provides a computer program product comprising a computer program which, when executed by a processor, implements the above-described signal processing method.
According to embodiments of the present disclosure, signal downscaling and fixed compensation unwrapping are achieved by compressing the signal for multiple differentiation of the wrapped phase signal and unwrapping of fixed transition values. And (3) carrying out signal expansion and eliminating an error trend term through overshoot integral and polynomial fitting, judging whether an initial unwrapping signal corresponding to the current iteration number meets the requirement or not according to the change of the frequency spectrum amplitude ratio of the high-frequency band, and realizing the unwrapping of the signal. Because multiple differentiation and integration are adopted, the phase demodulation correction can be realized through algorithm design on the premise of not adding additional devices, the problem that the phase can only be recovered from the original signal with the phase change not exceeding pi rad in the traditional phase unwrapping method is solved, and the dynamic range of phase demodulation is improved. By judging whether the initial unwinding signal meets the requirement or not through the current amplitude ratio, the influence on the accuracy of the phase unwinding signal result caused by the manual selection of an iteration threshold in the traditional phase unwinding method is reduced, the application range of the laser coherent vibration measurement is widened, and the displacement measurement accuracy of the laser coherent vibration measurement system is effectively improved.
Drawings
The foregoing and other objects, features and advantages of the disclosure will be more apparent from the following description of embodiments of the disclosure with reference to the accompanying drawings, in which:
fig. 1 schematically illustrates an application scenario of a signal processing method according to an embodiment of the present disclosure;
fig. 2 schematically illustrates a flow chart of a signal processing method according to an embodiment of the present disclosure;
fig. 3 schematically illustrates a flow chart of a signal processing method according to an embodiment of the disclosure;
fig. 4 (a) to (f) schematically show demodulation time domain waveform diagrams of a signal processing method according to an embodiment of the present disclosure, respectively;
fig. 5 (a) to (d) schematically show demodulation spectrograms of a signal processing method according to an embodiment of the present disclosure, respectively;
fig. 6 schematically shows a block diagram of a signal processing apparatus according to an embodiment of the present disclosure; and
fig. 7 schematically shows a block diagram of an electronic device adapted for a signal processing method according to an embodiment of the disclosure.
Detailed Description
Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. It should be understood that the description is only exemplary and is not intended to limit the scope of the present disclosure. In the following detailed description, for purposes of explanation, numerous specific details are set forth in order to provide a thorough understanding of the embodiments of the present disclosure. It may be evident, however, that one or more embodiments may be practiced without these specific details. In addition, in the following description, descriptions of well-known structures and techniques are omitted so as not to unnecessarily obscure the concepts of the present disclosure.
The terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the disclosure. The terms "comprises," "comprising," and/or the like, as used herein, specify the presence of stated features, steps, operations, and/or components, but do not preclude the presence or addition of one or more other features, steps, operations, or components.
All terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art unless otherwise defined. It should be noted that the terms used herein should be construed to have meanings consistent with the context of the present specification and should not be construed in an idealized or overly formal manner.
Where expressions like at least one of "A, B and C, etc. are used, the expressions should generally be interpreted in accordance with the meaning as commonly understood by those skilled in the art (e.g.," a system having at least one of A, B and C "shall include, but not be limited to, a system having a alone, B alone, C alone, a and B together, a and C together, B and C together, and/or A, B, C together, etc.).
In the implementation process of the present disclosure, it is found that the existing signal phase unwrapping method calculates the difference between adjacent values according to the feature of wrapping the arctangent calculation result, determines whether a jump occurs in the arctangent value near the sampling point according to whether the difference exceeds a certain threshold, and compensates the sampling point according to the jump count value, where the compensation value is generally an integer multiple of pi rad. The unwrapping method can realize correction of the phase demodulation result to a certain extent, but can only recover signals with the phase change of no more than pi rad between adjacent sampling points.
Another phase Unwrapping method is a Differential-Unwrapping-Integral (DUI) algorithm, where the amplitude of the signal is suppressed by Differential operation in a single Unwrapping, and by selecting an appropriate Differential order, a wrapped signal with a phase change greater than pi rad between adjacent sampling points can be unwrapped correctly. And then the original phase is restored without additional noise and additional device cost through multiple integral operation and polynomial fitting compensation algorithm. However, the stopping condition of the DUI algorithm is that the spearman coefficient of the algorithm result is greater than a certain threshold epsilon close to 1, and in the discrete digital signal processing process, it is difficult to select a proper threshold epsilon, and in practical application, the threshold selection error also causes an error of the unwrapping result, and affects the phase unwrapping effect.
The embodiment of the disclosure provides a signal processing method, which comprises the following steps: performing differential operation on a winding phase signal of an interference signal based on the current iteration times to obtain a differential signal of the winding phase signal, wherein the interference signal is obtained by heterodyne interference between a reference signal and a signal to be identified; based on the current iteration times, performing integral operation on the compensation signal of the differential signal to obtain an integral signal of the compensation signal; determining an initial unwrapping signal according to the integrated signal and a fitting signal of the integrated signal; determining a current amplitude ratio based on the initial unwind signal; and under the condition that the current amplitude ratio represents that the initial unwrapping signal meets the requirement, determining the signal to be identified according to the initial unwrapping signal corresponding to the previous iteration number.
Fig. 1 schematically illustrates an application scenario of a signal processing method according to an embodiment of the present disclosure.
As shown in fig. 1, an application scenario 100 according to this embodiment may include a laser source 101, a polarization splitting prism 102, a polarization splitting prism 103, a vibrating object 104, a polarization splitting prism 105, a polarization splitting prism 106, an acousto-optic frequency shifter 107, and a heterodyne coherence detector 108. The laser source 101 may be used to emit laser light, the polarization beam splitter prism 102 may be used to split the laser emitted by the laser source 101 into local oscillation light and measurement light, in fig. 1, the local oscillation light path is a dotted line, the measurement light path is a solid line, the polarization beam splitter prism 103 may be used to reflect light reflected by the measurement light to the vibration object, the vibration object 104 may vibrate, shift the frequency of the measurement light to obtain a signal to be identified, the polarization beam splitter prism 105 may be used to reflect the reference signal, the polarization beam splitter prism 106 may be used to reflect the local oscillation light, the acousto-optic frequency shifter 107 may be used to make the local oscillation light generate a fixed frequency shift to obtain the reference signal, and the heterodyne coherence detector 108 may be used to receive an interference signal generated by interference between the detection reference signal and the signal to be identified.
It should be understood that the number of laser sources, polarization splitting chills, vibrating objects, and acousto-optic frequency shifters in fig. 1 is merely illustrative. Any number of laser sources, polarization splitting chills, vibrating objects, and acousto-optic frequency shifters may be provided as desired.
The signal processing method of the disclosed embodiment will be described in detail below with reference to fig. 2 to 5 based on the scenario described in fig. 1.
Fig. 2 schematically shows a flow chart of a signal processing method according to an embodiment of the present disclosure.
As shown in fig. 2, the method 200 includes operations S210-S250.
In operation S210, a difference operation is performed on the wrapped phase signal of the interference signal based on the current iteration number, to obtain a difference signal of the wrapped phase signal, where the interference signal is obtained by heterodyne interference between the reference signal and the signal to be identified.
According to the embodiment of the disclosure, a laser beam emitted by a laser source can be divided into local oscillation light and measuring light through a polarization beam splitter, wherein the local oscillation light is subjected to fixed frequency shift after passing through an acousto-optic frequency shifter, so as to obtain a reference signal; the measuring light irradiates the surface of the vibration object, and Doppler frequency shift can be generated on the measuring light due to the vibration of the vibration object, so that a signal to be identified is obtained. The reference signal and the interference signal are coherent by heterodyne, and the interference signal can be obtained.
According to embodiments of the present disclosure, a wrapped-around phase signal of an interference signalCan be determined from the quadrature component of the interference signal. Performing differential operation with the frequency equal to the current iteration frequency n on the interference signal to obtain a differential signal D (t), wherein the differential signal D (t) is shown in a formula (1):
When the current iteration number is 0, the differential signal is 0 times of the difference of the winding phase signal, i.e. the winding phase signal is used as the differential signal.
In operation S220, the compensation signal of the differential signal is integrated based on the current iteration number, to obtain an integrated signal of the compensation signal.
According to the embodiment of the disclosure, the compensation signal u (t) is determined according to the jump and direction generated by the differential signal, and the signal is obtained after the differential signal is subjected to phase compensation. And (3) performing integral operation with the number of times equal to the current iteration number on the compensation signal to obtain an integral signal F (t), wherein the integral signal F (t) is shown in a formula (2):
F(t)＝I n [u(t)] (2)
in the case that the current iteration number is 0, the integration signal is 0 times of integration of the compensation signal, i.e., the compensation signal is taken as the integration signal.
In operation S230, an initial unwind signal is determined from the integrated signal and a fitting signal of the integrated signal.
According to the embodiment of the disclosure, when the integration operation is performed on the compensation signal of the differential signal, each time the integration is performed, the integration signal is increased by a constant C, and after n times of integration, the integration signal is increased by an error trend term, which is n times of polynomials. To eliminate the error trend term, the error trend term may be fitted by obtaining a fitting signal P (t) by a fitting operation, and then using the difference between the integrated signal and the fitting signal as an initial value Unwind signal
In operation S240, a current amplitude ratio is determined based on the initial unwind signal.
According to embodiments of the present disclosure, the initial unwind signal may beAnd converting the frequency-domain signal into a frequency-domain signal, and determining the current amplitude ratio according to the sum of the amplitudes of the high-frequency component spectrums in the frequency-domain signal and the sum of the amplitudes of all the component spectrums in the frequency-domain signal.
In operation S250, in the case that the current amplitude ratio represents that the initial unwrapped signal meets the requirement, the signal to be identified is determined according to the initial unwrapped signal corresponding to the previous iteration number.
According to the embodiment of the disclosure, in the case that the current amplitude ratio represents that the initial unwrapping signal is less distorted than the signal to be identified, the signal to be identified can be determined according to the initial unwrapping signal and the reference signal, so as to determine the vibration condition of the vibration object.
According to embodiments of the present disclosure, signal downscaling and fixed compensation unwrapping are achieved by compressing the signal for multiple differentiation of the wrapped phase signal and unwrapping of fixed transition values. And (3) carrying out signal expansion and eliminating an error trend term through overshoot integral and polynomial fitting, judging whether an initial unwrapping signal corresponding to the current iteration number meets the requirement or not according to the change of the frequency spectrum amplitude ratio of the high-frequency band, and realizing the unwrapping of the signal. Because multiple differentiation and integration are adopted, the phase demodulation correction can be realized through algorithm design on the premise of not adding additional devices, the problem that the phase can only be recovered from the original signal with the phase change not exceeding pi rad in the traditional phase unwrapping method is solved, and the dynamic range of phase demodulation is improved. By judging whether the initial unwinding signal meets the requirement or not through the current amplitude ratio, the influence on the accuracy of the phase unwinding signal result caused by the manual selection of an iteration threshold in the traditional phase unwinding method is reduced, the application range of the laser coherent vibration measurement is widened, and the displacement measurement accuracy of the laser coherent vibration measurement system is effectively improved.
According to an embodiment of the present disclosure, the signal processing method further includes: and carrying out orthogonal decomposition on the interference signal to obtain a first orthogonal component and a second orthogonal component, wherein the first orthogonal component is determined according to the interference signal and the sine of the local oscillator, and the second orthogonal component is determined according to the interference signal and the cosine of the local oscillator.
In accordance with embodiments of the present disclosure, the quadrature components I (t) and Q (t) of the interference signal may be determined using an IQ demodulation algorithm, but is not limited thereto, and other demodulation algorithms may be used to determine the quadrature components of the interference signal. Wherein the first quadrature component Q (t) is determined from the component of the interference signal and the sine of the local oscillator and the second quadrature component I (t) is determined from the other component of the interference signal and the cosine of the local oscillator
According to an embodiment of the present disclosure, the signal processing method further includes: determining the quadrature component of the interference signal using a demodulation algorithm; the wrapped-around phase signal of the interference signal is determined using an arctangent function based on the ratio of the quadrature components.
According to embodiments of the present disclosure, based on the ratio of the quadrature components, the winding phase signal may be determined according to equation (3)
According to an embodiment of the present disclosure, the first quadrature component Q (t) is determined from a component of the interference signal and the sine of the local oscillator, and the second quadrature component I (t) is determined from another component of the interference signal and the cosine of the local oscillator.
According to an embodiment of the present disclosure, the signal processing method further includes: and determining a compensation signal according to the differential signal and a count value, wherein the count value is determined according to the historical count value and the differential signal corresponding to the last iteration.
According to an embodiment of the present disclosure, according to the count value k n (t) and a differential signal D (t), the compensation signal u (t) can be determined according to equation (4):
u(t)＝D(t)+2π*k n (t) (4)
according to an embodiment of the present disclosure, determining a compensation signal from a differential signal and a count value includes: analyzing the difference value of adjacent sampling points of the differential signal to obtain the jump direction and amplitude; determining the current count value according to the direction and the amplitude of the jump; and summing the historical count value and the current count value to obtain the count value.
According to the embodiment of the disclosure, the difference value of adjacent sampling points of the differential signal is analyzed, if the signal amplitude between the adjacent sampling points is too large, the signal can be determined to jump, and the direction of signal jump can be determined according to the positive and negative of the difference value of the adjacent sampling points. According to the jump of the signal and the jump direction, the direction and the size of the phase compensation required by the iteration can be determined by comparing with the historical count value k n-1 (t) accumulating to determine the count value k n (t), wherein k n (t) an initial value of 0 when n is 0.
According to an embodiment of the present disclosure, determining an initial unwind signal from an integrated signal and a fitted signal of the integrated signal comprises: performing polynomial fitting on the integrated signal based on the current iteration times to obtain a fitted signal; and calculating the difference between the integrated signal and the fitting signal to obtain an initial unwrapping signal.
According to an embodiment of the present disclosure, based on the current iteration number n, a polynomial fit may be performed on the integrated signal F (t) n times, resulting in a fitted signal P (t). The fitting signal P (t) may be used to fit the error trend term due to the n-th order integration, and the initial unwind signal may be determined according to equation (5)
According to the embodiment of the disclosure, the influence of the error trend term in the integral signal is eliminated through polynomial fitting, so that the distortion of the signal can be reduced, and the accuracy of the determined signal to be identified is improved.
According to an embodiment of the present disclosure, determining a current amplitude ratio from an initial unwind signal includes: performing Fourier transform on the initial unwrapped signal to obtain a transformed signal; calculating the sum of the amplitudes of the components of the frequency of the transformed signal in a first predetermined frequency range, resulting in a first frequency sum, wherein the first predetermined frequency range may be determined from the frequency of the interference signal and the nyquist frequency of the initial unwrapped signal; calculating the sum of the amplitudes of the components of the transformed signal whose frequencies lie in a first predetermined frequency range, obtaining a second frequency sum, which can be based on the nyquist frequency of the initial unwrapped signal; and calculating the ratio of the first frequency sum to the second frequency sum to obtain the current amplitude ratio.
According to embodiments of the present disclosure, the time domain signal may be transformed by fourier transformConversion to frequency domain signalsThe frequency of the interference signal in the first predetermined frequency range may be the Nyquist frequency f Nyquist Three fourths of (f), the first predetermined frequency range is (0.75 × f) Nyquist ，f Nyquist ). The second predetermined frequency range may be a signal of all frequency ranges in the frequency domain signal, i.e. the second predetermined frequency range is (0, f Nyquist ). The amplitudes of the components in the first predetermined frequency range may be summed to obtain a first frequency sum; the amplitudes of the components in the second predetermined frequency range are summed to obtain a second frequency sum. Current amplitude ratio R when current iteration number is n n Can be determined according to equation (6):
wherein when the current iteration number n is 0, R 0 Can be determined according to equation (7):
according to an embodiment of the present disclosure, the signal processing method further includes: according to the current iteration times, determining a historical amplitude ratio corresponding to the last iteration; and under the condition that the current amplitude ratio is larger than the historical amplitude ratio, determining that the initial unwinding signal meets the requirement.
According to the embodiments of the present disclosure, since the high frequency component duty ratio after conversion into the frequency domain signal should be smaller and smaller as the number of iterations increases, the current amplitude ratio should be smaller and smaller, i.e., R n ＜R n-1 When the distortion is small compared with the signal to be identified from iteration to initial unwrapping, the current amplitude ratio reaches the minimum value, and the iteration is continued at this time, the current amplitude ratio becomes larger, namely R n ＞R n-1 And (3) representing that the initial unwrapping signal meets the requirement, and taking the initial unwrapping signal of the n-1 th iteration as a signal to be identified.
According to an embodiment of the present disclosure, the signal processing method further includes: under the condition that the current amplitude ratio is smaller than the historical amplitude ratio, determining that the initial unwinding signal does not meet the requirement; and accumulating the current iteration times by 1 under the condition that the current amplitude ratio represents that the initial unwrapping signal does not meet the requirement, wherein the initial value of the current iteration times is 0.
According to embodiments of the present disclosure, the historical amplitude ratio at the time when the current amplitude ratio was less than the last iteration, i.e., R n ＜R n-1 In the event that the initial unwind signal is determined to be unsatisfactory. And adding 1 to the current iteration number, and continuing to perform iterative computation.
Fig. 3 schematically shows a flow chart of a signal processing method according to an embodiment of the present disclosure.
As shown in fig. 3, the flow includes operations S301 to S310.
In operation S301, an interference signal is orthogonally decomposed.
In operation S302, a winding phase signal is determined using an arctangent function.
In operation S303, the current iteration number n is initialized to 0.
In operation S304, an n-order difference signal of the winding phase signal is calculated.
In operation S305, a compensation signal is calculated.
In operation S306, an n-order integrated signal of the compensation signal is calculated.
In operation S307, an n-order fitting signal of the integrated signal is calculated.
In operation S308, it is determined whether the current amplitude ratio is greater than the historical amplitude ratio. If so, operation S309 is performed; if not, operation S310 is performed.
In operation S309, an initial unwind signal corresponding to the previous iteration number is taken as a result.
In operation S310, n=n+1 is taken, and operation S304 is performed.
According to an embodiment of the present disclosure, an interference signal is subjected to quadrature decomposition, and a winding phase signal is determined using an arctangent function and a result of the quadrature decomposition. Initializing the current iteration number n to 0, calculating an n-order differential signal of the winding phase signal, and calculating a compensation signal according to the differential signal. And calculating an n-order integral signal of the compensation signal, calculating an n-order fitting signal of the integral signal, and calculating the difference between the integral signal and the fitting signal to obtain an initial unwrapping signal. Calculating a current amplitude ratio according to the initial unwinding signal, judging whether the value of the current amplitude ratio is larger than the value of the historical amplitude ratio corresponding to the previous iteration number, and if so, taking the initial unwinding signal corresponding to the previous iteration number as an unwinding result; if not, updating the value of n to n+1, and continuing to execute the differential operation.
Fig. 4 (a) to (f) schematically show demodulation time domain waveform diagrams of the signal processing method according to the embodiment of the present disclosure, respectively.
As shown in fig. 4, fig. 4 (a) is a waveform diagram of a signal to be identified, the gray waveform in fig. 4 (b) is a waveform diagram of a winding signal of an arctangent algorithm, the range of values is (-pi, pi), the black waveform is a waveform diagram of an initial unwinding signal of 0 th order, and the trend of the black waveform is consistent with that of the signal to be identified. Graph (c) is a waveform diagram of the 1 st order initial unwrap signal with fewer high frequency spurs than the 0 th order initial unwrap signal. The graphs (d), (e) and (f) respectively show waveform diagrams of the initial unwinding signals of the 2-order, the 3-order and the 4-order, the results of three iterations are similar, but the initial unwinding result of the 2-order is closest to the signal to be identified through the size relation of the first two wave peaks.
Fig. 5 (a) to (d) schematically show demodulation spectrograms of the signal processing method according to the embodiment of the present disclosure, respectively.
As shown in fig. 5, black in fig. 5 (a) represents the spectrum of the signal to be identified, and gray represents the spectrum of the wrapped signal of the arctangent algorithm. The graphs (b) and (c) are the spectrums of the 0-order and 1-order initial unwind signals, respectively, the solid line in the graph (d) is the spectrum of the 2-order initial unwind signal, the dotted line is the spectrum of the 3-order initial unwind signal, and the dot-dashed line is the spectrum of the 4-order initial unwind signal. As can be seen from fig. (a), fig. (b) and fig. (c), the spectrum of the initial unwind signal which is not completely unwound contains significant high-frequency noise, the high-frequency noise gradually decreases with increasing iteration number, and as can be seen from fig. (d), the high-frequency noise of the 2-order initial unwind signal is relatively lowest, and the iteration number is continuously increased, so that the high-frequency noise is increased, and the 2-order initial unwind signal can be taken as the best result of the signal to be identified.
According to an embodiment of the present disclosure, the phase initial unwrapping effect of the signal processing method and the comparison method of the present disclosure is verified through simulation calculation, and the experimental result is shown in table 1, where r (ori, DUI) represents a correlation coefficient between an initial unwrapped signal and a signal to be identified, r (inter, PF) represents a correlation coefficient between an integrated signal and a polynomial fitting signal after each iteration, where the correlation coefficient may be a spearman coefficient, and the closer the value is to 1, the stronger the positive correlation is, and the higher the similarity degree between the two is.The current amplitude ratio after each iteration is represented, and the smaller the current amplitude ratio is, the weaker the high-frequency component is represented.
TABLE 1
According to embodiments of the present disclosure, the conventional method comparison algorithm may select the DUI algorithm. Judging whether the initial unwrapping signal meets the requirement or not through the correlation coefficient in the DUI algorithm, if the threshold value is set to epsilon & gt 0.99999999, outputting the result after the third iteration finally, wherein the signal processing method is superior to the DUI algorithm in performance because the signal processing method is used for outputting the result after the second iteration of the bit according to the frequency spectrum judgment and the correlation coefficient between the initial unwrapping signal and the signal to be identified, and the initial unwrapping signal after the second iteration is closer to the signal to be identified.
Fig. 6 schematically shows a block diagram of a signal processing apparatus according to an embodiment of the present disclosure.
As shown in fig. 6, the signal processing apparatus 600 of this embodiment includes a signal differentiating module 610, a signal integrating module 620, an unwrapped signal determining module 630, a ratio determining module 640, and a signal to be identified determining module 650.
The signal differentiating module 610 is configured to perform a differentiating operation on a wrapped phase signal of an interference signal based on a current iteration number, to obtain a differential signal of the wrapped phase signal, where the interference signal is obtained by heterodyning interference between a reference signal and a signal to be identified.
The signal integration module 620 is configured to integrate the compensation signal of the differential signal based on the current iteration number, to obtain an integrated signal of the compensation signal.
The unwrap signal determining module 630 is configured to determine an initial unwrap signal based on the integrated signal and a fit of the integrated signal.
The ratio determination module 640 is configured to determine a current amplitude ratio based on the initial unwind signal.
The signal to be identified determining module 650 is configured to determine a signal to be identified according to the initial unwrapped signal corresponding to the previous iteration number if the current amplitude ratio characterizes the initial unwrapped signal to meet the requirement.
According to an embodiment of the present disclosure, the ratio determination module 640 includes a fourier transform unit, a first summing unit, a second summing unit, and a ratio calculation unit.
The Fourier transform unit is used for carrying out Fourier transform on the initial unwrapped signal to obtain a transformed signal.
The first summing unit is configured to calculate a sum of amplitudes of components of the transformed signal having frequencies within a first predetermined frequency range, where the first predetermined frequency range may be determined based on the frequency of the interference signal and a nyquist frequency of the initial unwrapped signal.
The second summing unit is configured to calculate a sum of amplitudes of components of the transformed signal having frequencies within a first predetermined frequency range, and obtain a second frequency sum, where the second predetermined frequency range may be based on a nyquist frequency of the initial unwrapped signal.
The ratio calculating unit is used for calculating the ratio of the first frequency addition and the second frequency addition to obtain the current amplitude ratio.
According to an embodiment of the present disclosure, the unwind signal determination module 630 includes a polynomial fitting unit and a signal difference calculation unit.
The polynomial fitting unit is used for performing polynomial fitting on the integral signal based on the current iteration times to obtain a fitting signal.
The signal difference calculation unit is used for calculating the difference between the integrated signal and the fitting signal to obtain an initial unwrapping signal.
According to an embodiment of the present disclosure, the signal processing apparatus 600 further includes a historical amplitude ratio calculation module and a first result determination module.
The historical amplitude ratio calculation module is used for determining the historical amplitude ratio corresponding to the last iteration according to the current iteration times.
The first result determining module is used for determining that the initial unwind signal meets the requirement under the condition that the current amplitude ratio is larger than the historical amplitude ratio.
According to an embodiment of the present disclosure, the signal processing apparatus 600 further includes a second result determination module and a number accumulation module.
The second result determining module is used for determining that the initial unwind signal does not meet the requirement under the condition that the current amplitude ratio is smaller than the historical amplitude ratio.
The frequency accumulation module is used for accumulating the current iteration frequency by 1 under the condition that the current amplitude ratio represents that the initial unwrapping signal does not meet the requirement, wherein the initial value of the current iteration frequency is 0.
According to an embodiment of the present disclosure, the signal processing apparatus 600 further comprises a compensation signal determination module.
The compensation signal determining module is used for determining a compensation signal according to the differential signal and a count value, wherein the count value is determined according to a historical count value corresponding to the last iteration and the differential signal.
According to an embodiment of the present disclosure, the signal processing apparatus 600 further includes a quadrature component determination module and a wrapping phase signal determination module.
The quadrature component determination module is used for determining quadrature components of the interference signal by using a demodulation algorithm.
The wrapping phase signal determining module is used for determining a wrapping phase signal of the interference signal by utilizing an arctangent function based on the ratio of the orthogonal components.
According to an embodiment of the present disclosure, the compensation signal determining module includes a difference analyzing unit, a current count value determining unit, and a count value summing unit.
The difference analysis unit is used for analyzing the difference value of adjacent sampling points of the differential signal to obtain the jump direction and amplitude.
The current count value determining unit is used for determining the current count value according to the direction and the amplitude of the jump.
And the count value summation unit is used for summing the historical count value and the current count value to obtain the count value.
According to an embodiment of the present disclosure, the orthogonal component determining module further includes an orthogonal decomposition unit.
The orthogonal decomposition unit is used for carrying out orthogonal decomposition on the interference signal to obtain a first orthogonal component and a second orthogonal component, wherein the first orthogonal component is determined according to the interference signal and the sine of the local oscillator, and the second orthogonal component is determined according to the interference signal and the cosine of the local oscillator.
Any of the signal differencing module 610, the signal integrating module 620, the unwrapped signal determining module 630, the ratio determining module 640, and the signal to be identified determining module 650 may be implemented in one module or any of the modules may be split into multiple modules according to embodiments of the present disclosure. Alternatively, at least some of the functionality of one or more of the modules may be combined with at least some of the functionality of other modules and implemented in one module.
According to embodiments of the present disclosure, at least one of the signal differencing module 610, the signal integrating module 620, the unwrapped signal determining module 630, the ratio determining module 640, and the signal to be identified determining module 650 may be implemented at least in part as hardware circuitry, such as a Field Programmable Gate Array (FPGA), a Programmable Logic Array (PLA), a system-on-chip, a system-on-substrate, a system-on-package, an Application Specific Integrated Circuit (ASIC), or in hardware or firmware in any other reasonable manner of integrating or packaging the circuitry, or in any one of or a suitable combination of three of software, hardware, and firmware. Alternatively, at least one of the signal differencing module 610, the signal integrating module 620, the unwrapped signal determining module 630, the ratio determining module 640, and the signal to be identified determining module 650 may be implemented at least in part as a computer program module that, when executed, performs the corresponding functions.
Fig. 7 schematically shows a block diagram of an electronic device adapted for a signal processing method according to an embodiment of the disclosure.
As shown in fig. 7, an electronic device 700 according to an embodiment of the present disclosure includes a processor 701 that can perform various appropriate actions and processes according to a program stored in a Read Only Memory (ROM) 702 or a program loaded from a storage section 708 into a Random Access Memory (RAM) 703. The processor 701 may include, for example, a general purpose microprocessor (e.g., a CPU), an instruction set processor and/or an associated chipset and/or a special purpose microprocessor (e.g., an Application Specific Integrated Circuit (ASIC)), or the like. The processor 701 may also include on-board memory for caching purposes. The processor 701 may comprise a single processing unit or a plurality of processing units for performing different actions of the method flows according to embodiments of the disclosure.
In the RAM 703, various programs and data necessary for the operation of the electronic apparatus 700 are stored. The processor 701, the ROM 702, and the RAM 703 are connected to each other through a bus 704. The processor 701 performs various operations of the method flow according to the embodiments of the present disclosure by executing programs in the ROM 702 and/or the RAM 703. Note that the program may be stored in one or more memories other than the ROM 702 and the RAM 703. The processor 701 may also perform various operations of the method flow according to embodiments of the present disclosure by executing programs stored in the one or more memories.
According to an embodiment of the present disclosure, the electronic device 700 may further include an input/output (I/O) interface 705, the input/output (I/O) interface 705 also being connected to the bus 704. The electronic device 700 may also include one or more of the following components connected to the input/output I/O interface 705: an input section 706 including a keyboard, a mouse, and the like; an output portion 707 including a Cathode Ray Tube (CRT), a Liquid Crystal Display (LCD), and the like, a speaker, and the like; a storage section 708 including a hard disk or the like; and a communication section 709 including a network interface card such as a LAN card, a modem, or the like. The communication section 709 performs communication processing via a network such as the internet. The drive 710 is also connected to the I/O interface 705 as needed. A removable medium 711 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, or the like is mounted on the drive 710 as necessary, so that a computer program read therefrom is mounted into the storage section 708 as necessary.
The present disclosure also provides a computer-readable storage medium that may be embodied in the apparatus/device/system described in the above embodiments; or may exist alone without being assembled into the apparatus/device/system. The computer-readable storage medium carries one or more programs which, when executed, implement methods in accordance with embodiments of the present disclosure.
According to embodiments of the present disclosure, the computer-readable storage medium may be a non-volatile computer-readable storage medium, which may include, for example, but is not limited to: a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the context of this disclosure, a computer-readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. For example, according to embodiments of the present disclosure, the computer-readable storage medium may include ROM 702 and/or RAM 703 and/or one or more memories other than ROM 702 and RAM 703 described above.
Embodiments of the present disclosure also include a computer program product comprising a computer program containing program code for performing the methods shown in the flowcharts. The program code, when executed in a computer system, causes the computer system to perform the methods provided by embodiments of the present disclosure.
The above-described functions defined in the system/apparatus of the embodiments of the present disclosure are performed when the computer program is executed by the processor 701. The systems, apparatus, modules, units, etc. described above may be implemented by computer program modules according to embodiments of the disclosure.
In one embodiment, the computer program may be based on a tangible storage medium such as an optical storage device, a magnetic storage device, or the like. In another embodiment, the computer program may also be transmitted, distributed over a network medium in the form of signals, downloaded and installed via the communication section 709, and/or installed from the removable medium 711. The computer program may include program code that may be transmitted using any appropriate network medium, including but not limited to: wireless, wired, etc., or any suitable combination of the foregoing.
In such an embodiment, the computer program may be downloaded and installed from a network via the communication portion 709, and/or installed from the removable medium 711. The above-described functions defined in the system of the embodiments of the present disclosure are performed when the computer program is executed by the processor 701. The systems, devices, apparatus, modules, units, etc. described above may be implemented by computer program modules according to embodiments of the disclosure.
According to embodiments of the present disclosure, program code for performing computer programs provided by embodiments of the present disclosure may be written in any combination of one or more programming languages, and in particular, such computer programs may be implemented in high-level procedural and/or object-oriented programming languages, and/or assembly/machine languages. Programming languages include, but are not limited to, such as Java, c++, python, "C" or similar programming languages. The program code may execute entirely on the user's computing device, partly on the user's device, partly on a remote computing device, or entirely on the remote computing device or server. In the case of remote computing devices, the remote computing device may be connected to the user computing device through any kind of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or may be connected to an external computing device (e.g., connected via the Internet using an Internet service provider).
The flowcharts and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams or flowchart illustration, and combinations of blocks in the block diagrams or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
Those skilled in the art will appreciate that the features recited in the various embodiments of the disclosure and/or in the claims may be provided in a variety of combinations and/or combinations, even if such combinations or combinations are not explicitly recited in the disclosure. In particular, the features recited in the various embodiments of the present disclosure and/or the claims may be variously combined and/or combined without departing from the spirit and teachings of the present disclosure. All such combinations and/or combinations fall within the scope of the present disclosure.
While the foregoing is directed to embodiments of the present disclosure, other and further details of the invention may be had by the present application, it is to be understood that the foregoing description is merely exemplary of the present disclosure and that no limitations are intended to the scope of the disclosure, except insofar as modifications, equivalents, improvements or modifications may be made without departing from the spirit and principles of the present disclosure.
Claims (10)
1. A signal processing method, comprising:
performing differential operation on a winding phase signal of an interference signal based on the current iteration times to obtain a differential signal of the winding phase signal, wherein the interference signal is obtained by heterodyne interference of a reference signal and a signal to be identified;
Based on the current iteration times, carrying out integral operation on the compensation signal of the differential signal to obtain an integral signal of the compensation signal;
determining an initial unwrapping signal based on the integrated signal and a fit of the integrated signal;
determining a current amplitude ratio based on the initial unwind signal;
and under the condition that the current amplitude ratio represents that the initial unwinding signal meets the requirement, determining the signal to be identified according to the initial unwinding signal corresponding to the previous iteration number.
2. The method of claim 1, wherein said determining a current amplitude ratio from said initial unwind signal comprises:
performing Fourier transform on the initial unwrapping signal to obtain a transformed signal;
calculating the sum of the amplitudes of the components of the frequency of the transformation signal in a first preset frequency range to obtain a first frequency sum, wherein the first preset frequency range can be determined according to the frequency of the interference signal and the Nyquist frequency of the initial unwrapping signal;
calculating the sum of the amplitudes of the components of the frequency of the transformation signal in a first preset frequency range to obtain a second frequency sum, wherein the second preset frequency range can be according to the Nyquist frequency of the initial unwrapped signal;
And calculating the ratio of the first frequency sum to the second frequency sum to obtain the current amplitude ratio.
3. The method of claim 1, wherein the determining an initial unwind signal from the integrated signal and a fitted signal of the integrated signal comprises:
performing polynomial fitting on the integral signal based on the current iteration times to obtain a fitting signal;
and calculating the difference between the integral signal and the fitting signal to obtain an initial unwrapping signal.
4. The method of claim 1, further comprising:
according to the current iteration times, determining a historical amplitude ratio corresponding to the last iteration;
and under the condition that the current amplitude ratio is larger than the historical amplitude ratio, determining that the initial unwind signal meets the requirement.
5. The method of claim 4, further comprising:
determining that the initial unwind signal does not meet a requirement if the current amplitude ratio is less than the historical amplitude ratio;
and accumulating the current iteration times by 1 under the condition that the current amplitude ratio represents that the initial unwrapping signal does not meet the requirement, wherein the initial value of the current iteration times is 0.
6. The method of claim 1, further comprising:
and determining the compensation signal according to the differential signal and a count value, wherein the count value is determined according to a historical count value corresponding to the last iteration and the differential signal.
7. The method of claim 1, further comprising:
determining a quadrature component of the interference signal using a demodulation algorithm;
a wrap-around phase signal of the interference signal is determined using an arctangent function based on the ratio of the quadrature components.
8. The method of claim 6, wherein said determining said compensation signal from said differential signal and a count value comprises:
analyzing the difference value of adjacent sampling points of the differential signals to obtain the jump direction and amplitude;
determining the count value according to the direction and the amplitude of the jump;
and summing the historical count value and the current count value to obtain the count value.
9. The method of claim 7, wherein said determining the quadrature component of the interference signal using a demodulation algorithm comprises:
and carrying out orthogonal decomposition on the interference signal to obtain a first orthogonal component and a second orthogonal component, wherein the first orthogonal component is determined according to the interference signal and the sine of the local oscillator, and the second orthogonal component is determined according to the interference signal and the cosine of the local oscillator.
10. A phase unwrapping apparatus comprising:
the signal difference module is used for carrying out difference operation on a winding phase signal of an interference signal based on the current iteration times to obtain a difference signal of the winding phase signal, wherein the interference signal is obtained by heterodyne interference of a reference signal and a signal to be identified;
the signal integration module is used for carrying out integration operation on the compensation signal of the differential signal based on the current iteration times to obtain an integrated signal of the compensation signal;
the unwrapping signal determining module is used for determining an initial unwrapping signal according to the integrated signal and a fitting signal of the integrated signal;
the ratio determining module is used for determining a current amplitude ratio based on the initial unwinding signal;
the signal to be identified determining module is configured to determine the signal to be identified according to the initial unwrapping signal corresponding to the previous iteration number when the current amplitude ratio represents that the initial unwrapping signal meets a requirement.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311380114.9A CN117433622A (en) | 2023-10-24 | 2023-10-24 | Signal processing method, device, equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311380114.9A CN117433622A (en) | 2023-10-24 | 2023-10-24 | Signal processing method, device, equipment and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117433622A true CN117433622A (en) | 2024-01-23 |
Family
ID=89554660
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311380114.9A Pending CN117433622A (en) | 2023-10-24 | 2023-10-24 | Signal processing method, device, equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117433622A (en) |
-
2023
- 2023-10-24 CN CN202311380114.9A patent/CN117433622A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6392289B2 (en) | Scaling fixed-point fast Fourier transforms in radar and sonar applications | |
| US8767193B2 (en) | Doppler tracking in presence of vehicle velocity uncertainty | |
| KR102458095B1 (en) | Phase correction method and device | |
| WO2017132670A1 (en) | Fixed-point high dynamic range fast fourier transforms | |
| JP5110531B2 (en) | FM-CW radar equipment | |
| CN117433622A (en) | Signal processing method, device, equipment and storage medium | |
| KR100817692B1 (en) | A method for estimating phase angle of time series data by discrete Fourier transform | |
| CN116701840A (en) | Method and system for optimizing and calculating cepstrum of mechanical vibration signal | |
| US20150134712A1 (en) | Signal processing device, signal processing method, and program | |
| US11181605B2 (en) | Method and device for estimating an angle of arrival of an incident radio signal | |
| JP2004150825A (en) | Spectral analysis apparatus and method | |
| George et al. | A novel and efficient hardware accelerator architecture for signal normalization | |
| CN115604066B (en) | Low-cost narrowband interference frequency estimation method based on discrete Fourier transform | |
| Sedlacek et al. | Iterative design of DSP cosine windows by placing spectrum zeros, and some new classes of DSP cosine windows | |
| Makdissi et al. | A signal approach analysis of the Ramsey pattern in cesium beam frequency standards | |
| WO2022107247A1 (en) | Radar signal processing device, radar device, and radar signal processing method | |
| Regimanu et al. | Dither filtering of real RLG signal using wavelet transforms | |
| Henry | An ultra-precise Fast Fourier Transform | |
| Zhu et al. | FouriDown: Factoring Down-Sampling into Shuffling and Superposing | |
| CN113129921B (en) | Method and apparatus for detecting frequency of fundamental tone in speech signal | |
| Glinchenko | Expanding the Dynamic Range of Digital Spectral Measurements | |
| Li et al. | On uncertainty principle of the local polynomial Fourier transform | |
| Belega et al. | Generalized multi-point interpolated DTFT frequency and damping factor estimators of real-valued damped sinusoids | |
| CN116820402A (en) | Variance calculation method, device, circuit and equipment for true random number generator | |
| CN113533238A (en) | Method and system for detecting sulfur hexafluoride decomposition gas based on absorption spectrum |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination |