APÉNDICE A

Tutorial de Matlab
Presentamos un tutorial de Matlab, una herramienta potentı́sima, casi estándar para cálculos en
muchas ramas de la Ingenierı́a, y de uso razonablemente simple. Haremos una descripción de los elementos
básicos de Matlab y remitiremos al estudiante interesado en mejorar sus conocimientos del mismo a cualquier
edición del manual de referencia del programa [21]. También es interesante el libro sobre Matlab de Higham
y Higham [16] y una buena referencia en castellano con ejemplos procedentes de problemas de Cálculo
Numérico es el de Quintela [23].

A.1. Conceptos básicos

Para arrancar Matlab, se procede como con cuaquier programa Windows, o sea, Inicio, Programas,
Matlab o Student Matlab caso de que utilicemos la versión educacional. Una vez arrancado aparece el cursor
con el sı́mbolo (>>) o (EDU >>), indicando que se pueden introducir órdenes. De hecho, en este tutorial,
cuando aparezca este sı́mbolo, se tiene que introducir por teclado la orden que aparece escrita a la derecha
del mismo.
La utilización más básica de Matlab es como calculadora 1 . Ası́, por ejemplo, para calcular cos(5) · 27.3 ,
se debe introducir 2 :

>>cos(5)*2^7.3
ans =
44.7013

Matlab mantiene en memoria el último resultado. Caso de que ese cálculo no se asigne a ninguna variable,
lo hace a una variable por defecto de nombre ans. Si queremos referirnos a ese resultado, lo haremos a través
de la variable ans, y si no se asigna ese nuevo cálculo a ninguna variable, volverá a ser asignado a ans.

>>log(ans)
ans =
3.8000

En este momento cabrı́a preguntarnos si tratamos con un logaritmo decimal o con uno neperiano (natural).
Para saberlo, pedimos ayuda acerca del comando log utilizando:

>>help log
LOG Natural logarithm.
LOG(X) is the natural logarithm of the elements of X.
Complex results are produced if X is not positive.

See also LOG2, LOG10, EXP, LOGM.

1 Funcionando de este modo, es similar a una calculadora programable, aunque bastante más versátil.
2 Los argumentos de las funciones trigonométricas siempre están en radianes.

381
382 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

Aunque en la explicación de las órdenes, los comandos aparezcan en mayúsculas, se deben usar en
minúsculas.
Por defecto, los resultados aparecen con 4 cifras decimales. Si se necesitara más precisión en los resultados,
se puede utilizar la orden format long repitiendo los cálculos:

>>format long

Para recuperar una orden y ejecutarla otra vez o modificarla se usan la flechas arriba y abajo del cursor
⇑, ⇓. Presionemos ⇑ hasta recuperar la orden:

>>cos(5)*2^7.3
ans =
44.70132670851334

Ejercicio A.1.1 Realizar la siguiente operación: 2.72.1 + log10 108.2.

2.1
Ejercicio A.1.2 Realizar la siguiente operación: e2.7 +log10 108.2
.

Si necesitamos referirnos a determinados cálculos, se asignan a variables y ası́ se pueden recuperar después.
Por ejemplo, podemos recuperar con ⇑ la orden cos(5) · 27.3 y asignar su valor a la variable x editando dicha
orden. Luego podremos utilizarla para otros cálculos.

>>x=cos(5)*2^7.3
x =
44.70132670851334
>>y=log(x)
y =
3.80000318145901

Ejercicio A.1.3 Realizar la siguiente operación: 2.72.1 + log10 108.2 y asignarla a la variable x.
2.1
Ejercicio A.1.4 Realizar la siguiente operación: e2.7 +log10 108.2
y asignarla a la variable t.

Si queremos saber cuánto vale una variable, no tenemos más que escribirla en la lı́nea de comandos y
pulsar Enter.

>>y
y =
3.80000318145901

Como es muy fácil recuperar órdenes previas podemos utilizar esta idea para simular los términos de una
sucesión recurrente. Por ejemplo, xn+1 = cos(xn )

>>x=0.2
x =
0.20000000000000
>>x=cos(x)
x =
0.98006657784124
>>x=cos(x)
x =
0.55696725280964
>>x=cos(x)
x =
0.84886216565827
>>x=cos(x)
x =
 Apéndice A: Tutorial de Matlab 383

