0% encontró este documento útil (0 votos)
43 vistas41 páginas

Todos Los Ejercicios (Barbieri)

Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como XLS, PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
43 vistas41 páginas

Todos Los Ejercicios (Barbieri)

Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como XLS, PDF, TXT o lee en línea desde Scribd

UTN FRBA

ANALICIS NUMERICO Y CALCULO AVANZADO

Profesores:

Ing. AnalisaPerez
Profesor Ing. Julian Barbieri

REFERENCIA:

Alumno:

Apellidos, Nombres:

Numero de Legajo:
El examen consta de un conjunto de actividades, a resolver, cada una
referenciada en una hoja del presente libro como EJ.1, EJ2, etc.
En estos espacio deben mostrarse la
memoria de calculo, que concluya en la/las soluciones que se requieran para
cada situación. La
hoja rotulada como "Memoria descriptiva" debe contener la redacción que de
respuesta a las consignas de cada actividad. La memoria descriptiva debe
seguir el orden propuesto para cada actividad. Considere a esta parte del
examen como las propuestas que efectúa un profesional de la ingeniería,
respecto de una problemática por la que fue consultado. Sean claros en la
redacción y contundente en los argumentos. La puntuación asignada a cada
actividad se corresponde a una solución adecuada, las desviaciones en este
proceso, o la falta de originalidad, traerán aparejadas una merma en la
puntuación.
El criterio en la asignación, de puntuaciones parciales, es
propio del cuerpo docente. En ninguna oportunidad se asignara puntuación
negativa por mala resolución.
Importante: Todas las soluciones deben ser realizadas usando
modelos de calculo desarrollados en planillas Excel.
Puntuación: EJ.1 .... 4 PUNTOS,
EJ.2.... 3 PUNTOS, EJ.3.... 3 PUNTOS. La condición de destacado requiere
un mínimo de dos EJ completos y correctos.
ejercicio 1

ejercicio 2

A P (0.524)= 0.51089186

B P(3)= 0.1486

C errores 0.01054443

0.00747999
ejercicio 3

f'(4) -1.4952381 no puede dar valor negativo una derivada


El crecimiento de bacterias en un organismo puede ser modelado con una función
exponencial de la forma , donde es el número de bacterias
(en millones) y t es el tiempo medidoN=N 0∗e
en horas
bt
a partir del N
momento en que se
detectaron bacterias (lo cual ocurre a tiempo t=0). Se realizo un experimento
para medir el numero de bacterias a distintos tiempos, el cual proveyo la siguiente
informacion:

t [ hs ] N [ millones ]
5 3
10 3.6
15 4.4
20 5.1
25 7.5

Encontrar N 0 y b aplicando el método de cuadrados mínimos a la ecuación


linealizada
¿Cuántas bacterias se espera encontrar a las 48 horas?
¿En qué momento se hallaran 2 millones de bacterias? ¿Cómo se interpreta el
valor del tiempo encontrado?

N=No∗e^Bt
t [ hs ] N [ millones ]

Ln(N)=ln⁡(No)∗ln⁡(e^Bt )
5 3
10 3.6
15 4.4

Ln(N)=ln⁡(No)+ b*t
20 5.1
25 7.5
V=ln N
A=b
U=x
B=ln ( No )

n x(i) y(i) U(i) V(i)


1 5 3 5 1.09861229
2 10 3.6 10 1.28093385
3 15 4.4 15 1.48160454
4 20 5.1 20 1.62924054
5 25 7.5 25 2.01490302

1375 7.50529424

75
562.897068

(■8(∑_(i=1)^5▒ 〖 Ui 〗 ^2
&∑_(i=1)^5▒Ui@∑_(i=1)^5▒Ui&∑_(i=1)^5▒
1))*(█(A@B))=(█(∑_(i=1)^5▒Ui∗Vi@∑_(i=1)^
5▒Vi))
(■8(∑_(i=1)^5▒ 〖 Ui 〗 ^2
&∑_(i=1)^5▒Ui@∑_(i=1)^5▒Ui&∑_(i=1)^5▒
1))*(█(A@B))=(█(∑_(i=1)^5▒Ui∗Vi@∑_(i=1)^
5▒Vi))

1375 75 A 562.89706765
x =
75 5 B 7.5052942353

FILA 1 1375A + 75B 562.897068


FILA 2 75A + 5B 7.50529424

1375A + 75B 562.897068


15*FILA 2 1125A + 75B 112.579414
F2 - F1 (-)250A + 0 -450.317654
A= 1.80127062
B= 0.31112856

