Términos Serie de fourier
Publicado por falmonacidc (1 intervención) el 08/10/2016 20:53:06
Hola, me encuentro trabajando con la serie de fourier para funciones simples, del tipo y=mx
busqué este algoritmo y hice unos pequeños ajustes.
Logro graficar la función original con la serie de fourier con el determinado número de armónicos de entrada.
Ahora quisiera que me entregara al final los términos de la serie.
Y también a partir de la serie generada lograr graficar el espectro de esta, según yo tendría que ingresar el algoritmo de la FFT, pero para eso necesito poder aislar en cada caso los términos de la serie o por lo menos la función que define la serie.
muchas gracias
Adjunto el algoritmo que llevo creado:
busqué este algoritmo y hice unos pequeños ajustes.
Logro graficar la función original con la serie de fourier con el determinado número de armónicos de entrada.
Ahora quisiera que me entregara al final los términos de la serie.
Y también a partir de la serie generada lograr graficar el espectro de esta, según yo tendría que ingresar el algoritmo de la FFT, pero para eso necesito poder aislar en cada caso los términos de la serie o por lo menos la función que define la serie.
muchas gracias
Adjunto el algoritmo que llevo creado:
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
function fourier
disp('Hola, función y(t)=mt')
m=sym('m')
t=sym('t')
m = input('Ingresar Pendiente: ')
func= m*t
T = input('Ingresar Periodo: ')
arm= input ('Ingrese N° Términos de la Serie: ')
disp('Favor esperar el gráfico')
syms t;
w = 2*pi/T;
n = 1:arm;
a0 = (2/T)*int(func,t, -T/2, T/2);
an = (2/T)*int(func*cos(n*w*t),t, -T/2, T/2);
bn = (2/T)*int(func*sin(n*w*t),t, -T/2, T/2);
f = a0/2 + dot(an,cos(n*w*t)) + dot (bn, sin(n*w*t));
ezplot(func, [-T/2 T/2])
hold on
grid on
ezplot (f, [-T/2 T/2])
end
Valora esta pregunta
0