PROGRAMACIÓN EN

COMPUTADORAS CON
MATLAB

José S. Orbegoso López

Programación en Computadora con MatLab – J.

E-mail:jose08052010@[Link]
Orbegoso
 PROGRAMA
•Introducción.
•Operaciones con matrices.
•Manipulación de vectores y
matrices.
•Funciones matriciales.
•Polinomios y procesado de
Programación en Computadora con MatLab – J.

señales.
•Funciones de funciones.
•Gráficos: 2D y 3D.
•Programación.
•Análisis de datos.
Orbegoso

•Análisis numérico
 INTRODUCCIÓN
¿Qué es Matlab?, MATrix LABoratory
MATLAB es un programa para realizar cálculos
numéricos con vectores y matrices. Como caso
particular puede también trabajar con números
escalares, tanto reales como complejos. Una de las
capacidades más atractivas es la de realizar una
Programación en Computadora con MatLab – J.

amplia variedad de gráficos en dos y tres

dimensiones. MATLAB.
Orbegoso
 Introducción_____________________________________
Los elementos básicos del Matlab, como
cualquier otro lenguaje de programación, son:
constantes, variables, operaciones, expresiones y
funciones.

Constante numéricas:
• Números enteros: 2 35 -48
• Números reales: 2. -35.2 48.45
• Máximo de 16 cifras significativas
• Utilizando la letra E a continuación de un
nº con punto decimal [2.2250e-308
Programación en Computadora con MatLab – J.

1.7e+308].
• Números complejos: 2+3i 4*j i,j=(-1)½
Operaciones aritméticas elementales:
Suma: + Multiplicación: * Exponenciación: ^
Resta: - División: /
Orbegoso

Primero exponenciaciones, luego divisiones y

o
 Introducción_____________________________________
Variables: es la etiqueta que identifica una
porción de memoria;
Matlab diferencia entre mayúsculas y minúsculas

Para ver las variables definidas en un instante

determinado se teclea:
>> who
o bien
>> whos
Programación en Computadora con MatLab – J. Orbegoso

Para eliminar alguna variable se ejecuta

>> clear variable1 variable2
Expresiones numéricas: son un conjunto de
números, funciones y variables previamente
definidas, relacionados todos ellos por operadores
aritméticos. Si una expresión es demasiado larga
se indica mediante ...
 Introducción_____________________________________

Formatos: por defecto matlab tiene formato corto

pero se puede elegir entre los siguientes formatos.
• >> format long (14 cifras significativas)
• >> format short (5 cifras significativas)
• >> format short e (notación exponencial)
• >> format long e (notación exponencial)
• >> format rat (aproximación racional)

Variables predefinidas en Matlab:

Programación en Computadora con MatLab – J. Orbegoso

i = (-1)½ pi =  Inf=  NaN= cálculos

indefinidos
eps = < nº que + otro nº=nº coma flotante 2.22e-
16
date = valor de la fecha actual
rand = genera números aleatorios [0 1]
realmin = <nº+ realmax= >nº+
 Introducción_____________________________________
Funciones de matlab:
nombre(argumento)
• sqrt (x) raiz cuadrada
• abs(x) módulo de x
• conj(z) conjugado de un complejo
• real(z), imag(z) parte real e imaginaria de z
respectivamente
• exp(x) calcula ex, siendo x real o complejo
• sin(x) asin(x) [-/2 /2] cos(x) acos(x) [0 ]
tan(x)
• atan(x) [-/2 /2] angle(z) log(x) (en base e)
Programación en Computadora con MatLab – J.

log10(x)
• rats(x) rem(x,y) resto de x/y round(x) sign(x)
Comandos de ayuda:
• help
• lookfor
Orbegoso

• what ficheros .m y .mat del directorio actual

 Introducción_____________________________________

• type nombre_fichero Muestra el contenido del

fichero
• delete nombre_fichero Borra el fichero
• cd cambia de directorio
• pwd indica el directorio actual
• which nombre_fichero india el directorio donde
esta
• ! Abre una ventana de MSDOS que se cierra
Programación en Computadora con MatLab – J. Orbegoso

cuando volvemos a Matlab

startup.m fichero de arranque al ejecutar matlab.

Para guardar en un fichero los comandos que se
ejecutan en una sesión se pone
>> diary nombre_fichero
...
>> diary off
 Introducción_____________________________________
