Tarea Métodos Numéricos

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 5

UNIVERSIDAD NACIONAL DE LOJA

INGENIERÍA ELECTROMECÁNICA

MÉTODOS NUMÉRICOS

Trabajo Grupal

“Ajuste de modelos no lineales mediante cambios de variables”

Autores:

● Ríos Víctor
● Barragán Alexander
● Fernández Danny
● Uyaguary Santiago
● Gutiérrez Steven

Docente:

Ing. Byron Agustín Solórzano Castillo. Ms. Sc.

2020
 Cuando un cuerpo sólido se mueve a través de un fluido viscoso se produce sobre aquel
una fuerza fr que se opone al movimiento (fuerza resistiva) y que depende de la
velocidad del cuerpo. Con vista a esclarecer esta dependencia, se efectúa un
experimento en el cual un cuerpo se hace mover a diferentes velocidades a través de un
medio viscoso y se mide la fuerza de fricción. Los resultados obtenidos se muestran en
la siguiente tabla. Halle el modelo del tipo fr = μ 𝑣 𝑘 que mejor se ajusta a estos datos.

Modelo lineal de la forma y = a 𝑥 𝑏


y = fr
a=μ
x=v
b=k
● aplicando logaritmos tenemos lo siguiente.

ln(fr) = ln(μ 𝑣 𝑘 )
ln(fr) = ln(μ) + k(lnv)
ln(y) = ln(μ) + k(lnv)
Y = ln(fr); X = ln(v); C1 = k; C2 = ln(μ)

● Obtenemos la siguiente ecuación.

Y= C1X + C2
G1 = X
G2 = 1
G1*G1 + G1*G2 = G1Y
G1*G2 + G2*G2 = G2Y

 Tabla de valores originales.

v fr Gráfica de datos
700
1 50 600
500
Fricción

2 153 400
300
3 293 200
100
4 464
0
0 1 2 3 4 5 6
5 663
velocidad
 Tabla de valores, con cambio de variable. Y = ln(fr); X = ln(v)

X Y

0 3.9120

0.6931 5.0304

1.0986 5.6802

1.3863 6.1399

1.6094 6.4968

𝐺1 ∙ 𝐺1 𝐺1 ∙ 𝐺2 𝐶1 𝐺1 ∙ 𝑦
[ ]∙[ ]=[ ]
𝐺1 ∙ 𝐺2 𝐺2 ∙ 𝐺2 𝐶2 𝐺2 ∙ 𝑦
𝐺1 ∙ 𝐺1 = 𝑋 ∗ 𝑋 = 𝑋 2
𝐺1 ∙ 𝐺2 = 𝑋 ∙ 1 = 𝑋
𝐺2 ∙ 𝐺2 = 1 ∙ 1 = 1
𝐺1 ∙ 𝑌 = 𝑋 ∙ 𝑌
𝐺2 ∙ 𝑌 = 1 ∙ 𝑌 = 𝑌
 Tabla de valores de las sumatorias de cada una de los términos.

G2*Y=1*Y
G1G1;( 𝑋 2 ) G1G2;(X) G2G2=1*1=1 G1*Y=X*Y
=Y
0 0 1 0 3.9120

0.48038 0.6931 1 3.4866 5.0304

1.20692 1.0986 1 6.2403 5.6802

1.92182 1.3863 1 8.5117 6.1199

2.59016 1.6094 1 10.4559 6.4968

28.6945
6.1993 4.7874 5 27.2593

 Obtenemos el siguiente conjunto de matrices.

6,1993 4,7874 𝐶1 28,6945


[ ]∙[ ]=[ ]
4,7874 5 𝐶2 27,2593

 De la matriz anteriormente encontrada despejamos 𝐶1 y 𝐶2 .

6.1993𝐶1 + 4.7875𝐶2 = 28.6945 → 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 1


4.7875𝐶1 + 5𝐶2 = 27.2593 → 𝐸𝑐𝑢𝑎𝑐𝑖ó𝑛 2
 Aplicando el método de la eliminación Gaussiana se obtiene.

6.1993𝐶1 + 4.7875𝐶2 = 28.6945


[ ]
4.7875𝐶1 + 5𝐶2 = 27.2593

𝐹1
𝐹2 = 𝐹2 + ( ) (−4.7875)
6.1993

6.1993𝐶1 + 4.7875𝐶2 = 28.6945
[ ]
0𝐶1 + 1.3028𝐶2 = 5.0996

 Despejando 𝐶2 se obtiene.
5.0996
𝐶2 =
1.3028
𝐶2 = 3,9142

 Reemplazando 𝐶2 en la ecuación 1 se obtiene el valor de 𝐶1 .

6.1993𝐶1 + 4.7875(3,9143) = 28.6945


28.6945 − 18.7392
𝐶1 =
6.1993
𝐶1 = 1,6059

𝐶1 = 1,6059
𝐶2 = 3,91424
𝐶1 = 𝑘
𝑘 = 1,6059
 De C2 necesitamos despejar 𝜇

𝑒(𝐶2) = (ln 𝜇)𝑒


𝜇 = 𝑒 𝐶2
𝜇 = 𝑒 3,91424
𝜇 = 50,111
ln 𝑦 = 𝑘 ln 𝑦 + ln 𝜇
ln 𝑦 = ln 𝑥 𝑘 + ln 𝜇
𝑒(ln 𝑦) = ln(𝜇 ∙ 𝑥 𝑘 ) 𝑒
𝑦 = 𝜇 ∙ 𝑥𝑘
𝑦 = 50,111 ∙ 𝑥 1,6059
𝑦 = 50,111 ∙ 𝑣1,6059

 El modelo del tipo fr = μ 𝑣 𝑘 que mejor se ajusta a estos datos.


𝑓𝑟 = 50,111 ∙ 𝑣1,6059
 Representación gráfica con la ecuación exponencial en Excel.

Gráfica y ecuación de un Modelo no lineal


700
y = 50,111x1,6059
600

500
Fricción

400

300

200

100

0
0 1 2 3 4 5 6
velocidad

- La línea de puntos rojos es el resultado de la ecuación que mejor se ajusta a este


conjunto de datos dados, se puede apreciar que tiene una gran aproximación a la
gráfica original de color azul.

También podría gustarte