MATLAB

TRABAJOS PRCTICOS
DE LGEBRA LINEAL
Para carreras de Ingeniera

F. C. E. y T.

U.N.S.E.

Mara Ins Morales de Barrionuevo

Trabajos Prcticos de lgebra Lineal con MATLAB

INTRODUCCION
Normalmente se requiere de modelos computacionales con el fin de resolver
problemas de ingeniera. Muchas veces puede ser til hacer un programa que utilice
matrices, complejos, y otras estructuras matemticas, fcil de escribir y revisar. MATLAB es
ideal para esto
El nombre MATLAB proviene de la contraccin de MATrix LABoratory. Es un sistema
de trabajo interactivo y una herramienta importante para cualquier tarea que requiera
clculos matriciales, ya sea que involucren ecuaciones, sistemas caractersticos, mnimos
cuadrados, etc. y la visualizacin grfica de los mismos. Se pueden resolver problemas
numricos relativamente complejos sin necesidad de escribir un programa para ello.
Tiene gran capacidad de expansin ya que permite que el usuario defina sus
propias funciones, es por ello que es considerado como un lenguaje de programacin
para clculos tcnicos y cientficos donde se encuentran implicados elevados clculos
matemticos y la visualizacin grfica.
MATLAB trabaja esencialmente con un solo tipo de objetos: matrices, con entradas
eventualmente complejas. Esto significa que todas las variables representan matrices, de
este modo los escalares quedan caracterizados por matrices 1x1 y los vectores, por
matrices de una fila o una columna.
Es importante destacar que el alcance y la potencia de MATLAB van ms all de lo
que pueda verse en esta serie de trabajos prcticos

Caractersticas bsicas:

Al ejecutarse MATLAB aparece una ventana, la ventana de orden, que es el lugar

fundamental donde se interacciona con MATLAB. Cuando esta ventana est activa
aparece a la izquierda el smbolo >> seguido del cursor, es el indicativo de peticin de
ordenes del MATLAB.

Operaciones aritmticas bsicas

En el cuadro de la pgina siguiente se ilustran las operaciones aritmticas bsicas,
el smbolo que emplea MATLAB y un ejemplo:

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 2 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Operacin

Smbolo

Ejemplo

Suma

3+6

Resta

12 7.5

Multiplicacin

23.2 5.12

Divisin

77 / 51

Potencia

5^ 2

El orden en que estas operaciones se evalan en una expresin est dada por la
siguiente regla de precedencia:
Las expresiones se evalan de izquierda a derecha, con la operacin de
potencia teniendo el orden de precedencia ms alto, seguido por multiplicacin y
divisin que tienen ambas igual precedencia y seguidas finalmente, por suma y resta
que tienen igual precedencia.
Se pueden emplear parntesis para alterar esta usual ordenacin, en cuyo caso la
evaluacin se inicia dentro del parntesis ms interno y procede hacia fuera.
Ejercicio:
Para observar como funcionan estas operaciones con MATLAB, resuelva:
a) 3^2 + 1 2 / 2 * 4
b) (3^2 + 1 2) / 2 * 4
c) 3^(2 + 1) 2 /( 2 * 4)

Variables:
En MATLAB como en cualquier otro lenguaje de programacin, y/o asistente
matemtico se utilizan variables. Las variables deben tener un nombre segn ciertas
reglas. Estas reglas son:

MATLAB distingue minsculas de maysculas, esto quiere decir que A y a representan
variables diferentes. Si se desea que su utilizacin sea indistinta se debe tipear la
sentencia casesen.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 3 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Las variables pueden contener hasta 19 caracteres y comenzar por una letra seguida
de cualquier nmero de letras, dgitos o guiones de subrayado
Ejemplo: suma, A, S3, CIRCULO, valor_propio, etc.

Los nombres de variables no pueden contener operadores ni puntos. (No es vlido
usar /, *, -, +, ...)
MATLAB cuenta con variables especiales, algunas de ellas son:
o

ans ( de answer respuesta-) es la variable que MATLAB crea

automticamente cuando hay una expresin a la que no se le asign
nombre de variable.
Ejemplo: 2 *5+1 .5
ans =
11.5000

da com o respuesta :

pi: Razn de una circunferencia y su dimetro

inf: Infinito. Por ejemplo:

a=1/0
Warning: Divide by zero.
a=
Inf

NaN: Magnitud no numrica. Por ejemplo:

b=0/0
Warning: Divide by zero.
b=
NaN

i y j : Se emplean para introducir nmeros complejos. Es decir: i = j =

Nota: los nmeros complejos se introducen siguiendo el convenio usual, vale decir que
z1 = 2 + 3i, se introduce como:
>> z1=2+3i

Ordenes de Propsito General

Para
recordar informacin: para poder visualizar cuales variables han sido ya
introducidas se teclea who

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 4 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Para conocer el contenido de una variable: simplemente se tipea su nombre y se oprime

enter.
Para eliminar una variable de la memoria se utiliza la instruccin clear seguida del
nombre de la variable( ej: clearA ). Si se escribe clear simplemente se borran todas las
variables.
Help: Si teclea help seguido de un comando, funcin o archivo de MATLAB, aparecer
una descripcin del comando.
Por ejemplo: help : dar informacin sobre como puede usar : en MATLAB
help rref dar informacin sobre el comando rref
Uso de flechas: Al emplear la flecha se desplazarn los comandos anteriores. Se
pueden utilizar las flechas , , , para localizar un comando y modificarlo y al oprimir
la tecla enter se ejecuta el comando modificado.
Comentarios: Si se inicia una lnea con el smbolo %, MATLAB interpretar como un
comentario.
Por ejemplo: % Este es un comentario
Supresin de vista en pantalla: Si se escribe una instruccin de MATLAB y no se desea ver
los resultados desplegados, se finaliza el comando con un ; (punto y coma). Esto es
esencial para evitar prdidas de tiempo al mostrar los resultados intermedios.
Cualquier tipo de clculo, grfico o impresin puede detenerse sin salir del programa con
CTRL-C.
Cuando se sale de MATLAB, se pierden todas las variables. Para evitarlo se puede utilizar
la instruccin save antes de salir. Esto hace que las variables se almacenen en el archivo
de disco matlab.mat. Al ingresar nuevamente a MATLAB se pueden recuperar todas las
variables con la instruccin load.

Ejercicios:
1. Realice las siguientes operaciones con MATLAB y verifique las respuestas.
a) 2.(9 3)+4

Rta.: 16

5 (1/2)2
b)

Rta: 2,7941
0,7 + 1
Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 5 _

Trabajos Prcticos de lgebra Lineal con MATLAB

c) (1 0,25)1/2 + (4/81)-1/2

d)

e)

256

(5 i )2 +
3i

(1/25) -2

2i

Rta: 5,3660

Rta: 0 +3i

( Nota: sqrt(x) calcula la raz

cuadrada de x )

Rta: -1.8780 - 8.3436i

2. Sean a = 4, b = -2, c = 2/3 calcule:

a) A = ab c
b) B = a(b c)

3. a) Ingrese la orden who para observar que variables ha introducido.

b) Guarde todas las variables, cierre el programa y a continuacin vuelva a abrirlo
c) Escriba el comando who y apriete enter. Qu ocurre?
d) Escriba el comando load y repita el paso (c)
e) Escriba el nombre de alguna variable y apriete enter para ver su contenido.
f) Guarde solo una de las variables, cierre el programa, vuelva a abrirlo y repita los
pasos c) y d).

ARCHIVOS M
MATLAB puede ejecutar una sucesin de instrucciones almacenadas en archivos de
disco. Estos archivos se denominan archivos-M debido a que su extensin deber ser m.
Los archivos-M pueden ser de comandos o de funciones.

