Optimization toolbox - capitolul II

Curs
9/10 (1 vot)
Domeniu: Automatică
Conține 20 fișiere: doc, m
Pagini : 22 în total
Cuvinte : 5423
Mărime: 1.02MB (arhivat)
Cost: Gratis
Profesor îndrumător / Prezentat Profesorului: Dan Selisteanu
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

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

Preview document

Optimization toolbox - capitolul II - Pagina 1
Optimization toolbox - capitolul II - Pagina 2
Optimization toolbox - capitolul II - Pagina 3
Optimization toolbox - capitolul II - Pagina 4
Optimization toolbox - capitolul II - Pagina 5
Optimization toolbox - capitolul II - Pagina 6
Optimization toolbox - capitolul II - Pagina 7
Optimization toolbox - capitolul II - Pagina 8
Optimization toolbox - capitolul II - Pagina 9
Optimization toolbox - capitolul II - Pagina 10
Optimization toolbox - capitolul II - Pagina 11
Optimization toolbox - capitolul II - Pagina 12
Optimization toolbox - capitolul II - Pagina 13
Optimization toolbox - capitolul II - Pagina 14
Optimization toolbox - capitolul II - Pagina 15
Optimization toolbox - capitolul II - Pagina 16
Optimization toolbox - capitolul II - Pagina 17
Optimization toolbox - capitolul II - Pagina 18
Optimization toolbox - capitolul II - Pagina 19
Optimization toolbox - capitolul II - Pagina 20
Optimization toolbox - capitolul II - Pagina 21
Optimization toolbox - capitolul II - Pagina 22

Conținut arhivă zip

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

Alții au mai descărcat și

Sisteme convenționale pentru reglarea proceselor continue

Capitolul 2 Sisteme Conventionale pentru Reglarea Proceselor Continue Rezumat: In acest capitol sunt tratate aspecte legate de metodologia...

Teoria Sistemelor

Reprezentarea Sistemelor Dinamice Liniare Multivariabile prin Matrice de Transfer 1. Matricea de transfer; legatura cu reprezentarile de tip...

Reprezentarea Informațiilor cu Obiecte

Informatiile pe care le reprezentam în memoria calculatorului sunt rareori atât de simple precum culorile sau literele. În general, dorim sa...

Sistemele Informatice

1.1. Contextul actual La sfârsitul secolului al XX-lea si începutul secolului al XXI-lea, clientii, concurenta si schimbarea au creat o noua lume a...

Cursuri Java

Cuvinte importante: - concepte fundamentale ale programarii orientate obiect in Java: incapsulare, mostenire, polimorfism; - crearea claselor de...

Aplicatii de retea în internet

Posta electronica (e - mail) Milioane de oameni sunt conectati într-un fel sau altul la reteaua Internet si pot trimite mesaje prin intermediul...

Optimizarea Conducerii Autovehiculelor

Titlul acestui subcapitol sugereaza utilizarea unor tehnici si a unor sisteme de conducere de tipul celor mentionate în primul capitol care sa...

Arhitectura modelului OSI(ISO)

ARHITECTURA MODELULUI OSI/ISO Modelul ISO/OSI (International Standards Organization / Open Systems Interconnection) este o arhitectura de retea...

Te-ar putea interesa și

Optimization Toolbox - Capitolul 1

Introducere Cursul îsi propune sã prezinte Optimization Toolbox utilizând MATLAB în versiunea 2. 0. Optimization Toolbox este o multime de...

Ai nevoie de altceva?