%Dla wskazanej przez prowadzącego częstotliwości %próbkowania fs i wybranej częstotliwości sygnału %prostokątnego (lub sinusoidalnego):
%1) wyliczyć częstotliwość (podstawową, tzn. po procesie %próbkowania) sygnału spróbkowanego
%2) wyznaczyć widmo sygnału spróbkowanego
%3) narysować widmo amplitudowe i sprawdzić czy jest %zgodne z wyliczeniami z pkt. 1) - uwaga - Bardzo ważny %jest opis %osi OX w zakresie (−fs/2:fs/2) i skala [dB]

N=1000;
fSampling=441;
F0=731.17;
A=2;
t=[0:1/fSampling:(N-1)/fSampling];
xt=A*sin(2*pi*F0*t);
%plot(t,xt)
Fn=fSampling/2;

fBasic = mod(F0,fSampling) %dodatnia czestotliwosc podstawowa
fBasicTrue = fBasic-fSampling %oczekiwana czestotliwosc podstawowa

X=fft(xt,1024); %widmo sygnału spróbkowanego
X1=fft(xt,512);
dB=20*log10(abs(X)); %obliczenie skali decybelowej
AmpX=fftshift(dB); %widmo amplitudowe, przesunięcie widma
AmpX1=fftshift(20*log10(abs(X1)));
fAxis=-Fn:fSampling/1024:Fn-fSampling/1024; %oś częstotliwości
fAxis1=-Fn:fSampling/512:Fn-fSampling/512;
plot(fAxis,AmpX,fAxis1,AmpX1)

%czestotliwosc podstawowa zgadza sie z oczekiwana