Least Squares (LS) is a deterministic parameter extraction algorithm that solves for the coefficient vector w by minimizing the cost function J(w) = ||d - Xw||², where d is the desired signal vector and X is the observation matrix. The closed-form solution, given by the Normal Equation w = (XᴴX)⁻¹Xᴴd, provides the optimal Wiener filter coefficients in a single batch computation without requiring iterative gradient descent steps.
Glossary
Least Squares (LS)

What is Least Squares (LS)?
Least Squares (LS) is a foundational batch estimation method that determines the optimal coefficient vector for a digital predistorter by minimizing the sum of squared errors between the observed output and the desired linear signal.
In digital predistortion, LS estimation is typically performed during offline training phases using captured input-output data from the power amplifier. While computationally intensive due to the required matrix inversion, LS provides the minimum Mean Squared Error (MSE) solution for stationary systems. Its numerical stability depends on the condition number of the autocorrelation matrix XᴴX, often necessitating regularization parameters or decomposition methods like QR decomposition (QRD) to mitigate ill-conditioning in wideband or highly correlated signal environments.
Key Characteristics of Least Squares
The Least Squares (LS) method provides a closed-form, batch solution for identifying optimal predistorter coefficients by minimizing the sum of squared errors. It serves as the theoretical bedrock for more complex adaptive algorithms.
The Normal Equation
The LS solution is derived analytically by setting the gradient of the cost function to zero, resulting in the Normal Equation: w = (X^H X)^-1 X^H y. This provides the optimal coefficient vector w in a single computation, assuming the input correlation matrix X^H X is invertible.
Batch Processing Paradigm
LS operates on a complete block of captured data. This offline training approach is ideal for initial parameter extraction and system characterization where real-time adaptation is not required. It processes the entire dataset at once to compute the global minimum of the error surface.
Computational Complexity
The primary cost is the matrix inversion of (X^H X), a square matrix of size M x M (where M is the number of model coefficients). This has a complexity of O(M³). For large models, this cubic growth becomes a bottleneck, motivating the use of recursive methods like Recursive Least Squares (RLS).
Numerical Stability & Regularization
Ill-conditioned input signals can make the matrix (X^H X) nearly singular, leading to unstable solutions. Tikhonov regularization is applied by adding a small penalty term: w = (X^H X + λI)^-1 X^H y. The regularization parameter λ improves the condition number and prevents overfitting.
Gauss-Markov Assumptions
Under the assumption of zero-mean, uncorrelated noise with constant variance, the LS estimator is the Best Linear Unbiased Estimator (BLUE). This means it achieves the minimum possible variance among all linear unbiased estimators, a critical property for precise behavioral modeling.
Least Squares vs. Adaptive Estimation Methods
Comparison of batch least squares against recursive and iterative adaptive algorithms for digital predistortion coefficient extraction.
| Feature | Least Squares (LS) | Recursive Least Squares (RLS) | Least Mean Squares (LMS) |
|---|---|---|---|
Estimation Mode | Batch (offline) | Recursive (online) | Iterative (online) |
Computational Complexity | O(N³) for matrix inversion | O(N²) per iteration | O(N) per iteration |
Convergence Rate | Instantaneous (one-shot) | Fast (typically < 50 samples) | Slow (signal-dependent) |
Numerical Stability | Moderate; requires regularization for ill-conditioned matrices | High with QR decomposition variant | High; inherently stable |
Memory Requirement | Stores full data matrix | Stores N×N correlation matrix | Stores coefficient vector only |
Tracks Time-Varying Systems | |||
Steady-State Misadjustment | 0% (optimal for stationary data) | 0.1-0.5% | 1-5% |
Suitable for FPGA Implementation |
Enabling Efficiency, Speed & Accuracy
Intelligent Analysis, Decision & Execution
We build AI systems for teams that need search across company data, workflow automation across tools, or AI features inside products and internal software.
Talk to Us
Search across company data
Give teams answers from docs, tickets, runbooks, and product data with sources and permissions.
Useful when people spend too long searching or get different answers from different systems.

Automate internal workflows
Use AI to route work, draft outputs, trigger actions, and keep approvals and logs in place.
Useful when repetitive work moves across multiple tools and teams.

