Código Implementado en Scilab
x=-%pi:0.1:%pi //Definir Dominio i=length(x)//Calculo de longitud; valores obtenidos de x y=cos(x);//Evaluación de lo los valores anteriore plot(x, y) funtion y=a0fserie(x),y=x^2,endfuntion a0=(1/(2*%pi))*intg(0,%pi,a0fserie) funcprot(0) //Calcular coeficientes: a_n for n=1:8 funtion y=f(x),y=x^2*cos(n*x),endfuntion an(n)=(1/%pi)*intg(0,%pi,f) end //Calcular coeficientes: b_n for n=1:8 funtion y=f(x),y=x^2*sin(x),endfuntion bn(n)=(1/%pi)*intg(0,%pi,f) end yfs=a0; for n=1:8 //Calcular todos los armonicos yfsp=a_n*cos(n*x)+bn*sin(n*x) yfs=yfs+yfsp end plot(x,yfs)
Resultado Obtenido
No hay comentarios:
Publicar un comentario