Matlab - espectro fourier

   
Vista:

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)
Valora esta pregunta
Me gusta: Está pregunta es útil y esta claraNo me gusta: Está pregunta no esta clara o no es útil
0
Responder
Imágen de perfil de JESUS

espectro fourier

Publicado por JESUS jesus.royeth@gmail.com (749 intervenciones) el 18/04/2015 18:30:06
Hola si tú quieres algo así : 11173326_10205002512035759_4317183127174786278_n entonces solo agrega estas dos líneas a tu código :
1
2
hold on
stem(-C)
Valora esta respuesta
Me gusta: Está respuesta es útil y esta claraNo me gusta: Está respuesta no esta clara o no es útil
1
Comentar