6 svar
116 visningar
Ygolopot 215
Postad: 22 apr 2021 10:31 Redigerad: 22 apr 2021 11:30

System av ODE, uppgift med komplexa egenvektorer och egenvärden

Hej, behöver hjälp med följande uppgift:

x'(t)y'(t)=11-23x(t)y(t), jag har löst ut egenvärden och egenvektorer: λ1=2-i, v1= 1+i2, λ2=2+i, v2= 1-i2

Lösningen på systemet blir då: r(t)=x(t)y(t)=C1e(2-i)t1+i2 + C2e(2+i)t1-i2 vilket sen kan skrivas om till:

r(t)=C1e2t(cost-isint)1+i2 + C2e2t(cost+sint)1-i2  som i sin tur kan skrivas om till:

r(t)=C1e2tcost+sint2cost + C1ie2tcost-sint-2sint +C2e2tcost+sint2cost + C2ie2tsint-cost2sint 

Där realdelen av lösningen blir:

r(t)=x(t)y(t)=C1e2tcost+sint2cost +C2e2tcost+sint2cost  =e2t(C1+C2)cost+(C1+C2)sint(C1+C2)2cost 

Men svaret i facit är:

x(t)y(t)=e2t C1cost+C2sint(C1+C2)cost +(C2-C1)sint

Det står inte i uppgiften att vi ska ha just realdelen av lösningen men i.o.m att svaret inte har i med i sig så antar jag att det är realdelen vi ska ha.

Kan någon se var jag gör fel?

Tack på förhand!

EDIT: Sorry! hade fel plats på C2 och C1 i facit, så justerade det

Torgil 21 – Fd. Medlem
Postad: 22 apr 2021 10:51

Observera att C1 och C2 är komplexa.  Medför att C1i, C2i är komplexa konstanter. Varför tar du realdelen av lösningen ovan?  Fel.

Jobba vidare med att skriva om lösningen du har innan "Där realdelen.." och tänk på att facits konstanter INTE är lika med dina

konstanter.

Torgil 21 – Fd. Medlem
Postad: 22 apr 2021 10:52

och det du påstår som realdelen är inte realdelen då Cx:na inte reella.

Ygolopot 215
Postad: 22 apr 2021 11:00

Tack för svaret Torgil! Jag ska räkna om nu baserat på den informationen du gav mig. Snabb fråga bara, jag hade ingen aning om att C1 och C2 var komplexa, varför är dom det?

Min enda tidigare erfarenhet av system likt dessa är från Linjär Algebra och då hade vi alltid en startvektor x(0) där:

x(0)=c1v1+...+cnvn, där ci(0)=ci, i=1, ..., n , och  v1, ..., vn  egenvektorer. Därför har jag alltid tänkt att C1 och C2 är två godtyckliga konstanter som vi får ut först efter att vi har något startvärde x(0) att ta hänsyn till.

Ygolopot 215
Postad: 22 apr 2021 11:41 Redigerad: 22 apr 2021 11:53

Jag har inte riktigt förstått varför jag kan göra såhär ännu, men om jag bortser från alla så får jag att:

r(t)=C1e2tcost+sint2cost + C1e2tcost-sint-2sint +C2e2tcost+sint2cost + C2e2tsint-cost2sint =

2C1cost + 2C2 sint2C1 cost-2C1 sint+2C2cost+2C2sint  ={Sätt D1= 2C1 och D2=2C2}=D1cost+D2sint(D1+D2)cost + (D2-D1)sint

Så det stämmer ju :)

Edit:

Fattar nu (till slut...) att om jag kallar den del av lösningen där jag har egenvektor v1 för x(t) och den del av lösningen där jag har egenvektor v2 y(t) så ska jag ta: Re[x(t)] = (1/2)(x(t)+y(t)) och lägga ihop med Im[x(t)] = (1/2i) (x(t)+y(t)) på så sätt går det ihop och jag ska alltså INTE skriva om konstanterna som D1 = 2C1 som jag gjort ovan.

x(t)=C1e2tcost+sint2cost + C1ie2tcost-sint-2sint  

y(t)=C2e2tcost+sint2cost + C2ie2tsint-cost2sint 

Torgil 21 – Fd. Medlem
Postad: 22 apr 2021 13:16 Redigerad: 22 apr 2021 14:05

Tänk ist Cx:Na som komplexa. Varför komplexa. Tänk så här: I övrigt är alla termer komplexa så finns ingen orsak till att Cx:na inte komplexa. I så fall lägger du in ett bivillkor som är omotiverat och allt ska kunna motiveras! Att de vore reella kan också ses som begränsning/bivillkor som minskar de möjliga utfallen.

Ygolopot 215
Postad: 22 apr 2021 16:03

Aha okej då fattar jag, det fins helt enkelt ingen anledning att avgränsa sig till reella tal från början. Om jag får lov att ställa en till fråga för jag har uppenbarligen fortfarande inte fattat hur det här går till och jag har verkligen försökt förstå och lusläst det kompendium vi följer men det går inte ihop för mig.

Två bilder, dom ligger direkt efter varandra i texten, men kunde inte få med dom i en bild. Historien innan de två bilderna är att vi har en matris A=3-24-1 och vi har beräknat egenvektorerna och egenvärdena: λ1=1-2i, v1=11+i, λ2=1+2i, v2= 11-i

Där vi kallar: x*(t) =eλ1tv1och y*(t) =eλ1tv2

 

Det står: "Two linearly independent real solutions can be chosen as real and imaginary parts of x*(t) (or y*(t))"

Ändå, beräknar man Re[x*(t)] = (1/2)(x*(t)+y*(t)) och Im[x*(t)] = (1/2i)(x*(t)-y*(t)), alltså verkar det som att man beräknar det som en kombination av x* och y*. Men sen i exemplet (precis innan answer) verkar det ju verkligen vara x* or y* man gör såhär med, för i exemplet så avslutar man ju med att ta just:

x(t) = C1Re[x(t)] + C2Im[x(t)] och behandlar dom båda som reella. Jag förstår inte hur det här hänger ihop med beskrivningen av vad t.ex Re[x*(t)] = (1/2)(x*(t)+y*(t)) är. 

Svara
Close