30.07.2013 Views

3. Egenværdier og egenvektorer. 1.4. Lineær afbildning Repetition

3. Egenværdier og egenvektorer. 1.4. Lineær afbildning Repetition

3. Egenværdier og egenvektorer. 1.4. Lineær afbildning Repetition

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Modellering <strong>og</strong> pr<strong>og</strong>rammering. <strong>Lineær</strong> algebra 31<br />

y ′ (t) = Ay(t) for t>0, y(0) = β<br />

Antag, at <strong>egenvektorer</strong>ne for A er lineært uafhængige.<br />

De udgør da en basis i R n ,<strong>og</strong>vikanskrive<br />

y(t) = ˜y1(t)v1+˜y2(t)v2+···+˜yn(t)vn = V ˜y(t)<br />

Det følger, at<br />

<strong>og</strong><br />

y ′ (t) = ˜y ′ 1 (t)v1+˜y ′ 2 (t)v2+···+˜y ′ n (t)vn = V ˜y ′ (t)<br />

˜y ′ (t) = V −1 AV ˜y(t) = V −1 V ΛV −1 V ˜y(t) = Λ˜y(t),<br />

som skal løses med begyndelsesbetingelsen ˜y(0) = ˜ β = V −1 β<br />

Simpelt problem: for j =1,2,...,n<br />

˜y ′ j (t) = λj˜yj(t) for t>0, ˜yj(0) = ˜ βj<br />

Løsning ˜yj(t) = ˜ βje λjt , j =1,2,...,n<br />

Modellering <strong>og</strong> pr<strong>og</strong>rammering. <strong>Lineær</strong> algebra 32<br />

Eksempel <strong>3.</strong>7. (compartment model)<br />

y ′ <br />

<br />

<br />

(t) =<br />

−0.18<br />

0.024<br />

0.15<br />

−0.27<br />

y(t) for t>0, y(0) =<br />

0<br />

1<br />

Egenløsninger for A:<br />

<br />

λ1 = −0.15 , v1 =<br />

0.981<br />

0.196<br />

λ2 = −0.30 , v2 =<br />

<br />

−0.781<br />

0.625<br />

Løsning:<br />

y1(t) = e −0.15t − e −0.30t<br />

y2(t) = 0.2e −0.15t +0.8e −0.30t<br />

y<br />

1<br />

0.8<br />

0.6<br />

0.4<br />

0.2<br />

0<br />

0 4 8 12 16 20 24<br />

t<br />

y 1<br />

y 2

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!