0.66083755111662
>>x=cos(x)
x =
0.78947843776687
>>x=cos(x)
x =
0.70421571334199
Ejercicio A.1.5 Repetir la operación anterior hasta que se estabilice el cuarto decimal de x de un paso al
siguiente.
Ejercicio A.1.6 Cambiar el formato para que otra vez se vean sólo cuatro decimales.
Ejercicio A.1.7 Empezando por x = 100 repetir la operación
x2 − 81
x=x−
2x
hasta que se converja en el cuarto decimal. ¿Qué relación hay entre el último x y 81?
Ejercicio A.1.8 Definir A como en vuestro documento de identidad o pasaporte. Empezando por x = 100
repetir la operación
x2 − A
x=x−
2x
hasta que se converja en el cuarto decimal. ¿A qué ha convergido la sucesión? 3
Es interesante comentar de este ejercicio que Matlab distingue entre letras mayúsculas y minúsculas en
los nombres de las variables.
A veces es bueno apagar y encender la calculadora para borrar todo y empezar de nuevo. Esto se hace
con la orden clear. Hay que tener cuidado al utilizarla, ya que borra todas las variables que estén en la
memoria sin pedir confirmación.
>>clear
>>x
??? Undefined function or variable ’x’.
Ejercicio A.1.9 Preguntar el valor de A igual que acabamos de preguntar x. ¿Tiene sentido el resultado?

A.2. Manejo de vectores

Para crear y almacenar en memoria un vector v que tenga como componentes v1 = 0, v2 = 2, v3 = 4,
v4 = 6 y v5 = 8 podemos hacerlo componente a componente:
>>v(1)=0
v =
0
>>v(2)=2
v =
0 2
>>v(3)=4
v =
0 2 4
>>v(4)=6
v =
0 2 4 6
>>v(5)=8
v =
0 2 4 6 8
3 Las calculadoras obtienen la raı́z cuadrada de un número mediante esta sucesión.
384 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

Se puede también definir este vector especificando su primer elemento, un incremento y el último elemen-
to. Matlab rellenará paso a paso sus componentes. Ası́, podemos definir igualmente el vector v como una
secuencia que empieza en 0, avanza de 2 en 2 y que termina en el 8:

>> v = [0:2:8]
v =
0 2 4 6 8
>> v
v =
0 2 4 6 8

Si ponemos ; al final de una lı́nea de comandos, cuando pulsemos la tecla Enter para ejecutarla, se
ejecutará pero no mostrará el resultado en pantalla (se anula el eco en pantalla). Esto es muy útil algunas
veces:

>> v = [0:2:8];
>> v
v =
0 2 4 6 8

Podemos construir el vector v editando directamente entre los corchetes las componentes del vector v:

>>v = [0 2 4 6 8];
>> v
v =
0 2 4 6 8

Es fácil acceder al contenido de una posición del vector, por ejemplo la primera.

>> v(1)
ans =
0

O modificarla:

>> v(1)=-3;
>> v
v =
-3 2 4 6 8

O hacer operaciones entre componentes, v2 · v53 :

>> v(2)*v(5)^3
ans =
1024

Ejercicio A.2.1 Calcular la suma de los elementos de v, elemento a elemento.

Para trasponer un vector o una matriz se usa el apóstrofo, que es el acento que está en la misma tecla
que el signo de interrogación “?”.

>> v’
ans =
-3
2
4
6
8
 Apéndice A: Tutorial de Matlab 385

Como hemos comentado, para recuperar una orden y ejecutarla otra vez o modificarla se usan la flechas
arriba y abajo del cursor ⇑, ⇓. Presionemos ⇑ hasta recuperar la orden:

>> v(1)=-3;

Modifiquémosla para dejar el valor original

>> v(1)=0;

Al definir ese vector v de 5 componentes, en realidad lo que definimos es una matriz fila de cinco columnas,
o sea, un matriz de 1 × 5. Esto se comprueba preguntando el tamaño de v con la sentencia size:

>>size(v)
ans =
1 5

que nos indica que v tiene una fila y 5 columnas.

Ejercicio A.2.2 Definir un nuevo vector que sea el traspuesto de v y aplicar a ese vector el comando size.
¿Es coherente el resultado?

A.3. Introducción al tratamiento de matrices

Haremos una introducción a la definición y manipulación de matrices. Se supone que se ha seguido la
sección anterior y que se dispone de los conocimientos básicos sobre la definición y manipulación de vectores
usando Matlab. La definición de una matriz es muy similar a la de un vector. Para definir una matriz, se
puede hacer dando sus filas separadas por un punto y coma (¡no olvidarse poner los espacios en blanco!):

