dobre1

%(17)Wczytanie sygnału i decymacja:
[x,fp]=wavread('s1.wav');
y=decimate(x,5);
subplot(211);
plot(y);

%2. %(18)Policzenie widma amplitudowego sygnału:
Nf=1024;
v=fft(x,Nf);
w=abs(v);
f=linspace(0,fp,Nf);
subplot(212);
plot(f,w);
ylabel('czestotliwosc[Hz]');
xlabe('modul widma');

%3. %(19)Filtr fir o zadanej długości i f granicznej
h=fir1( 51,[1400
plot(h);
freqz(h,1);
%y=filter(h,1,x);


%4. (20)Filtracja i wyznaczenie modułu widma sygnału po filtracji











%5. (21)Generacja sygnału o zadanej długości i częstotliwości
%fp=150; %próbkowanie
%t1=0:1/fp:2;
f1=5;
y1=sin(t1*2*pi*f1);
t2=2:1/fp:4;
f2=3.5;
fi2=pi/4;
y2=sin(t2*2*pi*f2+fi2);
t3=4:1/fp:8;
A3=1;
y3=2*rand(1,601)-A3; %rozkład normalny szumu
t4=8:1/fp:10;
f4=5;
fi4=pi/2;
y4=sin(2*pi*f4*t4+fi4);
y=[y1 y2 y3 y4];
t=[t1 t2 t3 t4];
plot(t,y);

%6. (22)Policzenie autokorelacji sygnału z 5tki z danym tau_max
[r,T]=xcorr(x,400);
plot(T,r);

%7. (23)Korelacja sygnału z 5tki z zadanym sygnalem sinusoidalnym
f7=5;
fi7=pi;
t7=0:1/fp:10;
y7=sin(2*pi*f7+fi7);
[r,T]=xcorr(y,y7,400);
plot(T,r);