Impact location detection from accelerometer data

Question

0 个投票

I have accelerometer data coming from a accelerometer sensor which is mounted on top of handle of a cricket bat, I have played some shots some are missed and some hit, now I want to detect which shots are missed shots and which shots are hit shots. So how I can detect impact location or collide between cricket bat and ball. I guess there should be some algorithm or procedures by which we can get clue of impact location. Do matlab provide any script for it or how I can get help regarding this.

5 个评论
显示 3更早的评论隐藏 3更早的评论

Mathieu NOE 2024-3-4

hello again

can you comment your graph and describes what you habe done to generate that data ? we needto make the correlation between the graph and the events

xubyr 2024-3-8

Hello Mathieu,

I have used a bno055 sensor accelerometer data to plot the graph on the browser and I am getting x,y, & z. I am attaching the image file of plot drawn. I am getting the buffer and converting them into json and saved in file.

I have also attached a data file in json format. This data is generated while playing cricket for one over ie 6 balls, it has hit and missed shots in it.

请先登录，再进行评论。

请先登录，再回答此问题。

Follow Question

Answer 1

Mathieu NOE 2024-3-8

编辑：Mathieu NOE 2024-3-8

在 MATLAB Online 中打开

0 个投票

hello

so this would be my suggestion

I hope my conclusion is not completely wrong !

At the end it seems I can clearly see 3 sharp and high amplitude peaks (above 5 g's ) corresponding to good hits; the other minor peaks are either missed hits or you shaking the baseball

code :

filename = "data1.csv";
data = importfile(filename);
% data = readtable('data1.csv') ; % 
% data format : 
% x,"y","z","timestamp"
% 0.10,"-1.69","9.85","14:58:30:267"
x = data.x;
y = data.y;
z = data.z;
timestamp = data.timestamp;
tmp = split(timestamp,':');
h = str2double(tmp(:,1)); % h column
m = str2double(tmp(:,2)); % m column
s = str2double(tmp(:,3)); % s column
ms = str2double(tmp(:,4)); % ms column
ts = h*3600 + m*60 + s + ms/1e3; % time in seconds 
ts = ts - ts(1);
% plot(diff(ts)) ; % Nb the non constant sampling interval
% resample the data on a constant time interval time vector
tt = linspace(ts(1),ts(end),numel(ts));
x = interp1(ts,x,tt);
y = interp1(ts,y,tt);
z = interp1(ts,z,tt);
dt = mean(diff(tt));
Fs = 1/dt;
figure(1)
plot(tt,x,tt,y,tt,z);
legend('x','y','z')
% some high pass filtering to remove DC / low frequency signals
flow = 1;
[b,a] = butter(2,flow*2/Fs,'high');
xf = filtfilt(b,a,x);
yf = filtfilt(b,a,y);
zf = filtfilt(b,a,z);
figure(2)
plot(tt,yf,tt,zf);
legend('y','z')
figure(3)
% combine x,y,z signals into one rms signal 
% see x direction contribution is neglectable , we can focus on y and z
% components only
rms_xyz = sqrt(xf.^2 + yf.^2 + zf.^2);
rms_yz = sqrt(yf.^2 + zf.^2);
plot(tt,rms_yz,tt,rms_xyz);
title('RMS (X) Y and Z acceleration');
xlabel('time (s)');
ylabel('amplitude (g)');
legend('rms(yz)','rms(xyz)')

3 个评论
显示 1更早的评论隐藏 1更早的评论

xubyr 2024-4-5

Hello,

Thanks for your humble response.

I am still away from the solution and could you please tell how to use the above code in matlab. I have also installed matlab desktop software. I am also trying python to draw some result. Kindly suggest what to do.

Mathieu NOE 2024-4-5

hello again

do you have access to matlab ? if yes , and if you are new to it , I would suggest to look at some tutorials :

MATLAB Onramp | Autoformations en ligne - MATLAB & Simulink (mathworks.com)

all the best

请先登录，再进行评论。

Impact location detection from accelerometer data

5 个评论
显示 3更早的评论隐藏 3更早的评论

回答（1 个）

3 个评论
显示 1更早的评论隐藏 1更早的评论

类别

产品

标签

Community Treasure Hunt

Impact location detection from accelerometer data

5 个评论 显示 3更早的评论 隐藏 3更早的评论

回答（1 个）

3 个评论 显示 1更早的评论 隐藏 1更早的评论

类别

产品

标签

另请参阅

Community Treasure Hunt

5 个评论
显示 3更早的评论隐藏 3更早的评论

3 个评论
显示 1更早的评论隐藏 1更早的评论