CN105635008B - A kind of adaptive equilibrium method, device and adaptive equalizer - Google Patents
A kind of adaptive equilibrium method, device and adaptive equalizer Download PDFInfo
- Publication number
- CN105635008B CN105635008B CN201610203866.1A CN201610203866A CN105635008B CN 105635008 B CN105635008 B CN 105635008B CN 201610203866 A CN201610203866 A CN 201610203866A CN 105635008 B CN105635008 B CN 105635008B
- Authority
- CN
- China
- Prior art keywords
- current time
- related coefficient
- vector
- calculates
- denominator
- 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
Classifications
-
- 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/03012—Arrangements for removing intersymbol interference operating in the time domain
- H04L25/03019—Arrangements for removing intersymbol interference operating in the time domain adaptive, i.e. capable of adjustment during data reception
-
- 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
-
- 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/03878—Line equalisers; line build-out devices
- H04L25/03885—Line equalisers; line build-out devices adaptive
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Filters That Use Time-Delay Elements (AREA)
- Cable Transmission Systems, Equalization Of Radio And Reduction Of Echo (AREA)
Abstract
The invention discloses a kind of adaptive equilibrium method, device and adaptive equalizers.This method includes：Error signal calculation step, related coefficient calculate step, decorrelation vector calculation step, step factor and calculate step and tap weights vector calculating step.
Description
Technical field
The present invention relates to channel equalization technique field more particularly to a kind of adaptive equilibrium method, device and it is adaptive
Weighing apparatus.
Background technique
Adaptive equalization technique can be used for the receiving end of digital communication system, its effect is to reduce intersymbol interference and noise
Caused compound influence.
A kind of existing adaptive equilibrium method is based on least mean-square error (LMS) algorithm, the process of signal processing
As shown in Figure 1.Its existing main problem is：Condition of the rate of convergence by filter tap input vector autocorrelation matrix
Number is affected.For this problem, the relevant least mean-square error of time solution (DLMS) algorithm is proposed in the prior art, but
It is that in filter tap values vector ω (n) renewal process each time, DLMS method will calculate the input signal at current time
With the value of the decorrelation factor mu (n) of the decorrelation coefficient a (n) and current time of the input signal of last moment, in Project Realization
When need to use a large amount of multiplier and the operand of algorithm is larger.
Summary of the invention
The purpose of the present invention is to provide a kind of adaptive equilibrium method, device and adaptive equalizers, existing to solve
Above-mentioned technical problem present in technology.
An embodiment provides a kind of adaptive equilibrium methods, including：Error signal calculation step：According to
The desired signal at current time, the filter tap weight vector at current time and the input signal at current time calculate current time
Error signal；Related coefficient calculates step：According to the phase relation at the molecular computing current time of the related coefficient of last moment
The denominator of several molecules and the related coefficient that current time is calculated according to the denominator of the related coefficient of last moment；Decorrelation vector
Calculate step：Worked as according to the calculating of the related coefficient of the input signal at current time, the input signal of last moment and current time
The decorrelation vector at preceding moment；Step factor calculates step：According to the error signal at current time, the related coefficient at current time
Molecule and the related coefficient at current time denominator calculate current time step factor；And tap weights vector calculates step
Suddenly：According to the decorrelation vector at current time, the step factor meter of the filter tap weight vector at current time and current time
Calculate the filter tap weight vector of subsequent time.
Another embodiment of the present invention provides a kind of adaptive equalization devices, including：Error signal calculation module：With
Work as according to the calculating of the desired signal at current time, the filter tap weight vector at current time and the input signal at current time
The error signal at preceding moment；Related coefficient computing module：For it is current according to the molecular computing of the related coefficient of last moment when
The denominator of the molecule of the related coefficient at quarter and the related coefficient according to the denominator of the related coefficient of last moment calculating current time；
Decorrelation vector calculation module：For according to the input signal at current time, the input signal of last moment and current time
The decorrelation vector at related coefficient calculating current time；Step factor computing module：For according to the error signal at current time,
The denominator of the related coefficient of the molecule and current time of the related coefficient at current time calculates the step factor at current time；And
Tap weights vector computing module：For according to the decorrelation vector at current time, the filter tap weight vector at current time and
The step factor at current time calculates the filter tap weight vector of subsequent time.
Another embodiment of the invention provides a kind of adaptive equalizer, including：It is according to the present invention adaptive equal
Weigh device；And transversal filter, it is connected with adaptive equalization device, the tap weight value of transversal filter is by described adaptive
Balancer is answered to be calculated.
Detailed description of the invention
By reading the following detailed description of the preferred embodiment, various other advantages and benefits are common for this field
Technical staff will become clear.The drawings are only for the purpose of illustrating a preferred embodiment, and is not considered as to the present invention
Limitation.And throughout the drawings, the same reference numbers will be used to refer to the same parts.Wherein in the accompanying drawings, reference number
Alphabetic flag later indicates that multiple identical components will omit its last alphabetic flag when referring to these components.Attached
In figure：
Fig. 1 show the signal flow graph of the existing adaptive equalization technique based on LMS；
Fig. 2 is the flow chart of one embodiment of adaptive equilibrium method of the invention；
Fig. 3 show the schematic diagram of one embodiment of adaptive equalization device of the invention.
In the accompanying drawings, same or similar element is referred to using same or similar label.
Specific embodiment
Carry out detailed description of the present invention illustrative embodiments with reference to the drawings.It should be appreciated that shown in attached drawing and
The embodiment of description is only exemplary, it is intended that is illustrated the principle and spirit of the invention, and is not limited model of the invention
It encloses.
Fig. 2 is the flow chart of one embodiment 200 of adaptive equilibrium method of the invention.Embodiment 200 may include
Following steps 201 to 205.
Step 201 is error signal calculation step：According to the desired signal at current time, the filter tap at current time
Weight vector and the input signal at current time calculate the error signal at current time.
In one embodiment of the invention, can first by the initial value of filter tap weight vector be set as ω (0)=
0。
The error signal at current time can be calculated by following equation (1)：
E (n)=d (n)-ωH(n)u(n) (1)
Wherein, e (n) indicates the error signal at current time, and d (n) indicates the desired signal at current time, ωH(n) it indicates
The conjugate transposition of the filter tap weight vector at current time, u (n) indicate the input signal at current time.
Step 202 is that related coefficient calculates step：According to the molecular computing current time of the related coefficient of last moment
The denominator of the molecule of related coefficient and the related coefficient according to the denominator of the related coefficient of last moment calculating current time.
In one embodiment of the invention, related coefficient α (n) can be counted as composed by molecule and denominator.
In one embodiment of the invention, therein when indicating related coefficient α (n) with α (n)=fz (n)/fm (n)
Molecule fZ(n) it can be calculated by following formula (2), denominator f thereinm(n) it can be calculated by following formula (3)
It arrives：
fZ(n)=fZ(n-1)+u(n)u(n-1)-u(n-M)u(n-M+1) (2)
fm(n)=fm(n-1)+|u(n)|2-|u(n-M-1)|2 (3)
Wherein, fZIt (n-1) is the molecule of the related coefficient of last moment, M is the sum of filter, fmIt (n-1) is upper a period of time
The denominator of the related coefficient at quarter.
Can use it can be seen from above formula (2) and (3) the last related coefficient being calculated molecule and
Denominator obtains the molecule and denominator of the related coefficient at current time by lesser operand.
Step 203 is decorrelation vector calculation step：According to the input signal at current time, the input signal of last moment
The decorrelation vector at current time is calculated with the related coefficient at current time.
In one embodiment of the invention, decorrelation vector z (n) can be calculated using following formula (4)：
Z (n)=u (n)-α (n) u (n-1) (4)
Step 204 is that step factor calculates step：According to the error signal at current time, the related coefficient at current time
The step factor at the denominator of molecule and the related coefficient at current time calculating current time.
In one embodiment of the invention, it can use the input signal u that following equation (5) first calculates current time
(n) with the conjugate transposition u of their ownH(n) vector product fu(n)：
fu(n)=uH(n) u (n)=fu(n-1)+u(n)u(n)-u(n-M)u(n-M) (5)
Then the step factor μ (n) at current time is calculated further according to following equation (6).
Wherein, μ is preset step parameter, can be positive real number, e*It (n) is the error signal e at current time
(n) conjugation.
The present invention also utilizes upper a period of time in material calculation factor mu (n) it can be seen from above formula (5) and (6)
The some variate-values being calculated are carved, such as：fu(n-1), some variate-values being calculated in step 203 are also used, this
Sample just considerably reduces the operand in step factor calculating process.
Step 205 is that tap weights vector calculates step：According to the decorrelation vector at current time, the filter at current time
The filter tap weight vector of the step factor at tap weights vector and current time calculating subsequent time.
In one embodiment of the invention, the filter tap power of subsequent time can be obtained by following formula (7)
Vector
Wherein,Indicate the filter tap weight vector at current time.
So far adaptive equilibrium method according to an embodiment of the present invention is described.The method of the present invention enough reduces channel equalization mistake
Operand in journey.
The present invention also provides adaptive equalization devices.Adaptive equalization device of the invention is shown with reference to Fig. 3, Fig. 3
One embodiment 300 schematic diagram.Device 300 may include：Error signal calculation module 301：For according to current time
Desired signal, the filter tap weight vector at current time and the input signal at current time calculate the error letter at current time
Number；Related coefficient computing module 302：For the related coefficient according to molecular computing current time of the related coefficient of last moment
Molecule and according to the denominator of the related coefficient of last moment calculate current time related coefficient denominator；Decorrelation is to meter
Calculate module 303：Based on the related coefficient according to the input signal at current time, the input signal of last moment and current time
Calculate the decorrelation vector at current time；Step factor computing module 304：For according to the error signal at current time, it is current when
The denominator of the related coefficient of the molecule and current time of the related coefficient at quarter calculates the step factor at current time；And tap weights
Vector calculation module 305：For according to the decorrelation vector at current time, the filter tap weight vector at current time and current
The step factor at moment calculates the filter tap weight vector of subsequent time.
In one embodiment of the invention, related coefficient computing module 302 may further include：Point of related coefficient
Sub- computing module, for according to formula
fZ(n)=fZ(n-1)+u (n) u (n-1)-u (n-M) u (n-M+1) calculates the molecule f of the related coefficient at current timeZ
(n), wherein fZIt (n-1) is the molecule of the related coefficient of last moment, u (n) is the input signal at current time, and M is filter
Sum；And the denominator computing module of related coefficient, for according to formula fm(n)=fm(n-1)+|u(n)|2-|u(n-M-1)
|2Calculate the denominator f of the related coefficient at current timem(n), wherein fmIt (n-1) is the denominator of the related coefficient of last moment.
In one embodiment of the invention, step factor computing module 304 may further include：For calculating
fu(n)=fu(n-1) module of+u (n) u (n)-u (n-M) u (n-M)；And for according to formulaCalculate the module of the step factor μ (n) at current time, wherein uH(n) turn for the conjugation of u (n)
It sets, μ is preset step parameter, and e* (n) is the conjugation of the error signal e (n) at current time.
So far adaptive equalization device according to an embodiment of the present invention is described.It is similar with above-mentioned method, present invention dress
Set the operand during enough reduction channel equalization.
Claims (3)
1. a kind of adaptive equilibrium method, characterized in that including：
Error signal calculation step：According to the desired signal at current time, the filter tap weight vector at current time and current
The input signal at moment calculates the error signal at current time：
E (n)=d (n)-ωH(n)u(n)
Wherein, e (n) indicates the error signal at current time, and d (n) indicates the desired signal at current time, ωH(n) indicate current
The conjugate transposition of the filter tap weight vector at moment, u (n) indicate the input signal at current time；
Related coefficient calculates step：According to the molecule of the related coefficient at the molecular computing current time of the related coefficient of last moment
And the denominator of the related coefficient at current time is calculated according to the denominator of the related coefficient of last moment；
Decorrelation vector calculation step：According to the input signal at current time, the input signal of last moment and current time
The decorrelation vector at related coefficient calculating current time：
Z (n)=u (n)-α (n) u (n-1)
Wherein, z (n) is decorrelation vector, and α (n) is related coefficient；
Step factor calculates step：According to the error signal at current time, the related coefficient at current time molecule and it is current when
The denominator of the related coefficient at quarter calculates the step factor at current time；And
Tap weights vector calculates step：According to the decorrelation vector at current time, the filter tap weight vector at current time and
The step factor at current time calculates the filter tap weight vector of subsequent time：
Wherein,Indicate the filter tap weight vector at current time, μ (n) is step factor；
The related coefficient calculates step：
According to formula：
fZ(n)=fZ(n-1)+u(n)u(n-1)-u(n-M)u(n-M+1)
Calculate the molecule f of the related coefficient at current timeZ(n), wherein fZIt (n-1) is the molecule of the related coefficient of last moment, u
It (n) is the input signal at current time, M is the sum of filter；And
According to formula fm(n)=fm(n-1)+|u(n)|2-|u(n-M-1)|2Calculate the denominator f of the related coefficient at current timem
(n), wherein fmIt (n-1) is the denominator of the related coefficient of last moment；
The step factor calculates step：
It calculates
fu(n)=fu(n-1)+u(n)u(n)-u(n-M)u(n-M)；And
According to formulaCalculate the step factor μ (n) at current time, wherein μ is preset step
Long parameter, e*It (n) is the conjugation of the error signal e (n) at current time, fu(n) for current time input signal u (n) and its from
Oneself conjugate transposition uH(n) vector product, fu(n-1) for last moment input signal and their own conjugate transposition to
Measure product.
2. a kind of adaptive equalization device, characterized in that including：
Error signal calculation module：For according to the desired signal at current time, the filter tap weight vector at current time and
The input signal at current time calculates the error signal at current time：
E (n)=d (n)-ωH(n)u(n)
Wherein, e (n) indicates the error signal at current time, and d (n) indicates the desired signal at current time, ωH(n) indicate current
The conjugate transposition of the filter tap weight vector at moment, u (n) indicate the input signal at current time；
Related coefficient computing module：For according to the related coefficient at molecular computing current time of the related coefficient of last moment
The denominator of molecule and the related coefficient according to the denominator of the related coefficient of last moment calculating current time；
Decorrelation vector calculation module：For according to the input signal at current time, the input signal of last moment and it is current when
The related coefficient at quarter calculates the decorrelation vector at current time：
Z (n)=u (n)-α (n) u (n-1)
Wherein, z (n) is decorrelation vector, and α (n) is related coefficient；
Step factor computing module：For according to the error signal at current time, the molecule of the related coefficient at current time and working as
The denominator of the related coefficient at preceding moment calculates the step factor at current time；And
Tap weights vector computing module：For decorrelation vector, the filter tap at current time according to current time weigh to
The step factor at amount and current time calculates the filter tap weight vector of subsequent time：
Wherein,Indicate the filter tap weight vector at current time, μ (n) is step factor；The related coefficient calculates mould
Block further comprises：
The molecular computing module of related coefficient, for according to formula fZ(n)=fZ(n-1)+u(n)u(n-1)-u(n-M)u(n-M+
1) the molecule f of the related coefficient at current time is calculatedZ(n), wherein fZIt (n-1) is the molecule of the related coefficient of last moment, u
It (n) is the input signal at current time, M is the sum of filter；And
The denominator computing module of related coefficient, for according to formula fm(n)=fm(n-1)+|u(n)|2-|u(n-M-1)|2Calculating is worked as
The denominator f of the related coefficient at preceding momentm(n), wherein fmIt (n-1) is the denominator of the related coefficient of last moment；
The step factor computing module further comprises：
For calculating
fu(n)=fu(n-1) module of+u (n) u (n)-u (n-M) u (n-M)；And
For according to formulaCalculate the module of the step factor μ (n) at current time；
Wherein, μ is preset step parameter, e*It (n) is the conjugation of the error signal e (n) at current time, fuIt (n) is current
The input signal u (n) at the moment and conjugate transposition u of their ownH(n) vector product, fu(n-1) believe for the input of last moment
Vector product number with the conjugate transposition of their own.
3. a kind of adaptive equalizer, is characterized in, including：
Adaptive equalization device according to claim 2；And
Transversal filter is connected with the adaptive equalization device, the tap weight value of the transversal filter by it is described from
Adaptive equalization device is calculated.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610203866.1A CN105635008B (en) | 2016-04-01 | 2016-04-01 | A kind of adaptive equilibrium method, device and adaptive equalizer |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610203866.1A CN105635008B (en) | 2016-04-01 | 2016-04-01 | A kind of adaptive equilibrium method, device and adaptive equalizer |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105635008A CN105635008A (en) | 2016-06-01 |
| CN105635008B true CN105635008B (en) | 2018-11-20 |
Family
ID=56049518
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610203866.1A Active CN105635008B (en) | 2016-04-01 | 2016-04-01 | A kind of adaptive equilibrium method, device and adaptive equalizer |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105635008B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107820677B (en) * | 2017-09-18 | 2021-03-19 | 深圳市汇顶科技股份有限公司 | Method and device for determining filter coefficient and terminal |
| CN114301529B (en) * | 2021-12-31 | 2024-04-02 | 杭州电子科技大学 | Volterra equalization method and system based on multi-symbol processing |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101478350A (en) * | 2009-01-20 | 2009-07-08 | 南京信息工程大学 | Combined super-exponential iteration blind equalization algorithm based on orthogonal wavelet transform |
| CN103780519A (en) * | 2014-01-07 | 2014-05-07 | 电子科技大学 | Channel equalization and frequency deviation estimation joint parallel method based on LMS |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2004091127A2 (en) * | 2003-03-31 | 2004-10-21 | University Of Florida | Blind equalizers using probability density matching and parzen windowing |
-
2016
- 2016-04-01 CN CN201610203866.1A patent/CN105635008B/en active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101478350A (en) * | 2009-01-20 | 2009-07-08 | 南京信息工程大学 | Combined super-exponential iteration blind equalization algorithm based on orthogonal wavelet transform |
| CN103780519A (en) * | 2014-01-07 | 2014-05-07 | 电子科技大学 | Channel equalization and frequency deviation estimation joint parallel method based on LMS |
Non-Patent Citations (2)
| Title |
|---|
| 无线通信系统中自适应信道均衡算法研究与实现;罗书建;《中国优秀硕士学位论文全文数据库 信息科技辑》;20160315(第3期);第I136-384页 * |
| 自适应均衡技术的研究;周杨，黄元峰;《国外电子元器件》;20081105;第72-74页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105635008A (en) | 2016-06-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN101902258B (en) | Method and device for acquiring digital pre-distortion processing parameter | |
| CN105635008B (en) | A kind of adaptive equilibrium method, device and adaptive equalizer | |
| JP2013542696A5 (en) | ||
| CN109274308A (en) | Parameters of hydraulic motor control system and parameters of hydraulic motor control method | |
| Butler et al. | Noniterative automatic equalization | |
| Niedźwiecki et al. | System identification based approach to dynamic weighing revisited | |
| CN103957177B (en) | A kind of multistage time domain blind balance method of ofdm signal | |
| JP6347190B2 (en) | Information processing apparatus, program, and method | |
| JP2008510323A5 (en) | ||
| CN105991484A (en) | Frequency offset estimation method and frequency offset estimation device | |
| CN111698026B (en) | Polarization state change monitoring method, device and receiver | |
| Zhao et al. | New variable step-sizes minimizing mean-square deviation for the LMS-type algorithms | |
| Mitra et al. | A block floating-point treatment to the LMS algorithm: efficient realization and a roundoff error analysis | |
| Cai et al. | New explicit multisymplectic scheme for the complex modified Korteweg-de Vries equation | |
| Zhang et al. | WDOP‐based Summation Inequality and its Application to Exponential Stability of Linear Delay Difference Systems | |
| Gray Jr | Numerical methods in digital real-time simulation | |
| CN108322410B (en) | Time domain equalizer and signal processing method thereof | |
| CN103297362A (en) | Error-updating-based system and method for channel estimation | |
| CN114577425B (en) | Transfer function identification method for multi-input multi-output vibration test system | |
| CN109728794A (en) | Filter coefficient computing device and method and its signal generating apparatus and method | |
| TWI270277B (en) | Transform domain adaptive filter, equalizer and wireless communication device using the same, and decision feedback equalization method | |
| CN104459316B (en) | The subharmonic measuring method and system of power system | |
| CN107547053B (en) | Signal digital predistortion processing method and device | |
| Preis | Least-squares time-domain deconvolution for transversal-filter equalisers | |
| EP2897334B1 (en) | Adaptive filtering method and system based on error sub-band |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |