CN107659527A - A kind of while co-channel full duplex communication for coordination phase noise reduction system and method - Google Patents

A kind of while co-channel full duplex communication for coordination phase noise reduction system and method Download PDF

Info

Publication number
CN107659527A
CN107659527A CN201710843738.8A CN201710843738A CN107659527A CN 107659527 A CN107659527 A CN 107659527A CN 201710843738 A CN201710843738 A CN 201710843738A CN 107659527 A CN107659527 A CN 107659527A
Authority
CN
China
Prior art keywords
mrow
msub
signal
mfrac
msup
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
CN201710843738.8A
Other languages
Chinese (zh)
Other versions
CN107659527B (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201710843738.8A priority Critical patent/CN107659527B/en
Publication of CN107659527A publication Critical patent/CN107659527A/en
Application granted granted Critical
Publication of CN107659527B publication Critical patent/CN107659527B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2691—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation involving interference determination or cancellation
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/022—Channel estimation of frequency response
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03006—Arrangements for removing intersymbol interference
    • H04L25/03178—Arrangements involving sequence estimation techniques
    • H04L25/03248—Arrangements for operating in conjunction with other apparatus
    • H04L25/0328—Arrangements for operating in conjunction with other apparatus with interference cancellation circuitry
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2692—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with preamble design, i.e. with negotiation of the synchronisation sequence with transmitter or sequence linked to the algorithm used at the receiver
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/26—Systems using multi-frequency codes
    • H04L27/2601—Multicarrier modulation systems
    • H04L27/2647—Arrangements specific to the receiver only
    • H04L27/2655—Synchronisation arrangements
    • H04L27/2689—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation
    • H04L27/2695—Link with other circuits, i.e. special connections between synchronisation arrangements and other circuits for achieving synchronisation with channel estimation, e.g. determination of delay spread, derivative or peak tracking
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L27/00—Modulated-carrier systems
    • H04L27/32—Carrier systems characterised by combinations of two or more of the types covered by groups H04L27/02, H04L27/10, H04L27/18 or H04L27/26
    • H04L27/34—Amplitude- and phase-modulated carrier systems, e.g. quadrature-amplitude modulated carrier systems
    • H04L27/38—Demodulator circuits; Receiver circuits
    • H04L27/3845—Demodulator circuits; Receiver circuits using non - coherent demodulation, i.e. not using a phase synchronous carrier
    • H04L27/3854—Demodulator circuits; Receiver circuits using non - coherent demodulation, i.e. not using a phase synchronous carrier using a non - coherent carrier, including systems with baseband correction for phase or frequency offset
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00—Arrangements affording multiple use of the transmission path
    • H04L5/14—Two-way operation using the same type of signal, i.e. duplex
    • H04L5/1461—Suppression of signals in the return path, i.e. bidirectional control circuits

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Power Engineering (AREA)
  • Radio Relay Systems (AREA)
  • Other Investigation Or Analysis Of Materials By Electrical Means (AREA)

Abstract

The invention discloses a kind of while co-channel full duplex communication for coordination phase noise reduction system and method, including source sub-system, trunk subsystem and terminal receiving subsystem；The output end of the source sub-system is connected by trunk subsystem with terminal receiving subsystem；The source sub-system, for producing echo signal and caused signal being handled, the signal that processing obtains is sent to trunk subsystem；The trunk subsystem, for receiving the signal from source sub-system, it is transformed into numeric field and carries out reception processing, digital suppression is carried out to the signal of reception processing in numeric field, to numeral suppress after signal be amplified and pilot tone insertion after, transmitting processing is carried out, and is transformed into analog domain, transmits a signal to receiving subsystem；Terminal receiving subsystem, for handling the signal from trunk subsystem, complete the reception of echo signal.Present invention reduces the bit error rate of signal transmission, it is suppressed that the influence for interference of mutually making an uproar.

Description

A kind of while co-channel full duplex communication for coordination phase noise reduction system and method
Technical field
The present invention relates to the cooperative transmission of trunk subsystem and phase noise reduction, more particularly to a kind of while complete with frequency The phase noise reduction system and method for duplexing communication for coordination.
Background technology
It is harmonic wave and intermodulation product caused by nonlinear interaction in non-linear equipment mutually to make an uproar, by oscillator it is non-ideal because The noise that element introduces forms random phase modulation to outgoing carrier, and in general, free oscillation device, which is mutually made an uproar, to be modeled as tieing up Receive process.Mutually make an uproar in multi-carrier transmission, such as using OFDM (orthogonal frequency division Multiplexing, OFDM) modulation system transmission system in, the inter-carrier interference (inter-carrier of signal can be caused Interference, ICI) and common phase error (common phase error, CPE), wherein, common phase error part Energy be much larger than inter-carrier interference part energy.Therefore, the influence mutually made an uproar can greatly be weakened by suppressing CPE；It is existing Have in mode, can estimate CPE and self-interference channel simultaneously, and remove in self-interference signal part CPE so as to Suppress the influence of the interference of mutually making an uproar of self-interference signal in relaying；Meanwhile power distribution rational on each subcarrier at relaying By the effective signal to noise ratio for improving receiving terminal, though prior art considers the power distribution under single carrier case, do not provide Optimal power distribution computational methods, it is unfavorable for reducing the bit error rate of signal transmission, is also unfavorable for suppressing the influence of interference of mutually making an uproar.
The content of the invention
It is an object of the invention to overcome the deficiencies of the prior art and provide a kind of while co-channel full duplex communication for coordination phase Position noise suppressing system and method.
The purpose of the present invention is achieved through the following technical solutions：A kind of while co-channel full duplex communication for coordination phase Position noise suppressing system, including source sub-system, trunk subsystem and terminal receiving subsystem；The output of the source sub-system End is connected by trunk subsystem with terminal receiving subsystem；
The source sub-system, for producing echo signal and caused signal being handled, obtained letter will be handled Number it is sent to trunk subsystem；Specifically, the source sub-system includes information source, and the output end of the information source passes sequentially through first Pilot tone insertion module, the first inverse Fourier transform module, first circulation prefix insertion module, the first D/A converter module and information source Emitter is connected, and the output signal of analog-to-digital conversion module is modulated to frequency band by information source emitter in oscillating fashion, and passes through letter Source antenna is sent to trunk subsystem.
The trunk subsystem, for receiving the signal from source sub-system, it is transformed into numeric field and carries out receiving area Reason, digital suppression is carried out to the signal of reception processing in numeric field, the signal after suppressing to numeral is amplified and pilot tone is inserted Afterwards, transmitting processing is carried out in next cycle, and is transformed into analog domain, transmit a signal to receiving subsystem；Specifically, it is described Trunk subsystem includes relay reception antenna, and the relay reception antenna, for receiving the signal from source sub-system, its is defeated Go out end to be connected with link receiver, reception signal is modulated to base band by link receiver in oscillating fashion, and link receiver is defeated The signal gone out presses down through the first analog-to-digital conversion module, first circulation prefix removal module, the first fourier transformation module, numeral successively Molding block, amplification module, the second pilot tone insertion module, the second inverse Fourier transform module, second circulation prefix insertion module, the After two D/A converter modules are handled, obtained signal is transferred to retransmitter, retransmitter turns the second modulus The output signal of mold changing block is modulated to frequency band in oscillating fashion, and sends a signal to reception subsystem by repeat transmitted antenna System；The output end of the second pilot tone insertion module is also connected with digital suppression module.
Terminal receiving subsystem, for handling the signal from trunk subsystem, complete the reception of echo signal. Specifically, described terminal receiving subsystem includes terminal antenna, and the terminal antenna receives the signal from retransmitter, The output end of terminal antenna is connected with terminal receiver, interrupts receiver with mode of oscillation by the signal modulation received to base band Afterwards, the second analog-to-digital conversion module, second circulation prefix removal module and the second fourier transformation module are transmitted a signal to successively Handled, the signal finally given by the output of the second fourier transformation module.
A kind of while co-channel full duplex communication for coordination phase noise inhibition method, comprises the following steps：
S1. source sub-system produces echo signal and caused signal is handled, and the signal that processing is obtained is sent To trunk subsystem；Specifically, the step S1 includes following sub-step：S101. hair is produced in m-th of communication cycle, information source It is P to penetrate powersSignalS102. N is uniformly inserted in signal caused by information sourcepIndividual pilot tone, the position collection of pilot tone It is combined into：
S103. the rear signal for inserting pilot tone is subjected to inverse Fourier transform, insertion cyclic prefix and digital-to-analogue conversion successively, obtained To analog signal to be launched and frequency domain is modulated to mode of oscillation, is sent to trunk subsystem.
S2. trunk subsystem receives the signal from source sub-system, is transformed into numeric field and carries out reception processing, in number Word domain carries out digital suppression to the signal of reception processing, the signal after suppressing to numeral be amplified and pilot tone insertion after, under A cycle carries out transmitting processing, and is transformed into analog domain, transmits a signal to receiving subsystem；Specifically, the step S2 Including following sub-step：S201. trunk subsystem is in oscillating fashion by the signal modulation received to base band, and base band is believed Number carry out analog-to-digital conversion, go after cyclic prefix and Fourier transformation；S202. the signal obtained to Fourier transformation carries out digital suppression System, the signal after numeral is suppressed is amplified to be inserted with pilot tone, and using the pilot signal of insertion as next cycle numeral The basis of suppression；S203. by after the signal delay a cycle obtained in step S202, through inverse Fourier transform to time domain, then To signal insertion cyclic prefix, digital-to-analogue conversion is carried out, frequency band is modulated to mode of oscillation, is sent to terminal receiving subsystem.
S3. terminal receiving subsystem is handled the signal from trunk subsystem, completes the reception of signal.Specifically Ground, the step S3 include following sub-step：Receiving terminal system with mode of oscillation by the signal modulation received to base band, and Successively through analog-to-digital conversion, go cyclic prefix and Fourier transformation after, the m+1 cycle has obtained the signal of final reception
Wherein, the step S202 includes following sub-step：
The signal received using m-th of cycle of trunk subsystem in pilot frequency locations, except above a cycle relays subsystem Unite in the signal of pilot frequency locations insertion, obtain the sub-carrier channels estimate of pilot frequency locations, using these estimates linearly to insert The mode of value obtains the estimation of other subcarrier upper signal channels, from trunk subsystem in the signal that m-th of cycle receives, subtracts The transmission signal and the product of estimation channel of upper a cycle delay, are suppressed with the self-interference numeral of real current period；
Calculate the optimal value of the sub-carrier power of transmission signal in trunk subsystem
By the power amplification of the signal of k-th of subcarrier after self-interference numeral suppression to Pk,Letter after being amplified Number；
I-th of position insertion power of signal after amplification is PtPilot tone, i ∈ D, using obtained signal as next The basis that cycle numeral suppresses.
Wherein, in trunk subsystem the sub-carrier power of transmission signal optimal valueCalculation procedure is as follows：
The first step, the value of parameter in setting system, including the total power constraint P at information sourceS, the total power constraint at relaying PR；Initialization tag value flag=0 and point=0 are set；
Second step, judge whether to meet condition：PS≤Np·PsOr PR≤Np·Pt, if so, terminate and reset be System value of consult volume；If it is not, then enter the 3rd step；
3rd step, calculate the signal power on information source useful signal subcarrier
4th step, calculateZero point PR0If PR0＞ PR-NpPt, by PR0Again it is entered as PR-NpPt, wherein：
εkIt is the mean-square value of the evaluated error of k-th of self-interference channel：
Wherein, A (i) represents sequenceNcI-th of value of point DFT transform, Δ f is three dB bandwidth of mutually making an uproar With the ratio of carrier wave frequency domain interval, ErrAnd EsrRepresent the energy of channel impulse response, Hsr[k] and Hrd[k] represents from information source respectively K-th of the value changed to relaying and from the DFT for the channel impulse response for being relayed to destination node；
5th step, outage threshold ε=10 are set-7, arrange parameter a1=0, b1=PR0, lef=a1+0.382·(b1-a1), Rig=a1+0.618·(b1-a1)；
6th step, judge whether to meet flag=0, if so, arrange parameter Pr=lef；If it is not, arrange parameter Pr=rig；
7th step, is proceeded as follows：
A, the power distribution P on each useful signal subcarrier at relaying is initialized0, P0It is a 1 × NcVector：
InitializationIteration stopping condition η is seto=10-2, ηi=10-5, t=4096, increment times mu=10 with And total constraint number m=2 (Nc-Np)；
B, set withFor the function of variableFor：
Wherein,
Eu[k] represents λ2Ps|Hrd[k]|2|Hsr[k]|2, FkRepresent the inter-carrier interference energy on k-th of subcarrier at relaying Amount, its estimate are
C, judge whether to meetIf it is not, into step F, if so, carrying out following steps：
(1) ifAccording to γFD[k] is on P0,kDerivativeCalculate functionOn P0,kLead Number, P0,kRepresent P0K-th of element value：
Then, according to γFD[k] is on P0,kSecond dervativeCalculate functionOn P0,kSecond order lead Number：
Median w is calculated again：
P is calculated according to median w0 The updated value increment of individual element；
As k ∈ D, by Δ PkIt is entered as 0.To P0It is updated, the P after renewal0Equal to the P before renewal0Plus 0.1 times Δ P, wherein, Δ P represents Δ PkThe vector of composition；
According to w pairs of medianIt is updated, after renewalBe：
(2) if after judging renewalWhether meetIf so, enter operation D, if it is not, return to step (1)；
D, t values are updated, after renewal t be equal to renewal before t be multiplied by mu；
E, judge whether the t after renewal meetsIf so, enter operation F, if it is not, returning to operation C；
F, median mw [k] is calculated：
According to median mw [k], bit error rate r is calculated：
Subsequently determine whether flag is 0, if so, by P0It is assigned to Pf and r is assigned to Bf, otherwise, then by P0It is assigned to Pg And r is assigned to Bg；
8th step, judge whether point is equal to 0, if so, point and flag are entered as into 1 and return the 6th again Step；If otherwise enter the 9th step；
9th step, judge whether to meet | Bf-Bg | ＜ ε, if so, returning to Pf as the optimal of relaying sub-carrier power distribution It is worth the vector of composition, calculating terminates；If it is not, into the tenth step；
Tenth step, compares Bf and Bg：
If Bf ＞ Bg, make a1=lef, lef=rig, rig=a1+0.618·(b1-a1), flag=1, return to the 6th step；
If Bf≤Bg, make b1=rig, rig=lef, lef=a1+0.382·(b1-a1), flag=0, return to the 6th step.
The beneficial effects of the invention are as follows：The present invention is in trunk subsystem, by the common phase error for suppressing self-interference Part, it is suppressed that the influence mutually made an uproar, and in ensuing amplification process, by calculating, the rational work(set on subcarrier Rate further suppress the influence mutually made an uproar to improve the signal to noise ratio of receiving terminal, so reducing the bit error rate of signal transmission, reduce Mutually make an uproar influence of the interference to full duplex communication for coordination.
Brief description of the drawings
Fig. 1 is the system principle diagram of the present invention；
Fig. 2 is flow chart of the method for the present invention；
Fig. 3 is that the CPE of self-interference signal suppresses the influence schematic diagram to error rate of system；
Fig. 4 is influence schematic diagram of the power distribution mode to error rate of system.
Embodiment
Technical scheme is described in further detail below in conjunction with the accompanying drawings, but protection scope of the present invention is not limited to It is as described below.
As shown in figure 1, a kind of while co-channel full duplex communication for coordination phase noise reduction system, including information source subsystem System, trunk subsystem and terminal receiving subsystem；The output end of the source sub-system is received by trunk subsystem and terminal Subsystem connects；
The source sub-system, for producing echo signal and caused signal being handled, obtained letter will be handled Number it is sent to trunk subsystem；Specifically, the source sub-system includes information source, and the output end of the information source passes sequentially through first Pilot tone insertion module, the first inverse Fourier transform module, first circulation prefix insertion module, the first D/A converter module and information source Emitter is connected, and the output signal of analog-to-digital conversion module is modulated to frequency band by information source emitter in oscillating fashion, and passes through letter Source antenna is sent to trunk subsystem.
The trunk subsystem, for receiving the signal from source sub-system, it is transformed into numeric field and carries out receiving area Reason, digital suppression is carried out to the signal of reception processing in numeric field, the signal after suppressing to numeral is amplified and pilot tone is inserted Afterwards, transmitting processing is carried out in next cycle, and is transformed into analog domain, transmit a signal to receiving subsystem；Specifically, it is described Trunk subsystem includes relay reception antenna, and the relay reception antenna, for receiving the signal from source sub-system, its is defeated Go out end to be connected with link receiver, reception signal is modulated to base band by link receiver in oscillating fashion, and link receiver is defeated The signal gone out presses down through the first analog-to-digital conversion module, first circulation prefix removal module, the first fourier transformation module, numeral successively Molding block, amplification module, the second pilot tone insertion module, the second inverse Fourier transform module, second circulation prefix insertion module, the After two D/A converter modules are handled, obtained signal is transferred to retransmitter, retransmitter turns the second modulus The output signal of mold changing block is modulated to frequency band in oscillating fashion, and sends a signal to reception subsystem by repeat transmitted antenna System；The output end of the second pilot tone insertion module is also connected with digital suppression module.
Terminal receiving subsystem, for handling the signal from trunk subsystem, complete the reception of echo signal. Specifically, described terminal receiving subsystem includes terminal antenna, and the terminal antenna receives the signal from retransmitter, The output end of terminal antenna is connected with terminal receiver, interrupts receiver with mode of oscillation by the signal modulation received to base band Afterwards, the second analog-to-digital conversion module, second circulation prefix removal module and the second fourier transformation module are transmitted a signal to successively Handled, the signal finally given by the output of the second fourier transformation module.
As shown in Fig. 2 a kind of while co-channel full duplex communication for coordination phase noise inhibition method, comprises the following steps：
S1. source sub-system produces echo signal and caused signal is handled, and the signal that processing is obtained is sent To trunk subsystem；Specifically, the step S1 includes following sub-step：S101. hair is produced in m-th of communication cycle, information source It is P to penetrate powersSignalS102. N is uniformly inserted in signal caused by information sourcepIndividual pilot tone, the position collection of pilot tone It is combined into：
S103. the rear signal for inserting pilot tone is subjected to inverse Fourier transform, insertion cyclic prefix and digital-to-analogue conversion successively, obtained To analog signal to be launched and frequency domain is modulated to mode of oscillation, is sent to trunk subsystem.
S2. trunk subsystem receives the signal from source sub-system, is transformed into numeric field and carries out reception processing, in number Word domain carries out digital suppression to the signal of reception processing, the signal after suppressing to numeral be amplified and pilot tone insertion after, under A cycle carries out transmitting processing, and is transformed into analog domain, transmits a signal to receiving subsystem；Specifically, the step S2 Including following sub-step：S201. trunk subsystem is in oscillating fashion by the signal modulation received to base band, and base band is believed Number carry out analog-to-digital conversion, go after cyclic prefix and Fourier transformation；S202. the signal obtained to Fourier transformation carries out digital suppression System, the signal after numeral is suppressed is amplified to be inserted with pilot tone, and using the pilot signal of insertion as next cycle numeral The basis of suppression；S203. by after the signal delay a cycle obtained in step S202, through inverse Fourier transform to time domain, then To signal insertion cyclic prefix, digital-to-analogue conversion is carried out, frequency band is modulated to mode of oscillation, is sent to terminal receiving subsystem.
S3. terminal receiving subsystem is handled the signal from trunk subsystem, completes the reception of signal.Specifically Ground, the step S3 include following sub-step：Receiving terminal system with mode of oscillation by the signal modulation received to base band, and Successively through analog-to-digital conversion, go cyclic prefix and Fourier transformation after, the m+1 cycle has obtained the signal of final reception
Wherein, the step S202 includes following sub-step：
The signal received using m-th of cycle of trunk subsystem in pilot frequency locations, except above a cycle relays subsystem Unite in the signal of pilot frequency locations insertion, obtain the sub-carrier channels estimate of pilot frequency locations, using these estimates linearly to insert The mode of value obtains the estimation of other subcarrier upper signal channels, from trunk subsystem in the signal that m-th of cycle receives, subtracts The transmission signal and the product of estimation channel of upper a cycle delay, are suppressed with the self-interference numeral of real current period；
Calculate the optimal value of the sub-carrier power of transmission signal in trunk subsystem
By the power amplification of the signal of k-th of subcarrier after self-interference numeral suppression to Pk,Letter after being amplified Number；
I-th of position insertion power of signal after amplification is PtPilot tone, i ∈ D, using obtained signal as next The basis that cycle numeral suppresses.
Wherein, in trunk subsystem the sub-carrier power of transmission signal optimal valueCalculation procedure is as follows：
The first step, the value of parameter in setting system, including the total power constraint P at information sourceS, the total power constraint at relaying PR；Initialization tag value flag=0 and point=0 are set；
Second step, judge whether to meet condition：PS≤Np·PsOr PR≤Np·Pt, if so, terminate and reset be System value of consult volume；If it is not, then enter the 3rd step；
3rd step, calculate the signal power on information source useful signal subcarrier
4th step, calculateZero point PR0If PR0＞ PR-NpPt, by PR0Again it is entered as PR-NpPt, wherein：
εkIt is the mean-square value of the evaluated error of k-th of self-interference channel：
Wherein, A (i) represents sequenceNcI-th of value of point DFT transform, Δ f is three dB bandwidth of mutually making an uproar With the ratio of carrier wave frequency domain interval, ErrAnd EsrRepresent the energy of channel impulse response, Hsr[k] and Hrd[k] represents from information source respectively K-th of the value changed to relaying and from the DFT for the channel impulse response for being relayed to destination node；
5th step, outage threshold ε=10 are set-7, arrange parameter a1=0, b1=PR0, lef=a1+0.382·(b1-a1), Rig=a1+0.618·(b1-a1)；
6th step, judge whether to meet flag=0, if so, arrange parameter Pr=lef；If it is not, arrange parameter Pr=rig；
7th step, is proceeded as follows：
A, the power distribution P on each useful signal subcarrier at relaying is initialized0, P0It is a 1 × NcVector：
InitializationIteration stopping condition η is seto=10-2, ηi=10-5, t=4096, increment times mu=10 with And total constraint number m=2 (Nc-Np)；
B, set withFor the function of variableFor：
Wherein,
Eu[k] represents λ2Ps|Hrd[k]|2|Hsr[k]|2, FkRepresent the inter-carrier interference energy on k-th of subcarrier at relaying Amount, its estimate are
C, judge whether to meetIf it is not, into step F, if so, carrying out following steps：
(1) ifAccording to γFD[k] is on P0,kDerivativeCalculate functionOn P0,kLead Number, P0,kRepresent P0K-th of element value：
Then, according to γFD[k] is on P0,kSecond dervativeCalculate functionOn P0,kSecond order lead Number：
Median w is calculated again：
P is calculated according to median w0 The updated value increment of individual element；
As k ∈ D, by Δ PkIt is entered as 0.To P0It is updated, the P after renewal0Equal to the P before renewal0Plus 0.1 times Δ P, wherein, Δ P represents Δ PkThe vector of composition；
According to w pairs of medianIt is updated, after renewalBe：
(2) if after judging renewalWhether meetIf so, enter operation D, if it is not, return to step (1)；
D, t values are updated, after renewal t be equal to renewal before t be multiplied by mu；
E, judge whether the t after renewal meetsIf so, enter operation F, if it is not, returning to operation C；
F, median mw [k] is calculated：
According to median mw [k], bit error rate r is calculated：
Subsequently determine whether flag is 0, if so, by P0It is assigned to Pf and r is assigned to Bf, otherwise, then by P0It is assigned to Pg And r is assigned to Bg；
8th step, judge whether point is equal to 0, if so, point and flag are entered as into 1 and return the 6th again Step；If otherwise enter the 9th step；
9th step, judge whether to meet | Bf-Bg | ＜ ε, if so, returning to Pf as the optimal of relaying sub-carrier power distribution It is worth the vector of composition, calculating terminates；If it is not, into the tenth step；
Tenth step, compares Bf and Bg：
If Bf ＞ Bg, make a1=lef, lef=rig, rig=a1+0.618·(b1-a1), flag=1, return to the 6th step；
If Bf≤Bg, make b1=rig, rig=lef, lef=a1+0.382·(b1-a1), flag=0, return to the 6th step.
By suppressing the CPE parts of self-interference in the present invention, it is suppressed that the influence mutually made an uproar, and in ensuing amplification process In, by calculating, the rational power set on subcarrierTo improve the signal to noise ratio of receiving terminal, further suppress The influence mutually made an uproar.
In embodiments herein, effect of the bit error rate with AF panel of mutually making an uproar is reduced for checking, carries out emulation experiment： Total carrier number is arranged to Nc=1057, sampling interval TsFor 3.3 × 10-8S, number of pilots NpFor 33, between inter-carrier frequency Every fcarrFor 15kHz, pilot power is set to 10dB at information source, and the pilot power at relaying is 15dB, therefrom the secondary signal sent It is 0 by tap number, Isosorbide-5-Nitrae and power delay profile (PDP) are 0dB, -5dB, -15dB self-interference channel hrr, information source is arrived The channel of relaying and it is set as from the channel parameter for being relayed to destination node identical, i.e. tap number 0,20,45 and corresponding Power delay profile is 0dB, -9dB, -20dB.Part of mutually making an uproar is modeled as Wiener-Hopf equation, and the three dB bandwidth mutually made an uproar at oscillator f3dB=80Hz.
As shown in figure 3, in the emulation experiment, the CPE of self-interference signal suppresses the influence schematic diagram to error rate of system, Self-interference signal CPE suppression is compared in figure and without error rate of system in the case of self-interference signal CPE two kinds of suppression Change.It can be seen that it will significantly reduce the BER of system by suppressing the CPE of self-interference signal at relaying.In information source work( When rate is 18dB, using the method for the suppression CPE in this patent, the bit error rate can be reduced an order of magnitude.
As shown in figure 4, influence schematic diagram of the power distribution mode to error rate of system, even power distribution is compared in figure With the influence of optimal power allocation both modes to error rate of system, even power allocative decision is that each height carries at relaying The optimal value that power on ripple is got when consistent.It can see from figure, the optimal power allocation method in this patent The bit error rate can be reduced to original half when source signal power is 25dB, system is improved compared to even power method Signal to noise ratio, it is suppressed that the influence for interference of mutually making an uproar.

Claims (10)

1. a kind of while co-channel full duplex communication for coordination phase noise reduction system, it is characterised in that：Including source sub-system, Trunk subsystem and terminal receiving subsystem；The output end of the source sub-system receives subsystem by trunk subsystem and terminal System connection；
The source sub-system, for producing echo signal and caused signal being handled, the signal that processing is obtained is sent out Give trunk subsystem；
The trunk subsystem, for receiving the signal from source sub-system, it is transformed into numeric field and carries out reception processing, Numeric field carries out digital suppression to the signal of reception processing, the signal after suppressing to numeral be amplified and pilot tone insertion after, Next cycle carries out transmitting processing, and is transformed into analog domain, transmits a signal to receiving subsystem；
Terminal receiving subsystem, for handling the signal from trunk subsystem, complete the reception of echo signal.
2. according to claim 1 a kind of while co-channel full duplex communication for coordination phase noise reduction system, its feature It is：The source sub-system includes information source, and the output end of the information source is passed sequentially through in the first pilot tone insertion module, first Fu Leaf inverse transform module, first circulation prefix insertion module, the first D/A converter module are connected with information source emitter, information source emitter The output signal of analog-to-digital conversion module is modulated to frequency band in oscillating fashion, and relaying subsystem is sent to by information source antenna System.
3. according to claim 1 a kind of while co-channel full duplex communication for coordination phase noise reduction system, its feature It is：The trunk subsystem includes relay reception antenna, the relay reception antenna, for receiving from source sub-system Signal, its output end are connected with link receiver, and reception signal is modulated to base band by link receiver in oscillating fashion, relaying The signal of receiver output is successively through the first analog-to-digital conversion module, first circulation prefix removal module, the first Fourier transformation mould Block, digital suppression module, amplification module, the second pilot tone insertion module, the second inverse Fourier transform module, second circulation prefix are inserted Enter module, after the second D/A converter module is handled, obtained signal is transferred to retransmitter, retransmitter is by The output signal of two analog-to-digital conversion modules is modulated to frequency band in oscillating fashion, and is sent a signal to by repeat transmitted antenna Receiving subsystem；The output end of the second pilot tone insertion module is also connected with digital suppression module.
4. according to claim 1 a kind of while co-channel full duplex communication for coordination phase noise reduction system, its feature It is：Described terminal receiving subsystem includes terminal antenna, and the terminal antenna receives the signal from retransmitter, eventually The output end of end antenna is connected with terminal receiver, and interruption receiver is with mode of oscillation by the signal modulation received to base band Afterwards, the second analog-to-digital conversion module, second circulation prefix removal module and the second fourier transformation module are transmitted a signal to successively Handled, the signal finally given by the output of the second fourier transformation module.
5. a kind of while co-channel full duplex communication for coordination phase noise inhibition method, it is characterised in that：Comprise the following steps：
S1. source sub-system produces echo signal and caused signal is handled, during the signal that processing obtains is sent to Stepchild's system；
S2. trunk subsystem receives the signal from source sub-system, is transformed into numeric field and carries out reception processing, in numeric field Digital suppression is carried out to the signal of reception processing, the signal after suppressing to numeral be amplified and pilot tone insertion after, next Cycle carries out transmitting processing, and is transformed into analog domain, transmits a signal to receiving subsystem；
S3. terminal receiving subsystem is handled the signal from trunk subsystem, completes the reception of signal.
6. according to claim 5 a kind of while co-channel full duplex communication for coordination phase noise inhibition method, its feature It is：The step S1 includes following sub-step：
S101. in m-th of communication cycle, it is P that information source, which produces transmission power,sSignal
S102. N is uniformly inserted in signal caused by information sourcepIndividual pilot tone, the location sets of pilot tone are：
<mrow> <mi>D</mi> <mo>=</mo> <mo>{</mo> <mi>i</mi> <mfrac> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> <mrow> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>|</mo> <mi>i</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>-</mo> <mn>1</mn> <mo>}</mo> <mo>,</mo> </mrow>
S103. the rear signal for inserting pilot tone is subjected to inverse Fourier transform, insertion cyclic prefix and digital-to-analogue conversion successively, treated The analog signal of transmitting is simultaneously modulated to frequency domain with mode of oscillation, is sent to trunk subsystem.
7. according to claim 5 a kind of while co-channel full duplex communication for coordination phase noise inhibition method, its feature It is：The step S2 includes following sub-step：
S201. trunk subsystem is in oscillating fashion by the signal modulation received to base band, and carries out modulus to baseband signal Change, go after cyclic prefix and Fourier transformation；
S202. the signal obtained to Fourier transformation carries out digital suppression, and the signal after numeral is suppressed is amplified and pilot tone Insertion, and the basis that the pilot signal of insertion is suppressed as next cycle numeral；
S203. by after the signal delay a cycle obtained in step S202, through inverse Fourier transform to time domain, then signal is inserted Enter cyclic prefix, carry out digital-to-analogue conversion, frequency band is modulated to mode of oscillation, is sent to terminal receiving subsystem.
8. according to claim 5 a kind of while co-channel full duplex communication for coordination phase noise inhibition method, its feature It is：The step S3 includes following sub-step：Receiving terminal system with mode of oscillation by the signal modulation received to base band, And successively through analog-to-digital conversion, go cyclic prefix and Fourier transformation after, the m+1 cycle has obtained the signal of final reception
9. according to claim 5 a kind of while co-channel full duplex communication for coordination phase noise inhibition method, its feature It is：The step S202 includes following sub-step：
The signal received using m-th of cycle of trunk subsystem in pilot frequency locations, except above a cycle trunk subsystem exists The signal of pilot frequency locations insertion, obtains the sub-carrier channels estimate of pilot frequency locations, using these estimates with linear interpolation Mode obtains the estimation of other subcarrier upper signal channels, from trunk subsystem in the signal that m-th of cycle receives, subtracts one The transmission signal of individual cycle delay and the product of estimation channel, to realize that the self-interference of current period numeral suppresses；
Calculate the optimal value of the sub-carrier power of transmission signal in trunk subsystem
The power amplification of the signal of k-th of subcarrier is arrived after self-interference numeral is suppressedSignal after being amplified；
I-th of position insertion power of signal after amplification is PtPilot tone, i ∈ D, using obtained signal as next cycle The basis that numeral suppresses.
10. according to claim 9 a kind of while co-channel full duplex communication for coordination phase noise inhibition method, its feature It is：The optimal value of the sub-carrier power of transmission signal in trunk subsystemCalculation procedure is as follows：
The first step, the value of parameter in setting system, including the total power constraint P at information sourceS, the total power constraint P at relayingR；If Put initialization tag value flag=0 and point=0；
Second step, judge whether to meet condition：PS≤Np·PsOr PR≤Np·Pt, if so, terminating and resetting system ginseng Value；If it is not, then enter the 3rd step；
3rd step, calculate the signal power on information source useful signal subcarrier
4th step, calculateZero point PR0If PR0＞ PR-NpPt, by PR0Again it is entered as PR- NpPt, wherein：
<mrow> <msub> <mi>L</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>x</mi> <mo>)</mo> </mrow> <mo>=</mo> <msqrt> <mfrac> <mrow> <mo>(</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> <mi>x</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>+</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> <mo>)</mo> <mo>(</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mi>x</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>+</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>s</mi> <mi>r</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> <mo>)</mo> </mrow> <mrow> <mo>(</mo> <msub> <mi>&amp;epsiv;</mi> <mi>k</mi> </msub> <mo>-</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> </msub> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>)</mo> <mo>(</mo> <mi>&amp;lambda;</mi> <mo>-</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <mo>)</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> </mrow> </mfrac> </msqrt> <mo>,</mo> </mrow>
<mrow> <mi>&amp;lambda;</mi> <mo>=</mo> <mfrac> <mn>1</mn> <msubsup> <mi>N</mi> <mi>c</mi> <mn>2</mn> </msubsup> </mfrac> <mo>{</mo> <mn>2</mn> <mfrac> <mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mn>4</mn> <mi>&amp;pi;</mi> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>/</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> </mrow> </msup> <mo>-</mo> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mn>4</mn> <mi>&amp;pi;</mi> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>/</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> </mrow> </msup> <mo>+</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> </mrow> <msup> <mrow> <mo>(</mo> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mn>4</mn> <mi>&amp;pi;</mi> <mi>&amp;Delta;</mi> <mi>f</mi> <mo>/</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> </mrow> </msup> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mn>2</mn> </msup> </mfrac> <mo>-</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>}</mo> <mo>,</mo> </mrow>
<mrow> <msub> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msub> <mo>=</mo> <mrow> <mo>(</mo> <mo>(</mo> <mrow> <mn>1</mn> <mo>-</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> </mrow> <mo>)</mo> <msub> <mi>P</mi> <mi>s</mi> </msub> <mo>+</mo> <mfrac> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mi>&amp;lambda;</mi> <mo>)</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </mfrac> <msub> <mi>P</mi> <mrow> <mi>s</mi> <mi>p</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>,</mo> </mrow>
εkIt is the mean-square value of the evaluated error of k-th of self-interference channel：
Wherein, A (i) represents sequenceNcI-th of value of point DFT transform, Δ f are mutually to make an uproar three dB bandwidth with carrying The ratio of ripple frequency domain interval, ErrAnd EsrRepresent the energy of channel impulse response, Hsr[k] and Hrd[k] represents from information source in respectively After and from be relayed to destination node channel impulse response DFT change k-th of value；
5th step, outage threshold ε=10 are set-7, arrange parameter a1=0, b1=PR0, lef=a1+0.382·(b1-a1), rig =a1+0.618·(b1-a1)；
6th step, judge whether to meet flag=0, if so, arrange parameter Pr=lef；If it is not, arrange parameter Pr=rig；
7th step, is proceeded as follows：
A, the power distribution P on each useful signal subcarrier at relaying is initialized0, P0It is a 1 × NcVector：
<mrow> <msub> <mi>P</mi> <mn>0</mn> </msub> <mo>=</mo> <mfrac> <msub> <mi>P</mi> <mi>r</mi> </msub> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> </mrow> </mfrac> <mo>&amp;lsqb;</mo> <mn>1</mn> <mo>,</mo> <mn>1</mn> <mo>,</mo> <mo>...</mo> <mo>,</mo> <mn>1</mn> <mo>&amp;rsqb;</mo> <mo>,</mo> </mrow>
InitializationIteration stopping condition η is seto=10-2, ηi=10-5, t=4096, increment times mu=10 and total Constraint number m=2 (Nc-Np)；
B, set withFor the function of variableFor：
<mrow> <mover> <mi>r</mi> <mo>^</mo> </mover> <mo>=</mo> <mfrac> <mi>t</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>k</mi> <mo>&amp;NotElement;</mo> <mi>D</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>e</mi> <mi>r</mi> <mi>f</mi> <mi>c</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </msqrt> <mo>)</mo> </mrow> <mo>-</mo> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>k</mi> <mo>&amp;NotElement;</mo> <mi>D</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>ln</mi> <mi> </mi> <msub> <mi>P</mi> <mi>k</mi> </msub> <mo>,</mo> </mrow>
Wherein,
<mrow> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> <mo>=</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>u</mi> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>P</mi> <mi>k</mi> </msub> <mo>+</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> <mo>)</mo> <mo>(</mo> <msub> <mi>&amp;epsiv;</mi> <mi>k</mi> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> </msub> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>+</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>s</mi> <mi>r</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> </mrow> <msub> <mi>P</mi> <mi>k</mi> </msub> </mfrac> <mo>)</mo> <mo>-</mo> <msub> <mi>E</mi> <mi>u</mi> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </mfrac> </mrow>
Eu[k] represents λ2Ps|Hrd[k]|2|Hsr[k]|2, FkThe inter-carrier interference energy on k-th of subcarrier at relaying is represented, its Estimate is
C, judge whether to meetIf it is not, into step F, if so, carrying out following steps：
(1) ifAccording to γFD[k] is on P0,kDerivativeCalculate functionOn P0,kDerivative, P0,kRepresent P0K-th of element value：
<mrow> <mfrac> <mrow> <mo>&amp;part;</mo> <mover> <mi>r</mi> <mo>^</mo> </mover> </mrow> <mrow> <mo>&amp;part;</mo> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mrow> </mfrac> <mo>=</mo> <mo>-</mo> <mfrac> <mi>t</mi> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mfrac> <mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </msup> <mo>&amp;part;</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> <mo>/</mo> <mo>&amp;part;</mo> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mrow> <mrow> <mn>2</mn> <msqrt> <mrow> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </msqrt> </mrow> </mfrac> <mo>-</mo> <mfrac> <mn>1</mn> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mfrac> <mo>,</mo> </mrow>
Then, according to γFD[k] is on P0,kSecond dervativeCalculate functionOn P0,kSecond dervative：
<mrow> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mover> <mi>r</mi> <mo>^</mo> </mover> </mrow> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </msup> <mi>t</mi> <mrow> <mo>(</mo> <mn>2</mn> <mo>+</mo> <msubsup> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msubsup> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> <mo>)</mo> </mrow> <mo>&amp;CenterDot;</mo> <msup> <mrow> <mo>(</mo> <mo>&amp;part;</mo> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> <mo>/</mo> <mo>&amp;part;</mo> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>)</mo> </mrow> <mn>2</mn> </msup> <mo>-</mo> <mn>2</mn> <mi>t</mi> <mo>&amp;CenterDot;</mo> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> <mo>/</mo> <mo>&amp;part;</mo> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mrow> <mrow> <mn>8</mn> <msqrt> <mrow> <msub> <mi>&amp;gamma;</mi> <mrow> <mi>F</mi> <mi>D</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </msqrt> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <mo>+</mo> <mfrac> <mn>1</mn> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mfrac> <mo>,</mo> </mrow>
Median w is calculated again：
<mrow> <mi>w</mi> <mo>=</mo> <msup> <mrow> <mo>(</mo> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>k</mi> <mo>&amp;NotElement;</mo> <mi>D</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mfrac> <mn>1</mn> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mover> <mi>r</mi> <mo>^</mo> </mover> <mo>/</mo> <mo>&amp;part;</mo> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>)</mo> </mrow> <mrow> <mo>-</mo> <mn>1</mn> </mrow> </msup> <munderover> <mi>&amp;Sigma;</mi> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>k</mi> <mo>&amp;NotElement;</mo> <mi>D</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mfrac> <mrow> <mo>&amp;part;</mo> <mover> <mi>r</mi> <mo>^</mo> </mover> <mo>/</mo> <mo>&amp;part;</mo> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mrow> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mover> <mi>r</mi> <mo>^</mo> </mover> <mo>/</mo> <mo>&amp;part;</mo> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>,</mo> </mrow>
P is calculated according to median w0 The updated value increment of individual element；
<mrow> <msub> <mi>&amp;Delta;P</mi> <mi>k</mi> </msub> <mo>=</mo> <mo>-</mo> <mfrac> <mrow> <mo>&amp;part;</mo> <mover> <mi>r</mi> <mo>^</mo> </mover> <mo>/</mo> <mo>&amp;part;</mo> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>+</mo> <mi>w</mi> </mrow> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mover> <mi>r</mi> <mo>^</mo> </mover> <mo>/</mo> <mo>&amp;part;</mo> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>,</mo> </mrow>
As k ∈ D, by Δ PkIt is entered as 0.To P0It is updated, the P after renewal0Equal to the P before renewal0Plus 0.1 times of Δ P, Wherein, Δ P represents Δ PkThe vector of composition；
According to w pairs of medianIt is updated, after renewalBe：
<mrow> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>k</mi> <mo>&amp;NotElement;</mo> <mi>D</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <msubsup> <mi>&amp;Delta;P</mi> <mi>k</mi> <mn>2</mn> </msubsup> <mo>&amp;CenterDot;</mo> <mfrac> <mrow> <msup> <mo>&amp;part;</mo> <mn>2</mn> </msup> <mover> <mi>r</mi> <mo>^</mo> </mover> </mrow> <mrow> <mo>&amp;part;</mo> <msubsup> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> <mn>2</mn> </msubsup> </mrow> </mfrac> </mrow>
(2) if after judging renewalWhether meetIf so, enter operation D, if it is not, return to step (1)；
D, t values are updated, after renewal t be equal to renewal before t be multiplied by mu；
E, judge whether the t after renewal meetsIf so, enter operation F, if it is not, returning to operation C；
F, median mw [k] is calculated：
<mrow> <mi>m</mi> <mi>w</mi> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> <mo>=</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mi>u</mi> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mrow> <mo>(</mo> <mi>&amp;lambda;</mi> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> <mo>+</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>r</mi> <mi>d</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> <mo>)</mo> <mo>(</mo> <msub> <mi>&amp;epsiv;</mi> <mi>k</mi> </msub> <mo>+</mo> <mfrac> <mrow> <msub> <mi>E</mi> <mrow> <mi>r</mi> <mi>r</mi> </mrow> </msub> <msub> <mi>F</mi> <mi>k</mi> </msub> <mo>+</mo> <msup> <mrow> <mo>|</mo> <mrow> <msub> <mi>H</mi> <mrow> <mi>s</mi> <mi>r</mi> </mrow> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> <mo>|</mo> </mrow> <mn>2</mn> </msup> <msub> <mover> <mi>P</mi> <mo>&amp;OverBar;</mo> </mover> <mi>s</mi> </msub> <mo>+</mo> <msub> <mi>E</mi> <mi>n</mi> </msub> </mrow> <msub> <mi>P</mi> <mrow> <mn>0</mn> <mo>,</mo> <mi>k</mi> </mrow> </msub> </mfrac> <mo>)</mo> <mo>-</mo> <msub> <mi>E</mi> <mi>u</mi> </msub> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </mfrac> </mrow>
According to median mw [k], bit error rate r is calculated：
<mrow> <mi>r</mi> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mn>2</mn> <mrow> <mo>(</mo> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <msub> <mi>N</mi> <mi>p</mi> </msub> <mo>)</mo> </mrow> </mrow> </mfrac> <munderover> <mo>&amp;Sigma;</mo> <mrow> <mi>k</mi> <mo>=</mo> <mn>0</mn> <mo>,</mo> <mi>k</mi> <mo>&amp;NotElement;</mo> <mi>D</mi> </mrow> <mrow> <msub> <mi>N</mi> <mi>c</mi> </msub> <mo>-</mo> <mn>1</mn> </mrow> </munderover> <mi>e</mi> <mi>r</mi> <mi>f</mi> <mi>c</mi> <mrow> <mo>(</mo> <msqrt> <mrow> <mi>m</mi> <mi>w</mi> <mo>&amp;lsqb;</mo> <mi>k</mi> <mo>&amp;rsqb;</mo> </mrow> </msqrt> <mo>)</mo> </mrow> </mrow>
Subsequently determine whether flag is 0, if so, by P0It is assigned to Pf and r is assigned to Bf, otherwise, then by P0It is assigned to Pg and incites somebody to action R is assigned to Bg；
8th step, judge whether point is equal to 0, if so, point and flag are entered as into 1 again and return to the 6th step；If Otherwise the 9th step is entered；
9th step, judge whether to meet | Bf-Bg | ＜ ε, if so, returning to optimal value groups of the Pf as relaying sub-carrier power distribution Into vector, calculating terminates；If it is not, into the tenth step；
Tenth step, compares Bf and Bg：
If Bf ＞ Bg, make a1=lef, lef=rig, rig=a1+0.618·(b1-a1), flag=1, return to the 6th step；
If Bf≤Bg, make b1=rig, rig=lef, lef=a1+0.382·(b1-a1), flag=0, return to the 6th step.
CN201710843738.8A 2017-09-19 2017-09-19 Phase noise suppression system and method for simultaneous same-frequency full duplex cooperative communication Active CN107659527B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710843738.8A CN107659527B (en) 2017-09-19 2017-09-19 Phase noise suppression system and method for simultaneous same-frequency full duplex cooperative communication

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710843738.8A CN107659527B (en) 2017-09-19 2017-09-19 Phase noise suppression system and method for simultaneous same-frequency full duplex cooperative communication

Publications (2)

Publication Number Publication Date
CN107659527A true CN107659527A (en) 2018-02-02
CN107659527B CN107659527B (en) 2020-09-08

Family

ID=61130666

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710843738.8A Active CN107659527B (en) 2017-09-19 2017-09-19 Phase noise suppression system and method for simultaneous same-frequency full duplex cooperative communication

Country Status (1)

Country Link
CN (1) CN107659527B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109104255A (en) * 2018-07-11 2018-12-28 东南大学 A kind of hardware corrected system of extensive broadband channel
CN109167744A (en) * 2018-11-06 2019-01-08 上海事凡物联网科技有限公司 A kind of phase noise combined estimation method
CN110602801A (en) * 2019-08-09 2019-12-20 北京展讯高科通信技术有限公司 Link configuration method and device
CN110932739A (en) * 2019-12-20 2020-03-27 成都大学 System and method for reducing error interference of communication and radar excitation signals
CN111082806A (en) * 2019-12-26 2020-04-28 北京化工大学 Method and system for eliminating noise
CN111211870A (en) * 2020-01-06 2020-05-29 华南理工大学 Iteration termination decision method in turbo receiving system
CN111726306A (en) * 2020-05-11 2020-09-29 北京大学 Full duplex system phase noise suppression method based on two-stage adaptive filtering

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1325198A (en) * 2000-05-22 2001-12-05 美国电报电话公司 Multi-input multi-output orthogonal frequency-division multiplexing system
US20070025461A1 (en) * 2005-07-29 2007-02-01 Samsung Electronics Co., Ltd. Phase noise compensation apparatus and an OFDM system having the apparatus and method thereof
CN101136731A (en) * 2007-08-09 2008-03-05 复旦大学 Method for eliminating phase noise using continuous transmission parameter signalling
CN102724027A (en) * 2012-06-15 2012-10-10 西安电子科技大学 Asynchronous space-time code coding/decoding system and method in full-duplex cooperative communication system
US20140198688A1 (en) * 2013-01-17 2014-07-17 Broadcom Corporation Method and Apparatus for Reducing Self Interference
CN104811213A (en) * 2014-01-28 2015-07-29 华为技术有限公司 Self-interference signal removing device and method
CN105187115A (en) * 2015-09-30 2015-12-23 西安电子科技大学 Orthogonal frequency division multiplexing (OFDM) co-frequency co-time full duplex relaying method
CN106101045A (en) * 2016-06-03 2016-11-09 北京邮电大学 A kind of OFDM full duplex based on phase noise Gauss albefaction polarization self-interference removing method
CN107154818A (en) * 2017-04-07 2017-09-12 西安电子科技大学 Co-channel full duplex bi-directional relaying transmission method while based on single carrier frequency domain equalization

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1325198A (en) * 2000-05-22 2001-12-05 美国电报电话公司 Multi-input multi-output orthogonal frequency-division multiplexing system
US20070025461A1 (en) * 2005-07-29 2007-02-01 Samsung Electronics Co., Ltd. Phase noise compensation apparatus and an OFDM system having the apparatus and method thereof
CN101136731A (en) * 2007-08-09 2008-03-05 复旦大学 Method for eliminating phase noise using continuous transmission parameter signalling
CN102724027A (en) * 2012-06-15 2012-10-10 西安电子科技大学 Asynchronous space-time code coding/decoding system and method in full-duplex cooperative communication system
US20140198688A1 (en) * 2013-01-17 2014-07-17 Broadcom Corporation Method and Apparatus for Reducing Self Interference
CN104811213A (en) * 2014-01-28 2015-07-29 华为技术有限公司 Self-interference signal removing device and method
CN105187115A (en) * 2015-09-30 2015-12-23 西安电子科技大学 Orthogonal frequency division multiplexing (OFDM) co-frequency co-time full duplex relaying method
CN106101045A (en) * 2016-06-03 2016-11-09 北京邮电大学 A kind of OFDM full duplex based on phase noise Gauss albefaction polarization self-interference removing method
CN107154818A (en) * 2017-04-07 2017-09-12 西安电子科技大学 Co-channel full duplex bi-directional relaying transmission method while based on single carrier frequency domain equalization

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
XINQUAN等: "Impacts of Phase Noise on Digital Self-Interference Cancellation in Full-Duplex Communications", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 *
张志亮: "同时同频全双工数字自干扰抑制关键技术", 《中国博士学位论文全文数据库信息科技辑》 *
文翔: "同时同频全双工自适应数字自干扰抑制技术研究与验证", 《中国优秀硕士学位论文全文数据库信息科技辑》 *
赵珺洁等: "OFDM系统原理及关键技术", 《通信与信息技术》 *

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109104255A (en) * 2018-07-11 2018-12-28 东南大学 A kind of hardware corrected system of extensive broadband channel
CN109104255B (en) * 2018-07-11 2021-01-26 东南大学 Hardware correction system of large-scale broadband channel
CN109167744A (en) * 2018-11-06 2019-01-08 上海事凡物联网科技有限公司 A kind of phase noise combined estimation method
CN109167744B (en) * 2018-11-06 2021-05-14 上海事凡物联网科技有限公司 Phase noise joint estimation method
CN110602801A (en) * 2019-08-09 2019-12-20 北京展讯高科通信技术有限公司 Link configuration method and device
CN110932739A (en) * 2019-12-20 2020-03-27 成都大学 System and method for reducing error interference of communication and radar excitation signals
CN110932739B (en) * 2019-12-20 2021-05-18 成都大学 System and method for reducing error interference of communication and radar excitation signals
CN111082806A (en) * 2019-12-26 2020-04-28 北京化工大学 Method and system for eliminating noise
CN111211870A (en) * 2020-01-06 2020-05-29 华南理工大学 Iteration termination decision method in turbo receiving system
CN111726306A (en) * 2020-05-11 2020-09-29 北京大学 Full duplex system phase noise suppression method based on two-stage adaptive filtering

Also Published As

Publication number Publication date
CN107659527B (en) 2020-09-08

Similar Documents

Publication Publication Date Title
CN107659527A (en) A kind of while co-channel full duplex communication for coordination phase noise reduction system and method
CN105978602B (en) The apparatus and method that a kind of while co-channel full duplex Nonlinear perturbations inhibit
CN104852752B (en) The system and method eliminated for the self-interference in efficient full-duplex communication
CN102487367B (en) Adaptive amplifying digital baseband pre-distortion method
CN105656830B (en) Ofdm signal method for inhibiting peak-to-average ratio based on distributed implementation
CN103338172A (en) Simultaneous common-frequency full duplex self-interference offset method in multi-path environment
CN103457638A (en) Restraining device and restraining method for burst impulse noise of power line communication channel
CN107026673A (en) The analog domain interference cancellation method and device of a kind of digital assistant
CN202085164U (en) Baseband digitalization apparatus for improving efficiency of power amplifier in adaptive ICS repeater
CN103166897B (en) The method of estimation of channel and IQI parameter in a kind of ofdm system
CN104935537A (en) Joint estimation compensation improving method of transmitter-receiver IQ imbalance and channel
CN102694562A (en) Method of adaptive interference cancellation by utilizing improved variable step size NLMS algorithm
CN103166879B (en) A kind of channel estimation methods and application the method ofdm system
CN115001913B (en) Full-duplex frequency domain self-interference elimination method based on digital assistance
CN107483095A (en) A kind of extensive nonopiate multi-access method assisted based on more relayings
CN105516044A (en) Orthogonal frequency division multiplexing (OFDM) system peak-to-average power ratio (PAPR) suppression method based on differential evolution
CN103384227A (en) Blind detection method for phase of partial transmit sequence of joint channel estimation
CN201947373U (en) High efficiency digital television transmitting device
CN103957178A (en) Multi-channel digital pre-distortion processing method and system
Brask et al. Deep digital signal recovery for LEO satellite communication in presence of system perturbations
CN108173611A (en) A kind of EVM test optimization methods based on OFDM system satellite repeaters
CN107171989A (en) Channel estimation methods based on DFT in visible light communication system
CN113852586A (en) Signal processing method and device based on 5G communication
Ibrahim et al. OFDM over wideband ionospheric HF channel: Channel modelling & optimal subcarrier power allocation
Singhal et al. Analysis of carrier aggregated OFDM signals in presence of dual band power amplifiers

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