Archivos-M de comandos
Estos archivos consisten en una sucesin de instrucciones de MATLAB y son guardados
como un simple archivo de texto (salvo que su extensin es m y no txt). Pueden ser usados
para introducir gran cantidad de datos o cuando se desea cambiar el valor de una o
ms variables y reevaluar una serie de rdenes.
Para elaborar un archivo-M se hace click en New del men File y se selecciona M-file.
Aparecer una nueva ventana para la edicin de texto donde se teclea la lista de

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 6 _

Trabajos Prcticos de lgebra Lineal con MATLAB

comandos de MATLAB. Para guardarlo como archivo-M en el disco, se hace clik en Save
as del men File y se escribe el nombre del archivo con la extensin m.
Escribiendo simplemente el nombre del archivo (sin la extensin m) en la lnea de orden
de MATLAB se ejecutarn, lnea por lnea las rdenes contenidas en el archivo.
Ejemplo: Un archivo-M para calcular el volumen de una esfera

% Volumen de una esfera

Lnea de comentario

r = 2.4;

Asigna a la variable r el valor del radio

de la esfera

V= (4/3)*pi*r^3;

Calcula el volumen y guarda el resultado

en la variable V

disp('El volumen de la esfera es:')

Imprime en pantalla el texto encerrado

entre comillas

Imprime el valor de V

Si este archivo es guardado con el nombre de esfera.m y luego se lo llama desde

la ventana de orden tecleando esfera, se obtiene:
El volumen de la esfera es:
V=
57.9058
Ejercicios:
a) Confeccione un archivoM que calcule el rea del tringulo cuya base y altura
miden respectivamente 12.34cm y 7.82cm. y gurdelo en el disco
b) Llmelo desde la ventana de orden para ver que obtiene.
c) Abra nuevamente el archivo haciendo click en Open M-file del men file y
luego de seleccionar el nombre del archivo haga click en Aceptar. Modifique
los valores de la base y altura del tringulo, guarde el archivo modificado y
ejectelo nuevamente desde la ventana de orden de MATLAB.

Archivos M de funcin
Estos archivos permiten que el usuario pueda definir sus propias funciones. Son similares a
los de comandos pues son archivos de texto que tienen una extensin m.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 7 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Los archivos-M de funcin, como los de comando admiten lneas de comentario que
comienzan con un % y no son activas al momento de ejecutarse.
La primera lnea activa (no comentario) tiene el siguiente formato obligatorio:
function[V1,V2,...,Vn]=nombre_archivo(a1,a2,...ak)
donde V1, V2, ...,Vn son las variables de salida y a1, a2,...,ak son los argumentos o
variables de entrada, que sern asignadas en el mismo orden al convocar la funcin.
Cuando se tipea help nombre_archivo aparece en pantalla el primer bloque de lneas de
comentario del archivo, en el cual usualmente se describe la funcin nombre_archivo.
Ejemplo:
El archivo esfera.m del ejemplo anterior puede ser modificado de manera tal que el
valor del radio de la esfera pueda ser ingresado directamente desde el teclado.
% vol-esf(r) calcula el volumen de una esfera de radio r
function V=vol_esf(r)
V = (4/3)pir^3;
disp('El volumen de la esfera es:')
Notas:
Este archivo debe guardarse con el nombre de vol_esf.m
Las lneas de comentario del archivo de funcin, constituyen el texto visualizado en
respuesta a la orden help vol_esf.

Ejercicio 1:
a) Guarde el archivo vol_esf en el disco y desde la ventana de orden de MATLAB
escriba V=vol_esf(5).
b) Halle el volumen para diferentes valores de r.
Ejercicio 2
Construya un archivo de funcin que calcule el valor de la funcin:
f(x) = 2x3 + 2x2 + 5x + 1

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 8 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Trabajo Prctico N1
Introduccin de Matrices. Operaciones
1- Los siguientes ejemplos definen diferentes formas de introducir matrices en MATLAB.
Pruebe y saque conclusiones.
(CUIDADO!! ... los espacios en blanco tienen su significado)
a) A=[1,2,3;4,5,6;7,8,9]
b) B=[1 2 3;4 5 6;7 8 9]
c) C=[1 2 3
4 5 6
7 8 9]
d) D=C
e) E=[3 0 5 6]
f) F=[1:9]
g) G=[1:9]
h) H =[0:2:12]
i) I=[-1,3 sqrt(3) 6.4]
j) K=[-1 2+i;i -1]
Para que se utiliza , ?
Para que se utiliza ; ?
Para que se utiliza : ?
Para que se utilizan los espacios en blanco?
Que ocurre cuando al final de una matriz se escribe ?Qu significa?
2- Existen funciones en MATLAB como por ejemplo: rand, magic, eye, etc. que
proporcionan una forma sencilla para crear matrices con las que se puede
experimentar. Introduzca las siguientes matrices y extraiga conclusiones como por
ejemplo que la instruccin rand(n) y rand(m,n) crean matrices nxn y mxn
respectivamente con entradas aleatoriamente generadas distribuidas entre 0 y 1.
Magic(n) crea una matriz cuadrada de orden n mgica (las filas y las columnas
suman la misma cantidad) con entradas enteras.
a) rand(4)
b) rand(2,3)
c) magic(3)
3-

d) eye(6)
e) eye(4,2)
f) ones(6)

g) ones(4,7)
h) zeros(4,5)
i) triu(C)

j) tril(A)
k) A(:)

Genere los siguientes vectores sin introducir explcitamente sus elementos.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 9 _

Trabajos Prcticos de lgebra Lineal con MATLAB

a) Un
vector
columna
cuyos
elementos
consecutivos entre el 15 y el 28
b) Un vector fila cuyos elementos sean:
1.5, 2.5, 3.5, 4.5, 5.5, 6.5, 7.5, 8.5 y 9.5

sean

nmeros

naturales

c) Un vector columna tal que su primer elemento sea 3, su ltimo elemento 2 y los
restantes se obtienen incrementando 0.5 al anterior.
4 Dadas las siguientes matrices:
M=

1
2
3

-3
1
-4

N=

-3
1
7

2
5
6

O=

3
1

-3
4

Investigue que operaciones son las siguientes:

a) M
b) M + N
c) M - N
d) O + M
e) M O
f) M N
g) M .N
h) M./N
i) 2 O
j) O^2
k) O.^2
l) 2/M
ll) 2./M
m) M/2
n) M./2
Siempre es posible realizar dichas operaciones? Qu ocurre cuando
posible?

no es

5 - Utilice MATLAB a fin de verificar los resultados que obtuvo al operar con matrices en
la gua de trabajos prcticos correspondiente
6 - Resuelva usando MATLAB la siguiente situacin:
En un bar se preparan 5 tipos de bebidas cuyas recetas se muestran a continuacin:

Bebida 1 (Cuba libre)

1 medida de ron
1 refresco de cola

Bebida 2 (Cuba campechana)

1 medida de ron
refresco de cola
agua mineral

Bebida 3 (Cuba doble)

2 medidas de ron
refresco de cola
de agua mineral

Bebida 4 (Shandi)

cerveza
limonada

Bebida 3 (Especialidad de la casa)

1 medida de ron
cerveza
limonada

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 10 _

Trabajos Prcticos de lgebra Lineal con MATLAB

