System and method for adaptive filter

Abstract

In one embodiment, a method for training an adaptive filter includes receiving, by a processor from a device, an input signal and a training reference signal and determining a correlation matrix in accordance with the input signal, the training reference signal, and a filter type. The method also includes determining a plurality of coefficients in accordance with the correlation matrix and adjusting the adaptive filter in accordance with the plurality of coefficients.

Description

TECHNICAL FIELD

The present invention relates to a system and method for filtering, and, in particular, to a system and method for adaptive filtering.

BACKGROUND

Adaptive filtering is used in a variety of situations, including power amplifier predistortion systems. Predistortion is a technique which improves the linearity of power amplifiers, for example in telecommunications systems. A power amplifier with nonlinearities causes interference on other radio channels. Predistortion circuits inversely model the power amplifier's gain and phase characteristics to produce a system that is more linear. Inverse distortion is introduced into the input of the power amplifier, cancelling nonlinearities in the amplifier. The characteristics of the adaptive filter may vary by the type of power amplifier or power amplifier sub-system architecture.

An adaptive equalizer provides feedback, for example to equalize the channel gain across frequency bandwidths to compensate for different gains at different frequencies. Also, adaptive filters may be used for interference calculations in other types of adaptive systems. An adaptive filter self-adjusts its transfer function based on an optimization algorithm from an error signal. In an example, an adaptive process uses a cost function, which is a criterion for optimum performance of the filter, as an input to an optimization algorithm. The algorithm determines how to modify the filter transfer function to minimize the cost of the next iteration.

SUMMARY

An embodiment method for training an adaptive filter includes receiving, by a processor from a device, an input signal and a training reference signal and determining a correlation matrix in accordance with the input signal, the training reference signal, and a filter type. The method also includes determining a plurality of coefficients in accordance with the correlation matrix and adjusting the adaptive filter in accordance with the plurality of coefficients.

An embodiment method for training an adaptive filter includes receiving, by a processor from a device, an input signal and a training reference signal and determining a tri-angle matrix in accordance with the input signal and the training reference signal. The method also includes storing the tri-angle matrix and determining a plurality of coefficients in accordance with the tri-angle matrix. Additionally, the method includes adjusting the adaptive filter in accordance with the plurality of coefficients.

An embodiment computer includes a processor and a computer readable storage medium storing programming for execution by the processor. The programming includes instructions to receive an input signal and a training reference signal and determine a correlation matrix in accordance with the input signal, the training reference signal, and a filter type. The programming also includes instructions to determine a plurality of coefficients in accordance with the correlation matrix and adjust an adaptive filter in accordance with the plurality of coefficients.

The foregoing has outlined rather broadly the features of an embodiment of the present invention in order that the detailed description of the invention that follows may be better understood. Additional features and advantages of embodiments of the invention will be described hereinafter, which form the subject of the claims of the invention. It should be appreciated by those skilled in the art that the conception and specific embodiments disclosed may be readily utilized as a basis for modifying or designing other structures or processes for carrying out the same purposes of the present invention. It should also be realized by those skilled in the art that such equivalent constructions do not depart from the spirit and scope of the invention as set forth in the appended claims.

BRIEF DESCRIPTION OF THE DRAWINGS

For a more complete understanding of the present invention, and the advantages thereof, reference is now made to the following descriptions taken in conjunction with the accompanying drawing, in which:

FIG. 1 illustrates an embodiment system for an adaptive filter;

FIG. 2 illustrates a flowchart for an embodiment method for adaptive filtering;

FIG. 3 illustrates another flowchart for an embodiment method of adaptive filtering; and

FIG. 4 illustrates a block diagram of an embodiment general-purpose computer system.

Corresponding numerals and symbols in the different figures generally refer to corresponding parts unless otherwise indicated. The figures are drawn to clearly illustrate the relevant aspects of the embodiments and are not necessarily drawn to scale.

DETAILED DESCRIPTION OF ILLUSTRATIVE EMBODIMENTS

It should be understood at the outset that although an illustrative implementation of one or more embodiments are provided below, the disclosed systems and/or methods may be implemented using any number of techniques, whether currently known or in existence. The disclosure should in no way be limited to the illustrative implementations, drawings, and techniques illustrated below, including the exemplary designs and implementations illustrated and described herein, but may be modified within the scope of the appended claims along with their full scope of equivalents.

