ANÁLISIS NUMÉRICO PARA

INGENIEROS CIVILES CI 708

2022-02
CONTENIDO
1. MÉTODOS NUMÉRICOS EN LA INGENIERÍA

2. MATLAB. TIPOS DE VARIABLES

3. MANEJO DE VECTORES. POLINOMIOS Y RAÍCES

4. MANEJO DE MATRICES. PROPIEDADES

5. REPRESENTACIÓN GRÁFICA EN MATLAB: 2D y 3D

6. CONCEPTOS BÁSICOS DE LA PROGRAMACIÓN: ALGORITMOS

7. CREACIÓN DE FUNCIONES Y APLICACIONES CON MATLAB

8. PROGRAMACIÓN EN MATLAB: IF, FOR , WHILE

9.INTRODUCCIÓN AL APPDESIGNER

10.TEORÍA DE ERRORES
MATLAB: Definición

Es un entorno computacional para realizar cálculos

numéricos, visualización gráfica, simulaciones, rutinas de
programación.
Es para la solución de problemas prácticos de ingeniería y
matemáticas.
MATLAB: Entorno

Ventana
de Comandos

Workspace

Ventana
de Historial
Carpeta actual
MATLAB: Tipos de Variables

Numéricos:
INTEGER (8, 16, 32, 64 bits)
SINGLE (32 bits)
DOUBLE (64 bits)
Lógico: valor “0” o “1”
Complejos: complex(a,b)
Caracteres: string
Matrices: 2, 3 o más dimensiones
MATLAB: Operadores

Operadores Lógicos Operadores Relacionales

OR MAYOR
QUE

AND IGUAL
QUE

MENOR
NOT
QUE
EJERCICIOS

>> 7 + 8/2
ans = 11
>> 27^(1/3)+32^0.2
ans = 5
>> (7+8)/2
ans = 7.5000
>> 27^1/3+32^0.2
ans =11
>> 4+5/3+1
ans = 6.6667
>>0.7854-(0.7854)^3/(1*2*3)+0.785^5/(1*2*3*4*5)...
-(0.785)^7/(1*2*3*4*5*6*7)
>> 5^3/2
ans = 0.7071
ans = 62.5000
MATLAB: Manejo de Vectores

>> V1 = [ 1 3/11 -10.90 sqrt(7) ]

V1 = 1.0000 0.2727 -10.9000 2.6458

>> V2= 0:1:5

V2 = 0 1 2 3 4 5

>> V3 = linspace( 0 , 1 , 5)
V3= 0 0.2500 0.5000 0.7500 1.0000

>> V1( 1 , 2 ) = pi;

V1 = 1.0000 3.1416 -10.9000 2.6458

>> V1(4) >> V2 = V1’;

ans = 2.6458
MATLAB: Polinomios y raices

𝑷 𝒙 = 𝟒𝒙𝟑 − 𝟐𝒙𝟐 + 𝟎. 𝟒𝒙 − 𝟏𝟎𝟎

Q 𝒙 = 𝒙𝟒 − 𝟕𝒙𝟑 + 𝟐. 𝟗𝒙 + 𝟏𝟎. 𝟖
MATLAB: Manejo de Matrices

>> M1 = [ 1 2 4 ; 5 0.1 -11 ]

M1 =
1.0000 2.0000 4.0000
5.0000 0.1000 -11.0000

>> M2 = ones( 3 , 3 );
>> M3 = zeros( 6 , 6 );
>> M4 = rand( 4 , 5 );

>> M4( 2 , 5 ) = 10;

>> M4( 1 : 4 , 2 ) = 1 + zeros( 4 , 1 );
MATLAB: Manejo de Matrices

2 −4 3.3
Z = 10 0.01 −0.4
3 2 −6

%Extraer todos los datos de la fila 2

z(2,:)

%Extraer todos los datos de la columna 3

z(:,3)

%Si queremos capturar el ultimo valor "-6"

z(end,end)

% Para extraer los dos primeros números de la fila 1

z(1,1:2)

% Para extraer el primer elemento y el ultimo de la columna 3

z([1 3],3)
MATLAB: Propiedades de las Matrices

Identidad: Producto de matrices

>> T1 = eye(5); >> T5 = T4 * T3;
Traspuesta: Producto vectorial
>> T1‘ >> T5 = T4 .* T3;
Inversa: División:
>> T2 = inv( T1 ); >> T6 = rand( 5 , 5 );
Suma y Resta: >> T7 = 4 ./ T6;
>> T3 = T1 + T2; Potencia:
Producto con escalar: >> T8 = T7.^2;
>> T4 = 3 * T1; >> T9 = power( T8,2/3); 2/3 );
MATLAB: Gráficos 2D

