Matemaattinen optimointi

Matemaattinen optimointi tarkoittaa määritellyn parhaan ratkaisun valintaa kaikkien mahdollisten ratkaisujen joukosta. Yleensä matemaattinen optimointi liittyy matemaattiseen malliin, kuten kohde-, hyöty- tai kustannusfunktioihin, joiden tuntemattomia parametrejä voidaan hakea tai sovittaa dataan . Kun kohdefunktio kuvaa ratkaisusta saatavaa hyötyä, jonka halutaan olevan mahdollisimman suuri, kutsutaan optimointitehtävää maksimoinniksi. Kun taas halutaan ratkaisusta koituvan mahdollisimman vähän kustannuksia tai haittaa, kutsutaan tehtävää minimoinniksi.

Formaalisti optimointi on sellaisen pisteen x∗ etsiminen ratkaisujoukosta A, jossa funktio f:A→ℝ saa joko pienimmän tai suurimman arvonsa eli minimin tai maksimin. Tätä pistettä x∗ kutsutaan vastaavasti minimipisteeksi tai maksimipisteeksi.

Jokaista maksimointiongelmaa vastaa tietty minimointiongelma, joka ratkaisee maksimointiongelman. Funktion f maksimointi on sama tehtävä kuin funktion −f minimointi: minxf(x)=−maxx−f(x) . Näin ollen matemaattisen optimointiteorian riittää tarkastella vain minimointiongelmaa.

Merkintätapa

Tehtävää, jonka tarkoituksena on etsiä f pienin arvo, kutsutaan minimointitehtäväksi. Tehtäväkuvaus merkitään formaalisti

minx∈ℝf(x)

jossa f on tehtävän kohdefunktio.

Vastaavasti kohdefunktion maksimointitehtävää merkitään

maxx∈ℝf(x)

Optimointitehtävässä

min⁡f(𝐱) , 𝐱∈S

joukkoa S kutsutaan käyväksi ratkaisujoukoksi tai vain käyväksi joukoksi. Käypä joukko voidaan määritellä esimerkiksi

S=[0,5]

Monen muuttujan tehtävässä käypä ratkaisujoukko voidaan määritellä esimerkiksi

S={(x,y)|x>0 ja y<=5−2x}

Funktiolla voi olla yksi tai useampi paikallinen optimipiste, lokaali optimi. Lokaali minimi määritellään pisteeksi x∗∈S, jolle on olemassa ε>0 siten, että f(x∗)≤f(x) kaikilla x∈S, joilla ‖x−x∗‖<ε. Toisin sanoen on olemassa jokin alue eli ympäristö, jossa piste x∗ antaa pienemmän funktion arvon kuin muut ympäröivät pisteet.

Optimointitehtävän globaali minimi on pienin likaaleista pisteistä, eli piste 𝐱∗, jolle pätee f(𝐱∗)≤f(𝐱) kaikilla 𝐱, jotka kuuluvat käypään alueeseen . Tätä pienempää funktion arvoa ei voida saavuttaa käyvässä joukossa.

Optimointitehtävätyyppejä

Lineaarinen optimointi

Lineaarinen optimointi tarkoittaa optimointia kun kohdefunktio ja käypää aluetta rajoittavat ehdot ovat lineaarisia. Lineaarista optimointia kutsutaan myös lineaariseksi ohjelmoinniksi. Yleinen lineaarinen tehtävä voidaan esittää muodossa

min⁡𝐜T𝐱
𝐀𝟏𝐱≤𝐛𝟏
𝐀𝟐𝐱≥𝐛𝟐
𝐀𝟑𝐱=𝐛𝟑

Tässä ≤ ja ≥ merkeillä tarkoitetaan, että jokaista alkiota verrataan riveittäin toisiinsa. Voidaan osoittaa, että kaikki lineaariset optimointitehtävät voidaan esittään ali- ja ylijäämämuuttujien (engl. slack variable) avulla ns. standardimuodossa

min⁡𝐜T𝐱
𝐀𝐱=𝐛
𝐱≥𝟎

jossa 𝐜∈ℝn×1, 𝐛∈ℝm×1 ja 𝐀∈ℝm×n. Eli aina kun tarkastellan yleistä lineaarista tehtävää, voidaan tarkastella pelkästään standarditehtävää menettämättä tuloksien yleispätevyyttä. Huomaa, että myös vektorit 𝐜 ja 𝐱 muuttuvat tehtävätyypin muunnoksessa.