An embodiment uses a flexible adaptation architecture to achieve convergence in different manners and adapts quickly for different types of training targets. A common hardware architecture may be used for various training algorithms, such as least squares (LS)-GIVEN, least squares QR decomposition (LS-QRD), QR decomposition recursive least squares (QRD-RLS) or (QRD), and least mean squares (LMS). Configurable parameters include the training algorithm, the filter type, the matrix size, the number of samples for the correlation matrix calculation, the number of rows within the correlation matrix (correlation U vectors or U vector) for basic tri-angle rotation matrix calculation, multiple level adaptations with forgetting factors λ₁, λ₂, and λ₃, coefficient update parameter μ, and regularization factor α.

FIG. 1 illustrates system 100 for adaptive filter training. Vector engine 104, rotation engine 106, and coefficient engine 110 are common to LS-GIVEN, LS-QRD, and QRD filter types. Correlation matrix engine 102 and tri-angle matrix engine 108 are programmable to facilitate LS-GIVEN, LS-QRD, QRD, and LMS architectures.

Correlation matrix engine 102 generates a U vector and a correlation matrix R_xd. In one example, the U vector is composed of the input signal x, which is the output of the pre-actuator model, when direct training is used, or the plant output when indirect training is used. The general formula for the training input U vector is a function of the input signal, and is given by:U=f(x₀,x₁, . . . ,x_i),where i represents the total number of training model parameters. In general, the k^th correlation matrix is defined by the recursive relation:R_xdk=R_xdk^bλ₁R_xdk−1.  (1)R_xdk−1 is the (k−1)^th correlation matrix, λ₁ is the forgetting factor of the matrix R_xd, and R_xdk^b is the basic correlation matrix at the k^th correlation matrix block. The correlation matrix R_xdk^b is composed of the correlation matrix U^H×U and the cross correlation vector U^H×d, where d is the training reference signal. That is:R_xdk^b=[U^H×UU^H×d]_M×(M+1).The correlation matrix of the U vector is renamed to:R_u=[U^H×U]_M×M+αI_M.Regularization factor α, which is programmable, handles the matrix ill condition. The U vector matrix is given by:

$U={\left[\begin{array}{ccccc}{u}_0^{k\times \left(N-1\right)+i}& u_1^{k\times \left(N-1\right)+i-1}& u_2^{k\times \left(N-1\right)+i-1}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-1}\\ u_0^{k\times \left(N-1\right)+i-1}& u_1^{k\times \left(N-1\right)+i-1}& u_2^{k\times \left(N-1\right)+i-1}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-1}\\ u_0^{k\times \left(N-1\right)+i-2}& u_1^{k\times \left(N-1\right)+i-2}& u_2^{k\times \left(N-1\right)+i-2}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-2}\\ u_0^{k\times \left(N-1\right)+i-3}& u_1^{k\times \left(N-1\right)+i-3}& u_2^{k\times \left(N-1\right)+i-3}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-3}\\ ⋮& ⋮& ⋮& \ddots & ⋮\\ u_0^{k\times \left(N-1\right)+i-N+1}& u_1^{k\times \left(N-1\right)+i-N+1}& u_2^{k\times \left(N-1\right)+i-N+1}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-N+1}\end{array}\right]}_{N\times M}.$ M is the number of coefficients, the superscript i is the i^th sample within the k^th correlation matrix, and N is the number of samples to be used for the correlation matrix calculation. The number of samples, N to be used for the R_xd matrix calculation may be programmable. The cross correlation vector of U and d is given by:R_d=[U^H×d]_N×1,where d is the training reference signal. In a direct learning algorithm, d is the error signal between the actuator input signal and the feedback signal. For indirect learning, d is the actuator output signal given by:

$d={\left[\begin{array}{c}{d}^{k\times \left(N-1\right)+i}\\ d^{k\times \left(N-1\right)+i-1}\\ ⋮\\ d^{k\times \left(N-1\right)+i-N+1}\end{array}\right]}_{N\times 1}.$ For QRD and LMS algorithms, the correlation matrix R_udk equals the U vector matrix and d vector, with R_udk given by:

${\left[\begin{array}{cccccc}{u}_0^{k\times \left(N-1\right)+i}& u_1^{k\times \left(N-1\right)+i}& u_2^{k\times \left(N-1\right)+i}& \dots & u_{M-1}^{k\times \left(N-1\right)+i}& d^{k\times \left(N-1\right)+i}\\ u_0^{k\times \left(N-1\right)+i-1}& u_1^{k\times \left(N-1\right)+i-1}& u_2^{k\times \left(N-1\right)+i-1}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-1}& d^{k\times \left(N-1\right)+i-1}\\ u_0^{k\times \left(N-1\right)+i-2}& u_1^{k\times \left(N-1\right)+i-2}& u_2^{k\times \left(N-1\right)+i-2}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-2}& d^{k\times \left(N-1\right)+i-2}\\ u_0^{k\times \left(N-1\right)+i-3}& u_1^{k\times \left(N-1\right)+i-3}& u_2^{k\times \left(N-1\right)+i-3}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-3}& d^{k\times \left(N-1\right)+i-3}\\ ⋮& ⋮& ⋮& \ddots & ⋮& ⋮\\ u_0^{k\times \left(N-1\right)+i-M+1}& u_1^{k\times \left(N-1\right)+i-M+1}& u_2^{k\times \left(N-1\right)+i-M+1}& \dots & u_{M-1}^{k\times \left(N-1\right)+i-M+1}& d^{k\times \left(N-1\right)+i-M+1}\end{array}\right]}_{M\times \left(M+1\right)}.$

Vector engine 104 calculates diagonal R elements for tri-angle matrix engine 108 and rotation matrix G_n for rotation engine 106. For the n^th rotation, the rotation matrix G_n is defined by:

$G_n=\left[\begin{array}{cc}{c}_n& s_n^{*}\\ -s_n& c_n\end{array}\right],$ where c_n is given by:

$c_n=\frac{r_{n-1}^i\left(n,n\right)}{\sqrt{{r_{n-1}^i\left(n,n\right)}^2+{u_{n-1}^i\left(n\right)}^2}},$ and s_n is given by:

$s_n=\frac{u_{n-1}^i\left(n\right)}{\sqrt{{r_{n-1}^i\left(n,n\right)}^2+{u_{n-1}^i\left(n\right)}^2}}.$ The diagonal elements of the R_ud matrix are then calculated by:

$r_n^i=\left[\begin{array}{cc}{c}_n& s_n^{*}\end{array}\right]\left[\begin{array}{c}{r}_{n-1}^i\left(n,n\right)\\ u_{n-1}^i\left(n\right)\end{array}\right]=\sqrt{{r_{n-1}^i\left(n,n\right)}^2+{u_{n-1}^i\left(n\right)}^2.}$

Rotation engine 106 calculates the rest of the elements in tri-angle matrix R_ud. The additional elements include r_n^k(i,j) and z_n^k(i). These additional elements may be obtained by:

$\left[\begin{array}{c}{x}_n\\ y_n\end{array}\right]=G_n\left[\begin{array}{c}{x}_{n-1}\\ y_{n-1}\end{array}\right].$

Tri-angle matrix engine 108 implements tri-angle matrix R_ud. The tri-angle matrix engine stores intermediate and final tri-angle matrix results. The general recursive j^th tri-angle matrix is given by:R_udj=R_udj^b+λ₃×R_udj−1,where R_udj^b is the basic tri-angle matrix. The basic tri-angle matrix is represented by:R_udj^b=λ₂×R_udj−1^b.The initial basic tri-angle matrix is configured based on the filter type. For LS-GIVEN, the initial tri-angle matrix equals the correlation matrix R_xdk. In LS-GIVEN's rotation, the rotation is performed on the full correlation matrix. In the LS-QRD and QRD, the basic tri-angle matrix is given by:

$R_{{\mathrm{ud}}_0}^b={\left[\begin{array}{cccccc}1& 0& 0& \dots & 0& 0\\ 0& 1& 0& \dots & 0& 0\\ 0& 0& 1& \dots & 0& 0\\ ⋮& ⋮& ⋮& \ddots & ⋮& ⋮\\ 0& 0& 0& \dots & 1& 0\end{array}\right]}_{M\times \left(M+1\right)}$ For QRD or LS-QRD, rotation is performed on a single row or a block row within the correlation matrix. After i rotations, the tri-angle matrix becomes:

$R_{{\mathrm{ud}}_j}={\left[\begin{array}{cccccc}{r}_{u_{11}}^i& r_{u_{12}}^i& r_{u_{13}}^i& \dots & r_{u_{1M}}^i& r_{d_1}^i\\ 0& r_{u_{22}}^i& r_{u_{23}}^i& \dots & r_{u_{2M}}^i& r_{d_2}^i\\ 0& 0& r_{u_{33}}^i& \dots & r_{u_{34}}^i& r_{d_3}^i\\ ⋮& ⋮& ⋮& \ddots & ⋮& ⋮\\ 0& 0& 0& \dots & r_{u_{\mathrm{MM}}}^i& r_{d_M}^i\end{array}\right]}_{M\times \left(M+1\right)}.$ For the LMS configuration, the tri-angle matrix is not used. Tri-angle matrix engine block 108 may be disabled, and the memory used for storing tri-angle matrix elements may be assigned to extend the size of the correlation matrix.

Coefficient engine 110 performs back substitution and updates the coefficients. The coefficient calculation and update are based on the matrix update configuration. This may be performed on a single U vector, a basic block, or several basic blocks based on the algorithm chosen and the training system requirements. Table 1 below illustrates various matrix update configurations.

	TABLE 1
	Matrix Update
	Correlation	Tri-angle
Algorithm	Matrix	Matrix	Description
LS-	Single	Basic Rotation	Matrix rotation and coefficient calculation are based on a basic
GIVEN's	basic block	Matrix	Correlation Matrix and basic Tri-angle Matrices
			λ1, and λ3 = 0; λ2 = 1
		Multiple basic	Coefficient calculation is based on a basic Correlation Matrix and
		Rotation	multiple basic Tri-angle Matrices
		Matrices	λ1 = 0, λ3 = (0, 1]; λ2 = 1
	Multiple	Basic Rotation	Matrix rotation and coefficient calculation are based on multiple
	basic blocks	Matrix	basic Correlation Matrices and a basic Tri-angle Matrix
			λ1 = (0, 1]; λ2 = 1; λ3 = 0
		Multiple Basic	Matrix rotation and coefficient calculation are based on multiple
		Rotation	basic Correlation Matrices and multiple basic Tri-angle Matrices
		Matrices	λ1 = (0, 1]; λ2 = 1; λ3 = (0, 1]
LS-QRD	Partial	Basic Rotation	Matrix rotation and coefficient calculation are based on several
	basic block	Matrix	rows within a basic Correlation Matrix or multiple Correlation
			Matrices and a basic Tri-angle Matrix
			λ1 = 0;, λ3 = 0; λ2 = (0, 1]
		Multiple Basic	Matrix rotation and coefficient calculation are based on several
		Rotation	rows within a basic Correlation Matrix or multiple Correlation
		Matrices	Matrices and multiple basic Tri-angle Matrices
			λ1, = 0; λ3 = (0, 1]; λ2 = (0, 1]
	Single	Basic Rotation	Matrix rotation and coefficient calculation are based on all rows of
	basic block	Matrix	a basic Correlation Matrix and a basic Tri-angle Matrix
			λ1, = 0; λ3 = 0; λ2 = (0, 1]
		Multiple Basic	Matrix rotation and coefficient calculation are based on all rows of
		Rotation	a basic Correlation Matrix and multiple basic Tri-angle Matrices
		Matrices	λ1, = 0; λ3 = (0, 1]; λ2 = (0, 1]
	Multiple	Basic Rotation	Matrix rotation and coefficient calculation are based on all rows of
	basic blocks	Matrix	multiple basic Correlation Matrices and a basic Tri-angle Matrix
			λ1, = (0, 1]; λ3 = 0; λ2 = (0, 1]
		Multiple Basic	Matrix rotation and coefficient calculation are based on all rows of
		Rotation	multiple basic Correlation Matrices and multiple basic Tri-angle
		Matrices	Matrices
			λ1, = (0,1]; λ3 = (0, 1]; A2 = (0, 1]
QRD	Single U	Basic Rotation	Matrix rotation and coefficient calculation are based on single
	vector	Matrix	rows of multiple basic Correlation Matrix and a basic Tri-angle
			Matrix
			λ1, = 0; λ3 = 0; λ2 = (0, 1]
		Multiple Basic	Matrix rotation and coefficient calculation are based on a single
		Rotation	row of a multiple basic Correlation Matrices and multiple basic Tri-
		Matrices	angle Matrices
			λ1, = 0; λ3 = (0, 1]; λ2 = (0, 1]
