8.

8 El volumen V de un líquido contenido en un tanque horizontal cilíndrico de radio r y longitud L

está relacionado con la profundidad del líquido h por

(
V = r 2∗acos ( r −hr )−(r −h∗√( 2∗r∗h−h )))∗L
2

Determine h para r = 2 m, L = 5 m y V = 8.5 m3

 Método grafico

>> syms h

>> r=2;L=5;V=8.5;

>> V=(r^2*acos((r-h)/r)-(r-h*sqrt(2*r*h-h^2 )))*L

V=

20*acos(1 - h/2) + 5*h*(- h^2 + 4*h)^(1/2) - 10

>> ezplot(V)

>> ezplot(V),grid on
>> ezplot(V,[0,2]),grid on

>> % el valor de h cuando V es 8.5 se encuentra aproximadamente entre 0.2 y 0.6

 Bisección
o Valor inicial: 0.2
o Valor final. 0.6
o Iteraciones: 18
o Valor obtenido para h=0.53015
 Falsa posición
o Valor inicial: 0.2
o Valor final: 0.6
o Iteraciones: 2
o Valor obtenido para h= 0.53041198
 Punto fijo
o f=@(x) 20*acos(1 - x./2) + 5.*x.*(- x.^2 + 4.*x)^(1/2) - 10-8.5;
o g=@(x) (20*acos(1 - x./2) - 10-8.5)./( -5.*(- x.^2 + 4.*x)^(1/2));
o Valor inicial: 0.2, 0.6
o No converge
 Newton-Raphson
o Valor inicial: 0.2
o Iteraciones: 4
o Valor obtenido para h=0.5301496
 Secante
o Valor inicial: 0.2
o Valor final: 0.6
o Iteraciones: 5
o Valor obtenido para h= 0.5301496
 Müller
o Valor x0: 0.6
o Valor x1: 0.7
o Valor x2: 0.8
o Iteraciones: 5
o Valor obtenido de h= 5301496
 Bairstow
o Polinomio de Taylor de 5to grado: (40*pi)/3 + 15*3^(1/2) + (20*3^(1/2)*(x - 3))/3 -
(35*3^(1/2)*(x - 3)^2)/9 - (40*3^(1/2)*(x - 3)^3)/27 - (55*3^(1/2)*(x - 3)^4)/81 -
(104*3^(1/2)*(x - 3)^5)/243 - 37/2
o r=1, s=1
o 5 iteraciones
o r_nueva= 0.3605, s_nueva= -0.0547
2
o Raíces del polinomio de la forma x +rx + s
 -0.47552963441456642388708508116322
 0.11502963441456642388708508116322
o No coincide con ningún valor anterior
 Polinomios de Lagrange
o Innecesario para este problema, seria buscar n puntos en la gráfica evaluados, lo
cual lo hace poco útil para este problema, aun asi, utilicé los puntos para
x_r=0:0.2:1 y los evalué en la función del problema.
 o La única raíz real que salió fue 0.419082847130 que está mas cercana a la
respuesta real.

8.20 La concentración de bacterias contaminantes c en un lago disminuye de acuerdo con la

ecuación
c=75 e – 1.5 t +20 e – 0.075t

Determine el tiempo que se requiere para que la concentración de bacterias se reduzca a 15

 Método grafico
o >> c = 75*exp(-1.5*x) + 20*exp(-0.075*x)-15;
o >> ezplot(c),grid on
o

o >> ezplot(c,[3.9,4.1,-.5,.5]),grid on
 o

 Bisección
o Valor inicial: 3.9
o Valor final. 4.1
o Iteraciones: 15
o Valor obtenido para t= 4.00163
 Falsa posición
o Valor inicial: 3.9
o Valor final. 4.1
o Iteraciones: 2
o Valor obtenido para t= 4.00166791
 Punto fijo
o f=@(x) 75*exp(-1.5*x) + 20*exp(-0.075*x)-15
o g=@(x) 75*exp(-1.5*x) + 20*exp(-0.075*x)-15+x;
o Valor inicial: 3.9
o Iteraciones: 14
o Valor obtenido para t= 4.00163
 Newton-Raphson
o Valor inicial: 3.9
o Iteraciones: 4
o Valor obtenido para t= 4.0016341
 Secante
 o Valor inicial: 3.9
o Valor final: 4.1
o Iteraciones: 4
o Valor obtenido para t= 4.0016341
 Müller
