Ottimizzazione con Picard

Messaggioda Danidila » 30/11/2020, 22:31

Salve a tutti, ho un problema con lo svolgimento di questo esercizio. Devo trovare il valore delle due variabili I1 I2 tale che la funzione sia minima, usando come da titolo il metodo di picard. La funzione è la seguente:
Pat(I1, I2) =

11*I1 + (5*I2)/2 + (I2 - 1/500)*(500*I1 + 1000*I2 - 5/2) + (I1 - 1/1000)*(10000*I1 + 500*I2 - 11) - 2/125

Ora, andando a fissare una delle due variabili (a zero) e calcolando la derivata della funzione con una variabile fissata, mi trovo le due seguenti funzione derivate:
20000*I1 - 11
2000*I2 - 5/2

Il dubbio che mi sorge è: essendo la funzione lineare non contrattiva in quanto non rispetta le ipotesi di Banach-Caccioppoli, è possibile comunque minimizzarla oppure continuerà a divergere comunque io scelga un punto per iniziare l’algoritmo iterativo?
Danidila
Starting Member
Starting Member
 
Messaggio: 4 di 8
Iscritto il: 08/06/2016, 14:44

Re: Ottimizzazione con Picard

Messaggioda feddy » 01/12/2020, 13:17

Danidila ha scritto:Salve a tutti, ho un problema con lo svolgimento di questo esercizio. Devo trovare il valore delle due variabili I1 I2 tale che la funzione sia minima, usando come da titolo il metodo di picard. La funzione è la seguente:
Pat(I1, I2) =

11*I1 + (5*I2)/2 + (I2 - 1/500)*(500*I1 + 1000*I2 - 5/2) + (I1 - 1/1000)*(10000*I1 + 500*I2 - 11) - 2/125



Per risolvere il problema, come intendi applicare l'algoritmo di Picard? Credo tu intenda applicarlo all'equazione $\nabla f(l_1,l_2)=(0,0)$. È questo ciò che intendi?
Avatar utente
feddy
Moderatore
Moderatore
 
Messaggio: 2775 di 5934
Iscritto il: 26/06/2016, 00:25
Località: SISSA

Re: Ottimizzazione con Picard

Messaggioda feddy » 01/12/2020, 13:25

Giustamente credo tu voglia capire se l'algoritmo converge. Non ho fatto i conti, ma non mi pare tu stia controllando in modo corretto. Per utilizzare Picard devi porre il problema in forma "equivalente" come nel caso unodimensionale $$x = G(x)$$, dove $x=(l_1,l_2)$ e $G(x)$ è il gradiente della tua funzione. Ciò che devi controllare è che la jacobiana della funzione di iterazione abbia norma minore di 1 in un intorno del punto fisso. Meglio ancora se è $|| JG(x)||<1$ ovunque, in quel caso converge per ogni guess iniziale. Prova a fare i conti. E per favore scrivi usando i tag appositi per leggere il post in LaTeX :-)
Avatar utente
feddy
Moderatore
Moderatore
 
Messaggio: 2776 di 5934
Iscritto il: 26/06/2016, 00:25
Località: SISSA


Torna a Analisi Numerica e Ricerca Operativa

Chi c’è in linea

Visitano il forum: Nessuno e 1 ospite