CN115112964A - Method and device for recognizing characteristic current code bits after noise addition - Google Patents
Method and device for recognizing characteristic current code bits after noise addition Download PDFInfo
- Publication number
- CN115112964A CN115112964A CN202110294524.6A CN202110294524A CN115112964A CN 115112964 A CN115112964 A CN 115112964A CN 202110294524 A CN202110294524 A CN 202110294524A CN 115112964 A CN115112964 A CN 115112964A
- Authority
- CN
- China
- Prior art keywords
- sequence
- code bit
- index
- num
- characteristic current
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 44
- 238000005070 sampling Methods 0.000 claims description 42
- 238000009499 grossing Methods 0.000 claims description 30
- 230000003321 amplification Effects 0.000 claims description 27
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 27
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 abstract description 4
- 230000005611 electricity Effects 0.000 abstract description 3
- 238000010586 diagram Methods 0.000 description 10
- 230000008569 process Effects 0.000 description 10
- 230000008859 change Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- -1 (i) -F Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R31/00—Arrangements for testing electric properties; Arrangements for locating electric faults; Arrangements for electrical testing characterised by what is being tested not provided for elsewhere
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R19/00—Arrangements for measuring currents or voltages or for indicating presence or sign thereof
- G01R19/0092—Arrangements for measuring currents or voltages or for indicating presence or sign thereof measuring current only
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Dc Digital Transmission (AREA)
Abstract
The invention relates to a method and a device for identifying a code bit of a noise-added characteristic current, which can identify the code bit of the characteristic current added with noises such as Gaussian white noise, pulse noise, user electricity consumption and the like. The cross-correlation operation is to perform cross-correlation operation on the noisy signal and a locally recorded sequence to obtain a new signal, so that code bit identification of the characteristic current can be realized without adopting a digital filter. Compared with the recognition method in the prior art, the code bit recognition method provided by the invention can recognize the code bit of the characteristic current after the noise is added without using a digital filter, thereby reducing the calculation amount required by the algorithm.
Description
Technical Field
The invention relates to the technical field of identification of a line topological structure of a station area, in particular to a method and a device for identifying a characteristic current code bit after noise is added.
Background
At present, the identification platform area circuit topological structure mainly has two large directions, one is a scheme which can directly identify and judge by directly utilizing a power grid acquisition system without hardware, and the other is a scheme which needs to add new hardware equipment to assist in completing circuit topological identification. In terms of hardware direction, there are various solutions for identifying the topology of the line in the cell, and one of them is to identify the topology of the line in the cell based on the characteristic current, see patent application CN110646690A, which describes the implementation of the hardware solution and system architecture for identifying the topology of the line in the cell.
For identifying the circuit topology based on the characteristic current, two methods are generally adopted for identifying the characteristic current after noise addition, one method is to convert a time domain signal into a frequency domain by using DFT or FFT, and judge and identify information transmitted by the characteristic current through the frequency domain signal and the amplitude on the characteristic frequency domain. The other is to reserve the frequency of the characteristic current through a digital filter and then use the autocorrelation or cross-correlation of the signal to judge the information conveyed by the identification characteristic current. In the technical scheme in the prior art, because a digital filter is required to identify the characteristic code bits, the calculation amount of an algorithm is large, the complexity is high, and the code bit identification speed is influenced.
Disclosure of Invention
Based on the foregoing situation in the prior art, an object of the present invention is to provide a method and an apparatus for recognizing code bits of a noise-added feature current, which can recognize code bits of a noise-added feature current without using a digital filter, and reduce the amount of computation required by an algorithm, compared with a recognition method in the prior art.
In order to achieve the above object, according to one aspect of the present invention, there is provided a method for identifying a code bit of a characteristic current after noise addition, comprising the steps of:
performing cross-correlation operation on the characteristic current subjected to noise addition to obtain a sequence C;
carrying out operation simplification processing on the sequence C to obtain a sequence D;
carrying out amplification smoothing treatment on the sequence D to obtain a sequence F;
and carrying out code bit recognition on the sequence F to obtain a code bit recognition result.
Further, the performing a cross-correlation operation on the noise-added characteristic current includes:
sampling the characteristic current with the duration of t after the noise is added by using a sampling frequency f to obtain a sequence A;
sampling characteristic currents with a period of T and code bits of 1 according to the same sampling frequency f to obtain a sequence a;
and performing cross-correlation operation on the sequences A and a according to the following formula to obtain a sequence C:
wherein C (i) is an element of the sequence C, n a Number of sampling points of sequence a, n a ＝T*f。
Further, the performing operation simplification processing on the sequence C includes:
for sequence C, start at point 1 and every n s Selecting the largest point as an element of the sequence D to form the sequence D;
wherein n is s The number of points occupied by a 50Hz sine wave in each sampling period, n s ＝f/50。
Further, the performing amplification smoothing processing on the sequence D includes:
constructing a sequence e, wherein the sequence e comprises T x f elements, each element is 1, and the number of the sequence e is n e ；
And performing cross-correlation operation on the sequence D and the sequence e according to the following formula to obtain a sequence F:
wherein F (i) is an element in the sequence F, and the number n of the sequence F F ＝n D -n e 。
Further, the code bit recognition on the sequence F includes:
acquiring the median of the sequence F, and recording as peak _ threshold;
let m be the index of sequence F, find the first occurrence of F (m +1) -F (m)<0, and F (m) -F (m-1)>0, and F (m)>The index of peak _ threshold, denoted as m ═ a; wherein m is 1,2,3 … … … n F ；
Finding out m when the value of F (i) is maximum from the index m ═ a to m ═ a + one _ cycle _ num, and marking the index m ═ begin _ index; wherein, one _ cycle _ num is T/0.02;
acquiring the total code number all _ num of the sequence F; wherein all _ num is T/T;
let j be 1 and let j be 1,
when begin _ index + (j-1) one _ cycle _ num is not more than n F When the temperature of the water is higher than the set temperature,
f (begin _ index + (j-1) one _ cycle _ num) > peak _ thermoshiftd, then the code bit of sequence F (i) at that point is 1, otherwise 0;
when begin _ index + (j-1) one _ cycle _ num>n F When the temperature of the water is higher than the set temperature,
F(n F )>peak thermally, the code bit of sequence f (i) is 1, otherwise it is 0;
if j is less than all + num, j is j +1, and the steps are repeatedly executed; otherwise, ending.
According to another aspect of the invention, a device for identifying code bits of a noisy characteristic current is provided, which comprises a cross-correlation operation module, an operation simplification processing module, an amplification smoothing processing module and a code bit identification module; wherein,
the cross-correlation operation module is used for performing cross-correlation operation on the characteristic current subjected to noise addition to obtain a sequence C;
the operation simplification processing module is used for carrying out operation simplification processing on the sequence C to obtain a sequence D;
the amplification smoothing processing module is used for carrying out amplification smoothing processing on the sequence D to obtain a sequence F;
and the code bit identification module is used for carrying out code bit identification on the sequence F to obtain a code bit identification result.
Further, the cross-correlation operation module performs cross-correlation operation on the noisy characteristic current, and includes:
sampling the characteristic current with the duration of t after the noise is added by using a sampling frequency f to obtain a sequence A;
sampling characteristic currents with a period of T and code bits of 1 according to the same sampling frequency f to obtain a sequence a;
and performing cross-correlation operation on the sequences A and a according to the following formula to obtain a sequence C:
wherein C (i) is an element of sequence C, n a Number of sampling points of sequence a, n a ＝T*f。
Further, the operation simplification processing module performs operation simplification processing on the sequence C, and includes:
for sequence C, start at point 1 and every n s Selecting the largest point as an element of the sequence D to form the sequence D;
wherein n is s The number of points occupied by a 50Hz sine wave in each sampling period, n s ＝f/50。
Further, the amplification smoothing processing module performs amplification smoothing processing on the sequence D, and includes:
constructing a sequence e, wherein the sequence e comprises T x f elements, each element is 1, and the number of the sequence e is n e ；
And performing cross-correlation operation on the sequence D and the sequence e according to the following formula to obtain a sequence F:
wherein F (i) is an element in the sequence F, and the number n of F (i) in the sequence F F ＝n D -n e 。
Further, the code bit identification module performs code bit identification on the sequence F, and includes:
acquiring the median of the sequence F, and recording as peak _ threshold;
let i be the index of sequence F, find the first occurrence of F (i +1) -F (i)<0, and F (i) -F (i-1)>0, and F (i)>The index of peak _ threshold is recorded as i ═ a; wherein i is 1,2,3 … … … n F ；
Finding out i with the maximum value of F (i) from the index i ═ a to i ═ a + one _ cycle _ num, and marking the index i ═ begin _ index; wherein, one _ cycle _ num is T/0.02;
acquiring the total code number all _ num of the sequence F; wherein all _ num is T/T;
let j be 1 and let j be 1,
when begin _ index + (j-1) one _ cycle _ num is not more than n F When the temperature of the water is higher than the set temperature,
f (begin _ index + (j-1) × one _ cycle _ num) > peak _ thermoshold, then the code bit of sequence F (i) there is 1, otherwise 0;
when begin _ index + (j-1) one _ cycle _ num>n F When the temperature of the water is higher than the set temperature,
F(n F )>peak thermally, the code bit of sequence f (i) is 1, otherwise it is 0;
if j is less than all + num, j is j +1, and the steps are repeatedly executed; otherwise, ending.
In summary, the present invention provides a method and an apparatus for identifying a code bit of a feature current after noise addition, which can identify a code bit of a feature current added with noise such as white gaussian noise, impulse noise, user power consumption, etc. according to the identification method provided by the present invention, the feature current after noise addition is processed through the steps of cross-correlation operation, simplified operation, amplified smoothing, code bit identification, etc., so as to identify a peak in a new signal, and further identify a code bit of the feature current after noise addition. The cross-correlation operation is to perform cross-correlation operation on the signal subjected to noise addition and a locally recorded sequence to obtain a new signal, so that code bit identification of the characteristic current can be realized without adopting a digital filter. Compared with the recognition method in the prior art, the code bit recognition method provided by the invention can recognize the code bit of the characteristic current after noise addition without using a digital filter, so that the calculation amount required by the algorithm is reduced, and the complexity of the method is reduced because the algorithm does not need a filter, and the algorithm can be still used under the condition that the digital filtering is not ideal, so that the algorithm is simpler to use and has a wider application range.
Drawings
FIG. 1 is a flow chart of a method for recognizing a code bit of a noise-added characteristic current according to the present invention;
FIG. 2 is a waveform diagram of an uncoded characteristic current, FIG. 2a is a waveform diagram of a characteristic current, and FIG. 2b is a code bit diagram of a characteristic current;
FIG. 3 is a waveform diagram of a characteristic current after noise addition;
FIG. 4 is a flowchart of an enlargement smoothing process;
FIG. 5 is a cross-correlation and magnification smoothing comparison graph;
FIG. 6 is a flow diagram of a code bit identification process;
fig. 7 is a block diagram showing the structure of the device for recognizing a code bit of a characteristic current after noise addition according to the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the accompanying drawings in conjunction with the following detailed description. It should be understood that the description is intended to be exemplary only, and is not intended to limit the scope of the present invention. Moreover, in the following description, descriptions of well-known structures and techniques are omitted so as to not unnecessarily obscure the concepts of the present invention.
The technical solution of the present invention will be described in detail below with reference to the accompanying drawings. According to an embodiment of the present invention, a method for identifying code bits of a noisy characteristic current is provided, where the method identifies peaks in an obtained new signal by performing cross-correlation operation, operation simplification processing, amplification smoothing processing, code bit identification, and other process steps on the noisy characteristic current, and further identifies code bits of the noisy characteristic current. The flow chart of the identification method of the embodiment is shown in fig. 1, and includes the steps of:
performing cross-correlation operation on the characteristic current subjected to noise addition to obtain a sequence C;
carrying out operation simplification processing on the sequence C to obtain a sequence D;
carrying out amplification smoothing treatment on the sequence D to obtain a sequence F;
and carrying out code bit recognition on the sequence F to obtain a code bit recognition result.
The above-described steps of the method will be specifically described below. The identification method is used for identifying the code bit of the characteristic current after noise addition, wherein the noise addition comprises noise generated by Gaussian white noise, pulse noise, household user electricity consumption and the like.
Firstly, cross-correlation operation is carried out on the characteristic current after noise addition to obtain a sequence C. The cross-correlation operation is to perform cross-correlation operation on the noisy signal and the locally recorded sequence to obtain a new signal. After receiving the current sequence, the receiver first performs cross-correlation operation on the current sequence, and fig. 2 shows a waveform diagram of an un-noise characteristic current, where a waveform in fig. 2a is a waveform diagram of the characteristic current, and fig. 2b is a code bit diagram of the characteristic current, and a signal in a period in which a code bit in fig. 2b is 1 is used to perform sliding point multiplication on the noise-added signal, so as to obtain a new signal. The step of cross-correlation may be performed as follows:
sampling the characteristic current with the duration of t after the noise is added by using a sampling frequency f to obtain a sequence A;
sampling characteristic currents with a period of T and code bits of 1 according to the same sampling frequency f to obtain a sequence a;
and performing cross-correlation operation on the sequences A and a according to the following formula to obtain a sequence C:
wherein C (i) is the i-th element of the sequence C, n a Number of sampling points of sequence a, n a T × f. i denotes the ith element of the sequence C, the total number of which is related to the number of elements of the sequence A, n C ＝n A –n a Wherein n is C Is the number of elements of the sequence C, n A Is the number of elements of sequence A, n A ＝t*f，n a T is the time period required to transmit one code bit.
As shown in fig. 2, the waveform of the characteristic current without noise is shown in fig. 2a, and the waveform of the characteristic current is shown in fig. 2a, and the waveform is formed on the alternating current of the power grid at the switching frequency of 833.3 Hz; FIG. 2b is a diagram of code bits for the characteristic current, with the switch turned off when the code bit is 0; when the code bit is 1, the characteristic current waveform in fig. 2a is formed, and the period of the code is 0.6 s.
Fig. 3 shows the characteristic current of the code bits plus gaussian white noise, impulse noise and user current noise, which is the current sequence sampled by the receiver and shown in fig. 3, with a sampling frequency of 5 kHz.
The current amplitude is relatively large due to the power consumption of a plurality of users, and the maximum amplitude is about 250A at the moment, so that the signal of 0.06A of the characteristic current is totally submerged in noise. The core idea of the method is that a filter is not used for filtering the influence of the 50Hz household electricity, and when the code bit of the characteristic current is 1, although the change trend of the code bit information of the sampled current after noise addition is not obvious due to the influence of noise, the trend is still contained in a current curve. When the characteristic current code bit is 1, from a single sampling point, the influence of noise is large, the overall change trend of a plurality of sampling points with the current code bit of 1 is more accurate than the characteristic of the single sampling point, so that the amplitude of the sequence with 1 can be selected, the number of the sequence points is the same as the number of the sampling points of the current code bit, and the sequence is subjected to cross-correlation operation with the original sequence, so that the change trend is amplified, and the noise interference is resisted. The method is specifically described below:
and (4) carrying out cross-correlation operation on the unfiltered signal to obtain a new signal, and carrying out amplification smoothing on the signal. Before this processing, simple arithmetic simplification processing can be performed. And simplifying the operation, namely selecting the maximum power in each sine wave period to represent the characteristic of each sine wave period, and obtaining the sequence D after processing.
After a sequence C is obtained through cross-correlation operation, operation simplification processing is carried out on the sequence C to obtain a sequence D, the operation simplification processing is that a part of signals to be processed are extracted according to rules to form new signals, and the operation simplification processing can be carried out by adopting the following steps:
for sequence C, start at point 1 and every n s Selecting the largest point as an element of the sequence D to form the sequence D;
wherein n is s Number of points occupied by 50Hz sine wave in each sampling period, n s ＝f/50。
And carrying out amplification smoothing treatment on the sequence D to obtain a sequence F. The amplification smoothing process is to amplify some rules hidden in the signal and smooth the comparatively jittered signal to obtain a new signal for the next identification step. Fig. 4 shows a flow chart of the steps of the amplification smoothing process, which mainly has two steps, the first step is to construct a new sequence e, the number of points of the sequence is the same as the number of samples whose code bits are 1, and each value of the sequence is 1; and secondly, performing sliding point multiplication by using the new sequence e and the sequence D to obtain a new sequence F. The purpose of the amplification smoothing processing step is to mine code bit information that cannot be reflected from the cross-correlation sequence, so that the hidden information is amplified and can be normally identified. A cross-correlation and amplification smoothing process comparison is shown in fig. 5. The specific steps of the amplification smoothing process may be as follows:
constructing a sequence e, wherein the sequence e comprises T x f elements, each element is 1, and the number of the sequence e is n e ；
And performing cross-correlation operation on the sequence D and the sequence e according to the following formula to obtain a sequence F:
wherein F (i) is an element in the sequence F, and the number n of F (i) in the sequence F F ＝n D -n e ，n D The number of elements representing sequence D may be n D And (3) calculating:
every n from sequence C s A sampling point (n) s F/50), from which the largest point is selected to form the new sequence D, n D ＝n C /n s 。
After the sequence F of recognizable peaks is obtained by the amplification smoothing process, the code bits of the sequence F may be recognized, which may refer to the code bit determination process shown in fig. 6. Firstly, some characteristic thresholds in the sequence F are obtained, then, the sequence index of the beginning of the code bit 1 is roughly positioned, then, the point with the maximum amplitude in the next code bit period is found according to the index, the point is taken as the middle point of the code bit 1, finally, the middle point of each code bit after the sequence is deduced according to the point, and the code bit condition is judged by utilizing the thresholds. The code bit identification can be specifically carried out by adopting the following steps:
acquiring the median of the sequence F, and recording as peak _ threshold;
let m be the index of sequence F, find the first occurrence of F (m +1) -F (m)<0, and F (m) -F (m-1)>0, and F (m)>The index of peak _ threshold, and the corresponding index m is recorded as m ═ a; wherein m is 1,2,3 … … … n F ；
Finding out an index value m corresponding to the maximum value of F (m) from the index m ═ a to m ═ a + one _ cycle _ num, and marking the index m ═ begin _ index; wherein, one _ cycle _ num is T/0.02;
acquiring the total code number all _ num of the sequence F; wherein all _ num is T/T;
let j be 1 and let j be 1,
when begin _ index + (j-1) one _ cycle _ num is not more than n F When the temperature of the water is higher than the set temperature,
f (begin _ index + (j-1) one _ cycle _ num) > peak _ thermoshiftd, then the code bit of sequence F (m) at that point is 1, otherwise 0;
when begin _ index + (j-1) one _ cycle _ num>n F When the temperature of the water is higher than the set temperature,
F(n F )>peak _ thermoshold, where the code bit of sequence f (m) is 1, otherwise 0;
if j is less than all + num, j is j +1, and the steps are repeatedly executed to search again; otherwise, the code bit identification process is finished.
According to another embodiment of the present invention, there is provided an apparatus for recognizing code bits of a noisy characteristic current, the apparatus being configured as shown in fig. 7 and including a cross-correlation operation module, an operation simplification processing module, an amplification smoothing processing module, and a code bit recognition module. The modules are connected in sequence to identify code bits of the characteristic current after noise is added. The functions of the respective modules are explained below.
And the cross-correlation operation module is used for performing cross-correlation operation on the characteristic current subjected to noise addition to obtain a sequence C: sampling the characteristic current with the duration of t after the noise is added by using a sampling frequency f to obtain a sequence A;
sampling characteristic currents with a period of T and code bits of 1 according to the same sampling frequency f to obtain a sequence a;
and performing cross-correlation operation on the sequences A and a according to the following formula to obtain a sequence C:
wherein C (i) is an element of the sequence C, n a Number of sampling points of sequence a, n a ＝T*f。
And the operation simplification processing module is used for carrying out operation simplification processing on the sequence C to obtain a sequence D: for sequence C, start at point 1 and every n s Selecting the largest point as an element of the sequence D to form the sequence D;
wherein n is s Number of points occupied by 50Hz sine wave in each sampling period, n s ＝50/f。
And the amplification smoothing processing module is used for carrying out amplification smoothing processing on the sequence D to obtain a sequence F: constructing a sequence e, wherein the sequence e comprises T x f elements, each element is 1, and the number of the sequence e is n e ；
And performing cross-correlation operation on the sequence D and the sequence e according to the following formula to obtain a sequence F:
wherein F (i) is an element in the sequence F, and the number n of F (i) in the sequence F F ＝n D -n e 。
And the code bit identification module is used for identifying the code bit of the sequence F to obtain a code bit identification result:
acquiring the median of the sequence F, and recording as peak _ threshold;
let m be the index of sequence F, find the first occurrence of F (m +1) -F (m)<0, and F (m) -F (m-1)>0, and F (m)>The index of peak _ threshold, denoted as m ═ a; wherein m is 1,2,3 … … n F ；
Finding out m when the value of F (m) is maximum from the index m ═ a to m ═ a + one _ cycle _ num, and marking the index m ═ begin _ index; wherein, one _ cycle _ num is T/0.02;
acquiring the total code number all _ num of the sequence F; wherein all _ num is T/T;
let j be 1 and let j be 1,
when begin _ index + (j-1) one _ cycle _ num is not more than n F When the temperature of the water is higher than the set temperature,
f (begin _ index + (j-1) one _ cycle _ num) > peak _ thermoshiftd, then the code bit of sequence F (i) at that point is 1, otherwise 0;
when begin _ index + (j-1) one _ cycle _ num>n F When the temperature of the water is higher than the set temperature,
F(n F )>peak thermally, the code bit of sequence f (m) is 1, otherwise 0;
if j is less than all + num, j is j +1, and the steps are repeatedly executed; otherwise, the code bit identification process is finished.
In summary, the present invention relates to a method and an apparatus for identifying a code bit of a noise-added characteristic current, which can identify a code bit of a characteristic current to which noise such as white gaussian noise, impulse noise, and user power consumption is added. The cross-correlation operation is to perform cross-correlation operation on the signal subjected to noise addition and a locally recorded sequence to obtain a new signal, so that code bit identification of the characteristic current can be realized without adopting a digital filter. Compared with the recognition method in the prior art, the code bit recognition method provided by the invention can recognize the code bit of the characteristic current after the noise is added without using a digital filter, thereby reducing the calculation amount required by the algorithm.
It is to be understood that the above-described embodiments of the present invention are merely illustrative of or explaining the principles of the invention and are not to be construed as limiting the invention. Therefore, any modification, equivalent replacement, improvement and the like made without departing from the spirit and scope of the present invention should be included in the protection scope of the present invention. Further, it is intended that the appended claims cover all such variations and modifications as fall within the scope and boundaries of the appended claims or the equivalents of such scope and boundaries.
Claims (10)
1. A method for recognizing a code bit of a characteristic current after noise addition is characterized by comprising the following steps:
performing cross-correlation operation on the characteristic current subjected to noise addition to obtain a sequence C;
carrying out operation simplification processing on the sequence C to obtain a sequence D;
carrying out amplification smoothing treatment on the sequence D to obtain a sequence F;
and carrying out code bit recognition on the sequence F to obtain a code bit recognition result.
2. The method of claim 1, wherein performing a cross-correlation operation on the noisy characteristic current comprises:
sampling the characteristic current with the duration of t after the noise is added by using a sampling frequency f to obtain a sequence A;
sampling characteristic currents with a period of T and code bits of 1 according to the same sampling frequency f to obtain a sequence a;
and performing cross-correlation operation on the sequences A and a according to the following formula to obtain a sequence C:
wherein C (i) is the ith element of the sequence C, n a Number of sampling points of sequence a, n a ＝T*f。
3. The method of claim 2, wherein the performing operation simplification processing on the sequence C comprises:
for sequence C, start at point 1 and every n s Selecting the largest point as an element of the sequence D to form the sequence D;
wherein n is s Number of sampling points occupied by 50Hz sine wave in each period, n s ＝f/50。
4. The method according to claim 3, wherein the performing amplification smoothing on the sequence D comprises:
constructing a sequence e, wherein the sequence e comprises T x f elements, each element is 1, and the number of the elements of the sequence e is n e ；
And performing cross-correlation operation on the sequence D and the sequence e according to the following formula to obtain a sequence F:
wherein F (i) is an element in the sequence F, and the number n of F (i) in the sequence F F ＝n D -n e ，n D Is the number of elements in sequence D.
5. The method of claim 4, wherein the code bit identification of the sequence F comprises:
acquiring the median of the sequence F, and recording as peak _ threshold;
let m be the index of sequence F, find the first occurrence of F (m +1) -F (m)<0, and F (m) -F (m-1)>0, and F (m)>The index of peak _ threshold, denoted as m ═ a; wherein m is 1,2,3 … … … n F ；
Finding out m when the value of F (m) is maximum from the index m ═ a to m ═ a + one _ cycle _ num, and marking the index m ═ begin _ index; wherein, one _ cycle _ num is T/0.02;
acquiring the total code number all _ num of the sequence F; wherein all _ num is T/T;
let j be 1 and let j be 1,
when begin _ index + (j-1) one _ cycle _ num is not more than n F When the utility model is used, the water is discharged,
f (begin _ index + (j-1) one _ cycle _ num) > peak _ thermoshiftd, then the code bit of sequence F (m) at that point is 1, otherwise 0;
when begin _ index + (j-1) one _ cycle _ num>n F When the utility model is used, the water is discharged,
F(n F )>peak thermally, the code bit of sequence f (m) is 1, otherwise 0;
if j is less than all + num, j is j +1, and the steps are repeatedly executed; otherwise, ending.
6. A recognition device of a characteristic current code bit after noise addition is characterized by comprising a cross-correlation operation module, an operation simplification processing module, an amplification smoothing processing module and a code bit recognition module; wherein,
the cross-correlation operation module is used for performing cross-correlation operation on the characteristic current subjected to noise addition to obtain a sequence C;
the operation simplification processing module is used for carrying out operation simplification processing on the sequence C to obtain a sequence D;
the amplification smoothing processing module is used for carrying out amplification smoothing processing on the sequence D to obtain a sequence F;
and the code bit identification module is used for carrying out code bit identification on the sequence F to obtain a code bit identification result.
7. The apparatus of claim 6, wherein the cross-correlation module performs a cross-correlation operation on the noisy characteristic current, and comprises:
sampling the characteristic current with the duration of t after the noise is added by using a sampling frequency f to obtain a sequence A;
sampling characteristic currents with a period of T and code bits of 1 according to the same sampling frequency f to obtain a sequence a;
and performing cross-correlation operation on the sequences A and a according to the following formula to obtain a sequence C:
wherein C (i) is an element of the sequence C, n a Number of sampling points of sequence a, n a ＝T*f。
8. The apparatus of claim 7, wherein the operation simplification processing module performs operation simplification processing on the sequence C, and comprises:
for sequence C, start at point 1 and every n s Selecting the largest point as an element of the sequence D to form the sequence D;
wherein n is s The number of points occupied by a 50Hz sine wave in each sampling period, n s ＝f/50。
9. The apparatus of claim 8, wherein the amplification smoothing module performs amplification smoothing on the sequence D, and comprises:
constructing a sequence e, wherein the sequence e comprises T x f elements, each element is 1, and the number of the sequence e is n e ；
And performing cross-correlation operation on the sequence D and the sequence e according to the following formula to obtain a sequence F:
wherein F (i) is an element in the sequence F, and the number n of F (i) in the sequence F F ＝n D -n e 。
10. The apparatus of claim 9, wherein the code bit identification module is configured to perform code bit identification on a sequence F, and comprises:
acquiring a median of the sequence F, and recording the median as peak _ threshold;
let i be the index of sequence F, find the first occurrence of F (m +1) -F (m)<0, and F (m) -F (m-1)>0, and F (m)>The index of peak _ threshold, denoted as m ═ a; wherein m is 1,2,3 … … … n F ；
Finding out m when the value of F (m) is maximum from the index m ═ a to m ═ a + one _ cycle _ num, and marking the index m ═ begin _ index; wherein, one _ cycle _ num is T/0.02;
acquiring the total code number all _ num of the sequence F; wherein all _ num is T/T;
let j be 1 and let j be 1,
when begin _ index + (j-1) one _ cycle _ num is not more than n F When the temperature of the water is higher than the set temperature,
f (begin _ index + (j-1) × one _ cycle _ num) > peak _ thermoshold, then the code bit of sequence F (m) there is 1, otherwise 0;
when begin _ index + (j-1) one _ cycle _ num>n F When the temperature of the water is higher than the set temperature,
F(n F )>peak thermally, the code bit of sequence f (m) is 1, otherwise 0;
if j is less than all + num, j is j +1, and the steps are repeatedly executed; otherwise, ending.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110294524.6A CN115112964B (en) | 2021-03-19 | 2021-03-19 | Method and device for recognizing characteristic current code bits after noise addition |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110294524.6A CN115112964B (en) | 2021-03-19 | 2021-03-19 | Method and device for recognizing characteristic current code bits after noise addition |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115112964A true CN115112964A (en) | 2022-09-27 |
| CN115112964B CN115112964B (en) | 2024-08-09 |
Family
ID=83322783
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110294524.6A Active CN115112964B (en) | 2021-03-19 | 2021-03-19 | Method and device for recognizing characteristic current code bits after noise addition |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115112964B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150316592A1 (en) * | 2013-01-21 | 2015-11-05 | Mitsubishi Electric Corporation | Electric device identification system, electric device, and electric device identification method |
| CN108847653A (en) * | 2018-07-12 | 2018-11-20 | 广东电网有限责任公司 | Method and device for identifying magnetizing inrush current of transformer based on waveform change rate |
| US10330812B1 (en) * | 2015-03-16 | 2019-06-25 | Elintrix | Identifying schematic and topological properties of an electrical grid through analysis of directional relationships of electrical currents |
| CN111600380A (en) * | 2020-03-28 | 2020-08-28 | 青岛鼎信通讯股份有限公司 | Variable load and constant current detection-based household variable relation identification method |
| CN112019241A (en) * | 2020-10-22 | 2020-12-01 | 中国电力科学研究院有限公司 | Method and system for identifying equipment association relation in power system |
| CN112213578A (en) * | 2020-09-23 | 2021-01-12 | 青岛鼎信通讯股份有限公司 | Medium voltage line variable relation identification method |
-
2021
- 2021-03-19 CN CN202110294524.6A patent/CN115112964B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20150316592A1 (en) * | 2013-01-21 | 2015-11-05 | Mitsubishi Electric Corporation | Electric device identification system, electric device, and electric device identification method |
| US10330812B1 (en) * | 2015-03-16 | 2019-06-25 | Elintrix | Identifying schematic and topological properties of an electrical grid through analysis of directional relationships of electrical currents |
| CN108847653A (en) * | 2018-07-12 | 2018-11-20 | 广东电网有限责任公司 | Method and device for identifying magnetizing inrush current of transformer based on waveform change rate |
| CN111600380A (en) * | 2020-03-28 | 2020-08-28 | 青岛鼎信通讯股份有限公司 | Variable load and constant current detection-based household variable relation identification method |
| CN112213578A (en) * | 2020-09-23 | 2021-01-12 | 青岛鼎信通讯股份有限公司 | Medium voltage line variable relation identification method |
| CN112019241A (en) * | 2020-10-22 | 2020-12-01 | 中国电力科学研究院有限公司 | Method and system for identifying equipment association relation in power system |
Non-Patent Citations (1)
| Title |
|---|
| 王丙元;林晓晶;李宝胜;: "基于数学形态学的工频通信信号识别技术", 中国民航大学学报, no. 01, 15 February 2008 (2008-02-15) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115112964B (en) | 2024-08-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109890043B (en) | Wireless signal noise reduction method based on generative countermeasure network | |
| CN106885971B (en) | Intelligent background noise reduction method for cable fault detection pointing instrument | |
| CN114093377B (en) | Splitting normalization method and device, audio feature extractor and chip | |
| CN110690931B (en) | Digital signal adaptive code rate estimation method and device based on multi-wavelet-base combination | |
| EP3206574A1 (en) | Frame based spike detection module | |
| CN112309417A (en) | Wind noise suppression audio signal processing method, device, system and readable medium | |
| CN113473522A (en) | Method for detecting synchronization information in cell search, related device and storage medium | |
| CN108039182B (en) | Voice activation detection method | |
| CN115112964A (en) | Method and device for recognizing characteristic current code bits after noise addition | |
| CN107465639B (en) | Multi-path delay synchronous decision demodulation method based on short-time discrete Fourier transform | |
| CN107507630A (en) | A kind of non-cooperation voice communication receives data dead time section recognition methods | |
| CN108737319B (en) | Real-time detection method and device for target OFDM signal | |
| CN101141143B (en) | Time-frequency RAKE receiving method and apparatus | |
| CN105959035B (en) | A kind of direct sequence signal intercepts and captures detection method | |
| Wu et al. | Joint time delay and frequency estimation via state-space realization | |
| CN116112039A (en) | Unmanned aerial vehicle frequency hopping signal rapid detection method based on FPGA | |
| CN115616280A (en) | Event detection method based on attention mechanism | |
| An et al. | Spectrum sensing based on KNN algorithm for 230 MHz power private networks | |
| CN114301175A (en) | Power distribution station area user transformation relation identification method and device based on injection signals | |
| CN111064686B (en) | Symbol timing synchronization method and device | |
| CN111131100B (en) | NBIOT system multi-cell iterative interference cancellation detection method | |
| Li et al. | Goertzel Algorithm Based on Matlab Platform System Simulation for Listening Dial Tone Recognition Research | |
| Gu et al. | Exploiting ResNeXt with Convolutional Shortcut for Signal Modulation Classification at Low SNRs | |
| CN114759680B (en) | Method and system for identifying topology of station area | |
| CN114299938B (en) | Intelligent voice recognition method and system based on deep learning |
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 |