clear clc close all for m=1:4 % DATOS if m==1 f=@(x)x; % funcion aa=0;bb=3; % intervalo n=12; % numero de iteraciones elseif m==2 f=@(x)(x-pi).^2; % funcion aa=0;bb=2*pi; % intervalo n=8; % numero de iteraciones elseif m==3 f=@(x)square(x); % funcion aa=0;bb=2*pi ; % intervalo n=12; % numero de iteraciones elseif m==4 f=@(x)tripuls(x);% funcion aa=-1;bb=1 ; % intervalo n=12; % numero de iteraciones end % DATOS % T=bb-aa; f1=@(x)f(x+T); f2=@(x)f(x-T); % a=zeros(1,n); b=zeros(1,n); xx=linspace(aa,bb); xx1=linspace(aa-T,bb-T); xx2=linspace(aa+T,bb+T); xxt=linspace(aa-T,bb+T,300); a0=quad(f,aa,bb)/T; s=a0*ones(1,length(xxt)); for k=1:n figure(m) clf hold on fun1=@(x)cos(2*pi*k*x/T).*f(x); fun2=@(x)sin(2*pi*k*x/T).*f(x); a(k)=2*quad(fun1,aa,bb)/T; b(k)=2*quad(fun2,aa,bb)/T; s1=a(k)*cos(2*pi*k*xxt/T)+b(k)*sin(2*pi*k*xxt/T); s=s+s1; plot(xx,f(xx),'linewidth',3) plot(xxt,s,'r','linewidth',3) plot(xxt,s1); legend('f(x)','Serie','Termino','location','Southeast') plot(xx1,f1(xx-T),'linewidth',3) plot(xx2,f2(xx+T),'linewidth',3) plot(xxt,s,'r','linewidth',3) if m==1 title(['f(x)=x en [0,3] k=',num2str(k)],'fontsize',16) plot([0 0],[0,3],'linewidth',3) plot([3 3],[0,3],'linewidth',3) axis([min(xxt) max(xxt) -1 max(f(xx))]) elseif m==2 title(['f(x)=(x-\pi)^2 en [0,2\pi] k=',num2str(k)],'fontsize',16) axis([min(xxt) max(xxt) -4 max(f(xx))]) elseif m==3 title(['Onda cuadrada en [0,2\pi] k=',num2str(k)],'fontsize',16) axis([min(xxt) max(xxt) -1.5 1.5]) elseif m==4 title(['Onda triangular en [-1,1] k=',num2str(k)],'fontsize',16) axis([min(xxt) max(xxt) -0.5 1]) end box on pause end end