Escuela Superior de Ingeniería y Arquitectura

Técnicas de Análisis Numérico

Problemario

Tabla de contenido
Tema. Regresión Lineal......................................................................................................................................1
Tema. Series de Taylor.......................................................................................................................................2
Tema. Polinomios de Interpolación de Newton................................................................................................3
Tema. Polinomios de Interpolación de Lagrange...............................................................................................3
Tema. Polinomios de Hermite............................................................................................................................3
Tema. Polinomios de Chevyshev........................................................................................................................4
Tema. Regresión Sinusoidal...............................................................................................................................4
Tema. Derivación Numérica...............................................................................................................................5
Tema. Integración Numerica - Regla del Trapecio y Reglas de Simpson...........................................................6
Tema. Integración Numerica - Cuadratura de Gauss e Integracion de Romberg..............................................6
Tema. Ecuaciones Diferenciales.........................................................................................................................7
Tema. Series de Fourier.....................................................................................................................................9
Tema. Transformada de Fourier......................................................................................................................10

Tema. Regresión Lineal

1. Con los siguientes datos use la regresión lineal para determinar la línea recta que mejor se ajuste:
y=a1 x +a o

x 1 2 4 7 10 12 15
y 3 5 7 11 14 17 21
2. El esfuerzo de cedencia, σ y, de muchos metales varía con el tamaño del grano. Frecuente, la relación entre el
tamaño del grano, d, y el esfuerzo de cedencia se modela con la ecuación de Hall-Petch:
y=a1 x +a o
1
x=
√d
(−12 )
σ y =σ 0+ k d
Los siguientes datos son los resultados de medir el tamaño promedio del grano y los esfuerzos de cedencia:
d (mm) 0.006 0.011 0.017 0.025 0.039 0.060 0.081 0.105
σy (MPa) 334 276 249 235 216 197 194 192

Determine las constantes σ0 y k de manera que la ecuación de Hall-Petch ajuste mejor a los datos. Grafique
los datos y la ecuación de Hall-Petch. Use la ecuación para estimar el esfuerzo de cedencia de una probeta
con un tamaño de grano de 0.05mm.

3. La resistencia R de un alambre de tungsteno está en función de la temperatura y se modela con la ecuación

R = R0 [1+α(T-T0)], donde R0 es la resistencia a la temperatura T 0, y α es el coeficiente de temperatura de la
resistencia. Determine R0 y α tal que usando una regresión lineal se ajuste mejor a los siguientes datos
usando T0 = 20 °C:
T (°C) 20 100 180 260 340 420 500
 R (Ω) 500 676 870 1060 1205 1410 1565
y=a1 x +a o

4. En un experimento para medir el coeficiente de expansión térmica, se coloca en un horno una barra de
acero inoxidable de 2.5 m de longitud. La temperatura de la barra se eleva desde 20 °C hasta 820 °C con
incrementos de 100 °C, y se mide la longitud de la barra cada vez. El cambio en la longitud ∆L se relaciona
con el cambio de temperatura ∆T mediante la función ∆L = αL 0∆T, donde α es el coeficiente de expansión
térmica, y L0 es la longitud inicial. Use la regresión lineal de mínimos cuadrados para determinar el
coeficiente de expansión térmica si se midieron los siguientes datos:
m
y=b x
ln ( y)=mln( x)+ ln(b)

∆T (°C) 100 200 300 400 500 600 700 800

∆L (mm) 4.2 8.9 17.3 14.8 23.5 28 30.8 34.2

5. La tasa de crecimiento de una bacteria puede modelarse con la siguiente ecuación: ln N t −ln N 0 =μ (t−t 0),
donde μ es la constante de crecimiento, N t y N0 son el número de bacterias al tiempo t y t 0 respectivamente.
Determine μ y N0 de manera que la ecuación obtenida ajuste a los siguientes datos usando t 0=0.
mx
y=b e
ln ( y)=m x +ln(b)

y=a1 x +a o
t(h) 0 2 4 6 8
N(células/ml) 35 1990 70800 2810000 141250000

