Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Grupo 1 Taller 2 Parcial 3

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

UNIVERSIDAD DE LAS FUERZAS ARMADAS-ESPE

DEPARTAMENTO DE CIENCIAS EXACTAS


ÁREA DE ANALISIS FUNCIONAL
MÉTODOS NUMÉRICOS

INTEGRANTES:  Hurtado Thomas


 Calle Fabian
 Naranjo Giovanny
 Ponce Evelyn
 Guaman Jonathan
 Zurita Carlos
NRC: 4142

FECHA: 2022-08-15

En un circuito de voltaje impreso 𝐸 que tiene la resistencia 𝑅, la inductancia 𝐿 y la


capacitancia 𝐶 en paralelo, la corriente 𝑖 satisface la ecuación diferencial
𝑑𝑖 𝑑 2 𝐸 1 𝑑𝐸 1
=𝐶 2 + + 𝐸
𝑑𝑡 𝑑𝑡 𝑅 𝑑𝑡 𝐿
Supongamos que C=0.3 faradios, R=1.4 ohms, L=1.7 henrios y que el voltaje está dado por
𝐸(𝑡) = 𝑒 −0.06𝜋𝑡 𝑠𝑒𝑛(2𝑡 − 𝜋)
Si 𝑖(0) = 0, calcule la corriente 𝑖 con los valores 𝑡 = 0.1𝑗, donde 𝑗 = 0,1, . . .100.
PLANTEAMIENTO:
Para resolver la siguiente EDO’s de valor inicial lo que procedemos a hacer es hallar 𝑓(𝑥, 𝑦)
𝑑2 𝐸 𝑑𝐸 𝑑𝑖
mediante el cálculo de , despeje 𝑑𝑡 y hacemos el reemplazo de los valores de C,R,L y
𝑑𝑡 2 𝑑𝑡
𝑖0 a continuación, calcularemos h mediante la diferencia entre 𝑡1 − 𝑡0 y aplicaremos la
fórmula de Euler (𝑦𝑘+1 = 𝑦𝑘 + ℎ𝑓(𝑥𝑘 , 𝑦𝑘 )) ,este proceso será iterativo en un intervalo [0,10]
RESOLUCIÓN:
𝒅𝒊
 RESOLVEMOS 𝒅𝒕

𝑑𝐸
= 0.18849𝑒 −0.06𝜋𝑡 𝑠𝑒𝑛 (2𝑡) − 2𝑒 −0.06𝜋𝑡 𝑐𝑜𝑠(2𝑡)
𝑑𝑡
𝑑2 𝐸
= 3.96446𝑒 −0.06𝜋 𝑡 𝑠𝑒𝑛(2𝑡) + 0.75398𝑒 −0.06𝜋𝑡 𝑐𝑜𝑠(2𝑡)
𝑑𝑡 2
𝑑𝑖 𝑑 2 𝐸 1 𝑑𝐸 1
=𝐶 2 + + 𝐸
𝑑𝑡 𝑑𝑡 𝑅 𝑑𝑡 𝐿
3
= (3.96446𝑒 −0.06𝜋 𝑡 𝑠𝑒𝑛(2𝑡) + 0.75398𝑒 −0.06𝜋𝑡 𝑐𝑜𝑠(2𝑡))
10
5
+ (0.18849𝑒 −0.06𝜋𝑡 𝑠𝑒𝑛 (2𝑡) − 2𝑒 −0.06𝜋𝑡 𝑐𝑜𝑠(2𝑡))
7
10
+ (𝑒 −0.06𝜋𝑡 𝑠𝑒𝑛(2𝑡 − 𝜋))
7
= 0.73573𝑒 −0.06𝜋𝑡 𝑠𝑒𝑛(2𝑡) − 1.20237𝑒 −0.06𝜋𝑡 𝑐𝑜𝑠(2𝑡) = 𝑓(𝑡, 𝑖)
 CALCULAMOS h
𝑆𝑖 𝑗 = 0 → 𝑡0 = 0
𝑆𝑖 𝑗 = 1 → 𝑡1 = 0.1
𝑡1 − 𝑡0 = 0.1 = ℎ
Calculamos 𝑖, 𝐶 = 0.3 𝑓𝑎𝑟𝑎𝑑𝑖𝑜𝑠, 𝑅 = 1.4 𝑜ℎ𝑚𝑠, 𝐿 = 1.7 ℎ𝑒𝑛𝑟𝑖𝑜𝑠, con ℎ = 0.1 para
𝑑𝑖
aproximar la solución de 𝑑𝑡 en el intervalo 0 ≤ 𝑡 ≤ 10

