Análisis de espectro basado en ventana FFT, responda antes del día 27.

Fs=1000;

T=1/Fs

Tp=0.1;

N=Tp*Fs;

w=100*pi;

n=1:N;

Xn=24*sin(w*n*T)+12*sin( 3*w*n* T+40*pi/180)+6*sin(5*w*n*T+60*pi/180)+3*sin(7*w*n*T+(80*pi) /180)+1.5*sin(9*w*n*T+120*pi/180);

Xn=Xn/max(abs(Xn));

Xk = fft(Xn,4096);

fk=Fs*(0:4095);

plot(fk,abs(Xk)/max(abs(Xk)));

xlabel('Hz');

ylabel('Magnitud');

title('ventana rectangular');

2)

Fs=1000;

T=1/Fs;

Tp=0.1;

N=Tp*Fs;

w=100*pi;

n=1:N;

Xn=24*sin(w*n*T)+12*sin( 3*w*n*T+40*pi/180)+6*sin(5*w*n*T+60*pi/180)+3*sin(7*w*n*T+(80*pi) /180)+1.5*sin(9*w*n*T+120*pi/180);

Xn=Xn/max(abs(Xn)),wn= hanning(N);

Xn1=Xn.*wn';

Xk=fft(Xn1,4096);

fk=Fs*(0:4095);

plot(fk,abs(Xk)/max(abs(Xk)));

xlabel('Hz');

ylabel('Amplitud');

title('Hanning Window');

El resto de cambios se basan en esto

.