CN115378467A - Power line noise sample extraction method based on diversity signal cancellation - Google Patents

Power line noise sample extraction method based on diversity signal cancellation Download PDF

Info

Publication number
CN115378467A
CN115378467A CN202210866696.0A CN202210866696A CN115378467A CN 115378467 A CN115378467 A CN 115378467A CN 202210866696 A CN202210866696 A CN 202210866696A CN 115378467 A CN115378467 A CN 115378467A
Authority
CN
China
Prior art keywords
noise
power line
signal
channel
diversity
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
Application number
CN202210866696.0A
Other languages
Chinese (zh)
Other versions
CN115378467B (en
Inventor
陈智雄
张志坤
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202210866696.0A priority Critical patent/CN115378467B/en
Publication of CN115378467A publication Critical patent/CN115378467A/en
Application granted granted Critical
Publication of CN115378467B publication Critical patent/CN115378467B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/54Systems for transmission via power distribution lines
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/02Details
    • H04B3/46Monitoring; Testing
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B3/00Line transmission systems
    • H04B3/54Systems for transmission via power distribution lines
    • H04B3/544Setting up communications; Call and signalling arrangements
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L1/00Arrangements for detecting or preventing errors in the information received
    • H04L1/0001Systems modifying transmission characteristics according to link quality, e.g. power backoff
    • H04L1/0006Systems modifying transmission characteristics according to link quality, e.g. power backoff by adapting the transmission format
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00Modulated-carrier systems
    • H04L27/26Systems using multi-frequency codes
    • H04L27/2601Multicarrier modulation systems
    • H04L27/2647Arrangements specific to the receiver only
    • H04L27/2655Synchronisation arrangements
    • H04L27/2668Details of algorithms
    • H04L27/2681Details of algorithms characterised by constraints
    • H04L27/2688Resistance to perturbation, e.g. noise, interference or fading
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE 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/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Power Engineering (AREA)
  • Quality & Reliability (AREA)
  • Radio Transmission System (AREA)

Abstract

The invention provides a power line noise sample extraction method based on diversity signal cancellation. The method comprises the following steps: diversity signals to be transmitted are respectively transmitted in a power line and a wireless channel in a diversity mode, and the power line is influenced by impulse noise; respectively carrying out discrete Fourier transform on diversity signals received by receiving ends of a power line and a wireless channel, and then obtaining noise data in the power line signals through cancellation operation; the impulse noise data is removed from the power line signal received at the receiving end by a non-linear function based on an optimal threshold estimate of the noise samples. The invention provides a power line noise sample extraction method based on diversity signal cancellation, which aims at a dual-interface communication architecture. The noise sample is extracted by utilizing the consistency of the power line and the wireless channel transmission signal, the limitation that only the power line is used for independently processing the impulse noise in the traditional hybrid communication is broken through, and the method can be applied to the prediction of the optimal threshold value.

Description

Power line noise sample extraction method based on diversity signal cancellation
Technical Field
The invention relates to the technical field of smart power grids, in particular to a power line noise sample extraction method based on diversity signal cancellation.
Background
The comprehensive perception and communication of information are the basis and guarantee for constructing the intelligent power grid. The smart grid is wide in coverage range, complex in topological structure and diversified in access equipment, and a single communication technology cannot meet the requirements of various scenes, so that hybrid communication becomes a key technology of the smart grid. The Power Line Communication (PLC) coverage is wide, and no extra wiring is needed; wireless Communication (WLC) is flexible in access and strong in expandability, so that PLC and WLC technologies are closely combined, the system performance of hybrid Communication is improved, and the Wireless Communication network has a wide application prospect in smart power grids and internet of things, and has become one of research hotspots.
However, the power line channel is highly susceptible to Impulse Noise (IN) due to sudden voltage changes caused by the switching on and off of the household appliance and power electronics switches. In particular, the burst noise is generated at a plurality of consecutive sampling points, and affects not only demodulation and decoding of signals but also synchronization performance of power line transmission. How to suppress the impulse noise of the power line channel in the dual-interface communication and further improve the coverage area of the communication system and the reliability of the access thereof is one of the key problems to be solved urgently.
At present, in a dual-interface communication system under the influence of impulse noise in the prior art, researchers mostly focus on technologies such as theoretical performance and signal combination. The Gamma distribution of the wireless fading is approximated to Log Normal (LogN) distribution, so that the mixed fading problem in the power line and wireless parallel communication system is converted into the LogN-LogN same distribution fading problem, and then the theoretical performances such as the system error rate and the like are analyzed. The performances of algorithms such as maximum ratio combination, selective combination and the like in a power line and a wireless parallel communication system are researched by adopting a Middleton Class A pulse noise model. However, these studies often only consider the impact of impulse noise on system performance, and do not address how impulse noise is eliminated. According to the scheme, a multi-stage orthogonal matching tracking recovery algorithm is provided for different types of impulse noise and narrow-band interference, and is used for recovering the impulse noise and the narrow-band interference jointly, but the complexity of parallel processing of double-interface signals is high.
Compared with a compressed sensing algorithm, the impulse noise suppression algorithm based on the nonlinear function obtains attention of researchers due to the advantages of high calculation speed, no need of null sub-carriers and the like. The most widely used non-linear functions include nulling, clipping/deep clipping, and combinations thereof. However, in an Orthogonal Frequency Division Multiplexing (OFDM) system, due to the existence of a Peak-to-Average Power Ratio (PAPR), if an accurate noise threshold cannot be obtained, part of the useful signal will be treated as noise, which greatly limits the performance of the nonlinear function.
In recent years, attention has been paid to the determination of an optimum threshold value for a nonlinear function. In some schemes, burst noise is discretized through an interleaving technology, then a moment estimation is used for obtaining a noise statistical parameter and calculating a noise threshold, but the scheme ignores the influence of channel fading on a system, and the interleaving technology can increase the transmission delay of the system. It is not beneficial to the application of the algorithm in the delay sensitive network. For deterministic signals in impulse noise, the Local Optimal Detection (LOD) algorithm has the best Detection performance at low signal-to-noise ratio, and approaches the theoretical Optimal solution. However, the LOD algorithm needs to know the Probability Density Function (PDF) of the noise accurately, which may not be obtained in the actual scene. Some schemes combine clipping and time/frequency channel equalization, and perform combined clipping before OFDM signal demodulation, and then perform time/frequency equalization, but this method does not significantly improve communication quality, and the threshold is generally obtained through experience.
In the above solutions in the prior art, impulse noise is processed separately in the PLC, and a dual-interface communication environment of a power line and a wireless is not involved.
Disclosure of Invention
The embodiment of the invention provides a power line noise sample extraction method based on diversity signal cancellation, so as to effectively extract noise samples of diversity signals in a power line channel.
In order to achieve the purpose, the invention adopts the following technical scheme.
A power line noise sample extraction method based on diversity signal cancellation comprises the following steps:
diversity signals to be transmitted are diversity-transmitted in a power line and a wireless channel, respectively, the power line being affected by impulse noise;
respectively carrying out discrete Fourier transform on the diversity signals received by the receiving ends of the power line channel and the wireless channel, and then obtaining noise data in the power line diversity signals through cancellation operation;
the impulse noise data is cancelled from the power line diversity signal received at the receiving end by a non-linear function based on an optimal threshold estimate of the noise samples.
Preferably, the diversity signals to be transmitted are diversity-transmitted in a power line and a wireless channel, respectively, the power line being affected by impulse noise, and the diversity signals include:
diversity signal to be transmitted
Figure BDA0003759573950000034
Respectively carrying out diversity transmission in a power line and a wireless channel, wherein the wireless channel is influenced by Rayleigh fading and white Gaussian noise, the power line channel is influenced by logarithmically normal fading and pulse noise, and the diversity signals are respectively subjected to Orthogonal Frequency Division Multiplexing (OFDM) modulation in the power line and a wireless communication module;
radio channel fading coefficient h W Satisfy Rayleigh distribution, then h W The probability density function PDF of (1) is:
Figure BDA0003759573950000031
wherein
Figure BDA0003759573950000032
Is the rayleigh fading variance;
power line channel fading coefficient h P Satisfy LogN distribution, then h P The probability density function PDF of (1) is:
Figure BDA0003759573950000033
in the formula mu P And
Figure BDA0003759573950000041
are each lnh P Normalizing the channel fading energy to obtain the mean value and the variance of the channel fading energy
Figure BDA0003759573950000042
Preferably, the obtaining of the impulse noise data in the power line diversity signal by performing discrete fourier transform on the diversity signals received by the receiving ends of the power line channel and the wireless channel, and then performing cancellation operation includes:
the OFDM modulated time domain symbol of the t-th transmission time slot is
Figure BDA0003759573950000043
Diversity signal to be transmitted by a transmitting end
Figure BDA0003759573950000044
Respectively in the channel fading coefficient of
Figure BDA00037595739500000431
The power line channel and the channel fading coefficient are
Figure BDA0003759573950000045
In a wireless channel, a power line channel and a receiving end of the wireless channel receive time domain OFDM sampling signals
Figure BDA0003759573950000046
And
Figure BDA0003759573950000047
expressed as:
Figure BDA0003759573950000048
Figure BDA0003759573950000049
wherein the symbols
Figure BDA00037595739500000410
Represents a convolution operation;
Figure BDA00037595739500000411
and
Figure BDA00037595739500000412
respectively representing time domain Gaussian noise vectors on a power line and a wireless channel, wherein elements in the vectors respectively satisfy the conditions that the mean value is zero and the variance is
Figure BDA00037595739500000413
And
Figure BDA00037595739500000414
(ii) a gaussian distribution of;
Figure BDA00037595739500000415
representing a time-domain impulse noise vector whose elements satisfy a mean of zero and a variance of
Figure BDA00037595739500000416
Is a Gaussian distribution of
Figure BDA00037595739500000417
Greater than 1;
will be provided with
Figure BDA00037595739500000418
And
Figure BDA00037595739500000419
after the cyclic prefix is removed, discrete Fourier Transform (DFT) is carried out to obtain equivalent frequency domain OFDM signals of the power line channel and the wireless channel
Figure BDA00037595739500000420
And
Figure BDA00037595739500000421
respectively expressed as:
Figure BDA00037595739500000422
Figure BDA00037595739500000423
wherein F represents the DFT operator and wherein,
Figure BDA00037595739500000424
denotes an OFDM frequency domain symbol vector transmitted by the transmitting end in the t-th transmission slot, N is the number of subcarriers, N] T Which represents the operation of transposition by means of a transposition operation,
Figure BDA00037595739500000425
and
Figure BDA00037595739500000426
equivalent channel coefficient matrices respectively representing a power line and a radio channel of a t-th transmission slot,
Figure BDA00037595739500000427
and
Figure BDA00037595739500000428
are all diagonal matrices, and
Figure BDA00037595739500000429
is reversible, make
Figure BDA00037595739500000430
The left-hand matrix Q can be obtained by:
Figure BDA0003759573950000051
wherein
Figure BDA0003759573950000052
Figure BDA0003759573950000053
To represent
Figure BDA0003759573950000054
The inverse matrix of (5) and (7) is subtracted to obtain the same signal part in the power line channel receiving end which is cancelled out of the wireless channel
Figure BDA0003759573950000055
The post frequency domain noise estimate vector Ψ t Comprises the following steps:
Figure BDA0003759573950000056
wherein
Figure BDA0003759573950000057
Is a frequency domain vector of the impulse noise of the time domain after DFT transformation,
Figure BDA0003759573950000058
representing a background noise interference term of the dual interface channel;
let F * Representing the inverse discrete Fourier transform operator, then for Ψ t Performing an inverse discrete Fourier transform F * Ψ t Obtaining a time-domain noise signal vector
Figure BDA0003759573950000059
Comprises the following steps:
Figure BDA00037595739500000510
wherein
Figure BDA00037595739500000511
And
Figure BDA00037595739500000512
subject to a Gaussian distribution with a mean value of 0, given a known Q, will
Figure BDA00037595739500000513
Approximate mean of 0 and variance of
Figure BDA00037595739500000514
A gaussian distribution of (a).
Preferably, the removing the impulse noise data from the power line diversity signal received by the receiving end by the nonlinear function based on the optimal threshold estimation of the noise sample comprises:
determining threshold T using a data sample based adaptive threshold estimation method t Using a zero-setting algorithm in a non-linear function to reduce the amplitude below a threshold T t Is replaced by zero, the reserved amplitude is greater than or equal to T t Of the noise signal
Figure BDA00037595739500000515
For the extracted noise time domain signal vector
Figure BDA00037595739500000516
M-th noise element of (1)
Figure BDA00037595739500000517
The zero-setting transformation is adopted:
Figure BDA00037595739500000518
where M represents the total number of samples in the tth transmission slot, |, represents the absolute value operation,
Figure BDA00037595739500000519
a nonlinear transformation result representing impulse noise at the mth sampling point of the tth transmission time slot;
order vector
Figure BDA00037595739500000520
To represent
Figure BDA00037595739500000521
Is a result of a non-linear transformation of
Figure BDA00037595739500000522
For power line signal
Figure BDA00037595739500000523
Eliminating impulse noise to obtain power line signal without impulse noise
Figure BDA00037595739500000524
Comprises the following steps:
Figure BDA0003759573950000061
like equations (5) and (6), F represents a discrete fourier transform operator.
Preferably, the method further comprises:
will be provided with
Figure BDA0003759573950000062
With signals received by the wireless port
Figure BDA0003759573950000063
Performing merging treatment to obtain
Figure BDA0003759573950000064
Estimated symbol obtained by maximum likelihood detection algorithm for kth subcarrier of t transmission time slot
Figure BDA0003759573950000065
Is composed of
Figure BDA0003759573950000066
Wherein | · | purple 2 Which represents the 2-norm operation of the signal,
Figure BDA0003759573950000067
represents an equivalent received signal obtained by combining the power line signal and the wireless signal on the kth subcarrier of the t-th time slot, wherein omega represents a modulation signal constellation point set,
Figure BDA0003759573950000068
a symbol on a constellation diagram is represented,
Figure BDA0003759573950000069
which represents the equivalent channel coefficient at the kth sub-carrier of the t-th transmission slot at the receiving end, N is the number of sub-carriers,
Figure BDA00037595739500000610
meaning that the best x is chosen so that the function f (x) on x is minimal.
Preferably, the threshold value T is determined by adopting an adaptive threshold value estimation method based on data samples t The method comprises the following steps:
minimizing the bit error rate: let Pr (y) t (T t ) Means that the threshold value T is used in one OFDM data frame in the tth transmission slot t The obtained OFDM bit error rate is expressed as:
Figure BDA00037595739500000611
wherein | · | charging 0 Representing a 0-norm operation (i.e. counting the number of non-zero elements),
Figure BDA00037595739500000612
indicating the symbol sent by the sending end on the kth subcarrier of the tth transmission slot, where N is the number of subcarriers of one OFDM data frame, the objective function is:
Figure BDA00037595739500000613
wherein
Figure BDA00037595739500000614
Indicating a selected noise threshold
Figure BDA00037595739500000615
So as to relate to T t Error rate Pr (y) t (T t ) Minimum threshold, s.t T. t ≧ 0 indicates satisfaction of the "noise threshold T t The selection range of the receiving terminal is greater than or equal to 0 ″, and the receiving terminal obtains the signal at the t-th transmission time slot
Figure BDA0003759573950000071
And a channel coefficient matrix
Figure BDA0003759573950000072
And
Figure BDA0003759573950000073
based on the construction/update of the real-time noise sample, the error rate problem of the minimum current time is converted into the average error rate problem of the minimum noise sample, and the optimal threshold value of the current time is obtained
Figure BDA0003759573950000074
The expression is as follows:
Figure BDA0003759573950000075
wherein L is D Representing the magnitude of the number of noise samples, pr (y) i (T t ) Means that the ith OFDM data frame in the noise sample is related to the noise threshold T t S.t is an abbreviation of subject to, indicating that the constraint condition is satisfied;
when the threshold value is updated in the t-th transmission time slot, the weight lambda of each time slot in the data sample is calculated according to the discount factor gamma i ,i=1,2,…,L D According to the threshold value T at the T-1 th moment t-1 Combining the weight λ i For data sampleThis D EN Calculating the error rate, and obtaining the threshold value T by gradient descent method t-1 Corresponding error rate gradient value
Figure BDA0003759573950000076
Using learning rate l R Updating the threshold value of the current time
Figure BDA0003759573950000077
Figure BDA0003759573950000078
The solution process of (2) is shown as equation (16):
Figure BDA0003759573950000079
wherein
Figure BDA00037595739500000710
Denotes the value of x with respect to f (x) 0 The value of the gradient of (a) is,
Figure BDA00037595739500000711
is indicated at a noise threshold T t-1 Then, on the kth subcarrier of the ith OFDM data frame in the noise sample, combining the power line signal and the wireless signal to obtain an equivalent received signal,
Figure BDA00037595739500000712
indicating the symbol transmitted by the transmitting end on the k subcarrier of the ith OFDM data frame among the noise samples,
Figure BDA00037595739500000713
indicating the number of errors in the ith OFDM data frame in the noise sample,
Figure BDA00037595739500000714
then the weight of the ith OFDM data frame in the noise sample in the data sample is represented, i.e.
Figure BDA00037595739500000715
L D N represents the total number of symbols in the data sample, hence
Figure BDA0003759573950000081
Is the average bit error rate after weighted summation based on data samples.
It can be seen from the technical solutions provided by the embodiments of the present invention that, the method of the present invention provides a power line noise sample extraction method based on diversity signal cancellation for a dual-interface communication architecture. The noise sample is extracted by utilizing the consistency of the power line channel and the wireless channel transmission signal, the limitation that only the power line is used for independently processing the impulse noise in the traditional hybrid communication is broken through, and the method can be applied to the prediction of the optimal threshold value.
Additional aspects and advantages of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a processing flow chart of a power line noise sample extraction method based on diversity signal cancellation according to an embodiment of the present invention;
fig. 2 is a schematic diagram illustrating a markov chain of a TSMG noise model according to an embodiment of the present invention;
fig. 3 is a schematic diagram of a frame of a power line and wireless dual-interface communication system for performing signal processing at a receiving end according to an embodiment of the present invention;
FIG. 4 shows a probability density according to an embodiment of the present invention
Figure BDA0003759573950000082
The statistical histogram of (1);
fig. 5 shows a noise data extraction effect based on dual-interface diversity signal cancellation according to an embodiment of the present invention; fig. 5 (a) is a waveform diagram of a power line channel time domain signal, fig. 5 (b) is a waveform diagram of a wireless channel time domain signal, and fig. 5 (c) is a comparison of an extracted signal with an actual IN;
FIG. 6 shows Bit Error Rate (BER) and η @ under the same SNR according to an embodiment of the present invention WCC And η SC The variation with the threshold value T is shown schematically.
Detailed Description
Reference will now be made in detail to embodiments of the present invention, examples of which are illustrated in the accompanying drawings, wherein like reference numerals refer to the same or similar elements or elements having the same or similar function throughout. The embodiments described below with reference to the accompanying drawings are exemplary only for explaining the present invention and are not construed as limiting the present invention.
As used herein, the singular forms "a", "an", "the" and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It will be further understood that the terms "comprises" and/or "comprising," when used in this specification, specify the presence of stated features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof. It will be understood that when an element is referred to as being "connected" or "coupled" to another element, it can be directly connected or coupled to the other element or intervening elements may also be present. Further, "connected" or "coupled" as used herein may include wirelessly connected or coupled. As used herein, the term "and/or" includes any and all combinations of one or more of the associated listed items.
It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. It will be further understood that terms, such as those defined in commonly used dictionaries, should be interpreted as having a meaning that is consistent with their meaning in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless expressly so defined herein.
For the convenience of understanding the embodiments of the present invention, the following description will be further explained by taking several specific embodiments as examples in conjunction with the drawings, and the embodiments are not to be construed as limiting the embodiments of the present invention.
In view of the analysis of dual interface communication and non-linear functions. The invention provides a Diversity Signal Cancellation and Adaptive Threshold Estimation (DSC-ATE) pulse noise suppression method based on double interfaces according to independence and difference of a power line and a wireless channel and consistency of Diversity signals. The method comprises the steps of firstly utilizing a wireless channel to assist a power line channel to extract a noise sample, then estimating an impulse noise threshold value in a self-adaptive mode through the noise sample, and finally separating impulse noise from a received signal.
The invention designs an algorithm based on diversity signal cancellation and adaptive threshold estimation to suppress burst noise in a power line channel in dual-interface communication. First, a wireless channel auxiliary power line channel is used for preprocessing in consideration of the consistency of diversity signals, and noise samples are extracted. In order to accurately estimate the noise threshold, the noise threshold is iteratively estimated using the data samples and an objective function minimization/maximization algorithm, thereby determining the location information of the impulse noise. Simulation results show that the method can effectively improve the suppression capability of the burst noise and the reliability of the communication system. The algorithm can be combined with a sparse theory, compressive sensing and the like for research aiming at other noise models such as narrow-band noise, mixed noise and the like, so that the elimination performance of the algorithm on other noises can be researched.
Fig. 1 shows a model of a dual-interface OFDM (Orthogonal Frequency Division Multiplexing) communication system used in the present invention. The transmitted signals are subjected to OFDM modulation in the power line and the wireless communication module respectively, and are transmitted in a power line and a wireless channel in a centralized mode. The wireless channel is affected by Rayleigh fading and Gaussian white noise, and the power line channel is affected by lognormal fading and impulse noise. The power line signal and the wireless signal at the receiving end simultaneously contain the same transmission information, and only the power line signal has impulse noise with larger amplitude, so that the impulse noise data in the power line signal can be obtained by carrying out cancellation operation on the signal after Discrete Fourier Transform (DFT). That is, impulse noise data can be obtained with the aid of diversity signals transmitted over a wireless channel. And then, the power line and the wireless signal are subjected to combination processing and decoding output. And storing the extracted pulse noise data, namely establishing a data sample space to provide data support for the algorithm or other algorithms.
Radio channel fading coefficient h W Satisfy Rayleigh distribution, then h W The probability density function PDF of (1) is:
Figure BDA0003759573950000111
wherein
Figure BDA0003759573950000112
Is the variance of rayleigh fading.
Power line channel fading coefficient h P Satisfy LogN distribution, then h P Has a probability density function PDF of
Figure BDA0003759573950000113
In the formula of P And
Figure BDA0003759573950000114
are each lnh P Mean and variance of. Normalizing the channel fading energy to obtain
Figure BDA0003759573950000115
Based on the correlation of burst noise in time, the invention establishesA Two-State Markov-Gaussian (TSMG) model. Fig. 2 is a schematic diagram illustrating a markov chain of a TSMG noise model according to an embodiment of the present invention. As shown in FIG. 2, in the TSMG, the statistical behavior of the power line noise is made to be the noise state s k E { B, I }, where B represents that the channel is only interfered by background noise and I represents that the channel is affected by impulse noise. The noise state generation process may use one s 1 ,s 2 ,…,s K A sequence characterization of. For this model, the K +1 th noise state s K+1 Can be represented by a first order Markov process:
Figure BDA0003759573950000116
wherein s is k And s k+1 Respectively representing the noise states, P(s), of the samples at different times K+1 ) Representing the K +1 th sample point, the noise being in state s K+1 The probability value of (2). s 1 Representing the initial noise state s without loss of generality 1 Randomly between a background noise state and an impulse noise state.
Thus, the state change process of the noise can be represented by the state transition probability Ps k s k+1 =P(s k+1 |s k ) And (4) showing. According to the state transition probability, respectively obtaining the probability P that the noise is in the state I and the state B I And P B Is composed of
Figure BDA0003759573950000117
Wherein P is BI Representing the transition probability, P, of state B to state I IB Is the transition probability of state I to state B. s k = I denotes at kth state (sample point), noise state s k An impulse noise state.
Figure BDA0003759573950000121
s k = B indicates the noise state s at the k-th state (sampling point) k Is a background noise state. Available parameters
Figure BDA0003759573950000122
To describe the size of the noisy memories,
Figure BDA0003759573950000123
indicating that the noise is memoryless, an
Figure BDA0003759573950000124
Indicating that the noise has memory, i.e. there is a temporal correlation.
In the following description, the power line channel is subjected to background noise
Figure BDA0003759573950000125
Is still subject to impulse noise
Figure BDA0003759573950000126
Interference of (2) requires a Markov sequence s in the TSMG model 1 ,s 2 ,…,s K Characterization, and background noise
Figure BDA0003759573950000127
And impulse noise
Figure BDA0003759573950000128
Is constrained by the gaussian distribution of the different parameters.
Fig. 3 is a schematic diagram of a frame of a power line and wireless dual-interface communication system for performing signal processing at a receiving end according to an embodiment of the present invention. Firstly, removing a power line and a Cyclic Prefix (CP) of a wireless time domain signal at a receiving end, and performing ideal channel estimation by using a pilot signal; because the two channels transmit the same information, under ideal channel estimation, the data part of the power line signal, which is the same as that of the wireless signal, can be removed through cancellation operation, so that a noise sample in the power line signal is reserved; then extracting the impulse noise in the noise sample through a nonlinear function, and eliminating the impulse noise from the power line signal; and finally, combining, decoding and outputting the diversity signals in the power line and the wireless channel.
The OFDM modulated time domain symbol of the t-th transmission time slot is
Figure BDA0003759573950000129
By a channel fading coefficient of
Figure BDA00037595739500001210
And
Figure BDA00037595739500001211
after the same diversity signals are transmitted by the power line and the wireless parallel channel, the time domain OFDM sampling signals received by the receiving end
Figure BDA00037595739500001212
And
Figure BDA00037595739500001213
can be expressed as:
Figure BDA00037595739500001214
Figure BDA00037595739500001215
wherein the symbols
Figure BDA00037595739500001216
Which represents a convolution operation, is a function of,
Figure BDA00037595739500001217
and
Figure BDA00037595739500001218
respectively representing time domain Gaussian noise vectors on a power line and a wireless channel, wherein elements in the vectors respectively satisfy that the mean value is zero and the variance is
Figure BDA00037595739500001219
And
Figure BDA00037595739500001220
(ii) a gaussian distribution of;
Figure BDA00037595739500001221
representing a time-domain impulse noise vector whose elements satisfy a mean of zero and a variance of
Figure BDA00037595739500001222
Is a Gaussian distribution of
Figure BDA00037595739500001223
Is greater than 1. Removing CP from the OFDM signal, and performing DFT to obtain equivalent frequency domain OFDM signal
Figure BDA00037595739500001224
And
Figure BDA00037595739500001225
respectively expressed as:
Figure BDA0003759573950000131
Figure BDA0003759573950000132
wherein F represents the DFT operator and wherein,
Figure BDA0003759573950000133
denotes an OFDM frequency domain symbol vector transmitted by the transmitting end in the t-th transmission slot, N is the number of subcarriers, N] T Representing a transpose operation.
Figure BDA0003759573950000134
And
Figure BDA0003759573950000135
the power line of the t-th transmission time slot and the wireless equivalent channel coefficient matrix are respectively represented, and because CP is utilized to resist frequency selective fading caused by multipath channels, the invention ignores the influence of intersymbol interference, so that
Figure BDA0003759573950000136
And
Figure BDA0003759573950000137
are all diagonal matrices, and
Figure BDA0003759573950000138
it is reversible. Order to
Figure BDA0003759573950000139
Left-hand multiplication matrix Q can be obtained
Figure BDA00037595739500001310
Wherein
Figure BDA00037595739500001311
Figure BDA00037595739500001312
Represent
Figure BDA00037595739500001313
The inverse matrix of (c). Subtracting the formula (8) and the formula (10) to obtain the same signal part in the power line channel receiving end which is cancelled out of the wireless channel
Figure BDA00037595739500001314
The resulting frequency domain noise estimate vector Ψ t Comprises the following steps:
Figure BDA00037595739500001315
wherein
Figure BDA00037595739500001316
Is a frequency domain vector of the time domain impulse noise after DFT transformation,
Figure BDA00037595739500001317
representing the background noise interference term of the dual interface channel.
To psi t Performing an inverse discrete Fourier transform F * Ψ t Obtaining a noisy time-domain signal vector
Figure BDA00037595739500001318
Comprises the following steps:
Figure BDA00037595739500001319
wherein the content of the first and second substances,
Figure BDA00037595739500001320
and
Figure BDA00037595739500001321
a gaussian distribution with a mean value of 0 is obeyed. Under the condition of known Q, the method can be used for
Figure BDA00037595739500001322
Approximated as a mean of 0 and a variance of
Figure BDA00037595739500001323
A gaussian distribution of (a). The present invention performs statistical analysis on noisy data samples, and fig. 4 is a probability density provided by an embodiment of the present invention
Figure BDA00037595739500001324
Is a statistical histogram of
Figure BDA00037595739500001325
Distribution of (2) and Gaussian distribution
Figure BDA00037595739500001326
The heights are overlapped.
Fig. 5 shows a Noise data extraction effect based on dual-interface diversity Signal cancellation according to an embodiment of the present invention, wherein a Signal-to-Noise Ratio (SNR) S NR =10dB, and the other system parameters are shown in table 1. Fig. 5 (a) and 5 (b) show time domain signals received by a power line and a wireless channel, respectively, and fig. 5 (c) shows a processed signal
Figure BDA0003759573950000141
And comparing the result with the actual impulse noise. As can be seen from fig. 5 (a) and 5 (b), the PAPR exists in the OFDM system, and it is difficult to directly acquire the position information of the impulse noise in the power line channel. After the diversity signal cancellation, as shown in FIG. 5 (c), the diversity signal is cancelled
Figure BDA0003759573950000142
Only background interference is left, the influence of the PAPR is minimal, and therefore, the position information of the impulse noise can be acquired using a nonlinear transformation. It is clear that the noise data in fig. 5 (c) can be saved to provide a big data support for the optimal threshold estimation algorithm.
The signal processing process based on the nonlinear function provided by the embodiment of the invention comprises the following steps: as can be seen from FIG. 5 (c), the amplitude is reduced to be less than the threshold T by using the zero-setting algorithm in the nonlinear function t Is replaced by zero, the reserved amplitude is greater than or equal to T t Is the simplest and most effective. Thus to the extracted noisy time-domain signal vector
Figure BDA0003759573950000143
M-th noise element of (1)
Figure BDA0003759573950000144
The zero-setting transformation is adopted:
Figure BDA0003759573950000145
where M represents the total number of samples in the tth transmission slot, | · | represents an absolute value operation,
Figure BDA0003759573950000146
representing the result of the nonlinear transformation of the impulse noise at the mth sampling point of the tth transmission slot.
Order vector
Figure BDA0003759573950000147
To represent
Figure BDA0003759573950000148
Is a result of a non-linear transformation of
Figure BDA0003759573950000149
For power line signal
Figure BDA00037595739500001410
Pulse noise elimination is carried out, and a power line signal without pulse noise can be obtained
Figure BDA00037595739500001411
Comprises the following steps:
Figure BDA00037595739500001412
like equations (8) and (9), F represents a discrete fourier transform operator.
The power line signal after pulse noise removal is then used
Figure BDA00037595739500001413
And wireless signals
Figure BDA00037595739500001414
Merging to obtain y t Finally, the estimated symbol obtained by the maximum likelihood detection algorithm is adopted for the kth subcarrier of the tth transmission time slot
Figure BDA00037595739500001415
Comprises the following steps:
Figure BDA00037595739500001416
wherein | · | purple 2 Which represents the 2-norm operation,
Figure BDA0003759573950000151
represents an equivalent received signal obtained by combining the power line signal and the wireless signal on the kth subcarrier of the t-th time slot, wherein omega represents a modulation signal constellation point set,
Figure BDA0003759573950000152
a symbol on a constellation diagram is represented,
Figure BDA0003759573950000153
which represents the equivalent channel coefficient at the kth sub-carrier of the t-th transmission slot at the receiving end, N is the number of sub-carriers,
Figure BDA0003759573950000154
meaning that the best x is chosen so that the function f (x) on x is minimal.
The nonlinear function adaptive threshold estimation process provided by the embodiment of the invention comprises the following steps:
adaptive threshold design based on traditional methods
How to accurately estimate the threshold T in equation (13) is crucial to the algorithm performance, wherein the conventional Weighted Combination Criterion (WCC) and Siegert Criterion (SC) balance the detection probability and false alarm probability of impulse noise to obtain the optimal threshold. Extensive research has been carried out due to the small number of parameters required for WCC and SC.
Definition P a For good detection probability of impulse noise, P b Is the false alarm probability, where the composite objective function is η WCC =P a -P b Then the optimum threshold T of the WCC criterion WCC Comprises the following steps:
Figure BDA0003759573950000155
the key step of the WCC criterion is to determine the optimal threshold T by targeting the composite function of the detection probability and the false alarm probability as shown in equation (16) WCC Eta is given by the formula (16) WCC The larger the false alarm probability and the smaller the missed detection probability. At a known power of impulse noise
Figure BDA0003759573950000156
Power of and background interference
Figure BDA0003759573950000157
Under the conditions of (1), P may be given separately a And P b Is then calculated from the target function eta WCC The optimum threshold T can be obtained by calculating the partial derivative of the threshold T WCC Comprises the following steps:
Figure BDA0003759573950000158
similarly, let the objective function η SC =P I P a +P B (1-P b ) The optimal threshold T of the SC criterion can be obtained SC
FIG. 6 shows Bit Error Rate (BER) and η according to the embodiment of the present invention WCC And η SC The variation with the threshold value T is shown schematically. It can be seen that there is an optimum threshold T with the goal of BER minimization * But the optimum threshold value T * And threshold T obtained by two criteria WCC 、T SC There is a certain deviation. This is because both criteria solve the threshold with the detection probability and the false alarm probability as the synthetic target, rather than directly using BER as the optimization target.
Adaptive threshold estimation based on data samples
Conventional threshold estimation algorithms require accurate noise parameters. Due to the non-stationarity of the communication environment, the acquisition of parameters is difficult, and a certain deviation exists between the threshold value acquired by the traditional algorithm and the real optimal solution.
Inspired by data-driven machine learning and transfer learning, obtaining the optimal threshold of a nonlinear function through a noise sample and gradient descent method is a new idea for processing impulse noise. Compared with the traditional algorithm, the optimization algorithm based on the noise sample can continuously adjust parameters according to the actual environment so as to achieve the purposes of directly optimizing the target function and reducing the deviation. Since there is no disclosed impulse noise database, the diversity cancellation algorithm proposed by the present invention can be used to construct a noise sample database and used for optimal threshold estimation. In order to ensure the real-time performance of the algorithm, the threshold value needs to be updated iteratively by using the latest noise sample so as to cope with the change of the environment or the parameter.
The optimal threshold estimation algorithm is closely related to the selection of the objective function, and a specific objective function needs to be established for different communication scenes and performance indexes. The invention mainly provides the following two forms:
1) Minimization of BER: in the OFDM communication system, the performance is most commonly evaluated using the BER, which is one of the indicators in the Quality of Service (QoS). The invention compares the symbol obtained by the decoding judgment of the formula (15) with the symbol sent by the sending end, thereby obtaining the error rate and enabling Pr (y) t (T t ) Means that the threshold value T is used at the T-th time t The obtained error rate of one OFDM frame can be expressed as:
Figure BDA0003759573950000161
wherein
Figure BDA0003759573950000162
Is indicated at a noise threshold T t On the kth subcarrier of the next tth transmission time slot, combining the power line signal and the wireless signal to obtain an equivalent receiving signal, | · | | | sweet wind 2 Representing 2-norm operation, | ·| luminance 0 Representing a 0-norm operation (i.e., counting the number of non-zero elements), omega represents the set of modulation signal constellation points,
Figure BDA0003759573950000171
a symbol on a constellation diagram is represented,
Figure BDA0003759573950000172
indicating the equivalent channel coefficient on the kth sub-carrier of the tth transmission slot at the receiving end,
Figure BDA0003759573950000173
indicating the symbol transmitted by the transmitting end on the kth subcarrier of the tth transmission slot. The same reasoning is that N is the number of subcarriers of one OFDM data frame,
Figure BDA0003759573950000174
meaning that the best x is chosen such that the function value f (x) for x is minimal. The objective function is then:
Figure BDA0003759573950000175
wherein Pr (y) t (T t ) Means that one OFDM frame uses the threshold T at the tth transmission slot t The obtained OFDM code error rate is calculated,
Figure BDA0003759573950000176
then indicates the selected noise threshold
Figure BDA0003759573950000177
So as to be related to T t Error rate Pr (y) t (T t ) ) minimum threshold.
2) Maximizing the transmission rate R t : the transmission rate and the channel quality are directly related. After the channel bandwidth is normalized, the equivalent signal-to-noise ratio can be improved by eliminating the impulse noise in the channel, and the transmission capability of the channel is indirectly improved. The noise threshold can thus be optimized based on short packet theory with the goal of maximizing the transmission rate of the finite length code, which is at the noise threshold T t Lower average transmission rate R t (T t ) Can be expressed as
Figure BDA0003759573950000178
Wherein
Figure BDA0003759573950000179
Indicating that the receiving end is at the noise threshold T t Equivalent signal-to-noise ratio of k sub-carrier of the next t transmission time slot, V represents channel dispersion, L represents coding length, and Q -1 (. Cndot.) is a complementary error function, and ε represents the bit error rate, then the objective function is
Figure BDA00037595739500001710
Wherein
Figure BDA00037595739500001711
Means that the optimum x is selected so that the function value f (x) with respect to x is maximum.
The characteristics of different objective functions are similar, a similar optimization design method can be adopted, and the optimization process with the minimum BER as the objective function is introduced.
In the t-th transmission time slot, the receiving end can only obtain signals
Figure BDA0003759573950000181
And a channel coefficient matrix
Figure BDA0003759573950000182
And
Figure BDA0003759573950000183
therefore, the bit error rate of the t-th transmission time slot cannot be directly obtained, and the threshold value cannot be optimized. Based on the construction/updating of the real-time noise sample, the invention converts the error rate problem of the minimized current time into the error rate problem of the minimized noise sample, thereby obtaining the optimal threshold value of the current time slot
Figure BDA0003759573950000184
The expression is
Figure BDA0003759573950000185
Wherein L is D Representing the magnitude of the number of noise samples. Under random channel conditions, L D The selection of (a) needs to take into account trade-offs in computational accuracy and complexity.
When the threshold value is updated in the t-th transmission time slot, the weight lambda of each time slot in the data sample is calculated according to the discount factor gamma i ,i=1,2,…,L D . And then according to the threshold value T of the previous time slot (i.e. T-1 time slot) t-1 Combining the weights λ i For data sample D EN Calculating the error rate, and then obtaining the threshold value T by a gradient descent method t-1 Corresponding error rate gradient value
Figure BDA0003759573950000186
Last using learning rate l R Updating the threshold value at the current moment
Figure BDA0003759573950000187
1 ). In the appendix are introduced T t The result of the solving process of (2) is shown in equation (23)
Figure BDA0003759573950000188
Wherein
Figure BDA0003759573950000189
Denotes the value of x with respect to f (x) 0 The value of the gradient of (a) is,
Figure BDA00037595739500001810
is indicated at a noise threshold T t-1 And combining the power line signal and the wireless signal on the kth subcarrier of the ith OFDM data frame in the noise sample to obtain an equivalent received signal.
Figure BDA00037595739500001811
RepresentIn the noise samples, the symbol sent by the transmitting end on the k subcarrier of the ith OFDM data frame.
Figure BDA00037595739500001812
Indicating the number of errors in the ith OFDM data frame in the noise sample.
Figure BDA00037595739500001813
Then the weight of the ith OFDM data frame in the noise sample in the data sample is represented, i.e.
Figure BDA00037595739500001814
L D N represents the total number of symbols in the data sample, hence
Figure BDA0003759573950000191
Is the average bit error rate after weighted summation based on data samples.
In addition, the present invention employs a queue storage (first-in first-out criterion) to save/update noise samples, thereby ensuring the validity of noise data. The specific steps of the algorithm are as follows:
Figure BDA0003759573950000192
through the above steps, a more accurate threshold value can be obtained. Reconstruction of impulse noise by equation (13)
Figure BDA0003759573950000193
And subtracting the reconstructed impulse noise from the frequency domain OFDM signal before demodulation by the formula (14), and then carrying out operations such as merging processing, decoding judgment and the like on the received signal from which the impulse noise is removed to obtain a data symbol.
The characteristics of the algorithm of the present invention are analyzed in terms of computational complexity, threshold accuracy, update rate, and the like.
1) And (3) complexity analysis: in the process of taking the minimum BER as a target, the complexity of the algorithm mainly relates to three operations of judgment, addition and subtraction and multiplicationAnd (4) calculating. The judgment and addition-subtraction are linear operations with the highest operation efficiency, mainly include formula (13) nonlinear transformation, (14) impulse noise removal and (18) decoding judgment, and the multiplication mainly includes
Figure BDA0003759573950000201
DFT, etc. Through analysis, in one iteration calculation, the complexity of linear operation quantity and multiplication operation is O (kN).
2) Threshold accuracy: the noise sample is directly from the physical environment where the communication is located, the communication quality at the current moment can be accurately reflected, and a more accurate threshold value is obtained on the premise that a noise model and parameters are not needed. The traditional algorithm needs to obtain statistical data and relevant models of the environment, and the phenomenon that the actual environment is not matched with the models can occur.
3) The updating rate is as follows: when the environment changes, it is required to quickly adapt to the dynamic change of the channel or the environment. By means of the learning rate and the discount factor, the DSC-ATE algorithm can dynamically adjust the threshold value according to environmental changes, and the deviation from the actual threshold value is reduced. The traditional algorithm is constrained by a statistical model, and when the model parameters are not changed, the threshold value cannot be adjusted quickly.
In summary, compared with the existing method, the embodiments of the present invention mainly contribute to the following:
1) Aiming at a dual-interface communication architecture, a power line noise sample extraction method based on diversity signal cancellation is innovatively provided. The noise sample is extracted by utilizing the consistency of the power line and the wireless transmission diversity signal, the limitation that only the power line is used for independently processing the impulse noise in the traditional hybrid communication is broken through, and the method can be applied to the prediction of the optimal threshold value.
2) Based on diversity transmission, a sample space of impulse noise is constructed, and the sample space comprises noise samples with low signal-to-noise ratio (abnormal communication); and by combining a nonlinear function, an optimal threshold estimation algorithm based on a noise sample is provided, so that the bit error rate minimization of a communication system can be realized, and the prior information of impulse noise is not required.
3) For the problem of real-time change of the optimal threshold in a non-stationary environment, the convergence speed of the algorithm is adjusted by introducing parameters such as a learning rate and a discount factor, so that the algorithm can better utilize noise samples, the robustness is improved, and the effective compromise of the algorithm between the robustness and the convergence speed is realized.
Those of ordinary skill in the art will understand that: the figures are schematic representations of one embodiment, and the blocks or processes shown in the figures are not necessarily required to practice the present invention.
From the above description of the embodiments, it is clear to those skilled in the art that the present invention can be implemented by software plus necessary general hardware platform. Based on such understanding, the technical solutions of the present invention may be embodied in the form of a software product, which may be stored in a storage medium, such as ROM/RAM, magnetic disk, optical disk, etc., and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute the method according to the embodiments or some parts of the embodiments.
The embodiments in the present specification are described in a progressive manner, and the same and similar parts among the embodiments are referred to each other, and each embodiment focuses on the differences from the other embodiments. In particular, apparatus or system embodiments, which are substantially similar to method embodiments, are described in relative ease, and reference may be made to some descriptions of method embodiments for related points. The above-described embodiments of the apparatus and system are merely illustrative, and the units described as separate parts may or may not be physically separate, and the parts displayed as units may or may not be physical units, may be located in one place, or may be distributed on a plurality of network units. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of the present embodiment. One of ordinary skill in the art can understand and implement it without inventive effort.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention should be subject to the protection scope of the claims.

