test-widmo

fs=1723;
f=3013;
fp=3446 %przynajmniej 2 razy wieksza
N=1000;
A=1;
t=0:1/N:(N-1)/fp;
x=A*sin(2*pi*f*t); % sygnal sinusoidalny
subplot(2,1,1)
plot(t,x)
title 'Sygnal'
subplot(2,1,2)
X=fft(x) %transformata
y=20*log10(abs(X)) %przejscie do skali decybeli
plot(fs,y);
title 'Widmo'