2.1.2.1

// Rappresento troncando a differenti ordini

// le serie di Fourier per un'onda quadrata e un'onda triangolare.

clear;clf;

function f=squared(x, n)

    f=0

    for i=1:2:n

        f=f+(4/%pi)*sin(i*x)/i

    end

endfunction

//

function f=triangular(x, n)

    f=%pi/2

    for i=1:2:n

        f=f-(4/%pi)*cos(i*x)/i^2

    end

endfunction

//

x=0:0.1:3*%pi

//

subplot(2,1,1)

plot(x,squared(x,1),x,squared(x,3),x,squared(x,5),x,squared(x,50),x,squared(x,100))

legend("n=1","n=2","n=3","n=50","n=100")

title("onda quadrata")

//

subplot(2,1,2)

plot(x,triangular(x,1),x,triangular(x,3),x,triangular(x,5),x,triangular(x,50),x,triangular(x,100))

legend("n=1","n=2","n=3","n=50","n=100")

title("onda triangolare")