USO DE SOFTWARE ESPECIALIZADO

CÁLCULO NUMÉRICO
INTRODUCCIÓN

IMPORTANCIA DEL MATLAB, OCTAVE Y SCILAB

COMANDOS BÁSICOS Y PRINCIPALES

MÉTODOS NUMÉRICOS
Ing. Patricio Pugarín Díaz, Mgs.
DEPARTAMENTO DE CIENCIAS EXACTAS - ESPE
 Semana 2

CONTENIDO

Título Programación para los métodos numérico

Duración 240 minutos

Información general Resaltar la importancia de la programación de los Métodos

Numéricos para resolver problemas de la Ingeniería.

Objetivo Conocer, aplicar los comandos básicos y de programación de

Matlab, Octave y Scilab.

1
 Semana 2

SUBMATRICES
Trabajemos con la matriz mágica A de orden (6x6)

Se quiere extraer de la matriz A la submatriz A1, que corresponden a las filas 3, 4 y todas las
columnas.

Se puede generar de 3 formas:

FORMA 1

FORMA 2

FORMA 3

2
 Semana 2

Matrices por Bloques

Las siguientes matrices se pueden considerar Bloques

Con estas matrices se puede formar La matriz por Bloques A

Ejercicio
Resolver el sistema lineal de ecuaciones
5𝑥1 − 2𝑥2 + 𝑥3 = 4
{2𝑥1 + 4𝑥2 − 3𝑥3 = −6
𝑥1 − 𝑥2 + 2𝑥3 = 1

Generamos las matrices A de coeficientes del sistema y B de términos independientes, en

Matlab y resolvemos el sistema:

3
 Semana 2

Ejercicio

9 8 7
Dada la matriz 𝑥 = (4 5 6). En no más de 2 líneas de ejecución generar la matriz
7 8 9
9 8 7
𝑦 = (8 5 8)
7 8 9

Ejercicio.-
1. Dada una matriz 𝑀 cuadrada aleatoria uniforme de orden 3, obtener su inversa, su
transpuesta y su diagonal. Transformarla en una matriz triangular inferior y en otra
superior y rotarla 90 grados. Obtener la suma de los elementos de la primera fila y la
suma de los elementos de la diagonal. Extraer la submatriz cuya diagonal son los
elementos 𝑎11 y 𝑎22 y extraer también la submatriz cuyos elementos de la diagonal
son 𝑎11 y 𝑎33 .

𝑖 2𝑖 3𝑖
2. Dada la matriz 𝑀 = (4𝑖 5𝑖 6𝑖 ), obtener su logaritmo neperiano elemento a
7𝑖 8𝑖 9𝑖
elemento y realice las operaciones matrices 𝑒 𝑀 y 𝐿𝑛(𝑀).

4
 Semana 2

GRÁFICAS DE FUNCIONES EN 𝑹𝟐
Matlab permite representar funciones de una y dos variables en:
 Coordenadas cartesianas
 Coordenadas polares
 Ecuaciones paramétricas

Gráfica de la función 𝑦 = 𝑠𝑒𝑛 𝑥 en el intervalo [0, 2pi]

FORMA 1

FORMA 2

5
 Semana 2

Añadir títulos y textos

axis → Permite cambiar la escala del eje x y del eje y
xlabel → Poner título en el eje x
ylabel → Poner título en el eje y
title → Coloca un título al gráfico
grid → Produce un mallado en el plano
text → Coloca un texto en una coordenada dada
En la pantalla gráfica de Matlab, que aparece cuando se hace un gráfico, tenemos numerosos
comandos que permiten modificar el aspecto de la gráfica, salvarla, imprimirla

Varias funciones en la misma gráfica

6
 Semana 2

GRÁFICAS DE FUNCIONES EN 𝑹𝟑

7
 Semana 2

Otros comandos que se pueden utilizar:

>> surf(z)
>> surfc(z)
>> meshc(z)
>> meshz(z)

Ejercicio de clase
Construir la gráfica de la función 𝑦 = 𝑒 𝑥 . 𝑠𝑒𝑛(𝑥). Adicionalmente, intuir el valor del cero
negativo de la función.

8
 Semana 2

Una primera intención sería:

Una segunda intención sería:

NOTA.-
Se puede seguir aproximando cada vez más

9
 Semana 2

Ejercicio de clase
Realizar el gráfico de la función
𝑺𝒆𝒏(√𝒙𝟐 + 𝒚𝟐 )
𝒛=
√𝒙𝟐 + 𝒚𝟐

10
Semana 2

11
 Semana 2

PROGRAMACIÓN

TIPOS DE 𝒎. 𝒇𝒊𝒍𝒆𝒔

1. ARCHIVOS DE INSTRUCCIONES (script).- Secuencia de instrucciones dentro de

𝑚. 𝑓𝑖𝑙𝑒𝑠 que persiguen una solución específica y única.

2. ARCHIVO DE FUNCIONES (𝒇𝒖𝒏𝒄𝒕𝒊𝒐𝒏).- Conjunto de instrucciones que se

pueden ejecutar como una función intrínseca de Matlab. Son nuevas funciones
definidas por el usuario.

Ejercicio
Hacer un programa que calcule el cuadrado del número 3

12
Semana 2

13
 Semana 2

𝑨𝒓𝒄𝒉𝒊𝒗𝒐 𝒅𝒆 𝒇𝒖𝒏𝒄𝒊ó𝒏

14
Semana 2

15
Semana 2

16
 Semana 2

ESTRUCTURAS DE CONTROL
Bucles for.- Repite un conjunto de instrucciones un determinado número de veces. Su
estructura es
𝒇𝒐𝒓 𝒊𝒏𝒅𝒆𝒙 = 𝒏1: 𝒏2
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
𝒆𝒏𝒅

If condición
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
Else
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
end

Bucle while.- Son similares a los bucles for; pero estos continúan hasta que se satisface algún
criterio. Su estructura es:

while criterio
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
end

Comando break.- Se usa para terminar un bloque prematuramente (sale del bucle). Su
estructura es:
if condición
𝑀𝑒𝑛𝑠𝑎𝑗𝑒 𝑑𝑒 𝑠𝑎𝑙𝑖𝑑𝑎
break
end

17
 Semana 2

Comado continue.- Similar al comando break. En lugar de salir del bucle, el programa salta
al siguiente paso. Su estructura es:

if condición
𝑀𝑒𝑛𝑠𝑎𝑗𝑒 𝑑𝑒 Continuar
continue
end

Comando switch case.- Se utiliza cuando existe una serie de rutas de programación, para una
variable dada, dependiendo de su valor. Permite elegir entre múltiples salidas. Su estructura
es:

N=input(‘ opción A=1, opción B=2, opción C=3’)

switch N
case 1
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
case 2
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
case 3
𝐿𝑖𝑛𝑒𝑎𝑠 𝑑𝑒 𝑒𝑗𝑒𝑐𝑢𝑐𝑖ó𝑛
otherwise
disp(‘no existe este caso’)
end

COMANDOS ADICIONALES IMPORTANTES

fix(x) Elimina la parte decimal de x
rem(x,y) Calcula el resto de dividir x para y
round(x) Toma la parte entera de x (aproximando)
floor(x) Toma el valor entero más próximo a la izquierda del número (no aproxima).
ceil(x) Toma el valor entero más próximo a la derecha del número (no aproxima).
sign(x) Determina el valor del signo de x. (-1 si es negativo, 1 si es positivo, 0 si x=0)

18
 Semana 2

Ejercicios para la clase. -

1. Calcular la suma de los elementos menores a 90 e indicar el número de elementos

mayores o iguales a 90, en el siguiente vector:
𝑣 = (60, 80, 95,38, 90, 96, 48, 102)
2. Construir un archivo de funciones que proporcione una matriz 𝐴 = (𝑎𝑖𝑗 ) de (mxn), m,
n son parámetros de entrada; siendo
1
𝑎𝑖𝑗 = .
𝑖+𝑗−1
Corresponde ésta a la matriz de Hilbert, mal condicionada.
3. Escribir la palabra ‘HOLA’, 10 veces, aplicando solo el bucle while.
4. Hacer un programa que calcule el logaritmo natural de un número positivo 𝑎,
ingresado por pantalla; 10 ingresos en total (Usar comando break).
5. Hacer un programa que calcule el logaritmo natural de un número positivo 𝑎,
ingresado por pantalla; 10 ingresos en total (Usar comando continue).

REFERENCIAS BIBLIOGRAFICAS

1. Sánchez Juan Miguel, Problemas de Cálculo Numérico para ingenieros con aplicaciones
Matlab, McGraw-Hill, Primera edición, 2005.
2. A. Quarteroni, F. Saleri, Cálculo Científico con Matlab y Octave. Springer-Verlag
Italia, milano 2006.
3. César Pérez López, MATLAB a través de ejemplos. IBERGARCETA
PUBLICACIONES, S.L., Madrid 2011.

19