6. En una prueba de tensión axial, un espécimen con forma de hueso es tensado en una maquina. Durante la
prueba se mide la fuerza aplicada al espécimen, F, y la longitud de la sección, L. El esfuerzo real, σ t y la
deformación real εt se definen como:
F L L
σt= ε t =ln
A 0 L0 y L0
Donde A0 y L0 son el área de la sección transversal y la longitud inicial. La curva esfuerzo-deformación en la
región posterior al esfuerzo de cedencia se modela como:
m
σ t = Kεt
Los siguientes valores de F y L se midieron en un experimento. Transforme la ecuación no lineal a su forma
lineal para determinar los valores de los coeficientes K y m que se ajustan mejor a los datos. El área de la
sección transversal y la longitud inicial son A 0=1.25x10-4m2 y L0=0.0125m.

F(kN) 24.6 29.3 31.5 33.3 34.8 35.7 36.6 37.5 38.8 39.6 40.4
L(mm) 12.58 12.82 12.91 12.95 13.05 13.21 13.35 13.49 14.0 14.21 14.48
8
m
y=b x
ln ( y)=mln( x)+ ln(b)

Tema. Series de Taylor

7. Dado f (x)= x3 −8
a. Encuentre el polinomio de segundo orden de Taylor P 2(x) alrededor de x0=0.
b. Encuentre el error verdadero cuando se usa P 2(0.5) para aproximar f(0.5).
c. Repita el inciso a. con x0=1.
d. Repita el inciso b. para el polinomio encontrado en el inciso c.

8. Dado f ( x )= √ x +1
a. Encuentre el polinomio de tercer orden de Taylor P 3(x) alrededor de x0=0.
b. Use P3(x) para aproximar √ 0.5 , √ 0.75 , √ 1.25 y √ 1.5
 c. Determine el error verdadero de las aproximaciones del inciso b.

9. Encuentre el polinomio de Taylor de tercer orden P 3(x) para la función f ( x )=( x−1 ) ln(x) alrededor de x0=1
y P3(x) para aproximar f(0.5).

10. En el caso de que x 0=0, el polinomio de Taylor se conoce como el polinomio de Maclaurin, y la serie de Taylor
se conoce como la serie de Maclaurin. Encuentre el polinomio de tercer orden de Maclaurin para la siguiente
función:
f ( x )=e x/ 2 sin ⁡¿

11. Use un polinomio de Taylor respecto a π/4 para aproximar sen 42° con una precisión de 10 -6.

12. Evalúe la siguiente función en x=0:

e x −1
f ( x )=
x
Pista: Expanda la serie de Taylor e x hasta cuarto orden del dividendo alrededor de x=0.

Tema. Polinomios de Interpolación de Newton

13. Dados los siguientes datos:

x 0.5 1 2 3 4 6
y 2 3.5 4.5 3.5 4 6

Escriba el polinomio de Newton y úselo para calcular el valor de y en x=3.5

14. Los siguientes datos:

x -2 1 4 -1 3 -4
y -1 2 59 4 24 -53
Caen en un polinomio. Determine el grado del polinomio la tabla de diferencias divididas.

15. Use los datos:

x 0 0.5 1 1.5 2
y -0.7854 0.6529 1.7390 2.2071 1.9425
Encuentre y en x=π/4 y x=π/2. Use el método de interpolación de Newton.

16. Los siguientes datos:

x -3 2 -1 3 1
y 0 5 -4 12 0
Caen en un polinomio. Determine el grado del polinomio usando la tabla de diferencias divididas del método de
Newton.

Tema. Polinomios de Interpolación de Lagrange

17. Encuentre la raíz del polinomio de Lagrange usando a) tres, y b) cuatro puntos más cercanos.

x 0 0.5 1 1.5 2 2.5 3

y 1.8421 2.4694 2.4921 1.9047 0.8509 −0.411 −1.5727
2

18. El factor de concentración de esfuerzos k es la relación entre el esfuerzo máximo σ max y el esfuerzo
promedio σ ave, k =σ max /σ ave. Para una placa de ancho D con dos muescas laterales de longitud d y cargada
con una fuerza axial F, el esfuerzo máximo se encuentra en el extremo de la muesca. El esfuerzo promedio
esta dado por σ ave =F/[t ( D−2 d ) ], donde t es el espesor de la placa. El factor de concentración de
 esfuerzos fue medido en cinco pruebas con placas con varias relaciones de d/D y los resultados se muestran
en la siguiente tabla:

d/D 0.05 0.15 0.25 0.35 0.45

k 3.84 4.63 4.4 3.63 2.22

Use el polinomio de Lagrange de tercer grado para modelar la relación entre k y d/D. Grafique los datos y el
polinomio encontrado. Use el modelo para predecir el factor de concentración de esfuerzos para d/D=0.21.