Tuloksia

  • Lineaarisen optimointitehtävän käypä alue on n dimensioisen avaruuden monitahokas (engl. polyhedra).
  • Oletetaan, että tehtävänä on minimoida lineaarista kohdefunktiota epätyhjän monitahokkaan sisällä (eli kyseessä on lineaarinen optimointitehtävä). Tällöin kohdefunktion optimaalinen arvo on −∞ tai on olemassa optimaalinen ratkaisu 𝐱∗. Huomaa, että optimipisteen 𝐱∗ yksikäsitteisyyttä ei ole taattu.
  • Jos lineaarisen optimointitehtävän ratkaisu 𝐱∗ on äärellinen, tulee yksi ratkaisu löytymään jostain rajoittavan monitahokkaan S kulmasta. Jos kaksi pistettä 𝐱1 ja 𝐱2 ovat optimaalisia ovat myös pisteet muotoa λ𝐱1+(1−λ)𝐱2 ,λ∈(0,1) optimaalisia. Tämän voi tulkita siten, että kahden kulman välillä oleva särmä on optimaalinen, jos kummatkin kulmat ovat optimaalisia. Todistus: Oletetaan, että annetut pisteet minimoivat kohdefunktion f. Voidaan osoittaa, että monitahokas on konveksi joukko, eli kaikilla 𝐱𝟏,𝐱𝟐∈S pätee λ𝐱1+(1−λ)𝐱2∈S ,λ∈(0,1). Koska annetut pisteet ovat optimaalisia, eli niitä pienempiä arvoja funktio ei voi monitahokkaassa saada, päteef(𝐱1)=f(𝐱2)=:z. Toisaalta f(λ𝐱1+(1−λ)𝐱2)=𝐜T(λ𝐱1+(1−λ)𝐱2)=𝐜Tλ𝐱1⏟=λf(𝐱1)+𝐜T(1−λ)𝐱2⏟=(1−λ)f(𝐱2), mikä voidaan kirjoittaa f(λ𝐱1+(1−λ)𝐱2)=λz+(1−λ)z=z . ◻. Helposti voidaan huomata, että tapaus voidaan yleistää koskemaan n:ää pistettä, eli jos pisteet 𝐱i,i=1,...n ovat optimaalisia, niin ovat myös niiden ns. konveksikombinaatiot ∑i=1nan𝐱n,∑i=1nan=1,an≥0 myös optimaalisia. Todistus on analoginen kahden pisteen tapauksen kanssa.

Esimerkkejä

Lineaarinen optimointitehtävä

min⁡ −x1+x2
x1+x2≤0
x1≥1

voidaan esittää standardimuodossa muunnoksella xi=xi+−xi−, jossa xi+ ja xi− ≥0. Kun lisätään vielä ali- ja ylijäämämuuttujat s1 ja s2 saadaan tehtäväksi

min⁡−x1++x1−+x2+−x2−
x1+−x1−+x2+−x2−+s1+0⋅s2=0
x1+−x1−+0⋅x2+−0⋅x2−+0⋅s1−s2=1
xi±,si≥1 ,i=1,2

Jos siis valitaan 𝐜=[−1,1,1,−1,0,0]T, 𝐛=[0,1]T, 𝐱=[x1+,x1−,x2+,x2−,s1,s2]T ja

𝐀=[1−11−1101−1000−1] ,

on tehtävä saatu standardimuotoon. Huomaa, että optimipisteessä vain toinen muuttujista xi+ ja xi− on nollasta poikkeava, joten ylimääräiset muuttujat eivät vaikuta rajoitusehtoihin.

Epälineaarinen optimointi

Epälineaarisella optimoinnilla tarkoitetaan sellaisen optimointitehtävän ratkaisemista, joka on muotoa

min⁡f(𝐱)
gi(𝐱)≤0,i=1,...,m
hj(𝐱)=0,j=1,...,l
x∈X⊆ℝn