>> A = [ 1 2 3; 3 4 5; 6 7 8]
A =
1 2 3
3 4 5
6 7 8

o definirla directamente fila a fila, que es más intuitivo:

>> A = [ 1 2 3
3 4 5
6 7 8]
A =
1 2 3
3 4 5
6 7 8

Se puede modificar alguno de los elementos de la matriz A, accediendo a cualquiera de sus posiciones,
por ejemplo:

>> A(2,2)=-9
A =
1 2 3
3 -9 5
6 7 8

Dejemos su valor original:

>> A(2,2)=4;

De igual modo, se la puede considerar como una fila de vectores columna:

386 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

>> B = [ [1 2 3]’ [2 4 7]’ [3 5 8]’]

B =
1 2 3
2 4 5
3 7 8
(Otra vez, es importante colocar los espacios en blanco.)
Ejercicio A.3.1 Sumar los elementos diagonales de la matriz A, refiriéndonos a ellos, elemento a elemento.
Podemos sumar o restar matrices para tener otras matrices.
>>C=A+B
C =
2 4 6
5 8 10
9 14 16
Ejercicio A.3.2 Definir la matriz D = 2B − A.
También podemos multiplicarlas.
>>C=A*B
C =
14 31 37
26 57 69
44 96 117
Ejercicio A.3.3 Definir la matriz D = B − A · B.
Ejercicio A.3.4 Definir la matriz C = AAt .
Podemos definir algunos tipos especiales de matrices, como por ejemplo una matriz de 3 × 3 que tenga
todos sus elementos nulos.
>>I=zeros(3)
I =
0 0 0
0 0 0
0 0 0
Podemos modificar sus elementos diagonales para tener la matriz identidad.
>>I(1,1)=1;
>>I(2,2)=1;
>>I(3,3)=1
I =
1 0 0
0 1 0
0 0 1
Ejercicio A.3.5 Definir la matriz D = B − B · A como D = B(I − A).
Otra forma de definir la matriz identidad es a través de la función diag, que recibe un vector que convierte
en diagonal de una matriz cuyos otros elementos son nulos.
>>J=diag([1 1 1])
J =
1 0 0
0 1 0
0 0 1
Ejercicio A.3.6 Definir una matriz D diagonal cuyos elementos sean −2, 1, 0.2 y −0.7.
Ejercicio A.3.7 Pedir ayuda de la función eye, y definir la matriz diagonal de 10 × 10.
 Apéndice A: Tutorial de Matlab 387

A.3.1. Definición de submatrices

La definición de “subvectores” o submatrices es muy fácil. Si v es
>> v = [0:2:8]
v =
0 2 4 6 8
Podemos definir un vector e cuyas componentes sean las tres primeras componentes del vector v poniendo
>> e=v(1:1:3)
e =
0 2 4
donde el primer uno indica que vamos a tomar el primer elemento de v. El segundo número es el incremento
de ı́ndices dentro de v y el último número marca el elemento final. Esta orden es equivalente a la siguiente
>> e=v(1:3)
e =
0 2 4
ya que cuando el incremento es la unidad, se puede suprimir.
Ejercicio A.3.8 Adivinar cuál va a ser el resultado de las dos órdenes siguientes
>> e=v(2:2:5)
>> e=v(1:3:5)
Como comentamos al principio, la notación usada por Matlab sigue en lo posible la notación estándar de
Álgebra Lineal que se asume conocida. Es muy sencillo multiplicar matrices y vectores, teniendo cuidado de
que las dimensiones sean las adecuadas.
>> A*v(1:3)
??? Error using == *
Inner matrix dimensions must agree.
>> A*v(1:3)’
ans =
16
28
46
Es importante acostumbrase a ver ese mensaje de error. Una vez que se empieza a trabajar con vectores
y matrices, es sencillo olvidar los tamaños de los objetos que se han ido creando.
Ejercicio A.3.9 Utilizando el comando size, razona sobre los problemas en lo que se refiere a dimensiones
en la multiplicación anterior.
Se pueden extraer columnas o filas de una matriz. Si queremos, por ejemplo, que C sea la tercera fila de
la matriz A:
>> C=A(3,:)
C =
6 7 8
O que C sea la segunda columna de la matriz B
>>C=B(:,2)
C =
2
4
7
388 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