Tema. Polinomios de Hermite

19. Use la interpolacion de Hermite para construir un polinomio de aproximación para los siguientes
datos:
a. b.
x f ( x) f ’ ( x) x f (x) f ’ ( x)
8.3 17.56492 3.116256 0.8 0.22363362 2.1691753
8.6 18.50515 3.151662 1.0 0.65809197 2.0466965
c. d. D
x f ( x) f ’ (x) x f ( x) f ’ ( x)
-0.5 -0.0247500 0.7510000 0.1 -0.62049958 3.58502082
-0.25 0.3349375 2.1890000 0.2 -0.28398668 3.14033271
0 1.1010000 4.0020000 0.3 0.00660095 2.66668043
0.4 0.24842440 2.16529366

20. Los datos del ejercicio anterior se generaron por medio de las siguientes funciones. Use los polinomios
construidos en el ejercicio anterior para el valor dado de x y con ellos aproxime f (x) y calcule el error real.
a. f ( x )=x ln x ; aproximar f ( 8.4)
b. f ( x )=sin ( e x −2 ) ; aproximar para f (0.9)
3 2 −1
c. f ( x )=x + 4.001 x + 4.002+1.101 ; aproximar para f ( )
3
d. f ( x )=x cos x−2 x 2+3 x−1 ; aproximar f (0.25)

21. Un automóvil realiza un recorrido por una carretera y se cronometra su recorrido en varios puntos. Los datos
recabados de las observaciones se incluyen en la tabla adjunta, donde el tiempo se indica en segundos, la
distancia en pies y la velocidad en pies por segundo.

Tiempo 0 3 5 8 13
Distancia 0 225 383 623 993
Velocidad 75 77 80 74 72

a. Use el polinomio de Hermite para predecir la posición del automóvil y su velocidad cuando t= 10 s.
b. Use la derivada del polinomio de Hermite para determinar si el automóvil rebasa el limite de
velocidad de 55 mi/h en la carretera. De ser así, ¿Cuál s la primera vez que la excede?
c. ¿Cuál es la velocidad máxima predecible del automóvil?

Tema. Polinomios de Chevyshev

~
22. Use los ceros de T 3 y construya un polinomio de interpolacion de segundo grado de Chevyshev paara las
siguientes funciones en el intervalo [-1,1].
a. f ( x )=e x b. f ( x )=sin x
c. f ( x )=ln (x +2) d. f ( x )=x 4

~
23. Use los ceros de T 4 y construya un polinomio interpolante de tercer grado para las funciones del ejercicio
anterior.
 Tema. Regresión Sinusoidal

24. La siguiente tabla muestra la temperatura diaria mensual de 1961 a 1990 en Halifax, Nueva Escocia en los
Estados Unidos de Norteamérica.

Mes Ene Feb Mar Abr May Jun Jul Ago Sep Oct Nov Dic
Temperatura - -6 -1.7 3.6 9.4 14.7 18.3 18.1 13.8 8.5 3.4 -3
5.8

Grafique y efectué el ajuste sinusoidal de los datos. Calcule el periodo de la función e interpreté el valor.
Determine la temperatura promedio para el 15 de Junio. Con base a los parámetros de la función, diga ¿cuál
es el valor de la temperatura promedio anual?, así como el rango de temperaturas esperadas. ¿Cuál sería el
efecto en la gráfica de función sinusoidal si el promedio de la temperatura anual se incrementara 2 oC?. La
temporada de siembra de algunos vegetales es la parte del año en la cual la temperatura permanece por
arriba de cierto valor. Para los frijoles, la temperatura debe estar por arriba de los 6 oC. ¿Cuál es la temporada
de siembra para ese tipo de frijoles en Halifax?.

25. Las olas más grandes del mundo se presentan en Minas Basin, Nueva Escocia. En la siguiente tabla se
muestran las olas medidas el 9 de abril del 2003.

Hr 00 02 03 06 07 08 11 14 15 17 20 21 22 23
Altur 3.3 4.7 6. 10.5 10.6 9.8 4. 3.7 5.1 8. 10.2 8.8 7.1 5.4
a 3 6 4

Efectue un ajuste con funciones sinusoidales. Si un barco camaronero de grandes dimensiones necesita un
tirante de 4 metros para flotar, ¿cuánto tiempo estará varado en Minas Basin? ¿Cuál es el tirante de agua a las
14:36?

