Maple algebra

Enviado por Programa Chuletas y clasificado en Matemáticas

Escrito el en español con un tamaño de 5,26 KB

 
a) Matriz asociada a f en la base B. Hallar el rango y el núcleo de f.
>
A:=matrix([[3, -2, -2], [-2, 3, -1/3], [-2, -1/3, 3]]);
El rango de f es:
>
rank(A);
El núcleo son las soluciones de AX=0.
>
vcero:=vector([0,0,0]);
>
linsolve(A,vcero);
Las ecuaciones paramétricas del núcleo son (x1,x2,x3)=(4/3 a,a,a)
El nucleo es un subespacio de dimension 1, que tiene por ecuaciones cartesianas (respecto de la base B):
b) subespacio conjugado del vector u.
Es el conjunto de vectores z =( ) que cumplen que f(u,z)=0.
>
matrix(1,3,[1,0,0])&*A&*matrix(3,1,[y1,y2,y3]);

>
evalm(%)[1,1]=0;
El subespacio conjugado de u tiene por ecuacion cartesiana:
''es un plano'' 0=
c) Subespacio conjugado del plano x+y+z=0.
Elijo una base del plano: { (1,-1,0) , (0,1,-1) } y calculo los subespacios conjugados de cada uno de los vectores de la base.
>
v1:=[1,-1,0];v2:=[0,1,-1];
ecu1:=evalm(matrix(1,3,v1)&*A&*matrix(3,1,[y1,y2,y3]))[1,1]=0;
ecu2:=evalm(matrix(1,3,v2)&*A&*matrix(3,1,[y1,y2,y3]))[1,1]=0;
 solve({ecu1,ecu2},{y1,y2,y3});

Hemos obtenido una recta: y1=4/3 a, y2=a, y3= a. .
d) Vectores autoconjugados.
Comprobamos si los vectores verifican.
>
matrix(1,3,[x1,x2,x3])&*A&*matrix(3,1,[x1,x2,x3])=0;
>
w1:=[4,3,3];w2:=[5,3,3];
>
evalm(matrix(1,3,w1)&*A&*matrix(3,1,w1));
>
evalm(matrix(1,3,w2)&*A&*matrix(3,1,w2));
Por tanto el vector w1=(4,3,3) es autoconjugado respecto de la forma bilineal dada pero w2 no.
e) Base B' de vectores conjugados. B'={e1 , e2, e3}



> e1:=[1,0,0];
>
evalm(matrix(1,3,e1)&*A&*matrix(3,1,e1))[1,1];
No es autoconjugado, nos sirve.
>
conju1:=evalm(matrix(1,3,e1)&*A&*matrix(3,1,[y1,y2,y3]))[1,1]=0;
e2 tiene que pertenecer al subespacio conjugado de e1.
>
solve(conju1);
>
e2:=[1,1/2,1];
Subespacio conjugado de e2:
>
conju2:=evalm(matrix(1,3,e2)&*A&*matrix(3,1,[y1,y2,y3]))[1,1]=0;
e3 tiene que pertenecer a conju1 y a conju2:
>
solve({conju1,conju2},{y1,y2,y3});
>
e3:=[4/3,1,1];
La base de vectores conjugados es
>
{e1=evalm(e1),e2=evalm(e2),e3=evalm(e3)};
f) Matriz D asociada a f respecto de la base B'
> evalm(e1&*A&*e1);
>
evalm(e2&*A&*e2);
>
evalm(e3&*A&*e3);
>
D1:=diag(3,5/12,0);
g) Hallar la matriz de cambio P tal que D=P^t AP
matriz de cambio de B' a B
P:=transpose(matrix([e1,e2,e3]));
>
evalm(transpose(P)&*A&*P);
A=Q^t D1 Q
>
Q:=inverse(P);
>
evalm(transpose(Q)&*D1&*Q);evalm(A);
h) Ya hemos diagonaliozado ''por congruencia '' la matriz A, y hemos obtenido una matriz diagonal D1:
>
evalm(D1);
Como el rango es 2 y la signatura tambien es (2,0), podemos concluir que la forma cuadratica asociada a f es SEMIDEFINIDA POSITIVA.
>
charpoly(A, lambda);
>
eigenvalues(A);
>
eigenvectors(A);

Entradas relacionadas: