Optimization Toolbox - Capitolul II

Imagine preview
(9/10 din 1 vot)

Acest curs prezinta Optimization Toolbox - Capitolul II.
Mai jos poate fi vizualizat un extras din document (aprox. 2 pagini).

Arhiva contine 20 fisiere doc, m de 22 de pagini (in total).

Profesor: Dan Selisteanu

Iti recomandam sa te uiti bine pe extras si pe imaginile oferite iar daca este ceea ce-ti trebuie pentru documentarea ta, il poti descarca.

Fratele cel mare te iubeste, acest download este gratuit. Yupyy!

Domeniu: Automatica

Extras din document

Cele mai mici patrate neliniare cu întreg. Modelul Jacobian de raritate.

Modelele de dimensiuni mari în Isqlin, Isqcurvefît sl fsolve pot fi folosite în problemele de mici dimensiuni pâna la probleme de dimensiuni medii fara a calcula Iacobinul în fun sau formarea modelului Jacobian de raritate. (Acest exemplu se apica de asemenea în care se folosesc frnincon sau fminunc fara a calcula Hessianul sau furnizarea modelului Hessian de raritate). Cât de mica este dimensiunea mica pâna la cea medie? Nu este disponibil nici un raspuns absolut, asa cum aceasta depinde de cantitatea virtuala a memoriei în configuratia sistemului calculator. Sa presupunem ca problema noastra are m ecuatii si n necunoscute. Daca comanda J=sparse(ones(m, n)) cauzeaza o eroare de depasire a memoriei calculatorului tau atunci aceasta problema este cu siguranta prea mare (depasita). Daca nu rezulta eroare, problema poate fi în continuare prea mare, dar putem sa aflam rulând si sa vedem daca MATLAB este capabil sa ruleze în cadrul cantitatii de memorie virtuala disponibila pe sistemul tau.

Sa presupunem ca avem o problema cu 10 ecuatii si doua necunoscute, ca de exemplu gaseste valoarea x care minimizeaza :

începând de la punctul x=[0.3, 0.4].

Pentru ca Isqnonlin presupune ca suma patratelor nu este formata explicit în functia folosita, functia introdusa Isqnonlin ar trebui în schimb sa calculeze vectorul functiei evaluate

, pentru k=l pâna la 10 (adica F ar trebui sa aiba k componente)

Pasul 1: Deschide un M-fisier myfun.m care calculeaza valorile functiei obiectiv

function f=myfun(x)

k=1:10;

f=2+2*k-exp(k*x(1))-exp(k*x(2));

Pasul2: Apeleaza programul celor mai mici patrate neliniare

x0=[0.3, 0.4]

[x, resnorm]=lsqnonlin('myfun',x0)

Întrucât Jacobianul nu este calculat în myfun.m (si parametrul Jacobian în options este "off" prin eroare) si nici un model de raritate Jacobian nu este format folosind parametrul JacobPattern în options, lsqnonlin va apela metoda de dimensiuni mari) eroarea pentru ldqnonlin, cu JacobPattern introdus în jstr=sparse(ones(10, 2)). Când programul de diferenta finita este apelat pentru prima data va detecta ca jstr este de fapt o matrice densa, adica nici un beneficiar al vitezei nu este dedus din magazinarea lui ca matrice rara, si ca de acum înainte va folosi jstr=ones(10, 2) (o matrice întreaga) pentru calculele optimizarii.

Dupa 24 de evaluari ale functiei acest exemplu da solutia:

x=

0.2578 0.2578

resnorm

resnorm=

124.3622

Majoritatea sistemelor calculator vor fi capabile sa lucreze cu probleme întregi mult mai mari, pâna la ecuatii si variabile de ordinul 100 .

Daca informatia exista, atunci se poate profita de existenta unei structuri de raritate în Jacobian (sau Hessian), metodele de mari dimensiuni vor fi întotdeauna rulate mai repede daca aceasta informatie este oferita.

Minimizarea neliniara cu gradient si Hessian cunoscut

Acest exemplu implica o rezolvare a unei probleme de minimizare neliniara cu o matrice H(x) Hessian tridiagonala, prima data calculata în mod explicit si apoi prin formarea structurii de raritate Hessian pentru problemele de diferenta finita. Problema este de a gasi x pentru a minimiza:

(1-7)

unde n=1000

Pasul1: Deschide un M-fisier brownfgh.m care calculeaza functia obiectiv, gradientul functiei obiectiv si matricea tridiagonala rara Hessian.

Acest fisier este mai degraba lung si nu este inclus aici. Poti vedea codul cu comanda

Fisiere in arhiva (20):

  • 1.Cele mai mici patrate neliniare cu intreg. Modelul Jacobian de raritate
  • 1.Cele mai mici patrate neliniare cu intreg. Modelul Jacobian de raritate
    • exemplul1.m
    • myfun.m
  • 10.Optimizarea obiectelor inline in schimbul M-fisierelor
  • 10.Optimizarea obiectelor inline in schimbul M-fisierelor
    • exemplul101.m
    • exemplul102.m
  • 2.Minimizarea neliniara cu gradient si Hessian cunoscut
  • 2.Minimizarea neliniara cu gradient si Hessian cunoscut
    • exemplul2.m
  • 3.Miminizarea neliniara cu gradient si model de raritate Hessian
  • 3.Miminizarea neliniara cu gradient si model de raritate Hessian
    • brownfg.m
    • exemplul3.m
  • 4.Minimizare neliniara cu restrictii de tip margini si preconditionar cu banda
  • 4.Minimizare neliniara cu restrictii de tip margini si preconditionar cu banda
    • exemplul40.m
    • exemplul41.m
    • exemplul42.m
  • 5.Minimizare neliniara cu Restrictii egalitati
  • 5.Minimizare neliniara cu Restrictii egalitati
    • exemplul5.m
  • 6.Minimizare patratica cu Restrictii margini simple
  • 6.Minimizare patratica cu Restrictii margini simple
    • exemplul61.m
    • exemplul62.m
    • exemplul63.m
  • 7.Cele mai mici patrate lineare cu restrictii cu margini simple
  • 7.Cele mai mici patrate lineare cu restrictii cu margini simple
    • exemplul71.m
    • exemplul72.m
  • 8.Programare lineara cu egalitati si inegalitati
  • 8.Programare lineara cu egalitati si inegalitati
    • exemplul8.m
  • 9.Programare lineara cu coloane pline in egalitati
  • 9.Programare lineara cu coloane pline in egalitati
    • exemplul9.m
  • Aplicatii.doc
  • Optimization Toolbox CAP II.doc

Alte informatii

Sunt si toate aplicatiile sub forma de fisiere cu extensia ".m". Aplicatiile meg, iar documentul propriu-zis poate fi aranjat in pagina dupa bunul vostru plac. A fost prezentat ca proiect la Facultatea de Automatica din Craiova