UNIVERSIDAD NACIONAL DE INGENIERIA 2017-2
FACULTAD DE INGENIERIA MECANICA 20/02/2018
DACIBHCC
EXAMEN DE APLAZADOS DE METODOS NUMERICOS (MB536)
SOLO SE PERMITE EL USO DE UNA HOJA DE FORMULARIO Y CALCULADORA CIENTIFICA
ESCRIBA CLARAMENTE SUS PROCEDIMIENTOS
PROHIBIDO EL USO DE CELULARES U OTROS EQUIPOS DE COMUNICACION ELECTRONICA
LAS MOCHILAS O MALETINES DEBERÁN DEJARSE EN LA PARTE INFERIOR DE LA PIZARRA
DURACION: 110 MINUTOS
Problema 1
Sea el sistema lineal, siendo a real:
a) (1.0 P) Para qué valores de a, la matriz es simétrica y definida positiva?
b) (1.0 P) Para qué valores de a, el sistema presenta solución única?
c) (1.0 P) Para qué valores de a, el sistema es convergente para Jacobi, según el criterio de
la diagonalmente estrictamente dominante?
d) (1.0 P) Para qué valores de a, el sistema es convergente para Jacobi, según el criterio
del Radio Espectral?
e) (1.0 P) Realice 2 iteraciones y muestre el error para a=2, partiendo de un vector inicial
nulo.
Problema 2
Sea la función:
a) (1.0 P) Localice las raíces de f, indicando intervalos de longitud unitaria
b) (2.5 P) Determine la menor raíz positiva, partiendo de un valor igual al punto medio del
intervalo obtenido en a) y realice 03 iteraciones de Newton-Raphson y estime el error
c) (1.5 P) Escriba un código Matlab, para resolver b) hasta tener un error menor a 10-12.
Problema 3
Se desea calcular el área de la región que delimitan las gráficas de f y g:
a) (1.0 P) Plantear la integral correspondiente
b) (2.0 P) Resolver mediante la cuadratura de Gauss, con N=3.
c) (0.5 P) Estime el error si la integral exacta es 12.92546881
d) (1.5 P) Escriba un código MATLAB, que lea N, entre 1 y 3 y muestra la Integral de
Gauss aproximada, para este problema.
UNIVERSIDAD NACIONAL DE INGENIERIA 2017-2
FACULTAD DE INGENIERIA MECANICA 20/02/2018
DACIBHCC
Problema 4
Sea el Péndulo de la figura, cuyo movimiento
es gobernado por la siguiente ecuación
diferencial ordinaria:
Considere g=9.8 m/s2 y L=1 m. Si en el
instante inicial el ángulo es 0.75 rad y la
velocidad angular es 0 rad/s:
a) (1.0 P) Reducir a una sistema de 2
EDO de primer orden
b) (2.5 P) Determine la posición y
velocidad angular para t=0.05, 0.1 seg, aplicando Taylor de segundo orden
c) (1.5 P) Escriba un programa que determine el tiempo para que el péndulo este en
posición vertical aplicando Taylor de orden 2.
Los Profesores
UNIVERSIDAD NACIONAL DE INGENIERIA 2017-2
FACULTAD DE INGENIERIA MECANICA 20/02/2018
DACIBHCC
Solucionario Examen de Aplazados de Métodos Numéricos
Problema 1
a) Por Silvester:
Intersectando:
b) Para solución única:
c) Para el criterio de la diagonal estrictamente dominante:
d) Para el criterio del radio espectral de Jacobi:
e) Aplicando Jacobi:
N xN yN Err
0 0 0 -----
1 0.125 0.1111 0.125
2 0.1111 0.1042 0.0139
UNIVERSIDAD NACIONAL DE INGENIERIA 2017-2
FACULTAD DE INGENIERIA MECANICA 20/02/2018
DACIBHCC
Problema 2
a)
Raíces en [-1,0] y [1,2]
b)
xn Error
1.500000000000000 -----------------
2.009894305760567 0.509894305760567
1.927253998159970 0.082640307600597
1.925693186393502 0.001560811766468
c)
% subsana2.m
s='((x-1)/2)^2-sin(x+1)'
ezplot(s) % raices en [-1,0] y [1,2]
grid
ds=diff(s)
f=inline(s)
df=inline(ds)
x=1.5, acum=[x NaN]; TOL=1e-12;
for i=1:10
xn=x-f(x)/df(x);
err=abs(xn-x);
acum=[acum; xn err];
x=xn;
if err<TOL
break
end
end
disp(acum)
UNIVERSIDAD NACIONAL DE INGENIERIA 2017-2
FACULTAD DE INGENIERIA MECANICA 20/02/2018
DACIBHCC
Problema 3
a)
dx
b) Cuadratura de Gauss
dt
I=12.67652495
c)
Err=0.2489
d)
syms t
s='10/(x^2+1)-(x^2-2*x+2)'
ezplot(s)
grid
Ie=double(int(s,-1,2))
ss=subs(s,(3*t+1)/2)
ff=inline(3/2*ss)
N=input('Ingrese N=')
switch N
case 1
I1=2*ff(0)
case 2
I2=ff(sqrt(3)/3)+ff(-sqrt(3)/3)
case 3
I3=5/9*ff(-
sqrt(3/5))+8/9*ff(0)+5/9*ff(sqrt(3/5))
% I3 = 12.676524953789279
otherwise
disp('N fuera de rango!!!')
end
Err=abs(Ie-I3)
UNIVERSIDAD NACIONAL DE INGENIERIA 2017-2
FACULTAD DE INGENIERIA MECANICA 20/02/2018
DACIBHCC
Problema 4
a)
=0.75
b)
Taylor de orden 2
t Ѳ v
0 0.75 0
0.05 0.7416 -0.3340
0.1 0.7167 -0.6620
c)
t(1)=0; u1(1)=0.75; u2(1)=0; h=0.05;
for i=1:20
t(i+1)=t(i)+h;
u1(i+1)=u1(i)+h*u2(i)+h^2/2*-9.8*sin(u1(i));
u2(i+1)=u2(i)+h*-9.8*sin(u1(i))+h^2/2*-9.8*cos(u1(i))*u2(i);
if u1(i+1)<=0
break
end
end
plot(t,u1,t,u2),grid
legend('u1(t)','u2(t)')
disp([t' u1' u2'])
