kol1pps

clear all;
close all;

%ZAD 17 wczytywanie i decymajca
[x,fp]=audioread('emergency005.wav'); %wczytujemy sygnal x i f probkowania fp
xpd=decimate(x,5); %5-cio krotna decymacja syg x
t17=linspace(0,fp,length(x));%tworzymy wektor czasu dla sygnalu.
subplot(331); plot(t17,x); title('17-decymacja'); xlabel('t'); ylabel('sygnal');

%ZAD18 widmo amplitudowe
%widmo amplitudowe to wykres transformaty Fouriera
Nf=2^12;
N18=Nf/2+1;
xptf=fft(x,Nf);
wx=abs(xptf);
f18=linspace(0,fp/2,N18);
subplot(332); plot(f18,wx(1:N18)); title('18-widmo amp'); xlabel('f'); ylabel('widmo');

%ZAD19 filtr o dlugosci 51 i fg=1400hz oraz wykres f-cji transmitancji
M=51;
fg=(1400/(fp/2)); %czest graniczna
h=fir1(M-1,fg); %robie filtr
v=fft(h,Nf); %ft z filtru=transmitancja
wh=abs(v); %modul fcji transmitancji
subplot(333); plot(f18,wh(1:N18));

%ZAD20 filtracja i modul widma syg po filt
y=filter(h,1,x); %robie filtracje
v=fft(y,Nf); %mam widmo
wy=abs(v);  %mod widma
subplot(334);plot(f18,wy(1:N18));

%ZAD21 dany sygnal
fp=150; %czestotliwosc probkowania

t1=0:1/fp:2; %pierwsze 2 sek
f1=5;
y1=sin(2*pi*f1*t1);

t2=2:1/fp:4; %2-4 sek
f2=3.5;
fi2=pi/4;
y2=sin(2*pi*f2*t2+fi2);

t3=4:1/fp:8; %od 4-8 sek
A3=1;
y3=2*rand(1,601)-A3;

t4=8:1/fp:10; %8-10 sek
f4=5;
fi4=pi/2;
y4=sin(2*pi*f4*t4+fi4);

y=[y1 y2 y3 y4]; %łączenie wykresów w jeden
t=[t1 t2 t3 t4];

subplot(335); plot(t,y); title('21-zmienny syg'); xlabel('t'); ylabel('sygnal');

%ZAD 22 autokorelacja syg z 21
tau_max=400;
[s,T]=xcorr(y,tau_max);
subplot(336); plot(T,s); title('22-autokorelacja'); xlabel('t'); ylabel('sygnal');

%ZAD23 korelacja z 2 sygnalem
f7=5;
t7=0:1/fp:10;
y7=sin(2*pi*f7*t7);
[ss,TT]=xcorr(y,y7,tau_max);
subplot(337), plot(TT,ss);  title('23-korelacja'); xlabel('t'); ylabel('sygnal');