test-widmo1

fs=441;
f=731.17; % czestotliwosc sygnalu
fp=882 % czestotliwosc probkowania musi być przynajmniej dwa razy wieksza,pozwala to uniknac alliasingu i teoretycznie w pelni odwzorowac dany sygnal

N=1000; % liczba probek
A=2;
t=(0:1/N:(N-1)/fp);
x=A*sin(2*pi*f*t); % uzyskany sygnal sinusoidalny
subplot(2,1,1)
plot(t,x)
title 'sygnal sinusoidalny'

subplot(2,1,2)
X=fft(x) %transformata Fouriera dla wskazanego sygnału

y=20*log10(abs(X)) % zamiana na skale decybelowa
plot(fs,y);
title 'Uzyskane widmo sygnalu'
xlabel('dB')