| Indice del artículo |
|---|
| Comparación entre métodos de integración numérica |
| Comparación entre el método de Euler y el trapezoidal |
| Comparativa práctica Euler-Trapecio con ANALOGIA.EXE |
| Todas las páginas |
Durante el desarrollo de este proyecto, la primera versión (que dio muy malos resultados) realizaba el análisis en frecuencia utilizando técnicas de análisis AC. La segunda versión hacía uso del algoritmo implícito de Euler, y la definitiva, del algoritmo del trapecio. Vamos a ver pues numéricamente y gráficamente (de forma práctica) las diferencias entre estos dos algoritmos y la superioridad del trapezoidal.
Sea el circuito:
que usaremos para las mediciones de exactitud y estabilidad, y donde Vr es la variable de interés.
Utilizando las aproximaciones a un condensador, ya estudiadas para cada algoritmo, obtenemos:
Euler Explícito: 
Euler Implícito: 
Trapezoidal: 
Se cumple en cualquier instante que: 
sustituyendo en las anteriores ecuaciones nos queda:
Euler Explícito: 
Euler Implícito: 
Trapezoidal: 
Sustituyendo RC por t, podemos escribir estas ecuaciones para producir relaciones iterativas:
Euler Explícito: 
Euler Implícito: 
Trapezoidal: 
Si se aplica una excitación en escalón de tensión E0, el valor inicial de Vr es Vr(0)=E0. Aplicando las ecuaciones anteriores n veces (con los incrementos siempre iguales), nos quedará:
Euler Explícito: 
Euler Implícito: 
Trapezoidal: 
La solución exacta para este circuito, sabemos que es:
y para t=n
Podemos ahora calcular las expresiones aproximadas y comparar con el valor exacto para ver cual produce mejores resultados. Si tomamos dos incrementos de ejemplo =0.1
y =, resulta:
| Salto | Exp.Euler
| Imp.Euler
| Trapecio
| Exacto
|
|---|---|---|---|---|
=0.1 | 0.90000 | 0.90909 | 0.90476 | 0.90483 |
| = | 0.0 | 0.5 | 0.333 | 0.368 |
Como vemos, para ambos incrementos el método del trapecio produce una mejor aproximación. La estabilidad de los métodos puede estimarse de forma aproximada buscando para qué valores de
la solución empieza a oscilar:
| Salto | Exp.Euler
| Imp.Euler
| Trapecio
|
|---|---|---|---|
| Estable para: | 0< | 0< | 0<<2 |
Para grande:  |  Divergente y oscilante |  Estable |  Estable pero oscilante |
Vemos que el método implícito de Euler y el del trapecio son estables, aunque el último presenta oscilaciones para >2 que se pueden eliminar disminuyendo el salto interno.
