Nnls lawson hanson method in linearized models
This presentation is the property of its rightful owner.
Sponsored Links
1 / 20

NNLS (Lawson-Hanson) method in linearized models PowerPoint PPT Presentation


  • 127 Views
  • Uploaded on
  • Presentation posted in: General

NNLS (Lawson-Hanson) method in linearized models. LSI & NNLS. LSI = Least square with linear equality constraints NNLS = nonnegative least square . Flowchart. Initial conditions. Sets Z and P Variables indexed in the set Z are held at value zero

Download Presentation

NNLS (Lawson-Hanson) method in linearized models

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript


Nnls lawson hanson method in linearized models

NNLS (Lawson-Hanson) method in linearized models


Lsi nnls

LSI & NNLS

  • LSI = Least square with linear equality constraints

  • NNLS = nonnegative least square


Flowchart

Flowchart


Initial conditions

Initial conditions

  • Sets Z and P

  • Variables indexed in the set Z are held at value zero

  • Variables indexed in the set P are free to take values different from zero

  • Initially and P:=NULL


Flowchart1

Flowchart


Stopping condition

Stopping condition

  • Start of the main loop

  • Dual vector

  • Stopping condition:

    set Z is empty or


Flowchart2

Flowchart


Manipulate indexes

Manipulate indexes

  • Based on dual vector, one parameter indexed in Z is chosen to be estimated

  • Index of this parameter is moved from set Z to set P


Flowchart3

Flowchart


Compute subproblem

Compute subproblem

  • Start of the inner loop

  • Subproblem

    where column j of Ep


Flowchart4

Flowchart


Nonnegativity conditions

Nonnegativity conditions

  • If z satisfies nonnegativity conditions then we set x:=z and jump to stopping condition

  • else continue


Flowchart5

Flowchart


Manipulating the solution

Manipulating the solution

  • x is moved towards z so that every parameter estimate stays positive. Indexes of estimates that are zero are moved from P to Z. The new subproblem is solved.


Testing the algorithm

Testing the algorithm

  • Ex. Values of polynomial

    are calculated at points x=1,2,3,4 with fixed p1 and p2.

  • Columns of E hold the values of polynomial y(x)=x and polynomial at points x=1,2,3,4.

  • Values of p1and p2 are estimated with NNLS.


Nnls lawson hanson method in linearized models

nnls_test 0.1 (c) 2003 by Turku PET Centre

Matrix E:

1 1

2 4

3 9

4 16

Vector f:

0.6 2.2 4.8 8.4

Result vector:0.1 0.5


Nnls lawson hanson method in linearized models

nnls_test 0.1 (c) 2003 by Turku PET Centre

Matrix E:

1 1 1

2 4 8

3 9 27

4 16 64

Vector f:

0.73 3.24 8.31 16.72

Result vector:0.1 0.5 0.13


Nnls lawson hanson method in linearized models

nnls_test 0.1 (c) 2003 by Turku PET Centre

Matrix E:

1 1 1 1

2 4 8 16

3 9 27 81

4 16 64 256

Vector f:

0.73 3.24 8.31 16.72

Result vector:0.1 0.5 0.13 0


Nnls lawson hanson method in linearized models

nnls_test 0.1 (c) 2003 by Turku PET Centre

Matrix E:

1 1 1

2 4 8

3 9 27

4 16 64

Vector f:

0.23 1.24 3.81 8.72

Result vector:0.1 7.26423e-16 0.13


Nnls lawson hanson method in linearized models

  • Kaisa Sederholm: Turku PET Centre Modelling report TPCMOD0020 2003-05-23


  • Login