LMS	Single U	N/A	Coefficient calculation and update is based on a single U vector
	vector		within a Single Basic Correlation Matrix (sample by sample update)
			λ1, = 0; λ3 = N/A; λ2 = N/A
	Single	N/A	Several or all U vectors within Single Basic Correlation Matrix are
	basic block		used for coefficient calculation and update (Block based update)
			λ1, = 0; λ3 = N/A; λ2 = N/A

For LS-GIVEN, a single basic block correlation matrix or multiple basic block correlation matrices may be used. With either a single basic block correlation matrix or multiple basic block matrices, the tri-angle matrix may be a basic rotation matrix or multiple basic rotation matrices. In a basic rotation matrix, the matrix rotation and coefficient calculation are based on the basic correlation matrix and the basic tri-angle matrix. In basic rotation, λ₁=0, λ₂=1, and λ₃=0. For multiple basic rotation matrices, calculations are based on the basic correlation matrix and multiple basic tri-angle matrices, where λ₁=0, λ₂=1, and λ₃=(0, 1]. When multiple basic blocks are used for the correlation matrix, a basic rotation matrix or multiple basic rotation matrix may be used for the tri-angle matrix. For a basic rotation matrix as the tri-angle matrix, matrix rotation and coefficient calculation are based on multiple basic correlation matrices, and the basic tri-angle matrix, where λ₁=(0, 1], λ₂=1, and λ₃=0. For multiple basic rotation matrix, the matrix rotation and coefficient calculation are based on multiple basic correlation matrices, and multiple basic tri-angle matrices, and λ₁=(0, 1], λ₂=1, and λ₃=(0, 1].

In an example with a LS-QRD filter, the correlation matrix may be a partial basic block matrix, a single basic block matrix, or multiple basic block matrices. The tri-angle matrix may be a basic rotation matrix or multiple basic rotation matrices. There may be any combination of correlation matrix and tri-angle matrix. When a partial basic block correlation matrix is used with a basic rotation tri-angle matrix, matrix rotation and coefficient calculation are based on several rows within a basic correlation matrix or multiple correlation matrices and the basic tri-angle matrix, where λ₁=0, λ₂=(0, 1], and λ₃=0. On the other hand, when a partial basic block correlation matrix is used with multiple basic rotation tri-angle matrices, matrix rotation and coefficient calculation are based on several rows within a basic correlation matrix or multiple correlation matrices and multiple basic tri-angle matrices, where λ₁=0, λ₂=(0, 1], and λ₃=(0, 1]. Also, when a single basic block correlation matrix is used with a basic rotation tri-angle matrix, the matrix rotation and coefficient calculation are based on all rows of the basic correlation matrix and the basic tri-angle matrix, where λ₁=0, λ₂=(0, 1], and λ₃=0. When a single basic block correlation matrix is used with multiple basic rotation matrices, the matrix rotation and coefficient calculation are based on all rows of a basic correlation matrix and multiple basic rotation matrices, where λ₁=0, λ₂=(0, 1], and λ₃=(0, 1]. Additionally, when multiple basic block correlation matrices and a basic rotation tri-angle matrix are used, the matrix rotation and coefficient calculations are based on all rows of multiple basic correlation matrices and a basic tri-angle matrix, where λ₁=(0, 1], λ₂=(0, 1], and λ₃=0. When multiple basic block correlation matrices are used with multiple basic rotation tri-angle matrices, matrix rotation and coefficient calculations are based on all rows of multiple basic correlation matrices and multiple basic tri-angle matrices, with λ₁=(0, 1], λ₂=(0, 1], and λ₃=(0, 1].

In another example, a QRD filter is used, and the correlation matrix is a single U vector matrix. The tri-angle matrix may be a basic rotation matrix or multiple basic rotation matrices. When a basic rotation matrix is use as the tri-angle matrix, the matrix rotation and coefficient calculations are based on single rows of multiple basic correlation matrices and a basic tri-angle matrix, where λ₁=0, λ₂=(0, 1], and λ₃=0. On the other hand, when the multiple basic rotation matrices are used, the matrix rotation and coefficient calculations are based on a single row of multiple basic correlation matrix and multiple basic tri-angle matrices, where λ₁=0, λ₂=(0, 1], and λ₃=(0, 1].

With an LMS filter, the correlation matrix may be a single U vector matrix or a single basic block matrix. There is no tri-angle matrix for LMS. When the correlation matrix is a single U vector, the coefficient calculation and update is based on a single U vector with a single basic correlation matrix for a sample by sample update, where λ₁=0, and λ₂ and λ₃ are not applicable. When a single basic block correlation matrix is used, some or all of the U vectors within a single basic correlation matrix are used for the calculation and update, with block based updating, where λ₁=0, and λ₂ and λ₃ are not applicable.

