Certamen 3 Mat 270 1° 2005

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

UNIVERSIDAD TECNICA FEDERICO SANTA MARIA

DEPARTAMENTO DE MATEMATICA

ANALISIS NUMERICO (MAT-270)


CERTAMEN No. 3
24 de Junio de 2005

R x −t2 /2
1. (25 puntos) Para resolver la ecuación 0 e√2π dt = 0.45 usando el método de Newton-
R x −t2 /2
Raphson con f (x) = 0 e√2π dt − 0.45 se requiere una fórmula de cuadratura que
permita evaluar f (xn ), para n = 0, 1, 2, ... .

(a) Escriba, explı́citamente la fórmula de Simpson compuesta para calcular el término


f (x0 ) en la primera iteración de Newton, con x0 = 0.5
(b) Determine el número de subintervalos para aproximar el valor f (0.5) con una
precisión de 10−5 .
5 3 2
√ −3t e−t /2 para t ∈ [0, ∞[. Pruebe que
2. (25 puntos) Considere la función g(t) = t +t2π
al usar un cambio de variables adecuado
R∞ y cuadraturas Gaussianas en dos puntos se
obtiene el valor exacto de la integral 0 g(t)dt.
3. (25 puntos) Considere la siguiente muestra de la función g(x), definida como el pro-
ducto de convolución de una función f (x) con la constante 1 en el intervalo [−1, 1]:

x -0.92388 -0.382683 0.382683 0.92388


g(x) 0.141093 0.70772 1.58517 3.56601

Con esa pequeña muestra se puede calcular la aproximación cúbica por mı́nimos
cuadrados de g(x) utilizando polinomios ortogonales de Chebyshev (para ”poder ver”
el producto de convolución). Calcular el coeficiente de Fourier correspondiente al poli-
nomio de Chebyshev de grado uno, esto es el coeficiente en T1 .
La fórmula de integración de Gauss-Tchebyshev:
Z 1 n   
f (x) π X 2k − 1 2π
√ dx = f cos π + 2n f (2n) (ξ).
−1 1 − x2 n k=1 2n 2 (2n)!

4. (25 puntos) Se considera la torsión de un sistema compuesto por dos discos de iguales
caracterı́sticas y separados por tubos iguales. Se trata de ver la evolución de los ángulos
de giro en el tiempo, problema modelado por el sistema de ecuaciones diferenciales
ordinarias:
I 00
2ϕ1 − ϕ2 = − ϕ
R 1
I
−ϕ1 + 2ϕ2 = − ϕ002
R
(a) Escribir lo anterior como un sistema de ecuaciones diferenciales de primer orden,
suponiendo que la rigidez R y la inercia I son iguales.
(b) Si se parte de una situación de reposo en la que las torsiones iniciales son (en
radianes):
(0) (0)
ϕ1 = −0.5 ϕ2 = 0.0
Efectuar dos pasos de un esquema de Euler implı́cito tomando un incremento
temporal de 1 [uds].

RAF/JFN/EHH/.

También podría gustarte