Tema. Derivación Numérica

26. Dados los siguientes datos para tres funciones encuentre la primera y segunda derivada en los puntos 1,
1.2472 y -0.1 respectivamente con 0(h2).

x f 1 ( x) x f 2 ( x) x f 3 ( x)
0.8000 −1.0880 0.8472 0.7494 −0.2000 1.2214
0.9000 −1.0710 0.9472 0.8118 −0.1000 1.1052
1.0000 −1.0000 1.0472 0.8660 0 1.0000
1.1000 −0.8690 1.1472 0.9116 0.1000 0.9048
1.2000 −0.6720 1.2472 0.9481 0.2000 0.8187

27. Los siguientes datos muestran la población de Nepal entre 1980 y 2005:

Año 1980 1985 1990 1995 2000 2005

Población (millones) 15 17 19.3 22 24.5 27.1

Calcule la tasa de crecimiento de la población en millones por año para el año 2005. Use la pendiente en
2005 y aplique la formula de diferencias centrales para extrapolar y predecir la población en el año 2010.

28. El voltaje E=E(t) en un circuito eléctrico obedece a la ecuación E ( t )=L ( dldt )+ RI (t), donde R es la
resistencia y L es la inductancia. Use L=0.05, R=2 y los valores de I (t) de la siguiente tabla para:
a. Encontrar I ’ (1.2) por diferenciación numérica y úsela para calcular E(1.2)
b. Compare sus resultados con I (t)=10e−t /10 sen( 2t) .
 t 1.0 1.1 1.2 1.3 1.4
I (t) 8.2277 7.2428 5.9908 4.5260 2.9122

29. Una viga de 10 metros de longitud con apoyos simples está sujeta a una carga triangular. La deflexión de la
viga está dada por la ecuación diferencial:
d 2 y −M (x)
2
=
dx EI
Donde y es la deflexión, x es la coordenada medida a lo largo de la viga, M(x) es el momento flexionante, y
EI =1.2 x 107 N −m2. Los siguientes datos se obtuvieron:

x(m) 0 1 2 3 4 5 6 7 8 9 10
y(mm) 0 -16 -30.6 -42.5 -50.7 -54.3 -52.6 -45.7 -33.9 -18.1 0

Usando los datos determine el momento flexionante M(x) en cada posición de x.

30. Dados los siguientes datos encuentre la primera derivada f’(x) en el punto x= 0.399.
x 0.389 0.399 0.400 0.401 0.402
f(x) 0.0630 0.06375 0.0649 0.0654 0.0658

31. Los siguientes datos corresponden a la velocidad de un vehículo durante una prueba de colisión:

t(ms) 0 10 20 30 40 50 60 70 80
v(mph) 30 29 27 24 18 12 5 1 0

Si el peso del vehículo es de 2000 lb, determine la fuerza instantánea F actuando en el vehículo en la colisión.
dv
La fuerza se puede calcular como F=m , y la masa del automóvil m es 2000/32.2 slug.
dt
Note que 1 ms = 10-3 s y 1 milla= 5280 pies.
Tema. Integración Numerica - Regla del Trapecio y Reglas de Simpson

32. Use la regla del trapecio de aplicación múltiple para evaluar:

π
4

∫ ln (1+tanx ) dx
0
Explique sus resultados.

33. La siguiente tabla muestra la potencia P aplicada a la transmisión de un automóvil en función de la velocidad
v. Si la masa del automóvil es m=2000 kg, determine el tiempo ∆t que tarda el automóvil para acelerar de 1
m/s a 6 m/s. Use la regla del trapecio para resolver la integral. Considere que:

6 m /s
∆ t=m ∫ ( Pv ) dv
1 m /s

La cual se deriva de la segunda ley de Newton F=m(dv/dt) y de la definición de potencia P=Fv.

v(m/s) 0 1. 1.8 2.4 3.5 4.4 5.1 6.0

0
P(kW) 0 4. 12.2 19.0 31. 40.1 43.8 43.2
7 8

34. La velocidad de un automóvil durante los primeros 6 segundos de una carrera está dada por:
t(s) 0 1 2 3 4 5 6
V(mi/h) 0 18 41 63 83 99 112
 Determine la distancia recorrida durante los primeros seis segundos.
a. Use la regla de Simpson que aplique.
b. Use la regla del trapecio de aplicación múltiple.

