Path: news.mathworks.com!not-for-mail
From: "Bruno Luong" <b.luong@fogale.findmycountry>
Newsgroups: comp.soft-sys.matlab
Subject: Re: inverse matrix inv pinv linfactor ginv
Date: Sun, 4 Oct 2009 14:05:01 +0000 (UTC)
Organization: FOGALE nanotech
Lines: 18
Message-ID: <haaa2d$61e$1@fred.mathworks.com>
References: <ha8r6f$pf7$1@fred.mathworks.com> <haa339$1tk$1@fred.mathworks.com> <haa48p$in7$1@fred.mathworks.com> <haa8su$ige$1@fred.mathworks.com>
Reply-To: "Bruno Luong" <b.luong@fogale.findmycountry>
NNTP-Posting-Host: webapp-03-blr.mathworks.com
Content-Type: text/plain; charset="ISO-8859-1"
Content-Transfer-Encoding: 8bit
X-Trace: fred.mathworks.com 1254665101 6190 172.30.248.38 (4 Oct 2009 14:05:01 GMT)
X-Complaints-To: news@mathworks.com
NNTP-Posting-Date: Sun, 4 Oct 2009 14:05:01 +0000 (UTC)
X-Newsreader: MATLAB Central Newsreader 390839
Xref: news.mathworks.com comp.soft-sys.matlab:574805


"Petros Mpogiatzis" <painter@geo.auth.gr> wrote in message <haa8su$ige$1@fred.mathworks.com>...

> 
> or redefine your initial system to include the L'L constraint and then to solve the redefined system
> x=Aext\bext

I don't see any obvious redefinition beside the obvious one:
Aext = A'*A + L'*L
bext = A'*b

Note that the Tikhonov regularization above has another equivalent formulation:

minimizing || A*x - b ||
such that ||L*x|| <= delta

delta is "some" number

Bruno