z=rand()
z=rand(5)
z=rand(5,1)
z=10*rand(5,1)
z=ceil(10*rand(5,1))
MATLAB: Gráficos 2D

• plot(x,y)
plot(x,y,'kp:')
• plot(x,y,S))

• plot(x1,y1,S1,x2,y2,S2,x3,y3,S3,...)
MATLAB: Gráficos 2D

2
𝑓 𝑥 = ⅇ 𝒔𝒆𝒏(π𝑥 3
−𝑥
൯ en el intervalo [-3,3]

PASO 1: x=linspace(-3,3,400);
PASO 2: y=exp(-x.^2).*sin(pi*x.^3);
PASO 3: plot(x,y);
MATLAB: Gráficos 2D
clear, clc
x=[ 0 2 5 8 ];
y = [ 4.1 6.6 1 8.9 ];
plot(x,y)
figure(2)
plot(x,y,'r*')
figure(3)
plot(x,y,'kp:')
figure(4)
plot(x,y)
hold on
plot(x,y,'r*')
MATLAB: Gráficos 2D

−𝑥 2 −5𝑥 2
𝑓 𝑥 = ⅇ , g(x)=ⅇ en el intervalo [-2,2]

(a) Con trazo continuo.

(b) La gráfica de 𝑓 𝑥 con asteriscos rojos y la de g 𝑥 con
cruces verdes

x=-2:0.1:2; (a)
y1=exp(-x.^2);
y2=exp(-5*x.^2);
plot(x,y1,x,y2)
figure(2)
plot(x,y1,’r*’,x,y2,’g+’) (b)
MATLAB: Gráficos 2D
X2 si X <0;
F(X) 1 si 0 <= X<1;
-X+2 si 1<=X

x= linspace(-2,3,1000);
y=(x.^2).*(x<0)+ 1*((0<=x)&(x<1))…
+(-x+2).*(1<=x);
plot(x,y,'.g');
grid on;
title('Funcion con discontinuidad');
MATLAB: Gráficos 2D

t = linspace( 0 , 1 , 200 );
x = sin(2*pi*3*t);
subplot(3 1 1), plot( t , x )
subplot(3 1 2), stem( t , x )
subplot(3 1 3), bar( t , x )
grid minor
xlabel(' Tiempo ')
MATLAB: Gráficos 2D
Dibujar las gráficas de las funciones y = sen(x), b = ln(a) en el intervalo [-5,5].
Utilizar la orden grid para estimar un intervalo en el que se encuentre el punto
de intersección. Con la orden axis para ver con mayor claridad la parte de la
figura en la que se encuentra el punto de intersección. Haciendo uso del zoom
comprobar que las coordenadas del punto de intersección son (2.218, 0.796)

x=-5:.1:5;
y=sin(x);
plot(x,y)
hold on
a=eps:0.1:5;
b=log(a);
plot(a,b)
grid on
axis([0,3,-2,2])
zoom on
MATLAB: Gráficos 3D

Arbitrario:
x=[ 0 2 5 8 ];
y = [ 4.1 6.6 1 8.9 ];
z = [ 0.1 -0.9 3 12 ];
plot3( x , y , z , 'b' , 'linewidth',3)
grid minor
xlabel('Eje X')
ylabel('Eje Y')
zlabel('Eje Z')
MATLAB: Gráficos 3D

𝑟Ԧ = (sen 𝑡 , cos 𝑡 , 𝑡 ) 0 ≤ 𝑡 ≤ 8𝜋

t=linspace(0,8*pi,200);
y1=sin(t);
y2=cos(t);
y3=t;
plot3(y1,y2,y3);
figure(2)
comet3(y1,y2,y3);
Curvas en el espacio
MATLAB: Gráficos 3D

x = 0 : 0.1 : 3 ; 𝒙𝟐 − 𝒚𝟐
𝒇 𝒙, 𝒚 =
𝒔𝒆𝒏 𝟑𝟎𝒙 + 𝟑
y=x;
[ X ,Y ] = meshgrid( x , y );
Z=(X.^ 2 - Y.^ 2)./(sin(30*X) + 3);
surf( X , Y , Z )
xlabel('Eje X')
ylabel('Eje Y')
zlabel('Eje Z’)

Superficies
MATLAB: Gráficos 3D

𝑥 = cos 𝑢 sen 𝑣 𝑢, 𝑣 ∈ 0.2𝜋 𝑥[0, 𝜋ሿ

𝑦 = 𝑠ⅇ𝑛 𝑢 sen 𝑣
z = cos 𝑢

u=linspace(0,2*pi,32);
v=linspace(0,pi,32);
[U V]=meshgrid(u,v);
X=cos(U).*sin(V);
Y=sin(U).*sin(V);
Z=cos(V);
surf(X,Y,Z) Ecuaciones Paramétricas