Tema 01 - Introduccion A Matlab

Tema 1
Introducción a Matlab
Prof. Neus. Garrido

Dra. Alicia Cordero
Dr. Juan R. Torregrosa
Máster en Ingenierı́a Matemática y Computación

Escuela Superior en Ingenierı́a y Tecnologı́a
N Garrido + A Cordero + JR Torregrosa Primavera 20’ 1 / 42

Contenido
1 ¿Qué es Matlab?
2 La interfaz gráfica
3 Instrucciones básicas
4 Operaciones con vectores y matrices
Vectores
Matrices
5 Funciones y scripts
6 Estructuras de control
Condicionales
Bucles
7 Representaciones gráficas
Plot
Surf
8 Cálculo simbólico
9 Ejercicios

1
¿Qué es Matlab?

2
La interfaz gráfica

Partes de la interfaz gráfica

Establecer el directorio de trabajo

Archivos de la carpeta actual

Workspace

Editor de funciones, scripts, texto, ...

Ventana de comandos

Historial de comandos

Aspecto que le vamos a dar

Aspecto que le vamos a dar

3
Instrucciones básicas

Comandos
help
save, clear, clc, who, load
Ejemplo 1. Información sobre el comando save

>> help save
Ejemplo 2. Información sobre la función ode45

>> help ode45

Tipos de archivos
.m, .mat, .fig
Ejemplo 3. Generación del archivo holaMundo.m

1. Abrir el editor de Matlab
2. Escribir disp(’Hola Mundo’);
3. Guardar como holaMundo.m
4. Ejecutar en la consola: >> holaMundo
Ejemplo 4. Generación del archivo Notas.mat

1. Introducir en la consola las notas de Modelado y Simulación Numérica: >> MSN=[5
7 8 4 10]
2. Introducir en la consola las notas de Métodos Numéricos: >> MN=[3 7 10 9 8]
3. Introducir en la consola las notas de Geometrı́a Diferencial Aplicada: >> GDA=[9 7
8 5 1]
4. Guardar todas las variables: >> save Notas

Funciones
sin, cos, log, log10, atan, atan2
Ejemplo 5. Cálculo del cos(π/4)

>> cos(pi/4)
Ejemplo 6. Cálculo del log(0)

>> log10(0)

Otros
ans, NaN, ...
π: pi, e: exp(1), i: 1i
Ejemplo 7. Cálculo del último elemento elevado al cuadrado

>> ans^2
π
Ejemplo 8. Cálculo de e−i 2
>> exp(-1i*pi/2)

Formatos numéricos de salida

format short: pi=3.1416
format long: pi=3.141592653589793
format short e: pi = 3.1416e+000
format long e: pi = 3.141592653589793e+000

Formatos numéricos de salida

format short: pi=3.1416
format long: pi=3.141592653589793
format short e: pi = 3.1416e+000
format long e: pi = 3.141592653589793e+000
Tipos numéricos
Booleanos: logical
Enteros sin signo de 2k bits: uintk, k = 3, 4, 5, 6
Enteros con signo de 2k bits: intk, k = 3, 4, 5, 6
Reales: double

4
Operaciones con vectores y matrices

Contenidos

Vectores
Matrices
9 Ejercicios
Vectores
¿Cómo se introducen los vectores?

u=[3,-1,5,6]
Entre corchetes, con comas o espacios entre los elementos

Vectores

u=[3,-1,5,6]

Ejemplo 9. Introduce el vector w = 4 π −1 0 .
>>w=[4 pi -1 0]

Vectores

u=[3,-1,5,6]

>>w=[4 pi -1 0]
Extracción de una componente

u(2)=-1
Entre paréntesis el elemento del que se quiere conocer el valor

Vectores

u=[3,-1,5,6]

>>w=[4 pi -1 0]

u(2)=-1
Ejemplo 10. Obtén la última componente del vector w

>> w(end)

Vectores

u=[3,-1,5,6]

>>w=[4 pi -1 0]

u(2)=-1
Ejemplo 10. Obtén la última componente del vector w

>> w(end)
Ejemplo 11. Obtén las componentes primera y tercera del vector w

>> w([1 3])

Vectores
Vectores con componentes equiespaciadas

