UNIVERSIDAD CATOLICA DE SANTA MARIA

FACULTAD DE CIENCIAS E INGENIERIASFISICASY FORMALES

PROGRAMA PROFECIONAL DE INGENIERIA MECANICA ELECTRICA
Y MECATRONICA
MECANICA COMPUTACIONAL II
INFORME DE FASE N1
NOMBRE
CODIGO
SECCI
ON
MEDINA VILLEGAS ARNULFO 200820168
A
ANDRE
1
SALINAS BARREDA EDISON 200880236
A
ERICK
1
ARENAS OVIEDO ALVARO
200820391
A
ALONSO
1
MONZN ARU DIEGO YAIR
200820156
A
1

ING. JUAN CARLOS CUADROS

AREQUIPA-PER 2009-10-8

TRABAJO DE FASE N1
1. Si se drena el agua desde un tanque cilndrico vertical por medio de abrir
una vlvula en la base, el liquido fluir rpidamente cuando el tanque este
lleno y despacio conforme se drene. La tasa a la que el nivel del agua
disminuye es:
dy
=k y
dt
donde es una constante que depende de la forma del agujero y del rea de
la seccin transversal del tanque y agujero de drenaje. La profundidad del
agua y se mide en metros y el tiempo t en minutos. Si k=0.06 determine
cunto tiempo se requiere para vaciar el tanque si el nivel del fluido se
encuentra en un inicio a 3 m.
2.1.1. Resuelva analticamente con por la metodologa de Euler. Utilice un paso de 0.5
minutos
2.1.2. Resuelva analticamente por la metodologa de RK2, bajo las mismas condiciones.
2.1.3. Haga una grafica de los resultados de cada uno de los mtodos aplicados

Resolucin de 2.1.1 y 2.1.2 adjuntado en el archivo de Excel

2.1.3
Grafica de Euler
3

2.5

1.5

0.5

-0.5
0

10

20

30

40

50

60

Grafica de RK2
3

2.5

1.5

0.5

10

20

30

40

50

60

VALOR REAL DE LA INTEGRACION:

0

1
dy
= dt

0.06 3 y 0
1
( 2 3 )=t
0.06
t=57,7

2. Implemente un programa en MATLAB que solucione el problema anterior. El

programa deber solicitar cual de los dos mtodos se emplear en la solucin, y
una vez hecha la eleccin por parte del usuario deber mostrar la respuesta al
problema en forma tabular y grfica. Adjuntar su diagrama de flujo.
Cdigo:
clc, clear all
q=menu('seleccione el metodo','metodo de Euler','Metodo de RK2');
switch q
case 1
F=inline('-0.06*sqrt(y)');
h=0.5;
y(1)=3;
t(1)=0;
i=1;
while y(i)>0
y(i+1)=y(i)+F(y(i))*h;
t(i+1)=t(i)+h;

fprintf('%5.5f %5.5f %5.5f %5.5f\n',t(i),y(i),t(i+1),y(i+1))

i=i+1;

