Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- 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
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement