PRÁCTICA 1

10/09/24 CALOR EXPLÍCITO

Dso
= 7 (xat) para (xt) ∈ Taib] ✗ [Tini, Thin]
2,4-2,/RAH 2)
Repasar ecuación
U /×, Tini) = No (x) para × Eta, b)
del calor!! I

Ula it) = ✗ (t) 24 2h

pares te [Tini, Trin]
2+-2×2 = 7
al bit) = PCH
Dominio de definición
✗ tu At de la solución UCxit)
Thin
}' Aquí conocemos la solución

α •
β Natación

• •
I ="n° de nodos iniciales"

Tini • • • •
N ="n° de pasos de tiempo"
a Δ! no [rig
" " " h
ti Xi XI
K

OBS: Se usarán siempre "mallas" (pasos uniformes.

OBS": Hay un test que necesitamos ejecutar para que el método

explícito funcione (EXPLICITheat. m). Se verá en la teoría.

NOTA IMPORTANTE Es recomendable

revisar el
Def.- 0/1) cuando h → 0" "O grande de h"
PDF de Landau
p.ej. 11h) = E. h en esta Carpeta

Si aplicamos una división usual, la transformación resultante

es directa. Es más, resulta ser lineal

111)-KE

pry:/ 11h1 = Kh + Kah'

Ahora la transformación también es casi lineal, pero no exacta
 Yonsideramos ahora 0/42) cuando h→ ♂ Para nuestro método
preferiremos un error
11h) - Kh' TRANSFORMACIÓN
CUADRÁTICA!!, con 0th) y así
pal 11%1=4%0 será más preciso

12/09/24 HAY QUE CORRER EL ARCHIVO prin EXPLICIT!!!

PREGUNTAS PRÁCTICA 1

Eje n°1
h
K K
h

a ½
XI-1 b Tini t EN_, Tfin
"X1 * " "
te tu
Sol/Ejecutando el código con N-500, I-1o se obtienen los

valores h= 1.1 E- 1 , K = 2 E- 3

2 /IR/lo k = HE-3 42 = 1.23 E- 2

Valor máximo (en la práctica nos referimos al máximo en los nodos)

OBS: Le dan valores de Ny de I para obtener los de Ky h

Eje n°2

Hay que ver en el explícito que se cumple la condición en

EXPLICIT heat. un 1h ↑ & Klatt)

Eje n°2

(N) K (I) h ERROR Ok/+ 01h4

2×10-3 1. 1×10-1
6. 18×10-5

10-1
102 NO
8.24×10-3
Lo que espero
103 10"
8. 24×10-4 ver
 (otras preguntas)

Eje n°1

Hay que dos valores constantes y ver si la solución deja

de moverse. Esto se comprueba yendo al OUTPUT y chequeando

que a partir de cierto modo se alcance dicho 'estado estacionario".

Eje ni 2

a soporte de No b

Solución con métodos numéricos

Solución del problema
continuo

(completar...?)

13/09/24

a (✗ it) = Xat", D= 1 , Taib/= [1,2] , [Tini,Tfin] = [0,1]

lo (x) = U (×,Tini) = U (Xo) = ✗ 3

LA) = 4 (art) = 4 (1H) = 1+7"

β A) = 4lb, +1 = 4/4+1 = 8+72

71×1 t) = 2¥-£ [DAH 27] = 217_22¥ = 2+-2

4E Aprox. de la sol. exacta en Cxi.tn)

OBS: El método es explícito siel "uti aparece solo en

un sitio, lo que permite despejar la incógnita,
Eje n°2
01k 1+0142)
⇔ (N) Error

5 100 2.83 E- 04
50 5000 3.59 E- 06 Se cumple
{= 2.09 E- 02 4=2.0004 E-04
4=2.5 E- 01 K 1. OIE-02
???

Ejc n°3

Sí. Véase con D (la 1ª)

Jon la última 2 no funciona. → Fallo de inestab.

Ejc n°4

Sí, es con continuidad ???

Eje n°5
No se hereda la No regularidad ???

17/09/24
VECTOR}
EXPLÍCITO VECTOR VECTOR 2

II."Ui-III.+ II.½) Ui + DE" ↳ { KIT (i. a)

l

4"" 4" ¥

4"-" =
I-I
44¥ I-I
, EI, LE III.+ EH) 4"+ DE,_Y {+ KI: "-
"

Hay que pensar el sistema rechorialmente!

 IMPLÍCITO Revisar
2+1
uti
"'_ (2-""+ I" - KF"" LEE {I-I
LII"-4:"-¥ Ri-½ Lin i-½ it ½ 14:"+ III. Wi. ,
- 4+1
EDIE 4:" = 4a" + KIE"
(i-z)
- ¥5"4.
2- ½
" Ir + ¥./Ray, + DE:) 4:"

19/09/24

crear matrices con alguna entrada específica

Yo mando Sparse: sirve para de ceros

p. ej/ A = sparse /Mil 15311, F- 4,61, 4,4)

