Series de Fourier (Raúl, Sofía, Jaime)

Trabajo realizado por estudiantes
Título	Series de Fourier
Asignatura	EDP
Curso	2023-24
Autores	Raúl Ortega Sofía Gómez Jaime Sáenz de Miera
Este artículo ha sido escrito por estudiantes como parte de su evaluación en la asignatura

2 Ej1

right

x=linspace(-1,1,200);
y=1/2*ones(1,200);
hold on
plot(x,y)
for n=1:4
    plot(x,cos(n*pi*x))
    plot(x,sin(n*pi*x))
end
plot(x,cos(5*pi*x))
title('Primeros 10 términos de la base trigonométrica');
xlabel('x');
ylabel('y');
legend('1/2', '$\cos(\pi x)$', '$\sin(\pi x)$', '$\cos(2\pi x)$', '$\sin(2\pi x)$', '$\cos(3\pi x)$', '$\sin(3\pi x)$', '$\cos(4\pi x)$', '$\sin(4\pi x)$','$\cos(5\pi x)$', 'Interpreter','latex');

3 Ej2

right

x1=0:10^(-3):1;
x2=-1:10^(-3):0;
figure(1)
hold on
plot(x1,x1.*(1-x1),'r',LineWidth=1.5)
plot(x2,x2.*(1+x2),'b',LineWidth=1.5)
title('Extension impar de $f(x)=x(1-x)$', 'interpreter','latex',FontSize=15)
xlabel('x')
ylabel('y')

right

a=0; b=1;                    
u_1=a:10^(-3):b;
N=length(u_1)-1;
h=(b-a)/N;
w=ones(N+1,1);              
w(1)=1/2; w(N+1)=1/2;
figure(2)
n=[1,5,10];
f_n=zeros(length(n),N+1);
for i=1:length(n)
    for k=1:n(i)
        g=(u_1.*(1-u_1).*sin(k*pi*u_1))'; 
        a_k=2*h*w'*g;
        f_n(i,:)=f_n(i,:)+a_k*sin(k*pi*u_1);
    end
    subplot(1,3,i)
    hold on 
    plot(u_1,f_n(i,:),'r',LineWidth=1.5)
    plot(u_1,u_1.*(1-u_1),'--b',LineWidth=1.5)
    title('n='+string(n(i)))
    legend('$f_n(x)$','$f(x)$','interpreter','latex','Location','northwest')
    grid('on')
end