principalMatrix.mws

diffgeometry[PrincipalMatrix] - Hlavní směry v tečné rovině plochy

Calling Sequence

PrincipalMatrix(X(u,v),u,v)

Parameters

X(u,v) - plocha daná vektorovou funkcí [x(u,v), y(u,v), z(u,v)] parametrů u, v.

u,v       - parametry plochy

Description

Examples

Nejprve je třeba nastavit cestu, kde máte uloženou knihovnu "diffgeometry", např

>    restart;

>    libname:=libname,"D:/Sarka/Maple/diffgeometry/libsarka":

>    with(diffgeometry):

Sféra je parametrizovaná hlavními křivkami, identicky platí F=M=0. Pro všechny rotační plochy jsou hlavní křivky rovnoběžkové kružnice a meridiány.

>    S:=[R*cos(u)*cos(v),R*sin(u)*cos(v),R*sin(v)];(PrincipalMatrix(S,u,v));

S := [R*cos(u)*cos(v), R*sin(u)*cos(v), R*sin(v)]

matrix([[x^2, -x*y, y^2], [R^2, 0, R^2*cos(v)^2], [-R, 0, -R*cos(v)^2]])

Hlavní  křivky helikoidu.

>    helicoid:=[u*cos(v),u*sin(v),v]:rovnice:=linalg[det](PrincipalMatrix(helicoid,u,v));

rovnice := -(-x^2+y^2+y^2*u^2)/(1+u^2)^(1/2)

Rovnici hlavních směrů řeší x = sqrt(u^2+1)*y  a x = - sqrt(u^2+1)*y  . po dosazení du/dt = x , dv/dt = y , řešíme diferenciální rovnici

>    x1:=solve(rovnice,x)[1];x2:=solve(rovnice,x)[2];

x1 := (1+u^2)^(1/2)*y

x2 := -(1+u^2)^(1/2)*y

>