o Diferentes iteraciones antes y después del valor ya obtenido con otros métodos
o i xf(i) Error aprox (i)
o
o 1 6.5131490 34.747
o 2 16.1420753 59.651
o 3 30.2699318 46.673
o 4 38.3482900 21.066
o 5 32.1048687 58.882
o 6 99.6873678 65.887
o 7 175.7506574 41.396
o 8 217.8521509 23.596
o 9 2093.3762036 89.034
o 10 202357.0822452 98.938
o 11 5161102.4888822 99.894
o 12 87421854181236481000000000000000000000000.0000000 100.000
o i xf(i) Error aprox (i)
o
o 1 1.4655716 166.108
o 2 5.0149231 70.776
o 3 21.9952300 77.200
o 4 24.9591010 11.875
o 5 17.9276842 80.818
o 6 69.5199615 69.671
o 7 113.1748818 37.742
o 8 134.6562431 18.229
o 9 498.1554887 75.748
o 10 5424.7979446 94.474
o 11 168965.3954863 96.931
o 12 1942101191594.6128000 100.000
o 13
183035133988030650000000000000000000000000000000000000000000000000
00000000000000000000000000000000.0000000 100.000
o Diverge…
 Bairstow
o Polinomio de Taylor de 5to grado: - 4.7460941455078125*x^5 +
15.8203388671875*x^4 - 42.18890625*x^3 + 84.43125*x^2 - 114.0*x + 80.0
o r=3.9, s=3.9
o 9 iteraciones
o r_nueva= 1.6326, s_nueva= -0.8316
2
o Raíces del polinomio de la forma x +rx + s
 o -2.0402059155016777677507667806902
o 0.40760591550167776775076678069025
o Se toma el mayor, t=0.40760591550167776775076678069025
 Polinomios de Lagrange
o X_r=3.8:.1:4.2.
o Las raíces reales son 3.9619617641 y 292.177117029, de las cuales, la primera
coincide con las anteriores.

7.20 En el análisis de sistemas de control, se desarrollan funciones de transferencia que relacionan

en forma matemática la dinámica de la entrada de un sistema con su salida. La función de
transferencia para un sistema de posicionamiento robotizado está dada por:
3 2
C( s) s + 12.5∗s +50.5∗s +66
G ( s )= = 4
N ( s ) s +19∗s 3 +122∗s2 +296∗s+192

(s 3 +12.5∗s 2+50.5∗s+66)/( s 4 +19∗s 3+ 122∗s2 +296∗s +192)

donde G(s) = ganancia del sistema, C(s) = salida del sistema, N(s) = entrada del sistema y s =
frecuencia compleja de la transformada de Laplace. Utilice una técnica numérica para obtener las
raíces del numerador y el denominador, y factorícelas en la forma siguiente:

( s +a1 ) ( s+a 2) ( s+ a3 )
G ( s )=
( s +b1 ) ( s+b 2) ( s+ b3 ) ( s +b 4 )

donde ai y bi = las raíces del numerador y el denominador, respectivamente.

PARA C (s )
 Método gráfico
o >> syms x
o >> cs=(x^3+12.5*x^2+50.5*x+66)
o >> ezplot(cs), grid on
 o >> ezplot(cs,[-6,-2.5]), grid on

 Bisección
o Valor inicial: -6
o Valor final. -5
o Iteraciones: 1
o Valor obtenido para s1=-5.5
o Valor inicial: -4.5
o Valor final. -3.5
o Iteraciones: 1
o Valor obtenido para s2=-4
o Valor inicial: -3.5
o Valor final. -2.5
o Iteraciones: 1
o Valor obtenido para s3=-3
 Falsa posición
o Valor inicial: -6
o Valor final. -5
o Iteraciones: 8
o Valor obtenido para s1= -5. 499671224443083042418097461154
o Valor inicial: -4.5
o Valor final. -3.5
o Iteraciones: 4
 o Valor obtenido para s2=-3.9999999999953034635181291186821
o Valor inicial: -3.4
o Valor final. -2.6
o Iteraciones: 2
o Valor obtenido para s3=-2.6266258607498087222647283856159
 Punto fijo
o f=@(x) x^3 + (25*x^2)/2 + (101*x)/2 + 66;
o g=@(x) -(x^3 + 12.5*x^2 + 66.0)/50.5;
o Valor inicial {x1,x2,x3}= {-6,-4,-3.5}
o Iteraciones para {x1,x2,x3}= {93,1,163}
%el valor inicial x2 si le ponemos -3.9 converge a -3, y si le ponemos -4.1 converge a -5.5
o Valor obtenido para {x1,x2,x3}={-5.50025,-4, -3.00038}
 Newton-Raphson
o Valor inicial: {-6,-4.5,-2.5}
o Iteraciones: {5,4,5}
o Valor obtenido={-5.5,-4, -3}
 Secante
o Valor inicial: {-6,-4.5,-3.5}
o Valor final: {-5,-3.5,-2.5}
o Iteraciones: {9,6,9}
o Valor obtenido={-5.5, -4, -3.0}
 Müller