a) Represente los datos en una matriz cuyas filas sean ingredientes en las
bebidas y cuyas columnas sean tipos de bebidas .
b) La experiencia demuestra que de lunes a jueves se consumen, en promedio
por da, 20 bebidas del tipo 1, 30 del tipo 2, 10 del tipo 3, 20 del tipo 4 y 30 del
tipo 5, mientras que tanto el viernes como el sbado se consumen 30 bebidas
del tipo 1, 40 del tipo 2, 20 del tipo 3, 40 del tipo 4 y 40 del tipo 5. Represente
estos datos en una matriz cuyas filas sean tipos de bebidas y cuyas
columnas sean das de la semana .
c) Use las matrices de los incisos anteriores para obtener una matriz en la que se
pueda leer la cantidad de ingredientes consumidos por da. Responda
entonces a las siguientes preguntas:
Cuntos refrescos de cola se consumen por semana?
Cuntas limonadas?
Cuntas botellas de ron, sabiendo que cada botella rinde 20 medidas?

Submatrices. Matrices por bloques y matrices aumentadas.

7 - Introduzca las siguientes matrices:
-1.3
A = -7
6

0 5 0
4 5
8 -1 3

B=

1
-3
-5

En cada uno de los siguientes casos explique lo que obtiene.

i) M = A(2,3)
iv) P = A(:,[1 3])
vii) T = [Q;ones(2)]

ii) N = A(3,:)
v) Q = A([2 3],[1 2])
viii) S = [P B;eye(3)]

iii) O = A(:,3)
vi) R= [A B]

8 - Introduzca las siguientes matrices:

A=

a)
b)
c)
d)

2 2 3 4 5
-6 -1 2 0 7
1 2 -1 3 4

B=

-1
2
5

Forme la matriz ampliada C = [A B]

Forme la matriz U, que consiste solo en la 2 y 4 columnas de A.
Forme la matriz T, que consiste solo en la 1 y 2 filas de A.
Forme la matriz V, que consiste solo en las filas 2 y 3 y las columnas 1 y 2 de A

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 11 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Trabajo Prctico N2
Determinante. Inversa de una matriz. Sistemas de Ecuaciones
Para calcular el determinante de una matriz MATLAB cuenta con la funcin det
det(A) calcula el determinante de la matriz cuadrada A.
1 - Aplique la funcin det a cada una de las siguientes matrices y explique la respuesta
de MATLAB en cada caso.

27
2
A=
17

0.32 4 0.5
10 1 3
0.23 5 6
B= 5 0 2
0.32 4 0
5 1 5

3 0 0

-1/2 5 -1
1+ 5i 3i
2 -1 1
C= 3 0 -1 D=
E =

3 2 + 4i
0 12 - 9

2 4 - 3

2 - Escriba help inv a fin de conocer lo que realiza esta funcin y ver su sintaxis. Para
cada una de las matrices del ejercicio anterior aplique inv y explique que obtiene
3 - Por anlisis de ejemplos estudie la posible validez de las siguientes propiedades:
a) Inv(I) = I
b) Si A es inversible, inv(inv(A)) = A
c) Si A y B son inversibles, AB es inversible e inv(AB) = inv(B)inv(A)

4 - Es posible encriptar un mensaje secreto usando una

matriz cuadrada inversible

cuyos elementos y los de su inversa son nmeros enteros, sta se llama matriz de
cdigo.
El procedimiento es el siguiente: se asigna a cada letra del abecedario un
nmero (por ejemplo, A=1, B=2, etc., y espacio = 28), se transforma el mensaje
asignando a cada letra el

nmero correspondiente,

los nmeros obtenidos se

ubican de izquierda a derecha en las filas de una matriz M (matriz de mensaje) de tal
modo que el nmero de elementos de cada fila sea igual al orden de la matriz de
cdigo, se multiplica M por la matriz de cdigo por derecha, se transcribe el mensaje
a una cadena de nmeros (que se lee de izquierda a derecha a lo largo de cada
fila) y se enva el mensaje.
El receptor del mensaje conoce la matriz de cdigo, arregla el mensaje encriptado
en una matriz ubicando los nmeros de izquierda a derecha en cada fila, de modo
tal que el nmero de columnas de dicha matriz coincida con el orden de la matriz de
cdigo, multiplica por derecha por la inversa de la matriz de cdigo y puede leer el
mensaje decodificado.
Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 12 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Actividad
Si C es la matriz de cdigo, decodifique el mensaje que se encuentra ms abajo:

C=

-3

-2

-5

-8 -9

-2

12

-6

11

Mensaje:
128, 107, 167, 55, -235, -206, -279, -39, 232, 188, 347, 183, 273, 194, 412, 274, 272, 221,
353, 126.

*****
MATLAB cuenta con dos funciones matriciales que resultan de gran utilidad a la hora
de resolver sistemas de ecuaciones lineales, ellas son rank y rref:
rank(A) proporciona el rango de la matriz A.
rref(A) devuelve la matriz escaln reducida por filas de A la cual nos permite
obtener un sistema de ecuaciones lineales ms simple equivalente al dado.
5 - Para cada uno de los sistemas siguientes, llame A a la matriz de coeficientes y B a la
matriz columna de trminos independientes. Forme la matriz ampliada y use la
funcin rref para encontrar la forma escaln reducida por filas. Muestre que cada
uno de estos sistemas tiene solucin nica y que la solucin est contenida en la
ltima columna de la forma reducida de la matriz ampliada. Use la notacin de
submatrices para asignar a la variable X la solucin.
x + 3y 3z = 1
4x

z = -1
5x 1/3y + 2/5z = -1
7/2x 3y + 2z = -2

x iy z = i
y iz = i
ix + 2y
=1+i

Calcule AX y compare con B

6 - Para cada uno de los siguientes sistemas introduzca la matriz ampliada y use la
funcin rref para encontrar la forma escaln reducida por filas. Concluya que
ninguno tiene solucin.
2x 3y = -2
2x + y = 1
3x + 2y = 1

x 2y + z 4w = 1
x + 3y + 7z + 2w = 2
x 12y 11z 16w = 5

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 13 _

Trabajos Prcticos de lgebra Lineal con MATLAB

7 - Las matrices siguientes son matrices ampliadas de sistemas de ecuaciones que

tienen ms de una solucin.

9 27 3 3 12
9 27 10 1 19

1 3 5 9 6

i
2
2 3 + i
0
2
1+ i
i

2 1 + i 1 2i 2 + i

para cada matriz use la funcin rref. Para obtener el conjunto solucin de estos
sistemas necesitar papel y lpiz.
Para cada sistema asigne la variable A a la matriz de coeficientes, B a la matriz de
trminos independientes, d una solucin particular del sistema y asgnele el nombre
X. Calcule AX y compare con B.
8 - Analice si los siguientes sistemas homogneos son determinados o indeterminados.
Resuelva el ejercicio empleando
a) la funcin rref
b) la funcin rank
c) Indique cual de las dos formas es la ms conveniente en este caso y explique los
alcances y limitaciones de ambas funciones.
x + 2y z + 3w = 0
2x + 4y 2z + 6w = 0
3x + 6y 3z + 9w = 0
x + 3y + z + 2w = 0

-4x + 2y
2w = 0
2x
3z + 2w = 0
x + 3y 4z + 3w = 0
-x
+ 4z 4w = 0

9 - Se pueden resolver sistemas cuadrados AX = B en los cuales la matriz de coeficientes

es inversible realizando X=inv(A)B (Teorema de Cramer).
Confeccione un archivo-M de funcin que resuelva, cuando sea posible, los
siguientes sistemas cuadrados.

1 1 3 1
2 1 1 1

a) A =
1 2 4 2

2 1
3 1

6
5
B=
6

8

b)

2 1+ i 3
A1 = 1 i 0 2 + 3i
0
i
2i