Claims (6)

1. A power line noise sample extraction method based on diversity signal cancellation is characterized by comprising the following steps:
diversity signals to be transmitted are diversity-transmitted in a power line and a wireless channel, respectively, the power line being affected by impulse noise;
respectively carrying out discrete Fourier transform on the diversity signals received by the receiving ends of the power line channel and the wireless channel, and then obtaining noise data in the power line diversity signals through cancellation operation;
the impulse noise data is cancelled from the power line diversity signal received at the receiving end by a non-linear function based on an optimal threshold estimate of the noise samples.
2. The method of claim 1, wherein the diversity signals to be transmitted are diversity transmitted in a power line and a wireless channel, respectively, the power line being affected by impulse noise, comprising:
diversity signal to be transmitted
Figure FDA0003759573940000014
Respectively carrying out diversity transmission in a power line and a wireless channel, wherein the wireless channel is influenced by Rayleigh fading and white Gaussian noise, the power line channel is influenced by logarithmically normal fading and pulse noise, and the diversity signals are respectively subjected to Orthogonal Frequency Division Multiplexing (OFDM) modulation in the power line and a wireless communication module;
radio channel fading coefficient h W Satisfy Rayleigh distribution, then h W The probability density function PDF of (1) is:
Figure FDA0003759573940000011
wherein
Figure FDA0003759573940000012
Is the rayleigh fading variance;
power line channel fading coefficient h P Satisfy LogN distribution, then h P The probability density function PDF of (1) is:
Figure FDA0003759573940000013
in the formula of P And
Figure FDA0003759573940000021
are each lnh P Normalizing the channel fading energy to obtain
Figure FDA0003759573940000022
3. The method of claim 1, wherein the obtaining impulse noise data in the power line diversity signal by performing discrete fourier transform on the diversity signals received by the receiving ends of the power line channel and the wireless channel and performing cancellation operation comprises:
the OFDM modulated time domain symbol of the t-th transmission time slot is
Figure FDA0003759573940000023
Diversity signal to be transmitted by a transmitting end
Figure FDA0003759573940000024
Respectively in the channel fading coefficient of
Figure FDA0003759573940000025
The power line channel and the channel fading coefficient are
Figure FDA0003759573940000026
In a wireless channel, a power line channel and a receiving end of the wireless channel receive time domain OFDM sampling signals
Figure FDA0003759573940000027
And
Figure FDA0003759573940000028
expressed as:
Figure FDA0003759573940000029
Figure FDA00037595739400000210
wherein the symbols
Figure FDA00037595739400000211
Represents a convolution operation;
Figure FDA00037595739400000212
and
Figure FDA00037595739400000213
respectively representing time domain Gaussian noise vectors on a power line and a wireless channel, wherein elements in the vectors respectively satisfy the conditions that the mean value is zero and the variance is
Figure FDA00037595739400000214
And
Figure FDA00037595739400000215
(ii) a gaussian distribution of;
Figure FDA00037595739400000216
representing a time-domain impulse noise vector whose elements satisfy a mean of zero and a variance of
Figure FDA00037595739400000217
Is gaussian distribution of (a) and
Figure FDA00037595739400000218
greater than 1;
will be provided with
Figure FDA00037595739400000219
And
Figure FDA00037595739400000220
after the cyclic prefix is removed, discrete Fourier Transform (DFT) is carried out to obtain equivalent frequency domain OFDM signals of the power line channel and the wireless channel
Figure FDA00037595739400000221
And
Figure FDA00037595739400000222
respectively expressed as:
Figure FDA00037595739400000223
Figure FDA00037595739400000224
wherein F represents the DFT operator and wherein,
Figure FDA00037595739400000225
denotes an OFDM frequency domain symbol vector transmitted by the transmitting end in the t-th transmission slot, N is the number of subcarriers, N] T Which represents the operation of transposition of the image,
Figure FDA00037595739400000226
and
Figure FDA00037595739400000227
equivalent channel coefficient matrices respectively representing the power line and the wireless channel of the tth transmission slot,
Figure FDA00037595739400000228
and
Figure FDA00037595739400000229
are all diagonal matrices, and
Figure FDA00037595739400000230
is reversible, make
Figure FDA00037595739400000231
The left-hand matrix Q can be obtained by:
Figure FDA0003759573940000031
wherein
Figure FDA0003759573940000032
Figure FDA0003759573940000033
To represent
Figure FDA0003759573940000034
The inverse matrix of (5) and (7) is subtracted to obtain the same signal part in the power line channel receiving end which is cancelled out of the wireless channel
Figure FDA0003759573940000035
The post frequency domain noise estimate vector Ψ t Comprises the following steps:
Figure FDA0003759573940000036
wherein
Figure FDA0003759573940000037
Is a frequency domain vector of the impulse noise of the time domain after DFT transformation,
Figure FDA0003759573940000038
representing a background noise interference term of the dual interface channel;
let F * Representing the inverse discrete Fourier transform operator, then for Ψ t Performing an inverse discrete Fourier transform F * Ψ t Obtaining a time-domain noise signal vector
Figure FDA00037595739400000319
Comprises the following steps:
Figure FDA0003759573940000039
wherein
Figure FDA00037595739400000310
And
Figure FDA00037595739400000311
subject to a Gaussian distribution with a mean value of 0, given a known Q, will
Figure FDA00037595739400000312
Approximate mean of 0 and variance of
Figure FDA00037595739400000313
A gaussian distribution of (a).
4. The method of claim 3, wherein the removing the impulse noise data from the power line diversity signal received from the receiving end by the nonlinear function based on the optimal threshold estimation of noise samples comprises:
determining threshold T using a data sample based adaptive threshold estimation method t Using a zero-setting algorithm in a non-linear function to reduce the amplitude below a threshold T t Is replaced by zero, the reserved amplitude is greater than or equal to T t Of the noise signal
Figure FDA00037595739400000314
For the extracted noise time domain signal vector
Figure FDA00037595739400000315
M-th noise element of (1)
Figure FDA00037595739400000316
The zero-setting transformation is adopted:
Figure FDA00037595739400000317
where M represents the total number of samples in the tth transmission slot, |, represents the absolute value operation,
Figure FDA00037595739400000318
a nonlinear transformation result representing impulse noise at the mth sampling point of the tth transmission time slot;
order vector
Figure FDA0003759573940000041
Represent
Figure FDA0003759573940000042
Is a result of a non-linear transformation of
Figure FDA0003759573940000043
For power line signal
Figure FDA0003759573940000044
Eliminating impulse noise to obtain power line signal without impulse noise
Figure FDA0003759573940000045
Comprises the following steps:
Figure FDA0003759573940000046
like equations (5) and (6), F represents a discrete fourier transform operator.
5. The method of claim 4, further comprising:
will be provided with
Figure FDA0003759573940000047
With signals received at the wireless port
Figure FDA0003759573940000048
Merging to obtain y t Estimating symbol obtained by maximum likelihood detection algorithm for kth subcarrier of t-th transmission time slot
Figure FDA0003759573940000049
Is composed of
Figure FDA00037595739400000410
Wherein | · | purple 2 Which represents the 2-norm operation,
Figure FDA00037595739400000411
represents an equivalent received signal obtained by combining the power line signal and the wireless signal on the kth subcarrier of the t-th time slot, wherein omega represents a modulation signal constellation point set,
Figure FDA00037595739400000412
a symbol on a constellation diagram is represented,
Figure FDA00037595739400000413
which represents the equivalent channel coefficient at the kth subcarrier of the t transmission slot at the receiving end, N is the number of subcarriers,
Figure FDA00037595739400000414
meaning that the best x is chosen such that the function value f (x) for x is minimal.
6. The method of claim 5 wherein the threshold T is determined using a data sample-based adaptive threshold estimation method t The method comprises the following steps:
minimizing the error rate: let Pr (y) t (T t ) Means that the threshold value T is used in one OFDM data frame in the tth transmission slot t The obtained OFDM bit error rate is expressed as:
Figure FDA00037595739400000415
wherein | · | purple 0 Representing a 0-norm operation (i.e. counting the number of non-zero elements),
Figure FDA00037595739400000416
indicating the symbol sent by the sending end on the kth subcarrier of the tth transmission slot, where N is the number of subcarriers of one OFDM data frame, the objective function is:
Figure FDA00037595739400000417
wherein
Figure FDA0003759573940000051
Indicating a selected noise threshold
Figure FDA0003759573940000052
So as to relate to T t Error rate Pr (y) t (T t ) Minimum threshold, s.tT) t ≧ 0 indicates satisfaction of the "noise threshold T t The selection range of the receiving terminal is greater than or equal to 0 ″, and the receiving terminal obtains the signal at the t-th transmission time slot
Figure FDA0003759573940000053
And a channel coefficient matrix
Figure FDA0003759573940000054
And
Figure FDA0003759573940000055
based on the construction/updating of the real-time noise sample, the error rate problem at the current moment is converted into the average error rate problem of the minimized noise sample, and the optimal threshold value at the current moment is obtained
Figure FDA0003759573940000056
The expression is as follows:
Figure FDA0003759573940000057
wherein L is D Representing the magnitude of the number of noise samples, pr (y) i (T t ) Means that the ith OFDM data frame in the noise sample is related to the noise threshold T t S.t is an abbreviation of subject to, indicating that the constraint condition is satisfied;
when the threshold value is updated in the t-th transmission time slot, the weight lambda of each time slot in the data sample is calculated according to the discount factor gamma i ,i=1,2,…,L D According to the threshold value T at the T-1 th moment t-1 Combining the weights λ i For data sample D EN Calculating the error rate, and obtaining the threshold value T by a gradient descent method t-1 Corresponding gradient value of error rate
Figure FDA00037595739400000516
Using learning rate
Figure FDA00037595739400000517
Updating the threshold value at the current moment
Figure FDA0003759573940000058
Figure FDA0003759573940000059
Is shown in equation (16):
Figure FDA00037595739400000510
wherein
Figure FDA00037595739400000511
Denotes the value of x with respect to f (x) 0 The value of the gradient of (a) is,
Figure FDA00037595739400000512
is indicated at a noise threshold T t-1 Then, on the kth subcarrier of the ith OFDM data frame in the noise sample, combining the power line signal and the wireless signal to obtain an equivalent received signal,
Figure FDA00037595739400000513
indicating the symbol transmitted by the transmitting end on the k subcarrier of the ith OFDM data frame among the noise samples,
Figure FDA00037595739400000514
indicating the number of errors in the ith OFDM data frame in the noise sample,
Figure FDA00037595739400000515
then the weight of the ith OFDM data frame in the noise sample in the data sample is represented, i.e.
Figure FDA0003759573940000061
L D N represents the total number of symbols in the data sample, hence
Figure FDA0003759573940000062
Is the average bit error rate after weighted summation based on data samples.
CN202210866696.0A 2022-07-22 2022-07-22 Power line noise sample extraction method based on diversity signal cancellation Active CN115378467B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210866696.0A CN115378467B (en) 2022-07-22 2022-07-22 Power line noise sample extraction method based on diversity signal cancellation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210866696.0A CN115378467B (en) 2022-07-22 2022-07-22 Power line noise sample extraction method based on diversity signal cancellation

