lab 3 -pcps

clc;
close all;
clear all;

[y,Fs]=audioread('C:\Users\macie\OneDrive\Desktop\Lab3\a-11.wav');
[p,f]=pwelch(y,[],[],16384,Fs);


%%%%%%%%%%%%%%%%%%%%%%%  ZAD 1  %%%%%%%%
t=(0:348480-1)/Fs;
figure(1);
plot(1:348480,y);
% xlim([232319 348480]);
% yobc=y(232319:1:348480);
% N=116159;
% f=fft(yobc);
N=length(y);
f=fft(y);
GWM=(abs(f).^2)/(Fs/2);
GWMdebprim=10*log10(GWM);
GWMdeb=GWMdebprim(1:N/2+1);
freq=0:Fs/N:Fs/2;
figure(2);
plot(freq,GWMdeb)

%%%%%%%%%%%%%%%%%%%%%%%  ZAD 2 %%%%%%%%%

%zerujacy
fp=250;
Tp=1/fp;
mnoznikpisku=365;
nom=[1 -2*cos(mnoznikpisku*Tp) 1];
denom=1;
y2a=filter(nom,denom,y);
figure(3);
plot(1:348480,y2a);
% sound(y2a,Fs);
f2a=fft(y2a);
GWM2a=(abs(f2a).^2)/(Fs/2);
GWMdebprim2a=10*log10(GWM2a);
fvtool(nom,denom);
GWMdeb2a=GWMdebprim2a(1:N/2+1);
figure(4);
plot(freq,GWMdeb2a)

%%%%%%% notch
close all;
plot(freq,GWMdeb2a)
wo=11500/(Fs/2);
bw=wo/35;
[btest,atest]=iirnotch(wo,bw);
fvtool(btest,atest);
y3=filter(btest,atest,y);
figure(6);
plot(1:348480,y3);

%%%%%%%sr ruchoma
close all;
y4=movmean(y,10);
plot(1:348480,y4);