1
B1 = 2
1 + i

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 14 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Trabajo Prctico N3
Espacios vectoriales. Combinaciones lineales. Subespacio generado.
Dependencia e independencia lineal. Base y dimensin
Los vectores de R n se introducen en MATLAB como matrices de una fila o una columna y
se puede operar con ellos de igual modo que con las matrices.
Es posible graficar vectores del plano empleando la funcin compass
compass(x,y) dibuja una flecha que emana del origen de coordenadas cartesianas y
cuyo extremo es el punto del plano (x,y).

Ejemplo:
compass(2,3)
Produce el siguiente grfico:
Se puede agregar un tercer argumento a la funcin compass, encerrado entre comillas,
que indica el color o tipo de lnea de la flecha, as por ejemplo compass(2,3,r) realizar
el mismo grfico anterior pero el color de la flecha ser rojo.
El siguiente es un cuadro que muestra los caracteres que deben ir encerrados entre
comillas para graficar con distintos tipos de colores y puntos o lneas.

Colores
Y
m
c
r
g
b
w
k

Tipos de punto o lnea

amarillo
violeta
celeste
rojo
verde
azul
blanco
negro

.
o
x
+

-.
-:

punto
crculo
cruz
ms
lnea llena
punto y trazo
lnea de trazos
lnea de puntos

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 15 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Algunas funciones importantes de MATLAB necesarias cuando se quieren realizar grficas:

La instruccin grid realiza un cuadriculado del grfico actual.
Pueden ponerse ttulos, comentarios en los ejes o en cualquier otra parte con los
siguientes comandos que tienen una cadena como argumento:
title
xlabel
ylabel
gtext
text

ttulo del grfico. Ej: title(Vectores)

comentario en el eje x. Ej: xlabel(Eje x)
comentario en el eje y. Ej: ylabel(Eje y)
texto posicionado con el mouse. Ej: gtext(v=(2,3))
texto posicionado mediante coordenadas Ej: text(2,3, v=(2,3))

hold se emplea cuando se desea ver simultneamente dos o ms grficos en una

misma figura, hold on congela la pantalla grfica actual de forma que las grficas
posteriores se superponen en ella. Escribiendo hold off se libera el hold.
1 - Ingrese el siguiente grupo de sentencias y observe el grfico obtenido
compass(2,3)
hold on
compass(1,-1)
compass(1/2,-3)
hold off
2 - Dados los vectores u = (2,4) y v = (-3,2)
a) Resuelva con MATLAB las siguientes operaciones:
u + v, u v,
2u + 5v, 1/2v
b) Grafique en una misma figura, los vectores dados y los obtenidos de resolver las
operaciones anteriores, de tal manera que se diferencien entre si ya sea por el color
o tipo de lnea. Ponga ttulos y nombre a los vectores.
3 - Escriba al vector s como combinacin lineal de los vectores s1 y s2. Utilice la funcin
rref para determinar los escalares.
3
2
1
s = -5
s1 = -1
s2 = 3
-3
1
5
4 - Determine en cada caso si el vector v pertenece al subespacio generado por el
conjunto dado. Resuelva el problema de dos maneras diferentes: empleando la
funcin rref y luego la funcin rank. Explique cual es la ms conveniente para este
caso en particular.
a)
3
1
0
1
v = -2
A=
2
,
-1
,
1
1
0
3
0
0
-1
2
0

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 16 _

Trabajos Prcticos de lgebra Lineal con MATLAB

b) v =

0 -3 -1

A=

1 -1 0

2 1 1 ,

1 2 1

5 - Determine si los vectores u, v y w generan al espacio vectorial V.

a) u =

-12 0 -3

v=

b) u =

1 0 -5 3

5 -4 15

v=

w = 1 -4

0 -1 4 2

w=

14

V = R1x3

1 -1 -1 5

V = R1x4

6 - Determine en cada caso si el conjunto de vectores {v1, v2, v3} es linealmente

independiente o linealmente dependiente, considerando los respectivos espacios
sobre el cuerpo de los nmeros reales.
a) v1=

b)

1 -1 2

1
v1 = 3
3

v2 =

v2 =

1 1 -2

-1
1
3

v3 =

v3 =

-2 3 1

-5
-7
3

c)
2
0
0.7
0

v1 =

v2 =

5
0

v3 =

5
-0.8
7
0.3

7 - Determine si los siguientes conjuntos constituyen una base del espacio vectorial
dado:
a)

[ 1 0 1 0 ] , [2 -1 1 0 ] , [0 1 1 1] , [3 0 -1 1]

b)

1
0
1

1
-1
-1

-1
0
1

V =R1x4

V = R3x1

8 - Dada la base B={ (1,2,3), (-1,0,7), (7,5,0)} de R3

Halle las coordenadas de v1, v2 y v3 respecto de la base dada. Realice el clculo
para los tres vectores en un solo procedimiento.
v1 = (1,2,3)

v2 = (0,-3,12)

v3 = (-7, 2,-8)

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 17 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Trabajo Prctico N4
Producto Interior. Norma. ngulo.
Conjuntos Ortogonales y Ortonormales

Si x e y son dos vectores fila o columna de la misma dimensin, dot(x,y) calcula el

producto escalar de los vectores x e y
1 - Halle el producto escalar de u y v en cada caso
a) u = (2,-3,0) v= (0,-1,8)
b) u = (0.2,-5,67,4) v= (1/9,2.75,0,-5)
c) u= (0.73,-1) v=(-2,0)

2 - La norma de un vector se puede calcular de diferentes maneras:

Proponga un vector v de cuatro componentes y calcule:
a)
b)
c)
d)

dot(v,v)^(1/2) Que operaciones realiza?

sqrt(dot(v,v)) Que operaciones realiza?
norm(v,2)
Compare los resultados y repita el ejercicio para vectores de ms componentes.

3 - Dados v= [ 1 -1 ] y w = [2 -3]
Calcule:
a) El versor de v
b) Un vector paralelo a v de mdulo 5
c) Un vector paralelo a w de mdulo 3
d) La distancia de v a w
e) La proyeccin de v sobre w. Utilizando el valor encontrado halle un vector
ortogonal a w
Realice una grfica de los vectores v y w y de cada uno de los vectores resultantes
de cada tem, todos en la misma figura para poder comparar.

4 - Calcule la medida del ngulo entre v y w para los vectores del ejercicio anterior.
(Ayuda: MATLAB cuenta con la funcin acos, tal que acos(x) da la medida del
ngulo cuyo coseno es x expresado en radianes).

5 - MATLAB cuenta en su librera con la funcin matricial trace que calcula la traza de
una matriz cuadrada. Para los siguientes pares de matrices calcule el producto
interior definido como A . B = Tr (Bt.A)

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 18 _

Trabajos Prcticos de lgebra Lineal con MATLAB

6 5
a) A =

0 1

3 5
B=

1 4

b)

3 1 1/ 2
A = 7
42 21
4 1 / 2 2 / 5

12 5 4
B = 1 / 7 6 20
2
4 10

6 - Calcule la norma de cada una de las matrices del ejercicio anterior.

7 - Determine cuales de los siguientes conjuntos son ortogonales y cuales ortonormales:
a) { (3 , -1, 4, 0), (1, 0, -3/4, 5), (-5, -15, 0, 1) }
b)

c)

-3
0
-5

0
1 ,
2
-3/4

2
-1
3
4/5

[-1/30 5/30 -2/30] , [2/5

0 -1/5 ] , [1/6 2/6 2/6]

8 - Para obtener bases ortonormales de un espacio vectorial, MATLAB cuenta con la

