detect_zero_vel.c 18 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863
  1. //#include "detect_zero_vel.h"
  2. //#include "hal_mt.h"
  3. //#include "hal_flash.h"
  4. //#define MAG_THRESHHOLD 1000.f
  5. //int front_zero_tmp = 0;
  6. //int back_zero_tmp = 0;
  7. //void start_cal_step(int16_t front[3], int16_t back[3], int16_t acc[3])
  8. //{
  9. // int16_t front_zero = 0;
  10. // int16_t back_zero = 0;
  11. // int16_t acc_zero = 0;
  12. //
  13. // detect_zero_vel( front, back, acc, &front_zero, &back_zero, &acc_zero);
  14. //}
  15. //void cal_step(int16_t front_zero, int16_t back_zero)
  16. //{
  17. // static int step_calm_down;
  18. //
  19. // if(front_zero || back_zero)
  20. // {
  21. // if(step_calm_down == 0)
  22. // {
  23. // mFlash.mStep.stepCur[0] ++;
  24. // }
  25. // step_calm_down = 10;
  26. // }
  27. //
  28. // if(step_calm_down > 0)
  29. // {
  30. // step_calm_down --;
  31. // }
  32. //}
  33. // void detect_zero_vel(int16_t front[3], int16_t back[3], int16_t acc[3],
  34. // int16_t *front_zero, int16_t *back_zero, int16_t *acc_zero)
  35. //{
  36. // static float front_mag_window[WINDOW_SIZE];
  37. //
  38. // static float back_mag_window[WINDOW_SIZE];
  39. //
  40. // static float acc_x_window[WINDOW_SIZE];
  41. // static float acc_y_window[WINDOW_SIZE];
  42. // static float acc_z_window[WINDOW_SIZE];
  43. //
  44. // static int last_front_zupt;
  45. //
  46. // static int last_back_zupt;
  47. //
  48. // static int front_zupt_wait;
  49. //
  50. // static int back_zupt_wait;
  51. //
  52. //
  53. // float front_val = sqrt((float) (front[0] * front[0] + front[1] * front[1] + front[2] * front[2]));
  54. //
  55. // float back_val = sqrt((float) (back[0] * back[0] + back[1] * back[1] + back[2] * back[2]));
  56. //
  57. //
  58. // //滑动窗口更新数据
  59. //
  60. // memcpy(front_mag_window, front_mag_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  61. //
  62. // memcpy(back_mag_window, back_mag_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  63. //
  64. // memcpy(acc_x_window, acc_x_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  65. //
  66. // memcpy(acc_y_window, acc_y_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  67. //
  68. // memcpy(acc_z_window, acc_z_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  69. //
  70. // front_mag_window[WINDOW_SIZE - 1] = front_val;
  71. //
  72. // back_mag_window[WINDOW_SIZE - 1] = back_val;
  73. //
  74. // acc_x_window[WINDOW_SIZE - 1] = acc[0] / 2048.f;
  75. // acc_y_window[WINDOW_SIZE - 1] = acc[1] / 2048.f;
  76. // acc_z_window[WINDOW_SIZE - 1] = acc[2] / 2048.f;
  77. // float front_max_val = front_val;
  78. // float front_min_val = front_val;
  79. //
  80. // uint16_t front_max_index = WINDOW_SIZE - 1;
  81. // uint16_t front_min_index = WINDOW_SIZE - 1;
  82. //
  83. // float back_max_val = back_val;
  84. // float back_min_val = back_val;
  85. //
  86. // uint16_t back_max_index = WINDOW_SIZE - 1;
  87. // uint16_t back_min_index = WINDOW_SIZE - 1;
  88. //
  89. // for(int i = 0 ; i < WINDOW_SIZE; i++ )
  90. // {
  91. // if(front_mag_window[i] > front_max_val)
  92. // {
  93. // front_max_val = front_mag_window[i];
  94. // front_max_index = i;
  95. // }
  96. //
  97. // if(front_mag_window[i] < front_min_val)
  98. // {
  99. // front_min_val = front_mag_window[i];
  100. // front_min_index = i;
  101. // }
  102. //
  103. // if(back_mag_window[i] > back_max_val)
  104. // {
  105. // back_max_val = back_mag_window[i];
  106. // back_max_index = i;
  107. // }
  108. //
  109. // if(back_mag_window[i] < back_min_val)
  110. // {
  111. // back_min_val = back_mag_window[i];
  112. // back_min_index = i;
  113. // }
  114. //
  115. // }
  116. //
  117. // /*
  118. // * 计算稳定的状态
  119. // */
  120. //
  121. // float acc_max_val_x = acc_x_window[WINDOW_SIZE - 1] ;
  122. // float acc_min_val_x = acc_x_window[WINDOW_SIZE - 1] ;
  123. //
  124. // float acc_max_val_y = acc_y_window[WINDOW_SIZE - 1] ;
  125. // float acc_min_val_y = acc_y_window[WINDOW_SIZE - 1] ;
  126. //
  127. // float acc_max_val_z = acc_z_window[WINDOW_SIZE - 1];
  128. // float acc_min_val_z = acc_z_window[WINDOW_SIZE - 1];
  129. //
  130. // for(int i = WINDOW_SIZE - 5 ; i < WINDOW_SIZE; i++ )
  131. // {
  132. // if(acc_x_window[i] < acc_min_val_x)
  133. // {
  134. // acc_min_val_x = acc_x_window[i];
  135. // }
  136. //
  137. // if(acc_x_window[i] > acc_max_val_x)
  138. // {
  139. // acc_max_val_x = acc_x_window[i];
  140. // }
  141. //
  142. // if(acc_y_window[i] < acc_min_val_y)
  143. // {
  144. // acc_min_val_y = acc_y_window[i];
  145. // }
  146. //
  147. // if(acc_y_window[i] > acc_max_val_y)
  148. // {
  149. // acc_max_val_y = acc_y_window[i];
  150. // }
  151. //
  152. // if(acc_z_window[i] < acc_min_val_z)
  153. // {
  154. // acc_min_val_z = acc_z_window[i];
  155. // }
  156. //
  157. // if(acc_z_window[i] > acc_max_val_z)
  158. // {
  159. // acc_max_val_z = acc_z_window[i];
  160. // }
  161. //
  162. // }
  163. //
  164. // if(front_max_index > front_min_index && front_max_val > front_min_val + MAG_THRESHHOLD)
  165. // {
  166. // front_zero_tmp = 1;
  167. //
  168. // front_zupt_wait = 10;
  169. // }
  170. // else if(front_max_index < front_min_index && front_max_val > front_min_val + MAG_THRESHHOLD)
  171. // {
  172. // front_zero_tmp = 0;
  173. // }
  174. // else if(front_zero_tmp == 1)
  175. // {
  176. // front_zero_tmp = 2;
  177. // }
  178. //
  179. //
  180. // /*
  181. // * 判断后脚往下压
  182. // */
  183. //
  184. // if(back_max_index > back_min_index && back_max_val > back_min_val + MAG_THRESHHOLD)
  185. // {
  186. // back_zero_tmp = 1;
  187. //
  188. // back_zupt_wait = 10;
  189. // }
  190. //
  191. // else if(back_max_index < back_min_index && back_max_val > back_min_val + MAG_THRESHHOLD)
  192. // {
  193. // back_zero_tmp = 0;
  194. // }
  195. // else if(back_zero_tmp == 1)
  196. // {
  197. // back_zero_tmp = 2;
  198. // }
  199. //
  200. // /*
  201. // * 判断仅前踮 或 后垫
  202. // */
  203. // if(back_zero_tmp == 2 && front_zero_tmp == 0)
  204. // {
  205. // back_zero_tmp = 0;
  206. // }
  207. //
  208. // if(front_zero_tmp == 2 && back_zero_tmp == 0)
  209. // {
  210. // front_zero_tmp = 0;
  211. // }
  212. //
  213. // //延迟等待稳定状态
  214. //// if(*front_zero == 0 && front_zupt_wait > 0 && acc_max_val - acc_min_val < 0.08f)
  215. //// {
  216. //// *front_zero = 1;
  217. //// }
  218. //
  219. // if(front_zero_tmp == 1)
  220. // {
  221. // *front_zero = 1;
  222. //
  223. // }
  224. // else if(front_zupt_wait > 0 && acc_max_val_x - acc_min_val_x < 0.07f && acc_max_val_y - acc_min_val_y < 0.07f && acc_max_val_z - acc_min_val_z < 0.07f)
  225. // {
  226. // *front_zero = 1;
  227. // front_zupt_wait = 20;
  228. // }
  229. // else
  230. // {
  231. // *front_zero = 0;
  232. // }
  233. //
  234. //
  235. //// if(*back_zero == 0 && back_zupt_wait > 0 && acc_max_val - acc_min_val < 0.08f)
  236. //// {
  237. //// *back_zero = 1;
  238. //// }
  239. //
  240. // if(back_zero_tmp == 1)
  241. // {
  242. // *back_zero = 1;
  243. // }
  244. // else if(back_zupt_wait > 0 && acc_max_val_x - acc_min_val_x < 0.07f && acc_max_val_y - acc_min_val_y < 0.07f && acc_max_val_z - acc_min_val_z < 0.07f )
  245. // {
  246. // *back_zero = 1;
  247. //
  248. // back_zupt_wait = 20;
  249. // }
  250. // else
  251. // {
  252. // *back_zero = 0;
  253. // }
  254. //
  255. //
  256. // for(int i = 0 ; i < WINDOW_SIZE - 5; i++ )
  257. // {
  258. // if(acc_x_window[i] < acc_min_val_x)
  259. // {
  260. // acc_min_val_x = acc_x_window[i];
  261. // }
  262. //
  263. // if(acc_x_window[i] > acc_max_val_x)
  264. // {
  265. // acc_max_val_x = acc_x_window[i];
  266. // }
  267. //
  268. // if(acc_y_window[i] < acc_min_val_y)
  269. // {
  270. // acc_min_val_y = acc_y_window[i];
  271. // }
  272. //
  273. // if(acc_y_window[i] > acc_max_val_y)
  274. // {
  275. // acc_max_val_y = acc_y_window[i];
  276. // }
  277. //
  278. // if(acc_z_window[i] < acc_min_val_z)
  279. // {
  280. // acc_min_val_z = acc_z_window[i];
  281. // }
  282. //
  283. // if(acc_z_window[i] > acc_max_val_z)
  284. // {
  285. // acc_max_val_z = acc_z_window[i];
  286. // }
  287. //
  288. // }
  289. //
  290. //
  291. // if(acc_max_val_x - acc_min_val_x < 0.05f && acc_max_val_y - acc_min_val_y < 0.05f && acc_max_val_z - acc_min_val_z < 0.05f)
  292. // {
  293. // *acc_zero = 1;
  294. // }
  295. // else
  296. // {
  297. // *acc_zero = 0;
  298. // }
  299. //
  300. //
  301. // if(*front_zero == 0)
  302. // {
  303. // *front_zero = *back_zero;
  304. // }
  305. //
  306. // if(front_zupt_wait > 0)
  307. // {
  308. // front_zupt_wait --;
  309. // }
  310. //
  311. // if(back_zupt_wait > 0)
  312. // {
  313. // back_zupt_wait --;
  314. // }
  315. //
  316. //
  317. // //利用加速度延续
  318. // if((last_front_zupt == 1|| last_back_zupt == 1)
  319. // && fabsf(acc_x_window[WINDOW_SIZE - 2] - acc_x_window[WINDOW_SIZE - 1]) < 0.05f
  320. // && fabsf(acc_y_window[WINDOW_SIZE - 2] - acc_y_window[WINDOW_SIZE - 1]) < 0.05f
  321. // && fabsf(acc_z_window[WINDOW_SIZE - 2] - acc_z_window[WINDOW_SIZE - 1]) < 0.05f)
  322. // {
  323. // *front_zero = 1;
  324. // *back_zero = 1;
  325. //
  326. // front_zero_tmp = 1;
  327. // back_zero_tmp = 1;
  328. //
  329. // front_zupt_wait = 15;
  330. // back_zupt_wait = 15;
  331. // }
  332. //
  333. //
  334. // last_front_zupt = front_zero_tmp;
  335. //
  336. // last_back_zupt = back_zero_tmp;
  337. //
  338. // /*
  339. // * 直接用这个来记步好了
  340. // */
  341. // cal_step(*front_zero, *back_zero );
  342. //}
  343. #include "detect_zero_vel.h"
  344. #include "hal_mt.h"
  345. #include "hal_flash.h"
  346. #include "hal_imu.h"
  347. #define MAG_THRESHHOLD 1500.f
  348. int front_zero_tmp = 0;
  349. int back_zero_tmp = 0;
  350. void start_cal_step(int16_t front[3], int16_t back[3], int16_t acc[3])
  351. {
  352. int16_t front_zero = 0;
  353. int16_t back_zero = 0;
  354. int16_t acc_zero = 0;
  355. detect_zero_vel( front, back, acc, &front_zero, &back_zero, &acc_zero);
  356. }
  357. void cal_step(int16_t front_zero, int16_t back_zero, float acc_x, float acc_y, float acc_z)
  358. {
  359. static int step_calm_down;
  360. static float acc_min[3];
  361. static float acc_max[3];
  362. if(!(front_zero || back_zero))
  363. {
  364. acc_max[2] = 1.0f;
  365. acc_min[2] = 1.0f;
  366. if(acc_x > acc_max[0])
  367. {
  368. acc_max[0] = acc_x;
  369. }
  370. if(acc_y > acc_max[1])
  371. {
  372. acc_max[1] = acc_y;
  373. }
  374. if(acc_z > acc_max[2])
  375. {
  376. acc_max[2] = acc_z;
  377. }
  378. if(acc_x < acc_min[0])
  379. {
  380. acc_min[0] = acc_x;
  381. }
  382. if(acc_y < acc_min[1])
  383. {
  384. acc_min[1] = acc_y;
  385. }
  386. if(acc_z < acc_min[2])
  387. {
  388. acc_min[2] = acc_z;
  389. }
  390. }
  391. if(front_zero || back_zero)
  392. {
  393. if(step_calm_down == 0 && (acc_max[0] - acc_min[0] > 0.8f || acc_max[1] - acc_min[1] > 0.8f || acc_max[2] - acc_min[2] > 0.8f))
  394. {
  395. mFlash.mStep.stepCur[0] ++;
  396. }
  397. //假设20ms相见为极限
  398. step_calm_down = 5;
  399. acc_max[0] = acc_x;
  400. acc_max[1] = acc_y;
  401. acc_max[2] = acc_z;
  402. acc_min[0] = acc_x;
  403. acc_min[1] = acc_y;
  404. acc_min[2] = acc_z;
  405. }
  406. if(step_calm_down > 0)
  407. {
  408. step_calm_down --;
  409. }
  410. }
  411. void setLongTimeUpTrend(int max_index, int min_index, float max_val, float min_val, float *longTime_maxVal, float *longTime_minVal)
  412. {
  413. if(max_index > min_index && max_val > min_val + 800.f)
  414. {
  415. *longTime_maxVal = *longTime_maxVal > max_val ? *longTime_maxVal : max_val;
  416. *longTime_minVal = *longTime_minVal < min_val ? *longTime_minVal : min_val;
  417. }
  418. else if(max_index < min_index && max_val > min_val + 800.f)
  419. {
  420. *longTime_maxVal = 0.0f;
  421. *longTime_minVal = 50000.0f;
  422. }
  423. }
  424. void setLongTimeDownTrend(int max_index, int min_index, float max_val, float min_val, float *longTime_maxVal, float *longTime_minVal)
  425. {
  426. if(max_index > min_index && max_val > min_val + 800.f)
  427. {
  428. *longTime_maxVal = 0.0f;
  429. *longTime_minVal = 50000.0f;
  430. }
  431. else if(max_index < min_index && max_val > min_val + 800.f)
  432. {
  433. *longTime_maxVal = *longTime_maxVal > max_val ? *longTime_maxVal : max_val;
  434. *longTime_minVal = *longTime_minVal < min_val ? *longTime_minVal : min_val;
  435. }
  436. }
  437. void detect_zero_vel(int16_t front[3], int16_t back[3], int16_t acc[3],
  438. int16_t *front_zero, int16_t *back_zero, int16_t *acc_zero)
  439. {
  440. static float front_mag_window[WINDOW_SIZE];
  441. static float back_mag_window[WINDOW_SIZE];
  442. static float acc_x_window[WINDOW_SIZE];
  443. static float acc_y_window[WINDOW_SIZE];
  444. static float acc_z_window[WINDOW_SIZE];
  445. static int last_front_zupt;
  446. static int last_back_zupt;
  447. static int front_zupt_wait;
  448. static int back_zupt_wait;
  449. static int front_calm_down_wait;
  450. static int acc_zero_count;
  451. static float back_mag_max_up = 0.0f;
  452. static float back_mag_min_up = 50000.f;//磁力计均方和大小小于32000
  453. static float back_mag_max_down = 0.0f;
  454. static float back_mag_min_down = 50000.f;//磁力计均方和大小小于32000
  455. float front_val = sqrt((float) (front[0] * front[0] + front[1] * front[1] + front[2] * front[2]));
  456. float back_val = sqrt((float) (back[0] * back[0] + back[1] * back[1] + back[2] * back[2]));
  457. //滑动窗口更新数据
  458. memcpy(front_mag_window, front_mag_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  459. memcpy(back_mag_window, back_mag_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  460. memcpy(acc_x_window, acc_x_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  461. memcpy(acc_y_window, acc_y_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  462. memcpy(acc_z_window, acc_z_window + 1, (WINDOW_SIZE - 1) * sizeof(float));
  463. front_mag_window[WINDOW_SIZE - 1] = front_val;
  464. back_mag_window[WINDOW_SIZE - 1] = back_val;
  465. acc_x_window[WINDOW_SIZE - 1] = acc[0] / 2048.f;
  466. acc_y_window[WINDOW_SIZE - 1] = acc[1] / 2048.f;
  467. acc_z_window[WINDOW_SIZE - 1] = acc[2] / 2048.f;
  468. float front_max_val = front_val;
  469. float front_min_val = front_val;
  470. uint16_t front_max_index = WINDOW_SIZE - 1;
  471. uint16_t front_min_index = WINDOW_SIZE - 1;
  472. float back_max_val = back_val;
  473. float back_min_val = back_val;
  474. uint16_t back_max_index = WINDOW_SIZE - 1;
  475. uint16_t back_min_index = WINDOW_SIZE - 1;
  476. for(int i = 0 ; i < WINDOW_SIZE; i++ )
  477. {
  478. if(front_mag_window[i] > front_max_val)
  479. {
  480. front_max_val = front_mag_window[i];
  481. front_max_index = i;
  482. }
  483. if(front_mag_window[i] < front_min_val)
  484. {
  485. front_min_val = front_mag_window[i];
  486. front_min_index = i;
  487. }
  488. if(back_mag_window[i] > back_max_val)
  489. {
  490. back_max_val = back_mag_window[i];
  491. back_max_index = i;
  492. }
  493. if(back_mag_window[i] < back_min_val)
  494. {
  495. back_min_val = back_mag_window[i];
  496. back_min_index = i;
  497. }
  498. }
  499. /*
  500. * 计算稳定的状态
  501. */
  502. float acc_max_val_x = acc_x_window[WINDOW_SIZE - 1] ;
  503. float acc_min_val_x = acc_x_window[WINDOW_SIZE - 1] ;
  504. float acc_max_val_y = acc_y_window[WINDOW_SIZE - 1] ;
  505. float acc_min_val_y = acc_y_window[WINDOW_SIZE - 1] ;
  506. float acc_max_val_z = acc_z_window[WINDOW_SIZE - 1];
  507. float acc_min_val_z = acc_z_window[WINDOW_SIZE - 1];
  508. for(int i = WINDOW_SIZE - 10 ; i < WINDOW_SIZE; i++ )
  509. {
  510. if(acc_x_window[i] < acc_min_val_x)
  511. {
  512. acc_min_val_x = acc_x_window[i];
  513. }
  514. if(acc_x_window[i] > acc_max_val_x)
  515. {
  516. acc_max_val_x = acc_x_window[i];
  517. }
  518. if(acc_y_window[i] < acc_min_val_y)
  519. {
  520. acc_min_val_y = acc_y_window[i];
  521. }
  522. if(acc_y_window[i] > acc_max_val_y)
  523. {
  524. acc_max_val_y = acc_y_window[i];
  525. }
  526. if(acc_z_window[i] < acc_min_val_z)
  527. {
  528. acc_min_val_z = acc_z_window[i];
  529. }
  530. if(acc_z_window[i] > acc_max_val_z)
  531. {
  532. acc_max_val_z = acc_z_window[i];
  533. }
  534. }
  535. if(front_max_index > front_min_index && front_max_val > front_min_val + MAG_THRESHHOLD)
  536. {
  537. front_zero_tmp = 1;
  538. }
  539. else if(front_max_index < front_min_index && front_max_val > front_min_val + MAG_THRESHHOLD)
  540. {
  541. front_zero_tmp = 0;
  542. }
  543. else if(front_zero_tmp == 1)
  544. {
  545. front_zero_tmp = 2;
  546. }
  547. //当触发了压力上升及平稳的时候,需要用加速度倒计时延续状态
  548. if(front_zero_tmp > 0)
  549. {
  550. front_zupt_wait = 20;
  551. }
  552. /*
  553. * 判断后脚往下压
  554. */
  555. if(back_max_index > back_min_index && back_max_val > back_min_val + MAG_THRESHHOLD)
  556. {
  557. back_zero_tmp = 1;
  558. }
  559. else if(back_max_index < back_min_index && back_max_val > back_min_val + MAG_THRESHHOLD)
  560. {
  561. back_zero_tmp = 0;
  562. }
  563. else if(back_zero_tmp == 1)
  564. {
  565. back_zero_tmp = 2;
  566. }
  567. //当触发了压力上升及平稳的时候,需要用加速度倒计时延续状态
  568. if(back_zero_tmp > 0)
  569. {
  570. back_zupt_wait = 20;
  571. }
  572. /*过滤一下类似于穿拖鞋,后鞋垫与磁力计传感器之间的距离在振荡*/
  573. if( back_zero_tmp == 2 && front_max_index < front_min_index && front_max_val > front_min_val + MAG_THRESHHOLD)
  574. {
  575. back_zero_tmp = 0;
  576. }
  577. if(front_zero_tmp == 2 && back_max_index < back_min_index && back_max_val > back_min_val + MAG_THRESHHOLD)
  578. {
  579. front_zero_tmp = 0;
  580. }
  581. /*
  582. * 规避先脚掌下地,后脚掌在下地;接着是后脚先离地,再前脚离地的情况
  583. */
  584. setLongTimeUpTrend(back_max_index, back_min_index, back_val, back_min_val, &back_mag_max_up, &back_mag_min_up);
  585. setLongTimeDownTrend(back_max_index, back_min_index, back_val, back_min_val, &back_mag_max_down, &back_mag_min_down);
  586. if(back_zero_tmp == 2 && back_mag_max_down > back_mag_min_down + 3000.f)
  587. {
  588. back_zero_tmp = 0;
  589. back_mag_max_down = 0.0f;
  590. back_mag_min_down = 50000.0f;
  591. }
  592. if((front_zupt_wait > 0 || back_zupt_wait > 0) && acc_max_val_x - acc_min_val_x < 0.05f && acc_max_val_y - acc_min_val_y < 0.05f && acc_max_val_z - acc_min_val_z < 0.05f)
  593. {
  594. *back_zero = 1;
  595. *front_zero = 1;
  596. }
  597. /*
  598. * 判断仅前踮 或 后垫, IMU往前移动,就不需要这个了
  599. */
  600. // if(back_zero_tmp == 2 && front_zero_tmp == 0)
  601. // {
  602. // back_zero_tmp = 0;
  603. // }
  604. //
  605. // if(front_zero_tmp == 2 && back_zero_tmp == 0)
  606. // {
  607. // front_zero_tmp = 0;
  608. // }
  609. //延迟等待稳定状态
  610. // if(*front_zero == 0 && front_zupt_wait > 0 && acc_max_val - acc_min_val < 0.08f)
  611. // {
  612. // *front_zero = 1;
  613. // }
  614. if(front_zero_tmp > 0)
  615. {
  616. /*
  617. * 如IMU放在前脚掌,就不用front_calm_down_wait == 0的判断了
  618. */
  619. *front_zero = 1;
  620. }
  621. else
  622. {
  623. *front_zero = 0;
  624. }
  625. // if(*back_zero == 0 && back_zupt_wait > 0 && acc_max_val - acc_min_val < 0.08f)
  626. // {
  627. // *back_zero = 1;
  628. // }
  629. if(back_zero_tmp > 0)
  630. {
  631. *back_zero = 1;
  632. }
  633. else
  634. {
  635. *back_zero = 0;
  636. }
  637. for(int i = 0 ; i < WINDOW_SIZE - 5; i++ )
  638. {
  639. if(acc_x_window[i] < acc_min_val_x)
  640. {
  641. acc_min_val_x = acc_x_window[i];
  642. }
  643. if(acc_x_window[i] > acc_max_val_x)
  644. {
  645. acc_max_val_x = acc_x_window[i];
  646. }
  647. if(acc_y_window[i] < acc_min_val_y)
  648. {
  649. acc_min_val_y = acc_y_window[i];
  650. }
  651. if(acc_y_window[i] > acc_max_val_y)
  652. {
  653. acc_max_val_y = acc_y_window[i];
  654. }
  655. if(acc_z_window[i] < acc_min_val_z)
  656. {
  657. acc_min_val_z = acc_z_window[i];
  658. }
  659. if(acc_z_window[i] > acc_max_val_z)
  660. {
  661. acc_max_val_z = acc_z_window[i];
  662. }
  663. }
  664. if(acc_max_val_x - acc_min_val_x < 0.03f && acc_max_val_y - acc_min_val_y < 0.03f && acc_max_val_z - acc_min_val_z < 0.03f )
  665. {
  666. *acc_zero = 1;
  667. acc_zero_count ++;
  668. }
  669. else
  670. {
  671. *acc_zero = 0;
  672. acc_zero_count = 0;
  673. }
  674. // /*
  675. // * 强制大于1秒的零速均视为触地
  676. // */
  677. if(acc_zero_count > 100)
  678. {
  679. *front_zero = 1;
  680. *back_zero = 1;
  681. front_zero_tmp = 1;
  682. back_zero_tmp = 1;
  683. }
  684. //利用加速度延续
  685. if((last_front_zupt == 1|| last_back_zupt == 1)
  686. && fabsf(acc_x_window[WINDOW_SIZE - 2] - acc_x_window[WINDOW_SIZE - 1]) < 0.05f
  687. && fabsf(acc_y_window[WINDOW_SIZE - 2] - acc_y_window[WINDOW_SIZE - 1]) < 0.05f
  688. && fabsf(acc_z_window[WINDOW_SIZE - 2] - acc_z_window[WINDOW_SIZE - 1]) < 0.05f)
  689. {
  690. *front_zero = 1;
  691. *back_zero = 1;
  692. }
  693. if(front_zupt_wait > 0)
  694. {
  695. front_zupt_wait --;
  696. }
  697. if(back_zupt_wait > 0)
  698. {
  699. back_zupt_wait --;
  700. }
  701. if(*front_zero || *back_zero)
  702. {
  703. front_zupt_wait = 20;
  704. back_zupt_wait = 20;
  705. }
  706. last_front_zupt = *front_zero ;
  707. last_back_zupt = *back_zero;
  708. /*
  709. * 直接用这个来记步好了
  710. */
  711. if(!IMU_GetGameMode())
  712. cal_step(*front_zero, *back_zero, acc_x_window[WINDOW_SIZE - 1], acc_y_window[WINDOW_SIZE - 1], acc_z_window[WINDOW_SIZE - 1]);
  713. }