O bien que D sea la submatriz cuadrada de orden dos inferior derecha de la matriz A.
>> D=A(2:3,2:3)
D =
4 5
7 8
Una vez que se es capaz de crear y manipular una matriz, se pueden realizar muchas operaciones estándar.
Por ejemplo, calcular su inversa. Hay que tener cuidado y no olvidar que las operaciones son cálculos numéri-
cos realizados por ordenador. En el ejemplo, A no es una matriz regular, y sin embargo Matlab devolverá su
inversa, pues los errores de redondeo durante su cálculo convierten en invertible a dicha matriz.
>> inv(A)
Warning: Matrix is close to singular or badly scaled.
Results may be inaccurate. RCOND = 4.565062e-18
ans =
1.0e+15 *
-2.7022 4.5036 -1.8014
5.4043 -9.0072 3.6029
-2.7022 4.5036 -1.8014
Con la matriz B sı́ que es posible calcular su inversa:
>>inv(B)
ans =
-3.0000 5.0000 -2.0000
-1.0000 -1.0000 1.0000
2.0000 -1.0000 0

Ejercicio A.3.10 Definir una matriz de nombre B1 como la inversa de B. Multiplicar B por B1 y razonar
la coherencia del resultado.

Hay que comentar que Matlab distingue entre mayúsculas y minúsculas. Este puede ser el origen de
algunas confusiones si se manejan algoritmos complejos.
>> inv(a)
??? Undefined function or variable a.

A.4. Cálculo de los autovalores

Hay dos versiones del comando que calcula aproximaciones a los autovalores de una matriz, una de
ellas da solamente los autovalores, mientras que la otra da además sus autovectores correspondientes. Si se
olvida cómo utilizarla, se obtiene esta información pidiendo ayuda sobre esta orden en la lı́nea de comandos
de Matlab.
>> eig(A)
ans =
14.0664
-1.0664
0.0000
>> [V,e] = eig(A)
V =
-0.2656 0.7444 -0.4082
-0.4912 0.1907 0.8165
-0.8295 -0.6399 -0.4082
e =
14.0664 0 0
 Apéndice A: Tutorial de Matlab 389

0 -1.0664 0
0 0 0.0000
>> diag(e)
ans =
14.0664
-1.0664
0.0000

Ejercicio A.4.1 Definir un vector w como la primera columna de la matriz V .

Ejercicio A.4.2 Calcular el producto Aw.

Ejercicio A.4.3 Calcular el producto 14.0664w. ¿Son parecidos? ¿Por qué sı́ o por qué no?

A.5. Resolución de sistemas lineales

También hay funciones para resolver sistemas lineales. Si Ax = b y queremos encontrar x, el modo
más directo es simplemente invertir A, y luego premultiplicar por la inversa ambos lados. Sin embargo,
hay medios mucho más eficientes y estables para resolver sistemas lineales (descomposición LU con pivote
o eliminación gaussiana, por ejemplo). Matlab dispone de comandos especiales que permiten realizar estas
operaciones. Si queremos resolver, por ejemplo, el sistema lineal Bx = v con:

>>v = [1 3 5]’
v =
1
3
5
>>B = [ [1 2 3]’ [2 4 7]’ [3 5 8]’];

utilizaremos eliminación gaussiana con retrosustitución

>> x = B\v
x =
2
1
-1
>> B*x
ans =
1
3
5

Ejercicio A.5.1 Definir una matriz B2 = BB t .

Ejercicio A.5.2 Calcular los autovalores de B2. ¿Qué tienen de especial?

Ejercicio A.5.3 Encontrar la solución del sistema lineal BB t x = v asignando esa solución al vector x.

Ejercicio A.5.4 Comprobar la solución obtenida realizando el cálculo BB T x − v.

Podemos crear una matriz aumentada a partir de B y del término independiente y reducirla hasta
convertir el sistema en uno equivalente triangular, efectuando las necesarias transformaciones elementales de
fila
390 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

>>BA=[B v]
BA =
1 2 3 1
2 4 5 3
3 7 8 5
>>BA(2,:)=BA(2,:)-2*BA(1,:)
BA =
1 2 3 1
0 0 -1 1
3 7 8 5
>>BA(3,:)=BA(3,:)-3*BA(1,:)
BA =
1 2 3 1
0 0 -1 1
0 1 -1 2
La segunda fila tiene el elemento diagonal nulo, ası́ que hay que realizar una permutación de filas,
premultiplicando por la identidad permutada:
>>IP=[1 0 0;0 0 1;0 1 0];
>>BA=IP*BA
BA =

