kramann.info
© Guido Kramann

Auswertung der Meßdaten zur Sprungantwort mit Scilab

(EN google-translate)

(PL google-translate)

Es wurden mit Hilfe von "BLUSCI7SPRUNG" (vergl. vorangehendes Unterkapitel) folgende Rohdaten gewonnen:

arr=[0, -124; 
68, -125; 
142, -117; 
201, -42; 
267, 48; 
333, 140; 
400, 120; 
467, 130; 
539, -27; 
572, 42];

Code 0-1: arrsprung4.sci - Rohdaten zur Sprungantwort.

Um die Daten leicht nach Scilab importieren zu können, schreibt das BlueJ-Projekt ein Scilab-Array heraus.
Die erste Spalte ist die Zeit in Millisekunden.
Die zweite Spalte ist die Position des grünen Kreises relativ zur Spiegel-/ Kamersabildmitte in Pixeln.
Um den Sprung über einen längeren Zeitraum aufnehmen zu können, wurde nicht bei Spiegelposition Null, sondern links davon begonnen.

1. Nach 400 Millisekunden sind die Daten unbrauchbar (Umgebungslicht dringt auf die Kamera):

arr=[0, -124; 
68, -125; 
142, -117; 
201, -42; 
267, 48; 
333, 140];

Code 0-2: Daten reduziert.

2. Der Sprungverlauf sollte bei Null losgehen, deshalb Offset von +124 auf die zweite Spalte addieren:

3. Der Rückgang um 1 in der zweiten Zeile macht keinen Sinn und wird auf Null gesetzt.

arr=[0, -124; 
68, -125; 
142, -117; 
201, -42; 
267, 48; 
333, 140; 
400, 120; 
467, 130; 
539, -27; 
572, 42]; 

A = arr(:,1);
B = arr(:,2);
B=B+124;
B(2)=0;
C = [A(1:6,1),B(1:6,1)]
plot(C(:,1),C(:,2))

Code 0-3: Die bisheringen Schritte als Scilab-Befehle umgesetzt.

Bild 0-1: Sprungantwort (Winkel) als Plot.

4. Nicht die Übertragung von PWM auf Winkel, sondern von PWM auf Winkelgeschwindigkeit wird identifiziert:

Dies erfordert eine numerische Ableitung des Winkels:

Q=diff(C(:,2))    //Differenzbildung phi(i+1)-phi(i), man erhält einen Wert weniger!
R = diff(C(:,1))  //wie Q, jedoch um die Zeitdifferenzen zu bilden.
RR=Q./R           //Bilden von delta phi(i) / delta t(i)
CC=C(2:6,1)       //Zeitachse neu bilden
CC=CC-68          //Offset abziehen: Erster Wert t0 soll Null sein.
plot(CC,RR)       //Testplot.

Code 0-4: Numerische zeitliche Ableitung des Winkels

Bild 0-2: Sprungantwort (Winkelgeschwindigkeit) als Plot.

Grundlage für die nachfolgende Systemidentifikation soll eine Darstellung in Sekunden sein:
Zeit in Sekunden (CA)
Winkelgeschwindigkeit in Pixel pro Sekunden (CB)

CA=CC(:,1)
CB=RR(:,1)
CA=CA.*0.001
CB=CB.*1000
ZZ=[CA,CB]

Code 0-5: Umrechnung in Sekunden.

ZZ  = [0.    ,   0.;         
       0.074 ,   94.594595;  
       0.133 ,   1271.1864;  
       0.199 ,   1363.6364;  
       0.265 ,   1393.9394];

Code 0-6: Grundlage für die nachfolgende Systemidentifikation.