FIG. 2 illustrates the relationship of the matrix and the coefficient update. The number of rows in the correlation matrix used for matrix rotation depends on the training algorithm. For LS-GIVEN, all rows within the correlation matrix are used for matrix rotation, while for LS-QRD or QRD, a single row or a block of rows may be used for matrix rotation, depending on the configuration. On the other hand, for LMS, the tri-angle matrix is not used, and coefficient update may be based on single or block U vectors. A variety of different step sizes may be used. Smaller steps are slower, but perform better than bigger steps. Multiple step sizes and multiple matrix sizes may be used adaptively.

In block 122, the correlation matrix is generated. λ₁ is forgetting factor used for the correlation matrix calculation as defined in Equation (1). When λ₁ is zero, the current correlation matrix R_xdi is calculated without basing itself on previous correlation matrix information. At a given time i, the correlation matrix, R_xdi, may be calculated independently or based on partial correlation matrix information from the previous iteration weighted by the forgetting factor.

Then, in block 124, λ₂ is the forgetting factor for individual tri-angle matrix, and λ₃ is the forgetting factor between tri-angle matrices. At a given time, the tri-angle matrix, R_xdi may be calculated independently or based on partial information of the previous triangle matrix R_xdi−1 weighted by the corresponding forgetting factor value.

Finally, in block 126, the coefficients are updated. In block 140, the coefficient is calculated from the tri-angle matrix. Note that μ is the coefficient forgetting factor. At a given time i (or sample), the coefficients may be calculated independently or based on weighted coefficients of the previous iteration. Depending on the configuration, the coefficient may not be applied to the actuator in the current iteration i.e. multiple iterations can happen before a coefficient update is applied.

FIG. 3 illustrates flowchart 150 for an adaptive filtering method. Initially, in step 152, U vectors are generated. In one example, direct training is used, and the U vector is composed of the input signal x from the actuator output. Alternatively, when indirect training is used, the U vector is composed of the output of the target plant. The general formula for the training input U vector is given by:U=f(x₀,x₁, . . . ,x_i),where i is an integer indicating the number of training model parameters.

Next, in step 154, the correlation matrix is generated. In general, the k^th correlation matrix is defined by the recursive relation:R_xdk=R_xdk^b+λ₁R_xdk−1.R_xdk−1 is the (k−1)^th correlation matrix, λ₁ is the forgetting factor of the matrix R_xd, and R_xdk^b is the basic correlation matrix at the k^th correlation matrix block. The correlation matrix R_xdk^b is composed of the correlation matrix U^H×U and the cross correlation vector U^H×d.

Then, in step 156, rotation matrix G_n is calculated. For the n^th rotation, the rotation matrix G_n is defined by:

$G_n=\left[\begin{array}{cc}{c}_n& s_n^{*}\\ -s_n& c_n\end{array}\right],$ where c_n is given by:

$c_n=\frac{r_{n-1}^i\left(n,n\right)}{\sqrt{{r_{n-1}^i\left(n,n\right)}^2+{u_{n-1}^i\left(n\right)}^2}},$ and s_n is given by:

$s_n=\frac{u_{n-1}^i\left(n\right)}{\sqrt{{r_{n-1}^i\left(n,n\right)}^2+{u_{n-1}^i\left(n\right)}^2}}.$

In step 158, the diagonal elements of the tri-angle matrix are calculated. For instance, the diagonal elements of the tri-angle matrix may be calculated according to the formula:

$r_n^i=\left[\begin{array}{cc}{c}_n& s_n^{*}\end{array}\right]\left[\begin{array}{c}{r}_{n-1}^i\left(n,n\right)\\ u_{n-1}^i\left(n\right)\end{array}\right]=\sqrt{{r_{n-1}^i\left(n,n\right)}^2+{u_{n-1}^i\left(n\right)}^2.}$

The remaining elements of the tri-angle matrix are calculated in step 160. The additional elements include r_n^k(i,j) and z_n^k(i). These additional elements may be calculated, for example by:

$\left[\begin{array}{c}{x}_n\\ y_n\end{array}\right]=G_n\left[\begin{array}{c}{x}_{n-1}\\ y_{n-1}\end{array}\right].$