1 2 3 1
0 1 -1 2
0 0 -1 1
Ahora ya es inmediato resolver este sistema por sustitución hacia atrás:
Ejercicio A.5.5 Aplicar lo anterior al problema 2.1.

Ejercicio A.5.6 Definir una matriz H de 3 × 3 a partir de las tres primeras columnas de la matriz BA.

Ejercicio A.5.7 Definir un vector h utilizando la última columna de BA.

Ejercicio A.5.8 Definir el vector z tal que Hz = h. ¿Es coherente el resultado?

Ejercicio A.5.9 Pedir ayuda de la función det.

Ejercicio A.5.10 Calcular el determinante de la matriz H.

A.6. Vectorización de operaciones

Ejercicio A.6.1 Borrar la memoria porque vamos a empezar operaciones nuevas reutilizando nombres de
variables ya usadas.

Con Matlab es sencillo crear vectores y matrices. La potencia de Matlab nace de la facilidad con la que
se pueden manipular estos vectores y matrices. Primero mostraremos cómo realizar operaciones sencillas,
sumar, restar y multiplicar. Luego las combinaremos para mostrar que se pueden realizar operaciones com-
plejar a partir de estas operaciones simples sin mucho esfuerzo. Primero definiremos dos vectores, los cuales
sumaremos y restaremos:
>> v = [1 2 3]’
v =
1
2
3
 Apéndice A: Tutorial de Matlab 391

>> b = [2 4 6]’
b =
2
4
6
>> v+b
ans =
3
6
9
>> v-b
ans =
-1
-2
-3
La multiplicación de vectores y matrices, igual que su suma, sigue las reglas estrictas del Álgebra Lineal.
En el ejemplo anterior, los vectores son ambos vectores columna con tres elementos. No se puede sumar un
vector fila con un vector columna. Se debe recordar que el número de columnas del primer operando debe
ser igual al número de filas del segundo.
>> v*b
Error using == *
Inner matrix dimensions must agree.
>> v*b’
ans =
2 4 6
4 8 12
6 12 18
>> v’*b
ans =
28
Matlab permite realizar las operaciones entre elementos de un vector o matriz de modo muy sencillo.
Supongamos que queremos multiplicar, por ejemplo, cada elemento del vector v con su correspondiente
elemento en el vector b. En otras palabras, supongamos que se quiere conocer v(1) ∗ b(1), v(2) ∗ b(2), y
v(3) ∗ b(3). Serı́a estupendo poder usar directamente el sı́mbolo “∗” pues en realidad estamos haciendo una
especie de multiplicación, pero como esta multiplicación tiene otro sentido, necesitamos algo diferente. Los
programadores que crearon Matlab decidieron usar el sı́mbolo “·∗” para realizar estas operaciones. De hecho,
un punto delante de cualquier sı́mbolo significa que las operaciones se realizan elemento a elemento.
>> v.*b
ans =
2
8
18
>> v./b
ans =
0.5000
0.5000
0.5000
Ejercicio A.6.2 Definir un vector tal que sus componentes sean las de v al cubo.
Una vez que hemos abierto la puerta a operaciones no lineales, ¿por qué no ir hasta el final? Si aplicamos
una función matemática predefinida a un vector, Matlab nos devolverá un vector del mismo tamaño en el
que cada elemento se obtiene aplicando la función al elemento correspondiente del vector original
392 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

>> sin(v)
ans =
0.8415
0.9093
0.1411
>> log(v)
ans =
0
0.6931
1.0986
Saber manejar hábilmente estas funciones vectoriales es una de las ventajas de Matlab. De este modo, se
pueden definir operaciones sencillas que se pueden realizar fácil y rápidamente. En el siguiente ejemplo, se
define un vector muy grande y lo manipulamos de este modo tan sencillo.
>> x = [0:0.1:100]
x =
Columns 1 through 7
0 0.1000 0.2000 0.3000 0.4000 0.5000 0.6000
......................
Columns 995 through 1001
99.4000 99.5000 99.6000 99.7000 99.8000 99.9000 100.0000
>> y = sin(x).*x./(1+cos(x));
Usando este tratamiento vectorial, se pueden generar gráficos de modo muy sencillo. Damos una muestra
de esto que luego completaremos.
>> plot(x,y)

Ejercicio A.6.3 Definir un vector t cuya primera componente sea −4, que tenga un incremento entre
componentes de 0.05 y termine en el punto 1.

