INFORME FINAL N°01: SIMULACION DE LA SERIE DE

FOURIER MEDIANTE EL SOFTWARE MATLAB

I. MOSTRAR Y ANALIZAR LOS GRAFICOS A continuación, se muestra las comparaciones de

OBTENIDOS EN LA SIMULACIÓN los gráficos aproximados con serie de Fourier
MEDIANTE MATLAB. con respecto a la onda triangular ideal (función
asignada).
A nuestro grupo se ha asignado la onda
triangular Serie de Fourier para n=1
.

Fig. 1. Gráfico de la onda triangular impar de amplitud Fig. 2. La grafica que se encuentra en la parte superior
10Vpp y periodo de 20ms. La gráfica fue obtenida (verde) muestra la aproximación en serie de Fourier con un
mediante la ejecución del algoritmo realizado en MatLab. armónico. La gráfica de la parte inferior muestra la
comparación entre la función aproximada (azul) y la onda
asignada (línea punteada).
En Fig.1 observamos el desarrollo de la onda
triangular impar durante un periodo, a priori
podemos decir que esta onda presenta un valor Podemos observar que la señal triangular impar
medio (Vdc) igual a cero, así mismo, al tratarse muestra un desarrollo casi similar al de una onda
de una función impar veremos que el desarrollo senoidal. Para una aproximación mediante la
en serie de Fourier presentará una expansión en serie de Fourier con un armónico, la cual se
senos, es decir bn tomara el valor de cero. visualiza en Fig. 2, podemos notar que la función
aproximada presenta semejanza con la onda
 triangular con la diferencia de que en los picos
aún falta una mejor aproximación.
Serie de Fourier para n=3

Fig. 3. La grafica que se encuentra en la parte superior (verde)

muestra la aproximación en serie de Fourier con tres armónicos. La
gráfica de la parte inferior muestra la comparación entre la función
aproximada (azul) y la onda asignada (línea punteada).

Para una aproximación mediante la serie de