end
grid on
plot(t,y)
grid on
case 2
F=inline('-0.06*sqrt(y)');
h=0.5;
y(1)=3;
t(1)=0;
i=1;
while y(i)>0
k1=F(y(i));
y1(i+1)=y(i)+F(y(i))*h;
if y1(i+1)<0
break
end
t(i+1)=t(i)+h;
k2=F(y1(i+1));
yn(i+1)=y(i)+(h/2)*(k1+k2);
fprintf('%5.5f %5.5f %5.5f %5.5f %5.5f %5.5f
%5.5f\n',t(i),y(i),k1,t(i+1),y1(i+1),k2,yn(i+1))
y(i+1)=yn(i+1);
i=i+1;
end
grid on
plot(t,y)
grid on
end

Algoritmo:
i. Elegir el mtodo a realizar
ii.
Si se elige el mtodo de Euler
i. Ingresar la funcin, h, y0, t0, i=1
ii. Mientras y(i)>0
1. y(i+1)=y(i)+f(y(i))*h
2. t(i+1)=t(i)+h;
3. i=i+1
iii.
Si se elige el mtodo de RK2
i. Ingresar la funcin, h, y0, t0, i=1
ii. Mientras y(i)>0
1. K1=f(y(i))
2. y(i+1)=y(i)+F(y(i))*h
3. t(i+1)=t(i)+h;
4. k2=F(y1(i+1));
5. yn(i+1)=y(i)+(h/2)*(k1+k2)
6. y(i+1)=yn(i+1)
7. i=i+1

y(i)>0
y ( i+1 ) = y ( i ) +

dy
( y ( i ) )h
dt

t ( i +1 ) =t ( i ) +h
i=i+1
t(i), y(i), t(i+1), y(i+1)
y(i)>0
dy
K1= ( y ( i ) )
dt
y1 (i+ 1 )= y (i )+

dy
( y ( i ) )h
dt

t ( i +1 ) =t ( i ) +h
K2=

dy
( y ( i+1 ) )
dt

h
y n (i+ 1 )= y (i )+ ( K1+ K2 )
2
y ( i+1 ) = y n (i+ 1 )
i=i+1
t(i),y(i),K1,t(i+1),y1(i+1),k2,yn(i+1)
INICION
h=0.5;y(1)=3;t(1)=0;i=1;q

q=1
q=2
FIN
Diagrama de flujo:

3. Es frecuente que en los anlisis avanzados de ingeniera surjan funciones de Bessel, como en
el estudio de los campos elctricos. Dichas funciones por lo general no son susceptibles de
evaluarse en forma directa y, por ello, no es raro que estn compiladas en tablas matemticas
estndar. Por ejemplo,

x
1.8
2
2.2
2.4
2.6

J1(x)
0.5815
0.5767
0.556
0.5202
0.4708

Estime J1(2.1) con el uso de un polinomio de interpolacin de Lagrange de grado mximo de

acuerdo a los datos proporcionados. Si el valor verdadero es 0.568292, haga el clculo del error
absoluto y el error relativo

x=2.1

f ( x n )=0.568292

Ea= ?
Er= ?
F ( x 2,1 )=

x2
x2.2
x2.4
x2.6
( 1.82
)( 1.82.2
)( 1.82.4
)( 1.82.6
) ( 0.5815) +

x1.8 x 2.2 x2.4 x2.6

( 21.8
)( 22.2 )( 22.4 )( 22.6 )( 0.5767 )+
x1.8
x2
x2.4
x2.6
( 2.21.8
)( 2.22
)( 2.22.4
)( 2.22.6
) ( 0.5560 ) +
x1.8
x2
x2.2
x2.6
( 2.41.8
)( 2.42
)( 2.42.2
)( 2.42.6
)( 0.5202 )+
x1.8
x2
x2.2
x2.4
( 2.61.8
)( 2.62
)( 2.62.2
)( 2.62.4
)( 0.4708 ) +
x=2.1
2.12 2.12.2 2.12.4
F ( x 2,1 )=
1.82 1.82.2 1.82.4

)(

)(

( 0.5815 ) +
)( 2.12.6
1.82.6 )

2.12.2 2.12.4 2.12.6

( 0.5767 )+
( 2.11.8
)(
21.8
22.2 )( 22.4 )( 22.6 )
2.12 2.12.4 2.12.6
( 0.5560 ) +
( 2.11.8
)(
2.21.8 2.22 )( 2.22.4 )( 2.22.6 )
2.12 2.12.2 2.12.6
( 0.5202 )+
( 2.11.8
2.41.8 )( 2.42 )( 2.42.2 )( 2.42.6 )

2.12 2.12.2 2.12.4

( 0.4708 ) +
( 2.11.8
)(
2.61.8 2.62 )( 2.62.2 )( 2.62.4 )
F ( x 2,1 )=0,571147

Ea=

F ( x 2,1 ) F (x n )
100
F x 2,1

0,5711470,568292
100
0,571147

Ea=

E a =0,499871

x
1.8
2
2.2
2.4
2.6

Er=0.0182( x1.8)(x 2)(x2.2)(x2.4)( x2.6)

F(x)
0.5815
0.5767
0.556
0.5202
0.4708

-0.024
-0.1035
-0.179
-0.247

Er=8.190010

-0.19875
-0.18875
-0.017

4. ados los datos de la siguiente tabla:

x
1
2
3
5
7
8

f(x)
3
6
19
99
291
444

0.016
0.03125

0.0182

1. Calcule f(4) con el uso de polinomios de interpolacin de Newton de rdenes de 1 a