Ejercicio A.6.4 Definir un vector y a partir de cada componente del vector t recién definido
2
y = 5e−t + sin(10t)

Ejercicio A.6.5 Dibujar la curva t, y.

A.7. Creación de gráficas

En esta sección presentamos los comandos básicos para crear representaciones gráficas de funciones.
Para mostrar el uso del comando plot, utilizaremos la función seno y su desarrollo limitado en torno al cero,
x − x3 /6. Para dibujar la gráfica, seleccionamos el paso del vector de muestreo x y sus valores primero y
último
>>h=0.1
>>xmin=-2;
>>xmax=2;
>>x=xmin:h:xmax;
>>yseno=sin(x);
>>ytaylor=x-x.^3/6;
Tras esto, tenemos en los vectores yseno e ytaylor los valores reales y los valores aproximados obtenidos
del desarrollo limitado. Para compararlos, dibujamos los valores exactos superpuestos con los aproximados
marcados por puntos verdes ‘o’.
El comando plot se utiliza para generar gráficas en Matlab. Admite una gran variedad de argumentos.
Aquı́ sólo utilizaremos el rango y el formato, y la posibilidad de representar dos curvas en la misma gráfica.
 Apéndice A: Tutorial de Matlab 393

>>plot(x,yseno,’go’,x,ytaylor);
La g se refiere al color verde (green), y la o significa que los puntos se van a marcar con un circulito. El
apóstrofo es el que está en la tecla de la interrogación de cierre.
Ejercicio A.7.1 En la ventana en la que aparece la figura, seleccionar Edit, Copy Figure. Abrir un nuevo
documento de Word y pegar la figura en ese documento.
También es buena idea representar la función error:
>>plot(x,abs(yseno-ytaylor),’mx’);
Para que al final del fichero con todas las órdenes aparezca en pantalla el gráfico, una vez que éste ya
ha sido ejecutado alguna vez, se utiliza la orden shg, que hace que la ventana del gráfico se convierta en la
activa. Usemos este comando para utilizar el comando de petición de ayuda help que es muy útil también
por sus referencias cruzadas a otros comandos.
>> help shg

SHG Show graph window.

SHG brings the current figure window forward.

Ejercicio A.7.2 Pedir ayuda de los comandos grid y plot.

Ejercicio A.7.3 Dibujar la curva t, y del ejercicio A.6.4 con cruces rojas y con una retı́cula incorporada
( grid).

También se puede copiar este gráfico al portapapeles desde la ventana del gráfico, para después pegarlo
en un documento Word por ejemplo, como ya vimos en el ejercicio A.7.1.

A.8. Conjuntos de órdenes

En esta sección explicaremos cómo reunir órdenes en ficheros ejecutables desde la lı́nea de comandos
de Matlab. Ello permite realizar operaciones más complejas, y facilita sus repeticiones.
Para empezar a trabajar sobre esta parte del tutorial, lo primero que haremos es ejecutar clear para
borrar las variables activas.
Como ejemplo, consideramos el fichero correspondiente al dibujo de las gráficas de la sección A.7. Para
ejecutar los comandos del fichero se debe especificar el intervalo entre los valores de las abscisas en el
muestreo. De este modo, se pueden construir infinidad de aproximaciones variando este parámetro.
Primero hay que crear el fichero. El editor más conveniente es el que trae incorporado el propio Matlab.
Este editor es muy simple y suficiente para este tipo de aplicaciones. A partir de la versión 5, viene incor-
porado al propio Matlab. Los ficheros ejecutables de Matlab, los M-files, deben tener la extensión “.m”. En
este ejemplo creamos un fichero de nombre tutorm.m. Para que Matlab ejecute los comandos en el fichero
sólamente hay que ejecutar el comando tutorm.
El fichero se puede guardar dónde se quiera pero se tiene que decirle a Matlab dónde está. Esto se hace
indicando la ruta del archivo en el path browser con el icono correspondiente, o desde el menú File con la
opción Set Path. Por defecto, si se guarda en el directorio ..\matlab\bin, Matlab lo encontrará 4 .
Una vez que el editor aparece en la pantalla (File, New, M-file) debemos ir escribiendo y/o copiando-
pegando los comandos necesarios. Se debe tener en cuenta que cuando una sentencia comienza por %, es un
comentario, y no se va a ejecutar. Por tanto, en este ejemplo, no es necesario reproducir esas lı́neas.
% file: tutorm.m
% Seno y desarrollo del seno.
%
% Para ejecutarlo tienes que fijar el paso
4 Si se utiliza Matlab en el centro de cálculo o laboratorio de una facultad o escuela, probablemente el usuario no tenga