N=1,36∗e^(1,80∗t)
A)
A=b 1.8012706165
A=LN (No)
No=e´B 1.3649646912

B) t 48 horas
N=N 0∗ebt No 1.36496469
b 1.80127062
86.4609896
n 4.83785E+37 3.5443E+37

N 2 MILLONES

C) N/NO 1.465239367
Ln(1,465) 0.3820186195 horas
Bt )
Dado los puntos (x;y) : (0;0),(1,571;1,000), (0,785;0,707),
(2,618;0,500), (3,142;0).
a) Determinar, por alguno de los métodos de interpolación vistos, la
imagen de x = 0,524, según un polinomio cubico.
b) Determinar, por alguno de los métodos de interpolación vistos, la
imagen de x = 3,000, según un polinomio cubico.
c) Considerando que la función generadora de los (x;y) datos es f(x) =
seno(x), determinar los errores cometidos en la determinación de la
imagen para las preguntas a) y b).
X Y
0 0 .--> Lo descarto, para encontrar la imágen de 3, necesito 4 puntos par
0.785 0.707
1.571 1
2.618 0.5
3.142 0 .--> Lo descarto, para encontrar la imágen de 0.524, necesito 4 puntos

Método de Newton: Diferencias divididas

A) Orden X Y 1era Dif. Divi 2da Dif. Divid 3era Dif. Dividida
0 0 0 0
1 0.785 0.707 0.90063694
2 1.571 1 0.37277354 -0.33600471
3 2.618 0.5 -0.47755492 -0.46389987 -0.04885224

a0 a1 a2 a3
0 0.90063694 -0.33600471 -0.04885224

Polinomio cúbico 0.524 3

P (x)= 0.90063694 (x) - 0.3360047 (x)(x-0.785)-0.0488522 (x)(x-0.785)(x-1.571)

CA
0.2637637

0.038349

-0.06024629 P (X)= ´-0,04885224(x^3)-0,22090885(x^2)-1,10415436(x)

0.07674686 P (0.524)= -0.00702875 -0.06065627 0.57857688 .=

1.10415436
-0.22090885
-0.04885224

C) Errores

e(0.524)= 0.51089186 - 0.01054443

e(3)= 0.1489-sen(3 0.00747999


la imágen de 3, necesito 4 puntos para la formación de un polinomio cúbico.

la imágen de 0.524, necesito 4 puntos para la formación de un polinomio cúbico.

B) Orden
3era Dif. Dividida 0
1
2
3

P3 (x)=

x-0.785)(x-1.571) P(3)=

0.51089186
X Y 1ra Dif. Divid 2da Dif. Divid 3era Dif. Dividida
0.785 0.707 0.707
1.571 1 0.37277354
2.618 0.5 -0.47755492 -0.46389987
3.142 0 -0.95419847 -0.30340137 0.0680944

a0 a1 a2 a3
0.707 0.37277354 -0.46389987 0.0680944
.

0.0681(x^3) - 0.8023(x^2) + 1.9692(x)-0.377

0.1486
Para el siguiente conjunto de valores (x;y): (-1;4), (4;8), (5;6); (8;9),
encuentre el valor de la derivada en x=4 usando el metodo de
diferencias centrado de orden 4, con h=1. Es necesario recurrir a los
conceptos vistos durante toda la cursada. Justifique la solucion
adoptada.
x y
-1 4 falta completa
4 8
5 6
8 9

M X Y 1 2 3
0 -1 4
0.8
1 4 8 -0.46666667
-2 0.13518519
2 5 6 0.75
1
3 8 9

F(x) = a0 + a1 (X-X0) +
F(x) = 4 + 0.8 (X-(1)) +

F(6) = 4 + 0.8 7 +
F(6) = 4.95925926

F(3) = 4 + 0.8 4 +
F(3) = 10.1481481

F(2) = 4 + 0.8 3 +
F(2) = 11.6333333

F'(x) = ( -F(x+2h)+8*f(x+h)-8(f(x-h) + F(x-2h) )/ 12h

F'(4) = (-f(6)+8*f(5)-8*f(3) + f(2)/12

F'(4) = -2.20925926
a0 a1 a2 a3
4 0.8 -0.46666667 0.13518519

a2 (X-X0) (X-X1) + a3 (X-X0) (X-X1)


-0.46666667 (X-(1)) (X-4) + 0.13518519 (X-(1)) (X-4)

-0.46666667 7 2 + 0.13518519 7 2

-0.46666667 4 -1 + 0.13518519 4 -1

-0.46666667 3 -2 + 0.13518519 3 -2
(X-X2)
(X-5)

-2

