espectro fourier
Publicado por osvaldo (2 intervenciones) el 17/04/2015 18:27:34
hola estoy realizando el espectro de la serie de fourier pero necesito graficarla simetricamente (-5 a 5 por ejemplo) pero solo grafico la parte positiva en el eje x como puedo graficar la parte simetrica aqui el codigo
fs=10000; %frecuencia de muestreo
p=input('ingrese periodo de la señal ');
n=input('Ingresa número de armonicos ');
t=linspace(-10,10,20);%vector tiempo
f1=0.5*rectpuls(t+0.39025,2.3485);%tramo de función por partes
f2=0.48305*rectpuls(t,1.5657);%tramo de función por partes
f3=f1+f2;%función por partes
w0=2*pi/p;
a0=(1/p)*sum(f3*(1/fs)); %Valor medio de la senal
F=a0; %Incializamos la serie con el valor medio de la senal
for n=1:10 %Numero de armonicos
an=(1/(n*pi))*(0.48305*sin(0.19554*pi*n)+(0.5)*sin(n*0.1176*pi)+0.5*sin(n*0.3129*pi)+0.48305*sin(n*0.1568*pi)); %Calculo de coeficiente an
bn=(1/(n*pi))*(0.48305*cos(0.19554*pi*n)-(0.5)*cos(n*0.1176*pi)+0.5*cos(n*0.3129*pi)-0.48305*cos(n*0.1568*pi)); %Calculo del coeficiente bn
%F1=sum((an*cos(n*w0*t))+bn*sin(n*w0*t)); %Serie trigonometrica de Fourier
C(n)=sqrt((an.^2)+(bn.^2));
end
stem(C)
fs=10000; %frecuencia de muestreo
p=input('ingrese periodo de la señal ');
n=input('Ingresa número de armonicos ');
t=linspace(-10,10,20);%vector tiempo
f1=0.5*rectpuls(t+0.39025,2.3485);%tramo de función por partes
f2=0.48305*rectpuls(t,1.5657);%tramo de función por partes
f3=f1+f2;%función por partes
w0=2*pi/p;
a0=(1/p)*sum(f3*(1/fs)); %Valor medio de la senal
F=a0; %Incializamos la serie con el valor medio de la senal
for n=1:10 %Numero de armonicos
an=(1/(n*pi))*(0.48305*sin(0.19554*pi*n)+(0.5)*sin(n*0.1176*pi)+0.5*sin(n*0.3129*pi)+0.48305*sin(n*0.1568*pi)); %Calculo de coeficiente an
bn=(1/(n*pi))*(0.48305*cos(0.19554*pi*n)-(0.5)*cos(n*0.1176*pi)+0.5*cos(n*0.3129*pi)-0.48305*cos(n*0.1568*pi)); %Calculo del coeficiente bn
%F1=sum((an*cos(n*w0*t))+bn*sin(n*w0*t)); %Serie trigonometrica de Fourier
C(n)=sqrt((an.^2)+(bn.^2));
end
stem(C)
Valora esta pregunta
0