Publications (2)

Publication Number Publication Date
CN115378467A true CN115378467A (en) 2022-11-22
CN115378467B CN115378467B (en) 2024-04-16

Family

ID=84061000

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210866696.0A Active CN115378467B (en) 2022-07-22 2022-07-22 Power line noise sample extraction method based on diversity signal cancellation

Country Status (1)

Country Link
CN (1) CN115378467B (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2519011C1 (en) * 2013-02-20 2014-06-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Омский государственный технический университет" Method of transmitting information over short-wave link using frequency-shift keyed signals
CN104852879A (en) * 2015-04-07 2015-08-19 湖南大学 Method and system for dynamically suppressing power line impulse noise
CN106877907A (en) * 2017-01-09 2017-06-20 北京邮电大学 A kind of anti-noise channel estimation methods in power line
CN106992800A (en) * 2017-03-16 2017-07-28 宁波大学 Electric line communication system impulse noise suppression method based on iteration self-adapting algorithm
CN110493151A (en) * 2019-07-04 2019-11-22 宁波大学 The impulse noise suppression method of sparse iteration covariance estimation in power line

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2519011C1 (en) * 2013-02-20 2014-06-10 Федеральное государственное бюджетное образовательное учреждение высшего профессионального образования "Омский государственный технический университет" Method of transmitting information over short-wave link using frequency-shift keyed signals
CN104852879A (en) * 2015-04-07 2015-08-19 湖南大学 Method and system for dynamically suppressing power line impulse noise
CN106877907A (en) * 2017-01-09 2017-06-20 北京邮电大学 A kind of anti-noise channel estimation methods in power line
CN106992800A (en) * 2017-03-16 2017-07-28 宁波大学 Electric line communication system impulse noise suppression method based on iteration self-adapting algorithm
CN110493151A (en) * 2019-07-04 2019-11-22 宁波大学 The impulse noise suppression method of sparse iteration covariance estimation in power line

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
陈智雄;韩东升;邱丽君;: "室内无线和电力线双媒质协作通信***性能研究", 中国电机工程学报, no. 09, 5 May 2017 (2017-05-05), pages 2589 - 2598 *

Also Published As

Publication number Publication date
CN115378467B (en) 2024-04-16

Similar Documents

Publication Publication Date Title
CN112637093B (en) Signal detection method based on model-driven deep learning
CN108900443B (en) A kind of underwater acoustic channel interference elimination method in underwater sound communication
WO2009099013A1 (en) Mobile communication system, reception device, and method
CN107332797B (en) Channel estimation method in power line OFDM communication system
CN112347871B (en) Interference signal modulation identification method for communication carrier monitoring system
WO2007137484A1 (en) A channel estimation method and the device thereof
CN113472706A (en) MIMO-OFDM system channel estimation method based on deep neural network
CN100571098C (en) The maximum likelihood detecting method of low complex degree and device in the communication system
CN110311876A (en) The implementation method of underwater sound OFDM receiver based on deep neural network
CN101018219A (en) Space frequency signal processing method
CN115250216A (en) Underwater sound OFDM combined channel estimation and signal detection method based on deep learning
Ru et al. Model-driven channel estimation for OFDM systems based on image super-resolution network
Arunkumar et al. Iterative sparse channel estimation and data detection for underwater acoustic communications using partial interval demodulation
Zhang et al. Efficient residual shrinkage CNN denoiser design for intelligent signal processing: Modulation recognition, detection, and decoding
Bai et al. Impulsive noise mitigation with interleaving based on MUSIC in power line communication
Gizzini et al. Adaptive channel estimation based on deep learning
CN111291511B (en) Soft Kalman filtering iteration time-varying channel estimation method based on historical information
CN112636855A (en) OFDM signal detection method
CN102025459A (en) Nonparametric estimation ICA-based MIMO-OFDM system blind deconvolution method
CN115378467B (en) Power line noise sample extraction method based on diversity signal cancellation
CN113709075B (en) Method for realizing underwater acoustic communication receiver by using underwater acoustic channel multipath effect
CN116094885A (en) Signal modulation mode identification method based on One2ThreeNet
CN113726710B (en) OFDM signal transmitting method, OFDM signal receiving method, OFDM signal transmitting device, OFDM signal receiving device, storage medium and electronic device
CN113067652B (en) Channel prediction method by using space characteristics of wide-band large-scale antenna array
Mahmood et al. A design of extreme learning machine based receiver for 2× 2 MIMO-OFDM system

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