-3
Construir un polinomio interpolador a partir de los siguientes valores (x;y), (1;150), (2;250
determinar la imagen de x=3,5. b) calcular la imagen de x=3,5 con una interpolacion lineal. C
imgenes obtenidas entre los puntos a y b.

X Y
1 150
2 250
3 200
4 350

Método de Newton: Diferencias divididas

A) Orden X Y 1era Dif. Div 2da Dif. Divi 3era Dif. Div.
0 1 150 150
1 2 250 100
2 3 200 -50 -75
3 4 350 150 100 58.3333333

a0 a1 a2 a3
150 100 -75 58.3333333

Polinomio cúbico 3.5

P3 (x)= 150 + 100*(x-1) -75*(x-1)*(x-2) + 58.3333 * (x-1)*(x-2)*(x-3)

P3 (x)= .-450 + 966.63 (x) - 424.98 (x)^2 + 58.33 (x)^3

Para la obtención de la imágen, reemplazo el valor de x = 3.5, en el polinomio obtenido


P3 (3.5)= 228.09875

B) Método de interpolación lineal: Interpolación de Lagrange

P1 (x) = y0*((x-x1)/(x0-x1)) + y1*((x-x0)/(x1-x0))

x0= 3
y0= 200

x1= 4
y1= 350

P1 (x)= 150x - 250

P1(x)= 275
res (x;y), (1;150), (2;250); (3;200); (4;350), a) para
na interpolacion lineal. Cual es la diferencia entre las
puntos a y b.
Considere la ecuación de estado de Van der Waals para un gas real,
dada por:

(P+a/V²)(V-b)=RT
Donde P es la presión medida en atmósferas, T es la temperatura medida
en grados Kelvin, V es el volumen molar del gas medido en lt/gr mol, R es
la constante universal de los gases 0.08205 lt atm /mol K, y a y b son
constantes particulares de cada gas. Calcule V a T=80 C (353 K), para
presiones de 10, 20, 30 y 100 atm, con un error menor a 10-6 , utlizando
el métodos de Bisección

Nota:
a b R T
He 0.0341 0.23700 0.08205 353
Método de Bisección

F(v) = (P+a/V²)(V-b) - R*T = 0

He a b P T R V e
0.0341 0.237 10 353 0.08205 3.1323584318 0.000001
0.0341 0.237 20 353 0.08205 1.6843129396 0.000001
0.0341 0.237 30 353 0.08205 1.2016950846 0.000001
0.0341 0.237 100 353 0.08205 0.5262807608 0.000001

Para P = 10

N n a b c F(a) F(b) Cumple


1 0 -1 6 2.5 -41.375832 28.671808842 Cumple
2 1 2.5 6 4.25 -6.3213031 28.671808842 Cumple
3 2 2.5 4.25 3.375 -6.3213031 11.1739261 Cumple
4 3 2.5 3.375 2.9375 -6.3213031 2.4257441992 Cumple
5 4 2.9375 3.375 3.15625 -1.9479781 2.4257441992 Cumple
6 5 2.9375 3.15625 3.046875 -1.9479781 0.2388427006 Cumple
7 6 3.046875 3.15625 3.1015625 -0.8545788 0.2388427006 Cumple
8 7 3.1015625 3.15625 3.1289063 -0.3078707 0.2388427006 Cumple
9 8 3.1289063 3.15625 3.1425781 -0.0345146 0.2388427006 Cumple
10 9 3.1289063 3.1425781 3.1357422 -0.0345146 0.1021638794 Cumple
11 10 3.1289063 3.1357422 3.1323242 -0.0345146 0.0338245877 Cumple
12 11 3.1323242 3.1357422 3.1340332 -0.000345 0.0338245877 Cumple
13 12 3.1323242 3.1340332 3.1331787 -0.000345 0.0167397773 Cumple
14 13 3.1323242 3.1331787 3.1327515 -0.000345 0.0081973739 Cumple
15 14 3.1323242 3.1327515 3.1325378 -0.000345 0.0039261727 Cumple
16 15 3.1323242 3.1325378 3.132431 -0.000345 0.0017905722 Cumple
17 16 3.1323242 3.132431 3.1323776 -0.000345 0.000722772 Cumple
18 17 3.1323242 3.1323776 3.1323509 -0.000345 0.0001888719 Cumple
19 18 3.1323509 3.1323776 3.1323643 -7.808E-05 0.0001888719 Cumple
20 19 3.1323509 3.1323643 3.1323576 -7.808E-05 5.539691E-05 Cumple
21 20 3.1323576 3.1323643 3.1323609 -1.134E-05 5.539691E-05 Cumple
22 21 3.1323576 3.1323609 3.1323593 -1.134E-05 2.202816E-05 Cumple
23 22 3.1323576 3.1323593 3.1323584 -1.134E-05 5.343778E-06 Cumple

