Stack Exchange Network

Stack Exchange network consists of 175 Q&A communities including Stack Overflow, the largest, most trusted online community for developers to learn, share their knowledge, and build their careers.

Visit Stack Exchange

Questions tagged [least-squares]

The tag has no usage guidance.

22
votes
3answers
3k views

FIR Filter design: Window vs Parks-McClellan and Least-Squares

Are there any advantages to use a window approach over Parks-McClellan (further abbreviated here as PMcC) or Least Squares algorithms for FIR filter design of a low pass filter? Assume with today's ...
15
votes
1answer
1k views

What Does Make an Error Surface Convex? Is It Determined by the Covarinace Matrix or the Hessian?

I am currently learning about least-squares (and other) estimations for regression, and from what I am also reading in some adaptive algorithm literatures, often times the phrase "... and since the ...
14
votes
1answer
700 views

How do I improve LMS / NLMS filter performance?

Are there ways to increase computational performance of a normalized least squares (NLMS) filter? Multidelay block frequency-domain (MDF) filters have been proposed to do this, but they also take away ...
8
votes
2answers
1k views

Is the Kalman Filter a Best Linear Unbiased Estimator (BLUE) for Heteroscedastic Noise?

According to the Gauss-Markov Theorem, a ordinary least squares estimator is BLUE if the noise entering a system is uncorrelated with zero mean and is homoscedastic (has a constant finite variance). I ...
6
votes
2answers
704 views

Zero Phase Filter: Determining Initial Conditions for Forward Backward Filtering

Is anybody familiar with Gustafson's algorithm for minimizing transients in forward backward filtering [1]? I'm trying to implement it and my first guess was to check Matlab's filtfilt.m, since they ...
5
votes
1answer
3k views

What's the Difference Between LMS and Gradient Descent Adaptation?

I found algorithms that seems the same to me, but they are described with different names (in field of adaptive filtering). For example: LMS - least-mean-squares seems to be GD - stochastic gradient ...
5
votes
1answer
245 views

Direct Correlation (DC) time delay estimation: variance keeps decreasing for increasing SNR?

