11 svar
278 visningar
andreaas10 behöver inte mer hjälp
andreaas10 26
Postad: 25 jun 2022 10:29

Optimering tre variabler och två bivillkor

Jag har försökt optimera 

f(x,y,z)=xyz med bivillkoren g(x,y,z)=x+y+z=1 och h(x,y,z)=x2+y2+z2=1.

f=(yz,xz,xy), g=(1,1,1) , h=(2x,2y,2z).

Jag försökte först använda determinanten det(f,g,h)=yz(2z2y)xz(2z2x)+xy(2y2x)=0, men det blir bara väldigt rörigt och jag kan inte förenkla determinanten med radoperationer på ett bra sätt.

Har också testat lösa ut från g. z=1-x-y och optimera

 f(x,y,1-x-y)=xy-x2y-xy2 med bivillkoret h(x,y,1-x-y)=2x2+2y2+2xy-2x-2y+1=1

med Lagrange multiplikator men de partiella derivatorna gör det för svårt.

D4NIEL 2932
Postad: 25 jun 2022 15:06 Redigerad: 25 jun 2022 15:07

Jag tycker att du kan börja med determinanten eftersom det blir så stiligt.

detf,g,h=2x2y+2x2z+2xy2-2xz2-2y2z+2yz2=-2(x-y)(x-z)(y-z)=0\det\left [\nabla f, \nabla g, \nabla h\right]=2 x^2 y+2 x^2 z+2 x y^2-2 x z^2-2 y^2 z+2 y z^2=-2(x-y)(x-z)(y-z)=0

Vi ser att determinanten är noll då x=y=tx=y=t eller x=z=ux=z=u eller y=z=vy=z=v. Om vi börjar med x=y=tx=y=t ger h(x,y,z)=1h(x,y,z)=1

2t2+(1-2t)2=12t^2+(1-2t)^2=1

Osv...

För metoden med Lagrange mulitplikator är det bara att kötta på, återigen kan en insikt om hur man kan skriva ekvationerna i faktorer vara till stor hjälp.

D4NIEL 2932
Postad: 25 jun 2022 15:22 Redigerad: 25 jun 2022 15:26

Jag kanske också ska nämna att de två restriktionerna utgör ett plan och en sfär. Dina extrempunkter ska alltså ligga på skärningslinjen, vilken kan parametriseras. Ditt problem reduceras då till att att optimera löpvariabeln vilket kan vara ett roligt tredje alternativ till dina föreslagna lösningsmetoder. Enklast är kanske att införa ett sfäriskt koordinatsystem med "z-axeln", dvs den axel mot vilken du mäter azimut,  i planets normalriktning.

tomast80 4245
Postad: 25 jun 2022 17:26

Funderar lite på denna uppgift, är den verkligen rätt formulerad/avskriven? Svaret blir lite skumt. Funderar på om något av bivillkoren skulle varit lite annorlunda. Blir ett mer intressant svar också om man istället söker minimum av funktionen.

andreaas10 26
Postad: 25 jun 2022 17:46
tomast80 skrev:

Funderar lite på denna uppgift, är den verkligen rätt formulerad/avskriven? Svaret blir lite skumt. Funderar på om något av bivillkoren skulle varit lite annorlunda. Blir ett mer intressant svar också om man istället söker minimum av funktionen.

Uppgiften var egentligen optimera f med bivillkoren g=1 och h≤1 men jag hittade punkten (1/3,1/3,1/3) där när jag undersökte området g=1 och h<1. Efter det hade jag bara området då g=h=1 att undersöka och det var det jag frågade om.

Vad får du för svar? Största värdet är 0 i de triviala punkterna (1,0,0) osv och minsta -4/27 i (2/3,2/3,-1/3), (2/3,-1/3,2/3) och (-1/3,2/3,2/3) med bivillkoren g=h=1.

Med billkoren h≤1=g finns största värdet 1/27 i (1/3,1/3,1/3).

tomast80 4245
Postad: 25 jun 2022 17:54

Då förstår jag. Jag får samma svar. En sista fundering, blir det verkligen en skärningslinje mellan klotet och planet? Det känns mer som de skär varandra i enstaka punkter, har jag fel? Tycker WolframAlpha indikerar det med.

tomast80 4245
Postad: 25 jun 2022 17:58
tomast80 skrev:

Då förstår jag. Jag får samma svar. En sista fundering, blir det verkligen en skärningslinje mellan klotet och planet? Det känns mer som de skär varandra i enstaka punkter, har jag fel? Tycker WolframAlpha indikerar det med.

Hade tydligen fel, det blir en ellips.

Smaragdalena 80504 – Avstängd
Postad: 25 jun 2022 18:01

Hur skulle man kunna skära ett klot med ett plan så att skärningsytan blir något annat än en cirkelskiva? Jag kan inte föreställa mig hur det skulle gå till.  (Jo, om de bara skär varandra i en enda punkt, men det är en tråkig lösning...)

andreaas10 26
Postad: 25 jun 2022 18:16
D4NIEL skrev:

Jag tycker att du kan börja med determinanten eftersom det blir så stiligt.

detf,g,h=2x2y+2x2z+2xy2-2xz2-2y2z+2yz2=-2(x-y)(x-z)(y-z)=0\det\left [\nabla f, \nabla g, \nabla h\right]=2 x^2 y+2 x^2 z+2 x y^2-2 x z^2-2 y^2 z+2 y z^2=-2(x-y)(x-z)(y-z)=0

Vi ser att determinanten är noll då x=y=tx=y=t eller x=z=ux=z=u eller y=z=vy=z=v. Om vi börjar med x=y=tx=y=t ger h(x,y,z)=1h(x,y,z)=1

2t2+(1-2t)2=12t^2+(1-2t)^2=1

Osv...

För metoden med Lagrange mulitplikator är det bara att kötta på, återigen kan en insikt om hur man kan skriva ekvationerna i faktorer vara till stor hjälp.

Tack, jag löste uppgiften med determinanten men hur tänkte du när du faktoriserade? Det ska stå -2x2y förresten.

tomast80 4245
Postad: 25 jun 2022 18:23
Smaragdalena skrev:

Hur skulle man kunna skära ett klot med ett plan så att skärningsytan blir något annat än en cirkelskiva? Jag kan inte föreställa mig hur det skulle gå till.  (Jo, om de bara skär varandra i en enda punkt, men det är en tråkig lösning...)

Ja, kom på det. Men såg lite skumt ut på WolframAlpha.

tomast80 4245
Postad: 25 jun 2022 18:26

D4NIEL 2932
Postad: 25 jun 2022 18:31 Redigerad: 25 jun 2022 18:48

Om du undrar hur man kommer på att man kan faktorisera så beror det förmodligen mest på lite metafakta.

0. Det var väldigt vad uttrycket ser symmetriskt och städat ut? Det här borde gå att faktorisera!

1. Ursprungsekvationerna är symmetriska i x och y, dvs om du byter namn på x till y och y till x får du exakt samma uttryck för f,gf,g och hh

En eventuell lösning har alltså egenskapen x=yx=y

2. "Alltså måste (x-y) vara en faktor i determinantens nollprodukt. Intressant. Vad får jag kvar om jag löser ut (x-y)"... osv

Om det känns krångligt kan du ju alltid chansa på att  lösa ut -2x-2x och notera att du kan kombinera för att få ut -2(x-y)...

Så här ser sfären och planet ut när man tittar ned på planet  (vi befinner oss i en punkt ovanför planet och tittar ned i riktningen (-1,-1,-1)(-1,-1,-1):

Svara
Close