Re: Minimizing the Frobenius norm



Peter,

So the thing that I am stumbling on is the fact that both Ti and Tj are
defined in terms of the unknown v:

Considering a single term in the inner summation, don't I really have
the following:

|| Ti * v - Tj * v||^2

I am not seeing how this turns into the familiar Ax - b = 0 type of
problem. If I knew either Ti or Tj then it easily fits into that form.
Treating (Ti - Tj) as the A matrix and b as a zero-filled vector
doesn't seem correct as the optimal value for x would always end up
being a zero-filled vector: Ax = 0

Am I missing the obvious?

P.S.

The v are linearly related to the T so it is indeed a linear
optimization problem.


the frobenius norm squared is the sum of the squares of the elements of the
matrix. hence you have here a typical "least squares" problem.
now it depends on the kind how your variables v_i enter the matrices T_i
what to do: should each of the elements of these matrices be an affine linear
function of the v_i, then you have a simple linear least squares problem and
can use one of the well established linear least squares solvers, for example
the one from LAPACK. (DGELSS)
Otherwise you have a nonlinear least squares problem and need a corresponding
code. there exist several good codes for this , for example ELSUNC is such ,
or , should the optimal sum of differences not be small,
port/n2f or port/n2g
anything in
http://www.netlib.org
or look here
http://plato.asu.edu/sub/nonlsq.html#lsqnres

hth
peter

.



Relevant Pages

  • Re: Minimizing the Frobenius norm
    ... The v are linearly related to the T so it is indeed a linear ... hence you have here a typical "least squares" problem. ... or, should the optimal sum of differences not be small, ... port/n2f or port/n2g ...
    (sci.math.num-analysis)
  • LMSE Estimation and Least Squares
    ... The solution of this Linear Mean Square Error (LMSE) error problem is ... Solution for this least squares problem is ... weight S by statistical considerations? ...
    (comp.dsp)
  • Re: Separable nonlinear least squares optimizer with simple bounds on linear and nonlinear variables
    ... solves a non-linear least squares problem, ... imposes simple bounds both on the linear and nonlinear variables, ... exploit the separability of my model. ...
    (sci.math.num-analysis)
  • Re: Linear Equations Similar to 1D Poisson
    ... how to deal with a linear system with matrix ... plus the constraint sum x=1 ... try to solve this harmless least squares problem directly, ... you write the coefficients of the linear equation on top of the matrix. ...
    (sci.math.num-analysis)
  • Re: Matrix / Estimation stuff
    ... > anything to do with the jacobian matrix above? ... linearise non-linear problems. ... linear least squares problem min|b-Ax| is to solve the normal ... the non-linear least squares problem minimises the linear model of the ...
    (sci.math.num-analysis)