Zde můžete vidět rozdíly mezi vybranou verzí a aktuální verzí dané stránky.
| Obě strany předchozí revize Předchozí verze Následující verze | Předchozí verze | ||
| 
                    2018:accel-roll-pitch [2019/01/14 21:14] Ondrej Vörös  | 
                
                    2018:accel-roll-pitch [2019/01/14 23:58] (aktuální) Ondrej Vörös  | 
            ||
|---|---|---|---|
| Řádek 15: | Řádek 15: | ||
| Derotačná matica bude počítaná pomocou Matlabu, pretože obsahuje goniometrické funkcie. | Derotačná matica bude počítaná pomocou Matlabu, pretože obsahuje goniometrické funkcie. | ||
| - | Program bude možno ovládať cez UART konzolu (vložiť derotačnú maticu, kompenzovať výpočet náklonu). | + | Program bude možné ovládať cez UART konzolu (vložiť derotačnú maticu, kompenzovať výpočet náklonu). | 
| {{https://www.novatel.com/assets/Web-Phase-2-2012/Solution-Pages/AttitudePlane.png }} | {{https://www.novatel.com/assets/Web-Phase-2-2012/Solution-Pages/AttitudePlane.png }} | ||
| Řádek 33: | Řádek 33: | ||
| ==== Výpočet náklonu ==== | ==== Výpočet náklonu ==== | ||
| - | Pre výpočet pozdĺžneho a priečneho náklonu roll a pitch bola vytvorená funkcia get%degree_tilt. Pri pomalom otáčaní | + | Pre výpočet pozdĺžneho a priečneho náklonu roll a pitch bola vytvorená funkcia ''get_degree_tilt''. Pri pomalom otáčaní | 
| pôsobí na akcelerometer hlavne gravitačná sila, zrýchlenie v jednotlivých osách akcelerometru je dané gravitačným | pôsobí na akcelerometer hlavne gravitačná sila, zrýchlenie v jednotlivých osách akcelerometru je dané gravitačným | ||
| zrýchlením. Skusmým otáčaním dosky s akcelerometrom boli zistené maximálne hodnoty zrýchlenia v jednotlivých osách, | zrýchlením. Skusmým otáčaním dosky s akcelerometrom boli zistené maximálne hodnoty zrýchlenia v jednotlivých osách, | ||
| Řádek 58: | Řádek 58: | ||
| </code> | </code> | ||
| - | {{2018:accel-roll-pitch:sur_system-na-doske.png?400 }} | + | Na nasledujúcom obrázku je zobrazený pravotočivý súradnicový systém na vývojovej doske. Podľa hodnôt získaných | 
| + | z akcelerometru má ten orientáciu osí podľa ľavotočivého súradnicového systému, preto pri výpočte náklonu sa | ||
| + | hodnotám zrýchlenia v ose x otáča znamienko (otáčame x-ovú os), keďže výpočet derotačnej matice podľa | ||
| + | [[https://www.st.com/resource/en/design_tip/dm00358510.pdf|DT0076]]  aj podľa | ||
| + | [[https://www.mathworks.com/help/phased/ref/rotx.html| funkcií MATLABU]] prebieha v pravotočivom súradnicovom | ||
| + | systéme. | ||
| + | |||
| + | {{ 2018:accel-roll-pitch:sur_system-na-doske.png?400 }} | ||
| + | |||
| + | Priečny náklon je definovaný ako rotácia okolo osi x (použité hodnoty zrýchlenia  v osi y) a pozdĺžny náklon | ||
| + | ako rotácia okolo osi y (použité hodnoty zrýchlenia  v osi x). Náklon v jednotlivých smeroch je indikovaný aj | ||
| + | pomocou štyroch LED diód umiestených  na vývojovej doske do kríža po stranách akceleremetru (funkcia ''tilt_light''). | ||
| + | |||
| + | <code c> | ||
| + | while (1) | ||
| + | { | ||
| + | if ((HAL_GetTick() - tick) >= per_meas) { // period of measurement | ||
| + | BSP_ACCELERO_GetXYZ(buffer); | ||
| + | accX = buffer[0]; | ||
| + | accY = buffer[1]; | ||
| + | accZ = buffer[2]; | ||
| + | |||
| + | tilt_light(accX, accY); | ||
| + | |||
| + | pitch = get_degree_tilt(-accX, max_xval); //pitch- acceleration around axis y --> accX , -accX --turning around axis X to get right-rotated coordinate system | ||
| + | roll = get_degree_tilt(accY, max_yval); // get roll- acceleration around axis x --> accY | ||
| + | |||
| + | printf("roll,pitch:\t%2.1f%c\t%2.1f%c\n", roll, deg_sign, pitch, deg_sign); | ||
| + | |||
| + | tick = HAL_GetTick(); | ||
| + | } | ||
| + | ... | ||
| + | } | ||
| + | </code> | ||
| ==== Výpočet derotačnej matice a kompenzácia náklonu ==== | ==== Výpočet derotačnej matice a kompenzácia náklonu ==== | ||
| + | Derotačná matica je inverzná alebo transponovaná rotačná matica. Rotačné matice pre rotáciu okolo jednotlivých osí | ||
| + | sú [[https://www.mathworks.com/help/phased/ref/rotx.html|definované]] takto: | ||
| + | Rx = [1 0 0; 0 cos(roll) -sin(roll); 0 sin(roll) cos(roll)];  % rotacna matica okolo osi x- counterclockwise | ||
| + | Ry = [cos(pitch) 0 sin(pitch); 0 1 0; -sin(pitch) 0 cos(pitch)]; % rotacna matica okolo osi y- counterclockwise | ||
| + | Rz = [cos(yaw) -sin(yaw) 0; sin(yaw) cos(yaw) 0; 0 0 1]; % rotacna matica okolo osi z- counterclockwise | ||
| - | ---- | + | Keďže pomocou akcelerometru nie je možné počítať uhol kurzu (otočenie okolo osi z- yaw), rotačnú maticu ''Rz'' | 
| - | ====== Záver ====== | + | nevyužijeme. Ako derotačnú maticu sme použili ''M = Rx*Ry'' (bola určená experimentálne, aby vektor zrýchlenia | 
| + | [1000 0 0] pri rotácii o 90° pozdĺžneho náklonu-pitch sa otočil do polohy [0 0 1000]- nasmerovanie osi z | ||
| + | do smeru gravitačného zrýchlenia). Rovnako, pri použití rotačnej matice M sa vektor zrýchlenia [0 -1000 0] | ||
| + | pri rotácii o 90° priečneho náklonu-roll otočí do polohy [0 0 1000]. | ||
| - | [[https://drive.google.com/open?id=1P0V9NppeeiWWQ49n4yWtXbpKTopOQkus|zdrojové súbory]] | + | Derotovanie dátového vektora (v našom prípade vektoru zrýchlenia) prebieha vynásobením vektoru s derotačnou maticou. | 
| + | Následne možno z derotovaného vektoru znovu vypočítať uhly náklonu. | ||
| + | accXr = M(1,1)*accX + M(1,2)*accY + M(1,3)*accZ | ||
| + | accYr = M(2,1)*accX + M(2,2)*accY + M(2,3)*accZ | ||
| + | accZr = M(3,1)*accX + M(3,2)*accY + M(3,3)*accZ | ||
| + | |||
| + | ==== Štruktúra programu ==== | ||
| + | Štruktúra programu je nasledovná, v nekonečnej smyčke prebieha výpočet a zobrazovanie náklonu v časovom intervale | ||
| + | jedna sekunda a tiež spracovanie znakov z kruhového zásobníka prijímaných sériovým rozhraním UART. K ovládaniu | ||
| + | programu cez sériové rozhranie slúžia dva príkazy:  | ||
| + | * ''COMPENSATE'' - spustenie kompenzovaného výpočtu náklonu (30 meraní s periódou 1 sekunda) | ||
| + | * ''MATRIX m11 m12 m13 m21 m22 m23 m31 m32 m33'' - zadanie derotačnej matice (zadanie prvkov matice po riadkoch) | ||
| + | |||
| + | ----- | ||
| + | |||
| + | ====== Záver ====== | ||
| + | |||
| + | Výsledkom riešenia projektu je meranie pozdĺžneho a priečneho náklonu. Tiež bola implementovaná derotácia dátového | ||
| + | vektoru zrýchlenia pre ľubovoľnú polohu akcelerometru voči referenčnej rovine. Kompenzácia dát derotáciou | ||
| + | dátového vektoru zrýchlenia podľa hodnôt náklonu akcelerometru vo zvolenej polohe nie je žiaľ presná (nepresnosť | ||
| + | kolíše v závislosti na zvolenej polohe akcelerometru) a nepodarilo sa mi prísť (kvôli nedostatku času) na to, | ||
| + | kde môže byť chyba.  | ||
| + | [[https://drive.google.com/open?id=16A5mXGF80x7sn3EjSTIg6NrUqWS896O5|zdrojové súbory]] | ||
| [[https://drive.google.com/open?id=17qJpzuhK15dT8O1v8BeUGPA5j-lbc9JG|videoukážka]] | [[https://drive.google.com/open?id=17qJpzuhK15dT8O1v8BeUGPA5j-lbc9JG|videoukážka]] | ||
| ---- | ---- | ||
| + | |||
| ====== Zdroje ====== | ====== Zdroje ====== | ||
| - | DT0076- Compensating for accelerometer installation error: zeroing pitch and roll for a reference orientation: | + | https://www.st.com/resource/en/design_tip/dm00358510.pdf\\ | 
| - | [[https://www.st.com/resource/en/design_tip/dm00358510.pdf]] | + | https://www.zive.cz/clanky/pojdme-programovat-elektroniku-jak-vlastne-funguje-akcelerometr-a-gyroskop-nejen-ve-vasem-telefonu/sc-3-a-194858/default.aspx#articleStart\\ | 
| - | https://www.zive.cz/clanky/pojdme-programovat-elektroniku-jak-vlastne-funguje-akcelerometr-a-gyroskop-nejen-ve-vasem-telefonu/sc-3-a-194858/default.aspx#articleStart | + | |
| https://www.mathworks.com/help/phased/ref/rotx.html | https://www.mathworks.com/help/phased/ref/rotx.html | ||