funcin orth.
orth(A) da una base ortonormal del subespacio generado por las columnas de la
matriz A.
Realice lo siguiente:
a) Introduzca la matriz A=[1 2;0 4]
b) Escriba B=orth(A)
c) Verifique que las columnas de B generan al subespacio columna de A y que son
linealmente independientes. Es decir, que forman una base.
d) Realice el producto interior de los vectores columnas de B y calcule sus normas a
fin de verificar que dicha base es ortonormal.
9 - Determine una base ortonormal para los siguientes subespacios:
a) El subespacio generado por v1= ( -2,3,-1), v2=(0,-2,2).
b) El subespacio generado por v1 =(-2,5,0,4), v2 = (1,1,2,-1), v3 = (1,0,0,0)
c) W = { (x,y,z) R3 / x + 2y z = 0 }

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 19 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Trabajo Prctico N 5
Recta y Plano

 Grficos con MATLAB en el plano

Introduzca los vectores fila u=[3 0 5 4], v= [0 0 1 2] y a continuacin escriba plot(u,v),
oprima enter y observe que ocurre.
En este caso el comando plot representa una biyeccin entre los elementos de u y
los elementos del u. En los ejes aparece una porcin de escala necesaria (autoescalado) para que se puedan representar las componentes de los vectores dados.
El comando plot grafica entonces en este caso los puntos (3,0), (0,0), (5,1) y (4,2) y
dibuja un segmento de recta entre los consecutivos como lo muestra la figura:

Proponga vectores con gran diferencia de valores (positivos y negativos) y

observe la acomadacin de las escalas.

 Grfica de funciones
Se desea graficar la funcin seno en el intervalo [-2, 2]. Para ello se escriben las
siguientes instrucciones:
x=-2*pi:0.1:2*pi;
y=sin(x);
plot(x,y)
El vector x es una particin del dominio con incremento de 0.1, mientras que y es un
vector con los valores que toma el seno en los puntos de esa particin.
Este conjunto de instrucciones produce el siguiente grfico:
Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 20 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Pueden tambin realizarse grficas de curvas definidas paramtricamente. Por

ejemplo:
t=0:2*pi/100:2*pi;
x=sin(2*t);
y=sin(3*t);
plot(x,y)
produce el grfico:

Veamos a modo de ejemplo como graficar la elipse cuya ecuacin est dada por:

x2 y2
+
=1
9
4
Esta elipse tiene por centro a C = (0,0), la longitud del semieje mayor es 3, del semieje
menor es 2, su eje focal se encuentra sobre el eje X, sus vrtices son A = (3, 0), A = (-3, 0).
Para realizar la grfica, se dan valores a x y stos deben incluir las absisas de los vrtices
de modo que pueda verse en el grfico la elipse completa y no una porcin de ella, Por
lo tanto se tomarn los elementos del vector x en el intervalo [-3, 3]. Por otra parte, el
vector y debe estar expresado en funcin de x.
Entonces, despejando y se tiene que:

y= 4

4 2
x
9

Por lo tanto las instrucciones siguientes producirn el grfico de la elipse dada.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 21 _

Trabajos Prcticos de lgebra Lineal con MATLAB

x=-3:.1:3;
y1=sqrt(4-4/9*x.^2);
y2=sqrt(4-4/9*x.^2);
plot(x,y1)
hold on
plot(x,y2)
axis equal

2
1.5
1
0.5
0
-0.5
-1
-1.5
-2
-3

-2

-1

De manera similar se puede graficar cualquier cnica

Ejercicio:
Realice el grfico de las siguientes rectas de R2 de modo tal que en cada eje
aparezcan escritos sus nombres X e Y, el ttulo del grfico sea Recta, sobre la recta
se vea escrita su ecuacin y la grfica est cuadriculada. De la observacin de la figura
identifique algunos puntos pertenecientes a cada recta. (Nota: tenga en cuenta que
para graficar una recta son necesarios solo dos puntos).
a) y = 2x 3

b) y = 4x

c)

x=t
y= 3 t

d) (x,y) = t(1,2)

(Tenga en cuenta en el punto (d) que la ecuacin vectorial debe expresarse en forma
paramtrica o cartesiana)

Existen dos formas en las que se puede realizar una graficacin simultanea de
curvas. Esto es importante cuando se quieren comparar los grficos de dos o ms
funciones:
Ejemplo:
x=[0 3];
y1=2x;
y2=2x-1
plot(x,y1,x,y2)
Produce el grfico:

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 22 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Otra forma es usando el comando hold,

Ejemplo:
x=[0 3];
y1= 2x;
plot(x,y1)
hold on
y2=-3x+4;
plot(x,y2)
y3= (-1/2)x
plot(x,y3)
hold off
Escriba la lista precedente de instrucciones y observe lo que obtiene.
Por defecto los ejes se autoescalan. Si se desea cambiar la escala se usa el
comando axis
Si c=[Xmin Xmax Ymin Ymax] es un vector de cuatro elementos donde Xmin, Xmax son los
valores mnimo y mximo respectivamente que se toman sobre el eje X y Ymin, Ymax, son
los valores mnimo y mximo que se toman sobre el eje Y, axis(c) establece el escalado
de los ejes a los lmites prescritos
axis sin argumento congela el escalado actual para los grficos siguientes y
escribiendo axis nuevamente, se vuelve al autoescalado.
axis(equal) produce separaciones de igual medida en los dos ejes.
Ejercicios:

1 - D las ecuaciones de cuatro rectas paralelas tal que una de ellas contenga al
origen. Grafquelas en una misma figura de forma que:
a) El ttulo de la figura sea Rectas paralelas
b) Los ejes estn identificados con x e y
c) La grfica aparezca cuadriculada.
d) En cada recta aparezca su ecuacin
e) La escala en ambos ejes sea la misma.

2 - Realice la grfica de los siguientes pares de rectas y

determine su punto de
interseccin a partir de la figura. Luego calcule, usando MATLAB, mediante sistemas
de ecuaciones la interseccin de dichas rectas y compare resultados.

a) L1: y = 2x 5
L2: y = 1/3 x 5/3
b) L1: y = -7x
L2: y = x

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 23 _

Trabajos Prcticos de lgebra Lineal con MATLAB

c) L1:

x4= t
y
= -t

L2:

x 1 = -1/2 t
y =-t

3 - Las siguientes rectas son ortogonales:

L2: y = -1/2x
L1: y = 2x 5
a) Realice la grfica de ambas rectas en la misma figura, cuadricule y observe si
stas se ven perpendiculares.
b) Escriba ahora la instruccin axis(equal) y vea que ocurre.
c) Explique la diferencia entre los items a) y b).

4 - Escriba la ecuacin de la recta que contiene al origen y es ortogonal a la recta de

ecuacin:
x 2 = 3t
y = -t
Como en el caso anterior verifique lo obtenido (tenga en cuenta que en este
caso deber graficar dos rectas).

 Grficos con MATLAB en el espacio

La instruccin mesh crea en un grfico tridimensional una superficie de malla que
est definida por las terceras componentes de los puntos sobre una cuadrcula
rectangular en el plano XY.
Para realizar el grfico de una funcin se debe definir en primer lugar los vectores x e y
que dan la particin del rectngulo (es decir se define el dominio) con la funcin
meshgrid.
Ejemplo 1:
[x y]=meshgrid(-4:0.2:4);

% genera una rejilla rectangular de puntos

% igualmente espaciados entre -4 y 4 en cada eje.

z=y.^2/16-x.^2/9;
mesh(x,y,z)
Este conjunto de instrucciones genera la grfica del paraboloide hiperblico

z=

y 2 x2

sobre el rectngulo [ -4,4] x [-4,4] que se muestra en la figura siguiente:

16
9

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 24 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Un grfico de superficie de la misma funcin se obtiene con la funcin surf en lugar de