𝑑𝑖
= 0.73573𝑒 −0.06𝜋𝑡 𝑠𝑒𝑛(2𝑡) − 1.20237𝑒 −0.06𝜋𝑡 𝑐𝑜𝑠(2𝑡); 𝑖(0) = 0
𝑑𝑡
Resolución de la ecuación diferencial por medio del método de Euler, utilizando la
herramienta Octave/Mathlab
Código implementado:

Datos que se obtuvieron:


METODO DE EULER
j 𝑡𝑗 𝑖𝑗+1
0 0 0
1 0,1 -0,1202
2 0,2 -0,2215
3 0,3 -0,3006
4 0,4 -0,3551
5 0,5 -0,3839
6 0,6 -0,3866
7 0,7 -0,3643
8 0,8 -0,3187
9 0,9 -0,2524
10 1 -0,1689
11 1,1 -0,0721
12 1,2 0,0338
13 1,3 0,1441
14 1,4 0,2545
15 1,5 0,3604
16 1,6 0,4580
17 1,7 0,5436
18 1,8 0,6143
19 1,9 0,6679
20 2 0,7029
21 2,1 0,7186
22 2,2 0,7151
23 2,3 0,6933
24 2,4 0,6547
25 2,5 0,6013
26 2,6 0,5360
27 2,7 0,4617
28 2,8 0,3816
29 2,9 0,2992
30 3 0,2178
31 3,1 0,1406
32 3,2 0,0703
33 3,3 0,0097
34 3,4 -0,0393
35 3,5 -0,0752
36 3,6 -0,0971
37 3,7 -0,1046
38 3,8 -0,0979
39 3,9 -0,0779
40 4 -0,0458
41 4,1 -0,0033
42 4,2 0,0475
43 4,3 0,1042
44 4,4 0,1645
45 4,5 0,2259
46 4,6 0,2858
47 4,7 0,3419
48 4,8 0,3922
49 4,9 0,4349
50 5 0,4686
51 5,1 0,4924
52 5,2 0,5055
53 5,3 0,5080
54 5,4 0,5000
55 5,5 0,4824
56 5,6 0,4561
57 5,7 0,4226
58 5,8 0,3833
59 5,9 0,3401
60 6 0,2949
61 6,1 0,2494
62 6,2 0,2055
63 6,3 0,1648
64 6,4 0,1289
65 6,5 0,0990
66 6,6 0,0761
67 6,7 0,0607
68 6,8 0,0532
69 6,9 0,0537
70 7 0,0618
71 7,1 0,0769
72 7,2 0,0981
73 7,3 0,1244
74 7,4 0,1546
75 7,5 0,1873
76 7,6 0,2212
77 7,7 0,2548
78 7,8 0,2869
79 7,9 0,3162
80 8 0,3416
81 8,1 0,3624
82 8,2 0,3779
83 8,3 0,3876
84 8,4 0,3915
85 8,5 0,3894
86 8,6 0,3818
87 8,7 0,3692
88 8,8 0,3522
89 8,9 0,3317
90 9 0,3086
91 9,1 0,2839
92 9,2 0,2587
93 9,3 0,2339
94 9,4 0,2106
95 9,5 0,1896
96 9,6 0,1716
97 9,7 0,1572
98 9,8 0,1469
99 9,9 0,1410
100 10 0,1394

Es por ello que i(10) es igual a 0.1394 Amperios

CONCLUSIÓN:
En resumen, para resolver nuestra ecuación diferencial, lo que realizamos fue implementar
una derivación de nuestra función ya dada, por la cual obtendremos una segunda derivada
de la misma, para posteriormente reemplazar los valores de las resistencias, capacitores e
inductores, propuestos en nuestro ejercicio, y por medio de ello obtener valores de t en un
determinado intervalo, con el fin de obtener por medio del uso de la herramienta de Octave-
Matlab, los valores de resolución de nuestro ejercicio.
Cabe recalcar que se empleó el método de Euler ya que es una herramienta numérica que
empleamos para aproximar los valores para las soluciones de ecuaciones diferenciales.

También podría gustarte