jossa f,gi,hj:ℝn↦ℝ, f on mielivaltainen kohdefunktio, funktiot gi epäyhtälörajoitukset, hj yhtälörajoitukset ja X käypä joukko. Vaikka yhtälörajoitukset voidaan esittää epäyhtälörajoitusten avulla (hi≤0 ja −hi≤0), tämä ei käytännössä aina ole mahdollista, sillä jotkin ratkaisumenetelmät olettavat tehtävän rajoitusten olevan optimipisteessä lineaarisesti riippumattomat.

Epälineaarinen ja lineaarinen tehtävä (LP) eroavat toisistaan seuraavasti. Epälineaarisessa tehtävässä käypä alue voi olla mielivaltainen, kun LP-tehtävän käypä joukko on aina monitahokas. Toiseksi LP-tehtävän lineaarinen kohdefunktio on aina myös konveksi, joten tehtävälle löydetään globaali optimi hyödyntämällä konveksin optimoinnin tuloksia. Epälineaarinen tehtävä on lineaarista tehtävää yleisempi, joten kaikki tulokset, jotka on johdettu epälineaariselle tehtävälle, pätevät myös lineaariselle tehtävälle.

Välttämättömät optimaalisuusehdot

Kun jokin optimointitehtävän käyvän alueen piste on optimaalinen, se täyttää välttämättömät optimaalisuusehdot. Lause ei päde toisin päin, eli tehtävälle voi olla mahdollista löytää optimaalisuusehdot täyttävä piste, joka ei välttämättä ole optimipiste.

Rajoittamattoman tehtävän välttämätön optimaalisuusehto optimointitehtävän ratkaisupisteelle 𝐱∗ on

∇f(𝐱∗)=0

Yleiselle rajoitetulle tehtävälle voidaan johtaa Karush-Kuhn-Tuckerin välttämättömät optimaalisuusehdot (välttämättömät KKT-ehdot), ja hieman yleisemmät Fritz Johnin välttämättömät optimaalisuusehdot (välttämättömät FJ-ehdot).

Fritz Johnin välttämättömät optimaalisuusehdot

Olkoon käypä joukko epätyhjä ja avoin sekä 𝐱 jonkin epälineaarisen tehtävän lokaali optimi. Tällöin on olemassa vektori (u0,𝐮,𝐯) siten, että

u0∇f(𝐱)+∑i=1mui∇gi(𝐱)+∑i=1lvi∇hi(𝐱)=0
uigi(𝐱)=0,i=1,...,m
u0≥0,ui≥0,i=1,...,m

jossa 𝐮 ja 𝐯 ovat m ja l vektoreita joiden i:nnet komponintit ovat ui ja vi. Luonnollisesti on oletettava, että gradientit ovat olemassa. FJ-ehdot saadaan toteutumaan mielivaltaiselle käyvälle pisteelle lisäämällä sopivia merkityksettömiä rajoituksia.

Karush-Kuhn-Tuckerin välttämättömät optimaalisuusehdot

Fritz Johnin ehdot redusoituvat tietyin ehdoin KKT-ehdoiksi. KKT-ehdot ovat FJ-ehtoja hyödyllisempiä, sillä ne karsivat pois turhia kandidaattipisteitä. Käytännössä näiden ehtojen on taattava, että FJ-ehtojen u0 on nollasta poikkeava, jolloin jakamalla sillä puolittain saadaan ns. KKT-ehdot.

Olkoon 𝐱 jonkin epälineaarisen tehtävän lokaali optimi, jolla on sopivat rajoitusehdot. Tällöin on olemassa vektori (𝐮,𝐯) siten, että

∇f(𝐱)+∑i=1mui∇gi(𝐱)+∑i=1lvi∇hi(𝐱)=0
uigi(𝐱)=0,i=1,...,m
ui≥0,i=1,...,m

jossa 𝐮 ja 𝐯 ovat m ja l vektoreita joiden i:nnet komponentit ovat ui ja vi.

Esimerkkejä

Lähteet

  1. 1 2 3 4 5 Rantala, Tapio: Optimointitehtävien ratkaiseminen. 3. painos. Helsinki: CSC, 2004. ISBN 952-9821-95-6 Teoksen verkkoversio (PDF) Viitattu 23.4.2026. (suomeksi)
  2. ↑ George B. Dantzig: Origins of the simplex method, s. 141–151. New York, NY, USA: Association for Computing Machinery, 1.6.1990. ISBN 978-0-201-50814-7 Teoksen verkkoversio Viitattu 23.4.2026.

Kirjallisuutta