V2/2001 Harj. 10 ohjeita
../H/ohjeita10.html      4.4.01

Harjoitustehtäviin 10 liittyviä ohjeita ja ehdotuksia

Huomaa, että jos olet kopioinut v201.mpl:n itsellesi, siellä ei ole ehkä uusimpia funktioita (kuten SVD), joten lue näin:
read("/p/edu/mat-1.414/maple/v201.mpl"):

Teht. 2 ja muut SD-hommat

Voi olla opettavaisempaa ja helpompaakin tehdä gradientit ilman grad-komentoa, suoraan osittaisderivoimalla. Alla on näin tehty ja saatu gradienttifunktio (jatkossa kätevämpi kuin lauseke).

Tässä on toimiva SD-askel komento komennolta. Hiukan piti kokeilla ja sen tuloksena lisätä mm. evalm. Tästä pitäisi olla helppo jatkaa. (Pari yllättävää kaatumista tapahtui. Muistakaa, että CTR-S tallettaa, kannattaa viljellä usein.)

Tässä pidetään rinnakkain

jonoa x0 ja Vectoria x0v
      x1             x1v     
          ....
Funktiota f voidaan soveltaa suoraan jonoon, kuten x0, x1 ...
Vektorit x0v, x1v ... soveltuvat vektorioperaatioihin, kuinkas muuten.

Jos vektorilasku ei mene komponenteille, täytyy komentaa vanhakantaiseen tyyliin evalm();

f:=(x,y)->4*x^2-4*x*y+2*y^2; g:=[D[1](f),D[2](f)];
x0:=1,1: x0v:=;
u:=-Vector(g(x0));
jana:=v2l(evalm(x0v+t*u));
phi:=simplify(f(op(jana)));
plot(phi,t=0..1);
t1:=solve(diff(phi,t)=0,t);
x1v:=x0v+t1*u;x1:=op(v2l(x1v));
f(x1),subs(t=t1,phi);