Next, in step 164, the tri-angle matrix is implemented. The tri-angle matrix engine stores the intermediate and final tri-angle matrix results. The initial basic tri-angle matrix is configured based on the filter type. For LS-GIVEN, the initial tri-angle matrix equals the correlation matrix R_xdk. On the other hand, LS-GIVEN's rotation is performed on the full correlation matrix. In another example with QRD or LS-QRD, rotation is performed on a single row or a block row within the correlation matrix. For the LMS configuration, the tri-angle matrix is not used. Tri-angle matrix engine block 108 may be disabled, and the memory used for storing tri-angle matrix elements may be assigned to extend the size of the correlation matrix.

Finally, in step 166, back substitution and coefficient updates are performed. The coefficient calculation and update are based on the matrix update configuration. This may be performed on a single U vector, a basic block, or several basic blocks based on the algorithm chosen and the training system requirements.

Embodiments may have a variety of physical implementations. For example, a field programmable gate array (FPGA), application specific integrated circuit (ASIC), digital signal processor, or general purpose computer may be used. These are all considered as examples of the general category of processor.

FIG. 4 illustrates a block diagram of processing system 270 that may be used for implementing the devices and methods disclosed herein. Specific devices may utilize all of the components shown, or only a subset of the components, and levels of integration may vary from device to device. Furthermore, a device may contain multiple instances of a component, such as multiple processing units, processors, memories, transmitters, receivers, etc. The processing system may comprise a processing unit equipped with one or more input devices, such as a microphone, mouse, touchscreen, keypad, keyboard, and the like. Also, processing system 270 may be equipped with one or more output devices, such as a speaker, a printer, a display, and the like. The processing unit may include central processing unit (CPU) 274, memory 276, mass storage device 278, video adapter 280, and I/O interface 288 connected to a bus.

The bus may be one or more of any type of several bus architectures including a memory bus or memory controller, a peripheral bus, video bus, or the like. CPU 274 may comprise any type of electronic data processor. Memory 276 may comprise any type of system memory such as static random access memory (SRAM), dynamic random access memory (DRAM), synchronous DRAM (SDRAM), read-only memory (ROM), a combination thereof, or the like. In an embodiment, the memory may include ROM for use at boot-up, and DRAM for program and data storage for use while executing programs.

Mass storage device 278 may comprise any type of storage device configured to store data, programs, and other information and to make the data, programs, and other information accessible via the bus. Mass storage device 278 may comprise, for example, one or more of a solid state drive, hard disk drive, a magnetic disk drive, an optical disk drive, or the like.

Video adaptor 280 and I/O interface 288 provide interfaces to couple external input and output devices to the processing unit. As illustrated, examples of input and output devices include the display coupled to the video adapter and the mouse/keyboard/printer coupled to the I/O interface. Other devices may be coupled to the processing unit, and additional or fewer interface cards may be utilized. For example, a serial interface card (not pictured) may be used to provide a serial interface for a printer.

The processing unit also includes one or more network interface 284, which may comprise wired links, such as an Ethernet cable or the like, and/or wireless links to access nodes or different networks. Network interface 284 allows the processing unit to communicate with remote units via the networks. For example, the network interface may provide wireless communication via one or more transmitters/transmit antennas and one or more receivers/receive antennas. In an embodiment, the processing unit is coupled to a local-area network or a wide-area network for data processing and communications with remote devices, such as other processing units, the Internet, remote storage facilities, or the like.

While several embodiments have been provided in the present disclosure, it should be understood that the disclosed systems and methods might be embodied in many other specific forms without departing from the spirit or scope of the present disclosure. The present examples are to be considered as illustrative and not restrictive, and the intention is not to be limited to the details given herein. For example, the various elements or components may be combined or integrated in another system or certain features may be omitted, or not implemented.

In addition, techniques, systems, subsystems, and methods described and illustrated in the various embodiments as discrete or separate may be combined or integrated with other systems, modules, techniques, or methods without departing from the scope of the present disclosure. Other items shown or discussed as coupled or directly coupled or communicating with each other may be indirectly coupled or communicating through some interface, device, or intermediate component whether electrically, mechanically, or otherwise. Other examples of changes, substitutions, and alterations are ascertainable by one skilled in the art and could be made without departing from the spirit and scope disclosed herein.

Claims