Para P = 20

N n a b c F(a) F(b) Cumple


1 0 -1 6 2.5 -53.746 86.301808842 Cumple
2 1 -1 2.5 0.75 -53.746 16.308696928 Cumple
3 2 0.75 2.5 1.625 -18.673 16.308696928 Cumple
4 3 1.625 2.5 2.0625 -1.1857 16.308696928 Cumple
5 4 1.625 2.0625 1.84375 -1.1857 7.560983503 Cumple
6 5 1.625 1.84375 1.734375 -1.1857 3.187467535 Cumple
7 6 1.625 1.734375 1.6796875 -1.1857 1.0008245765 Cumple
8 7 1.6796875 1.734375 1.7070313 -0.0925 1.0008245765 Cumple
9 8 1.6796875 1.7070313 1.6933594 -0.0925 0.4541777549 Cumple
10 9 1.6796875 1.6933594 1.6865234 -0.0925 0.1808565738 Cumple
11 10 1.6796875 1.6865234 1.6831055 -0.0925 0.0441965526 Cumple
12 11 1.6831055 1.6865234 1.6848145 -0.0241 0.0441965526 Cumple
13 12 1.6831055 1.6848145 1.68396 -0.0241 0.0100316072 Cumple
14 13 1.68396 1.6848145 1.6843872 -0.0071 0.0100316072 Cumple
15 14 1.68396 1.6843872 1.6841736 -0.0071 0.0014903746 Cumple
16 15 1.6841736 1.6843872 1.6842804 -0.0028 0.0014903746 Cumple
17 16 1.6842804 1.6843872 1.6843338 -0.0006 0.0014903746 Cumple
18 17 1.6842804 1.6843338 1.6843071 -0.0006 0.0004227206 Cumple
19 18 1.6843071 1.6843338 1.6843204 -0.0001 0.0004227206 Cumple
1 20 19 1.6843071 1.6843204 1.6843138 -0.0001 0.0001558071 Cumple
21 20 1.6843071 1.6843138 1.6843104 -0.0001 2.23504E-05 Cumple
22 21 1.6843104 1.6843138 1.6843121 -4E-05 2.23504E-05 Cumple
23 22 1.6843121 1.6843138 1.6843129 -1E-05 2.23504E-05 Cumple

Para P = 30

N n a b c F(a) F(b) Cumple


1 0 -1 6 2.5 -66.115832 143.93180884 Cumple
2 1 -1 2.5 0.75 -66.115832 38.938696928 Cumple
3 2 0.75 2.5 1.625 -13.542551 38.938696928 Cumple
4 3 0.75 1.625 1.1875 -13.542551 12.69427409 Cumple
5 4 1.1875 1.625 1.40625 -0.4256653 12.69427409 Cumple
6 5 1.1875 1.40625 1.296875 -0.4256653 6.1340121428 Cumple
7 6 1.1875 1.296875 1.2421875 -0.4256653 2.854088831 Cumple
8 7 1.1875 1.2421875 1.2148438 -0.4256653 1.2141890195 Cumple
9 8 1.1875 1.2148438 1.2011719 -0.4256653 0.3942559731 Cumple
10 9 1.2011719 1.2148438 1.2080078 -0.0157062 0.3942559731 Cumple
11 10 1.2011719 1.2080078 1.2045898 -0.0157062 0.1892745381 Cumple
12 11 1.2011719 1.2045898 1.2028809 -0.0157062 0.0867840991 Cumple
13 12 1.2011719 1.2028809 1.2020264 -0.0157062 0.0355389496 Cumple
14 13 1.2011719 1.2020264 1.2015991 -0.0157062 0.0099163924 Cumple
15 14 1.2015991 1.2020264 1.2018127 -0.0028949 0.0099163924 Cumple
16 15 1.2015991 1.2018127 1.2017059 -0.0028949 0.0035107549 Cumple
17 16 1.2015991 1.2017059 1.2016525 -0.0028949 0.0003079365 Cumple
18 17 1.2016525 1.2017059 1.2016792 -0.0012935 0.0003079365 Cumple
19 18 1.2016792 1.2017059 1.2016926 -0.0004928 0.0003079365 Cumple
20 19 1.2016926 1.2017059 1.2016993 -9.242E-05 0.0003079365 Cumple
21 20 1.2016926 1.2016993 1.2016959 -9.242E-05 0.0001077603 Cumple
22 21 1.2016926 1.2016959 1.2016943 -9.242E-05 7.672233E-06 Cumple
23 22 1.2016943 1.2016959 1.2016951 -4.237E-05 7.672233E-06 Cumple