diary [Link]
clear
Ejercicio 1.1 Calcular el valor de la expresión
42.1768  234
J  10
2  10247
Ejercicio 1.2 Calcular el valor de la expresión
9.8 * 1014  5.876 * 10  5
H
9.987 * 10 5  10 6
Y escribir el resultado en al menos 2 formatos
Ejercicio 1.3 Calcular 3sen(32º15' )
I 7
Programación en Computadora con MatLab – J. Orbegoso

42.13
Ejercicio 1.4 Según Hill y Lounasmaa, la ecuación de
la curva de inversión del helio  21.0  5.44T  0.132T 2
P es
Donde la presión vienen dada en atmósfera y la
temperatura T en grados Kelvin. Calcular el valor de
la presión a una temperatura de 293 ºK. Calcular el
valor de la temperatura para una presión de 1N/m2.
Nota: 1N/m2=9.265*10-6atm
 Introducción_____________________________________

diary off
dir
type [Link]
Programación en Computadora con MatLab – J. Orbegoso
 VECTORES Y MATRICES
Las matrices son el tipo
fundamental de dato en Matlab.
» A=[1 3 5; 6 9 2; 4 8 7] » A^2+3*A

A = ans =

1 3 5 42 79 61
86 142 68
Programación en Computadora con MatLab – J.

6 9 2
4 8 7 92 164 106

» det(A) MATrix LABoratory

ans = -- datos son matrices
-- reglas del álgebra lineal
Orbegoso

5
 Vectores y matrices_______________________________
Los vectores pueden ser vectores fila o vectores columna
Vectores fila; los elementos de una misma fila están
separados por blancos o comas,
>> v =[2 3 4]
Vectores columna: los elementos de una misma columna
están separados por intro o por caracteres punto y coma
(;).
>> w =[2;3;4;7;9;8]
La dimensión de un vector se obtiene por el comando
length(nombre del vector) sub_w=w(i:k:j)
Programación en Computadora con MatLab – J.

Generación de vectores:
• Especificando el incremento de sus componentes v=a:h:b;
• Especificando su dimensión linspace(a,b,n) si se omite n
toma 100 por defecto; el incremento es k=(b-a)/(n-1)
• Con componentes logarítmicamente espaciadas
logspace(a,b,n) genera un vector fila de n puntos
logarítmicamente espaciados entre 10a y 10b. Si se omite el
Orbegoso

valor de n se toma 50
 Vectores y matrices_______________________________
Operaciones con escalares:
v+k adición o suma
v-k sustracción o resta
v*k multiplicación
v/k divide por k cada elemento de v
k./v divide k por cada elemento de v
v.^k potenciación cada componente de v esta elevado a k
k.^v potenciación k elevado cada componente de v
Operaciones entre vectores:
v+w adición o suma
v-w sustracción o resta
Programación en Computadora con MatLab – J.

v.*w multiplicación cada elemento de v por el

correspondiente de w
v./w divide cada elemento de v por el correspondiente de w
v.^w potenciación cada componente de v esta elevado al
correspondiente de w
Producto escalar de vectores v*w calcula el producto
Orbegoso

escalar de v por w
 Vectores y matrices_______________________________

Funciones de matlab específicas para

vectores:

sum(v) suma
prod(v) producto
v` transposición de vectores (filas  columnas)
dot(v,w) producto escalar de vectores
cross(v,w) producto vectorial de vectores
Programación en Computadora con MatLab – J. Orbegoso

[y,k]=max(v) valor máximo de las componentes de

un vector k indica la posición, lo mismo para
min(v) valor mínimo
 Vectores y matrices_______________________________

Ejemplo: Cálculo de errores relativos

Supongamos que para resolver una ecuación
diferencial ordinaria hemos utilizado.
-Un método analítico mediante el cual sabemos
que su solución en el intervalo [0,1] es
y(x)=x2+cos(x).
-Un método numérico para aproximar la solución
en el intervalo [0,1] con parámetro de
discretización
Ángulo0.2
(radianes) Solución
aproximada
Programación en Computadora con MatLab – J.

0 1.0030
0.2 1.0234
0.4 1.0825
0.6 1.1869
0.8 1.3342
Orbegoso

1 1.5415
 Vectores y matrices_______________________________

Para definir una matriz no hace falta establecer de antemano su tamaño

(de hecho, se puede definir un tamaño y cambiarlo posteriormente).
MATLAB determina el número de filas y de
columnas en función del número de elementos que se proporcionan (o
se utilizan). Las matrices se definen por filas; los elementos de una
misma fila están separados por blancos o comas, mientras que las
filas están separadas por pulsaciones intro o por caracteres punto y
coma (;). Por ejemplo, el
siguiente comando define una matriz A de dimensión (3x3):
» A=[1 2 3; 4 5 6; 7 8 9]
La respuesta del programa es la siguiente:
A=
123
Programación en Computadora con MatLab – J.

456
789

matriz traspuesta: En MATLAB el apóstrofo (') es el símbolo de

trasposición matricial.

matriz inversa: la inversa de A se calcula con la función inv( )

Orbegoso

B=inv(A).
 Vectores y matrices_______________________________
En MATLAB se accede a los elementos de un vector poniendo el índice
entre paréntesis (por ejemplo x(3) ó x(i)). Los elementos de las matrices
se acceden poniendo los dos índices entre paréntesis, separados por una
coma (por ejemplo A(1,2) ó A(i,j)). Las matrices se almacenan por
columnas (aunque se introduzcan por filas, como se ha dicho antes),
y teniendo en cuenta esto puede accederse a cualquier elemento de una
matriz con un sólo subíndice. Por ejemplo, si A es una matriz (3x3) se
obtiene el mismo valor escribiendo A(1,2) que escribiendo A(4).
Los operadores matriciales de MATLAB son los siguientes:

+ adición o suma
Programación en Computadora con MatLab – J. Orbegoso

– sustracción o resta
* multiplicación
' adjunta (transpuesta o transpuesta conjugada)
^ potenciación
\ división-izquierda
/ división-derecha
.* producto elemento a elemento
./ y .\ división elemento a elemento
.^ elevar a una potencia elemento a elemento
 Vectores y matrices_______________________________
diag(A) Obtencion de la diagonal de una matriz.
sum(diag(A)) calcula la traza de la matriz A
diag(A,k) busca la k-ésima diagonal.

norm(A) norma de una matriz, calcula el máximo de

los valores absolutos de los elementos de A

Generación de matrices:
Programación en Computadora con MatLab – J. Orbegoso

•Generación de una matriz de ceros, zeros(n,m)

•Generación de una matriz de unos, ones(n,m)
•Inicialización de una matriz identidad eye(n,m)
•Generación de una matriz de elementos aleatorios
rand(n,m)
• Matrices con diagonal dada diag(v), diag(v,k)
• [X Y] columnas, [X; Y] filas
 Vectores y matrices_______________________________

Desde la versión 5 de matlab se admiten variables

subindicadas multidimensionalmente
a=ones(2,2,3)
a(:,:,1)= 1 1 a(:,:,2)= 1 1 a(:,:,3)= 1 1
11 11 11

Matrices dispersas o huecas son aquellas que tienen

gran cantidad de elementos nulos
Programación en Computadora con MatLab – J. Orbegoso

>> a= sparse(i,j,c,m,n)
m indica filas, n columnas, c vector que contiene los
elementos no nulos i,j son dos vectores que indican
la posición de cada elemento de c
para visualizar la matriz entera >>full(a)
[i2,j2,c2]=find(a)
 Vectores y matrices_______________________________

Ejemplo1: Tomemos la siguiente matriz

 12  4 0 0 0 
 
 7 3 0 0 8 
 0 0 0 0 0 
 
 0 0  13 11 0 
 
 0 0 2 7  4 
Programación en Computadora con MatLab – J. Orbegoso

Y trabajemos con ella como una matriz dispersa,

para ello definimos el vector de elementos no
nulos, el vector definido por las filas y el vector
definido por las columnas
 Vectores y matrices_______________________________
m=[12,-4,7,3,-8,-13,11,2,7,-4];
f=[1,1,2,2,2,4,4,5,5,5];
c=[1,2,1,2,5,3,4,3,4,5];
a=sparse(f,c,m,5,5)
full(a)
m1=[12,7,-4,3,-13,2,11,7,-8,-4];
f1=[1,2,1,2,4,5,4,5,2,5];
c1=[1,1,2,2,3,3,4,4,5,5];
Programación en Computadora con MatLab – J. Orbegoso

b=sparse(f1,c1,m1,5,5];
full(b)
a=sparse(f,c,m)
s=a+b
p=a*b
[f2,c2,m2]=find(p)
e=full(sparse(f2,c2,m2))
 Vectores y matrices_______________________________
Ejercicio 2.1 Dados los vectores definidos por
x=(1,4.5,7.8) y=(sen(90),cos(45),0)
Realizar los cálculos siguientes:
x+y;el producto escalar de x e y; calcular el
ángulo que forman ambos vectores.
Ejercicio 2.2 Para un laboratorio se compran los
materiales especificados en la tabla siguiente.
Ref. Precio Cantida
arti. d
Programación en Computadora con MatLab – J. Orbegoso

1520 1146 200

1621 3450 250
1428 6225 150
1429 7100 150

Utilizar vectores1628 8500

y el producto 100
de vectores para
calcular el coste de cada producto y el total a pagar
(la tabla de precios no incluye IVA, debiendo por
tanto aplicarse un 16% de incremento.
 Vectores y matrices_______________________________
Ejercicio 2.3 Dadas las matrices a y b calcular:
a+b, a+0, BI. Comprobar que el producto de
matrices no es conmutativo. Elegir la
submatriz de a formada por la primera
columna y la tercera fila y la submatriz de b
formada por la segunda
 1 2 3 y tercera
 4 2 columna
1  y
calcular el producto    
a  2 3 4  b  2 0 4 
 4 2 5  1  2  5
   
Programación en Computadora con MatLab – J. Orbegoso

Ejercicio 2.4 Utilizar matrices para construir una

tabla que contenga: En la 1º columna la variable
grados celsius en el intervalo [0 100] con un paso
de 2. En la 2º columna su valor en grados
fahrenheit y en la 3º en grados kelvin. ¿Cómo
9 * celsius
guardarías en un 32
fichero la tabla anterior?
5
 Vectores y matrices_______________________________

Ejercicio 2.5 Considérese la matriz tridiagonal

definida por:
i, si i  j 
 j, si i  j  1 
 
aij  1 i , j n
 j, si i  j - i 
0, en otro caso 
Introducir la matriz en el espacio de trabajo como
Programación en Computadora con MatLab – J.

matriz dispersa para n=10

Orbegoso
 POLINOMIOS
Los polinomios se representan en matlab por
un vector fila de dimensión n+1 siendo n el
grado del polinomio.

Dado un polinomio
x3+2x
se representa por
>> pol1=[1 0 2 0]
Programación en Computadora con MatLab – J.

para el cálculo de las raíces de un polinomio

existe el comando roots.
>>raices=roots(pol1) (da un vector
columna, aunque pol1 es un vector fila)
Orbegoso
 Polinomios______________________________________
Un polinomio puede ser reconstruido a partir de
sus raíces con el comando poly
>> p=poly(raices) (da un vector fila) **
en caso de que el argumento de poly fuera una
matriz obtendríamos como resultado el polinomio
característico de la matriz. Así mismo si
queremos calcular los autovalores de la matriz
bastaría con calcular las raíces del polinomio
característico.
Programación en Computadora con MatLab – J. Orbegoso

Ejemplo1:
pol2=[2 4 0 1]; % definición del polinomio
2x3+4x2+1
raices=roots(pol2) % cálculo de sus raíces
pol2_n=poly(raices) % reconstrucción del
polinomio
real(pol2_n) % ¿qué pasó?
 Polinomios______________________________________
Para calcular el valor de un polinomio p en un
punto dado x basta con utilizar el comando
polyval
>>y=polyval(p,x)
Ejemplo3:
p=[1 –1 –1 1] % definición del polinomio x3-x2-
x+1
polyval(p,2)
Ejemplo4:
x=-2:0.1:2;
Programación en Computadora con MatLab – J.

Y=polyval(p,x);
Para multiplicar y dividir polinomios tenemos los
comandos especiales conv(p1,p2) y
deconv(p1,p2)
Ejemplo5:
Orbegoso

p1=[1, -2, 1]; p2=[1,1]; p3=conv(p1,p2)

 Polinomios______________________________________
Para conocer el resto de la división de polinomios
basta con poner
>>[p4,r] = deconv(p3,p2)
El comando residue, permite el cálculo del
desarrollo en suma de fracciones simples del
cociente p1/p2.
p2 debe tener raíces reales
El formato del comando es: >>[r,p,k]
Programación en Computadora con MatLab – J. Orbegoso

=residue(p1,p2)
donde:
r= vector columna con los numeradores
p= vector columna con las raices del denominador
pk=
1( xvector
) r (1)con los coeficientes
fila r (n ) del polinomio

independiente.    k(x)
p2( x ) x  p(1) x  p(n )
 Polinomios______________________________________

>>[p1,p2]=residue(r,p,k) hace la operación inversa

Ejemplo6:
Descomponer en fracciones simples el cociente
x3  x2 1
p1=[1 1 0 1] x 3  3x 2  4
p2=[1 –3 0 4]
[r,p,k]=residue(p1,p2)
Programación en Computadora con MatLab – J. Orbegoso

rats(r)
x3  x2 1 35 13 1
3 2
  2
 1
x  3 x  4 9( x  2) 3( x  2) 9( x  1)
[pol1,pol2]=residue(r,p,k)
 Polinomios______________________________________

Para calcular la derivada de un polinomio tenemos

el comando,
>>polyder(p)
Ejemplo7:
Dado el polinomio x3+6x2+1 su derivada es

p=[1, 6,0,1];
d=polyder(p) % es decir 3x2+12x
Programación en Computadora con MatLab – J.
Orbegoso
 Polinomios______________________________________

Ejercicio3.1 Consideremos el polinomio p(x)=x-1.

Calcular p(x)3 e identificar el polinomio obtenido y
calcular sus raíces
Ejercicio3.2 Según Hill y Lounasmaa, la ecuación
de la curva de inversión del helio es P=-21+5.44T-
0.132T2 donde P viene dada en atmósferas y T en
grados Kelvin.
a) Calcular el valor de la presión a una
temperatura de –268.25 ºC.
Programación en Computadora con MatLab – J.

b) Calcular el valor de la temperatura para una

presión de 3N/m2. Nota 1N/m2=9.265e-6 atm.
Ejercicio 3.3 Calcular la solución general de la
EDO 4
x  2x  1
I 
y -y +2y’’’-2y’’+y’-y=0
(5) (4)
dx
Orbegoso

Ejercicio3.4 Calcular al integral x  1

 Polinomios______________________________________
Ejercicio3.1
p=[1 –1];
q=conv(p,p)
q=conv(q,p)
raices=roots(q)
format short
real(raices)

Ejercicio3.2
a) presion=[-0.132 4.55 –21.0];
Programación en Computadora con MatLab – J.

presion0=polyval(presion,273.15-268.25)
b) presion1=3*9.265e-006;
presion(3)=presion(3)-presion1
temp=roots(presion)
Orbegoso
 Polinomios______________________________________
Ejercicio3.3
format long
p=[1 –1 2 –2 1 –1];
raices=roots(p)

der_p=polyder(p)
polyval(der_p,round(raices(2)))

La solución exacta es
y(x)=c1ex+c2sen(x)+c3cos(x)+c4xsen(x)+c5xcox(x)
Programación en Computadora con MatLab – J.

dsolve(‘D5y-D4y+2*D3y-2*D2y+Dy-y=0’,’x’)
Orbegoso
 Polinomios______________________________________
Ejercicio3.4
numerador=[1 0 0 2 1];
denominador=[1 –1];
[cociente,resto]=deconv(numerador,denominador)
4 3 2
 3 4  x x x
I   x  x 2  x  3   dx     3 x  4 ln( x  1)
 x  1 4 3 2

int(‘(x^4+2*x+1)/(x-1)’)
Programación en Computadora con MatLab – J.
Orbegoso
 ECUACIONES LINEALES
Dado un sistema algebraico de ecuaciones
lineales de la forma Ax=b. Lo resolveremos por
métodos clásicos y con funciones propias de
matlab
Algunos comandos cuyos argumentos son
matrices y resultan útiles para la resolución de
Programación en Computadora con MatLab – J. Orbegoso

sistemas.

>>det(A) %determinante de una matriz cuadrada

>>inv(A) %inversa de una matriz cuadrada
>>rank(A) %rango de una matriz el orden del
mayor menor con determinante no nulo
Ejemplo1:
x=ones(4,4);v=[2 2 2 2];y=x+diag(v)
 Ecuaciones lineales________________________________

Para ver si un sistema es compatible

Ax=b compatible >>rank(A)-rank([A,b])
Número de condición de una matriz
>>cond(A)
este número indica la sensibilidad de la solución de
un problema con respecto
Programación en Computadora con MatLab – J. Orbegoso

a cambios relativos en los datos de entrada. (es la

relación entre
el mayor y el menor valor singular de la matriz,
debe ser próximo a 1)
 Ecuaciones lineales________________________________

Ejemplos de resolución de un sistema Ax=b:

Regla de Cramer
Ejemplo2: Dado un sistema Ax=b

A=[1,2,3;2,3,4;4,2,5]
b=[4;5;1];
%comprobamos que la matriz no es singular
Programación en Computadora con MatLab – J. Orbegoso

det(A)
D1=A;D1(:,1)=b
D2=A;D2(:,2)=b
D3=A;D3(:,3)=b
x=[det(D1);det(D2);det(D3)]/det(A)
 Ecuaciones lineales________________________________
La solución por Cramer es tediosa y poco eficaz.
Una primera posibilidad más cómoda.
Si la matriz es cuadrada y su determinante es
distinto de cero.
Ainv=inv(A)
Y el sistema se resolvería como
x=Ainv*b
Otra solución sería utilizar la división matricial
X=A\b
Programación en Computadora con MatLab – J. Orbegoso

A\b, produce un resultado aunque la matriz no

sea cuadrada y aunque no sea invertible.
Por ejemplo si el sistema esta sobredeterminado
(más ecuaciones que incógnitas) se calcula la
solución que minimiza la norma 2 de Ax-b
(solución de mínimos cuadrados)
 Ecuaciones lineales________________________________
Con matrices huecas
Ejemplo3: Consideremos una matriz tridiagonal de
dimensión 32 con la diagonal principal
idénticamente igual a 4 y las diagonales superior
e inferior igual a 2, y estudiamos la variación en el
número de operaciones a realizar según tratemos
la matriz como llena o dispersa. Al resolver el
sistema “Ax=b”, con b=d=[1:32]’;
Programación en Computadora con MatLab – J. Orbegoso

%se define la diagonal inferior

diag=sparse(2:32,1:31,2*ones(1,31),32,32)
%se define la diagonal principal
a=sparse(1:32,1:32,4*ones(1:32));
%se construye la matriz tridiagonal dispersa
a=a+diag+diag’;
%se construye la matriz llena
b=full(a)
 Ecuaciones lineales________________________________

%se define el segundo miembro

d=[1:32]’;

Autovalores y autovectores de una matriz.

>>eig(A) (vector columna)

>>[V,D]=eig(A)
Programación en Computadora con MatLab – J.
Orbegoso
 Ecuaciones lineales________________________________
Ejercicio4.1
Se considera la matriz A. Calcular el determinante
de la matriz A. Resolver el sistema siendo b un
vector columna igual a (1 2 3)
 6 18 2 
 
A=  7  2  4 
 4 10  6 
 
Programación en Computadora con MatLab – J. Orbegoso

Sustituir en la matriz A el elemento de la fila 1

columna 3 por –10 y hacer lo mismo.
Ejercicio 4.2 Dada la matriz M, comprobar el
teorema de Cayley-Hamilton: Toda matriz es raíz
de su polinomio característico
 1 1  2
 
M=   1  2 0 
 3 0 1 

 Ecuaciones lineales________________________________
Ejercicio4.3
Se considera la matriz A. Calcular los autovectores
y autovalores 1 2 0 
 
A=  2 5  1
 4 10  1
 

Ejercicio4.1
A=[6 18 2;7 –2 –4;4 10 –6];
d=det(A)
Programación en Computadora con MatLab – J.

b=[1;2;3];
sol=A\b
A(1,3)=-10
d=det(A)
sol=A\b
Orbegoso

A*sol-b
 Ecuaciones lineales________________________________

Ejercicio 4.2
M=[1 1 –2; -1 –2 0; 3 0 1];
poly(M)
M^3+4*M+13*eye(3)

Ejercicio 4.3
A=[1 2 0; 2 5 -1;4 10 -1];
[V,D]=eig(A)
diag(D)
Programación en Computadora con MatLab – J.
Orbegoso
 GRÁFICOS: 2D Y 3D
Funciones gráficas 2D elementales:
MATLAB dispone de 4 funciones básicas para crear
gráficos 2-D. Estas se diferencian principalmente
por el tipo de escala que utilizan en los ejes Estas
cuatro funciones son las siguientes:

• plot() crea un gráfico a partir de vectores y/o

columnas de matrices, con escalas lineales sobre
Programación en Computadora con MatLab – J.

ambos ejes.
• loglog() ídem con escala logarítmica en ambos
ejes.
• semilogx() ídem con escala lineal en el eje de
ordenadas y logarítmica en el eje de abscisas.
Orbegoso

• semilogy() ídem con escala lineal en el eje de

 GRÁFICOS: 2D Y 3D______________________________
Existen funciones orientadas a añadir títulos al
gráfico, a los ejes, a dibujar una cuadrícula auxiliar,
a introducir texto, etc.
• title('título') añade un título al dibujo
• xlabel('tal') añade una etiqueta al eje de abscisas.
Con xlabel off desaparece
• ylabel('cual') idem al eje de ordenadas. Con
ylabel off desaparece
• text(x,y,'texto') introduce 'texto' en el lugar
especificado por las coordenadas x e y. Si x e y son
Programación en Computadora con MatLab – J.

vectores, el texto se repite por cada par de

elementos.
•gtext('texto') introduce texto con ayuda del ratón:
legend() define rótulos para las distintas líneas o
ejes utilizados en la figura.
Orbegoso

• grid activa una cuadrícula en el dibujo.

 GRÁFICOS: 2D Y 3D______________________________

plot es la función clave de todos los gráficos 2-D

en MATLAB. Ya se ha dicho que el elemento
básico de los gráficos bidimensionales es el
vector. Se utilizan también cadenas de 1, 2 ó 3
caracteres para indicar colores y tipos de línea. La
función plot(), no hace otra cosa que dibujar
vectores. Ejemplo1:
» x=[1 3 2 4 5 3]
Programación en Computadora con MatLab – J. Orbegoso

x=
132453
» plot(x)

Figura 5.1. Gráfico del vector x=[1 3 2 4 5 3].

 GRÁFICOS: 2D Y 3D______________________________

FUNCIÓN PLOT
Ejemplo2
x=0:pi/90:2*pi;
y=sin(x).*cos(x);
plot(x,y)

grid on
Programación en Computadora con MatLab – J. Orbegoso

grid off
xlabel(‘eje x (en radianes)’)
ylabel(‘eje y’)
title(‘y=sen(x)*cos(x)’)
 GRÁFICOS: 2D Y 3D______________________________

FUNCIÓN PLOT
Es posible incluir en el título o en la etiqueta de los
ejes el valor de una variable numérica. Ya que el
argumento de los comandos title, xlabel e ylabel
es una variable carácter, es preciso transformar
las variables numéricas
Programación en Computadora con MatLab – J. Orbegoso

int2str(n) convierte el valor de la variable entera n

en carácter
num2str(x) convierte el valor de la variable real o
compleja x en carácter
 GRÁFICOS: 2D Y 3D______________________________

Texto sobre la gráfica

gtext(’texto’)
text(x,y,’texto a imprimir’)

Calcular las coordenadas de puntos sobre la curva

ginput(n)
[x,y]=ginput(n)
Programación en Computadora con MatLab – J. Orbegoso

Elección del trazo y color de la curva

plot(x,y,’opcion’)

hold on
hold off
 GRÁFICOS: 2D Y 3D______________________________
y yellow . point - solid
m magenta o circle : dotted
c cyan x x-mark -. dashdot
r red + plus -- dashed
g green * star
b blue s square
w white d diamond Opciones de plot
Programación en Computadora con MatLab – J. Orbegoso

k black v triangle (down)

^ triangle (up)
< triangle (left)
> triangle (right)
p pentagram
h hexagram
 GRÁFICOS: 2D Y 3D______________________________
Ejemplo3: Calcular gráficamente las soluciones de
la ecuación
2 x  cos(2 x )
teta=0:pi/360:pi/4;
0.4
f1=(2*teta-cos(2*teta))/2;
2
f2=0.4*ones(size(f1));
figure
plot(teta,f1,'g--',teta,f2,'r')
Programación en Computadora con MatLab – J. Orbegoso

axis square
xlabel('Angulo (radianes)')
gtext('2x-cos(2x))/2')
text(0.2,0.43,'y=0.4')
[teta0,y0]=ginput(1)
title(['Raiz aproximada=',num2str(teta0)])
 GRÁFICOS: 2D Y 3D______________________________

Elección de la escala de los ejes

axis([x0 x1 y0 y1])
axis auto: devuelve la escala a la de defecto
axis off: desactiva los etiquetados de los ejes
desapareciendo los ejes sus etiquetas y la grid.
axis on: lo activa de nuevo
axis xy: sistema de coordenas cartesianas origen
Programación en Computadora con MatLab – J. Orbegoso

en el ángulo inferior izquierdo, eje ox de izqda. A

dcha. y oy de abajo a arriba.
axis equal: los mismos factores de escala para
los dos ejes
axis square: cierra con un cuadrado la región
delimitada por los ejes de coordenadas actuales.
 GRÁFICOS: 2D Y 3D______________________________

Impresión de gráficas

print -dps % PostScript for black and white

printers
-dpsc % PostScript for color printers

-deps % Encapsulated PostScript

-depsc % Encapsulated Color PostScript
Programación en Computadora con MatLab – J. Orbegoso

print -djpeg<nn> % JPEG imagen, nn nivel de

calidad Ejemplo. print -djpeg90 figura1 (nn
75 por defecto)
 GRÁFICOS: 2D Y 3D______________________________
Ejemplo4: plot(vector,Matriz)
x=0:pi/180:2*pi;
y=sin(x);
z=cos(x);
plot(x,y,x,z)
A=[y' z']
plot(x,A)
Función eval se utiliza con funciones definidas
Programación en Computadora con MatLab – J. Orbegoso

con un carácter. y=eval(’caracter’)

Ejemplo5:
f=‘sin(x)-2*cos(x)’;
x=0:pi/90:2*pi;
y=eval(f);
plot(x,y)
axis([0 6 0 2.4]);gtext(‘sen(x)-2cos(x)’)
 GRÁFICOS: 2D Y 3D______________________________
Función fplot se utiliza con funciones definidas
con un carácter. fplot=(f,[0 2*pi ymin
ymax])
Ejemplo6:
f=‘sin(x)-2*cos(x)’;
fplot(f,[0 2*pi],’g—’)

Esta función puede utilizarse también en la

forma:
Programación en Computadora con MatLab – J. Orbegoso

[x,y]=fplot(f,[0 2*pi ymin ymax])

y en este caso se devuelven los vectores x e y,
pero no se dibuja nada.

Llamar una nueva figura figure o referirnos a

una figura ya echa figure(n)
 GRÁFICOS: 2D Y 3D______________________________

FUNCIÓN SUBPLOT

Una ventana gráfica se puede dividir en m

particiones horizontales y n verticales, con
objeto de representar múltiples gráficos en ella.
Cada una de estas subventanas tiene sus propios
ejes, aunque otras propiedades son comunes a
toda la figura. La forma general de este comando
es:
Programación en Computadora con MatLab – J.

subplot(m,n,i) donde m y n son el número de

subdivisiones en filas y columnas, e i es la
subdivisión que se convierte en activa. Las
subdiviones se numeran consecutivamente
empezando por las de la primera fila, siguiendo
Orbegoso

por las de la segunda, etc.

 GRÁFICOS: 2D Y 3D______________________________

Ejemplo7:
subplot(121)
f=‘sin(x)-2*cos(2*x)’;
fplot(f,[0 2*pi])
legend(‘sen(x)-2cos(2x)’)
subplot(122)
fplot(‘sin’,[0 4*pi],’r’)
legend(‘sen(x)’)
Programación en Computadora con MatLab – J.
Orbegoso
 GRÁFICOS: 2D Y 3D______________________________
Otras funciones gráficas 2-D

• bar() crea diagramas de barras.

• barh() diagramas de barras horizontales.
• bar3() diagramas de barras con aspecto 3-D.
• bar3h() diagramas de barras horizontales con aspecto 3-D.
• pie() gráficos con forma de “tarta”.
• pie3() gráficos con forma de “tarta” y aspecto 3-D.
• area() similar plot(), pero rellenando en ordenadas de 0 a y.
• stairs() función análoga a bar() sin líneas internas.
• errorbar() representa sobre una gráfica –mediante barras–
valores de errores.
Programación en Computadora con MatLab – J.

• compass() dibuja los elementos de un vector complejo como un

conjunto de vectores partiendo de un origen común.
• feather() dibuja los elementos de un vector complejo como un
conjunto de vectores partiendo de orígenes uniformemente
espaciados sobre el eje de abscisas.
• hist() dibuja histogramas de un vector.
Orbegoso
 GRÁFICOS: 2D Y 3D______________________________

Representación de polígonos
función especial para dibujar polígonos planos,
rellenándolos de un determinado color.
La forma general es la siguiente:
» fill(x,y,c)
· Si c es un carácter de color
('r','g','b','c','m','y','w','k'), o un vector de valores
[r g b], el polígono se rellena de modo uniforme
Programación en Computadora con MatLab – J.

con el color especificado.

· Si c es un vector de la misma dimensión que x e
y, sus elementos se trasforman de acuerdo con
un mapa de colores determinado, y el llenado del
polígono –no uniforme en este caso–se obtiene
Orbegoso

interpolando entre los colores de los vértices.

Este comando con matrices:
 GRÁFICOS: 2D Y 3D______________________________

» fill(A,B,C)
donde A y B son matrices del mismo tamaño. En
este caso se dibuja un polígono por cada par de
columnas de dichas matrices. C puede ser un
vector fila de colores uniformes para cada
polígono, o una matriz del mismo tamaño que las
anteriores para obtener colores de relleno por
interpolación.
Programación en Computadora con MatLab – J.

Ejemplo8:
x=[1,2,1,0];
y=[0,1,2,1];
figure
fill(x,y,’r’)
Orbegoso

title(‘rombo’)
 GRÁFICOS: 2D Y 3D______________________________
Funciones gráficas 3D elementales:
La función plot3 es análoga a su homóloga
bidimensional plot. Su forma más sencilla es la
siguiente:
» plot3(x,y,z)
Ejemplo9:
teta=0:pi/80:8*pi;
x=1+2*cos(teta);
y=1+2*sin(teta);
Programación en Computadora con MatLab – J.

z=4*teta;
plot3(x,y,z)
axis([-1 3 -1 3 0 120]);
xlabel('eje x')
ylabel('eje y')
Orbegoso

zlabel('eje z')
 GRÁFICOS: 2D Y 3D______________________________
Representación gráfica de superficies.
mesh(x,y,Z),
Creación de una malla [X, Y]=meshgrid(x,y)
Gráfica de la malla construida sobre la superficie z
mesh(X,Y,Z), meshz(X,Y,Z)
Además hace una proyección sobre el plano z=0,
meshc(X,Y,Z), líneas de contorno en el plano z=0
Ejemplo10:
Programación en Computadora con MatLab – J. Orbegoso

x=[Link];y=[0:50];
%Obtenemos la malla del dominio
[X Y]=meshgrid(x,y);
length(x),length(y)
size(X), size(Y)
Z=X.^2-Y.^2;
figure(1);mesh(X,Y,Z)
figure(2);meshz(X,Y,Z)
 GRÁFICOS: 2D Y 3D______________________________
Lo mismo con surf(X,Y,Z), surfc(X,Y,Z), surfl(X,Y,Z)
Una forma distinta de representar funciones tridimensionales es por
medio de isolíneas o curvas de nivel. Con contour(x,y,Z) y con
contour3(X,Y,Z) generamos las líenas de nivel de una superficie.
Existen etiquetas especiales, primero necesitamos saber los valores
del contorno
cs=contour(Z) y luego ponemos clabel(cs) o clabel(cs,v)
pcolor(Z) dibuja una proyección con sombras de color sobre el
plano, la gama de colores está en consonancia con las variaciones
de la matriz Z.
La función surf y pcolor tiene diversas posibilidades referentes a la
forma en que son representadas las facetas o polígonos coloreados.
Programación en Computadora con MatLab – J.

Las tres posibilidades son las siguientes:

• shading flat: determina sombreado con color constante para cada
polígono. Este sombreado se llama plano o flat.
• shading interp: establece que el sombreado se calculará por
interpolación de colores entre los vértices de cada faceta. Se llama
también sombreado de Gouraud
• shading faceted: consiste en sombreado constante con líneas
Orbegoso

negras superpuestas. Esta es la opción por defecto

 GRÁFICOS: 2D Y 3D______________________________
Cuando se desea dibujar una figura con un determinado mapa de colores se
establece una correspondencia (o un mapping) entre los valores de la función
y los colores del mapa de colores.
» caxis([cmin, cmax]) %escala el mapa de colores
colormap(opcion). Distintas escalas de colores
hsv - Hue-saturation-value color map.
hot - Black-red-yellow-white color map.
gray - Linear gray-scale color map.
bone - Gray-scale with tinge of blue color map.
copper - Linear copper-tone color map.
pink - Pastel shades of pink color map.
Programación en Computadora con MatLab – J. Orbegoso

white - All white color map.

flag - Alternating red, white, blue, and black color map.
lines - Color map with the line colors.
colorcube - Enhanced color-cube color map.
vga - Windows colormap for 16 colors.
jet - Variant of HSV.
prism - Prism color map.
cool - Shades of cyan and magenta color map.
autumn - Shades of red and yellow color map.
spring - Shades of magenta and yellow color map.
winter - Shades of blue and green color map.
summer - Shades of green and yellow color map.
 GRÁFICOS: 2D Y 3D______________________________

colorbar(‘horiz’), colorbar(‘vertical’) %barra con la

escala de colores.
brighten ajusta el brillo de color del mapa de
colores

[X,Y,Z]=sphere(n) dibuja una representación de la

esfera unidad con n puntos de discretización

[X,Y,Z]=cylinder(rad,n) dibuja una representación de

un cilindro unidad con n puntos igualmente espaciados
Programación en Computadora con MatLab – J.

cuya sección viene dada por la curva cuyos radios se

guardan en el vector rad. (1,1) y 20 valores por
defecto.

fill3(x,y,z,c) %polígonos tridimensionales

Orbegoso
 GRÁFICOS: 2D Y 3D______________________________

Manipulación de gráficos
view: view(azimut, elev), view([xd,yd,zd]).
view(2)
view(3)
rotate(h,d,a) o rotate(h,d,a,o) h es el objeto, d es
un vector que indica la dirección y a un ángulo, o
el origen de rotación

En el dibujo de funciones tridimensionales, a veces

Programación en Computadora con MatLab – J.

también son útiles los NaNs. Cuando una parte de

los elementos de la matriz de valores Z son NaNs,
esa parte de la superficie no se dibuja, permitiendo
ver el resto de la superficie.
Orbegoso
 GRÁFICOS: 2D Y 3D______________________________

hidden off desactiva le eliminación de líneas

escondidas
hidden on situación previa
Programación en Computadora con MatLab – J. Orbegoso
 GRÁFICOS: 2D Y 3D______________________________

Transformación de coordenadas
[ang,rad]=cart2pol(x,y) %De cartesianas a polares
[ang,rad,z]=cart2pol(x,y,z) %De cartesianas a cilindricas

[x,y]=pol2cart(ang,rad) %De polares a cartesianas

[x,y,z]=pol2cart(ang,rad,z) %De cilindricas a cartesianas

[angx,angz,rad]=cart2sph(x,y,z)
Programación en Computadora con MatLab – J. Orbegoso

%De cartesianas a esfericas

[x,y,z]=aph2cart(angx,angz,rad) %De esfericas a cartesianas

Ejemplo11:
%cilindricas
[ang,rad,z]=cart2pol(sqrt(3),1,2)
%esfericas
[ang1,ang2,rad1]=cart2sph(sqrt(3),1,2)
 GRÁFICOS: 2D Y 3D______________________________
Creación de películas
Para preparar pequeñas películas o movies se pueden utilizar
las funciones movie, moviein y getframe. Una película se
compone de varias imágenes, denominadas frames. La
función getframe devuelve un vector columna con la
información necesaria para reproducir la imagen que se
acaba de representar en la figura o ventana gráfica activa,
por ejemplo con la función plot. El tamaño de este vector
columna depende del tamaño de la ventana, pero no de la
complejidad del dibujo. La función moviein(n) reserva
memoria para almacenar n frames. Una vez creada la
Programación en Computadora con MatLab – J.

película se puede representar el número de veces que se

desee con el comando movie.
Ejemplo12:
for j=1:10
x=0:0.01:2*pi;
plot(x,sin(j*x)/2)
Orbegoso

M(j)=getframe;
 GRÁFICOS: 2D Y 3D______________________________

Ejercicio5.1: Consideremos la ecuación de Van

der Waals. Se considera el benceno para el cual
a=18.78 atml2/mol2, b=0.1208 l/mol. Representar
sobre una misma gráfica las dos subgráficas
correspondientes a:
Isotermas de 100, 200, 300 y 400 ºC
Isobaras de 25, 35, 45 y 55 atm
Cada curva debe ir con trazo diferenciado, con el
texto que indique la isolínea que se ha
Programación en Computadora con MatLab – J.

representado, así como el título de la gráfica y la

a R=0.0821 y
etiqueta de los ejes.
(P  2 )(V  b ) RT
V=[2:100]; V
Orbegoso
 GRÁFICOS: 2D Y 3D______________________________
Ejercicio5.1:
a=18.78;
b=0.1208;
R=0.0821;
subplot(1,2,1)
T=[Link];
V=[2:100];
fac1=R./(V-b);
fac2=a./V.^2;
P=zeros(4,length(V));
P(1,:)=T(1)*fac1-fac2;
P(2,:)=T(2)*fac1-fac2;
Programación en Computadora con MatLab – J.

P(3,:)=T(3)*fac1-fac2;
P(4,:)=T(4)*fac1-fac2;
plot(V,P(1,:),'-',V,P(2,:),'--',V,P(3,:),':',V,P(4,:),'-.')
title('Ecuacion de Van der Waals: Isotermas')
xlabel('Volumen, ltr.')
ylabel('Presion, atm.')
Orbegoso

axis([0, 50,0,15])
legend('T=100ºC','T=200ºC','T=300ºC','T=400ºC')
 GRÁFICOS: 2D Y 3D______________________________
subplot(1,2,2)
P=[Link];
T=zeros(4,length(V));
fac1=(V-b)/R;
T(1,:)=((P(1)+fac2).*fac1)-273.15;
T(2,:)=((P(2)+fac2).*fac1)-273.15;
T(3,:)=((P(3)+fac2).*fac1)-273.15;
T(4,:)=((P(4)+fac2).*fac1)-273.15;
plot(V,T(1,:),'-',V,T(2,:),'--',V,T(3,:),':',V,T(4,:),'-.')
title('Ecuacion de Van der Waals: Isobaras')
xlabel('Volumen, ltr.')
Programación en Computadora con MatLab – J.

ylabel('Temperatura, ºC')
legend('P=25 atm.','P=35 atm.','P=45 atm.','P=55 atm.')
Orbegoso
Programación en Computadora con MatLab – J.
Orbegoso

GRÁFICOS: 2D Y 3D______________________________
 GRÁFICOS: 2D Y 3D______________________________

Ejercicio5.2 Dada la función f(x,y)=xy, obtener

sobre una ventana gráfica las representaciones
siguientes:
-La superficie definida por la función sobre el
dominio [-10,10]*[-10,10].
-Las líneas de contorno sobre la superficie
-La proyección de las líneas de contorno sobre
el domino de definición
-La proyección de las líneas de contorno sobre
Programación en Computadora con MatLab – J.

el plano xy correspondientes a los valores –4,-1,

1 y 4.
Orbegoso
 GRÁFICOS: 2D Y 3D______________________________
Ejercicio5.2
x=[-10:0.5:10]; y=x;
[X,Y]=meshgrid(x,y);
Z=X.*Y;
subplot(221)
mesh(X,Y,Z)
legend('z=xy')
xlabel('eje x')
Programación en Computadora con MatLab – J. Orbegoso

ylabel('eje y')
zlabel('eje z')
title('superficie z=xy')
subplot(222)
contour3(Z)
grid off
xlabel('eje x')
ylabel('eje y')
 GRÁFICOS: 2D Y 3D______________________________
zlabel('eje z')
title('lineas de contorno 3D')
subplot(223)
cs=contour(Z);
contour(x,y,Z)
grid off
clabel(cs)
xlabel('eje x')
Programación en Computadora con MatLab – J. Orbegoso

ylabel('eje y')
title('Proyeccion de las lineas de contorno')
subplot(224)
contour(Z,[-4,-1,1,4])
grid off
xlabel('eje x')
ylabel('eje y')
title('Proyeccion de las lineas de contorno')
Programación en Computadora con MatLab – J.
Orbegoso

GRÁFICOS: 2D Y 3D______________________________
 GRÁFICOS: 2D Y 3D______________________________

Ejercicio5.3 Representar la superficie de

revolución obtenida al girar la curva y=x2+1
alrededor del eje ox
x=[0:0.1:1];
Programación en Computadora con MatLab – J. Orbegoso
 GRÁFICOS: 2D Y 3D______________________________
Ejercicio5.3
x=[0:0.1:1]; %puntos de discretizacion del eje ox
rad=x.^2+1; %vector de radios
n=length(rad); %numero de radios
cylinder(rad,n) %representacion del cilindro
xlabel('eje x')
ylabel('eje y')
zlabel('eje z')
[X,Y,Z]=cylinder(rad,n);
Programación en Computadora con MatLab – J. Orbegoso

h=surf(X,Y,Z); %calculo del objeto

rotate(h,[0,1,0],90)
%Al rotar desaparecen las etiquetas de los ejes
xlabel('eje x')
ylabel('eje y')
zlabel('eje z')
view(15,15) %cambiamos el punto de observacion
grid off
legend('f(x)=x^2+1')
Programación en Computadora con MatLab – J. Orbegoso

GRÁFICOS: 2D Y 3D______________________________
 PROGRAMACIÓN
Todos los ficheros de comandos matlab deben de
llevar la extensión .m hay 2 tipos
Ficheros de función la primera línea es ejecutable
y comienza con la palabra function
Ficheros de programa no constituyen funciones y
se construye mediante una secuencia de
Programación en Computadora con MatLab – J. Orbegoso

comandos se ejecuta tecleando el nombre sin

extensión
Ejemplo1:
function temp_c=convert(temp_f);
%CREAR UNA FUNCION
% PASA DE GRADOS ºF A ºC
temp_c=5/9*temp_f-5/9*32;
 PROGRAMACIÓN________________________________

Comandos de entrada salida

input: nos permite introducir datos
variable=input(‘mensaje a pantalla’);

disp: muestra un texto por pantalla

disp(‘El algoritmo no ha convergido’)

menu: genera un menú que permite al usuario

Programación en Computadora con MatLab – J. Orbegoso

elegir entre distintas opciones

opcion=menu( ‘titulo del
mensaje’,’opcion1’,...’opcionp’ )

error: informa de la existencia de un error y

detiene la ejecución del programa devolviendo el
control al teclado
 PROGRAMACIÓN________________________________

Programación de funciones
La primera línea es ejecutable y empieza por la
palabra function de la forma
function arg_salida=nombre_funcion(arg_entrada)

Después cuantos comandos sean necesarios

incluidos comentarios como si se tratase de un
fichero programa
Programación en Computadora con MatLab – J. Orbegoso

El fichero se debe guardar con nombre_funcion.m

Para devolver el control al programa desde

cualquier punto de una función basta con escribir la
sentencia return

Ejemplo2:Construir una función que defina para un

gas ideal la temperatura dependiendo de la presión
 PROGRAMACIÓN________________________________

function t=gases(p,v,n)
%t=gases(p,v,n)
%funcion que considra la ley de los gases ideales
%Argumentos de entrada:
%presion p (atmosferas)
%volumen v (litros)
%número de moles n
%Argumentos de salida:
%temperatura t (grados kelvin)
%R=0.0821 [Link]/[Link]
Programación en Computadora con MatLab – J.

R=0.0821
t=p*v/(n*R);

Lo grabo como gases.m

Orbegoso

temp=gases(20,10,10)
 PROGRAMACIÓN________________________________

-El m-fichero solo puede ejecutarse a continuación

del >> de Matlab o como una línea de otro fichero,
una función puede ser llamada desde cualquier
expresión.
-El m-fichero no admite argumentos de entrada,
simplemente trabaja con datos existentes en el
espacio de trabajo.
-Las variables de una función son locales de la
Programación en Computadora con MatLab – J. Orbegoso

función y las de un fichero-m son globales.

Si queremos que una variable sea compartida por
varias funciones a de definirse en todas ellas como
global
global variable
echo escribe cada comando del fichero sobre la
pantalla
pause detiene la ejecución hasta pulsar una tecla
 PROGRAMACIÓN________________________________

Testear el número de argumentos

nargin devuelve el número de argumentos de

entrada con los que el usuario ha llamado a la
función
nargout devuelve el número de argumentos de
salida con los que el usuario ha llamado a la función
nargchk chequea si el número de argumentos de
entrada calculados con nargin esta entre el valor
máximo y mínimo previsto, si no da error.
Programación en Computadora con MatLab – J.

nargchk(mínimo,máximo,número_entrada)
Orbegoso
 PROGRAMACIÓN________________________________

Bucles. Bucles anidados

for k=n1:incre:n2
end

for k=vector_columna
end

Con break se rompe la ejecución

Estructuras de control condicionadas
Programación en Computadora con MatLab – J.

if
if if
if elseif whil
else elseif e
end else
end end end
Orbegoso

end
 PROGRAMACIÓN________________________________
Operadores relacionales y lógicos
Menor: < Menor o igual:<= Mayor: > Mayor o
igual: >=
Igual: == Distinto: ~=
O: | Y: & no: ~
El resultado de realizar operaciones relacionales o
lógicas será un 1 si es verdadera o un 0 si es falsa
Otras funciones relacionales o lógicas:
xor(x,y): “or” exclusivo 0 si x o y son no nulos y 1 en
otro caso.
any(x): si x es un vector devuelve un 1 si alguna
Programación en Computadora con MatLab – J.

componente de x es no nula. Si es una matriz

devuelve un vector fila con un 1 para cada columna
de la matriz x que tenga alguna de sus filas no nulas
y 0 en otro caso.
all(x): Si es un vector devuelve un 1 si todas sus
Orbegoso

componentes son no nulas. Si es una matriz

 PROGRAMACIÓN________________________________
isempty(x): devuelve un 1 si la matriz x es vacía y 0
en caso contrario.
isequal(x1,x2,…,xn): Devuelve un 1 si todas las
matrices son idénticas y un 0 en caso contrario.
Estructuras de control predefinidas:
Ejemplo3:
q=[];v=[1,-1,2,-1];
for i=1:length(v)
Programación en Computadora con MatLab – J. Orbegoso

if(v(i)>=0)
q=q[q,1];
else
q=[q,0];
end
end
q=v>=0
 PROGRAMACIÓN________________________________
r=v==q
s=(v>=0)&(v<2)
y=v(abs(v)>=2)
A=[2,4,6;1,3,5];
[i,j]=find(A<3)
Desde Matlab 5 se ha incorporado la estructura
switch que permite realizar bifurcaciones en un
programa atendiendo al resultado de una expersión
Programación en Computadora con MatLab – J. Orbegoso

switch expresion
case expresion_0
comandos_0 que deben ejecutarse
case expresion_1
comandos_1 que deben ejecutarse
case expresion_2
comandos_2 que deben ejecutarse
otherwise
comandos escoba que deben ejecutarse
end
 PROGRAMACIÓN________________________________

Lectura y escritura en ficheros externos

La lectura y escritura de información en ficheros

externos se lleva a cabo esencialmente con los
comandos
fread y fscanf para lectura
fprintf y fwrite para escritura
El procedimiento general en todos los casos es:
Programación en Computadora con MatLab – J. Orbegoso

-Abrir el fichero del que se desea leer o en el cual

deseamos escribir.
-Colocar el puntero de lectura o escritura en la
posición deseada.
-Leer o escribir las variables.
-Cerrar el fichero.
 PROGRAMACIÓN________________________________
-Abrir el fichero el comando es fopen y su estructura
ident=fopen(‘nombre de fichero’)
En ident se guarda un número de identificación si es
-1 es que el fichero no se ha podido abrir.
-Cerrar ficheros el comando es fclose y su estructura
fclose(ident)
fclose(‘all’)
verif=fclose(‘all’)
-Posicionamiento del puntero.
1.-Posicionar el puntero al inicio del archivo con
Programación en Computadora con MatLab – J.

número de identificación ident:

frewind(ident)
2.-Posicionar el puntero dentro de un archivo
test=fseek(ident,posi,’origen’)
Orbegoso
 PROGRAMACIÓN________________________________
Esta sentencia coloca el puntero del archivo ident en
la posición indicada en posi:
.-Si posi>0 mueve posi bytes hacia delante.
.-Si posi<0 mueve posi bytes hacia atrás.
.-Si posi=0 el puntero no se mueve.
La variable carácter origen indica desde donde se
empieza a mover el puntero.
.-’bof’: principio del archivo.
.-’cof’: posición actual del archivo.
.-’eof’: final del archivo.
Programación en Computadora con MatLab – J.

La ejecución del comando fseek devuelve la variable

test que valdrá -1 si algo falla.
Para conocer la posición del puntero:
posi=ftell(ident)
Lectura de datos debemos distinguir entre ficheros
Orbegoso

formateados y no formateados.
 PROGRAMACIÓN________________________________
Lectura de datos formateados:
[datos,contador]=fscanf(ident,’formato’,cuantos)
1.-Lee datos del archivo ident.
2.-Los datos leídos se guardan en datos.
3.-cuantos, indica cuántos datos vamos a leer.
.-un escalar k
.-un vector[n,m], en este caso se leeran los
datos necesarios para rellenar una matriz de orden
n x m columna a columna.
.-inf todos los datos del archivo.
Programación en Computadora con MatLab – J.

4.-La variable contador indica cuántos se han leido

con éxito.
5.-La variable formato indica el formato de lectura.
%d: decimales
Orbegoso

%e: notación
exponencial
 PROGRAMACIÓN________________________________

%f: not. Pto. fijo

%g: no considera los ceros no
significativos.
%s: variable carácter
Lectura de datos no formateados:
datos=fread(ident,cuantos,’precision’)
Programación en Computadora con MatLab – J.

Escritura en fichero
Escritura de datos formateados:
contador=fprintf(ident,’formato’,datos,controles)
\n línea nueva
Orbegoso

\t avanza hasta la siguiente posición de

 PROGRAMACIÓN________________________________

Escritura de datos no formateados:

contador=fwrite(ident,datos,’precision’)

Ejercicio6.1:
Utilizando estructuras de control construir una
función que nos calcule las raíces de una ecuación
Programación en Computadora con MatLab – J. Orbegoso

de segundo grado
ax2+bx+c=0
 PROGRAMACIÓN________________________________
Ejercicio6.1:
function raiz=sole2(a,b,c)
%raiz=sole2(a,b,c)
%solucion de la ecuacion de segundo grado
%ax^2+bx+c=0, a~=0
%
if (nargin ~=3)
error(‘El numero de argumentos de entrada debe
ser 3’)
Programación en Computadora con MatLab – J. Orbegoso

end
discri=b^2-4*a*c;
raiz=[];
if(discri==0)
raiz=-b/(2*a);
disp([‘Raiz doble=‘,num2str(raiz)])
return
 PROGRAMACIÓN________________________________
raiz(1)=(-b+sqrt(discri))/(2*a);
raiz(2)=(-b-sqrt(discri))/(2*a);
disp([‘Raices reales simples=‘,num2str(raiz)])
return
else
raiz(1)=(-b+sqrt(-discri)*i)/(2*a);
raiz(2)=(-b-sqrt(-discri)*i)/(2*a);
disp([‘Raices complejas=‘,num2str(raiz)])
return
Programación en Computadora con MatLab – J. Orbegoso

end

Esta función se ejecuta:

raiz=sole2(2,3,1);
raiz=sole2(1,0,1);
 PROGRAMACIÓN________________________________
Ejercicio6.2:
Problema para lectura de datos
Nombre Cálculo Álgebra Física Estadística
Fernando Gómez Pereira 3 6 5 7
Susana Rodríguez Pérez 7 4 3 1
Carlos Leis Álvarez 8 9 7 9
Arturo Gómez Álvarez 5 4 5 9
Silvia Tais Álvarez 10 9 10 9
Andrea Gallego Nimes 3 3 2 4
Alicia Caballero Leis 6 8 8 5
Antonio Fraga Gómez 5 7 6 5
Programación en Computadora con MatLab – J. Orbegoso

Beatriz Machado Gómez 4 3 5 4

Laura Tobío Manzanal 7 8 5 9
Juan Rico Fraga 4 7 5 5
Andrés Pena Gómez 6 8 5 6
Luis Blanco Villa 8 6 6 4
Sandra Puentes Gallego9 9 7 9
Isolina Prieto Gómez 5 5 6 6
Teresa Sieiro Gon 4 2 5 3
Ricardo López Amigo 8 6 2 9
 PROGRAMACIÓN________________________________
Ejercicio6.2:
ident=fopen('datos.m');
frewind(ident)
variables=[];
for i=1:5
variable=fscanf(ident,'%s',1);
long_v(i)=length(variables)+length(variable);
variables=[variables variable];
if(i==1)
Programación en Computadora con MatLab – J. Orbegoso

fprintf(1,'\t%s\t\t',variable)
else
fprintf(1,'\t%s\t',variable)
end
end
nombres=[];
califica=[];
 PROGRAMACIÓN________________________________
for i=1:17
nombre=fscanf(ident,'%s,%c',1);
apellido1=fscanf(ident,'%s,%c',1);
apellido2=fscanf(ident,'%s,%c',1);
l=3*(i-1)
long_n(l+1)=length(nombres)+length(nombre);
long_n(l+2)=long_n(l+1)+length(apellido1);
long_n(l+3)=long_n(l+2)+length(apellido2);
nombres=[nombres nombre apellido1 apellido2];
Programación en Computadora con MatLab – J. Orbegoso

califica=[califica; fscanf(ident,'%i',4)'];
fprintf(1,'\n%9s %s %9s',nombre, apellido1,
apellido2)
fprintf(1,'\t%i\t\t%i\t\t%i\t\t%i',califica(i,:))
end
fclose(ident)
 ANÁLISIS DE DATOS
Para realizar con Matlab análisis estadísticos de
un conjunto de datos estos deben ser
almacenados utilizando matrices. Cada columna
de la matriz representará una variable medida y
cada fila los valores que toman las variables
consideradas en un determinado punto de
Programación en Computadora con MatLab – J. Orbegoso

medida.
Veamos algunos comandos interesantes:
Cálculo del mínimo, máximo, media y mediana.
[y,k]=min(x)
[y,k]=max(x)
m=mean(y)
me=median(y) (la mediana de un conjunto de
números ordenados en magnitud es o el valor
 ANÁLISIS DE DATOS______________________________

Desviación típica.
s=std(y)
Nos da la dispersión o variación de los datos para
dar una idea de cuan esparcidos están.
  y (i )  y 2 
std ( y )   i 
 N 
 
Coeficientes de correlación y matriz de
covarianza
Programación en Computadora con MatLab – J.

r=corrcoef([x,y])
Nos da el grado de relación entre x e y.

N
x ( i ) y ( i )
r  i 1

 N
x ( i ) 2
 N
y ( i ) 2

Orbegoso

i 1 i 1
 ANÁLISIS DE DATOS______________________________

Ejemplo1:
R=corrcoef(califica)
 1.0000 0.6438 0.4563 0.4628 
 
 0.6438 1.0000 0.6164 0.6089 
 0.4563 0.6164 1.0000 0.3332 
 
 0.4628 0.6089 0.3332 1.0000 

R(2,3) = R(3,2) ; R(4,3) poca relación entre física
y estadística
Programación en Computadora con MatLab – J.

Para calcular la covarianza (numerador de la

expresión r=...)


N
s=cov([x,y]) x ( i ) y ( i )
s i 1
N1
Orbegoso
 ANÁLISIS DE DATOS______________________________
En el caso de matrices S=cov(X) nos daría la
matriz de covarianza y su diagonal se llama
vector de varianzas = desviaciones típicas al
cuadrado.
Gráficos estadísticos.
Histograma:
>>hist(y)
Ejemplo2:
y=rand(40,1);
hist(y)
Programación en Computadora con MatLab – J.

Variantes: hist(y,n) (con n subintervalos de clases

regulares); hist(y,x); [n,p]=hist(y)
Gráfico de barras:
bar(y)
Perfil de muestra:
Orbegoso

stairs(y)
 ANÁLISIS DE DATOS______________________________
Gráfico de errores:
errorbar(x,y,e)
Histograma angular:
rose similar a hist pero realiza un histograma
angular los valores de la muestra es de ángulos
en radianes.

Curvas de regresión
p=polyfit(x,y,n)
Programación en Computadora con MatLab – J. Orbegoso

Calcula el polinomio de regresión de grado n; es

decir el polinomio p de grado 2
n que minimiza.
N
(ver ejercicio 7.2)
 p( x )  y 
i 1
i i
 ANÁLISIS DE DATOS______________________________
Interpolación uni y bidimensional
.-1D
vector_y=interp1(x,y,vector_x,opcion)
opcion:
-’linear’: interpolación lineal
-’cubic’ :interpolación cúbica
-’spline’: interpolación spline cúbica (ptos de
interpolación igualmente espaciados).
Programación en Computadora con MatLab – J. Orbegoso

.-2D
matriz_Z=interp2(X,Y,Z,matriz_X,matriz_Y,opcion
)
opcion:
-’bilinear’: interpolación lineal
-’bicubic’ :interpolación cúbica
-’nearest’
 ANÁLISIS DE DATOS______________________________
Ejercicio7.1:
La tabla siguiente recoge el peso de 30
estudiantes. Construir sobre una misma ventana
las 4 figuras siguientes:
1.-Un histograma de frecuencias con el ox peso y
oy frecuencia de valores.
2.-El polígono de frecuencias, curva obtenida
entre los puntos definidos por las marcas de
clase y la frecuencia.
Programación en Computadora con MatLab – J. Orbegoso

3.-El perfil de muestra mediante stairs

4.-El perfil de muestra mediante stem.
Calcular el máximo y mínimo peso y la media de
pesos. 71 82 65 75 77 91 59 84 89
81
73 91 82 75 96 85 69 76 81
92
 ANÁLISIS DE DATOS______________________________
Ejercicio7.1:
peso=[71 82 65 75 77 91 59 84 89 81 ...
73 91 82 75 96 85 69 76 81 92 ...
84 79 77 95 81 79 84 85 76 82];
figure
subplot(221)
hist(peso)
title(‘Histograma’)
xlabel(‘Peso’)
Programación en Computadora con MatLab – J. Orbegoso

ylabel(‘Frecuencia’)
subplot(222)
[n,p]=hist(peso);
plot(p,n)
title(‘poligono de frecuencias’)
xlabel(‘Peso’)
ylabel(‘Frecuencia’)
subplot(223)
 ANÁLISIS DE DATOS______________________________

stairs(peso)
title(‘Stairs’)
ylabel(‘Peso’)
xlabel(‘Individuo’)
subplot(224)
stem(peso)
title(‘Stem’)
ylabel(‘Peso’)
Programación en Computadora con MatLab – J. Orbegoso

xlabel(‘Individuo’)
maximo=max(peso)
minimo=min(peso)
media=mean(peso)
 ANÁLISIS DE DATOS______________________________

Ejercicio7.2
Se ha medido experimentalmente la
conductividad eléctrica del acero a distintas
temperaturas, recogiéndose los siguientes
valores:
Obtener los polinomios de regresión de primero y
segundo grado. Calcular para los polinomios
obtenidos cuál sería el valor estimado de la
Programación en Computadora con MatLab – J. Orbegoso

conductividad eléctrica
T(ºC) a K(cm)
los 600 -1
y a los 1000ºC.
100 51813
300 28571
500 17483
700 11696
900 9116
 ANÁLISIS DE DATOS______________________________

Ejercicio7.2
temp=[Link];
conduc=[51813 28571 17483 11696 9116];
pol1=polyfit(temp,conduc,1)
temp1=[Link];
conduc1=polyval(pol1,temp1)
conduc1e=polyval(pol1,temp);
plot(temp1,conduc1,temp,conduc1e,’o’,temp,conduc,’
Programación en Computadora con MatLab – J. Orbegoso

*’)
legend(‘recta de regresion’,’conductividades
estimadas’,...
‘conductividades medidas’)
xlabel(‘temperatura ºC’)
ylabel(‘Conductividad electrica (ohm cm)^{-1}’)
cond1_600=polyval(pol1,600)
cond1_1000=polyval(pol1,1000)
 ANÁLISIS DE DATOS______________________________

pol2=polyfit(temp,conduc,2)
temp2=[Link];
conduc2=polyval(pol2,temp2)
conduc2e=polyval(pol2,temp);
figure
plot(temp2,conduc2,temp,conduc2e,'o',temp,conduc,'
*')
legend('recta de regresion','conductividades
Programación en Computadora con MatLab – J. Orbegoso

estimadas',...
'conductividades medidas')
xlabel('temperatura ºC')
ylabel('Conductividad electrica (ohm cm)^{-1}')
cond2_600=polyval(pol2,600)
cond2_1000=polyval(pol2,1000)
corre1=corrcoef([conduc',conduc1e'])
corre2=corrcoef([conduc',conduc2e'])
 ANÁLISIS DE DATOS______________________________

Ejercicio7.3:
[x,y]=meshgrid(-[Link]);
z=peaks(x,y);
Dada esta función ver la diferencia entre la
interpolación nearest, bilinear y bicubic
Programación en Computadora con MatLab – J. Orbegoso
 ANÁLISIS DE DATOS______________________________

Ejercicio7.3:

[x,y]=meshgrid(-[Link]);
z=peaks(x,y);
[xi,yi]=meshgrid(-3:0.25:3);
zi1=interp2(x,y,z,xi,yi,'nearest');
zi2=interp2(x,y,z,xi,yi,'bilinear');
zi3=interp2(x,y,z,xi,yi,'bicubic');
figure
mesh(xi,yi,zi1)
Programación en Computadora con MatLab – J.

figure
mesh(xi,yi,zi2)
figure
mesh(xi,yi,zi3)
Orbegoso
 ANÁLISIS DE DATOS______________________________

Ejercicio7.4:
Para el ejemplo del ejercicio 6.2 calcular:
.- La nota media obtenida por los alumnos en
cada una de las disciplinas.
.- La nota máxima en estadística y los alumnos
que la han obtenido
Programación en Computadora con MatLab – J. Orbegoso
 ANÁLISIS DE DATOS______________________________

Ejercicio7.4:
mean(califica);
for i=1:4
materia=variables(long_v(i)+1:long_v(i+1));
disp([‘La nota media de ‘,materia, ...
‘es ...’,num2str(ans(i))])
end
yM=max(califica(:,4)); %calculo de la calificacion
maxima
Programación en Computadora con MatLab – J.

disp([‘La nota máxima es’, num2str(yM),’ obtenida

por : ‘])
%identificación alumno
iM=find(califica(:,4)==yM);
%se extrae su nombre y apellidos
Orbegoso
 ANÁLISIS DE DATOS______________________________
Ejercicio7.4:
for i=1:length(iM)
n=3*(iM(i)-1);
if n==0
p=[0 long_n(1:3)];
else
p=long_n(n:n+3);
end
alumno=[nombres(p(1)+1:p(2)) ‘ ‘
nombres(p(2)+1:p(3)) ...
Programación en Computadora con MatLab – J.

‘ ‘ nombres(p(3)+1p(4))];
disp(alumno)
end
Orbegoso
 ANÁLISIS NUMÉRICO
En este capítulo nos introduciremos al análisis
numérico de problemas básicos relacionados con
el estudio de funciones de una y varias variables,
la aproximación de extremos locales o el cálculo
de ceros de una función.

El cálculo de los mínimos locales de una función

de una variable se realiza utilizando la sentencia.
Programación en Computadora con MatLab – J.

min=fminbnd(‘funcion’,a,b,opciones)
Ejemplo1:
Calcular un mínimo local de la función f(x)=3x4-
4x3 en el intervalo [-1 2]
Orbegoso

fminbnd(‘3*x^4-4*x^3’,-1,2)
 ANÁLISIS NUMÉRICO_____________________________
Para calcular un máximo de f en el intervalo [a,b] es
lo mismo que calcular un mínimo sobre –f
fminbnd('-(3*x^4-4*x^3)',-1,2)
fminbnd('-(3*x^4-4*x^3)‘,0,2)

El cálculo de mínimos locales para funciones de

varias variables puede llevarse a cabo ejecutando
Programación en Computadora con MatLab – J. Orbegoso

min=fminsearch(‘funcion’,x0)
La diferencia es que en lugar de dar un intervalo
damos un vector x0 que indica el pto. Entorno al
cual deseamos minimizar la función.
Ejemplo2:
Minimizar la función f(x)=sen(xy) en el entorno de
[0,0]
 ANÁLISIS NUMÉRICO_____________________________
Para calcular los ceros de un función tenemos el
comando fzero
raiz=fzero(‘funcion’,x0)
Calcula la raíz de la función f(x)=0 a partir del
iterante inicial x0
Ejemplo3:
Calcular una solución de la ecuación sen(x)-2cos(2x)
+x2=2-2.
Programación en Computadora con MatLab – J. Orbegoso

Esto es f(x)=sen(x)-2cos(2x)+x2-2+2 =0 .
Función continua para la elección del iterante
buscamos un intervalo donde la función cambie de
signo, [0,10]

x=0;eval(‘sin(x)-2*cos(2*x)+x^2-pi^2+2’)
x=10;eval(‘sin(x)-2*cos(2*x)+x^2-pi^2+2’)
 ANÁLISIS NUMÉRICO_____________________________

En el intervalo [0,10] hay un cambio de signo esto

implica que existe al menos un cero elegiremos
x0=5
x=fzero(‘sin(x)-2*cos(2*x)+x^2-pi^2+2’,5)
eval(‘sin(x)-2*cos(2*x)+x^2-pi^2+2’)
Programación en Computadora con MatLab – J. Orbegoso
 ANÁLISIS NUMÉRICO_____________________________
Métodos de integración numérica:
Podemos calcular el área bajo curvas F(x) por
integración numérica.
.-Integración de funciones unidimensionales
quad
quadl
.-Integración doble
dblquad
Ejemplo4:
Programación en Computadora con MatLab – J. Orbegoso

ia=quad(‘sin(x)+1’,0,2*pi); %(6.2832)
xmin=pi;xmax=2*pi;
ymin=0;ymax=pi;
result=dblquad('y*sin(x)
+x*cos(y)',xmin,xmax,ymin,ymax)
%-9.8698
 ANÁLISIS NUMÉRICO_____________________________

int('sin(x)^2*cos(x)^2',0,4*pi) %1/2*pi=1.5708
ia=quad('sin(x).^2.*cos(x).^2',0,4*pi) %
1.0051e-030
ial=quadl('sin(x).^2.*cos(x).^2',0,4*pi) %1.5708
Programación en Computadora con MatLab – J.
Orbegoso
 ANÁLISIS NUMÉRICO_____________________________

Métodos de resolución de ecuaciones

diferenciales
Resolución de problemas de valores iniciales para
ecuaciones diferenciales ordinarias (ODEs)
[T,Y]=solver(‘F’,tspan,y0)
.- solver algoritmo de resolución de ODEs, ode45,
Programación en Computadora con MatLab – J.

ode23, ode113, ode15s,ode23s.

.-F string conteniendo el nombre del fichero ODE.
.-tspan vector de tiempos [t0 tfinal] de
integración.
Orbegoso

.-y0 vector columna de condiciones iniciales en t0

 ANÁLISIS NUMÉRICO_____________________________

Solvers ode45 Ec. Dif. No rígidas. Orden

ode23 medio
ode113 Ec. Dif. No rígidas. Orden
bajo
ode15s
Ec. Dif. No rígidas. Orden
ode23s
variab
Ec. Dif. rígidas. Orden variab
Ec. Dif. rígidas. Orden bajo
Programación en Computadora con MatLab – J.

Opcione odeset Crear/modificar opciones

s odeget Obtener opciones
Salidas odeplot Dibujar series temporales
odephas Dibujar fases
2 bidimensionales
Orbegoso

odephas Dibujar fases

 ANÁLISIS NUMÉRICO_____________________________

Ejemplo5:
La ecuación de van der Pol

 
y 1   1  y y 1  y 1 0
2
1

 0
Reescribimos el sistema
Programación en Computadora con MatLab – J.

y 1  y 2
 
y 2   1  y 12 y 2  y 1
Escribimos el fichero ODE
Orbegoso
 ANÁLISIS NUMÉRICO_____________________________
function dy=vdp1(t,y)
dy=[y(2); (1-y(1)^2)*y(2)-y(1)];
Llamamos a solver
[T,Y]=ode45(‘vdp1’,[0 20],[2;0]);
plot(T,Y(:,1),'-',T,Y(:,2),'--')
title('Solucion de la ecuacion de Van der Pol, \
Programación en Computadora con MatLab – J. Orbegoso

mu=1')
xlabel('Tiempo T')
ylabel('Solucion Y')
legend('Y1','Y2')
Programación en Computadora con MatLab – J. Orbegoso

ANÁLISIS NUMÉRICO_____________________________
 AYUDA
La ayuda es lo mas importante de matlab,
help nombre_comando
help nombre_toolbox
Algunas importantes son
help graph2d
help graph3d
Programación en Computadora con MatLab – J. Orbegoso

help specgraph
help save
Podemos ver ejemplos hechos con matlab
poniendo
demo
Para saber más: [Link]
[Link]
Esta última hace referencia a los compiladores de c de matlab, los mex files