permiso de escritura en ese directorio y no pueda guardar ahı́ sus ficheros. En este caso, se pueden guardar en la carpeta que
se desee que después se incorpora a la ruta de búsqueda (path), bien con el comando path o con el icono correspondiente.
394 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

% h : intervalo entre las x

%
% La rutina genera tres vectores, x con las abscisas, yseno con
% el seno evaluado en esas abscisas, e ytaylor con el desarrollo
% hasta el termino cubico del seno en torno al cero.
%
xmin=-2;
xmax=2;
x=xmin:h:xmax;
yseno=sin(x);
ytaylor=x-x.^3/6;

Una vez que se hayan introducido las sentencias, se salva el fichero volviendo a la ventana con la lı́nea
de comando y se teclea en esta lı́nea el nombre del fichero quitando .m. En este caso tutorm.

>>tutorm
??? Undefined function or variable ’h’.
Error in ==> C:\MATLAB\bin\tut.m
On line 13 ==> x=xmin:h:xmax;

Si se trata de llamar al fichero sin haber definido primero la variable h, aparecerá un mensaje de error.
Se deben definir todas las variables que no se definen en el propio fichero y que éste utiliza.

>>h = 0.1;
>>tutorm
>>plot(x,yseno,’rx’,x,ytaylor)

Una vez ejecutada esta instrucción se deberá ver una gráfica como la de la Figura A.1.
1

0.8

0.6

0.4

0.2

−0.2

−0.4

−0.6

−0.8

−1
−2 −1.5 −1 −0.5 0 0.5 1 1.5 2

Figura A.1: Gráfica correspondiente al ejemplo tutorm.m.

Al teclear tutorm en la lı́nea de comandos, Matlab buscará en los directorios indicados en el path un
fichero llamado tutorm.m. Una vez que lo encuentre lo leerá y ejecutará los comandos como si se hubiesen
tecleado uno detrás de otro en la lı́nea de comandos. Si se desea ejecutar el programa otra vez pero con
un paso diferente, hay que tener cuidado. El programa sobreescribirá los vectores x, yseno e yatylor. Si se
quieren guardar estos vectores hay que especificarlo, almacenándolos en nuevas variables.

>>xp = x;
>>ysenop = yseno;
>>ytaylorp = ytaylor;

Ahora podemos seleccionar un nuevo paso h y volver a ejecutar tutor.

 Apéndice A: Tutorial de Matlab 395

>>h = 0.01;
>>tutorm

Tenemos dos aproximaciones; la primera con un paso h de 0.1 que se almacena en los vectores xp, ysenop
e ytaylorp y la segunda relativa a un paso de 0.01 que guardamos en los vectores x, yseno e ytaylor.

Ejercicio A.8.1 Calcular la dimensión que tienen que tener los vectores x y xp y confirmar el resultado
utilizando la orden size.

Ejercicio A.8.2 Crear y ejecutar desde Matlab un fichero que se llame BAIP.m con la secuencia de coman-
dos siguiente

v = [1 3 5]’;
B = [ [1 2 3]’ [2 4 7]’ [3 5 8]’];
BA=[B v]
BA(2,:)=BA(2,:)-2*BA(1,:)
BA(3,:)=BA(3,:)-3*BA(1,:)
IP=[1 0 0;0 0 1;0 1 0];
BA=IP*BA

Ejercicio A.8.3 Pedir ayuda del comando pause e incorporarlo entre algunas lı́neas del ejercicio anterior
para ver todos los pasos de la secuencia de comandos.

Ejercicio A.8.4 Crear y ejecutar desde Matlab un fichero que se llame CURVATY.m con una secuencia de
comandos que realicen las operaciones siguientes:
1. Borrar todas las variables activas de la memoria.
2. Definir un vector t cuya primera componente sea −4, que tenga un incremento entre componentes de
0.05 y termine en el punto 1.
3. Definir un vector y a partir de cada componente del vector t recién definido como:
2
y = 5e−t + sin(10t)

4. Dibujar la curva (t, y) con cruces rojas y con una retı́cula ( grid) incorporada.

A.9. Matlab y números complejos

Matlab entiende la aritmética compleja y es perfectamente posible trabajar con números complejos.
Podemos multiplicar dos números complejos como:

>>(2+3*i)*(3-7*i)
ans =
27.0000 - 5.0000i
√
Podemos usar indistintamente el sı́mbolo i, o j para −1. El segundo se usa mucho al tratar con los
desarrollos de Fourier tanto en sus versiones continuas como discretas.

>>(2+3*j)*(3-7*j)
ans =
27.0000 - 5.0000i

Podemos utilizar también la fórmula de Euler para definir la exponencial compleja de modo directo,
ejt = cos(t) + j sin(t). La usaremos a menudo en el capı́tulo de aproximación al tratar del desarrollo en serie
de Fourier y de la DFT.

>>exp(j*5.3)
ans =
0.5544 - 0.8323i
396 Problemas de Cálculo Numérico para ingenieros con aplicaciones Matlab

Ejercicio A.9.1 Calcular el valor del seno y del coseno de 5.3 para comprobar que la fórmula de Euler
está bien usada.

Ejercicio A.9.2 Multiplicar exp(jπ/6) por el complejo 1 − j. Comprobar el resultado realizando manual-
mente las operaciones correspondientes.

Podemos también generar un vector de números complejos de modo similar a como lo hemos venido
haciendo. Con las siguientes lı́neas, dibujamos la parte real de la exponencial correspondiente, una función
de periodo 3. Para ello usamos la función real, que asocia al número complejo su parte real.

% complejos.m
clear;
t=-1:0.001:2;
w0=2*pi/3;
a=0.5+0.3*j;
fhat=a*exp(j*w0*t);
plot(t,fhat);
axis([-1.5 2.5 -0.8 0.8]);
shg;

Al empaquetarlas en un fichero y ejecutarlas, obtenemos la gráfica A.2 y el siguiente resultado

>>complejos
Warning: Imaginary parts of complex X and/or Y arguments ignored.
> In complejos.m at line 6

Ejercicio A.9.3 Modificar complejos.m para representar la función suma de dos exponenciales de periodos
2 y 3 y factores 1 − j y 2 − 3j, respectivamente.

Ejercicio A.9.4 Pidiendo ayuda sobre la orden real, obtener la parte imaginaria de f hat en complejos.m.

0.8

0.6

0.4

0.2

−0.2

−0.4

−0.6

−0.8
−1.5 −1 −0.5 0 0.5 1 1.5 2 2.5

Figura A.2: Gráfica resultado del código complejos.m.

A.10. Matemáticas simbólicas con Matlab

Matlab dispone de herramientas para cálculo simbólico. Para ello es necesario instalar el Symbolic
Math Toolbox, que es una especie de versión reducida de Maple, un programa de cálculo simbólico muy
conocido. Aquı́ podemos usar esta caja de herramientas para resolver integrales y calcular determinantes de
modo simbólico entre otras cosas. Lo primero que tenemos que hacer es definir una variable como susceptible
de ser utilizada en cálculo simbólico:
 Apéndice A: Tutorial de Matlab 397

>>syms x

Ahora podemos definir una función que dependa de x y cuya integral queramos calcular:

>>f=cos(x)^2;
>>int(f)
ans=
1/2*cos(x)*sin(x)+1/2*x

Podemos también definir una matriz que dependa de x y de una nueva variable y:

>>syms y
>> A=[x y x-y
2 x^2 y
-x -y 0]
A =
[ x, y, x-y]
[ 2, x^2, y]
[ -x, -y, 0]
Y podemos calcular su determinante de modo simbólico:
>> det(A)
ans =
-2*y*x+2*y^2+x^4-x^3*y

Ejercicio A.10.1 Calcular de modo simbólico la inversa de la matriz A.

Podemos evaluar este determinante para valores reales de x e y asignando valores a esas variables y
utilizando después la orden eval :

>> x=2.41
x =
2.4100
>> y=-3.2
y =
-3.2000
>> eval(det(A))
ans =
114.4301

En el momento en que hemos asignado valores a las variables, éstas dejan de ser sı́mbolos. Si queremos
que vuelvan a serlo tenemos que hacerlo de modo explı́cito

>>syms x
2
Ejercicio A.10.2 Definir una función f como e−x .

Ejercicio A.10.3 Pedir ayuda de la función diff y calcular la derivada de f . Evaluar esta derivada para
x = −3.327.

Ejercicio A.10.4 Pedir ayuda de la función limit y calcular el lı́mite de f cuando x → ∞.