4. Elija los puntos base para obtener una buena exactitud.

2. Estime el error para cada prediccin segn la ecuacin de error Rn

F(4)=???
x=4

a) Grado 1 :
F ( 4 )=F ( x 0 ) + F ( x 1 , x 0 ) ( xx 0 )

x
3
5

F(x)
19
99

40

F ( 4 )=19+ 40(43)
F ( 4 )=59

Error (Rn) :
x
3
5
7

F(x)
19
99
291

40
96

Rn =14( xx 0 )( x x 1 )

Rn =14=14
b) Grado 2 :
F ( 4 )=F ( x 0 ) + F ( x 1 , x 0 ) ( x x 0 ) + F ( x 2, x 1 , x 0 ) ( x x 0 ) (x x 1)

x
2
3
5

F(x)
6
19
99

13
40

14

F ( 4 )=6+13 ( 42 )+ 9 ( 42 ) ( 43)
F ( 4 )=50

x
2
3
5
7

F(x)
6
19
99
291

Error (Rn) :
13
40
96

9
14

Rn =1( xx 0 )( x x 1)( xx 2)

Rn =2=2
c) Grado 3 :
F ( 4 )=F ( x 0 ) + F ( x 1 , x 0 ) ( x x 0 ) + F ( x 2, x 1 , x 0 ) ( x x 0 )( xx 1 ) + F ( x 3 , x 2 x 1 , x 0 )

( x x 0 ) ( x x1 ) ( xx 2)
x
1
2
3
5

F(x)
3
6
19
99

3
13
40

5
9

F ( 4 )=3+3 ( 41 ) +5 ( 41 )( 42 ) +1 ( 41 ) (42)( 43)

F ( 4 )=48

x
1
2
3
5
7

F(x)
3
6
19
99
291

Error (Rn) :
3
13
40
96

5
9
14

1
1

Rn =0
d) Grado 4 :
x 4 , x 3 , x 2 x 1 , x0
F ( 4 )=F ( x 0 ) + F ( x 1 , x 0 ) ( x x 0 ) + F ( x 2, x 1 , x 0 ) ( x x 0 )( xx 1 ) + F ( x 3 , x 2 x 1 , x 0 ) ( x x 0 ) ( x x 1) ( xx 2 ) + F
)

x
1
2
3
5
7

( x x 0 ) ( x x1 ) ( xx 2 ) ( xx 3 )
F(x)
3
6
19
99
291

3
13
40
96

5
9
14

1
1

F ( 4 )=3+3 ( 41 ) +5 ( 41 )( 42 ) +1 ( 41 ) (42) ( 43 ) +0( 41)(42)(43)(45)

F ( 4 )=48

x
1
2
3
5
7
8

F(x)
3
6
19
99
291
444

Error (Rn) :
3
13
40
96
153

Rn =0

5
9
14
19

1
1
1

0
0

5. Suponga que esta diseando un tanque esfrico para almacenar agua para un poblado
pequeo del pas. El volumen de lquido que puede contener el tanque se calcula con:

V = h 2

(3Rh)
3

donde V=volumen(m3), h=profundidad del agua en el tanque (m), y R=radio del tanque (m)
Resuelva por el mtodo de la Falsa Posicin hasta que el error relativo se menor o igual que 0.5e4.
Intervalo: [1:3]

a
1.00000
1.89773
2.01679
2.02617
2.02685
2.02690
2.02691

f(a)
21.62242
-3.21471
-0.25564
-0.01871
-0.00136
-0.00010
-0.00001

f(b)

xr

f(xr)

3.0000
3.0000
3.0000
3.0000
3.0000
3.0000
3.0000

26.54867
26.54867
26.54867
26.54867
26.54867
26.54867
26.54867

1.89773
2.01679
2.02617
2.02685
2.0269
2.02691
2.02691

-3.21571
-0.25564
-0.01871
-0.00136
-0.0001
-0.00001
0

100
5.90317
0.4628
0.03385
0.00246
0.00018
0.00001

Resuelva por el mtodo de Newton Raphson de 2 Orden, hasta que el error relativo se menor o
igual que 0.5e-4.

xr
1.9
2.02682
2.02691
2.02691

e
100.0000
0
6.25715
0.00417
0.00000