I'm trying to reproduce the results from this paper "Discrete time techniques for Time Delay Estimation" doi:10.1109/78.193195, for both the correlation and Least Squares. I've generated a random (...
4
votes
2answers
3k views

Estimate the transfer function of an unknown system?

Suppose you have a system, H, that you want to estimate its transfer function. You have a finite number of complex input samples, x, and noisy complex (magnitude and phase) output samples, y: In ...
4
votes
6answers
336 views

The Least Norm Solution of Under Determined Linear System

Suppose I have a matrix $$A= \begin{pmatrix} 1 & 0 & 1 & 0\\ 0 & 1 & 1& 0\\ \end{pmatrix} $$ where the variables are channel information like assume $X_1$, $X_2$, $X_3$ ...
4
votes
3answers
6k views

How to remove the periodic oscillations from a signal

The task that I have is to remove the annual and semiannual oscillation from a set of temperature measurements, taken over several years, by means of least squares method. I found the method ...
4
votes
1answer
1k views

Constant Modulus Algorithm and the gradient operation

CMA is a blind channel equalization algorithm with the details presented above. I am facing 3 difficulties and shall appreciate help Q1: Does $H$ and the bar over $\bar{y_k}$ represent the Transpose ...
3
votes
3answers
241 views

LMS filter weight update

I have a general question regarding Least mean squares adaptive filters. Using the example of noise cancellation, I understand that if you have a set of reference signals (S) and corrupted signals (S+...
3
votes
2answers
10k views

Difference between Equiripple & Least Squares design for FIR digital filters

For an efficient and optimized digital FIR filter design, there are two methods available broadly, Equiripple filter design & Least Squares filter design. A general method for designing a filter ...
3
votes
2answers
254 views

How do you combine 2 least squares estimates?

Say you want to compute the least squares estimate of $w$ from a data-set: $$ \begin{bmatrix}d_1 \\d_2 \\\vdots\\d_N \end{bmatrix} =\begin{bmatrix} x_1 \\x_2 \\ \vdots \\x_N\end{bmatrix}w + \begin{...
3
votes
3answers
372 views

Force Linear Phase for a FIR filter synthesized using Berchin's FDLS?

As a follow-on to this post, how would you force linear phase for a FIR filter you synthesized using Berchin's FDLS?
3
votes
1answer
2k views

Use MATLAB to Restore a Signal from a Given Degraded Signal Using Tikhonov Regularization

Anyone could share how to develop an application in MATLAB to restore the signal from a given degraded signal using Tikhonov regularization i.e restoring the signal $f$ via solving $$ \min || g - f ...
3
votes
1answer
196 views

Least Angle Regression (LARS) without Matrix Inversion

Sorry if this is too damned long. I did what I could to abbreviate it. The question is about Least Angle Regression (LARS). I'm new to numerical work with matrices. I believe I have a way to ...
3
votes
1answer
1k views

Python: Least Squares Support Vector Machine (LS-SVM)

I'm looking for a Python package for a LS-SVM or a way to tune a normal SVM from scikit-learn to a Least-Squares Support Vector Machine for a classification problem. The goal of a SVM is to maximize ...
3
votes
0answers
193 views

Least Squares Approximation for FIR Filter Design

Following this paper , I am trying to make a least-squares algorithm in MATLAB, but for type I (I know about firls()). ...
2
votes
4answers
268 views

Sequential Form of the Least Squares for Linear Least Squares Model

I'm currently working on a project in which I need to find the tilt of a surface. Let's assume I'm only concerned with a single dimension tilt (i.e. slope) to begin. I currently have the ability to ...
2
votes
2answers
158 views

Least Square Error Estimation doubt when can we write $(A^TA)^{ -1} = A^{-1}(A^T)^{-1}$?

I am new to linear algebra and have this simple question... in least sqaure estimation...the best estimation of the equation $Ax = b$ is $x_{Estimated} = A (A^t A )^{-1} A^t b$...the projection of $b$ ...
2
votes
2answers
235 views

Doubt on Weighted Least Square Estimation

This is a page from the book linear algebra,geodesy and gps by Gilbert Strang.... the page explains about the justification of the inverse of the of the co variance matrix of measurement vector $b$ in ...
2
votes
1answer
100 views

What is the technique to find variance of estimation error

Given an $n$-vector $y$ (responses) and a design matrix $X$, I wish to fit them with a simple linear regression model $$y=X\beta+e,$$ or, $y_t = x_t'\beta_0 + e_t$ where $e\sim\mathcal{N}(0, \sigma^...
2
votes
1answer
286 views

How to Use the DFT (FFT) to Solve a Least Squares Regularization Problem (Inverse Problem)?

Let $X$ and $K$ be an image and a Point Spread Function (PSF), respectively. The blur image $B$ is obtained as follows $$B = X * K$$ I want to solve the following general regularization problem $$...
2
votes
2answers
132 views

Finding the Best Gaussian Smoothing Kernel to Minimize the Discrepancy Between Two Images

Suppose we have two grayscale images, $A$ and $B$. $A$ and $B$ very strongly resemble each other, such that the mean of the absolute difference $\lvert A - B\rvert$ is fairly small. Suppose further ...
2
votes
1answer
444 views

What Is the Difference between RLS, LMS and Wiener Filter? When Is One Preferred Over Another?

I'm dealing with a channel equalization problem where the channel is modeled as a WSS process. I understand LMS utilises a Wiener-like approach, ie it converges to the optimal (wiener) solution. I ...
2
votes
0answers
285 views

Motion artifact cancellation from ECG signal using 3-axis accelerometer data

I have a ECG signal measured as a test subject is performing 3 different activities - sitting, walking(2mph) and jogging(5mph) and also the corresponding 3-axis accelerometer(X,Y and Z) signals. I can ...
1
vote
3answers
130 views

Complex Least Squares Approximation

In the case of frequency domain FIR filter design, the error function given by : $$E(\omega)=H(e^{j\omega})-D(e^{j\omega}) \tag{1}$$ is a linear function with respect to the unknown filter ...
1
vote
2answers
98 views

Matlab : Stuck in implementing Least Squares estimator for this model

The formula to estimate $\mathbf{h}$ is then $$\hat{\mathbf{h}} = (X^T X)^{-1} X^T \vec{y}\tag{2}$$ I think this can be implemented in Matlab using ...
1
vote
2answers
56 views

What is the definition of Linear Predictive Coefficients when the optimal value is not unique?

Linear predictive coefficients of signal y are defined as the best $k$ coefficients $a_i, i = 1, \ldots, k$, that will approximate $y_n$ by $-\sum_{i=1}^k{a_iy_{n-i}}$. (Best approximation is that ...
1
vote
1answer
81 views

mpiir_l2(): Results are not IIR but FIR filters - why and how to avoid this?

I'm using mpiir_l2() from user Matt L's PhD thesis to design IIR filters. I set the number of numerator and denominator coefficients both to the same value (between ...
1
vote
1answer
216 views

Least Squares Linear Phase FIR Filter Design

In explication ''the geometric interpretation of least squares'' Typically, the number of frequency constraints is much greater than the number of design variables (filter coefficients). In these ...
1
vote
2answers
105 views

Solving a Linear Mean-Square Estimation the easy way

I have an exercise which is quite trivial. However I got stuck and I'm not sure if this the end-result. I assume there has to be a way to get this result much quicker. Given are two randomly ...
1
vote
1answer
99 views

Obtaining Correct (Least Squares Sense) Affine Transform Parameters Between Two Images

I have two images that I want to compute the affine motion model parameters. The model that I use is $$x' = a_1x+a_2y+a_3$$ $$y' = a_4x+a_5y+a_6$$ To calculate those 6 parameters, I picked 6 points (...
1
vote
2answers
271 views

Least Squares Fitting to Inverse Exponential Function

I have a time series of measurements that resembles the shape of an exponential function. The samples are a bit noisy and sometimes there is a weak sine like ripple signal ontop of it. Simplified the ...
1
vote
1answer
146 views

RLS Algorithm (Memoryless)

I have been studying the adaptive filters lately and now when I am at RLS (Recursive Least Squar) Algorithm I came across the term used in the weighting function of the RLS called forgetting factor ($\...
1
vote
2answers
285 views

Conceptual Question on Least Squares Estimation Method

What is meant by Least square plus Gaussian method of estimating the parameters of an unknown system? I am aware of least squares in general but some works refer it as least square+gaussian method. ...
1
vote
1answer
56 views

Adaptive LMS Algorithim MATLAB

I'm having some trouble implementing my LMS Adaptive Filter in MATLAB to seperate wideband and narrowband signals from a voice signal. I'm using a delayed version of my input as a reference as well ...
1
vote
2answers
584 views

Least squares and auto and cross correlation

I am trying to understand why auto and cross correlation helps find the best fit line in least squares. I have an equation as stated here: $Ax=b$ -- I have not exact solution, so I use the least ...
1
vote
4answers
223 views

Best signal processing method for this application?

I have a signal I'm getting from a sensor (shown below in the photo of an excel plot) and I need to process it to get it to look something like the red line I've overdrawn on it. Sort of a moving ...
1
vote
1answer
71 views

Equivalence of ML and DFT Peak Finding for Single Tone Estimation

My understanding is Maximum Likelihood and DFT Peak Finding for a single tone produce the same results assuming the ML is restricted to the same frequencies as the DFT. I was wondering if there was ...
1
vote
1answer
96 views

Sequential Non Linear Least Squares Problem

I have the the following non-linear function, $$s(x;A_k,\mu_k,\sigma_k)=\sum_{k=1}^2 A_k \exp\left(\frac{-(x-\mu_k)^2}{\sigma_k^2}\right)$$ with unknown (but deterministic) parameters $A_k,\mu_k,\...
1
vote
0answers
64 views

Solve for transfer function coefficients embedded in a non-linear system

Given a complex input signal $x$ and real input signal $v$, a 4th order (for simplicity) transfer function $H(z)$ is first applied to $v$ to obtain $w$, which in the time domain is represented by the ...
1
vote
0answers
35 views

Digital Filter Design Accuracy [duplicate]

Having only dealt with digital filter scarcely, the question dawned to me when I used the firls function in matlab to design an equalizer with a certain gain response. In general, can we prescribe ...
1
vote
0answers
54 views

MIL given for RLS equations vs the Woodbury derivation [closed]

Can any one help me in deriving the matrix inversion lemma rule for RLS algorithm? I don't know how to start with. Many books have just stated but they haven't derived it.
1
vote
0answers
116 views

Instability problem in NLMS adaptive algorithm

I have been used NLMS algorithm to equalize 4x4 MIMO signals, but the bit-error-rate (BER) after equalization is unstable with iterations. I don't know if it is the normal behavior of the adaptive ...
1
vote
0answers
57 views

Solving an array signal processing estimation problem based on the Rayleigh quotient

The Rayleigh quotient for a covariance matrix $\mathbf{C}$ and a non-zero steering vector $\mathbf{a}$ is given by $$ R(\mathbf{C},\mathbf{a}) := \frac{\mathbf{a}^H\mathbf{C}\mathbf{a}}{\mathbf{a}^H\...
1
vote
0answers
81 views

How to fit a function (find a lower bound) to (for) the plots?

I am dealing with a complicated optimization problem with two parameters $N, K$. Setting $K =$ constant and finding the minimum for different $N$'s lead to a set of function values. I can do this for $...
0
votes
2answers
91 views

Solving LASSO ($ {L}_{1} $ Regularized Least Squares) with Gradient Descent

To the best of my knowledge, state of the art methods for optimizing the LASSO objective function include the LARS algorithm and proximal gradient methods. I was wondering however, if the LASSO ...
0
votes
2answers
337 views

Why Is Non Linear Least Squares Method from MATLAB and Alglib Gives Different Results on the Same Data?

i'm trying to rewrite my Matalab prototype for some DSP to C++ and encountering a displeasing problem. I'm trying to fit data to a function $y = a * (\pi / 2 + arctg(b * x))$. In Matlab it works well ...