PAULA ANDREA MURCIA DIAZ
VALENTINA VARGAS ROCHA
TALLER 2 METODOS NUMERICOS
4. Por un canal trapezoidal fluye agua tratada procedente de una mina a
una tasa de Q = 20 m3/s. La profundidad crtica y para dicho canal
satisface la ecuacin
Donde g = 9.81m/s2, Ac = rea de la seccin transversal (m2), y B =
ancho del canal en la superficie (m). Para este caso, el ancho y el rea
de la seccin transversal se relacionan con la profundidad y por medio
de:
Resuelva para la profundidad crtica con el uso de los mtodos a)
biseccin, y b) falsa posicin. En los incisos a), Haga elecciones iniciales
de xl = 0.5 y xu = 2.5, y ejecute iteraciones hasta que el error
aproximado caiga por debajo del 1% o el nmero de interacciones
supere a 10.
Ecuacion:
400
f ( y) 1 (3 y )
9.81(3 y y 2 / 2) 3
GRAFICA
Grafica 1. Grafica de la funcin
INTERSECCION
Grafica 2. Corte en el eje x
METODO DE LA FALSA POSICION
Muestra de Clculos
Xi=0,5 y xu=0,25
f ( xu)( xi xu)
x r xu
f ( xi) f ( xu)
0.812841(0.5 2.5)
xr 2.5 2.450889
32.292154 0.812841
f (0.5) f (2.450889) 32.292154(0.799686) 25.823592
Como -25,823592 cumple con el teorema de Bolzano Si f(xl)f(xr) < 0,
entonces la raz se encuentra dentro del subintervalo inferior o
izquierdo. Por lo tanto, haga xu = xr
Entonces xU= 2.4508989 y la segunda iteracin sera:
0,799686(0.5 2.450889)
x r 2.450889 2.403748
32.292152 0.799686
%ERROR
2.403748 2.450889
a 100% 1.961%
2.403748
f (0.5) f (2.403748) 32.292152(0.785941) 25.379737
Nuevamente al ser el nmero anterior menor cero, xu = 2.403748; y
como se muestra en la siguiente tabla, xu es siempre xr
Tabla 1. Mtodo falsa posicin
IT XI Xr Xu F(XI) F(Xr) F(XI)*F(X ERROR
r) %
1 0,50000 2.45089 2500000 - 0.79968 -
00 31 0 32.292152 63 25.82359
2
2 0,50000 2.40374 2.45089 - 0.78594 - 1.961
00 85 31 32.292152 13 25.37973
7
3 0,50000 2.35851 2.40374 - 0.77161 - 1.918
00 51 85 32.292152 64 24.91715
3
4 0,50000 2.31514 2.35851 - 0.75672 - 1.873
00 25 51 32.292152 55 24.43629
5
5 0,50000 2.27358 2.31514 - 0.74128 - 1.828
00 09 25 32.292152 73 23.93776
2
6 0,50000 2.23378 2.27358 - 0.72532 -23.42229 1.782
00 08 09 32.292152 45
7 0,50000 2.19569 2.23378 - 0.70886 - 1.735
00 33 08 32.292152 44 22.89075
7
8 0,50000 2.15926 2.19569 - 0.69193 - 1.687
00 97 33 32.292152 84 22.34417
9
9 0,50000 2.12446 2.15926 - 0.67458 - 1.638
00 16 97 32.292152 19 21.78370
1
10 0,50000 2.09122 2.12446 - 0.65683 - 1.590
00 11 16 32.292152 44 21.21059
5
13 0,50000 1.51462 1.51463 - 0.00021 -0.007028 0
7 00 80 53 32.292152 76
COMPROBACION
clear
clc
funcprot(0)
function y=f(x)
y=1-((400)/((9.8)*((3*x)+(x^2/2))^3))*(3+x);
endfunction
function xr=falsaposicion(xai, xbi, tol)
i=1;
ea(1)=100;
if f(xai)*f(xbi) < 0
xi(1)=xai;
xu(1)=xbi;
xr(1)=xi(1)-f(xi(1))*(xu(1)-xi(1))/(f(xu(1))-f(xi(1)));
printf('It.\t\t Xi\t\t Xr\t\t Xu\t Error \n');
printf('%2d \t %11.7f \t %11.7f \t %11.7f \n',i,xi(i),xr(i),xu(i));
while abs(ea(i))>=tol,
if f(xi(i))*f(xr(i))< 0
xi(i+1)=xi(i);
xu(i+1)=xr(i);
end
if f(xi(i))*f(xr(i))> 0
xi(1)=xr(i);
xu(1)=xu(i);
end
xr(i+1)=xi(i+1)-((f(xi(i+1))*(xu(i+1)-xi(i+1)))/(f(xu(i+1))-f(xi(i+1))));
ea(i+1)=abs((xr(i+1)-xr(i))/(xr(i+1)))*100;
printf('%2d \t %11.7f \t %11.7f \t %11.7f \t %7.3f \n',
i+1,xi(i+1),xr(i+1),xu(i+1),ea(i+1));
i=i+1;
end
else
printf('No existe una raiz en ese intervalo');
end
endfunction
MTODO DE LA BISECCIN
IT xi f(xi) xu f(xu) xr f(xr) ERROR
2.500000 1.500000
1 0.5 -322.921 0.8128 -0.03199
0 0
2.500000 2.000000
2 1.5 -0.03199 0.8128 0.60140 25.000
0 0
2.000000 1.750000
3 1.5 -0.03199 0.6014 0.3782 14.286
0 0
1.750000 1.625000
4 1.5 -0.03199 0.3782 0.2061 7.692
0 0
1.625000 1.562500
5 1.5 -0.03199 0.2061 0.09703 4.000
0 0
1.562500 1.531250
6 1.5 -0.03199 0.09703 0.03527 2.041
0 0
1.531250 1.515625
7 1.5 -0.03199 0.03527 0.0023 1.031
0 0
1.515625 1.507812
8 1.5 -0.03199 0.0023 -0.01462
0 5 0.518
Se puede concluir que los tres mtodos mostraron resultados
satisfactorios con respecto a la raz ya que todos concuerdan; adems,
para este caso en especfico, el mtodo de la falsa posicin requiri
muchas ms iteraciones para poder llegar a la raz (1.51
aproximadamente). El numero obtenido que corresponde al corte de la
grfica con el eje x es igual a la profundidad critica la cual satisface la
ecuacin dada.