u=0:.1:0.7→u=0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
u=4:7→u=4 5 6 7
u=linspace(-3,3,4)→u=-3 -1 1 3

Vectores
Vectores con componentes equiespaciadas

u=0:.1:0.7→u=0 0.1 0.2 0.3 0.4 0.5 0.6 0.7
u=4:7→u=4 5 6 7
u=linspace(-3,3,4)→u=-3 -1 1 3
Ejemplo 12. Obtén un vector de 300 nodos equiespaciados entre los valores 0 y 5
>> x=linspace(0,5,300);
Ejemplo 13. Calcula 34 + 44 + 54 + 64 + 74 .

>> x=3:7;
>> y=x.^4;
>> z=sum(y);

Vectores
Operaciones con vectores

Mismo número de componentes
u+v, u-v, u*v, u.*v, u^2, u.^2, u.\v
cross(u,v), dot(u,v), norm(u), length(u)
T
Ejemplo 14. Sean w = 4 π −1 0 yt= i 0 2 1−i , calcula
T
w + t : w+t.’
w · t: w*t

z = w1 t1 w2 t2 w3 t3 w4 t4 : w.*t.’

Contenidos

Vectores
Matrices
9 Ejercicios
Matrices
¿Cómo se introducen las matrices?

A=[1,-1,3;0,2,-1;1,5,3]
Entre corchetes, con comas o espacios entre los elementos de una fila, y con punto y
coma para separar una fila de otra.

Matrices
¿Cómo se introducen las matrices?

A=[1,-1,3;0,2,-1;1,5,3]
Entre corchetes, con comas o espacios entre los elementos de una fila, y con punto y
coma para separar una fila de otra.
 
1 0 0
Ejemplo 15. Introduce la matriz B =  1 2 −1 
1 5 −1
>> B=[1 0 0;1 2 -1;1 5 1]

Matrices
Extracción de elementos
Una componente: A(2,3)=-1
Varias componentes: A([2,3],2)=2 5
Una fila: A(3,:)=1 5 3
Una columna: A(:,2)=-1 2 5

Matrices
Una fila: A(3,:)=1 5 3
Ejemplo 16. Obtén la primera fila de la matriz B

>> B(1,:)

Matrices
Una fila: A(3,:)=1 5 3
Ejemplo 16. Obtén la primera fila de la matriz B

>> B(1,:)
Ejemplo 17. Obtén un vector b compuesto por los elementos de la diagonal

principal de B
>> b=[B(1,1) B(2,2) B(3,3)]

Matrices
Operaciones con matrices

Respetando las dimensiones de las operaciones
A+B, A*B, A.*B
det(A), inv(A), size(A)
zeros(n,m), ones(n,m), rand(n,m), eye(n)

Matrices
Operaciones con matrices

Respetando las dimensiones de las operaciones
A+B, A*B, A.*B
det(A), inv(A), size(A)
zeros(n,m), ones(n,m), rand(n,m), eye(n)
   
1 0 0 8 1 6
Ejemplo 18. Sean las matrices C =  1 2 −1  y D =  3 5 7 , obtén
1 5 −1 4 9 2
las siguientes operaciones
C + D: C+D
F tal que fij = cij · dij : F=C.*D
|D|: det(D)
C −1 : inv(C)

5
Funciones y scripts

Funciones y scripts
Ambos ejecutan automática y secuencialmente una serie de lı́neas de comandos

Funciones y scripts
Funciones vs scripts
Funciones Scripts
Tiene parámetros de entrada No tiene parámetros de entrada
Útil con pocos y muchos parámetros Útil con pocos parámetros

Funciones y scripts
Funciones Scripts
raizcuadradafunction.m raizcuadradascript.m

Funciones y scripts
Funciones
Tiene parámetros de entrada
Scripts
Útil con pocos y muchos parámetros No tiene parámetros de entrada
raizcuadradafunction.m Útil con pocos parámetros
function y=raizcuadradafunction(x) raizcuadradascript.m
y=sqrt(x);

Funciones y scripts
Funciones Scripts
raizcuadradafunction.m raizcuadradascript.m
function y=raizcuadradafunction(x) x=5;
y=sqrt(x); y=sqrt(x);

