Dobre

%ZAD 17 wczytywanie i decymajca
[x,fp]=audioread('cartoon008.wav'); %wczytuje plik
y=decimate(x,5);
subplot(331); plot(x); title('17-decymacja'); xlabel('t'); ylabel('A');
fp
%ZAD18 widmo amplitudowe
Nf=1024;
v=fft(x,Nf); %fft dla syg x i liczbie probek=Nf
f=linspace(0,fp,Nf); %wektor o Nf wartosciach podzielonych od 0 do fp
subplot(332); plot(f,v); title('18-widmo amp'); xlabel('czestotliwosc'); ylabel('widmo');
vv=abs(v);
subplot(333); plot(f,vv); title('18A-modul widma'); xlabel('czestotliwosc'); ylabel('mod widma');

%ZAD19 filtr o dlugosci 51 i fg=1400hz
h=fir1(51,[1400/(fp/2)],'low');
q=fft(h,Nf);
q1=abs(q);
subplot(334);plot(f,q1); title('19-mod widma?'); xlabel('czestotliwosc'); ylabel('mod widma');
figure(2); freqz(h,1); title('19');

%ZAD20
y=filter(h,1,x);
rr=fft(y,Nf);
ww=abs(rr);
figure(1);
subplot(335); plot(f,ww); title('20-mod widma'); xlabel('czestotliwosc'); ylabel('mod widma');

%ZAD21 dany sygnal
fp=150; %czestotliwosc probkowania
t1=0:1/fp:2;
f1=5;
y1=sin(2*pi*f1*t1);

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

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

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];
subplot(336); plot(t,y); title('21-zmienny syg'); xlabel('t'); ylabel('A');
figure(3); plot(t,y); title('21-zmienny syg'); xlabel('t'); ylabel('A');

%ZAD22 autokorelacja syg z 21
tau_max=400;
[s,T]=xcorr(y,tau_max);
figure(1);
subplot(337); plot(T,s); title('22-autokorelacja'); xlabel('t'); ylabel('A');

%ZAD23 korelacja z 2 sygnalem
f7=5;
fi7=0; %mozna pominac
t7=0:1/fp:10;
y7=sin(2*pi*f7*t7+fi7);
[ss,TT]=xcorr(y,y7,tau_max);
subplot(338), plot(TT,ss); title('22-korelacja'); xlabel('t'); ylabel('A');