mesh. Por lo que:
surf(x,y,z)
producir el grfico siguiente:

Y si se le agregan las sentencias:

shading interp
colormap(copper)
se produce:

Del igual modo se pueden graficar otras cudricas.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 25 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Ejercicio: Escriba las siguientes instrucciones en MATLAB

[x y] = meshgrid(-2:0.1:2);
z= -x + y 1;
surf(x,y,z)
Qu figura obtiene?
Nota: Los comandos axis, hold, grid, title, xlabel, ylabel, zlabel, etc. tienen en este caso la
misma aplicacin.
Ejercicios:
5 - Realice el grfico de los siguientes planos, de tal modo que en cada eje est escrito
su nombre( X, Y y Z) y el grfico est cuadriculado.
a) 2x +3y z =2
b) x + y - z =0
c) x z =5
6 - Se desea realizar la grfica del plano x = 3
Tenga en cuenta que la informacin que nos brinda esta ecuacin es que mientras y
y z varan en todos los reales x toma siempre el valor 3. Luego el conjunto de
instrucciones que se deben escribir varan ligeramente de las dadas anteriormente.
En este caso se definirn primero los valores para y y z ya que stas no dependen de
ninguna variable, y se escribir:
[y z]=meshgrid(-2:0.2:3);
x=0y+0z+3;
mesh(x,y,z)
Por qu se escribe x=0y+0z+3 y no x=3 simplemente? Haga la prueba usando esta
instruccin y vea que ocurre.
7 - Realice la grfica de los siguientes planos:
a) El plano paralelo al plano XY que contenga al punto (0,1,3)
b) El plano paralelo al eje X que contenga al punto (-1,1,5)
c) El plano YZ
d) Un plano que contenga al eje X

8 - En cada uno de los siguientes casos grafique los planos dados en una misma figura
(use hold on y hold off) y observe su posicin relativa. Resuelva tambin en cada

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 26 _

Trabajos Prcticos de lgebra Lineal con MATLAB

caso, usando MATLAB, el sistema formado por las ecuaciones de dichos planos.
Compare la posicin relativa con la compatibilidad o incompatibilidad del sistema,
con el hecho de ser determinado o indeterminado y con el rango de la matriz de
coeficientes. Extraiga conclusiones.
a)

1 :2x + 3y 2z = 0

2: 5x y + 4 z = 0

b)

1 :x + y z + 10 = 0

2 :2x + 2y 2z = 0

c)

1: x + 3y 2z = 0

2: -x + 5y + 3z 1=0

3: 8y + z 1 = 0

d)

1: 5x 2y + z 3 = 0

2: -x + z = 0

3: 2x 2z 20 = 0

Realice los grficos sobre el rectngulo [-4,4] x [-4.4]

 Grfica de rectas en el espacio

Como pudo observar, con el comando mesh se realiza la grfica de una superficie
de malla en el espacio tridimensional, de modo que podemos, mediante este comando,
ver la recta pero solo como la interseccin de dos planos.
Para graficar una recta en el espacio se emplea el comando plot3, que es una
extensin del comando plot a grficos en el espacio tridimensional.
Ejemplo:
Se desea graficar la recta que contiene al punto (1,0,3) y es paralela al vector (-1,1,2).
La ecuacin vectorial de esta recta es (x,y,z) = (1,0,3) + t(-1,1,2)
Y las ecuaciones paramtricas:
x=1t
y=t
z = 3 + 2t
Se puede realizar la grfica empleando las ecuaciones paramtricas del siguiente modo:
t = -3:0.1:3;
x =1-t; y = t; z=3+2t;
plot3(x,y,z)

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 27 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Si se agrega la instruccin grid puede verse el grfico cuadriculado

Puede realizarse tambin la grfica empleando las ecuaciones cartesianas de la recta,

en este caso son :
z-3
-x + 1 = y =
2
Es importante tener en cuenta que para graficar una recta con plot3 a partir de sus
ecuaciones cartesianas, stas deben estar expresadas de tal forma que dos de las
variables, estn en funcin de una misma (por ejemplo y y z en funcin de x). En este
caso las ecuaciones cartesianas quedaran expresadas:
y = -x + 1
z = -2x + 5
Escriba hold on para graficar en la misma figura y as poder verificar si efectivamente se
obtiene la misma recta, luego escriba las siguientes instrucciones:
x =-3:0.1:3;
y=-x+1;
z=-2x+5;
plot3(x,y,z,b)
hold off
El argumento b en plot3 significa que la recta se graficar en color azul.

9 - Realice el grfico de las siguientes rectas:

a)
b)
c)
d)

La recta que contiene al origen y es paralela al vector (-1,2,6)

La recta que contiene al punto (-2,0,2) y es paralela al eje X.
La recta que contiene al punto (1,1,1) y es ortogonal al plano XY.
La recta cuyas ecuaciones cartesianas son:
2x + 3y + z = 0
3x + 5y 2z +3 = 0
(Sugerencia: pase primero estas ecuaciones a la forma paramtrica)
Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 28 _

Trabajos Prcticos de lgebra Lineal con MATLAB

10 - En cada uno de los siguientes casos realice los grficos de la recta y el plano dados,
en una misma figura y analice su posicin relativa. Luego resuelva el sistema
formado por las ecuaciones de la recta y el plano, empleando MATLAB, y extraiga
conclusiones (Sugerencia: en este caso es aconsejable que trabaje con las
ecuaciones cartesianas de la recta)
x5
a) L:

y5

z1

=
2

: 2x 3y + 5z = 0

=
-3

a) L: (x,y,z) =(1,-1,0) + t(1,1,0)

: z = 3
b)

x1=t
y
= -t
z 2 = 5t

: 5x z 3 = 0

Trabajo Prctico N 6
Transformaciones lineales.
1 - Halle la matriz asociada a la transformacin lineal dada respecto de las bases B y B.
Tenga en cuenta que primero deber trabajar con papel y lpiz hasta formar
los
sistemas de ecuaciones, que tendrn la misma matriz de coeficientes pero variarn
en sus trminos independientes.
f: R3 R2 tal que f(x,y,z) = (x,z)
B=

(-1,2,0), (1,0,3), (0,1,3)

B =

(1,2), (-1,1)

Observe la matriz ampliada que ha formado con las tres columnas de trminos
independientes. Qu vectores estn ubicados en las columnas de la matriz de
coeficientes? Qu vectores estn ubicados en las columnas de trminos
independientes? Luego de aplicar rref a la matriz ampliada qu matriz representa la
submatriz formada por las tres ltimas columnas?
De estas observaciones, confeccione una regla que le permita en lo sucesivo
trabajar menos con papel y lpiz.

1 - En cada uno de los siguientes casos halle la matriz asociada a la transformacin

lineal dada respecto de las bases B y B.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 29 _

Trabajos Prcticos de lgebra Lineal con MATLAB

a) f: R2 R3 tal que f(x , y) = (x, x+y, 2y)

B=

(2,3), (-1,5)

B =

(1,2,3), (-3,-3,-3), (2,0,1)

b) f : R3 R3 tal que
f( -2,3,1) = (0,1,2) , f(1,1,1) = (-3, 5,2) ,

f(0, 0, 1) = (1,0,1)

B=

(1,2,0), (0,1,2),(0,0,3)

(-2,3,1), (1,1,1), (0,0,1)

B =

2 - Construya un archivo de funcin tal que, dada la matriz asociada a una

transformacin lineal respecto de las bases connicas, calcule la imagen de un
vector.

3 - Sea A la matriz asociada a una transformacin lineal f: R2 R3 respecto de las bases