Funciones y scripts
Ejemplo 19. Genera el archivo e2g.m que resuelve polinomios de segundo grado
del tipo ax2 + bx + c = 0. Como parámetros de entrada tienes que introducir los
coeficientes a, b y c, y como parámetros de salida las soluciones.
2. Escribir la función
function [sol]=e2g(a,b,c)
sol1=(-b+sqrt(b^2-4*a*c))/(2*a);
sol2=(-b-sqrt(b^2-4*a*c))/(2*a);
sol=[sol1 sol2];

Funciones y scripts
Ejemplo 19. Genera el archivo e2g.m que resuelve polinomios de segundo grado
del tipo ax2 + bx + c = 0. Como parámetros de entrada tienes que introducir los
coeficientes a, b y c, y como parámetros de salida las soluciones.
2. Escribir la función
function [sol]=e2g(a,b,c)
sol1=(-b+sqrt(b^2-4*a*c))/(2*a);
sol2=(-b-sqrt(b^2-4*a*c))/(2*a);
sol=[sol1 sol2];
Ejemplo 20. Ejecuta el programa e2g para obtener la solución de los siguientes
polinomios de segundo grado.
2x2 − 5 = 0: >> s1=e2g(2,0,-5)
x2 + 1 = 0: >> s2=e2g(1,0,1)
x2 − i2x + π = 0: >> s3=e2g(1,-2*1i,pi)

6
Estructuras de control

Contenidos
Condicionales
Bucles
9 Ejercicios
Condicionales
if
if condición
instrucciones
end

Condicionales
if
if condición
instrucciones
end
if condición
intrucciones
else
instrucciones
end

Condicionales
if
if condición
instrucciones
end
if condición
intrucciones
else
instrucciones
end
if condición
intrucciones
else if condición
instrucciones
else
instrucciones
end
end
Condicionales
if Ejemplo 21. if
if nota<5
if condición disp(’Suspenso’)
instrucciones else
end disp(’Aprobado’)
end
if condición
intrucciones
else
instrucciones
end
if condición
intrucciones
else if condición
instrucciones
else
instrucciones
end
end
Condicionales
if Ejemplo 21. if
if nota<5
instrucciones else
end
if condición
intrucciones Ejemplo 22. if
else if nota<5
instrucciones disp(’Suspenso’)
end else if nota<7
disp(’Aprobado’)
else if nota<9
if condición disp(’Notable’)
intrucciones else
else if condición disp(’Sobresaliente’)
end
instrucciones end
else end
instrucciones
end
end
Condicionales
if Ejemplo 21. if
if nota<5
instrucciones else
end
if condición
intrucciones Ejemplo 22. if
else if nota<5
instrucciones disp(’Suspenso’)
end else if nota<7
disp(’Aprobado’)
else if nota<9
if condición disp(’Notable’)
intrucciones else
else if condición disp(’Sobresaliente’)
end
instrucciones end
else end
instrucciones
end Alineación en Matlab
end
Ctrl+A + Ctrl+I
Condicionales
switch
switch expresión
case valor1
instrucciones
case valor2
instrucciones
end

Condicionales
switch
switch expresión
case valor1
instrucciones
case valor2
instrucciones
end
switch expresión
case valor1
instrucciones
case valor2
instrucciones
otherwise
instrucciones
end

Condicionales
switch
switch expresión
case valor1 Ejemplo 23. switch
instrucciones
case valor2 Raı́ces de f (x) = x2 + c para
instrucciones c ∈ {−1, 0, 1}
end switch c
case -1
disp(’Raı́ces reales’)
switch expresión
case 0
case valor1
disp(’Raı́z real’)
instrucciones
otherwise
case valor2
disp(’Raı́ces complejas’)
instrucciones
end
otherwise
instrucciones
end

Contenidos
Condicionales
Bucles
9 Ejercicios
Bucles
while
while condición
instrucciones
end

Bucles
while
while condición
instrucciones
end
Uso de condiciones
Incrementales
n=0;
while n<10
instrucciones
n=n+1;
end

Bucles
while
while condición
instrucciones Ejemplo 24. while
end iter=0;
Uso de condiciones while iter<=20&&fx>1e-6
fx=f(x); dfx=df(x);
Incrementales x=x-fx/dfx;
n=0; iter=iter+1;
while n<10 end
instrucciones
n=n+1;
end