o Valor x0: {-6,-4.5,-3.95}
o Valor x1: {-5.9,-4.1,-3.9}
o Valor x2: {-5.8,-4.001,-3.7}
o Iteraciones: {7,3,6}
o Valor obtenido={-5.5,- 4.0010230,-3}
 Bairstow
o x^3 + (25*x^2)/2 + (101*x)/2 + 66
o r=1, s=-3
o 11 iteraciones
o {-3.0000, -4.0000, -5.5000}
 Polinomios de Lagrange
o X_r=[-6:2]
o F(x-r)= [-3,1,0,0,7,27,66,130,225]
o La única raíz real que salió fue -5.7884263870400484748292301684442

PARA N (s)
 Método gráfico
o >> syms x
o >> ns=x^4+19*x^3+122*x^2+296*x+192
o >> ezplot(ns), grid on
 o >> ezplot(ns,[-10,0]), grid on

 Bisección
o Valor inicial: {-9,-7,-5,-2}
o Valor final. {-7,-5,-3,0}
o Iteraciones: {1,1,1,1}
o Valor obtenido={-8,-6,-4,-1}
 Falsa posición
o Valor inicial: {-9,-7,-5,-2}
o Valor final. {-7,-5,-3,0}
o Iteraciones: {11,3,2,3}
o Valor obtenido={-7.997412, -5.99894, -3.1731025, 0.5773195876}
 Punto fijo
o f=@(x) x^4+19*x^3+122*x^2+296*x+192;
o g=@(x) -(x^4+19*x^3+122*x^2+192)/296;
o Valor inicial {x1,x2,x3,x4}= {-9,-7,-5,-2}
o Iteraciones para {x1,x2,x3,x4}= {indef, 142,149,31}
o Valor obtenido para {x1,x2,x3,x4}={divergente, -6.00005, -5.99995, -1.00001}
 Newton-Raphson
o Valor inicial: {-9,-6.5,-4.5,-2}
o Iteraciones: {5,4,5,11}
o Valor obtenido={ -8.0000000,-6.00000, -4.000, -1.00}
 Secante
 o Valor inicial: {-8.5,-6.5,-5,-2}
o Valor final. {-7.5,-5.5,-3,0}
o Iteraciones: {9,6,7,9}
o Valor obtenido={-8, -6, -4,-1}
 Müller
o Valor x0: {-6,-4.5,-3.95}
o Valor x1: {-5.9,-4.1,-3.9}
o Valor x2: {-5.8,-4.001,-3.7}
o Iteraciones: {7,3,6}
o Valor obtenido={-5.5,- 4.0010230,-3}
 Bairstow
o x^3 + (25*x^2)/2 + (101*x)/2 + 66
o r=1, s=-3
o 11 iteraciones
o {-3.0000, -4.0000, -5.5000}
 Polinomios de Lagrange
o X_r=[-6:2]
o F(x-r)= [-3,1,0,0,7,27,66,130,225]
o La única raíz real que salió fue -5.7884263870400484748292301684442

Problema resuelto:

( s+ 5.5 )( s+ 4 ) ( s +3 )
G ( s )=
( s+ 8 ) ( s+ 6 ) ( s+ 4 )( s+1 )

18.22 Emplee la porción de la tabla de vapor que se da para el H2O supercalentada a 200 MPa,
para a) encontrar la entropía correspondiente s para un volumen específico v de 0.108 m3 /kg con
interpolación lineal.

valores del vector x_r: [0.10377,0.11144,0.1254]

valores del vector f(x_r): [6.4147,6.5453,6.7664]

valor de evaluación: 0.108

La funcion de polinomios f=(547231660*X)/1659021 - (327265*((100000*X)/767 - 10377/767)*(X -

627/5000))/698 + (641470*((100000*X)/767 - 11144/767)*(X - 627/5000))/2163 +
(169160*((100000*X)/2163 - 3459/721)*(X - 1393/12500))/349 - 30762738233/1106014000

el valor de la funcion evaluada es f(0.108)=14.2975

ingrese el grado del polinomio, n=5

Se necesitan 6 puntos

ingrese los puntos

x0= 1

y0= 3

x1= 2

y1= 6

x2= 3

y2= 19

x3= 5

y3= 99

x4= 7

y4= 291

x5= 8

y5= 444

La matriz de diferencias divididas es:

3 0 0 0 0 0

6 3 0 0 0 0

19 13 5 0 0 0

99 40 9 1 0 0

291 96 14 1 0 0

444 153 19 1 0 0
El polinomio de newton es

3 2

x -x -x+4

ingrese el valor de x a interpolar,x=4

el valor interpolado es 48.00