codigosmemo
-
Upload
joel-gomez-reyes -
Category
Documents
-
view
212 -
download
0
description
Transcript of codigosmemo
Codigo 1
%Ejemplo del calculo de una serie de fourier de una funcion continua en un%periodo.clc clear allclose all%Definimos la funcin y la graficamossyms x y=x;a=-2;b=2;L=(b-a)/2;t=a:(b-a)/200:b;g=subs(y,t);subplot(1,3,1)plot(t,g)title('Funcion principal')%Encontramos el valor de a0a0=1/(L)*int(y,a,b);h=a0/2;%Calculamos la serie de fourier con un ciclo for para obtener 20 terminos y en cada ciclo calculado%el valor de ak y bkfor k=1:20w=k*pi/L;ak=1/(L)*int(y*cos(w*x),a,b);bk=1/(L)*int(y*sin(w*x),a,b);h=h+ak*cos(w*x)+bk*sin(w*x);end%Graficamos la expansion en serie de fourierm=subs(h,t);subplot(1,3,2)plot(t,m)title('Serie de Fourier de la funcin')%Graficamos la serie de Fourier mostrando tres periodosc=a-2*L:(b-a)/200:b+2*L;d=subs(h,c);subplot(1,3,3)plot(c,d)title('Tres periodos de la funcin')
cdigo 2%Ejemplo del calculo de una serie de fourier de una funcion a trozoa en un%periodo. La solucin involucra el considerar dos funciones distintasclc clear allclose all%Definimos la funcin y la graficamossyms x ry=x;%r=0;r=subs(r,0); %Hacemos uso del subs, ya que si solo usamos cero, matlab no identifica a%r como funcin sino como variable doublea=-2;c=0;b=2;L=(b-a)/2;%Grafica de la primera parte de la funcint=a:(b-a)/200:c;subplot(1,3,1)plot(t,r)%Grafica de la segunda parte de la funcinhold onl=c:(b-a)/200:b;g=subs(y,t);subplot(1,3,1)plot(l,g)title('Funcion principal')%Encontramos el valor de a0a0=1/L*int(r,a,c)+1/L*int(y,c,b);h=a0/2;%Calculamos la serie de fourier con un ciclo for para obtener 20 terminos y en cada ciclo calculado%el valor de ak y bkfor k=1:20w=k*pi/L;ak=1/(L)*int(r*cos(w*x),a,c)+1/(L)*int(y*cos(w*x),c,b);bk=1/(L)*int(r*sin(w*x),a,c)+1/(L)*int(y*sin(w*x),c,b);h=h+ak*cos(w*x)+bk*sin(w*x);end%Graficamos la expansion en serie de fourierp=a:(b-a)/200:b;m=subs(h,p);subplot(1,3,2)plot(p,m)title('Serie de Fourier de la funcin')%Graficamos la serie de Fourier mostrando tres periodosc=a-2*L:(b-a)/200:b+2*L;d=subs(h,c);subplot(1,3,3)plot(c,d)title('Tres periodos de la funcin')