Bucles
for simple
for ı́ndices=vector
instrucciones
end

Bucles
Ejemplo 25. for simple

for simple
v=[4 7 1 9]; v2=[];
for k=1:length(v)
instrucciones
v2(k)=v(k)*2;
end
end

Bucles

for simple
v=[4 7 1 9]; v2=[];
for k=1:length(v)
instrucciones
v2(k)=v(k)*2;
end
end
for doble
for ı́ndices1=vector1
instrucciones
end
end

Bucles

for simple
v=[4 7 1 9]; v2=[];
for k=1:length(v)
instrucciones
v2(k)=v(k)*2;
end
end
Ejemplo 26. for doble

for doble
A=rand(3,4); A2=[];
for fil=1:size(A,1)
for col=1:size(A,2)
instrucciones
A2(fil,col)=A(fil,col)*2;
end
end
end
end

7
Representaciones gráficas

Contenidos
Plot
Surf
9 Ejercicios
Plot
plot(x,y)
Representa en abscisas los valores de x y en ordenadas los valores de y, unidos por lı́neas.
Los vectores tienen que tener las mismas dimensiones.

Plot
plot(x,y)
Opciones
Añadir tı́tulo a los ejes: xlabel(texto), ylabel(texto)
Añadir tı́tulo a la figura: title(texto)
Los tı́tulos permiten texto en LATEX
Guardar figura como *.fig
Añadir mallado a los ejes: grid
Incluir leyenda: legend(texto1,texto2)

Plot
plot(x,y)
Opciones
Añadir tı́tulo a los ejes: xlabel(texto), ylabel(texto)
Añadir tı́tulo a la figura: title(texto)
Los tı́tulos permiten texto en LATEX
Guardar figura como *.fig
Añadir mallado a los ejes: grid
Incluir leyenda: legend(texto1,texto2)
Ejemplo 27. plot

1
x
t=linspace(0,2*pi); 0.5 y
x=cos(t);
0
y=sin(t);
plot(t/2/pi,[x;y]); -0.5
xlabel(’t’), legend(’x’,’y’)
-1
0 0.2 0.4 0.6 0.8 1
t

Contenidos
Plot
Surf
9 Ejercicios
Surf
mesh(x,y,z)
Abscisas: vector x. Ordenadas: vector y. Cota: matriz z.
La matriz tiene length(x) filas y length(y) columnas.

Surf
mesh(x,y,z)
Opciones
Interpolar el resultado: shading interp
Visualizar la barra de colores: colorbar
Cambiar la paleta de colores: colormap mapa de colores

Surf
mesh(x,y,z)
Opciones
Interpolar el resultado: shading interp
Visualizar la barra de colores: colorbar
Cambiar la paleta de colores: colormap mapa de colores
Ejemplo 28. surf

|i\sqrt{c}|
rc=linspace(-2,2,20); ic=linspace(-3,3,40); 2
C=RC+1i*IC; zC=1i*sqrt(C); 1.5
surf(rc,ic,abs(zC)) 1
shading interp 0.5
xlabel(’\Re\{c\}’) 0
3
ylabel(’\Im \{c \}’) 2

1
0 1
2
-1 0
-2 -1
={c} -3 -2 <{c}

8
Cálculo simbólico

Cálculo simbólico
Matlab también puede trabajar con cálculo simbólico, aunque no es tan potente como
Mathematica
Es necesario definir qué variables serán simbólicas
21/12/17 17:23 MATLAB Command Window 1 of 1
Matlab Mathematica
>> syms x c In[3]:= fx_, c_  x ^ 2  c
>> fx=x^2+c
fx = Out[3]= c  x2
x^2 + c
>> v=x+fx In[4]:= vx_, c_  x  fx, c
v =
x^2 + x + c Out[4]= c  x  x2
Mx_, c_  x  fx, c ^ 2  fv, c  fx, c  Simplify
>> fv=v^2+c
fv = In[5]:=
c  x2 
c + (x^2 + x + c)^2
2
>> M=simplify(x-fx^2/(fv-fx))
M = Out[5]= x
(x^3 + x^2 + c*x - c)/(x^2 + 2*x + c) v2  x2
>> solve(M==x,x) In[6]:= SolveMx, c  x, x
x    c , x    c , x   c , x   c 
ans =
(-c)^(1/2)
Out[6]=
-(-c)^(1/2)
>>
>>
>>
>>
>>
9
Ejercicios

