codigosmemo

3
Codigo 1 %Ejemplo del calculo de una serie de fourier de una funcion continua en un %periodo. clc clear all close all %Definimos la función y la graficamos syms 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 a0 a0=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 bk for k=1:20 w=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 fourier m=subs(h,t); subplot(1,3,2) plot(t,m) title('Serie de Fourier de la función') %Graficamos la serie de Fourier mostrando tres periodos c=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 función') código 2 %Ejemplo del calculo de una serie de fourier de una funcion a trozoa en un

description

matlab

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')