cannicas en ambos espacios.
1 0
A = -2 1
6 8
Use el archivo creado en el ejercicio anterior a fin de calcular la imagen de los
siguientes vectores:
a) f(0,-1)
b) f(3,-4)
c) f((0,0)

 Modificacin de figuras del plano por medio de transformaciones lineales.

Si se escriben siguientes instrucciones de MATLAB se obtiene el grfico de un
tringulo:
x=[2 4 6 2];y=[5 7 2 5];
plot (x,y)
grid

Recuerde que plot grafica los puntos de una biyeccin entre x e y, en este caso (2,5),
(4,7), (6,2) y (2,5), y los une con un segmento de recta. Como puede observarse en la

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 30 _

Trabajos Prcticos de lgebra Lineal con MATLAB

grfica estos puntos son justamente los vrtices de un tringulo. Por lo tanto si se desea
dibujar un tringulo conociendo sus vrtices se forman los vectores x e y de tal manera
que los elementos de x sean las absisas y los elementos de y sean las ordenadas de cada
punto. Observe que el primer elemento debe repetirse al final tanto en x como en y,
porque si esto no ocurriera plot no unira el ltimo vrtice con el primero. Pruebe
realizando la grfica para los vectores x e y dados pero sin el ltimo elemento.

Ejercicio:
Dibuje las siguientes figuras, respetando el orden en el que estn dados los vrtices:
a) El tringulo de vrtices (-1,3), (0,0), (3,-2)
b) El paralelogramo de vrtices (1,-2), (2,1), (6,1), (5,-2)
c) El polgono de vrtices (-3,0), (-1,2), (2,2), (4,-2), (0,-5)
d) Modifique el orden dado a los vrtices en el tem c) y observe la figura que
obtiene. Extraiga una conclusin.

Dada una figura en el plano se desea ver que efectos produce sobre ella
determinada transformacin lineal:

Ejemplo:
La transformacin lineal f: R2 R2 tal que f(x, y) = (3x, y) produce una expansin a
lo largo del eje X.
La matriz asociada a esta transformacin lineal respecto de la base cannica de R2
es:
3 0
A =
0 1
Ahora bien, se desea aplicar esta transformacin lineal al tringulo de vrtices
P=(0,0), Q = (2,5) y R =(6,3).
El procedimiento a seguir es el siguiente:
1 - Se introducen los vectores P, Q y R. Es conveniente que stos sean
expresados como vectores columnas, ya que posteriormente debern multiplicarse
por la matriz A.
2- Se construye una matriz T cuyas columnas son P, Q y R. Esta matriz tiene
como primera fila a las absisas de los vrtices y como segunda fila a las ordenadas.
3 - Se forman los vectores x de las absisas e y de las ordenadas de los vrtices
y se realiza la grfica del tringulo.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 31 _

Trabajos Prcticos de lgebra Lineal con MATLAB

4- Se realiza el producto AT a fin de obtener las imgenes.

5-Se construyen los vectores x1 de las absisas e y1 de las ordenadas de las
imgenes y se realiza la grfica, en una misma figura, del transformado del tringulo.
Las instrucciones correspondientes son:
1

P=[0 0];Q=[2 5];R=[6 3];

T=[P Q R P]

x=T(1, : ); y=T(2, : ) ;
plot(x,y)

A=[3 0
0 1] ;
Im=AT;

hold on
x1=Im(1,: ); y1=Im(2,: )
plot(x1,y1,g)

Realice a continuacin para el mismo tringulo y siempre en la misma figura del

MATLAB una expansin a lo largo del eje Y con un factor igual a 2.
La matriz asociada a esta transformacin lineal respecto de la base cannica es:
1
0
M=
0
2

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 32 _

Trabajos Prcticos de lgebra Lineal con MATLAB

De la misma forma en que se aplica una transformacin lineal a una figura del plano,
puede tambin aplicarse a curvas planas, el procedimiento es similar solo que se trabaja
con una mayor cantidad de puntos.
Suponga que se desea aplicar una rotacin un ngulo =

a la porcin de la parbola

y = 4 x 2 + 2 correspondiente al intervalo [-2,2] y teniendo en cuenta que la matriz de

cos
rotacin es

sen

sen
cos

, el siguiente conjunto de instrucciones producirn lo

deseado.
x=-2:.2:2;
y=4*x.^2+2;
plot(x,y)
hold on

dibuja la parbola original

m=[cos(pi/4) sin(pi/4);-sin(pi/4) cos(pi/4)];

puntos=[x;y];

Introduce la matriz asociada

forma una matriz cuya primera fila est compuesta por las absisas
y la segunda por las ordenadas de los puntos de la parbola

puntos_imag=m*puntos;
x1=puntos_imag(1,:);
y1=puntos_imag(2,:);

Calcula la imagen de dichos puntos

Forma el vector fila de las absisas de los puntos imagen
Forma el vector fila de las ordenadas de los puntos imagen

plot(x1,y1,'r')

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 33 _

Trabajos Prcticos de lgebra Lineal con MATLAB

En el siguiente cuadro figuran las matrices asociadas respecto de la base cannica de

R2 de algunas transformaciones lineales:

Expansin por un factor k, a lo largo del eje X

si k >1
Contraccin por un factor k, a lo largo del eje X si 0 < k < 1

Expansin por un factor k, a lo largo del eje Y si k > 1

Contraccin por un factor k, a lo largo del eje Y si 0 < k < 1

Expansin por un factor k, en ambas direcciones si k > 1

Contraccin por un factor k, en ambas direcciones si 0< k < 1

Reflexin respecto del eje X

0

-1

-1

0
Reflexin respecto del eje Y

cos

-sen
Rotacin con centro en el origen de coordenadas y ngulo

sen

cos

Ejercicios:
4 - Construya un archivo-M para que dibuje el cuadriltero de vrtices (0,1), (2,4), (4,4),
(6,1) en ese orden y luego aplique a ste las siguientes transformaciones lineales, de
modo que todos los grficos se vean en una misma figura, recuerde que para
diferenciarlos puede usar distintos tipos de lneas y/o colores.
a) Expansin por un factor 5 a lo largo del eje Y.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 34 _

Trabajos Prcticos de lgebra Lineal con MATLAB

b) Contraccin por un factor 1/2 en ambas direcciones.

c) Reflexin respecto del eje X

5 - Modifique el archivo anterior para que dibuje la figura de vrtices (2,-2), (2,7), (4, 5),
(2,3) en ese orden y luego aplique a esta las siguientes transformaciones lineales
precediendo del mismo modo que en el ejercicio anterior.
a) Expansin por un factor 4 en ambas direcciones.
b) Reflexin respecto del eje Y.
c) Rotacin en el sentido positivo con centro en el origen y un ngulo de 60 (Tener
en cuenta que MATLAB trabaja con ngulos expresados en radianes por lo que
60 = /3 que en el lenguaje de MATLAB se escribe pi /3).

6 - Dibuje una curva cualquiera y aplique la o las transformaciones lineales que ms le

agraden a fin de experimentar.

 Modificacin de figuras y cuerpos en el espacio por medio de transformaciones

lineales
Al igual que en el plano, se puede ver en el espacio tridimensional, los efectos que
produce una determinada transformacin lineal sobre una figura o un cuerpo.

Para realizar la grfica de figuras o cuerpos en el espacio usaremos el comando

plot3 que funciona de la siguiente manera:
Dados los vectores x=[x1 x2 x3...xn ] , y=[y1 y2 y3...yn] , z=[z1 z2 z3...zn] plot3(x,y,z)
dibuja los puntos (x1,y1,z1), (x2,y2,z2), (x3,y3,z3),..., (xn,yn,zn) y los une con un segmento de
recta. De esto se deduce que el vector x est formado por las primeras componentes de
tales puntos, el vector y por las segundas componentes y el vector z por las terceras.
Si se desea aplicar a estas figuras o cuerpos una transformacin y luego poder
visualizar su efecto se procede en forma anloga a lo realizado en el plano.