Add AI to products and internal tools
Build assistants, guided actions, or decision support into the software your team or customers already use.
Useful when AI needs to be part of the product, not a separate tool.
Frequently Asked Questions
Clear, technically precise answers to the most common questions about the Least Squares method for digital predistortion coefficient estimation.
The Least Squares (LS) method is a batch estimation technique that finds the optimal coefficient vector by minimizing the sum of squared errors between the observed power amplifier output and the desired linear signal. It operates by constructing an overdetermined system of linear equations from a block of captured input-output data, then solving the Normal Equation to obtain the coefficient vector that minimizes the Euclidean norm of the error vector. The solution is given by w = (X^H X)^{-1} X^H y, where X is the regression matrix formed from basis functions of the input signal, y is the vector of observed outputs, and w is the predistorter coefficient vector. This closed-form solution provides the optimal estimate in a single computation step, making it ideal for offline training and initial coefficient extraction.
Related Terms
Explore the fundamental algorithms and mathematical concepts that underpin Least Squares estimation for digital predistortion coefficient extraction.
Normal Equation
The closed-form analytical solution to the linear least squares problem. It is derived by setting the gradient of the cost function to zero, yielding the optimal coefficient vector w = (X^H X)^-1 X^H y. While conceptually elegant, direct computation via matrix inversion is computationally expensive (O(n^3)) and numerically unstable for ill-conditioned matrices, making it unsuitable for real-time DPD but essential for theoretical understanding.
QR Decomposition (QRD)
A numerically robust method for solving least squares problems without explicitly computing the matrix inverse. It factors the data matrix X into an orthogonal matrix Q and an upper triangular matrix R (X = QR). The solution is then found via back-substitution on Rw = Q^H y. This approach, often implemented using Givens rotations, dramatically improves numerical stability and is the foundation for high-performance RLS implementations like QR-RLS.
Recursive Least Squares (RLS)
An adaptive algorithm that recursively computes the least squares solution, converging significantly faster than gradient-based methods like LMS. It updates the inverse of the input correlation matrix using the matrix inversion lemma, achieving a convergence rate that is an order of magnitude faster. The trade-off is a high computational cost of O(n^2) per iteration, which drives the need for efficient implementations like QR-RLS for wideband DPD applications.
Least Mean Squares (LMS)
A foundational stochastic gradient descent algorithm that offers extreme computational simplicity at O(n) per iteration. It updates coefficients by moving in the negative direction of the instantaneous gradient of the squared error. While its convergence is slow and highly dependent on the eigenvalue spread of the input signal, its low complexity makes it a viable option for resource-constrained FPGA implementations where tracking speed is not the primary bottleneck.
Regularization Parameter
A critical scalar value (often denoted by λ or δ) added to the diagonal of the correlation matrix before inversion. This technique, known as Tikhonov regularization or ridge regression, serves two vital purposes in DPD:
- Prevents numerical instability when the data matrix is ill-conditioned or rank-deficient.
- Mitigates overfitting to measurement noise, ensuring the extracted model generalizes well to new signal statistics.
Forgetting Factor
A scalar parameter (0 << λ ≤ 1) used in recursive algorithms like RLS to exponentially weight recent data more heavily than past data. A value closer to 1 provides better steady-state accuracy for stationary systems, while a smaller value enables faster tracking of time-varying amplifier characteristics caused by thermal drift or bias changes. This mechanism is essential for maintaining linearity in dynamic operating environments.

About the author
Prasad Kumkar
CEO & MD, Inference Systems
Prasad Kumkar is the CEO & MD of Inference Systems and writes about AI systems architecture, LLM infrastructure, model serving, evaluation, and production deployment. Over 5+ years, he has worked across computer vision models, L5 autonomous vehicle systems, and LLM research, with a focus on taking complex AI ideas into real-world engineering systems.
His work and writing cover AI systems, large language models, AI agents, multimodal systems, autonomous systems, inference optimization, RAG, evaluation, and production AI engineering.
Partnered with leading AI, data, and software stack.
How We Work
Custom AI workflows for your Business
One-fit-all AI don't work for modern businesses. At Inferensys, we aim to understand your business & custom requirements; which we use to define most efficient agentic workflows, the data, and the tools for your business.
01
Review the use case
We understand the task, the users, and where AI can actually help.
Read more02
Pick the right approach
We define what needs search, automation, or product integration.
Read more03
Build the first useful version
We implement the part that proves the value first.
Read more04
Improve from there
We add the checks and visibility needed to keep it useful.
Read moreThe first call is a practical review of your use case and the right next step.
Talk to Us