Métodos Numéricos
Métodos Numéricos
Métodos Numéricos
2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
Problema 3
Al golpear una pelota de tenis con una raqueta, la pelota cambia de forma: primero se comprime y
después se expande. Sea la compresión de la pelota, con , y sea la fuerza
ejercida por la raqueta sobre la pelota. La energía trasferida es proporcional al área bajo la curva
. Supongamos que es la fuerza durante la compresión y la fuerza durante la
expansión. La energía es transferida a la pelota durante la compresión y transferida por la pelota
durante la expansión, de modo que la energía perdida por la pelota en la colisión (debido al
rozamiento) es proporcional a . Por lo tanto, el porcentaje de energía perdida en la
Universidad Nacional de Ingeniería P.A. 2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
colisión viene dada por: Se han efectuado medidas en varios golpes de tenis con
los resultados de la tabla.
(libras) 0 25 50 90 160
(libras) 0 23 46 78 160
a) (2.0 ptos) Aplicar el método de Simpson (1/3) compuesto, para estimar el porcentaje de
energía perdida por la pelota en la colisión con la raqueta. Comente su respuesta.
b) (2.0 ptos) Usando los mismos puntos de a), estimar el porcentaje de energía perdida por la
pelota por la colisión con la raqueta utilizando el método del trapecio compuesto.
c) (1.0 pto) Implementar un script en MATLAB que resuelva a).
Problema 4
Sea el sistema masa-resorte mostrado en la figura:
El cual es gobernado por el siguiente conjunto de ecuaciones
diferenciales:
d 2 y1 dy
m1 2
c1 1 k1 y1 0
dt dt
dy dy
2
c2 2 1 k 2 y2 y1 0
d y2
m2
dt dt
2
dt
Se sabe que c1=4, c2=1, k1=2, k2=4, unidades en el Sistema Internacional.
Si al aplicar un paso de Euler se obtuvo los siguientes valores:
T y1 y1’ y2 y2’
0 1 1 2 -1
Los Profesores
Universidad Nacional de Ingeniería P.A. 2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
Solucionario
Pregunta 1:
<1
Por lo tanto converge.
d)
A=[ 8 9
4 6]
b=[2 10]'
D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1);
T=inv(D)*(L+U); c=inv(D-L)*b;
x=[0 0 0]';
xa=x;
for i=1:300000
x=T*x+c;
er=norm(x-xa,2);
xa=x;
if er<10e-14
break
end
end
disp(x)
disp(i)
Universidad Nacional de Ingeniería P.A. 2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
Solucionario
Pregunta 1:
<1
Por lo tanto converge.
d)
A=[ 8 9
4 6]
b=[2 10]'
D=diag(diag(A)); L=-tril(A,-1); U=-triu(A,1);
T=inv(D)*(L+U); c=inv(D-L)*b;
x=[0 0 0]';
xa=x;
for i=1:300000
x=T*x+c;
er=norm(x-xa,2);
xa=x;
if er<10e-14
break
end
end
disp(x)
disp(i)
Universidad Nacional de Ingeniería P.A. 2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
Pregunta 2:
Universidad Nacional de Ingeniería P.A. 2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
Pregunta 3:
a)
Regla de Simpson:
b)
0, 4
f c ( x)dx 24,5
0
0,1
I1 f c (0) 2 f c (0,1) 2 f c (0,2) 2 f c (0,3) f c (0,4) 24,5
2
0, 4
f e ( x)dx 24,5
0
0,1
I 2 f e (0) 2 f e (0,1) 2 f e (0,2) 2 f e (0,3) f e (0,4) 22,7
2
c)
fcpar=fc(2:2:length(x)-1);
fcimpar=fc(3:2:length(x)-1);
fepar=fe(2:2:length(x)-1);
feimpar=fe(3:2:length(x)-1);
I=100*(fc(1)-fe(1)+4*sum(fcpar-fepar)+2*sum(fcimpar-feimpar)+(fc(length(x))-
fe(length(x))))./(fc(1)+4*sum(fcpar)+2*sum(fcimpar)+fc(length(x)))
Pregunta 4
dy1
v1
dt
dv1 cv k y
1 1 1 1
dt m1
dy2
v2
dt
dv2 c v v k 2 y2 y1
2 2 1
dt m2
b) Algoritmo de Euler:
y1N 1 y1N h * v1
N
c1v1 k1 y1
N N
v1N 1 v1N h *
m1
y2N 1 y2N h * v2
N
v N 1
2 v h*
N
2
c2 v2 v1
N N
k y
2 2
N
y1
N
m2
c) Algoritmo de Taylor 2:
h2
y1N 1 y1N h * v1
N
* v1 ' N
2
c v k1 y1 c v ' N k1 y1 ' N
N N
h2
v1N 1 v1N h * 1 1 * 1 1
m1 2 m1
h2
y2N 1 y2N h * v2
N
* v2 ' N
2
v2N 1 v2 h *
N
N
N N
c2 v2 v1 k 2 y2 y1
N
h2
* 2 2
c v ' N v1 ' N k 2 y2 ' N y1 ' N
m2 2 m2
Universidad Nacional de Ingeniería P.A. 2015-1
Facultad de Ingeniería Mecánica 21/07/15
DACIBAHCC
t y1 y1’ y2 y2’
0 1 1 2 -1
e) Comando MATLAB:
>> [y,z]=dsolve('D2y=-2*Dy-y','D2z=-1/5*(Dz-Dy)-4/5*(z-y)','y(0)=1','Dy(0)=1','z(0)=2','Dz(0)=-1')