Para P = 100

N n a b c F(a) F(b) Cumple


1 0 -1 6 2.5 -152.70583 547.34180884 Cumple
2 1 -1 2.5 0.75 -152.70583 197.34869693 Cumple
3 2 -1 0.75 -0.125 -152.70583 22.3674492 Cumple
4 3 -0.125 0.75 0.3125 -65.953679 22.3674492 Cumple
5 4 0.3125 0.75 0.53125 -21.387287 22.3674492 Cumple
6 5 0.3125 0.53125 0.421875 -21.387287 0.4969027308 Cumple
7 6 0.421875 0.53125 0.4765625 -10.440729 0.4969027308 Cumple
8 7 0.4765625 0.53125 0.5039063 -4.9714306 0.4969027308 Cumple
9 8 0.5039063 0.53125 0.5175781 -2.2371812 0.4969027308 Cumple
10 9 0.5175781 0.53125 0.5244141 -0.870122 0.4969027308 Cumple
11 10 0.5244141 0.53125 0.527832 -0.1866057 0.4969027308 Cumple
12 11 0.5244141 0.527832 0.526123 -0.1866057 0.1551494475 Cumple
13 12 0.526123 0.527832 0.5269775 -0.0157279 0.1551494475 Cumple
14 13 0.526123 0.5269775 0.5265503 -0.0157279 0.0697108361 Cumple
15 14 0.526123 0.5265503 0.5263367 -0.0157279 0.0269914859 Cumple
16 15 0.526123 0.5263367 0.5262299 -0.0157279 0.0056317997 Cumple
17 16 0.5262299 0.5263367 0.5262833 -0.005048 0.0056317997 Cumple
18 17 0.5262299 0.5262833 0.5262566 -0.005048 0.0002918769 Cumple
19 18 0.5262566 0.5262833 0.5262699 -0.0023781 0.0002918769 Cumple
20 19 0.5262699 0.5262833 0.5262766 -0.0010431 0.0002918769 Cumple
21 20 0.5262766 0.5262833 0.5262799 -0.0003756 0.0002918769 Cumple
22 21 0.5262799 0.5262833 0.5262816 -4.187E-05 0.0002918769 Cumple
23 22 0.5262799 0.5262816 0.5262808 -4.187E-05 0.0001250043 Cumple
N = Ln(b-a)-ln(error)/ln(2) = 22.738923 Con esto concluimos, que el número de

F(c) Error
-6.3213031 3.5
11.173926 1.75
2.4257442 0.875
-1.9479781 0.4375
0.2388427 0.21875
-0.8545788 0.109375
-0.3078707 0.0546875
-0.0345146 0.0273438
0.1021639 0.0136719
0.0338246 0.0068359
-0.000345 0.003418
0.0167398 0.001709
0.0081974 0.0008545
0.0039262 0.0004272
0.0017906 0.0002136
0.0007228 0.0001068
0.0001889 5.341E-05
-7.808E-05 2.67E-05
5.54E-05 1.335E-05
-1.134E-05 6.676E-06
2.203E-05 3.338E-06
5.344E-06 1.669E-06
-3E-06 8.345E-07

F(c) Error
16.308697 3.5
-18.672551 1.75
-1.1857259 0.875
7.5609835 0.4375
3.1874675 0.21875
1.0008246 0.109375
-0.0924631 0.0546875
0.4541778 0.0273438
0.1808566 0.0136719
0.0441966 0.0068359
-0.0241333 0.003418
0.0100316 0.001709
-0.0070509 0.0008545
0.0014904 0.0004272
-0.0027802 0.0002136
-0.0006449 0.0001068
0.0004227 5.341E-05
-0.0001111 2.67E-05
0.0001558 1.335E-05
2.235E-05 6.676E-06
-4.438E-05 3.338E-06
-1.101E-05 1.669E-06
5.668E-06 8.345E-07

F(c) Error
38.938697 3.5
-13.542551 1.75
12.694274 0.875
-0.4256653 0.4375
6.1340121 0.21875
2.8540888 0.109375
1.214189 0.0546875
0.394256 0.0273438
-0.0157062 0.0136719
0.1892745 0.0068359
0.0867841 0.003418
0.0355389 0.001709
0.0099164 0.0008545
-0.0028949 0.0004272
0.0035108 0.0002136
0.0003079 0.0001068
-0.0012935 5.341E-05
-0.0004928 2.67E-05
-9.242E-05 1.335E-05
0.0001078 6.676E-06
7.672E-06 3.338E-06
-4.237E-05 1.669E-06
-1.735E-05 8.345E-07

