function czPodst = dzwiek(filename)
clear all;
close all;
filename = 'mbi04jeden.wav';
[y, fs] = wavread(filename); %wektor scieżki y, fs - czestprobkowania
t=(0:length(y)-1)/fs; % wyskalowane w sekundach
ti =length(y)/fs % czas trwania wav
figure(1);
plot(t, y);
%biore 2 "e", odczytuje z wykresu ok. 17500-22000

ye1 = y(17500:22000);

[r,lag] = xcorr(ye1,1000, 'coeff'); %1000 jako ograniczenie do mneijszej lcizby powtórzeń

%fPowtarzania = 1/0.01;
if(nargout == 0)
	%podglad sygnalu
	figure(2);
	t2 = (1:length(ye1))/fs; %otrzymuje w czasie;
	plot(t2, ye1);

	%rysujemy autokoralacje
	figure(3);
	lag = lag/fs;
	plot(lag,r);



f_podst = mod(fs, 1000);  % ponieważ widmo po spróbkowaniu jest okresowe z okresem Fs
if f_podst > Fs/2
	f_podst = -(Fs-f_podst);	% jeżeli czest. wyliczona jest większa niż Fs/2 to jest to częstotliwość ujemna
	f_podst = abs(f_podst);

%jeszcze musza ogarnac czy dobrze wyskalowalem 2 plot i wyznacyzc f powtarzania
end
endfunction