CN117201250B - Phase generation carrier demodulation method and device, electronic equipment and storage medium - Google Patents
Phase generation carrier demodulation method and device, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN117201250B CN117201250B CN202311466287.2A CN202311466287A CN117201250B CN 117201250 B CN117201250 B CN 117201250B CN 202311466287 A CN202311466287 A CN 202311466287A CN 117201250 B CN117201250 B CN 117201250B
- Authority
- CN
- China
- Prior art keywords
- signal
- ellipse
- orthogonal
- component
- components
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 41
- 238000012937 correction Methods 0.000 claims abstract description 24
- 238000012545 processing Methods 0.000 claims abstract description 15
- 230000006870 function Effects 0.000 claims description 31
- 239000013598 vector Substances 0.000 claims description 27
- 238000004590 computer program Methods 0.000 claims description 14
- 238000001914 filtration Methods 0.000 claims description 14
- 238000013507 mapping Methods 0.000 claims description 12
- 230000008859 change Effects 0.000 claims description 8
- 230000036962 time dependent Effects 0.000 claims description 5
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000004891 communication Methods 0.000 abstract description 2
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000003071 parasitic effect Effects 0.000 description 6
- 238000005259 measurement Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 4
- 230000003287 optical effect Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000002093 peripheral effect Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02D—CLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
- Y02D30/00—Reducing energy consumption in communication networks
- Y02D30/70—Reducing energy consumption in communication networks in wireless communication networks
Landscapes
- Digital Transmission Methods That Use Modulated Carrier Waves (AREA)
Abstract
The embodiment of the application provides a phase generation carrier demodulation method, a phase generation carrier demodulation device, electronic equipment and a storage medium, and relates to the technical field of communication. The method comprises the following steps: determining an orthogonal component corresponding to an interference signal, wherein the interference signal comprises a high-frequency carrier signal and a signal to be detected carried on the high-frequency carrier signal; performing ellipse fitting on the orthogonal components to obtain ellipse parameters; nonlinear correction is carried out on the orthogonal components based on the ellipse parameters, so that mutually orthogonal signal components to be processed are obtained; and performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be detected. The embodiment of the application solves the technical problem that the reliability of the existing phase generation carrier demodulation technology is poor, and improves the demodulation accuracy.
Description
Technical Field
The present invention relates to the field of communications technologies, and in particular, to a method and apparatus for demodulating a phase-generating carrier, an electronic device, and a storage medium.
Background
In a modem system, an accurate carrier frequency and phase are required during demodulation. However, when the signal undergoes frequency and phase changes during transmission, such as doppler effect, carrier drift, phase jitter, etc., conventional demodulation techniques are susceptible to interference, resulting in inaccurate demodulation results or lost data.
In contrast, the phase generating carrier demodulation technology transfers the phase information to be detected to the sideband of the high-frequency carrier and the higher harmonic thereof by introducing the high-frequency carrier signal, has the advantages of low-frequency interference resistance, high sensitivity, large dynamic range, phase fading resistance and the like, and is widely applied to the phase demodulation of an interference type optical fiber sensor, a self-mixing interferometer and a sine phase modulation interferometer. At present, the phase generation carrier demodulation technology mainly comprises a differential cross multiplication algorithm (pgc-dcm) and an arctangent algorithm (pgc-arctan), and the problem that the traditional demodulation technology is easy to interfere, so that a demodulation result is inaccurate or data is lost is solved. However, the operation result of the differential cross multiplication algorithm has a linear relation with the phase to be measured, but the demodulation result is affected by the fluctuation of the laser light intensity, the phase delay and the modulation depth. The arctangent algorithm removes the effect of light intensity fluctuations, but is still affected by phase delay and modulation depth, and introduces nonlinear errors. Therefore, the existing phase generation carrier demodulation technique is poor in reliability.
Disclosure of Invention
The embodiment of the application provides a phase generation carrier demodulation method, a device, electronic equipment and a storage medium, which are used for solving the technical problem of poor reliability of the existing phase generation carrier demodulation technology and improving the demodulation accuracy.
According to an aspect of the embodiments of the present application, there is provided a phase generation carrier demodulation method, including:
determining an orthogonal component corresponding to an interference signal, wherein the interference signal comprises a high-frequency carrier signal and a signal to be detected carried on the high-frequency carrier signal;
performing ellipse fitting on the orthogonal components to obtain ellipse parameters;
nonlinear correction is carried out on the orthogonal components based on the ellipse parameters, so that mutually orthogonal signal components to be processed are obtained;
and performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be detected.
In one possible implementation, the determining the orthogonal component corresponding to the interference signal includes:
and mixing the interference signal with a single frequency multiplication carrier signal and a double frequency multiplication carrier signal respectively, and then carrying out low-pass filtering to obtain a pair of orthogonal components containing the signal to be detected, wherein the orthogonal components contain sine components and cosine components which are mutually orthogonal.
In one possible implementation, the sinusoidal component includes a first amplitude and a corresponding first bias that characterize the sinusoidal component over time; the cosine component includes a second amplitude and a corresponding second bias characterizing a time-dependent change of the cosine component;
the nonlinear correction is performed on the orthogonal components based on the ellipse parameters to obtain mutually orthogonal signal components to be processed, including:
determining mapping relations between each ellipse parameter and the first amplitude, the first bias, the second amplitude and the second bias;
and based on an ellipse parameter obtained by ellipse fitting and the mapping relation, obtaining the first amplitude, the first offset, the second amplitude and the second offset, so that the sine component and the cosine component are corrected to obtain a signal component to be processed.
In one possible implementation manner, the performing ellipse fitting on the orthogonal components to obtain ellipse parameters includes:
and carrying out ellipse fitting on the orthogonal components at any time, and carrying out nonlinear solving by utilizing a pre-established ellipse fitting model to output ellipse parameters corresponding to the time when an observation function fitted in the ellipse fitting model is minimized, wherein the observation function is used for indicating deviation between the orthogonal components and a fitted ellipse equation.
In one possible implementation manner, the non-linear solving by using the pre-established ellipse fitting model includes:
and for solving any ellipse parameter, iteratively predicting the current optimal state vector for representing the ellipse parameter by using a prediction model in the ellipse fitting model to obtain the next optimal state vector, and taking the current optimal state vector as the output of the ellipse fitting model until the difference value of the optimal state vectors obtained in two adjacent times reaches a preset stable condition.
In one possible implementation manner, the performing nonlinear correction on the orthogonal component based on the ellipse parameter to obtain mutually orthogonal signal components to be processed further includes:
and inputting the ellipse parameters into the observation function to correct the signal components which are not mutually orthogonal, so as to obtain mutually orthogonal signal components to be processed.
In one possible implementation manner, the performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be measured includes:
dividing the mutually orthogonal signal components to be processed to obtain tangent signals;
and carrying out high-pass filtering after carrying out arc tangent operation on the tangent signal so as to demodulate a signal to be detected and obtain phase information in the signal to be detected.
According to another aspect of the embodiments of the present application, there is provided a phase generating carrier demodulation apparatus, including:
the system comprises a quadrature component determining module, a quadrature component determining module and a quadrature component determining module, wherein the quadrature component determining module is used for determining a quadrature component corresponding to an interference signal, and the interference signal comprises a high-frequency carrier signal and a signal to be detected carried on the high-frequency carrier signal;
the ellipse fitting module is used for performing ellipse fitting on the orthogonal components to obtain ellipse parameters;
the nonlinear correction module is used for carrying out nonlinear correction on the orthogonal components based on the elliptic parameters to obtain mutually orthogonal signal components to be processed;
and the arctangent demodulation module is used for carrying out arctangent demodulation processing on the signal component to be processed so as to obtain phase information in the signal to be detected.
According to another aspect of the embodiments of the present application, there is provided an electronic device including a memory, a processor and a computer program stored on the memory, the processor executing the computer program to implement the steps of the phase-generating carrier demodulation method described in the above embodiments.
According to still another aspect of the embodiments of the present application, there is provided a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the phase-generating carrier demodulation method described in the above embodiments.
The beneficial effects that technical scheme that this application embodiment provided brought are:
according to the phase generation carrier demodulation method, the orthogonal component corresponding to the interference signal is determined, the interference signal comprises the high-frequency carrier signal and the signal to be detected borne on the high-frequency carrier signal, ellipse fitting is conducted on the orthogonal component to obtain ellipse parameters, nonlinear correction is conducted on the orthogonal component based on the ellipse parameters to obtain mutually orthogonal signal components to be processed, and therefore arc tangent demodulation processing is conducted on the signal components to be processed to obtain phase information in the signal to be detected, influences such as phase delay and modulation depth fluctuation caused by interference noise and parasitic amplitude are eliminated, nonlinear errors are eliminated, the technical problem that reliability of an existing phase generation carrier demodulation technology is poor is solved, and demodulation accuracy is improved.
Drawings
Fig. 1 is a flow chart of a phase generating carrier demodulation method according to an embodiment of the present application;
fig. 2 is a schematic structural diagram of a phase generating carrier demodulation device according to an embodiment of the present application;
fig. 3 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
The following description of the embodiments of the present application will be made clearly and fully with reference to the accompanying drawings, in which it is evident that the embodiments described are only some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the present disclosure, are within the scope of the present disclosure.
Example 1
Referring to fig. 1, a flowchart of a phase generating carrier demodulation method according to an embodiment of the present application is provided, where the method includes steps S101 to S104.
S101, determining orthogonal components corresponding to interference signals, wherein the interference signals comprise high-frequency carrier signals and signals to be detected loaded on the high-frequency carrier signals.
In some embodiments, the determining the quadrature component corresponding to the interference signal includes:
and mixing the interference signal with a single frequency multiplication carrier signal and a double frequency multiplication carrier signal respectively, and then carrying out low-pass filtering to obtain a pair of orthogonal components containing the signal to be detected, wherein the orthogonal components contain sine components and cosine components which are mutually orthogonal.
It should be noted that, by introducing a high-frequency carrier signal in a high frequency band outside the bandwidth of the signal to be measured, the signal to be measured is modulated onto the side band of the high-frequency carrier signal to form an interference signal. In this embodiment, the interference signal is mixed with the single frequency multiplication carrier signal and the double frequency multiplication carrier signal respectively by using the mixer to output two mixed signals, and then the mixed signals are subjected to low-pass filtering by the low-pass filter to obtain a pair of mutually orthogonal sine components and cosine components.
Specifically, in the PGC system of the external modulation system, the interference signal outputted by the interferometer is expressed asWherein S (t) is an interference signal, A and B are respectively the DC bias and AC amplitude of the interference signal, C is the phase modulation depth, omega c For the frequency of the high-frequency modulation signal, t is the sampling time point, delta theta is the carrier phase delay caused by the factors of optical path propagation, circuit transmission, digital-to-analog conversion and the like, and +.>For the phase to be measured (i.e. the phase information of the signal to be measured) of the interferometer->Where d (t) is the displacement of the interferometer to be measured and λ is the wavelength of the laser. Further, the sinusoidal components are expressed asWherein P (t) is a sinusoidal component, x LPF [·]For low-pass filtering operation, J 1 (C) And J 2 (C) For the first and second Bessel functions of the first class, the value of C ranges from 1.5 to 3.5rad and the corresponding sum value is between 0 and 1 in the research process. And cosine component is expressed asWhere Q (t) is a sinusoidal component.
S102, carrying out ellipse fitting on the orthogonal components to obtain ellipse parameters.
In the method, the signals of the high-frequency carrier signals subjected to low-pass filtering are discrete point signals, the Euclidean-based ellipse fitting method is used for accurately estimating geometric parameters (namely ellipse parameters) of ellipses, and the influence of phase delay and modulation depth can be eliminated simultaneously by carrying out ellipse fitting on orthogonal signals, so that nonlinear errors are reduced, and the accuracy of demodulation is improved.
S103, carrying out nonlinear correction on the orthogonal components based on the ellipse parameters to obtain mutually orthogonal signal components to be processed.
S104, performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be detected.
In the method, elliptical parameters obtained by elliptical fitting are input into a nonlinear model, and non-orthogonal signals are corrected, so that the problem that the maximum signals of the two channels are non-orthogonal is solved. Therefore, through ellipse fitting and nonlinear correction, accurate control of an input signal of an arctangent algorithm is realized, phase information to be detected for eliminating influences such as carrier phase delay and modulation depth fluctuation can be obtained, the phase of a signal can be quickly and adaptively tracked and recovered under nonlinear errors caused by interference noise and parasitic amplitude modulation in a phase generation carrier demodulation technology, the carrier phase in the signal can be accurately tracked and extracted, the reliability and robustness are improved, and the demodulation process is more stable and accurate.
According to the phase generation carrier demodulation method provided by the embodiment, through determining the quadrature component corresponding to the interference signal, the interference signal comprises a high-frequency carrier signal and a signal to be detected borne on the high-frequency carrier signal, elliptical fitting is carried out on the quadrature component to obtain elliptical parameters, nonlinear correction is carried out on the quadrature component based on the elliptical parameters to obtain mutually orthogonal signal components to be processed, and thus arctangent demodulation processing is carried out on the signal components to be processed to obtain phase information in the signal to be detected, the effects of phase delay, modulation depth fluctuation and the like caused by interference noise and parasitic amplitude are eliminated, nonlinear errors are eliminated, the technical problem that the reliability of an existing phase generation carrier demodulation technology is poor is solved, and demodulation accuracy is improved.
In some embodiments, the performing ellipse fitting on the orthogonal components to obtain ellipse parameters includes:
and carrying out ellipse fitting on the orthogonal components at any time, and carrying out nonlinear solving by utilizing a pre-established ellipse fitting model to output ellipse parameters corresponding to the time when an observation function fitted in the ellipse fitting model is minimized, wherein the observation function is used for indicating deviation between the orthogonal components and a fitted ellipse equation.
For example, based on ellipse fitting of geometric distance, 5 parameters can be uniquely expressed for ellipse on plane, namely center (X, Y), principal axes a, b, angle alpha, ellipse equation expressed as aX 2 +fXY+bY 2 +cx+dy+e=0, wherein a, b, c, d, e, f are all elliptical parameters. From the acquired plurality of measurement points Pi (Xi, yi), the fitted objective function (i.e., observation function) is expressed asn is the number of samples, i.e. the objective function F is used to quantify the deviation between all the measurement points Pi and the elliptic equation, thus solving a set of elliptic parameters a, b, c, d, e such that the objective function F reaches a minimum value so that all the measurement points are as far as possible on the ellipse. For this purposeTo minimize the objective function F, the first-order differentiation of the objective function F with respect to the ellipse parameters a, b, c, d, e is made 0, i.e.:to obtain corresponding ellipse parameters.
In this embodiment, for ellipse fitting, the data points corresponding to the orthogonal components fall on the ellipse when the orthogonal components satisfy the above-mentioned ellipse equation, and since the two orthogonal components are orthogonal to each other in this application, that is, XY is 0 in the ellipse equation, that is, the ellipse parameter f is not discussed, this application only describes the ellipse parameters a, b, c, d, e, so that the expression is aP (t) 2 +bQ(t) 2 +cp (t) +dq (t) +e=0. And then, carrying out nonlinear solution by using an ellipse fitting model, wherein Gaussian-Newton iteration is carried out on the ellipse fitting model to improve the model accuracy. Specifically, taking the discretized orthogonal component as an input, the observation function fitted in the ellipse fitting model (i.e., the objective function described above) is expressed as:
wherein Ln is an observation function of ellipse parameters, P n As discretized sinusoidal components, Q n The state vector x is the discretized cosine component n =[a n b n c n d n e n ],a n ,b n ,c n ,d n ,e n Is the discretized ellipse parameter.
Further, the observation function is used for indicating deviation between the orthogonal component and the fitted elliptic equation, namely, the observation function characterizes fitting condition for one orthogonal component, so that a group of elliptic parameters a, b, c, d and e capable of minimizing the observation function are solved, and data points corresponding to the orthogonal component fall on the ellipse. Therefore, in the embodiment, by performing ellipse fitting on discrete orthogonal components, the ellipse parameters are accurately estimated, and the demodulation accuracy is improved.
In some embodiments, the non-linear solution using a pre-established ellipse fitting model includes:
and for solving any ellipse parameter, iteratively predicting the current optimal state vector for representing the ellipse parameter by using a prediction model in the ellipse fitting model to obtain the next optimal state vector, and taking the current optimal state vector as the output of the ellipse fitting model until the difference value of the optimal state vectors obtained in two adjacent times reaches a preset stable condition.
In this embodiment, the orthogonal components are ellipse fitted based on a gaussian-newton iterative algorithm, that is, the ellipse fitting model is solved based on the gaussian-newton iterative algorithm, so as to solve the ellipse parameters corresponding to the minimization of the observation function. Specifically, a predictive model based on an ellipse fitting model of Gaussian-Newton iteration is used to iteratively calculate a state vector x for the nth time n =[a n b n c n d n e n ]Is the optimal state vector of (a)Then is +.>Predict the a priori state vector n+1th time +.>And optimal state vector->And so on. Wherein, illustratively, during the iterative recursion from the nth to the n+1th, if the optimal state vector +.>And (4) optimal state vector->The absolute value of the difference between them is atWhen the preset difference value is within the preset difference value, the difference value is considered to reach the preset stable condition, and the optimal state vector of the nth time is +.>Can be regarded as a stable matrix so that the optimal state vector +.>The corresponding elliptic function serves as an elliptic parameter for minimizing the observation function and serves as an output of the elliptic fitting model for nonlinear solution. Therefore, in the embodiment, through carrying out Gaussian-Newton iteration on the elliptic fitting model, the accuracy of the model can be effectively improved, and the stability and accuracy of the demodulation process are further improved.
For example, for a nonlinear optimization method based on Gaussian-Newton hair, the general nonlinear least squares fitting problem is min||X-F (a) || 2 Where a is a parameter vector, X is a set of measurements, and F represents a nonlinear function with respect to a; a is equal to Rqa/in, rqa/in is the inner long radius of a nonlinear curve, and R≡Rq is an optimization parameter; x ε RqX/in, rqX/in is the inside long radius of the measured value, R++qX ε Rp is a known vector. The fitting problem can then be understood as assuming that there are q parameter vectors a associated with p (> q) measurements X, minimizing the error of both. Thus, the predictive model is denoted as a k+1 =a k +λAJF T (X-F(a k )),a k 、a k+1 For any value, JF is a Jacobian matrix, lambda is a step size, A is a scaling factor, and Deltaa is a difference value. When a= (JF) T JF) -1 When the Gaussian-Newton iteration is carried out, the A is brought into the prediction model, and a prediction equation set a is obtained k+1 =a k +λΔa and jfΔa=x-F (a k ) At this timei=1, 2,..p, j=1, 2,..q. Therefore, the above-described prediction equation set is solved by singular values, and Δa is solved for iteration.
In some embodiments, the sinusoidal component includes a first amplitude and a corresponding first bias used to characterize the sinusoidal component over time; the cosine component includes a second amplitude and a corresponding second bias characterizing a time-dependent change of the cosine component;
the nonlinear correction is performed on the orthogonal components based on the ellipse parameters to obtain mutually orthogonal signal components to be processed, including:
determining mapping relations between each ellipse parameter and the first amplitude, the first bias, the second amplitude and the second bias;
and based on an ellipse parameter obtained by ellipse fitting and the mapping relation, obtaining the first amplitude, the first offset, the second amplitude and the second offset, so that the sine component and the cosine component are corrected to obtain a signal component to be processed.
In this embodiment, the carrier phase delay and the phase modulation depth are dynamically changed due to the effects of carrier phase delay, modulation depth fluctuation and the like caused by low signal-to-noise ratio and high spurious amplitude, so that the quadrature component obtained by demodulating the phase generating carrier contains amplitude fluctuation and bias which change with time. At this time, the expression of the magnitudes and bias of the sine component and the cosine component is: wherein X (t) is a first amplitude, X 0 (t) is a first bias, Y (t) is a second amplitude, Y 0 (t) is a second bias. In the existing arctangent algorithm, the quadrature component is directly divided and arctangent operated, so that the phase information in the signal to be detected obtained by demodulation contains nonlinear errors caused by interference noise and parasitic amplitude modulation, and therefore, the embodiment measures and corrects the amplitude and bias of the quadrature component to correct the quadrature component, thereby eliminating the nonlinear errors of the demodulated signal to be detected and improving the demodulation accuracy.
Specifically, an ellipse is determinedMapping relation between parameters a, b, c, d, e and the first amplitude, the first bias, the second amplitude and the second bias, namelyThen, elliptical fitting is carried out according to the observation model Ln to obtain an optimal state vector +.>The amplitude and bias of the discrete orthogonal component are obtained by combining the mapping relation and expressed as +.>Wherein (1)>And (5) obtaining ellipse parameters for iterative calculation. Substituting the magnitudes and offsets into expressions of the magnitudes and offsets of the sine component and the cosine component to correct the sine component and the cosine component.
In some embodiments, the nonlinear correction is performed on the orthogonal component based on the ellipse parameter to obtain mutually orthogonal signal components to be processed, and the method further includes:
and inputting the ellipse parameters into the observation function to correct the signal components which are not mutually orthogonal, so as to obtain mutually orthogonal signal components to be processed.
In this embodiment, since the observation model Ln is a nonlinear model, the ellipse parameters obtained by ellipse fitting are brought into the nonlinear model, so that each time the orthogonal component is sampled and detected, the signal components detected as not being orthogonal to each other are corrected, and the situation that the two-channel maximum signals for the sine component and the cosine component are not orthogonal is avoided.
In some embodiments, the performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be measured includes:
dividing the mutually orthogonal signal components to be processed to obtain tangent signals;
and carrying out high-pass filtering after carrying out arc tangent operation on the tangent signal so as to demodulate a signal to be detected and obtain phase information in the signal to be detected.
In the present embodiment, the signal components to be processed which are orthogonal to each other are subjected to the inverse tangent operation after division operation based on the correction of the amplitude and offset of the orthogonal components, that is
Then, the demodulation of the signal to be tested is realized through high-pass filtering.
Therefore, the embodiment of the application eliminates the influences of carrier phase delay, modulation depth fluctuation and the like of the signal input into the arc tangent demodulation processing through elliptical fitting and nonlinear correction, namely, the influence of nonlinear error is not required to be considered, the accurate control of the input signal is realized, the phase of the signal can be quickly and adaptively tracked and recovered under the nonlinear error caused by interference noise and parasitic amplitude modulation, the carrier phase in the signal is accurately tracked and extracted, and the stability and reliability of the demodulation process are greatly improved.
Example two
Referring to fig. 2, a schematic structural diagram of a phase generating carrier demodulation apparatus according to an embodiment of the present application, where the apparatus 200 includes:
a quadrature component determining module 201, configured to determine a quadrature component corresponding to an interference signal, where the interference signal includes a high-frequency carrier signal and a signal to be detected carried on the high-frequency carrier signal;
an ellipse fitting module 202, configured to perform ellipse fitting on the orthogonal components to obtain ellipse parameters;
a nonlinear correction module 203, configured to perform nonlinear correction on the orthogonal component based on the ellipse parameter, so as to obtain mutually orthogonal signal components to be processed;
and the arctangent demodulation module 204 is configured to perform arctangent demodulation processing on the signal component to be processed, so as to obtain phase information in the signal to be detected.
In some embodiments, the orthogonal component determination module 201 includes:
the frequency mixing low-pass filtering unit is used for carrying out low-pass filtering after mixing the interference signal with the single frequency multiplication carrier signal and the double frequency multiplication carrier signal respectively to obtain a pair of orthogonal components containing the signal to be detected, wherein the orthogonal components contain sine components and cosine components which are mutually orthogonal.
In some embodiments, the sinusoidal component includes a first amplitude and a corresponding first bias used to characterize the sinusoidal component over time; the cosine component includes a second amplitude and a corresponding second bias characterizing a time-dependent change of the cosine component;
the nonlinear correction module 203 includes:
the mapping relation acquisition unit is used for determining the mapping relation between each ellipse parameter and the first amplitude, the first bias, the second amplitude and the second bias;
and the component correction unit is used for obtaining the first amplitude, the first bias, the second amplitude and the second bias based on the elliptic parameter obtained by elliptic fitting and the mapping relation, so as to correct the sine component and the cosine component to obtain a signal component to be processed.
In some embodiments, the ellipse fitting module 202 includes:
and the nonlinear solving unit is used for carrying out ellipse fitting on the orthogonal components at any time and carrying out nonlinear solving by utilizing a pre-established ellipse fitting model so as to output ellipse parameters corresponding to the time when an observation function fitted in the ellipse fitting model is minimized, wherein the observation function is used for indicating the deviation between the orthogonal components and a fitted ellipse equation.
In some embodiments, the nonlinear solving unit includes:
and the iteration unit is used for solving any ellipse parameter, predicting the current optimal state vector used for representing the ellipse parameter by using the prediction model in the ellipse fitting model in an iteration manner so as to obtain the next optimal state vector, and taking the current optimal state vector as the output of the ellipse fitting model until the difference value of the optimal state vectors obtained in two adjacent times reaches a preset stable condition.
In some embodiments, the nonlinear correction module 203 further comprises:
and the non-orthogonal correction unit is used for inputting the ellipse parameters into the observation function to correct the non-orthogonal signal components so as to obtain the mutually orthogonal signal components to be processed.
In some embodiments, the arctangent demodulation module 204 includes:
the division operation unit is used for carrying out division operation on the mutually orthogonal signal components to be processed to obtain tangent signals;
and the demodulation processing unit is used for carrying out high-pass filtering after carrying out arctangent operation on the tangent signal so as to demodulate the signal to be detected and obtain the phase information in the signal to be detected.
The apparatus of the embodiments of the present application may perform the method provided by the embodiments of the present application, and implementation principles of the method are similar, and actions performed by each module in the apparatus of each embodiment of the present application correspond to steps in the method of each embodiment of the present application, and detailed functional descriptions of each module of the apparatus may be referred to in the corresponding method shown in the foregoing, which is not repeated herein.
Example III
The embodiment of the application provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory, wherein the processor executes the computer program to realize the steps of a phase generation carrier demodulation method, and compared with the related art, the method can realize the steps of: the method and the device realize the elimination of the influences of phase delay, modulation depth fluctuation and the like caused by interference noise and parasitic amplitude, thereby eliminating nonlinear errors, solving the technical problem of poor reliability of the existing phase generation carrier demodulation technology and improving the demodulation accuracy.
In an alternative embodiment, an electronic device is provided, as shown in fig. 3, the electronic device 300 shown in fig. 3 includes: a processor 301 and a memory 303. Wherein the processor 301 is coupled to the memory 303, such as via a bus 302. Optionally, the electronic device 300 may further comprise a transceiver 304, the transceiver 304 may be used for data interaction between the electronic device and other electronic devices, such as transmission of data and/or reception of data, etc. It should be noted that, in practical applications, the transceiver 304 is not limited to one, and the structure of the electronic device 300 is not limited to the embodiment of the present application.
The processor 301 may be a CPU (Central Processing Unit ), general purpose processor, DSP (Digital Signal Processor, data signal processor), ASIC (Application Specific Integrated Circuit ), FPGA (Field Programmable Gate Array, field programmable gate array) or other programmable logic device, transistor logic device, hardware components, or any combination thereof. Which may implement or perform the various exemplary logic blocks, modules, and circuits described in connection with this disclosure. Processor 301 may also be a combination that implements computing functionality, e.g., comprising one or more microprocessor combinations, a combination of a DSP and a microprocessor, etc.
Bus 302 may include a path to transfer information between the components. Bus 302 may be a PCI (Peripheral Component Interconnect, peripheral component interconnect Standard) bus or an EISA (Extended Industry Standard Architecture ) bus, or the like. Bus 302 may be divided into an address bus, a data bus, a control bus, and the like. For ease of illustration, only one thick line is shown in fig. 3, but not only one bus or one type of bus.
The Memory 303 may be a ROM (Read Only Memory) or other type of static storage device that can store static information and instructions, a RAM (Random Access Memory ) or other type of dynamic storage device that can store information and instructions, an EEPROM (Electrically Erasable Programmable Read Only Memory ), a CD-ROM (Compact Disc Read Only Memory, compact disc Read Only Memory) or other optical disk storage, optical disk storage (including compact discs, laser discs, optical discs, digital versatile discs, blu-ray discs, etc.), magnetic disk storage media, other magnetic storage devices, or any other medium that can be used to carry or store a computer program and that can be Read by a computer, without limitation.
The memory 303 is used for storing a computer program for executing the embodiments of the present application, and is controlled to be executed by the processor 301. The processor 301 is arranged to execute a computer program stored in the memory 303 for carrying out the steps shown in the previous method embodiments.
Embodiments of the present application provide a computer readable storage medium having a computer program stored thereon, where the computer program, when executed by a processor, may implement the steps and corresponding content of the foregoing method embodiments.
In the description of embodiments of the present invention, the terms "first," "second," "third," "fourth" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defining "a first", "a second", "a third" and a fourth "may explicitly or implicitly include one or more such feature. In the description of the present invention, unless otherwise indicated, the meaning of "a plurality" is two or more.
The foregoing is merely an optional implementation manner of the implementation scenario of the application, and it should be noted that, for those skilled in the art, other similar implementation manners based on the technical ideas of the application are adopted without departing from the technical ideas of the application, and also belong to the protection scope of the embodiments of the application.
Claims (7)
1. A phase-generating carrier demodulation method, comprising:
determining a quadrature component corresponding to an interference signal, the interference signal comprising a high frequency carrier signal and a signal to be measured carried on the high frequency carrier signal, comprising: mixing the interference signal with a single frequency multiplication carrier signal and a double frequency multiplication carrier signal respectively, and then carrying out low-pass filtering to obtain a pair of orthogonal components containing the signal to be detected, wherein the orthogonal components contain sine components and cosine components which are mutually orthogonal; wherein the sinusoidal component includes a first amplitude and a corresponding first bias for characterizing a change in the sinusoidal component over time; the cosine component includes a second amplitude and a corresponding second bias characterizing a time-dependent change of the cosine component;
and carrying out ellipse fitting on the orthogonal components to obtain ellipse parameters, wherein the ellipse parameters comprise: carrying out ellipse fitting on orthogonal components at any time, and carrying out nonlinear solving by utilizing a pre-established ellipse fitting model to output ellipse parameters corresponding to the time when an observation function fitted in the ellipse fitting model is minimized, wherein the observation function is used for indicating deviation between the orthogonal components and a fitted ellipse equation;
non-linear correction is carried out on the orthogonal components based on the ellipse parameters to obtain mutually orthogonal signal components to be processed, and the method comprises the following steps: determining mapping relations between each ellipse parameter and the first amplitude, the first bias, the second amplitude and the second bias; based on an ellipse parameter obtained by ellipse fitting and the mapping relation, the first amplitude, the first bias, the second amplitude and the second bias are obtained, so that the sine component and the cosine component are corrected to obtain a signal component to be processed;
and performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be detected.
2. The phase-generating carrier demodulation method according to claim 1, wherein the non-linear solution using a pre-established elliptic fitting model comprises:
and for solving any ellipse parameter, iteratively predicting the current optimal state vector for representing the ellipse parameter by using a prediction model in the ellipse fitting model to obtain the next optimal state vector, and taking the current optimal state vector as the output of the ellipse fitting model until the difference value of the optimal state vectors obtained in two adjacent times reaches a preset stable condition.
3. The phase generating carrier demodulation method according to claim 2, wherein the nonlinear correction is performed on the orthogonal components based on the ellipse parameters to obtain mutually orthogonal signal components to be processed, further comprising:
and inputting the ellipse parameters into the observation function to correct the signal components which are not mutually orthogonal, so as to obtain mutually orthogonal signal components to be processed.
4. The phase generating carrier demodulation method according to claim 1, wherein the performing arctangent demodulation processing on the signal component to be processed to obtain phase information in the signal to be detected includes:
dividing the mutually orthogonal signal components to be processed to obtain tangent signals;
and carrying out high-pass filtering after carrying out arc tangent operation on the tangent signal so as to demodulate a signal to be detected and obtain phase information in the signal to be detected.
5. A phase generating carrier demodulation apparatus, comprising:
the orthogonal component determining module is configured to determine an orthogonal component corresponding to an interference signal, where the interference signal includes a high-frequency carrier signal and a signal to be measured carried on the high-frequency carrier signal, and includes: mixing the interference signal with a single frequency multiplication carrier signal and a double frequency multiplication carrier signal respectively, and then carrying out low-pass filtering to obtain a pair of orthogonal components containing the signal to be detected, wherein the orthogonal components contain sine components and cosine components which are mutually orthogonal; wherein the sinusoidal component includes a first amplitude and a corresponding first bias for characterizing a change in the sinusoidal component over time; the cosine component includes a second amplitude and a corresponding second bias characterizing a time-dependent change of the cosine component;
the ellipse fitting module is configured to perform ellipse fitting on the orthogonal component to obtain an ellipse parameter, and includes: carrying out ellipse fitting on orthogonal components at any time, and carrying out nonlinear solving by utilizing a pre-established ellipse fitting model to output ellipse parameters corresponding to the time when an observation function fitted in the ellipse fitting model is minimized, wherein the observation function is used for indicating deviation between the orthogonal components and a fitted ellipse equation;
the nonlinear correction module is configured to perform nonlinear correction on the orthogonal component based on the ellipse parameter to obtain mutually orthogonal signal components to be processed, and includes: determining mapping relations between each ellipse parameter and the first amplitude, the first bias, the second amplitude and the second bias; based on an ellipse parameter obtained by ellipse fitting and the mapping relation, the first amplitude, the first bias, the second amplitude and the second bias are obtained, so that the sine component and the cosine component are corrected to obtain a signal component to be processed;
and the arctangent demodulation module is used for carrying out arctangent demodulation processing on the signal component to be processed so as to obtain phase information in the signal to be detected.
6. An electronic device comprising a memory, a processor and a computer program stored on the memory, characterized in that the processor executes the computer program to carry out the steps of the phase-generating carrier demodulation method of any one of claims 1-4.
7. A computer readable storage medium having stored thereon a computer program, characterized in that the computer program when executed by a processor realizes the steps of the phase generating carrier demodulation method according to any one of claims 1-4.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311466287.2A CN117201250B (en) | 2023-11-07 | 2023-11-07 | Phase generation carrier demodulation method and device, electronic equipment and storage medium |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202311466287.2A CN117201250B (en) | 2023-11-07 | 2023-11-07 | Phase generation carrier demodulation method and device, electronic equipment and storage medium |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117201250A CN117201250A (en) | 2023-12-08 |
CN117201250B true CN117201250B (en) | 2024-01-23 |
Family
ID=88998291
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202311466287.2A Active CN117201250B (en) | 2023-11-07 | 2023-11-07 | Phase generation carrier demodulation method and device, electronic equipment and storage medium |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117201250B (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002060114A2 (en) * | 2001-01-09 | 2002-08-01 | Qualcomm Incorporated | Filter for multicarrier signals |
CN109000782A (en) * | 2018-09-27 | 2018-12-14 | 哈尔滨工程大学 | A kind of ellipse fitting non-linear error calibration method based on Kalman filtering |
CN110836638A (en) * | 2019-10-28 | 2020-02-25 | 哈尔滨工业大学 | Carrier phase delay and associated amplitude modulation elimination method in phase generation carrier arc tangent |
CN111711589A (en) * | 2020-06-15 | 2020-09-25 | 中国人民解放军海军航空大学 | Continuous phase modulation and demodulation method based on elliptical spherical wave signal |
IT202000019267A1 (en) * | 2020-08-05 | 2022-02-05 | Scuola Superiore Santanna | PHASE GENERATION CARRIER DEMODULATION METHOD WITH DISTORTION CORRECTION, USING A MULTIPLE HARMONIC MIXING |
CN115834312A (en) * | 2022-11-18 | 2023-03-21 | 西安石油大学 | PGC-Arctan-based modulation method and device for eliminating carrier phase delay and compensating modulation depth demodulation |
CN116683997A (en) * | 2023-05-30 | 2023-09-01 | 上海大学 | Method for demodulating low-frequency signal by improving phase generation carrier algorithm |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050027489A1 (en) * | 2002-08-14 | 2005-02-03 | Yale University | Phase extraction between coupled atom interferometers using ellipse-specific fitting |
JP2010213107A (en) * | 2009-03-11 | 2010-09-24 | Fujitsu Semiconductor Ltd | Communication apparatus |
US10476731B2 (en) * | 2015-12-22 | 2019-11-12 | Intel Corporation | Joint noncoherent demodulation and carrier frequency offset correction based on non-linear filtering |
-
2023
- 2023-11-07 CN CN202311466287.2A patent/CN117201250B/en active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2002060114A2 (en) * | 2001-01-09 | 2002-08-01 | Qualcomm Incorporated | Filter for multicarrier signals |
CN109000782A (en) * | 2018-09-27 | 2018-12-14 | 哈尔滨工程大学 | A kind of ellipse fitting non-linear error calibration method based on Kalman filtering |
CN110836638A (en) * | 2019-10-28 | 2020-02-25 | 哈尔滨工业大学 | Carrier phase delay and associated amplitude modulation elimination method in phase generation carrier arc tangent |
CN111711589A (en) * | 2020-06-15 | 2020-09-25 | 中国人民解放军海军航空大学 | Continuous phase modulation and demodulation method based on elliptical spherical wave signal |
IT202000019267A1 (en) * | 2020-08-05 | 2022-02-05 | Scuola Superiore Santanna | PHASE GENERATION CARRIER DEMODULATION METHOD WITH DISTORTION CORRECTION, USING A MULTIPLE HARMONIC MIXING |
CN115834312A (en) * | 2022-11-18 | 2023-03-21 | 西安石油大学 | PGC-Arctan-based modulation method and device for eliminating carrier phase delay and compensating modulation depth demodulation |
CN116683997A (en) * | 2023-05-30 | 2023-09-01 | 上海大学 | Method for demodulating low-frequency signal by improving phase generation carrier algorithm |
Non-Patent Citations (1)
Title |
---|
参数估计误差对光纤水听器PGC反正切解调性能的影响;王凯;施清平;蒋佳佳;戴之光;田长栋;段发阶;张敏;廖延彪;;光电子.激光(第10期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN117201250A (en) | 2023-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110307780B (en) | PGC phase demodulation error real-time compensation method based on iterative computation | |
CN110411486B (en) | PGC-DCDM demodulation method insensitive to phase delay and modulation depth | |
CN1649334B (en) | Apparatus and method for adjusting quadrature modulator, and acommunication apparatus | |
CN109855719B (en) | High-stability interference type optical fiber hydrophone signal demodulation method | |
JP3737819B2 (en) | Carrier quadrature error detection method and quadrature modulation apparatus for quadrature modulator | |
CN106291105A (en) | A kind of sweep generator based on digital zero intermediate frequency | |
CN117201250B (en) | Phase generation carrier demodulation method and device, electronic equipment and storage medium | |
US7068983B2 (en) | Method for detecting quadrature modulator carrier leak adjusting point by geometrical analysis/calculation method, carrier leak adjusting method, and quadrature modulation apparatus | |
JP6415825B2 (en) | Method and apparatus for determining a DC offset in an IQ modulator | |
CN109521269B (en) | Amplitude modulation signal digital frequency measurement method | |
CN113922894B (en) | Multi-channel phase calibration method, system, medium and terminal of quadrature modulation digital receiver | |
CN117109719A (en) | Vibration velocity demodulation system and demodulation method of laser Doppler vibration meter | |
WO2022110832A1 (en) | Azimuthal electromagnetic wave logging while drilling signal processing method and apparatus, and storage medium | |
CN115834312A (en) | PGC-Arctan-based modulation method and device for eliminating carrier phase delay and compensating modulation depth demodulation | |
CN110060700B (en) | Short sequence audio analysis method based on parameter spectrum estimation | |
CN118857110A (en) | Nonlinear error correction method and device for laser measurement system | |
CN111130660B (en) | Phase noise measuring probe circuit | |
CN109510661B (en) | Method and device for measuring IQ delay difference in optical transmitter and optical transmitter | |
CN114630477A (en) | Stepping phase comparison method, device, equipment and storage medium | |
JP5783539B2 (en) | Identification apparatus and identification method | |
CN114184837B (en) | Instantaneous frequency measurement method based on Cordic algorithm | |
RU2799110C1 (en) | Method for measuring signal to noise ratio of signals with angular manipulation | |
CN116538909A (en) | Sinusoidal phase modulation interferometer PGC demodulation method, system, device and storage medium | |
CN116953596A (en) | Electrical parameter calculation method and device for non-whole period | |
CN115993087A (en) | Real-time phase delay compensation method in single-frequency laser interferometry |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right | ||
PE01 | Entry into force of the registration of the contract for pledge of patent right |
Denomination of invention: A phase generation carrier demodulation method, device, electronic device, and storage medium Granted publication date: 20240123 Pledgee: Guanggu Branch of Wuhan Rural Commercial Bank Co.,Ltd. Pledgor: WUHAN CPCTECH CO.,LTD. Registration number: Y2024980015042 |