Harj. 9 AV teht. 1
20.11.2005 HA
1.
> | restart: |
> | with(LinearAlgebra):with(linalg): |
Warning, the previous binding of the name GramSchmidt has been removed and it now has an assigned value
Warning, the protected names norm and trace have been redefined and unprotected
> | A:=<<1,9>|<-1,-5>>; |
Harj. 9 AV tehtävässä 3 saatiin HY:n yleinen ratkaisu , otetaan siitä (eräs) perusmatriisi latomalla y1(t) ja y2(t)
sarakkeiksi. Näin saadaan:
> | Y:=exp(-2*t)*<<1,3>|<t,3*t-1>>; |
Luennolla esitetyn määritelmän mukaan:
> | Y0:=eval(subs(t=0,Y));Y0I:=MatrixInverse(Y0); |
> | E:=t->exp(-2*t)*(<<1,3>|<t,3*t-1>>.Y0I); |
> | E(t); |
Tarkistetaan ja todetaan, että oikein meni:
> | Matrix(exponential(A,t)); |
Sitten lasketaan:
> | y0:=<1,-1>; # Alkuehto |
> | yh:=E(t).y0; # (HY):n osuus |
> | yh:=exp(-2*t)*((<<1,3>|<t,3*t-1>>.Y0I).y0); # Sievenettynä |
(EHY)-osuus:
> | EHYintgroitava:=map(simplify,E(t-s).<s,1>); |
> | yp:=map(int,EHYintgroitava,s=0..t); # (EHY):n erikoinen |
> | y:=yh+yp; y:=simplify(%); # Tässä koko ratkaisu. |
> |
Tarkistus:
> | subs(t=0,y); |
> | vasen:=map(diff,y,t); |
> | oikea:=A.y+<t,1>; |
> |