Quiero escribir la serie de Fourier para \ $ x (t) = 1 \ $ para \ $ | t | \ leq 1 \ $ y cero de otra manera (que se expande periódicamente a otros intervalos con este período de 1) y compárelo gráficamente con \ $ x (t) \ $, para valores diferentes para \ $ N \ $, donde \ $ x (t) = \ sum _ {- N} ^ {N} a_k e ^ {jk \ Omega_0 t} \ $.
Aquí está el código que escribí, básicamente, el intervalo debería estar entre -2 y 2, pero recibo un error todo el tiempo, espero que alguien pueda corregirme o encontrar una mejor manera de hacerlo.
n=linspace(-2,2,5);
x=abs(n)<=1;
y=symsum(exp(-1i*k*pi/4)*(sin(k*pi/4)/(k*pi))*exp(1i*k*pi*n/2),k,-1,1);
y2=symsum(exp(-1i*k*pi/4)*(sin(k*pi/4)/(k*pi))*exp(1i*k*pi*n/2),k,-3,3);
y3=symsum(exp(-1i*k*pi/4)*(sin(k*pi/4)/(k*pi))*exp(1i*k*pi*n/2),k,-7,7);
y4=symsum(exp(-1i*k*pi/4)*(sin(k*pi/4)/(k*pi))*exp(1i*k*pi*n/2),k,-19,19);
y5=symsum(exp(-1i*k*pi/4)*(sin(k*pi/4)/(k*pi))*exp(1i*k*pi*n/2),k,-79,79);
Recibo el siguiente error:
??? Undefined function or variable 'k'. Error in ==> AS at 3 y=symsum(exp(-1i*k*pi/4)*(sin(k*pi/4)/(k*pi))*exp(1i*k*pi*n/2),k,-1,1);