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
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