CN107104740A - A kind of signal processing method and system - Google Patents
A kind of signal processing method and system Download PDFInfo
- Publication number
- CN107104740A CN107104740A CN201710029044.0A CN201710029044A CN107104740A CN 107104740 A CN107104740 A CN 107104740A CN 201710029044 A CN201710029044 A CN 201710029044A CN 107104740 A CN107104740 A CN 107104740A
- Authority
- CN
- China
- Prior art keywords
- function
- signal
- msup
- output signal
- mrow
- 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 13
- 238000012545 processing Methods 0.000 claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 19
- 230000006870 function Effects 0.000 claims description 112
- 230000009131 signaling function Effects 0.000 claims description 25
- 238000001914 filtration Methods 0.000 claims description 15
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000003287 optical effect Effects 0.000 claims description 6
- 239000004576 sand Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 abstract description 6
- 238000002474 experimental method Methods 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 239000013307 optical fiber Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/60—Receivers
- H04B10/66—Non-coherent receivers, e.g. using direct detection
- H04B10/69—Electrical arrangements in the receiver
- H04B10/693—Arrangements for optimizing the preamplifier in the receiver
- H04B10/6931—Automatic gain control of the preamplifier
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B10/00—Transmission systems employing electromagnetic waves other than radio-waves, e.g. infrared, visible or ultraviolet light, or employing corpuscular radiation, e.g. quantum communication
- H04B10/60—Receivers
- H04B10/66—Non-coherent receivers, e.g. using direct detection
- H04B10/69—Electrical arrangements in the receiver
- H04B10/697—Arrangements for reducing noise and distortion
Landscapes
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
The invention discloses a kind of signal processing method and system.This method includes setting up the function of initial signal, and is demodulated according to the function pair initial signal of initial signal, to obtain the function of output signal;The function of output signal includes amplitude factor;Amplitude factor is estimated, and the function of output signal calculated according to the estimated result of amplitude factor, to obtain signal to be detected.The system includes first processing module and Second processing module.A kind of signal processing method and system that the present invention is provided, can real-time be estimated the amplitude factor because of interference of light intensity and modulation depth generation, amplitude correction be carried out to output signal, to obtain more accurately signal to be detected.
Description
Technical Field
The invention relates to the technical field of sensors, in particular to a signal processing method and a signal processing system.
Background
A phase-generated carrier (PGC) demodulation algorithm may be used to demodulate the signal of the interferometric optical fiber sensor. PGC demodulation algorithms are classified into differential-and-cross-multiplication (DCM) and Arctan (Arctan). The demodulation result of the DCM algorithm is affected by the light interference intensity and the modulation depth, and when the interference intensity and the modulation depth fluctuate, the demodulation result generates amplitude distortion. The arctan algorithm can eliminate the influence of the interference intensity on the demodulation result by dividing, but is still influenced by the modulation depth, and when the modulation depth fluctuates, the demodulation result generates nonlinearity to cause serious harmonic distortion.
For the influence of the optical interference intensity, the traditional DCM algorithm is implemented by adding a complex Automatic Gain Control (AGC) circuit, and the effect thereof has a great dependence on the performance of devices in the circuit and a limited improvement on the performance. In a digital PGC system, the compensation of the light interference intensity can be more flexible, and when the influence of the low-pass filter and the light interference intensity is researched, the modulation depth is assumed to be unchanged, and a trigonometric function is used for calculating a correction coefficient to eliminate the influence. In order to solve harmonic distortion possibly caused by an arc tangent algorithm, a method is provided for calculating a correction parameter by utilizing differential self-multiplication and demodulating a signal by combining an anti-cotangent method.
In summary, both of the current algorithms can eliminate the influence of interference intensity; the DCM algorithm has no scheme for solving the influence of modulation depth on a demodulation signal; the arctan algorithm currently has a relevant scheme to address the effects of modulation depth. However, one of the main problems of the arctangent method is the non-monotonicity of the tangent function, which affects the arctangent demodulation result.
Disclosure of Invention
The invention aims to solve the problem that the fluctuation of light interference intensity and modulation depth influences the signal demodulation result in the prior sensor technology, and provides a signal processing method and a signal processing system, which can carry out real-time and rapid estimation on an amplitude factor generated by the light interference intensity and the modulation depth and carry out amplitude correction on an output signal so as to obtain a more accurate signal to be detected.
To achieve the above object, in one aspect, the present invention provides a signal processing method. Establishing a function of an initial signal, and demodulating the initial signal according to the function of the initial signal to obtain a function of an output signal; the function of the output signal contains an amplitude factor; and estimating the amplitude factor, and calculating a function of the output signal according to the estimation result of the amplitude factor to acquire the signal to be detected.
Preferably, the demodulating the initial signal according to the function of the initial signal to obtain the function of the output signal specifically includes: mixing the modulation signal and the initial signal according to the function of the initial signal, and obtaining a low-frequency signal function through low-pass filtering; a function of the low frequency signal is calculated to obtain a function of the output signal.
Preferably, the function of calculating the low frequency signal function to obtain the output signal specifically comprises: calculating a low-frequency signal function by adopting a differential cross multiplication algorithm to obtain a first calculation result; and calculating the first calculation result by adopting a difference algorithm, an integral algorithm and high-pass filtering to obtain a function of the output signal.
Preferably, estimating the amplitude factor, and calculating a function of the output signal according to an estimation result of the amplitude factor to obtain the signal to be detected specifically includes: deducing a positive elliptic function according to the low-frequency signal function, and fitting the positive elliptic function by using a least square method so as to estimate related parameters of the positive elliptic function; and estimating the amplitude factor according to the estimation result of the related parameters, and acquiring the signal to be detected according to the estimation result of the amplitude factor and the function of the output signal.
Preferably, the function of the output signal is:
V0=B2GHJ1(C)J2(C)φs(t)
wherein, V0To output a signal, B2GHJ1(C)J2(C) Is the amplitude factor, B is the optical interference intensity, C is the modulation depth, G and H are the amplitude of the modulation signal, J1(C) And J2(C) Are the magnitudes of the Bessel function of order 1 and 2, respectively, phisAnd (t) is a signal to be detected.
Preferably, the positive elliptic function is:
wherein, I1fAnd I2fIs a low-pass filtered low-frequency signal, B is the optical interference intensity, C is the modulation depth, G and H are the amplitude of the modulation signal, J1(C) And J2(C) Are the 1 st order and 2 nd order bessel function vector values, respectively.
In another aspect, the present invention provides a signal processing system. The system includes a first processing module and a second processing module. The first processing module is used for establishing a function of the initial signal and demodulating the initial signal according to the function of the initial signal to obtain a function of an output signal; the function of the output signal contains an amplitude factor. The second processing module is used for estimating the amplitude factor and calculating the function of the output signal according to the estimation result of the amplitude factor so as to obtain the signal to be detected.
Preferably, the first processing module is specifically configured to: establishing a function of an initial signal, mixing a modulation signal and the initial signal according to the function of the initial signal, and obtaining a low-frequency signal function through low-pass filtering; calculating a low-frequency signal function to obtain a function of an output signal; the function of the output signal contains an amplitude factor.
Preferably, the first processing module is specifically configured to: establishing a function of an initial signal, mixing a modulation signal and the initial signal according to the function of the initial signal, and obtaining a low-frequency signal function through filtering; calculating a low-frequency signal function by adopting a differential cross multiplication algorithm to obtain a first calculation result; and calculating the first calculation result by adopting a difference algorithm, an integral algorithm and high-pass filtering to obtain a function of the output signal.
Preferably, the second processing module is specifically configured to: deducing a positive elliptic function according to the low-frequency signal function, and fitting the positive elliptic function by using a least square method so as to estimate related parameters of the positive elliptic function; and estimating the amplitude factor according to the estimation result of the related parameter, and acquiring the signal to be detected according to the estimation result of the amplitude factor and the function of the output signal.
The signal processing method and the signal processing system provided by the invention can be used for quickly estimating the amplitude factor generated by the light interference intensity and the modulation depth in real time and correcting the amplitude of the output signal so as to obtain a more accurate signal to be detected.
Drawings
Fig. 1 is a schematic flow chart of a signal processing method according to an embodiment of the present invention;
FIG. 2 is a graph comparing the estimated value of the amplitude factor with the actual value in FIG. 1 during the experiment;
FIG. 3 is a graph comparing the estimated value of the amplitude factor of FIG. 1 with the actual value in another experiment;
fig. 4 is a schematic structural diagram of a signal processing system according to an embodiment of the present invention.
Detailed Description
The technical solution of the present invention is further described in detail by the accompanying drawings and examples.
Fig. 1 is a schematic flowchart of a signal processing method according to an embodiment of the present invention. As shown in fig. 1, the method comprises steps 101-102:
step 101, establishing a function of the initial signal, and demodulating the initial signal according to the function of the initial signal to obtain a function of the output signal, wherein the function of the output signal comprises an amplitude factor.
Specifically, a function of an initial signal is established, a modulation signal and the initial signal are adopted for frequency mixing according to the function of the initial signal, and a low-frequency signal function is obtained through low-pass filtering; calculating a low-frequency signal function by adopting a differential cross multiplication algorithm to obtain a first calculation result; and calculating the first calculation result by adopting a difference algorithm, an integral algorithm and high-pass filtering to obtain a function of the output signal.
In an optical fiber sensor adopting a phase-generated carrier (PGC) demodulation algorithm, a function of an initial signal after photoelectric conversion is:
I=A+Bcos[Ccos(ω0t)+φ(t)](1)
wherein, I is initial signal, A is interference light intensity direct current quantity, B is light interference intensity, C is modulationDepth of formation, omega0For modulating the frequency, phi (t) is the phase difference caused by the signal acting on the probe of the fiber-optic sensor, including the phase difference phi caused by the signal to be detecteds(t) phase difference due to noisen(t)。
The Bessel (Bessel) function is used to expand equation (1) to:
wherein, I is output signal, A is interference light intensity direct current quantity, B is light interference intensity, C is modulation depth, omega0Is the angular frequency, phi (t) is the phase difference, J0(C)、J2k(C) And J2k+1(C) Respectively, 0 th order, 2k order and 2k +1 th order bessel function vector values, k being 0,1,2, …, k ∈ N.
According to the formula (2), the amplitude G and the angular frequency omega are adopted0Sum amplitude of H and angular frequency of 2 omega0The modulation signals are respectively mixed with the initial signals, and the low-frequency signal functions are obtained after the modulation signals pass through low-pass filters:
I1f=-BGJ1(C)sinφ(t) (3)
I2f=-BHJ2(C)cosφ(t) (4)
in order to overcome the blanking and distortion phenomena of the signal caused by the fluctuation of the external interference signal, differential-and-cross-multiplication (DCM) is adopted. The low-frequency signals from the low-pass filter pass through a differential circuit respectively, and the function of the differentiated signals is as follows:
I1d=-BGJ1(C)φ'(t)cosφ(t) (5)
I2d=-BHJ2(C)φ'(t)sinφ(t) (6)
and (3) adopting a cross multiplication algorithm, namely multiplying two ends of the formulas (3) and (6) respectively, and multiplying two ends of the formulas (4) and (5) respectively to obtain:
I1e=B2GHJ1(C)J2(C)φ'(t)sin2φ(t) (7)
I2e=B2GHJ1(C)J2(C)φ'(t)cos2φ(t) (8)
the two paths of signals represented by the formulas (7) and (8) are subjected to differential operation through an integral operational amplifier to obtain:
V'=B2GHJ1(C)J2(C)φ'(t) (9)
and then the integral operation is carried out by an integral operational amplifier to obtain:
V=B2GHJ1(C)J2(C)φ(t)=B2GHJ1(C)J2(C)[φs(t)+φn(t)](10)
in the formula (10), phis(t) is the phase difference, phi, caused by the signal to be detectedn(t) is a phase difference due to noise. Noise is a phase disturbance term caused by external environment, and is generally a slowly-varying signal according to the following formula, so that noise can be filtered by high-pass filtering. The signal represented by formula (10) is passed through a high-pass filter to filter out noise terms, and the function of the output signal is obtained as follows:
V0=B2GHJ1(C)J2(C)φs(t) (11)
wherein, V0To output a signal, B2GHJ1(C)J2(C) Is the amplitude factor, B is the optical interference intensity, C is the modulation depth, G and H are the amplitude of the modulation signal, J1(C) And J2(C) Are the magnitudes of the Bessel function of order 1 and 2, respectively, phisAnd (t) is a signal to be detected.
In formula (11), B2GHJ1(C)J2(C) As an amplitude factor, since the intensity of light interference B and the modulation depth C are unknown and are being sensedSlow fluctuations occur during the operation of the demodulator, and therefore the amplitude factor needs to be estimated in real time during the demodulation process.
Step 102: and estimating the amplitude factor, and calculating a function of the output signal according to the estimation result of the amplitude factor to acquire the signal to be detected.
Specifically, a positive elliptic function is deduced according to a low-frequency signal function, and the positive elliptic function is fitted by utilizing minimum two multiplication so as to estimate related parameters of the positive elliptic function; and estimating the amplitude factor according to the estimation result of the related parameters, and acquiring the signal to be detected according to the estimation result of the amplitude factor and the function of the output signal.
According to the formulas (3) and (4) and the trigonometric function formula sin2φ(n)+cos2φ (n) ≡ 1 may be extrapolated to the following equation:
the function expressed by the formula (12) is a positive elliptic function, and the elliptic curve expressed by the formula (12) is fitted by using the least square method (LS), so that the relevant parameters can be estimatedAnd B2G2And then an estimated magnitude of the amplitude factor is obtained. The standard positive elliptic function can be expressed as:
x2+a1y2+a2=0 (13)
wherein x is I1f,y=I2f,a2=B2G2。
Equation (13) can be expressed in a matrix as:
YA=X(14)
wherein,n is the data length used to make the minimum two-times estimate.
From equation (14), it can be deduced that the fitting result is:
A=(YTY)-1YTX (15)
the square expression of the amplitude factor can be derived in conjunction with equation (11) as:
B4G2H2J1 2(C)J2 2(C)=a2 2/a1(16)
according to the characteristics of the first class of Bessel functions, i.e. the function value fluctuates around the point 0 with the change of the independent variable, B2J1(C)J2(C) The value of (c) may be positive or negative. In order to reduce the dependence of the output result on the Bessel function, the modulation depth needs to be properly selected, so that J is enabled1(C)J2(C) Maximum, for example, the value of C can be limited to [0,3.83 ]]Within the interval of B2J1(C)J2(C) Take a positive value.
And (3) carrying out amplitude correction on the formula (11) by using the formula (17) to obtain a signal to be detected as follows:
the following experiment further proves the beneficial effects of the signal processing method provided by the invention.
In the experiment, the sampling frequency of the signal is 200KHZ, the amplitude is G-1, and the angular frequency is omega020kHz, amplitude H1, angular frequency 2 ω0The 40kHz modulated signal is mixed with the original signal separately. The settings of the low-pass filter are: an equal ripple method is adopted to design an FIR filter, the passband frequency is 8kHz, the stopband frequency is 12kHz, the stopband attenuation is 80dB, and the order of the filter is 127. The data length N used for ellipse fitting is 2000, i.e. the amplitude factor estimation update interval is 10 ms.
Assuming that the signal to be detected is a single frequency signal of 1kHz, the amplitude of the phase change caused is 0.8rad, and the estimation result of the amplitude factor is shown in fig. 2, and the relative error between the estimation result and the actual value is 0.5%.
Assuming that the signal to be detected is gaussian noise below 2kHz, the phase change amplitude caused by the signal to be detected is 1rad, and the estimation result of the amplitude factor is shown in fig. 3, with a relative error of 0.93%. It can be seen that the present invention can estimate the amplitude factor due to the intensity of light interference and modulation depth for different types of signals to be detected.
Therefore, the signal processing method provided by the embodiment of the invention can quickly estimate the amplitude factor generated by the light interference intensity and the modulation depth in real time, and perform amplitude correction on the output signal to obtain a more accurate signal to be detected.
Fig. 4 is a signal processing system according to an embodiment of the present invention, and as shown in fig. 4, the system 200 includes a first processing module 201 and a second processing module 202:
the first processing module 201 is configured to establish a function of an initial signal, and demodulate the initial signal according to the function of the initial signal to obtain a function of an output signal; the function of the output signal contains an amplitude factor.
Specifically, the first processing module is configured to establish a function of an initial signal, mix frequency with the initial signal using a modulation signal according to the function of the initial signal, and obtain a low-frequency signal function through filtering. Calculating a low-frequency signal function by adopting a differential cross multiplication algorithm to obtain a first calculation result; and calculating the first calculation result by adopting a difference algorithm, an integral algorithm and high-pass filtering to obtain a function of the output signal.
The second processing module 202 is configured to estimate the amplitude factor, and calculate a function of the output signal according to an estimation result of the amplitude factor to obtain a signal to be detected.
Specifically, the second processing module is configured to derive a positive elliptic function according to the low-frequency signal function, and fit the positive elliptic function by using a least square method to estimate relevant parameters of the positive elliptic function; and estimating the amplitude factor according to the estimation result of the related parameter, and acquiring the signal to be detected according to the estimation result of the amplitude factor and the function of the output signal.
The embodiment of the system provided by the present invention corresponds to the embodiment of the method provided by the present invention, and is not described herein again.
The signal processing system provided by the embodiment of the invention can quickly estimate the amplitude factor generated by the light interference intensity and the modulation depth in real time, and correct the amplitude of the output signal to obtain a more accurate signal to be detected.
The above embodiments are intended to illustrate the objects, aspects and advantages of the present invention in further detail, and it should be understood that the above embodiments are only exemplary embodiments of the present invention and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (10)
1. A signal processing method, comprising the steps of:
establishing a function of an initial signal, and demodulating the initial signal according to the function of the initial signal to obtain a function of an output signal; the function of the output signal comprises an amplitude factor;
and estimating the amplitude factor, and calculating a function of the output signal according to the estimation result of the amplitude factor to acquire a signal to be detected.
2. The method according to claim 1, wherein the demodulating the initial signal as a function of the initial signal to obtain the function of the output signal specifically comprises:
mixing a modulation signal with the initial signal according to the function of the initial signal, and obtaining a low-frequency signal function through low-pass filtering; the low frequency signal function is calculated to obtain a function of the output signal.
3. The method according to claim 2, wherein said calculating the function of the low frequency signal to obtain the output signal comprises in particular:
calculating the low-frequency signal function by adopting a differential cross multiplication algorithm to obtain a first calculation result; and calculating the first calculation result by adopting a difference algorithm, an integral algorithm and high-pass filtering to obtain a function of the output signal.
4. The method according to claim 2, wherein the estimating the amplitude factor and calculating the function of the output signal according to the estimation result of the amplitude factor to obtain the signal to be detected specifically comprises:
deriving a positive elliptic function according to the low-frequency signal function, and fitting the positive elliptic function by using a least square method so as to estimate related parameters of the positive elliptic function; and estimating the amplitude factor according to the estimation result of the related parameter, and acquiring a signal to be detected according to the estimation result of the amplitude factor and the function of the output signal.
5. The method of claim 2, wherein the function of the output signal is:
V0=B2GHJ1(C)J2(C)φs(t)
wherein, V0To output a signal, B2GHJ1(C)J2(C) Is the amplitude factor, B is the optical interference intensity, C is the modulation depth, G and H are the amplitude of the modulation signal, J1(C) And J2(C) Are the magnitudes of the Bessel function of order 1 and 2, respectively, phisAnd (t) is a signal to be detected.
6. The method of claim 4, wherein the positive elliptic function is:
<mrow> <msup> <msub> <mi>I</mi> <mrow> <mn>1</mn> <mi>f</mi> </mrow> </msub> <mn>2</mn> </msup> <mo>+</mo> <mfrac> <mrow> <msup> <mi>G</mi> <mn>2</mn> </msup> <msup> <msub> <mi>J</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>C</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msup> <mi>H</mi> <mn>2</mn> </msup> <msup> <msub> <mi>J</mi> <mn>2</mn> </msub> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>C</mi> <mo>)</mo> </mrow> </mrow> </mfrac> <msup> <msub> <mi>I</mi> <mrow> <mn>2</mn> <mi>f</mi> </mrow> </msub> <mn>2</mn> </msup> <mo>-</mo> <msup> <mi>B</mi> <mn>2</mn> </msup> <msup> <mi>G</mi> <mn>2</mn> </msup> <msup> <msub> <mi>J</mi> <mn>1</mn> </msub> <mn>2</mn> </msup> <mrow> <mo>(</mo> <mi>C</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>0</mn> </mrow>
wherein, I1fAnd I2fIs a low-pass filtered low-frequency signal, B is the optical interference intensity, C is the modulation depth, G and H are the amplitude of the modulation signal, J1(C) And J2(C) Are the 1 st order and 2 nd order bessel function vector values, respectively.
7. A signal processing system, comprising:
the first processing module is used for establishing a function of an initial signal and demodulating the initial signal according to the function of the initial signal so as to obtain a function of an output signal; the function of the output signal comprises an amplitude factor;
and the second processing module is used for estimating the amplitude factor and calculating the function of the output signal according to the estimation result of the amplitude factor so as to acquire the signal to be detected.
8. The system of claim 7, wherein the first processing module is specifically configured to:
establishing a function of an initial signal, mixing a modulation signal and the initial signal according to the function of the initial signal, and obtaining a low-frequency signal function through low-pass filtering; calculating the low-frequency signal function to obtain a function of an output signal; the function of the output signal comprises an amplitude factor.
9. The system of claim 8, wherein the first processing module is specifically configured to:
establishing a function of an initial signal, mixing a modulation signal and the initial signal according to the function of the initial signal, and obtaining a low-frequency signal function through filtering; calculating the low-frequency signal function by adopting a differential cross multiplication algorithm to obtain a first calculation result; and calculating the first calculation result by adopting a difference algorithm, an integral algorithm and high-pass filtering to obtain a function of the output signal.
10. The system of claim 8, wherein the second processing module is specifically configured to:
deriving a positive elliptic function according to the low-frequency signal function, and fitting the positive elliptic function by using a least square method so as to estimate related parameters of the positive elliptic function; and estimating the amplitude factor according to the estimation result of the related parameter, and acquiring a signal to be detected according to the estimation result of the amplitude factor and the function of the output signal.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710029044.0A CN107104740A (en) | 2017-01-16 | 2017-01-16 | A kind of signal processing method and system |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710029044.0A CN107104740A (en) | 2017-01-16 | 2017-01-16 | A kind of signal processing method and system |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107104740A true CN107104740A (en) | 2017-08-29 |
Family
ID=59676488
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710029044.0A Pending CN107104740A (en) | 2017-01-16 | 2017-01-16 | A kind of signal processing method and system |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107104740A (en) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102434783A (en) * | 2011-10-09 | 2012-05-02 | 中国计量学院 | Device for detecting leakage of optical fiber underwater long-distance transmission pipeline based on phase generator carrier demodulation |
CN203115505U (en) * | 2013-01-16 | 2013-08-07 | 上海健驰物联网科技有限公司 | Optical fiber leak detection device for underwater long-distance pipeline |
CN105157733A (en) * | 2015-06-02 | 2015-12-16 | 哈尔滨工程大学 | Improved phase generated carrier (PGC) demodulation method |
CN105486331A (en) * | 2015-12-01 | 2016-04-13 | 哈尔滨工程大学 | High-precision optical signal phase demodulation system and demodulation method |
CN105486225A (en) * | 2015-12-01 | 2016-04-13 | 哈尔滨工程大学 | Phase demodulation device and demodulation method capable of suppressing light-intensity fluctuation noise |
CN205808544U (en) * | 2016-06-02 | 2016-12-14 | 长沙湘计海盾科技有限公司 | A kind of fibre optic hydrophone acceleration sensitivity test system |
-
2017
- 2017-01-16 CN CN201710029044.0A patent/CN107104740A/en active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102434783A (en) * | 2011-10-09 | 2012-05-02 | 中国计量学院 | Device for detecting leakage of optical fiber underwater long-distance transmission pipeline based on phase generator carrier demodulation |
CN203115505U (en) * | 2013-01-16 | 2013-08-07 | 上海健驰物联网科技有限公司 | Optical fiber leak detection device for underwater long-distance pipeline |
CN105157733A (en) * | 2015-06-02 | 2015-12-16 | 哈尔滨工程大学 | Improved phase generated carrier (PGC) demodulation method |
CN105486331A (en) * | 2015-12-01 | 2016-04-13 | 哈尔滨工程大学 | High-precision optical signal phase demodulation system and demodulation method |
CN105486225A (en) * | 2015-12-01 | 2016-04-13 | 哈尔滨工程大学 | Phase demodulation device and demodulation method capable of suppressing light-intensity fluctuation noise |
CN205808544U (en) * | 2016-06-02 | 2016-12-14 | 长沙湘计海盾科技有限公司 | A kind of fibre optic hydrophone acceleration sensitivity test system |
Non-Patent Citations (1)
Title |
---|
张爱玲等: "干涉型光纤传感器PGC解调算法的研究", 《光电技术应用》 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105791179B (en) | Sampling frequency offset compensation device and method | |
WO2015069359A3 (en) | System and method for gyroscope zero-rate-offset drift reduction through demodulation phase error correction | |
CN101005480B (en) | Demodulation circuit and demodulation method | |
JP4030482B2 (en) | I / Q demodulation circuit | |
CN109862545B (en) | Frequency offset compensation method and device of Bluetooth signal, computer equipment and storage medium | |
US8681896B1 (en) | Transmitter I/Q and carrier leak calibration | |
JP3371876B2 (en) | Demodulator with automatic quadrature control function | |
JP2017112606A5 (en) | ||
JPH09130337A (en) | Digital radio receiver | |
JP5935631B2 (en) | Compensation device and wireless communication device | |
CN107104740A (en) | A kind of signal processing method and system | |
JP4821379B2 (en) | Demodulator, distance measuring device, and data receiving device | |
TW201631939A (en) | Maximum likelihood sequence detection in the phase domain | |
JP5598243B2 (en) | FM detector, signal interpolation method, and program | |
CN111212006B (en) | Signal demodulation system and demodulation method | |
WO2016062033A1 (en) | Frequency compensation processing method and device | |
CN108111451B (en) | Frequency offset estimation method, system and controller | |
TW201442472A (en) | Receiver with inphase-quadrature imbalance compensation and inphase-quadrature imbalance compensation method thereof | |
US20170180182A1 (en) | Joint noncoherent demodulation and carrier frequency offset correction based on non-linear filtering | |
CN106969787A (en) | Signal processing method and system | |
TWI528352B (en) | Method and apparatus for multiplex signal decoding | |
CN108270705B (en) | Demodulation device and demodulation method for frequency modulation signal | |
JP3522347B2 (en) | Transmission system and receiver | |
JP2000295304A (en) | Demodulator | |
JP2000244592A (en) | Demodulator |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170829 |