35. Evalúe las siguientes integrales usando la regla de Simpson 1/3 con 2, 4 y 6 segmentos. Compare con la
solución analítica y explique sus resultados.
3π
1 2

∫ cos ⁡(2 cos−1¿ x )dx y ∫ cos3 xdx ¿

−1 π
2

1
36. Evalúe la integral ∫ 225 x (1−x ) √ (1−x ) dx usando la integración de Romberg O(h8).
2 2 3

Tema. Integración Numerica - Cuadratura de Gauss e Integracion de Romberg

1
37. Evalúe la integral ∫ 225 x (1−x ) √ (1−x ) dx usando la integración de Romberg O(h8).
2 2 3

38. Evalúe la integral con la cuadratura de Gauss-Legendre usando dos y cuatro nodos
π
ln x
∫ x 2−2 x+2
dx
1
39. Una viga de aluminio con forma de cono truncado de longitud L = 2 m, esta empotrada en el extremo base
del cono y se aplica una fuerza axial en el otro extremo P= 20 kN. El área de la sección transversal de la viga
varia con la longitud de la viga y se puede calcular como función de x:
R−r 2
A=π ( R− x)
L

Donde R= 50 mm y r=20 mm. La deformación normal, ε, a lo largo de la viga está dada por:
du P
ε= =
dx EA
Donde E= 70GPa es el modulo de elasticidad de la viga. El desplazamiento δ en el punto x= xP a lo largo de la
viga se puede calcular como:
xP
P
δ =∫ dx
0
EA
a) Determine el desplazamiento δ en el extremo de la viga (xP=L). Calcule el desplazamiento usando la
cuadratura de Gauss-Legendre y con alguna de las funciones de MATLAB.
b) Escriba un programa de Python que calcule y grafique el desplazamiento δ como una función de x.

40. Evaluate

∫ sen
√x
x
dx
0
with Romberg integration O(h8). Hint: use transformation of variable to eliminate the indeterminacy at x = 0.

41. Un automóvil recorre una pista de carreras en 84 segundos. Su velocidad en cada intervalo de 6 segundos se
determina mediante una pistola de radar y está dada, en pies/s, desde el principio del recorrido, por los
datos de la tabla siguiente:
Tiempo 0 6 12 18 24 30 36 42 48 54 60 66 72 78 84
Velocidad 124 134 148 156 147 133 121 109 99 85 78 89 104 116 123
¿Qué longitud tiene la pista?
42. Por medio de la integración de Romberg calcule las siguientes integrales con O(h 6):
2π
a) ∫ x senx dx
0
3
b) ∫ x √1+ x 2
0
1
c) ∫ senπx dx
0

Tema. Ecuaciones Diferenciales

43. Convierta la siguiente ecuación diferencial a una ecuación diferencial de primer orden de la forma
y ' =F ( x , y ) y resuélvala por el método del punto medio en el intervalo de 0 a 4 con h=0.5
'
ln y + y=senx y ( 0 )=0

44. Considere la siguiente EDO de primer orden

dy
=x+ y ; de x=0 a x=2.4 con y ( 0 ) =2
dx
 Resuelva con el método de Heun usando h=0.8
 Resuelva con el método clásico de cuarto orden de Runge-Kutta usando h=0.8
La solución analítica de la EDO es y=3 e x −x−1.

45. Considere la siguiente EDO de primer orden

dy 2
=2 y−x ; de x =1a x=2.8 con y ( 1 )=2
dx
 Resuelva con el método del punto medio usando h=0.6
 Resuelva con el método clásico de cuarto orden de Runge-Kutta usando h=0.6
1 1 2 3 2 x−2 1
La solución analítica de la EDO es y= x + x+ e +
2 2 4 4

46. Considere la siguiente EDO de primer orden

dy 2 3 y
=t − ; de t=1 a t=2.2con y ( 1 )=1
dt t
 Resuelva con el método explicito de Heun usando h=0.4
 Resuelva con el método del punto medio usando h=0.4
 Resuelva con el método clásico de cuarto orden de Runge-Kutta usando h=0.4
1 5 3
La solución analítica de la EDO es y= ( +t )
6 t3

47. Considere la siguiente EDO de primer orden

dy 2 2 y
=t − ;de t=1a t=4 con y ( 1 )=2
dt t
 Resuelva con el método del punto medio usando h=1
 Resuelva con el método clásico de cuarto orden de Runge-Kutta usando h=1
