fp=1000;                  % częstotliwość próbkowania
N=100;                    % ilość próbek
t=0:1/fp:(N-1)/fp;        % oś czasu
x=1.2*sin(2*pi*87*t+0.3); % zmienna
subplot(121);
stem(t,x);                % wykres
xlabel('czas[s]');        % oznaczenie osi x
ylabel('sinus');          % oznaczenie osi y
Nf=128;                   % długość transformaty Fouriera (liczba próbek)
N21=Nf/2+1;               % zmienna pomocnicza
v=fft(x,Nf);              %
wx=abs(v);                %
f=linspace(0,fp/2,N21);      % rysowanie funkcji liniowo rozłożonych punktów
subplot(122);
stem(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

fp=128;                  % częstotliwość próbkowania
N=128;                    % ilość próbek
t=0:1/fp:(N-1)/fp;        % oś czasu
x=1.2*sin(2*pi*37*t+0.3); % zmienna
subplot(121);
stem(t,x);                % wykres
xlabel('czas[s]');        % oznaczenie osi x
ylabel('sinus');          % oznaczenie osi y
Nf=128;                   % długość transformaty Fouriera (liczba próbek)
N21=Nf/2+1;               % zmienna pomocnicza
v=fft(x,Nf);              %
wx=abs(v);                %
f=linspace(0,fp/2,N21);      % rysowanie funkcji liniowo rozłożonych punktów
subplot(122);
stem(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

fp=128;                  % częstotliwość próbkowania
N=128;                    % ilość próbek
t=0:1/fp:(N-1)/fp;        % oś czasu
x=1.2*sin(2*pi*37*t+0.3); % zmienna
subplot(121);
stem(t,x);                % wykres
xlabel('czas[s]');        % oznaczenie osi x
ylabel('sinus');          % oznaczenie osi y
Nf=128;                   % długość transformaty Fouriera (liczba próbek)
N21=Nf/2+1;               % zmienna pomocnicza
v=fft(x,Nf);              %
wx=abs(v);                %
f=linspace(0,fp/2,N21);      % rysowanie funkcji liniowo rozłożonych punktów
subplot(122);
stem(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

clear;
[a,fa]=wavread('mbi02becz.wav');
Na=length(a);
dr=4;
N=floor(Na/dr);
for i=1:N
  x(i)=a((i-1)*dr+1);
  end;
fp=fa/dr;
t=0:1/fp:(N-1)/fp;
subplot(321);
plot(t,x);
xlabel('czas[s]');
ylabel('wav');

Nf=2^15;                   % długość transformaty Fouriera (liczba próbek)
N21=Nf/2+1;               % zmienna pomocnicza
v=fft(x,Nf);              %
wx=abs(v);                %
f=linspace(0,fp/2,N21);      % rysowanie funkcji liniowo rozłożonych punktów
subplot(322);
plot(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

return;

fp=128;                  % częstotliwość próbkowania
N=128;                    % ilość próbek
t=0:1/fp:(N-1)/fp;        % oś czasu
x=1.2*sin(2*pi*37*t+0.3); % zmienna
subplot(221);             % liczba wierszy, kolumn, odwołanie do pola
stem(t,x);                % wykres
xlabel('czas[s]');        % oznaczenie osi x
ylabel('sinus');          % oznaczenie osi y

Nf=128;                   % długość transformaty Fouriera (liczba próbek)
N21=Nf/2+1;               % zmienna pomocnicza
v=fft(x,Nf);              %
wx=abs(v);                %
f=linspace(0,fp/2,N21);      % rysowanie funkcji liniowo rozłożonych punktów
subplot(222);
stem(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');

y=randn(1,N);
subplot(223);
plot(t,y);
xlabel('czas[s]');
ylabel('szum');

Nf=128;                   % długość transformaty Fouriera (liczba próbek)
N21=Nf/2+1;               % zmienna pomocnicza
v=fft(y,Nf);              %
wx=abs(v);                %
f=linspace(0,fp/2,N21);      % rysowanie funkcji liniowo rozłożonych punktów
subplot(224);
plot(f,wx(1:N21));
xlabel('czestotliwosc[Hz]');
ylabel('modul widma');