Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.
Nástěnka
❗22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
❗04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
❗23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.
Nejste přihlášen(a). Přihlásit
Stránky: 1
Dobrý den,
Prosím o radu, jak nejlépe najít minimum následující cílové funkce:
Pokoušel jsem se tuto funkci minimalizovat klasicky pomocí Newton-Raphsonovy metody. Nejprve jsem tedy vytvořil její gradient a optimalizoval ho tak, aby se rovnal nulovému vektoru.
Pokud zvolím MM1=konst. a vykreslím závislost cílové funkce a její derivace dle EE1 na EE1, získám následující průběhy:
na levo: cílová funkce pro MM1=1
na pravo derivace cílové funkce dle EE1 pro MM1=1
Červeným puntíkem je označen bod kdy derivace cílové funkce dle EE1 se rovná 0 (tedy bod minima cílové funkce)
Pokud bude počáteční odhad v blízkosti červeného bodu, bude Newton-Raphsonova metoda správně směřovat do minima cílové funkce. Pokud dám počáteční odhad například 0.8 bude Newthon-Raphsonova metoda směřovat nesprávně.
Prosím o radu, jaká jiná optimalizační metoda by šla použít? Případně, pokud bych chtěl použít Newtonovu metodu, jak tento nešvar obejít?
Offline
Dobrý deň.
Ja by som to skúsil takto.
Funkcia má tvar . Keďže tie prvé dve veci sú nezáporné, táto funkcia bude mať minimum v bode, v ktorom sú "a" a "b" nulové. Samozrejme, ak taký bod existuje.
Napíšem si rovnice:
a druhá rovnica bude vyzerať podobne, proste chcem tú druhú zátvorku nulovú.
Skúsim to vyriešiť, a ak to má riešenie a patrí do definičného oboru tej funkcie, tak to je minimum.
Ak tie rovnice nemajú riešenie, tak by som použil štandard, teda gradient nulový ...
Niečo tam výjde:
http://www.wolframalpha.com/input/?i=47 … 3%5Ex-1%29
Nula nie je v definičnom obore, takže tam ostane to druhé číslo.
Offline
Ahoj ↑ kafe_arabica:
Pokud bys našel bod, ve kterém by byly obě závorky nulové, byla by to dost velká náhoda (obě ty funkce uvnitř jsou na kladné poloose monotonní). Zdá se mi, že ty numera pochází z nějakých reálných dat a proto musíme předpokládat obecnou situaci, kdy bod, který hledáš neexistuje. Ten tvůj výpočet z wolframu je irelevantní vzhledem k dané situaci - nikde nepředepisuješ nulovost těch výrazů.
↑ eso.rimrak:
Napadá mě akorát vzít tu derivaci a zkusit najít nějakou vhodnou zhlazovací funkci, kterou to přenásobíš. Ta zhlazovací funkce by měla splňovat přinejmenším tohle:
1) je kladná
2) blízko nuly bude dostatečně malá
To by mělo vyřešit problém toho, že to protíná tu osu moc "strmě". K tomu, aby ses trefil do správné části grafu samozřejmě ale nějakou apriorní informaci mít musíš.
Offline
↑ Bati: Ahoj.
No, to máš pravdu. Ja som to skúsil a vyšlo to. Určite to takto nejde vždy. Aj som tam napísal, "ak to má riešenie ...".
Za pokus to stojí, keďže výpočet trvá pol minúty.
Po vyjadrení autora, zistíme či treba počítať so všeobecnou situáciou alebo nie.
Offline
↑ kafe_arabica:
No tak, už jenom to, že tam má nějaký parametr naznačuje, že to nechce pro situaci . A nevím teda, co ti vyšlo, ale pro to nevychází, že by obě závorky mohly být nulové.
Offline
↑ Bati:
To číslo z wolframu, 3.07... je hodnota . K tomu treba dopočítať (sústava). A v tom bode je minimum.
Offline
Děkuji za odpovědi,
↑ kafe_arabica:
Tento způsob bohužel nepůjde použít. Uvedená cílová funkce je pouze ukázková. Úlohu potřebuji řešit obecně pro MM1, MM2, ..., MMn a EE1, EE2, ..., EEn, přičemž v závorkách se mohou EE a MM různě střídat. Splněno je vždy pouze to, že EE je exponentem ve jmenovateli a MM je v čitateli.
↑ Bati:
Pokusím se najít zhlazovací funci, ale problém je jednak strmost v bodě minima, ale také pomalý pokles směrem k vyšším hodnotám EE1. Jakmile se v nějaké iteraci dostane EE1 za ten "kopeček" přestává to fungovat
Pokusím se najít nějakou vhodnou funkci, kterou přičtu k cílové funkci, tak aby se minimum zachovalo a zároveň se funkce na pravé straně zvednula, případně zcela zahladila ten kopeček...
Offline
↑ eso.rimrak:
Chápu, to by šlo nejspíš zachránit tak, že by ta zhlazovací funkce rostla dostatečně rychle pro velká EE1.
Offline
Stránky: 1