data=load('dowling.dat');
realacc=load('dowling.acc');

SF=512;
order=2;
i=1:length(data);
"filtered_data3 = filtering(data, 512, 3, order);"
"filtered_data20 = filtering(data, 512, 20, order);"

vel3=diff(filtered_data3)*512;
vel20=diff(filtered_data20)*512;

acc3=diff(vel3)*512;
acc20=diff(vel20)*512;

%extend the calculated acceleration by two points to 
%match accelerometer data
ndataacc=length(acc3);
acc3(ndataacc+1)=acc3(ndataacc);
acc3(ndataacc+2)=acc3(ndataacc+1);
acc20(ndataacc+1)=acc20(ndataacc);
acc20(ndataacc+2)=acc20(ndataacc+1);

m=length(acc3);

%plot graphs
i=1:m;
figure(1)
"plot(i,realacc,'-',i,acc3,'k.');"
ylabel('angular acceleration rad/s^2')
xlabel('data point');
figure(2)
"plot(i,realacc,'-',i,acc20,'k.');"
ylabel('angular acceleration rad/s^2')
xlabel('data point');

difference3=abs(realacc-acc3);
difference20=abs(realacc-acc20);

mean3=mean(difference3)
mean20=mean(difference20)