9 t3
La solución analítica de la EDO es y= +
5t 5
2

48. Considere el siguiente sistema de dos EDO:

dx dy de t=0 a t =2 con x ( 0 )=1 , y ( 0 )=2
=x+ y = y−x
dt dt
 Resuelva con el método de punto medio h=0.5
 Resuelva con el método de Heun usando h=0.5
La solución analítica del sistema es x=e t ¿
49. Considere el siguiente sistema de dos EDO:
dx dy de t=0 a t =1.2con x ( 0 )=1 , y ( 0 )=0.5
=xt− y = yt+ x
dt dt

 Resuelva con el método de Ralston usando h=0.4

 Resuelva con el método del punto medio usando h=0.4

[ ] [ ]
1 1
2
t 1 t 1
 La solución analítica del sistema es x=e cos ( t )− sin ( t ) , y =−e 2 −sin ( t ) − cos ( t )
2 2

50. Considere la siguiente EDO de segundo orden:

d2 y
dx 2
dy
−1.8 +0.8 y =0 ; de x=0 a x=3 con y ( 0 )=1 ,
dx
dy
dx |
x=0
=0.5
 Resuelva usando el método del punto medio con h=1
 Resuelva con el método clásico de Runge-Kutta de cuarto orden h=1

51. Considere la siguiente EDO de segundo orden:

d2 y
dx 2
3x 2 1 /2
=1.5(e − y ) ;de x=0 a x=1.5 con y ( 0 ) =0 ,
dy
dx | x=0
=1
 Resuelva usando el método de Heun con h=0.5
 Resuelva con el método clásico de Runge-Kutta de cuarto orden h=0.5
52. Considere un tanque de agua cilíndrico como se muestra en la figura. El tanque
se llena desde la parte superior, y el agua fluye fuera del tanque a través de un
tubo conectado en el fondo. La razón de cambio de la altura h del agua está dada
por la siguiente ecuación:
dh
ρA tanque =K 1 +K 2 cos(Ct )−ρA tubo √ 2 gh
dt

π
Para el tanque dado, A tanque=3.13 m2 , Atubo=0.06 m2 , C= , K1=300 kg/s, K2=200 kg/s,
12
ρ=1000 kg/m3, y g=9.81 m/s2. Determine y grafique la altura del agua como una
función del tiempo para 0≤t≤150 s, si en t=0 h=3 m.
 Use también un programa de python.

53. Un pequeño cohete en reposo con un peso inicial de 3000 lbs (incluyendo 2400 lbs de
combustible) se lanza verticalmente. El cohete consume combustible a razón de 80
lbs/s que provee un empuje constante T de 8000 lbs. El peso instantáneo del cohete
es w(t)=3000-80t. La resistencia D que experimenta el cohete está dada por

D=0.005 g ( )
dy 2lbs, donde y es la distancia en pies, y g=32.2 pies/s 2. Usando la
dt
ecuación de movimiento de la segunda ley de Newton:
2
wd y
=T −w−D
g dt 2
Determine y grafique la posición, velocidad, y aceleración del cohete (tres figuras separadas) como función
del tiempo desde t=0, cuando el cohete comienza a moverse hacia arriba desde el reposo, hasta t=3.
Reduzca la EDO de segundo orden a un sistema de ecuaciones de primer orden. Use el método de Runge-
Kutta de cuarto orden para un h=0.05 s.

54. Una masa de 2 lb m está sostenida verticalmente por un resorte con una constante k=1.2 lb/in y un
amortiguador que ejerce una fuerza F D sobre la masa la cual depende del desplazamiento y y la velocidad de

la masa dy/dt de acuerdo con F D =μ ( y −1 )

2
( dydt ) donde µ=4 lb-s/ft . Aplicando la segunda ley de Newton
3

se obtiene la siguiente EDO:

 ( )
2
d y dy
+ μ ( y −1 )
2
m 2
+ky=0
dt dt
Re-escriba la EDO de segundo orden como un sistema de dos EDO de primer orden. Determine el
dy dy
desplazamiento y(t) y la velocidad para 0≤t≤20 s con las condiciones iniciales y(0)=0.75 pies y ↑ =0.
dt dt t =0
Escriba un script de Python usando el método de Runge-Kutta con h=0.05 s.

Tema. Series de Fourier