1. A method for training an adaptive filter, the method comprising: selecting a filter type from the group consisting of least squares (LS)-GIVEN, least squares QR decomposition (LS-QRD), QR decomposition recursive least squares (QRD), and least mean squares (LMS);receiving, by a processor from a device, an input signal x and a training reference signal d;determining, by the processor, a correlation matrix Rxd in accordance with the input signal x, the training reference signal d, and the selected filter type;determining, by the processor, a plurality of coefficients for the adaptive filter in accordance with the correlation matrix Rxd;adjusting, by the processor, the adaptive filter in accordance with the plurality of coefficients; andpredistorting an input of a power amplifier utilizing the adjusted adaptive filter.

2. The method of claim 1, further comprising determining a U vector, wherein determining the correlation matrix Rxd further comprises determining the correlation matrix Rxd in accordance with the U vector.

3. The method of claim 2, wherein determining the U vector comprises determining the U vector in accordance with the input signal x and the filter type.

4. The method of claim 2, wherein determining the U vector comprises determining the U vector in accordance with the training reference signal d and the filter type.

5. The method of claim 2, further comprising determining a size of the U vector.

6. The method of claim 1, further comprising determining a size of the correlation matrix Rxd.

7. The method of claim 1, further comprising: determining a first forgetting factor λ1 in accordance with the filter type;determining a second forgetting factor λ2 in accordance with the filter type; anddetermining a third forgetting factor λ3 in accordance with the filter type, wherein determining the correlation matrix Rxd further comprises determining the correlation matrix Rxd in accordance with the first forgetting factor λ1, the second forgetting factor λ2, and the third forgetting factor λ3.

8. The method of claim 1, further comprising determining a regularization factor α, wherein determining the correlation matrix Rxd further comprises determining the correlation matrix Rxd in accordance with the regularization factor α.

9. The method of claim 1, wherein the correlation matrix Rxd consists of a basic block correlation matrix.

10. The method of claim 1, wherein the correlation matrix Rxd comprises: a first basic block correlation matrix; anda second basic block correlation matrix.

11. The method of claim 1, wherein the correlation matrix Rxd comprises a partial basic block correlation matrix.

12. The method of claim 1, wherein the correlation matrix Rxd comprises a U vector.

13. A method for training an adaptive filter, the method comprising: selecting a filter type from the group consisting of least squares (LS)-GIVEN, least squares QR decomposition (LS-QRD), and QR decomposition recursive least squares (QRD);receiving, by a processor from a device, an input signal x and a training reference signal d;determining, by the processor, a tri-angle matrix Rud in accordance with the input signal x, the training reference signal d, and the selected filter type;storing, by the processor, the tri-angle matrix Rud;determining, by the processor, a plurality of coefficients for the adaptive filter in accordance with the tri-angle matrix Rud;adjusting, by the processor, the adaptive filter in accordance with the plurality of coefficients; andpredistorting an input of a power amplifier utilizing the adjusted adaptive filter.

14. The method of claim 13, wherein determining the tri-angle matrix Rud comprises: determining a first plurality of diagonal elements of the tri-angle matrix Rud; anddetermining a second plurality of non-diagonal elements of the tri-angle matrix Rud.

15. The method of claim 14, wherein determining the second plurality of non-diagonal elements comprises determining a rotation matrix.

16. The method of claim 13, further comprising determining a size of the tri-angle matrix Rud.

17. The method of claim 13, further comprising determining a coefficient forgetting factor μ, wherein determining the plurality of coefficients for the adaptive filter comprises determining the plurality of coefficients in accordance with the coefficient forgetting factor μ.

18. The method of claim 13, wherein the tri-angle matrix Rud consists of a basic rotation matrix.

19. The method of claim 13, wherein the tri-angle matrix Rud comprises: a first basic tri-angle matrix Rudjb; anda previous tri-angle matrix Rudj−1.

20. A power amplifier predistortion system comprising: a processor; anda computer readable storage medium storing programming for execution by the processor, the programming including instructions to: selecting a filter type from the group consisting of least squares (LS)-GIVEN, least squares QR decomposition (LS-QRD), QR decomposition recursive least squares (QRD), and least mean squares (LMS);receive an input signal x and a training reference signal d,determine a correlation matrix Rxd in accordance with the input signal x, the training reference signal d, and the selected filter type,determine a plurality of coefficients for an adaptive filter in accordance with the correlation matrix Rxd,adjust the adaptive filter in accordance with the plurality of coefficients; andpredistort an input of a power amplifier utilizing the adjusted adaptive filter.