F(c) Error
197.3487 3.5
22.367449 1.75
-65.953679 0.875
-21.387287 0.4375
0.4969027 0.21875
-10.440729 0.109375
-4.9714306 0.0546875
-2.2371812 0.0273438
-0.870122 0.0136719
-0.1866057 0.0068359
0.1551494 0.003418
-0.0157279 0.001709
0.0697108 0.0008545
0.0269915 0.0004272
0.0056318 0.0002136
-0.005048 0.0001068
0.0002919 5.341E-05
-0.0023781 2.67E-05
-0.0010431 1.335E-05
-0.0003756 6.676E-06
-4.187E-05 3.338E-06
0.000125 1.669E-06
4.157E-05 8.345E-07
ncluimos, que el número de iteraciones será "23".
Resolver el ej.1 usando el metodo de punto fijo o Newton.

P a
F(V)={(P+a/V²)*(V-b)}-RT = (P*V)+(a/V)-(P*b)-(a*b/V^2)-RT He 10 0.0341
F'(V)=P-(a/V^2)+(2*a*b)/V^3
F''(V)=(-2*a/V^3)-(6*a*b/V^4) 20 0.341
30 0.341
100 0.341

P= 10

N Vo F(Vo) F´(Vo) F´´(Vo) Verifico que el Vo tomado, cumpla con las condiciones de conv
1 4 8.67436989 9.9981213 0.00087621
2 3.13240001 0.00041271 9.99705054 (F(Vo)*F´´(Vo))/F´(Vo) = 0.0007602
3 3.13235873 1.46017E-12 9.99705047

P= 20

N Vo F(vo) F´(Vo) F´´(Vo) Verifico que el Vo tomado, cumpla con las condiciones de conv
1 2 6.31137958 19.9934954 0.00549436
2 1.68432836 0.00031386 19.9913627 (F(Vo)*F´´(Vo))/F´(Vo) = 0.00173441
3 1.68431266 1.01252E-12 19.9913626

P= 30

N Vo F(Vo) F´(Vo) F´´(Vo) Verifico que el Vo tomado, cumpla con las condiciones de conv
1 2 23.9413796 29.9934954 0.00549436
2 1.20178094 0.00255723 29.9857019 (F(Vo)*F´´(Vo))/F´(Vo) = 0.0043857
3 1.20169566 5.83533E-11 29.9857005

P= 100

N Vo F(Vo) F´(Vo) F´´(Vo) Verifico que el Vo tomado, cumpla con las condiciones de conv
1 4 347.34437 99.9981213 0.00087621
2 0.52649104 0.02106735 99.987735 (F(Vo)*F´´(Vo))/F´(Vo) = 0.00304352
3 0.52628034 -3.6386E-09 99.9877696
b R T V
0.237 0.08205 353 3.13235873
0.237 0.08205 353 1.68431266
0.237 0.08205 353 1.20169566
0.237 0.08205 353 0.52628034

cumpla con las condiciones de convergencia

-------> Cumple!

cumpla con las condiciones de convergencia

-------> Cumple!

cumpla con las condiciones de convergencia

--------> Cumple!

cumpla con las condiciones de convergencia

--------> Cumple!
La relación del espacio en función del tiempo, para dos partículas,
están dadas por las siguientes funciones:
Xp1(x) = 3*t^3 – 2*t^2 + t – 4
Xp2(x) = 4*t^2 + t + 6

a) Usando el método de bisección encuentre el valor de velocidad, para


cada partícula, en una posición de encuentro; considerando una
precisión del orden de una centésima.
b) Justifique el número de iteraciones requeridas, por el método, para
encontrar la respuesta del punto anterior.
Nota: la resolución del punto a, implica la construcción de una tabla
donde se pueda ver claramente la evolución de las variables
significativas, considerar la evolución del error

Método de Bisección

Xp1 (x) = 3*t^3 – 2*t^2 + t – 4


Xp2 (x) = 4*t^2 + t + 6

F(x) = Xp1 - Xp2 = 0 Raices de F(x) ---------> X1 =


3 * t^3 - 6* t^2 - 10 = 0 X2 =
X3 =

A)

N n a b c (t) h(a) h(b) CUMPLE


