[DywanTv]

Siemka. Piszę apke w #matlab mam. Mam takie dane, które reprezentują przemieszczenie. Chciałbym policzyć średnią prędkość w kolejnych powtórzeniach ruchu, czyli potrzebuję wykryć początek wzniesienia, pik no i najlpeiej koniec kiedy opada. Macie jakieś pomysły jak mógłbym to obliczać?
Do głowy przychodzi mi algorytm Pan-Tompkinsa, ale czy nie jest to overkill?
#pytanie #programowanie #matematyka

[motoinzyniere]

@DywanTv: Pochodna Ci powie kiedy zmienia kierunek, druga pochodna kiedy się przegina. Maxima znajdziesz przy pomocy findpeaks.

[DywanTv]

@motoinzyniere: Pierwsza pochodna chyba powie o tym jak szybko wzrasta albo maleje, a druga o tym kiedy zmienia kierunek - czyli się przegina - nie tak będzie?

[jalop]

@motoinzyniere: @DywanTv: możesz jeszcze lokalnie spróbować przez islocalmax

[DywanTv]

@motoinzyniere: @jalop: Właściwie nawet samo findpeaks(y) i findpeaks(-y) z odpowiednimi parametrami MinPeakProminence działa całkiem dobrze. Niebieski sygnał to -y, ale punkty uzyskuję całkiem dokładnie. Dzięki za pomoc :)

[motoinzyniere]

@DywanTv: No...o to mi chodziło.

[heheszkant]

@DywanTv: Dane z akcelerometru telefonu? Jak załatwiłeś orientację osi pionowej względem osi akcelerometru, przez zwracanie uwagi na odpowiednie ustawienie urządzenia względem pionu czy obróciłeś se układ współrzędnych na podstawie przyspieszenia ziemskiego już w matlabie?
Przysiad czy CGBP kogoś z długimi rękami czy sumo?

[DywanTv]

@heheszkant: Dane na podstawie śledzenia toru ruchu sztangi na nagranym wcześniej filmiku. Ciągi klasyczne Clarence'a Kennedy'ego xD https://youtu.be/QdTc2Gv52uI?t=102

Wykorzystuję to, że nagrane filmiki są zazwyczaj idealnie poziomo ze względu na budowę telefonów czy aparatów, które mają płaską obudowę - ustawia się je też na płaskim boksie czy podłodze, więc trafiłeś:

zwracanie uwagi na odpowiednie ustawienie urządzenia względem pionu

Większa dokładność nie jest mi potrzebna

[DywanTv]

@heheszkant: Tamten wykres wyżej nie jest w odpowiedniej skali. Tutaj masz wykres przemieszczenia [cm] w zależności od czasu [s]. Jak widać dokładność takiego śledzenia nie jest bardzo duża, ale to raczej projekt do nauki programowania, a nie jakiś komercyjny projekt. Fajnie, że widać, że 3 i 4 powtórzenia są wyraźnie wolniejsze, że zbocze wstępujące jest mniej strome, widać nawet początkowe szarpnięcie sztangi. Może później przepiszę to na Pythona w ramach nauki,

[heheszkant]

@DywanTv: Zakres ruchu mi nie pasował do zwykłego ciągu, ale nie odjąłem połowy wysokości talerza xD
Rozumiem, że śledzisz ruch przyśrodkowego krańca lewego rękawa sztangi i skalujesz wymiary według najbliższego talerza?
Perspektywa trochę średnia pod względem zniekształceń. Wiesz, skąd się wziął drift zera w danych?
Dorób sobie jeszcze drugie badane kryterium, czyli minimalne generowane przyspieszenie w trakcie fazy koncentrycznej i na jakiej wysokości przypada, bo tam jest generowana najmniejsza siła.

[DywanTv]

@heheszkant: Tak dokładnie tak. Chociaż tu użytkownik określa, który punkt śledzić. Po prostu środek talerza niemal zawsze wyraźnie się odznacza, dzięki czemu jest prosty do śledzenia. Wymiary mam na podstawie najbliższego talerza. Drift zera prawdopodobnie wziął się stąd, że obliczam ruch w osi X i Y, a nie tylko Y, przez co ruch "do przodu" i "do tyłu" wpływa na wyniki. Wystarczy, że sztanga skończy w nieco innym miejscu np. przed