Ejercicios:

8 - Dado el tringulo de vrtices (2,3,-1),

(5,0,-2), (4,-2,0), aplique las siguientes

transformaciones lineales de modo que todas puedan verse en una misma figura,
diferencindose entre si por distintos colores o tipos de lneas.

a) Simetra respecto del plano XY definida por: f: R3 R3 / f(x,y,z) = (x,y,-z)

b) Simetra respecto del origen de coordenadas definida por

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 35 _

Trabajos Prcticos de lgebra Lineal con MATLAB

f: R3 R3 / f(x,y,z) = (-x,-y,-z)
c) Simetra respecto del eje Z definida por

f: R3 R3 / f(x,y,z) = (-x,-y,z)

9 - Dada la figura de vrtices (2,0,0), (1,5,3), (0,-3,4) y (-6,5,3) en ese orden. Halle su
imagen respecto de la transformacin lineal f y a la figura resultante aplique la
transformacin lineal g. Siendo
f: R3 R3 / f(x,y,z) = (2x,2y,2z) y
g: R3 R3 / g(x,y,z) = (x+y,y,x+y+z)
Grfico de cuerpos
Se desea dibujar un tetraedro de vrtices P = (2,4,0), Q = (-1,6,0), R = (0,2,0) y S = (0,4,5)
Se realizar el grfico de este cuerpo con el comando plot3.
Observe los siguientes grficos y ejecute las instrucciones para poder comprender el
modo de confeccionar los vectores x, y y z
S
S

R
P

R
Q

Figura 1

Q
Figura 2

Intente unir los vrtices P, Q, R y S de la Figura 1 formando un tetraedro como el de la

Figura 2 pero sin levantar el lpiz (se le permite dibujar una lnea sobre otra).
Escriba el orden que sigui para unir estos puntos (por ejemplo: P, Q, R, P, S, R, Q, S).
Forme los vectores x, y y z con las primeras, segundas y terceras componentes
respectivamente de estos puntos en el orden establecido
Los pasos siguientes para realizar la grfica son ya conocidos.
Ejercicios:

10 - Aplique al tetraedro del ejemplo anterior las siguientes transformaciones lineales de

modo que todas se puedan ver en una misma figura.
a) Una rotacin alrededor del eje Z un ngulo de 60
b) Una rotacin alrededor del eje X un ngulo de 30
c) Una rotacin alrededor del eje Y un ngulo de 45
Las matrices asociadas a estas transformaciones lineales se describen en el
cuadro siguiente.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 36 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Las siguientes son matrices asociadas de algunas transformaciones lineales en el espacio

respecto de la base cannica de R3

1
0
0

0
cos
sen

0
- sen
cos

Rotacin alrededor del eje X un ngulo

cos
0
sen

0
1
0

-sen
0
cos

Rotacin alrededor del eje Y un ngulo

cos
sen
0

-sen
cos
0

0
0
1

Rotacin alrededor del eje Z un ngulo

k
0
0

0
k
0

0
0
k

1
0
0

0
1
0

0
0
0

Proyeccin sobre el plano XY

0
0
0

0
1
0

0
0
1

Proyeccin sobre el plano YZ

1
0
0

0
-1
0

0
0
-1

Reflexin respecto del eje X

1
0
0

0
-1
0

0
0
1

Reflexin respecto del plano XZ

Expansin por un factor k, si k > 1

Contraccin por un factor k , si 0 < k < 1

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 37 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Trabajo Prctico N 7
Valores y vectores propios. Diagonalizacin.
Formas cuadrticas

Las funciones matriciales que se describen a continuacin le sern tiles para

resolver los problemas planteados en este Trabajo Prctico:
eig: eig(A) produce un vector columna cuyos elementos son los valores propios
(eigenvalores) de una matriz cuadrada A.
[V D]=eig(A) produce una matriz V cuyas columnas son los vectores propios de A y una
matriz D diagonal en la cual los elementos de la diagonal son los valores propios de A
poly: si A es una matriz cuadrada de orden n , poly(A) es un vector fila, tal que sus n + 1
elementos son los coeficientes del polinomio caracterstico de A ordenados en forma
decreciente.

1 - Halle los valores propios de las siguientes matrices consideradas sobre el cuerpo de
los nmeros reales.
a)

c)

-3
-3

10
-8

1 0
0
3 10 15
-2 -6 -9

b)

d)

2 -2 3
10 -4 5
5 -4 6

2
0
0
0

-1/2
-2
-1
-2

0 1
0 1
3 2
0 1

2 - Para cada una de las matrices del ejercicio anterior encuentre su polinomio
caracterstico y verifique que los valores propios obtenidos son efectivamente las
races de dicho polinomio.

3 - Verifique la posible validez de las siguientes proposiciones mediante el anlisis de

ejemplos. En todos los casos A es una matriz cuadrada.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 38 _

Trabajos Prcticos de lgebra Lineal con MATLAB

a) Si A es una matriz diagonal entonces sus valores propios son los elementos de la
diagonal.
b) 0 es un valor propio de A si y solo si A no es inversible
c) A y At tienen el mismo polinomio caracterstico.
d) A es raz de su polinomio caracterstico. (Teorema de Cayley-Hamilton)

4 - En cada uno de los siguientes casos indique si la matriz A dada es diagonalizable. Si

as lo fuera encuentre la matriz C que la diagonaliza y calcule C-1 A C.

a)

-1 -4
0 3
0 -2

-1
1
0

b)

1
-2
2

0
-1
0

0
0
-1

5 - Para cada una de las siguientes matrices simtricas, halle una matriz ortogonal que la
diagonalice:

A=

12

12

1
1
B=
6
8

1
5
0
0

0 0
8 3
3

C=

3
2

1
2

3
2

4
9
2
7
2

9
2

1
2
7
2

6 - En los siguientes casos:

i) Realice la grfica correspondientes a las ecuaciones cuadrticas dadas
ii) Escriba la ecuacin en la forma Av.v = d (donde A es una matriz simtrica ) y
elimine el trmino que contiene a xy rotando los ejes un ngulo .
iii) Escriba la ecuacin en trminos de las nuevas variables e identifique la cnica
que se obtiene.
iv) Verifique los resultados comparndolos con la grfica realizada.
a) 3x2 2xy 5 = 0
b) 3x2 6xy + 5y2 = 36
c) x2 + 2xy y2 = 0

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 39 _

Trabajos Prcticos de lgebra Lineal con MATLAB

Bibliografa
The Math Works, Inc. (1997). MATLAB, Versin 5, Gua del Usuario.
The Math Works, Inc. (1996). MATLAB, Edicin de Estudiante, Versin 4, Gua del Usuario.
(1 Ed).
Nakamura, S. (1997). Anlisis Numrico y Visualizacin Grfica con MATLAB (1 Ed.).
Mexico: Ed. Prentice-Hall Hispanoamericana, S.A.
Nakos, G. y

Joyber, D. (1999). lgebra Lineal con aplicaciones. Mxico: International

Thomson Editores
Poole, David. (2003). lgebra Lineal. Mxico: International Thomson Editores.
Gerber, H. (1992). lgebra Lineal. Mxico: Grupo Editorial Iberomrica.
Noble B., Daniel J. (1989). lgebra Lineal Aplicada (3 Ed.). Mxico: Ed. Prentice Hall.

Fac. de Ciencias Exactas y Tecnologas - UNSE

_ Pgina 40 _