55. Encuentre la serie de Fourier de la función f(x), y grafique la suma parcial de los primeros diez términos.
a. f ( x )=−1 (−2< x <0 ) , f ( x )=1 ( 0< x <2 ) ,T =4
b. f ( x )=0 (−2< x <0 ) , f ( x )=4 ( 0< x <2 ) , T =4
c. f ( x )=x 2 (−1< x<1 ) , T =2
d. f ( x )=π x 3 /2, T =2
e. f ( x )=sin πx ( 0< x< 1 ) , T =1

f. f ( x )=cos πx ( −12 < x< 12 ) , T =1

g. f ( x )=| x|(−1< x <1 ) ,T =2

h. f ( x )= { 1−x si 0< x <1

1+ x si−1< x <0
, T=2

i. f ( x )=1−x 2 (−1< x <1 ) ,T =2

j. f ( x )=0 (−2< x <0 ) , f ( x )=x ( 0< x <2 ) ,T =4
k. f ( x )=−x (−1< x <0 ) , f ( x )=x ( 0< x<1 ) , f ( x )=1 , T =4

56. Grafique las sumas parciales {SN(x)} de la serie trigonométrica dada. Experimente con diferentes valores de n
y con graficas en intervalos distintos del eje x. Utilice sus gráficas para formular una expresión de forma
cerrada para una función definida por 0 < x < T que esté representada por las series.

( )
∞
−1 4 1 nπ nπ
a) f ( x )= + 2 ∑ 2 1−cos cos x
4 π n=1 n 2 2

[ ]
∞
−π (−1 )n−1 1−2 (−1 )n
b) f ( t )= +∑ cos nx + sin nx
4 n=1 n2 π n

57. Escriba la serie de Fourier compleja de f, trace algunos puntos de espectro de frecuencia.
a. f ( x )=2 x ( 0≤ x <3 ) ,T =3
b. f ( x )=x 2 ( 0 ≤ x <2 ) ,T =2

c. {
f ( x )= 0 ,( 0 ≤ x <1), T =4
1 ,(1 ≤ x <4 )
d. f ( x )=1−x , ( 0 ≤ x<6 ) ,T =6
e. {
f ( x )= −1 ,(0 ≤ x< 2) , T =4
2 ,(2≤ x <4)
−x
f. f ( x )=e ( 0 ≤ x <5 ) ,T =5
g. {
f ( x )= x ,(0 ≤ x <1) , T =2
2−x ,(1 ≤ x <2)
h. Sea f la funcion periodica, parte de cuya grafica se muestra en la siguiente figura. Encuentre la serie
de Fourier compleja de f y trace algunos puntos de su espectro de amplitud.
 Tema. Transformada de Fourier

58. Utilice Matlab para obtener el espectro de frecuencias de la función dada por un número finito de puntos y
despliegue los espectros de línea de Re e Im en función de la frecuencia.
a.
t 0 0.25 0.5 0.75 1.0 1.25 1.5 1.75 2.0
f (t) 0 0.25 0.5 0.75 1.0 1.25 1.5 1.75 2.0
b.
t 0 0.25 0.5 0.75 1.0 1.25 1.5 1.75 2.0
f (t) 0 0.4375 0.75 0.9375 1.0 0.9375 0.75 0.4375 2.0

59. Calcule la transformada de Fourier del pulso exponencial dado por

{
at
f ( x )= e−at,(t ≤ 0) , (a> 0)
e ,(t >0)

60. Determine la transformada de Fourier del pulso “encendido-apagado”

61. Un pulso triangular se define como sigue. Dibuje la funcion y determine la transformada de Fourier. ¿Cuál es
la relación entre el pulso y el del ejercicio anterior.

{ ( )t + A ,(−T ≤ t ≤ 0)
A
f ( t )= T
( −AT ) t + A ,(0<t ≤ T )
62. Determine la transformada de Fourier de las siguientes funciones. Dibuje la función h( t)=f (t)−g (t)y
determine su transformada de Fourier.

{
f ( t )= 2 K ,(|t|≤ 2)
0 ,(|t |>2)

{
g ( t ) = K ,(|t|≤ 1)
0 ,(|t|>1)

63. Calcule la transformada de Fourier del pulso “apagado-encendido-apagado” f (t)definido por

{
0(t ←2)
−1(−2 ≤ t ←1)
f ( t )= 1(−1 ≤ t ≤ 1)
−1(1<t ≤ 2)
0(t >2)