MATLAB para Economistas

1
MATLAB para Economistas

• José Luis Hueso
• Matemática Aplicada
• Universidad Politécnica de Valencia

2
Itinerario

• 1ª Etapa Invertir en MATLAB
• 2ª Etapa MATLAB funciona
• 3ª Etapa MATLAB marca la diferencia

3
Contenido

• Mandatos básicos
• Variables
• Constantes
• Funciones
• Vectores
• Gráficas
• Matrices
• Polinomios
• Ecuaciones diferenciales

• Ajuste mínimo cuadrático
• Splines
• La cesta de la compra
• Interés compuesto
• Venta a plazos
• Modelos económicos

4
Invertir en MATLAB

• Mandatos básicos
• Variables
• Vectores
• Gráficas
• Matrices
• Funciones
• Polinomios

• La cesta de la compra
• El IPC
• Interés compuesto
• El IGBM
• Matriz de transición
• La Cuenta Naranja

5
Mandatos básicos

• help, help tema
• dir
• diary fichero
• Comentarios
• Edición de líneas de comando
• Cortar y Pegar

6
Variables

• Asignar
• a 3, b 4
• Listar
• ans
• who
• whos
• Eliminar
• clear b

• Guardar
• save fichero
• Recuperar
• load fichero

7
Constantes y cadenas

• Número p pi
• Unidad imaginaria i,j
• Precisión eps
• Infinito Inf
• Indeterminación NaN
• Cadenas de caracteres
• 'Esto es una cadena'
• Esto no

• Formatos numéricos
• format long
• format short
• format rat
• format bank

8
Vectores

• Edición
• u 11 9 2001
• v 1,2,3
• w 123
• w 123

• Componentes
• u(1), u(2), u(3)
• V. Progresivos
• 1731
• linspace(1,29,5)
• Norma
• norm(v)

9
Operaciones con vectores

• Suma de vectores uv
• Producto por escalar 2u
• Producto escalar dot(u,v)
• Producto matricial uw
• Producto elemental u.v
• Potencia elemental u2, u.v

10
Cesta de la compra

• Productos
• Pan, Leche, Carne, Pescado
• Precio unitario
• p 80, 100, 2000, 1000
• Cantidad
• c 5, 12, 0.750, 2.5
• Importe
• dot(p,c)

11
Más operaciones con vectores

• Suma de componentes
• sum(1 1 1 1 1 1 1)
• Suma acumulada
• cumsum (1 1 1 1 1 1 1)
• Producto de componentes
• prod(1 2 3 4 5 6)
• Producto acumulado
• cumprod()

12
IPC Septiembre 2001

• Intermensual
• Im 0 0.3 0.4 0.5 0.4 0.3 ... 0.2 0.2 0
• Incremento desde principio de año
• Imn 10.01Im
• Ima cumprod(Imn)
• Índice mes a mes (inicio de año 133.4)
• Imm 133.4Ima

13
Interés compuesto

• Tipo de interés anual (1998-2001)
• i 1.020 1.025 1.035 1.030
• Tipo acumulado
• ic cumprod(i)
• Valor de la inversión
• C Coic

14
Gráficas

• Gráfico de barras
• any 19982001
• bar(any,i)
• Títulos
• title('Tipo de interés anual')
• xlabel('Año')
• ylabel('Tipo')

15
Matrices

• Edición
• A 1,23,4
• B -1 -2-3 -4
• Elemento A(2,1)
• Fila A(2,)
• Columna A(,1)

• Bloques
• M A,BB,A
• Submatrices
• M41 M(13,24)
• fil 1,2,4
• col 1,3,4
• M32 M(fil,col)

16
Operaciones con matrices

• Suma y resta -
• Producto .
• Potencia .
• Cociente izq. / ./
• Cociente der. \ .\
• Transpuesta ' .'

• Determinante
• det(A)
• Inversa
• inv(A)
• Rango
• rank(A)

17
Matrices usuales

• Identidad de orden n eye(n)
• Nula de tamaño mn zeros(m,n)
• Matriz de mn unos ones(m,n)
• Matriz mn aleatoria rand(m,n)
• Cuadrado mágico magic(n)

18
Matriz de transición

• Tres plataformas de TV compiten en el mercado
  C, con el 60 Q, con el 20 y D, con el resto.
• De una temporada a otra, C pierde el 10 de sus
  clientes a favor de Q y el 35 a favor de D.
• Al mismo tiempo, el 5 de los clientes de Q
  pasan a C y el 10 a D.
• El 10 de los clientes de D pasan a Q y el 5 a
  C.
• Cómo se distribuye el mercado la temporada
  siguiente?

19
Matriz de transición
20
Matriz de transición
21
Funciones

• sin
• cos
• tan

• help elfun
• sqrt
• exp
• log
• log10
• ezplot (x31)/x

22
Funciones estadísticas

• Densidad y distribución
• normpdf
• normcdf
• logncdf
• Funciones inversas
• norminv

• Estadística descriptiva
• mean
• std
• median
• var

23
Funciones de fecha

• Texto
• d date
• Vector
• c clock
• Conversiones
• datestr(t)
• datestr(t,8)
• fix(datevec(t))
• datenum('1-jan-0000')

• Valor numérico
• t now
• f floor(t)
• h rem(t,1)
• Más
• help timefun

24
Importar datos de una Hoja de Cálculo

• Nombrar el rango a importar datos
• Posición inicial del rango fila, columna
• Guardar el fichero como .wk1 mihoja
• Leer los datos desde MATLAB
• ffila-1 ccolumna-1
• Awk1read('mihoja',f,c,'datos')

25
Exportar una matriz a una Hoja de Cálculo

• Amagic(5)
• wk1write('Cuadradomagico',A,4,2)

Nombre de fichero (.wk1)
Matriz a exportar
Filas y columnas de margen
26
Importar de un fichero ASCII

• load fichero.txt
• Lee filas de datos numéricos separados por
  espacios.
• Admite comentarios precedidos por .
• Genera una variable llamada "fichero".

27

• Gráfico de líneas
• load igbm.txt ascii
• plot(igbm(,2)),hold
• plot(igbm(,2),'ro')
• Títulos
• title('IGBM del 3/9 al 26/10')
• xlabel('Sesión')
• ylabel('Índice')
• gtext('11 de Septiembre')

28
Polinomios Cuenta Naranja

• Interés nominal anual
• r 0.0393
• Interés mensual
• rm i/12
• i (1rm)
• T.A.E
• tae i12-1

29
Polinomios Cuenta Naranja

• Movimientos
• p 1 3 2 4 1 5 0 2... 1 2 3 5 -7
• Saldo final
• polyval(p,i)

30
F I Nde la primera parte