close all
clear all;
fp=48000;
load dane_Lud3;
tk=100;
t=0:1/fp:tk-1/fp;
z1=yt(:,2);

%filtracja

N   = 100;          % FIR filter order
Fp  = 8e3;          % 8 kHz passband-edge frequency
Fs  = 48e3;         % 48 kHz sampling frequency
Rp  = 0.01;         % Corresponds to 0.01 dB peak-to-peak ripple
Rst = 1e-4;         % Corresponds to 80 dB stopband attenuation

lowpassFilt = dsp.LowpassFilter('DesignForMinimumOrder',false, ...
    'FilterOrder',N,'PassbandFrequency',Fp,'SampleRate',Fs,...
    'PassbandRipple',Rp, 'StopbandAttenuation',80);


%z = lowpassFilt(z1);
%plot(z1); hold
%plot(z)

z=z1;

% d�ugo�� transformaty i rozmiar bufora
L = 2^23; 

% szybka transformata Fouriera
Z = abs(fft(z,L))/(L/2);
% generowanie "osi cz�stotliwo�ci"
f = 0:fp/L:fp/2-fp/L; %fp/L - rozdzielczo�� cz�stotliwo�ciowa,fp/2- pasmo analizy

% generowanie wykres�w
figure
plot(f,Z(1:L/2),'LineWidth',2);
xlabel('cz�stotliwo�� [Hz]'); ylabel('przyspieszenie [m/s^2]'); hold on
xlim([0,100])


f_obr=8.68607
orders=f/f_obr;

figure
plot(orders,Z(1:L/2),'LineWidth',2); grid
xlabel('rz�d'); ylabel('przyspieszenie [m/s^2]'); hold on
xlim([0,6])

BPO=3.07
BPI=4.92
BS=2.06
FT=0.38