1 0 1.5 3.5 2.5 -13.375 45.125 CUMPLE
2 1 2.5 3.5 3 -0.625 45.125 CUMPLE
3 2 2.5 3 2.75 -0.625 17 CUMPLE
4 3 2.5 2.75 2.625 -0.625 7.015625 CUMPLE
5 4 2.5 2.625 2.5625 -0.625 2.9199219 CUMPLE
6 5 2.5 2.5625 2.53125 -0.625 1.0808105 CUMPLE
7 6 2.5 2.53125 2.515625 -0.625 0.2115173 CUMPLE
8 7 2.515625 2.53125 2.5234375 -0.210804 0.2115173 CUMPLE

B)
N= ln(b-a)-ln(error)/ln(2) = 7.3370034 El número de iteraciones será "8"
2.523462 La raíz a estudiar estará en el X1= 2.5234.
Imag El intervalo que tomaré será [1.5;3.5]
Imag

h(c) Error Vp1 = p1 dX/dt -------> Vp1 = 9t^2-4t+1 = 48.215881


-0.625 1 Vp2 = p2 dX/dt -------> Vp2 = 8t^2+1 = 21.1875
17 0.25
7.015625 0.0625
2.9199219 0.015625
1.0808105 0.0039063
0.2115173 0.0009766
-0.210804 0.0002441
-0.0006633 6.104E-05
La distancia, d, recorrida por un cuerpo en caída libre sobre cualquier planeta
sigue una ley cuadrática con el tiempo del tipo, d = R t2, donde R es la mitad de la
aceleración gravitatoria en ese planeta.
Se midió la distancia recorrida en una serie de instantes y se obtuvo la tabla
que abajo se muestra.

a) Determinar la aceleración de la gravedad en el lugar.


b) Cuanto tiempo tardara en llegar al piso un objeto que cae desde 500 mt
de altura, con velocidad inicial de 0.
c) Cual es la distancia recorrida luego de transcurrido 7,5 y 14 segundos,
desde que se inicio el movimiento.

t [s] d [m]
2 8
4 31
6 71
8 130
10 200

Modelo Aplicado

y= b +ax
y = ax^b
Log y = Log a + b*Log x

d= Rt^2
Ln d = Ln R + 2 Ln t

Referencia

y = log d
bo = log R
a0 = 2
x = log t

Log t Log d (logt)^2 logd*logt n Bo =


0.30103 0.90309 0.0906191 0.2718572 1
0.60206 1.4913617 0.3624762 0.8978892 2
0.7781513 1.8512583 0.6055194 1.440559 3 R = 10^Bo =
0.90309 2.1139434 0.8155715 1.9090811 4
1 2.30103 1 2.30103 5 Entonces
∑ 3.5843312 8.6606834 2.8741862 6.8204165
d = 1.96*t^2

A)

a = 2*R

a= 3.9231565

B)

d = 1.96*t^2 d= 500

500 = 1.96*t
t= (500/1.96)^1 15.965488

C)

t= 7.5
d(t=7.5)= 110.33878

t= 14
d(t=14)= 384.46933
0.2926056

1.9615782
Para los siguientes puntos: (1,3) , (3,5) y (5,6)
a) Encontrar la ecuación de la recta que mejor ajusta usando el método
de cuadrados mínimos.
b) Graficar
c) Calcular la distancia de cada punto a la recta
d) Calcular el desvío total como la suma de los cuadrados de las
distancias de cada punto a la recta.
e) Encuentre la ecuacion del polinomio cuadratici completo que mejor
ajusta a los puntos datos.
f) Calcular el desvío total como la suma de los cuadrados de las
distancias de cada punto a la cuadratica, compare los resltados con los
del ajuste lineal, que puede decir, según la teoria vista en clase.
g) Encotrar la ecuacion del polinomio de grado 3 completo que mejor
ajusta a los puntos datos.

X Y
1 3
3 5
5 6

A) Obtención de la Recta

Y = ax + b

X Y n X^2 X*Y
1 3 1 1 3
3 5 2 9 15
5 6 3 25 30
∑ 9 14 3 35 48

Fórmula Nª1

∑ X*Y = a * ∑ X^2 + b * ∑ X

Fórmula Nª2

∑Y = a*∑X+b*n

Despejamos alguna de las dos variables, la que más nos convenga, y la reemplazamos en la segunda fórmula.

a= 0.75 Y´ = 0.75X + 2.42 EC. DE LA RECTA


b= 2.4166667

C) Calcular la distancia de cada punto a la recta


X Y Y´ Y´-Y
1 3 3.1666667 0.1666667
3 5 4.6666667 0.3333333
5 6 6.1666667 0.1666667

D)