full (A) filas colum. dimensión matriz

vector que
O -40
contiene nuestros
6 O O O
valores

O O O

dr = 2. * ones / ???) → 2 2 2 2

da = Sparse (1: 4, 1: 4, du, 4,4)

full (LA/ = Matriz diagonal de doses

udv = -1.* ones/1,3)

Udr = - 1,-11-1

Uda = sparse 1113, 2: 4, udr, 4,4)

0 -1 O °

0 -1 o

◦ o o →

0 0 0 0
 A- sparse (51: 4,113, 2:47, [1: 4, 214, 1:3], [du, Udr, Udr/, 4,4/;

full/A) . ²:­
)

☒ ESQUEMA FTFS.
KILI-1
¡"-Li" + c: Hii-Hi" = fin
\ h
14h EN-1

Reformulando! Véase guía de programas

01/10/24

Veamos cómo plantean problemas para comprobar programas:

TRANSPORTE
pujo
an
at + C 2£ =L condición alfa

Fijamos c: Cant) = 4 + cosa sin + →

Fijamos next): Mx,t 1=+7+76

Fijamos No (x): U/×,Tini) = ✗ " + Tinio

LCH = Wait) = 97+76 (que tenga Sentido con alxit)!!!)

Se calcula la 7 en función de los datos anteriores usando la EDP

7kt) = Ett c 227=675+14 + cosxtsint) 7×6

Ya podemos meter el problema en nuestro programa

para comprobarlo!!!
 ☒ PRÁCTICA 3.

PARTE Δ (walk 10)

Si ✓ = [3 ⇒ ✗= linspace (a, b, Nii)

Si no, en ✓ marcamos los puntos en los que queremos acumula

Pj/ Ktaib] se acumulen en los extremos

5/12/2024

El error que obtenemos al ejecutar nuestro algoritmo es la

suma de EMI-Chi

ERROR DISCREI + ERROR REDONDEO

(método numérico) (precisión máquina del ordenador)

Gemonmente, hay un lo (sub." pequeño) tal que pare valoro

de h inferiores, el peso del error de Redondeo es dominante

y los resultados se descontrolen.
ELN)
EIN) = CN-P En:(¼)" (E) P
b- a
h= N-s P → p-GEMIEIN))
N El N) log (MIN)
N- A+ b-a
h

Ni-10 ECM)

N = 32 ELNa)
Ejercicio

Storm- Liouville
(a.b) = (-1/1)

- (pa'l' + que 1
U 1×1 = 3A = /⅓

N 'A) = ⅓ %
Util =-1 , 4111=1
-113×43 x- ⅔ + ✗ =p
Pix) = 9- (x) = 7%

N ELN) p

1. 4393 . só'

32 1. 3294. 10" 0.0683

100 1. 0443.15' 0.2118
316 7. 6772. 152 0.2674

1000 5. 4640 - 102 0.2444

3162 3. 8134. 152 0.3091

2. 6574.152 0.3178
10000

El orden de convergencia sólo se observe si la solución y

suficientemente regular.

Pruébese a hacer la malla con ✓mesh-o.

☒ EXACTITUD EN LAS FÓRMULAS DE CUADRATURA

(Am""/ja = hi fuí E) 9- (4E)/Wilde se calcula exactamente

grado 1 grados al aproximar...

Componente (jak) de
con Poncelet o con trapecios si 7=0
la matriz elemental
de masa
con Gauss de 2 puntos si qe Pa
componente (jak) de
la matriz elemental
de rigidez

' (wj)'e) p/4:( )) (War)'(E)DE se calcula exactamente

(Art") ju-Ei
0 al aproximar...
grado o grado o

con Poncelet o trapecios si fe R

con Gauss de 2 puntos si pelPs

componente j-ésima
del 2º miembro
elemental

1
Wj (E) L (Y:( )) de te calcule exactamente
bj'"'-hi
al aproximar...
O
grados

con Poncelet o con trapecios si LE Po

con Gauss de 2 puntos si felpa

Casos de exactitud del MEF Lagrange Pe:

1.) HE Pr y los integrales se calcula exactamente.

2.) PETE , q ⇒ y los integrales se calculen exactamente

↳ -a"-1 [Formulación diferencia]

OJO!!! ⇒ En el punto 2 hablamos solo de exactitud en

los nodos. Obviamente si 4¢ Ps (pres., si es polinomio de

grado mayor o cualquier elm función sucia no polinómica a

grado ≤ 1) entonces a ≠ Uh ∈ In"