Ejercicio 1
Ejercicio 1.
1. Genera un script que exporte a Excel en formato tabla los resultados de las
siguientes funciones en el intervalo x ∈ [−5, 5], tomando 11 valores en cada
intervalo. Muestra por pantalla los gráficos resultantes.
f1 (x) = x2 − x
f2 (x) = cos2 (x)
f3 (x) = exp(−x)
f4 (x) = sin(x) + cos(x)
2. Genera un informe con los valores numéricos y las representaciones gráficas.

Ejercicio 2
Ejercicio 2.
1. Genera un programa (una function que muestre por pantalla la derivada analı́tica y
la derivada aproximada de una función. La derivada aproximada viene dada por la
expresión
f (x + h) − f (x)
f 0 (x) ≈ ,
h
donde h = xi − xi−1 . Como parámetros de entrada tendrás que dar la función y los
valores de x.
2. Modifica dicho programa para que exporte los valores numéricos de la derivada
analı́tica y la derivada aproximada para las siguientes funciones y los siguientes
intervalos, tomando siempre 11 nodos:
f1 (x) = x3 − 3x + 2, x ∈ [−1, 1],
2
√ (x) − ln(x + 5), x ∈ [−3, 0],
f2 (x) = sin
f3 (x) = 3x + exp(−x/2), x ∈ [−1, 2].
3. Genera un informe con los valores numéricos y las representaciones gráficas.

Referencias
A. Cordero, J. L. Hueso, E. Martı́nez, J. R. Torregrosa. Métodos Numéricos con

Matlab. Universidad Politécnica de Valencia, 2005. ISBN 84-9705-854-2.
F. I. Chicharro. Introducción a SciLab. UNIR TV, 2017.
http://tv.unir.net/videos/22053/0/Scilab.


Tema 01 - Introduccion A Matlab

Cargado por

Copyright:

Formatos disponibles

Tema 01 - Introduccion A Matlab

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Tema 01 - Introduccion A Matlab

Cargado por

Copyright:

Formatos disponibles

Tema 1

Prof. Neus. Garrido

Máster en Ingenierı́a Matemática y Computación

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 1 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 2 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 3 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 4 / 42

Partes de la interfaz gráfica

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

Establecer el directorio de trabajo

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

Archivos de la carpeta actual

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

Editor de funciones, scripts, texto, ...

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

Aspecto que le vamos a dar

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

Aspecto que le vamos a dar

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 5 / 42

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 6 / 42

Ejemplo 1. Información sobre el comando save

Ejemplo 2. Información sobre la función ode45

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 7 / 42

Ejemplo 3. Generación del archivo holaMundo.m

Ejemplo 4. Generación del archivo Notas.mat

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 8 / 42

Ejemplo 5. Cálculo del cos(π/4)

Ejemplo 6. Cálculo del log(0)

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 9 / 42

Ejemplo 7. Cálculo del último elemento elevado al cuadrado

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 10 / 42

Formatos numéricos de salida

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 11 / 42

Formatos numéricos de salida

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 11 / 42

Operaciones con vectores y matrices

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 12 / 42

4 Operaciones con vectores y matrices

¿Cómo se introducen los vectores?

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 14 / 42

¿Cómo se introducen los vectores?

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 14 / 42

¿Cómo se introducen los vectores?

Extracción de una componente

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 14 / 42

¿Cómo se introducen los vectores?

Extracción de una componente

Ejemplo 10. Obtén la última componente del vector w

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 14 / 42

¿Cómo se introducen los vectores?

Extracción de una componente

Ejemplo 10. Obtén la última componente del vector w

Ejemplo 11. Obtén las componentes primera y tercera del vector w

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 14 / 42

Vectores con componentes equiespaciadas

N Garrido + A Cordero + JR Torregrosa Primavera 20’ 15 / 42

C=RC+1iIC; zC=1isqrt(C); 1.5