6. Implemente un programa en MATLAB que solucione el problema anterior. El programa deber

solicitar cual de los dos mtodos se emplear en la solucin, y una vez hecha la eleccin por
parte del usuario deber mostrar la respuesta al problema en forma tabular y grfica. Adjuntar
su diagrama de flujo.

Cdigo:
clc,clear all
q=menu('elija el metodo','Falsa Posicion','Newton Raphson 2do orden');
switch q
case 1
a=1;
b=3;
e=0.5*10^-4;
x0=0;
er=100;
F=inline('(pi*x^2*(9-x)/3)-30');
while F(a)*F(b)>0
a=input('ingrese un nuevo primer valor: ');
b=input('ingrese un nuevo segundo valor: ');
end
while e<er
xr=a-((F(a)*(b-a))/(F(b)-F(a)));
er=abs((xr-x0)/xr)*100;
x0=xr;
fprintf('%5.5f %5.5f %5.5f %5.5f %5.5f %5.5f
%5.5f\n',a,F(a),b,F(b),xr,F(xr),er)
if F(xr)*F(a)<0
b=xr;
elseif F(xr)*F(a)>0
a=xr;
elseif F(xr)*F(a)==0
break
end
end
case 2
F=inline('(pi*x^2*(9-x)/3)-30');
DF=inline('pi*(6*x-x^2)');
D2F=inline('pi*(6-2*x)');
er=100;
e=0.5*10^-4;
con=inf;
while con>1
x0=1.9;
con=abs(F(x0)*D2F(x0)/(DF(x0)^2));
end
while e<er
x1=x0-(DF(x0)/D2F(x0))+((sqrt((DF(x0)^2)(2*D2F(x0)*F(x0))))/D2F(x0));
x2=x0-(DF(x0)/D2F(x0))-((sqrt((DF(x0)^2)(2*D2F(x0)*F(x0))))/D2F(x0));
e1=abs((x1-x0)/x1)*100;
e2=abs((x2-x0)/x2)*100;
if e1<e2
x0=x1;
er=e1;
else
x0=x2;
er=e2;
end
fprintf('%5.5f %5.5f\n',x0,er)
end
end

Algoritmo:

Algoritmo de la falsa posicin:

i.
ii.

Ingresar f(x), a, b, es, x0=0, er=100

Comparar
1.
Si f(a)*f(b)<0
a.
Mientras es<er
f ( a )( ba )
x r =a
i.
f ( b ) f ( a )

er =

ii.
iii.

x r x 0
100
xr

Comparar
1. Si f(xr)*f(a)<0
a.
b=xr
2. Si f(xr)*f(a)>0
a.
a=xr
3. Si f(xt)*f(a)=0
a.
Xr=0
Si f(a)*f(b)>0
Volver a pedir un intervalo

2.
a.

Algoritmo de Newton Raphson de segundo orden

i. Ingresar f(x), f(x), f(x), es, con=inf, er=inf

ii. Mientras con>1
1. Ingresar x0
'
''
f ( x 0 ) f ( x 0 )
'
2
2. con=
f (x )

iii. Mientras es<er

1.

x 1= x0

f ' ( x0 )
''

f ( x0 )
f ' ( x0 )

2.

x 1= x0

3.

e 1=

x r x 0
100
xr

4.

e 2=

xr x0
100
xr

f ' ' ( x0 )

5. Comparar
a. Si e1<e2
i.
x0=x1
ii.
er=e1
b. Sino
i.
x0=x2
ii.
er=e2

'

2
( x 0 ) 2 f '' ( x 0 ) f ( x 0 )
''
f ( x0 )
'

( x 0 ) 2 f '' ( x 0 ) f ( x 0 )
f ' ' ( x0 )

Diagrama de flujo
Diagrama de la falsa posicin:

INICIO
f(x), a ,b, es, er=100, x0=0
f(a)*f(b)>0
a, b
es<er
f ( a )( ba )
x r =a
f ( b ) f ( a )

er =
f(a)*f(xr)>0
f(a)*f(xr)<0
A, f(a), b, f(b), xr, f(xr) ,er

a=xr
b=xr
xr=0
FIN

x r x 0
100
xr

Diagrama de flujo de Newton Raphson de Segundo orden

Inicio
f(x), f(x), f(x), x0, es, er=100, con=inf

con>1
X0

con=

'

''

f ( x 0 ) f (x 0 )
'

f ( x0 )

es<er
'

x 1= x0

x 1= x0

f ( x0 )
''

f ( x0 )
f ' ( x0 )
''

f ( x0 )

e 1=

x r x 0
100
xr

e 2=

e1<e2
er=e1
x0=x1
er=e2
x0=x2

FIN

xr x0
100
xr

'

2
( x 0 ) 2 f '' ( x 0 ) f ( x 0 )
''
f ( x0 )

'

( x 0 ) 2 f '' ( x 0 ) f ( x 0 )
''
f ( x0 )

Diagrama del programa de seleccin de Mtodo

INICIO

q
q=1
Realizar los procedimientos de la falsa posicin
Realizar los procedimientos de Newthon-Rhanposon de segundo orden
q=2
Fin

7. Suponga que la fuerza hacia arriba de la resistencia del aire sobre un objeto que cae es
proporcional al cuadrado de la velocidad. Para este caso, la velocidad se calcula con:
donde =coeficiente de arrastre de segundo orden.

v ( t )=

gm
tanh
Ca

( gCam t )

9.8 m/ s2 , m=68.1 kg, y ca=0.25

Si

Kg
m , use integracin analtica para determinar qu

tan lejos cae el objeto en 10 segundos.

Haga lo mismo, pero evalu la integral empleando la regla de Simpson 1/3 compuesta. Pruebe
con diferentes hasta obtener tres dgitos significativos de exactitud.

Integracin Analtica:
10

9.868.1
tanh
0.25

t ) dt =333.9262
( 9.80.25
68.1

Integracin por Simpson 1/3

Cuarta derivada
F4=67228/2318805*(1tanh(7/1362*1362^(1/2)*t)^2)^2*tanh(7/1362*1362^(1/2)*t)*1362^(1/2)33614/2318805*tanh(7/1362*1362^(1/2)*t)^3*(1tanh(7/1362*1362^(1/2)*t)^2)*1362^(1/2)
M4=max x [0,10 ] f 4

Para valor mximo de la cuarta derivada

t=2.2
F4 =0.2732
a=0

b=10
ba 4
h M4
180

0.0005

ba
2n

n=10

h=

x0
x1
x2
x3
x4
x5
x6
x7
x8
x9
x10

10

S = f ( t ) dt
0

0.00000
1.00000
2.00000
3.00000
4.00000
5.00000
6.00000
7.00000
8.00000
9.00000
10.00000

2n1

h
f ( x0 )+ 4 f ( xi ) +2
3
i =1

10

S = f ( x ) dx
0

h 0.4356

f(x0)
f(x1)
f(x2)
f(x3)
f(x4)
f(x5)
f(x6)
f(x7)
f(x8)
f(x9)
f(x10)

2n 2

i=2

f ( xi ) + f (x 2n )

h 1

0
9.6841
18.711
26.5902
33.0832
38.1846
42.0446
44.883
46.9266
48.3755
49.3918

1
[ 0+ 4 ( 9,1841+ 26,5002+38.1846+ 44,8830+48,3755 ) + 2 ( 18,711+33,0832+ 42.0446+46,
3
S334,064

CONCLUSIONES:

Como se ve en los programas realizados, el comando syms es de

mucha utilidad al momento de realizar mtodos de interpolacin, porque
permiten trabajar con las funciones creadas con una gran libertad ya que, a
diferencia del comando inline, se define al inicio una variable
El comando inline tiene un mejor uso en los programas y
mtodos que requieren una fcil evaluacin, ya que define la
funcin no solo en una variable, sino en varias

Los clculos de errores relativos son tiles si se quiere encontrar el

grado de erro de un resultado, sin tener el verdadero resultado en el cual
basarse y as, dar una idea de cunto es el error del mtodo

En los mtodos de interpolacin, como se puede apreciar en esta

prctica, mientras haiga ms valores, el resultado ser ms exacto al
resultado real.

Como se pudo observar en el los mtodos para encontrar races, el

mtodo de Newton de segundo orden lo encuentra ms rpido que la falsa
posicin, esto se debe, a que Newton de segundo orden trabaja con la
segunda derivada.