Fourier con tres armónicos, la cual se visualiza
en Fig. 3, podemos notar que la función
aproximada se va asemejando mucho más a la
onda triangular asignada. Cuando superponemos Fig. 6. La grafica que se encuentra en la parte superior
(verde) muestra la aproximación en serie de Fourier con
ambas gráficas podemos observar que la quince armónicos. La gráfica de la parte inferior muestra la
aproximación en lo picos se va mejorando comparación entre la función aproximada (azul) y la onda
asignada (línea punteada).
cuando se aumenta los armónicos.
Para una aproximación mediante la serie de
Serie de Fourier para n=5 Fourier con quince armónicos, podemos notar
que superponiendo la función aproximada sobre
Fig. 4. La grafica que se encuentra en la parte superior (verde)
la función ideal ya casi no se distingue diferencia
muestra la aproximación en serie de Fourier con cinco armónicos. alguna, la mínima diferencia que existe es en la
La gráfica de la parte inferior muestra la comparación entre la zona de los picos de la función triangular, en
función aproximada (azul) y la onda asignada (línea punteada).
pretty(An);
disp('Bn =');
Para una aproximación mediante la serie de pretty(Bn);
Fourier con cinco armónicos, la cual se visualiza %Aproximacion de la Función
en Fig. 4, podemos notar que la función disp('Aproximacion de la funcion');
aproximada va mostrando las pendientes que la disp('mediante serie de Fourier');
función triangular posee, así mismo, la Fou1=0;
aproximación en los picos se sigue mejorando. for i=1:N
Fou1=Fou1+subs((An*cos(n*wo*x)),n,i);
end
Serie de Fourier para n=7 Fou2=0;
for i=1:N
Fou2=Fou2+subs((Bn*sin(n*wo*x)),n,i);
Fig. 5. La grafica que se encuentra en la parte superior
(verde) muestra la aproximación en serie de Fourier con siete end
armónicos. La gráfica de la parte inferior muestra la Serie=0.5*Ao+Fou1+Fou2;
comparación entre la función aproximada (azul) y la onda disp('Serie de Fourier =');
asignada (línea punteada).
disp(Serie);
x=-(2*T):1/Fs:(2*T);
Para una aproximación mediante la serie de figure(3);
Fourier con siete armónicos, podemos notar que subplot(2,1,1);
hL=ezplot(Serie,x);
la función aproximada ya casi presenta la misma set(hL,'Color','g');
pendiente que la onda triangular. grid on;
Serie de Fourier para n=15 xlabel('\bf TIEMPO (seg.)','FontSize',15);
ylabel('\bf AMPLITUD (V)','FontSize',15);
title('\bf APROXIMACION DE "f(t)" CON SERIE
DE FOURIER','FontSize',15);
legend('Función aproximada'); for i=1:length(y)
subplot(2,1,2); Bn=Bn+int(y(i)*sin(n*wo*x),x,A(i),A(i+1));
ezplot(Serie,x);hold on end
grid on; Bn=simplify(2*Bn/T);
xlabel('\bf TIEMPO (seg.)','FontSize',15); disp('___________________________');
ylabel('\bf AMPLITUD (V)','FontSize',15); disp('___________________________');
title('\bf COMPARACIÓN','FontSize',15); syms n x
,15); disp('RESULTADOS');
title('\bf ONDA TRIANGULAR IMPAR disp('Los Coeficientes de Fourier');
"f(t)"','FontSize',15); disp('Son los siguientes:');
legend('V vs. t'); disp('Ao =');
set(gca,'FontSize',10) pretty(Ao);
%Espectro de frecuencias disp('An =');
L=length(f); pretty(An);
freq=Fs*(0:(L/2))/L; disp('Bn =');
Y=fft(f); pretty(Bn);
P2=abs(Y/L); %Aproximacion de la Función
P1=P2(1:(L-1)/2+1)*2; disp('Aproximacion de la funcion');
p=0:1:10; disp('mediante serie de Fourier');
x2=abs(8*(Vpp/2)./(p.^2.*pi^2).*sin(p.*(pi/2)));
figure(2); CON SERIE DE FOURIER','FontSize',15);
subplot(2,1,1); legend('Función aproximada');
plot(freq,P1,'r','linewidth',2.2) subplot(2,1,2);
grid on ezplot(Serie,x);hold on
xlabel('\bf FRECUENCIA (Hz)'); grid on;
ylabel('\bf |Cn|'); xlabel('\bf TIEMPO (seg.)','FontSize',15);
title('\bf ESPECTRO DE FRECUENCIA ylabel('\bf AMPLITUD (V)','FontSize',15);
DISCRETA'); title('\bf COMPARACIÓN','FontSize',15);
legend('|F(f)| vs. f'); plot(t,f,'--k');
subplot(2,1,2); disp('___________________________');
stem(p,x2,'linewidth',2) disp('___________________________');
grid on disp('¿Desea Cambiar parametros de la señal?');
xlabel('\bf VALOR DE "n"'); disp('Solo ingrese números');
ylabel('\bf AMPLITUD'); disp('Si=1');
%Coeficiente de Fourier disp('No=2');
syms n x res=input('\n Respuesta (SI(1)/NO(2)) = ');
A=[-T/2 -T/4 T/4 T/2]; if res==2
y=[(-2*Vpp*x/T)-Vpp 2*Vpp*x/T disp('Fin de programa')
-2*Vpp*x/T+Vpp]; else
y=sym(y); disp ('syntaxis incorrecta, Fin de Programa')
Per=max(A)-min(A); end
wo=2*pi/Per; end
Ao=0; elseif res==2
for i=1:length(y) disp('Fin de programa')
Ao=Ao+int(y(i),x,A(i),A(i+1)); else
end disp('syntaxis incorrecta, Fin de Programa')
Ao=simplify(2*Ao/(T)); end

An=0;
for i=1:length(y) II. DIFERENCIAS DE VALOR TEORICO
An=An+int(y(i)*cos(n*wo*x),x,A(i),A(i+1)); CON VALOR EXPERIMENTAL
end OBTENIDO POR SIMULACIÓN
An=simplify(2*An/T);
Bn=0;
 A continuación realizaremos los cálculos
matemáticos. Consideramos la señal triangular
Para el desarrollo en serie de Fourier, derivamos
impar: la serie:
∞
nπ π
f ' ( t ) =∑ b n
∞
n=1
(
( )) m
cos n t
m
2 2
n π π
f ' ' (t )=∑ (−b sen ( n t ))
n 2
n=1 m m

Entonces tenemos:
Definimos la onda mediante funciones m
n2 π 2 1 −2 k 4k m 4k m 2k
singulares: −b n 2 = ∫
m m −m m ( δ ( t+m ) +
m
δ t+ −
2 m
δ t− + δ ( t
2 m ( ) ( )
−2 k 4k m 4k m 2k
f ( t )=
m
μ−2 ( t +m )+
m
μ−2 t+ −
2 m ( )
μ−2 t− + μ−2 (t2−m
2
n
( )
)
m−b n π 2=2 ksen ( nπ )−4 ksen n π −4 ksen n π +2 ksen
( ) ( )
2 2
8k π
Se tiene que f(t) es una función periódica impar b= n sen ( n )
con periodo 2m, por lo tanto, su desarrollo en n π 2
2 2

serie de Fourier consta de términos de seno

solamente, es decir:
8k π 1 3π 1 5π
∞
f ( t )=∑ ( bn sen ( n w0 t ) )
f ( t )= 2 2
n π
sen
m ( )
t − 2 sen
3 m
t + 2 sen
5 m ( )
t …… ( )
n =1

Como podemos en el desarrollo de la serie de

Fourier solo estarán presente los armónicos
a 0=0 impares. Para nuestro caso tenemos:
a n=0 k=5
m=10ms

Donde: w 0=2 π /T =π /m. Entonces:

Para el cálculo de los coeficientes de Fourier 40 π (1)

desarrollamos lo siguiente:
b n= 2 2
n π
sen n
2 ( )
−2 k 4k m 4k m 2k 40 ) 1 1
f ' ( t )=
m
μ−1 ( t +m )+
m
μ−1 t+ −
2 m ( )
μ−1 t− + f ( tμ)−1
2 m ( )
= ( t−m
π 2
3 (
sen ( 100 πt )− 2 sen ( 300 πt ) + 2 sen ( 500 πt ) … …
5
(2)

−2 k 4k m 4k m 2 k Ahora procederemos a mostrar los resultados que

f ' ' (t )=
m
δ ( t+ m ) +
m
δ t+ −
2 m( )
δ t− +
2 m ( )
δ ( t −m )
nos arroja la ejecución del programa y lo
compararemos con lo obtenido mediante cálculos 40
matemáticos. f ( t )= sen ( 100 πt )
π2
Según algoritmo:

Como podemos ver ambos métodos nos arroja el

mismo resultado.

 Para n=3
Según cálculos matemáticos:

40 1
Fig. 8. Coeficientes de Fourier obtenidos mediante el
algoritmo en MATLAB.
f ( t )=
π 2 (
sen ( 100 πt )− sen ( 300 πt )
9 )
Según algoritmo:

Como podemos observar en Fig. 8, el coeficiente

de Fourier “bn” obtenido es:

b n=
−40
2 2
n π ( )( cos(n π2 )−1)
sen n
π
2
El cual es equivalente al obtenido mediante
cálculos (1), con la diferencia de que lo obtenido
mediante el algoritmo presenta el termino  Para n=5

π Según cálculos matemáticos:

( ( ) )
− cos n
2
−1
40 1 1
Donde podemos ver que se puede simplificar
f ( t )=
π 2 (
sen ( 100 πt )− sen ( 300 πt ) + sen ( 500 πt )
9 25 )
π Según algoritmo:
( ( ) )
− cos n
2
−1 =1

Para el caso de las aproximaciones mediante la

serie de Fourier, procederemos a comparar lo
obtenido mediante cálculos matemáticos con lo
obtenido mediante el algoritmo.

 Para n=1
De lo anteriormente mostrado para las
Según cálculos matemáticos:
aproximaciones mediante series de Fourier,
podemos ver que los resultados tanto obtenidos
 mediante cálculos como los obtenidos por el 40
algoritmo son iguales.
f ( t )= ( sen (100 πt ) )
π2
Evaluaremos en t1=8ms y t2=5ms:
III. RESPUESTAS A PREGUNTAS
f ( t 1 )=2.38
1. De lo observado en el computador y del
desarrollo teórico, aproxime el ancho de f ( t 2 )=4.05
banda de la señal que se le asigno.
Aproximación Onda
El ancho de banda absoluto es f2 – f1 , donde t Error
n=1 Triangular
el espectro es cero fuera del intervalo f 1<f< f2
8 ms 2.38 2 0.38
a lo largo del eje positivo de frecuencia. 5 ms 4.05 5 0.95
Para nuestro caso, de acuerdo con lo
desarrollado en el aparatado I. y observando
Fig. 9 (Espectro discreto de amplitud)  Para n=3
notamos que el intervalo que cumple con la
40 1
definición anteriormente mencionada es: f ( t )=
π 2 (
sen ( 100 πt )− sen ( 300 πt )
9 )
0< f <250 Hz
Evaluaremos en t1=8ms y t2=5ms:
0< n<5 armónicos
f ( t 1 )=1.95
Por lo tanto, el ancho de banda de la señal
asignada es 250 Hz. f ( t 2 )=4.50

2. ¿Qué diferencias observa entre la función

ideal y la obtenida por la serie de Fourier? Aproximación Onda
t Error
n=3 Triangular
Gráficamente podemos observar, de acuerdo
8 ms 1.95 2 0.05
con lo desarrollado en el apartado I., que la
5 ms 4.50 5 0.5
función obtenida por la serie de Fourier
presenta una buena aproximación a la función
ideal a partir del armónico siete, pero que la  Para n=5
aproximación en los picos de la onda
triangular no sigue siendo muy buena (ver 40 1 1
Fig. 7). Debido a que la función que se está
f ( t )=
π 2 (
sen ( 100 πt )− sen ( 300 πt ) + sen ( 500 πt )
9 25 )
desarrollando en serie de Fourier no presenta
señales de variación rápida (discontinuidades Evaluaremos en t1=8ms y t2=5ms:
como lo podría presentar la onda cuadrada, f ( t 1 )=1.95
por ejemplo) , no se presenta el fenómeno de
Gibbs.

A continuación, evaluaremos los valores

obtenidos con las aproximaciones de la serie
de Fourier con respecto a los valores de la
onda triangular asignada:

 Para n=1