CN103929212B - A kind of broadband receiver nonlinear blind identification and compensation method - Google Patents
A kind of broadband receiver nonlinear blind identification and compensation method Download PDFInfo
- Publication number
- CN103929212B CN103929212B CN201410160564.1A CN201410160564A CN103929212B CN 103929212 B CN103929212 B CN 103929212B CN 201410160564 A CN201410160564 A CN 201410160564A CN 103929212 B CN103929212 B CN 103929212B
- Authority
- CN
- China
- Prior art keywords
- signal
- distortion
- component
- compensation
- receiver
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Noise Elimination (AREA)
Abstract
This invention one broadband receiver nonlinear blind identification and compensation method, relate to radio monitoring field, particularly relate to broadband receiver nonlinear blind identification and compensation method.This invention is by compensating distorted signal, obtain mixed signal, calculate the cost function of distortion component short-time energy information in mixed signal again, calculate the compensating parameter making distortion component minimum, then the compensating parameter in compensation model is updated, repeat this step until distortion component is preferably suppressed, thus there is the SFDR of raising receiver and the effect of the ability of extraction weak signal.
Description
Technical field
The invention belongs to radio monitoring field, particularly relate to broadband receiver nonlinear blind identification and compensation method.
Background technology
Relative to transmitter system strong nonlinearity that may be present, the nonlinear characteristic of receiver is more weak.After receiver receives signal,
Through analog devices such as amplifier, wave filter, frequency mixers, signal can produce non-linear distortion through these analog devices.ADC
(analog-digital converter), as a kind of mixed processing device, is also the nonlinear important sources of receiver link, and signal passes through ADC
Can produce distortion in time domain, frequency domain produces harmonic distortion, intermodulation distortion etc..Reducing receiver one of them method nonlinear is choosing
With speed, conversion accuracy is higher, the ADC that SFDR is bigger.But, the price of high performance ADC is very
Expensive.Now, research to the linearization technique of receiver both at home and abroad is less.In engineering, in order to ensure the linearity of receiver,
One is the thought that can use back-off, will input signal power limited below 1dB compression point, only utilize receiver
Linear gain region.The shortcoming of Power back is to have had a strong impact on the dynamic range of receiving of receiver, sacrifices efficiency of transmission simultaneously.
Two is to use feedback technique, is extracted by the distortion component in nonlinear output signal, feeds back to input, the most again with input
Signal carries out anti-phase cancellation thus obtains the output signal after suppressing nonlinear component, but the method also exists stablizing of reponse system
Sex chromosome mosaicism, signal lag limits the maximum bandwidth allowing transmission signal present in transmission system simultaneously.Three is to use predistortion skill
Art, is i.e. that input signal carries out certain conversion in advance so that its conversion characteristics is just reciprocal with the nonlinear characteristic of receiver, from
And the non-linearity of cancellation receiver.But carrying out self-adapted pre-distortion before radio-frequency front-end, actual realization is upper exists bigger difficulty.
Current a kind of relatively conventional modeling compensation method is to look for table method.The method is to use the excitation of specific input signal to be compensated
System, by all possible output code value and it relative to the Difference Storage of preferable input code in look-up table, linearize in reality
During, using the instant amplitude value of output signal as index, search and obtain the corresponding margin of error, then by it from output signal
Deduct, reach to compensate the purpose of systematic error.But owing to frequency and the amplitude of input signal have bigger scope, compensate simultaneously
System be likely to be of memory effect, so for ensure compensate precision, need generate multidimensional look-up table various possible to index
Input condition.The most so, to ensure precision, then great memory capacity it is necessarily required to, simultaneously need to address faster and locate
Reason speed, the real-time that guarantee compensates.Therefore, the hardware implementation complexity of the method limits its application greatly.
Summary of the invention
It is an object of the invention to just weak point Curve guide impeller to background technology a kind of broadband receiver nonlinear blind identification and benefit
Compensation method, thus reach to improve the purpose of the ability of the SFDR of receiver and extraction weak signal.
The solution of the present invention is a kind of broadband receiver nonlinear blind identification and compensation method, and the method includes:
Step 1: receiver receives analog signal, this signal processes through receiver and distortion occurs, obtains distorted signal, and distortion is believed
Number it is made up of linear component and distortion component;
Step 2: build compensation model, and predesigned compensation parameter according to distorted signal, compensation model output compensates signal and believes distortion
Number compensate, obtain mixed signal;
Step 3: use FFT (Fast Fourier Transform (FFT)) to ask for the power spectral density function of mixed signal, at power spectral density function
In, judge linear component and the frequency of distortion component and bandwidth according to default received signal power thresholding;
Step 4: according to receiving signal and the frequency of distortion component and bandwidth, build the spectral response of comb filter, then carry out
Inverse FFT obtains filter coefficient, is filtered with compensation signal after filter coefficient windowing, calculates the humorous of linear component
Ripple item and intermodulation item;
Step 5: calculated the short-time energy of distortion component by the harmonic term of linear component and intermodulation item, reconstruct cost function,
Obtain the compensating parameter when cost function obtains minimum of a value;
Step 6: use adaptive algorithm that the compensating parameter preset in step 2 is replaced with the compensating parameter obtained in step 5, weight
Multiple above-mentioned steps, until this adaptive algorithm convergence during Ji Suaning, reaches to suppress the effect of distortion component.
The compensation model wherein built in step 2 is Volterra compensation model.
One broadband receiver nonlinear blind identification of the present invention and compensation method, by compensating distorted signal, obtain mixing letter
Number, then calculate the short-time energy cost function of distortion component in mixed signal, calculate the compensating parameter making distortion component minimum,
Then update the compensating parameter in compensation model, repeat this step until distortion component is preferably suppressed, thus there is raising and connect
The SFDR of receipts machine and the effect of the ability of extraction weak signal.
Accompanying drawing explanation
Fig. 1 is the flow chart of the method that inventive receiver blind discrimination compensates;
Fig. 2 is that the power threshold of linear component judges schematic diagram.
Detailed description of the invention
Distortion is there is after step 1: receiver receives analog signal, the process through receiver radio frequency front end of this signal and AD sampling,
Obtaining distorted signal, distorted signal is made up of linear component and distortion component.
Step 2: structure Volterra compensation model, the initial value ω of preset model compensating parameter ω0If, the rank of Volterra model
Number D=3, then compensation model output compensates signal is p (k)=v (k) ω0,
Wherein v (k)=[y2(k),y(k)y(k-1)…,y3(k),y2(k)y(k-1)…y3(k-1)],
ω0=[ω1,ω2,ω3,ω4,ω5,ω6,ω7]T=[1,0,0,0,0,0,0]TInitial value for model parameter, it may be assumed that
p(k)=ω1y2(k)+ω2y(k)y(k-1)+ω3y2(k-1)+ω4y3(k)
+ω5y2(k)y(k-1)+ω6y(k)y2(k-1)+ω7y3(k-1)
Then obtaining mixed signal s (k)=y (k)-p (k), wherein y (k) is distorted signal, and p (k) is for compensating signal.
Step 3: ask for the power spectral density of mixed signal s (k), and judge linear component and the frequency of distortion component and bandwidth:
Mixed signal s (k) being made FFT computing, is compensated the power spectral density function of signal s (k), previously given judgement is linear
Power threshold P of componentthreshold, according to linear component judgment criterion P > and Pthreshold, then can get linear component frequency in frequency spectrum
Rate f1,f2…fNWith bandwidth l1,l2…lN, frequency and bandwidth further according to linear component obtain harmonic component
2f1,3f1,…2fN,3fN... and intermodulation component f1+f2,…2fN-1+fN... frequency and bandwidth, as shown in Figure 2.
Step 4: according to frequency and the bandwidth of distortion component, structure comb filter calculates distortion component:
Construct the spectral response of comb filter according to harmonic wave and the frequency of intermodulation component and bandwidth, spectral response is obtained as IFFT
To filter coefficient, use FIR filter to be filtered after filter coefficient windowing with mixed signal, calculate linear component
Harmonic term and intermodulation item, if FIR filter vector is g=[g0,g1…gB-1]T, dimension is B × 1, and the most finally extract is mixed
The nonlinear component closing signal is expressed as:
sf(k,ω)=gTs(k)=gT[y(k)-VT(k)ω]
What wherein y (k) represented is the distorted signal after broadband receiver AD samples, and V (k) is to characterize to compensate mission nonlinear memory
The matrix of part.ω is the compensating parameter that the present invention wants identification,
y(k)=[y(k-B+1),y(k-B+2)…y(k)]T, V (k)=[v (k-B+1), v (k-B+2) ... v (k)],
v(k)=[y2(k),y(k)y(k-1)…,y3(k),y2(k)y(k-1)…y3(k-1)]T。
Step 5: harmonic term and intermodulation item according to linear component construct cost function J (ω), utilize adaptive algorithm to compensation
Model parameter is updated.
The core of the present invention is in the case of ensureing that linear component energy is constant, makes nonlinear component energy minimum, due to
sf(k,ω)=gTNot including linear component part in s (k), therefore making cost function is:
J(ω)=E[s2 f(k,ω)]
Recycling LMS adaptive algorithm (conventional is LMS or RLS algorithm) updates compensation model parameter ω so that cost
The value of function is minimum.
LMS adaptive approach:
Cost function J (ω) is launched and seeks gradient, can obtain:
According to Wiener filtering steepest descent method thought: Can obtain:
ω(n+1)=ω(n)-μ{E[V(n)ggTVT(n)]ω(n)-E[V(n)ggTy(n)]}
Instantaneous estimation value is utilized to approach assembly average, the LMS Adaptive Identification algorithm in available blind discrimination compensation system:
ω(n+1)=ω(n)-μ{V(n)ggTVT(n)ω(n)-V(n)ggTy(n)}。
Step 6: the ω substitution compensation model that adaptive algorithm is tried to achieve, thus mixed signal s (k) after being compensated:
The ω tried to achieve after adaptive algorithm being restrained substitutes in step 2, is compensated model output signal p (k), thus is mixed
Close signal s (k)=y (k)-p (k).
After blind discrimination compensates output, the distortion component in signal is effectively suppressed, and weak signal can be separated from distortion component
Out, improve the dynamic range of receiver, improve the nonlinear characteristic of its front end.
Claims (2)
1. the identification of broadband receiver nonlinear blind and a compensation method, the method includes:
Step 1: receiver receives analog signal, this signal processes through receiver and distortion occurs, obtains distorted signal, and distortion is believed
Number it is made up of linear component and distortion component;
Step 2: build compensation model, and predesigned compensation parameter according to distorted signal, compensation model output compensates signal and believes distortion
Number compensate, obtain mixed signal;
Step 3: use Fast Fourier Transform (FFT) to ask for the power spectral density function of mixed signal, in power spectral density function, root
Linear component and the frequency of distortion component and bandwidth is judged according to default received signal power thresholding;
Step 4: according to receiving signal and the frequency of distortion component and bandwidth, build the spectral response of comb filter, then carry out
Inverse FFT obtains filter coefficient, is filtered with mixed signal after filter coefficient windowing, calculates the humorous of linear component
Ripple item and intermodulation item;
Step 5: calculated the short-time energy of distortion component by the harmonic term of linear component and intermodulation item, reconstruct cost function,
Obtain the compensating parameter when cost function obtains minimum of a value;
Step 6: use adaptive algorithm that the compensating parameter preset in step 2 is replaced with the compensating parameter obtained in step 5, weight
Multiple above-mentioned steps, until this adaptive algorithm convergence during Ji Suaning, reaches to suppress the effect of distortion component.
2. a kind of broadband receiver nonlinear blind identification as claimed in claim 1 and compensation method, it is characterised in that in step 2
The compensation model built is Volterra compensation model.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410160564.1A CN103929212B (en) | 2014-04-21 | 2014-04-21 | A kind of broadband receiver nonlinear blind identification and compensation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410160564.1A CN103929212B (en) | 2014-04-21 | 2014-04-21 | A kind of broadband receiver nonlinear blind identification and compensation method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103929212A CN103929212A (en) | 2014-07-16 |
CN103929212B true CN103929212B (en) | 2016-08-24 |
Family
ID=51147308
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410160564.1A Active CN103929212B (en) | 2014-04-21 | 2014-04-21 | A kind of broadband receiver nonlinear blind identification and compensation method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103929212B (en) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107065524B (en) * | 2017-01-20 | 2020-05-12 | 东南大学 | Servo system frequency domain identification system and identification method |
CN111294295A (en) * | 2020-01-15 | 2020-06-16 | 中国科学院海洋研究所 | Multi-channel FIR filter ocean underwater sound channel blind identification algorithm based on subspace |
CN112187266B (en) * | 2020-09-29 | 2024-03-08 | 中国电子科技集团公司第三十六研究所 | Nonlinear correction method and device of analog-to-digital converter and electronic equipment |
CN115941402B (en) * | 2022-11-18 | 2024-07-26 | 西安电子科技大学 | Self-adaptive nonlinear equalization method for broadband receiver |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1612560A (en) * | 2003-10-28 | 2005-05-04 | 三星电子株式会社 | Receiver for compensating nonlinearly distorted multicarrier signals |
CN101188463A (en) * | 2007-12-20 | 2008-05-28 | 北京北方烽火科技有限公司 | A power calibration and real time correction device for CDMA receiver and its implementation method |
CN102257734A (en) * | 2008-12-19 | 2011-11-23 | 艾利森电话股份有限公司 | Interference cancellation in an OFDM receiver |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8139630B2 (en) * | 2005-10-03 | 2012-03-20 | Clariphy Communications, Inc. | High-speed receiver architecture |
US8467438B2 (en) * | 2010-08-02 | 2013-06-18 | Bassel F. Beidas | System and method for iterative nonlinear compensation for intermodulation distortion in multicarrier communication systems |
US8873971B2 (en) * | 2010-10-11 | 2014-10-28 | Nec Laboratories America, Inc. | Nonlinear compensation using an enhanced backpropagation method with subbanding |
-
2014
- 2014-04-21 CN CN201410160564.1A patent/CN103929212B/en active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1612560A (en) * | 2003-10-28 | 2005-05-04 | 三星电子株式会社 | Receiver for compensating nonlinearly distorted multicarrier signals |
CN101188463A (en) * | 2007-12-20 | 2008-05-28 | 北京北方烽火科技有限公司 | A power calibration and real time correction device for CDMA receiver and its implementation method |
CN102257734A (en) * | 2008-12-19 | 2011-11-23 | 艾利森电话股份有限公司 | Interference cancellation in an OFDM receiver |
Non-Patent Citations (2)
Title |
---|
大动态宽带干涉式雷电测向接收机新技术;彭亮,马洪,张辉;《武汉大学学报•信息科学版》;20111130;第36卷(第11期);1358-1362 * |
宽带接收机前端非线性的盲辨识与数字后补偿方法;彭隽,马洪等;《通信学报》;20110430;第32卷(第4期);119-125 * |
Also Published As
Publication number | Publication date |
---|---|
CN103929212A (en) | 2014-07-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3089414A1 (en) | Digital pre-distortion parameter obtaining method and pre-distortion system | |
US9231733B2 (en) | Method and system for updating multi-frequency-band pre-distortion coefficient lookup table | |
CN101272373B (en) | Self-adapting analog quadrature modulation disbalance compensation method and device | |
CN103929212B (en) | A kind of broadband receiver nonlinear blind identification and compensation method | |
CN103715992B (en) | Based on the power amplifier predistortion device and method simplifying Volterra progression | |
US20130329833A1 (en) | Modeling Digital Predistorter | |
CN104796364B (en) | A kind of pre-distortion parameters acquiring method and pre-distortion system | |
WO2017095869A1 (en) | Digital compensator | |
JP2009213113A (en) | Apparatus and method for identifying reverse characteristic of nonlinear system, power amplifier, and predistorter of power amplifier | |
CN105656434B (en) | Power amplifier digital predistortion apparatus and method based on modification piecewise linear function | |
CN103856429A (en) | Adaptive predistortion system and method based on hybrid indirect learning algorithm | |
CN104363191A (en) | Digital pre-distortion method of frequency hopping communication system | |
US20070133713A1 (en) | Power amplifier pre-distorter training | |
CN105635009A (en) | Self-adaptive MIMO pre-distortion method for hybrid compensation of multi-branch crosstalk and IQ imbalance | |
KR101386239B1 (en) | Predistorter for compensating of nonlinear distortion and method for the same | |
CN112019221B (en) | Signal processing method, device and storage medium | |
Liu et al. | On the robustness of look-up table digital predistortion in the presence of loop delay error | |
CN104092633A (en) | Multi-wavelet feedback blind equalization method for nonlinear satellite channel | |
CN101710887B (en) | Method and system for processing digital predistortion | |
CN105897189A (en) | Digital pre-distorter structure and method for look-up table | |
CN110971197B (en) | Digital predistortion device and method | |
Dvorak et al. | Adaptive-order polynomial methods for power amplifier model estimation | |
CN102983820A (en) | Nonlinear injection type linearization system and digital predistortion method | |
CN203800956U (en) | Digital pre-distortion system of off-line estimation pre-distortion coefficient | |
KR20090125597A (en) | Method for extracting nonlinear model parameter of wideband signal using narrowband signal, apparatus and method for digital predistortering its using |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right |
Effective date of registration: 20200630 Address after: 610000 No. 1301-1304, 13th floor, No. 2 Building, Group B, 200 Tianquan Road, Chengdu High-tech Zone, Sichuan Province Patentee after: CHENGDU DIANKE HUIAN TECHNOLOGY Co.,Ltd. Address before: 611731 Chengdu province high tech Zone (West) West source Avenue, No. 2006 Patentee before: UNIVERSITY OF ELECTRONIC SCIENCE AND TECHNOLOGY |
|
TR01 | Transfer of patent right |