Question

求大神翻译一下matlab代码。

m = mobiledev
m.AccelerationSensorEnabled = 1
m.Logging = 1
pause(15)
m.Logging = 0
[a, t] = accellog(m)
plot(t, a)
legend('X', 'Y', 'Z')
xlabel('Relative time (s)')
ylabel('Acceleration (m/s^2)')
x = a(:,1)
y = a(:,2)
z = a(:,3)
mag = sqrt(sum(x.^2 + y.^2 + z.^2, 2))
plot(t, mag)
xlabel('Time (s)')/
ylabel('Acceleration (m/s^2)')
magNoG = mag - mean(mag)
plot(t, magNoG)
xlabel('Time (s)')
ylabel('Acceleration (m/s^2)')
minPeakHeight = std(magNoG)
[pks, locs] = findpeaks(magNoG, 'MINPEAKHEIGHT', minPeakHeight)
numSteps = numel(pks)
hold on
plot(t(locs), pks, 'r', 'Marker', 'v', 'LineStyle', 'none')
title('Counting Steps')
xlabel('Time (s)')
ylabel('Acceleration Magnitude, No Gravity (m/s^2)')
hold off
m.AccelerationSensorEnabled = 0
clear m

Answer 1

m = mobiledev
获得结构体mobiledev句柄
m.AccelerationSensorEnabled = 1
m.Logging = 1
改变结构体中AccelerationSensorEnabled Logging 为1
pause(15)暂停15秒
m.Logging = 0
改变结构体中Logging 为0
[a, t] = accellog(m)
调用accellog函数，结果为a,t
plot(t, a)
legend('X', 'Y', 'Z')
xlabel('Relative time (s)')
ylabel('Acceleration (m/s^2)')
画图

x = a(:,1)
y = a(:,2)
z = a(:,3)
mag = sqrt(sum(x.^2 + y.^2 + z.^2, 2))
计算mag的值
plot(t, mag)
xlabel('Time (s)')/
ylabel('Acceleration (m/s^2)')
画图

magNoG = mag - mean(mag)
mag的值减去mag平均值
plot(t, magNoG)
xlabel('Time (s)')
ylabel('Acceleration (m/s^2)')
画图
minPeakHeight = std(magNoG)
[pks, locs] = findpeaks(magNoG, 'MINPEAKHEIGHT', minPeakHeight)
numSteps = numel(pks)
按顺序是三次调用函数
hold on
保持当前绘图figure
plot(t(locs), pks, 'r', 'Marker', 'v', 'LineStyle', 'none')
title('Counting Steps')
xlabel('Time (s)')
ylabel('Acceleration Magnitude, No Gravity (m/s^2)')
绘图
hold off
取消当前绘图figure
m.AccelerationSensorEnabled = 0 将m中的AccelerationSensorEnabled 变为0
clear m
清除m
子函数没有，只能这样