PRACTICA No.
2
“DESARROLLO DE LA SERIE DE FURIER"
I. OBJETIVO
Haciendo uso de octave u otro software, verificar la serie Trigonométrica y Exponencial de
Fourier y desarrollar los ejercicios propuestos en el cuestionario.
II. EQUIPOS Y MATERIALES
Matlab
Octave
PC Pentium II – actual
III. INTRODUCCION TEORICA
Matlab
Matlab es una herramienta interactiva basada en matrices para cálculos científicos y de
ingeniería (de hecho, el termino matlab procede de matrix laboratory). Desde el punto de vista
del control, matlab se puede considerar un entorno matemático de simulación que puede
utilizarse para modelar y analizar sistemas. Permitir a el estudio de sistemas continuos,
discretos, lineales y no lineales, mediante descripción interna y externa, en el dominio temporal
y frecuencia.
Octave
Octave o GNU Octave es un programa y lenguaje de programación para realizar cálculos numéricos.
Como su nombre indica, Octave es parte del proyecto GNU. Es considerado el
equivalente libre de MATLAB. Entre varias características que comparten, se puede destacar que
ambos ofrecen un intérprete, permitiendo ejecutar órdenes en modo interactivo. Nótese que Octave
no es un sistema de álgebra computacional, como lo es Maxima, sino que está orientado al análisis
numérico.
IV. PROCEDIMIENTO
1. Desarrolle la serie trigonométrica de Fourier de la función:
𝑨 ; 𝒆𝒏 𝟎 ≤ 𝒕 < 𝝅
𝒇(𝒕) = {
−𝑨 ; 𝒆𝒏 𝝅 ≤ 𝒕 ≤ 𝟐𝝅
Grafique la serie de Fourier f(t).
SOLUCION:
La función f(t) es una función impar cuya serie trigonométrica de Fourier es:
4𝐴 1 1
𝑓(𝑡) = ( ) [𝑠𝑒𝑛(𝜔𝑡) + 𝑠𝑒𝑛 (3𝜔𝑡) + 𝑠𝑒𝑛(5𝜔𝑡) … ]
𝜋 3 5
Programando para mostrar la gráfica de la serie de Fourier.
Fs=1000;
t=(1:100)/Fs;
w=2*pi*10;
f=(8/pi)*(sin(w*t)+(1/3)*sin(3*w*t)+(1/5)*sin(5*w*t)+(1/7)*sin(7*w*t)+(1/9)*sin(9*w*t);
plot(t,f)
grid
2. Desarrolle la serie trigonométrica de Fourier, para:
𝝅 𝝅
𝑨 ; 𝒆𝒏 −≤𝒕<
𝟐 𝟐
𝒇(𝒕) = { 𝝅 𝟑𝝅
−𝑨 ; 𝒆𝒏 ≤ 𝒕 ≤
𝟐 𝟐
SOLUCION:
La función f(t) es una función impar cuya serie trigonométrica de Fourier es:
4𝐴 1 1 1 1
𝑓(𝑡) = ( ) [𝑐𝑜𝑠(𝜔𝑡) − 𝑠𝑒𝑛 (𝜔𝑡) + 𝑐𝑜𝑠(5𝜔𝑡) − 𝑐𝑜𝑠(7𝜔𝑡) + 𝑐𝑜𝑠(9𝜔𝑡) … ]
𝜋 3 5 7 9
Cuyo programa es:
Fs=1000;
t=(1:100)/Fs;
w=2*pi*10;
f=(8/pi)*(cos(w*t) -(1/3)*cos(3*w*t) +(1/5)*cos(5*w*t)-(1/7)*cos(7*w*t)
+(1/9)*cos(9*w*t) -(1/11)*cos(11*w*t) +(1/13)*cos(13*w*t) );
plot(t,f)
grid
3. De acuerdo al problema 2, la expresión general de la serie trigonométrica de Fourier
para función f(t) par, está dado por:
∞
𝟒𝑨 𝒏𝝅
𝒇(𝒕) = ( ) ∑ 𝒔𝒆𝒏 ( ) 𝐜𝐨𝐬(𝒏𝝎𝒕)
𝝅 𝟐
𝒏=𝟎
Desarrolle mediante la instrucción de control de flujo:
SOLUCION:
Fs=100;
t=(-200:200)/Fs;
w= 2*pi;
A=2;
f=0
for n=1:1000;
f=f+(4*A/(n*pi))*(sin(n*0.5*pi))*cos(n*w*t);
end;
plot(t,f)
xlabel('t(seg)')
ylabel('AMPLITUD')
title('FUNCION PAR ONDA CUADRADA')
grid
V. CUESTIONARIO
1. Dada la expresión de la serie de Fourier trigonométrica, desarrolle la gráfica de f(t).
Usando el criterio de problema 3. Dada la serie:
∞
𝑨 𝟏
𝒇(𝒕) = − ∑ 𝐬𝐞𝐧(𝒏𝝎𝒕)
𝟐 𝐧
𝒏=𝟎
𝑺𝒊 𝒇(𝒕) = 𝑨𝒕, 𝒕 ∈ (𝟎, 𝟏)
2. Desarrolle la serie exponencial de Fourier, si 𝒇(𝒕) = 𝑨𝒔𝒆𝒏(𝝅𝒕) en el intervalo
[0,1].Grafique la S.E.F.
𝟒𝑨
3. Programe en Octave la siguiente serie trigonométrica 𝒇(𝒕) = ∑∞
𝒏=𝟏 (𝒏𝝅)𝟐 𝒄𝒐𝒔(𝒏𝝎𝒕) con
n=impar de la onda triangular.
4. Grafique la serie exponencial de Fourier de la función 𝒇(𝒕) = 𝑨𝒆−𝟐𝒕 en t [0.1].