Desvio Tot. = 0.1666667

E)

X Y X^2 X^3 X^4 X^5 X^6 X*Y


1 3 1 1 1 1 1 3
3 5 9 27 81 243 729 15
5 6 25 125 625 3125 15625 30
∑ 9 14 35 153 707 3369 16355 48

N ∑X ∑(X^2) a0 ∑Y
∑X ∑(X^2) ∑(X^3) * a1 .= ∑(Y*X)
∑(X^2) ∑(X^3) ∑(X^4) a2 ∑(Y*X^2)

3 9 35 a0 14
9 35 153 * a1 .= 48
35 153 707 a2 198

a0 = 1.62
a1 = 1.5 --------> Y´´= 1.62+1.5x+(-0.13)*x^2
a2 = -0.13

F)

X Y Y´´ Y´´-Y (Y´´-Y)^2


1 3 2.99 0.01 1E-04
3 5 4.95 0.05 0.0025
5 6 5.87 0.13 0.0169

Desvio total = 0.0195

Podemos concluir que el polinomio de grado "2", ajusta de manera más eficiente a los puntos en cuestión. Siendo más pr

G)

No se puede formar un polinomio de grado "3", ya que el grado debe ser menor a la cantidad de puntos dados.
B) Graficar

Y´ = 0.75X + 2.42 Puntos brindados por el enunciado

X Y X Y
0 2.4166667 1 3
1 3.1666667 3 5
2 3.9166667 5 6
3 4.6666667
4 5.4166667
5 6.1666667
6 6.9166667
7 7.6666667
8 8.4166667
9 9.1666667
10 9.9166667

en la segunda fórmula.
YX^2 YX^3 N
3 3 1
45 135 2
150 750 3
198 888

untos en cuestión. Siendo más preciso que el ajuste lineal.

ntidad de puntos dados.


ados por el enunciado
a- Dados los puntos (x;y): (1;-3), (3;1), (4;10), (11/2;8) construir
el polinomio cubico interpolador por tramos natural.
b-
Determine la imagen esperada para un x=2.

c- Encuentre el/los valores de x que anulan el


polinomio interpolador, justifique la elección del procedimiento
utilizado.
d- Determinar el valor de la derivada en
x=7/2, usando el método de diferencias centradas O de orden
dos, con h=1.
e- Cual es el valor del error cometido en
la determinación del valor para la derivada del punto anterior.
X Y h0 h1 h2 d0
0 1 -3 2 1 1.5 2
1 3 1
2 4 10
3 5.5 8
Por polinomio cúbico natural > m0
m3

A 42 = 0 + 6 * m1 + 1
B -62 = 1 * m1 + 5 * m2

A m1 = (42 - m2) / 6
B -62 = 7 - m2/6 + 5 * m2
B -69 = 4.833333 m2
B -14.27586 = m2

A m1 = 9.37931

S0,0 = -3 S1,0 = 1 S2,0 =


S0,1 = -1.126437 S1,1 = 8.252874 S2,1 =
S0,2 = 0 S1,2 = 4.689655 S2,2 =
S0,3 = 0.781609 S1,3 = -3.942529 S2,3 =

S0(x) = -3 + -1.126437 (X-X0) + 0 (X-X0)^2 +


S1(x) = 1 + 8.252874 (X-X1) + 4.689655 (X-X1)^2 +
S2(x) = 10 + 5.804598 (X-X2) + -7.137931 (X-X2)^2 +

S0(2) = -3.344828 Esta OK


S1(2) = 1.37931

S0(x) = -3 + -1.126437 (X-X0) + 0 (X-X0)^2 +


S0(1) = -3
S0(3) = 0.996

i a b c So(a) So(b) Verifica So(c)


0 1 3 2 -3 0.996 si -3.344828
1 2 3 2.5 -3.344828 0.996 si -2.051724
2 2.5 3 2.75 -2.051724 0.996 si -0.782328
3 2.75 3 2.875 -0.782328 0.996 si 0.04014

Si X=2.875 se anula el polinomio

D
x h
3.5 1

S´x = (S0(4,5)-S0(2,5))/2*1

S0(4,5) = 26.56
S0(2,5) = -2.05

S´x = 14.305
d1 d2 U1 U2
9 -1.333333 42 -62

= 0 m1 = ?
= 0 m2 = ?

* m2
+ 0

10
5.804598
-7.137931
1.586207

0.781609 (X-X0)^3 Esta OK


-3.942529 (X-X1)^3 Esta OK
1.586207 (X-X2)^3 Esta OK

0.781609 (X-X0)^3

También podría gustarte