Practica2 matII 2022 23
Practica2 matII 2022 23
Practica2 matII 2022 23
Practica 2 Matlab
1. Objetivos
- Estudiar y resolver sistemas de ecuaciones lineales con parámetros.
En primer indicamos que el parametro m queremos que sea considerado como una
variable. Para a continuación definir las matrices del sistema y termino independiente.
syms m
A=[1 m 1; 1 1 m; m 1 1]
b=[m+2 ; -2m-2; m]
det(A)
ans= 2 − 3m + m3
Luego resolvemos dicho determinante para obtener los valores de m para los cuales
dicho determinante es nulo.
solve(ans)= -2,1,1
caso1= A\b
m/(m-1)
(m+2)/(m-1)
-2*(m+1)/(m-1)
ans=
caso2=subs([A b],m,-2)
Obteniendo la matriz
ans =
1 −2 1 0
1 1 −2 2
−2 1 1 −2
rref(caso2) nos saldrá una matriz donde debemos de indicar que el sistema es com-
patible indeterminado.
caso3=subs([A b],m,1)
Obteniendo la matriz
ans =
1 1 1 3
1 1 1 −4
1 1 1 1
3. Espacios vectoriales
Veremos como podemos usar matlab para el estudio de espacios y subespacios vecto-
riales.
U1 =< (3, 8, 11, 2, −11), (0, 2, 2, 2, −5), (3, 4, 7, −2, −1), (2, 2, 4, −2, 1) >
U2 = {(x, y, z, w, t) ∈ R5 / x − y = 0, y + w = 0}
Se pide:
b) ¿Pertenece el vector −
→
u = (−3, 1, 0, 1, 1) al subespacio U1 ? ¿y el vector −
→
v =
(3, 2, 5, −4, 4) ?
Apartado a
En primer lugar generamos una matriz con los vectores que generan al subespacio U1 .
Luego calculamos los vectores que son linealmente independientes con la instrucción
rref (U1 )
1,0000 0 1,0000 −2,0000 3,0000
0 1,0000 1,0000 1,0000 −2,5000
0 0 0 0 0
0 0 0 0 0
Como las únicas filas no nulas se corresponden con las filas 1 y 2 quiere decir que
dichos vectores son los generadores y una base de U1 .
BU1=ans([1,2],:)
BU1 =
[ 1, 0, 1, -2, 3]
[ 0, 1, 1, 1, -5/2]
syms x y z w t
M=[BU1;[x y z w t]]
M=
[ 1, 0, 1, -2, 3]
[ 0, 1, 1, 1, -5/2]
[ x, y, z, w, t]
Ahora para calcular las ecuaciones implı́citas hallamos los siguientes determinantes:
d1=det(M(:,[1,2,3]))
d1 =z-y-x
d2=det(M(:,[1,2,4]))
d2 =w-y+2*x
d3=det(M(:,[1,2,5]))
d3 =t+5/2*y-3*x
EU1=[d1,d2,d3]
Apartado b
Para resolver el apartado segundo podemos proceder de dos maneras distintas. Pode-
mos hallar el rango de la matriz formada por los vectors BU1 y el nuevo vector o sustituir
en las ecuaciones implı́citas las coordenadas del vector y comprobar si las satisface.
C=[BU1;[-3 1 0 1 1]]
C=
[ 1, 0, 1, -2, 3]
[ 0, 1, 1, 1, -5/2]
[ -3, 1, 0, 1, 1]
rank(C)=3
subs(EU1,{’x’,’y’,’z’,’w’,’t’},{-3,1,0,1,1}) o subs(EU1,{x,y,z,w,t},{-3,1,0,1,1})
Apartado c
s=solve(’x-y=0’,’y+w=0’,’x,y,z,w,t’)
s=
t: [1x1 sym]
w: [1x1 sym]
x: [1x1 sym]
y: [1x1 sym]
z: [1x1 sym]
s=[s.x,s.y,s.z,s.w,s.t]
s=
[ -w, -w, z, w, t]
v1=subs(s,{z,w,t},{1,0,0})
v1 = 0 0 1 0 0
v2=subs(s,{z,w,t},{0,1,0})
v2 = -1 -1 0 1 0
v3=subs(s,{z,w,t},{0,0,1})
v3 = 0 0 0 0 1
BU2=[v1;v2;v3]
BU2 =
0 0 1 0 0
−1 −1 0 1 0
0 0 0 0 1