Signalverarbeitung die zweite

Hallöchen alle zusammen,

kurz vor Weihnachten hat mich "Kvasir" auf einen besseren Programmierstil aufmerksam gemacht.
Jetzt bin ich dabei diesen Stil für ein bestehendes Programm anzuwenden.

In meinem alten Programm suche ich im ersten Teil nach Punkten die meine x-Achse schneiden und zwar immer von negativ -> positiv und danach von positiv -> negativ. Im weiteren Text sage ich dazu immer Wertepaar, weil ich glaube dass dies einfacher ist.

Problem 1:

Mein Programm registriert natürlich nicht nur die Wertepaare die ich haben will, sondern auch hier und da kleine Peaks die ca. 12ms lang sind. Diese würde ich gern herausschneiden das zum Schluss nur noch meine richtigen Wertepaare im Arry stehen.
Die Lösung für dieses Problem ist einfach ich nehme einfach den letzten Wert meine Wertepaares und ziehe den ersten davon ab. Die Differenz sofern kleiner als z.B.: 12ms muss ein falsches Wertepaar sein und wird entfernt.

Frage:

Geht es bei diesem Programmierstil das ich diese Bedingung gleich mit in meine Wertepaarsuche integriere oder ist es besser dass ich erst alle Wertepaare aufnehme und hinterher untersuche.

Problem 2:

Ich habe festgestellt dass der 2te Punkt meines ersten Wertepaares nicht gefunden wird. Er steht in der Zeile 20 meiner Beispieldatei. Der erste liegt in der Zeile 16. (Edit: Sorry ab der Zeile 1900 tritt ein ähnliches Problem auf.)

Frage:

Was passiert dort? Das alte Programm findet diesen Punkt. Erst dachte ich der Abstand ist zu klein aber dann würde das NulldurchgangsVi
auch die andern zu kleinen Wertepaare nicht oder nur teilweise finden was nicht der Fall ist. Oder liegt es vielleicht daran dass nur ein Punkt kurzzeitig negativ ist bevor der Graphverlauf wieder steigt.

Ihr findet hier das Programm und die Datendatei. Vielleicht solltet ihr euch die letzte Spalte der Datendatei einmal graphisch anzeigen lasse.
Diese Werte sind übrigens die Differenz der Werte zweier Drucksensoren. Es wurde deshalb gemacht damit man nicht danach suchen muss wo die Schnittpunkte der beiden Graphen liegen.

Gruß Gerrit

Anzeige

Zu 1:
Für die Diskriminierung ungültiger Nulldurchgänge gibt es verschiedne Strategien, aber ganz frei von Willkür ist das immer nicht. Man kann immer Beispiele finden, in denen es nicht so wunschgemäß funktioniert, z.B. wenn mehrere Peaks dicht hintereinander folgen.
Ich würde die Entfernung auf jeden Fall "Online" vornehmen, und zwar so: Ein Übergang, z.B. -/+, ist nur dann gültig, wenn z.B die nachfolgenden 10 Samples alle positiv bleiben. Beispiel:
---++-++++--+++++++++++------ (nur der rote Wert wird gezählt)
Das läßt sich mit einem Zähler machen, wenn Du damit nicht zurechtkommst, helfe ich Dir gern weiter.

Zu 2:
Kann das nicht nachvollziehen, bei mir kommt genau das raus was ich erwarte. Ich hasse allerdings Kommas als Dezimalzeichen in meinem System, und um Deine Komma-Daten lesen zu können, mußte ich den Formatstring ändern.
   

Habe Dein VI etwas abgespeckt. Das Nulldurchgang-Vi ist wird hier nicht gebraucht. Wenn Du z.B in den State -+ eintrittst, sind die Werte ja immer erst mal negativ, und Du suchst den ersten positiven Wert. Da genügt ein einfache Abfrage >0.

Hallo Lucki

ich finde den ersten Weg auch besser. Die Idee kam mir auch schon nur die Umsetzung ist eine andere Sache. Bin eben kein Programmierer. (Ein Einfaches subtrahieren bekomme aber ich mittlerweile hin.)
Und die LV Vorlesung die dieses Semester an unsere FH angeboten wurde ging nur über 6 x 2h.
Von daher nehme ich dein Angebot mit dem Zähler an.

Gerrit

Habe mal was gemacht. Aber wie gesagt, Aussortierung "falscher Nulldurchgänge" ist nie frei von Willkür. Ich habe es hier so gemacht, das ein Nulldurchgang nur dann zählt, wenn die nachfolgenden N Werte alle das gleiche Vorzeichen haben. Ebenso gut hätte man es so machen können, daß ein Nulldurchgang nur zählt, wenn es in den vorangegangenen Werten keinen Vorzeichenwechsel gab. (das entspricht dem ersten Fall, wenn man das Array von hinten anfangen würde zu untersuchen).
Weitere Möglichkeiten, um "falsche" Nulldurchgänge zu eliminieren:
[list]
[*]Hysterese (Schmitt.Trigger-Verhalten)<>
[*]Totzeit nach jedem gefundenem Durchgang<>
[*]Tiefpassfilterung<>
[st]Habe festgestellt, ab N=5 bleibt die Anzahl Nulldurchgänge sehr schön stabil auf 11

Danke für deine Hilfe Lucki,

jetzt habe ich schon mal einen Ansatz wie ein Zähler programmiert wird.

Ich werde jetzt versuchen dass das Programm als ersten nach einem ND von minus -> plus sucht. Damit dürfte dann auch der erste Wert (19) verschwinden. (Sofern du das noch nicht gemacht hast!)

Nur eine Sache ist merkwürdig! Ich habe mal die Indexe überprüft - wenn das Programm eine ND n->p findet merkt es sich nicht den ND sondern 3 Werte vorher und von p->n 2 Werte vorher. Mal sehen wie das zustande kommt! Ein wenig Arbeit muss mir schließlich auch noch bleiben.

Ich danke nochmals
Gerrit

Zitat:jetzt habe ich schon mal einen Ansatz wie ein Zähler programmiert wird.
Ich werde jetzt versuchen dass das Programm als ersten nach einem ND von minus -> plus sucht. Damit dürfte dann auch der erste Wert (19) verschwinden. (Sofern du das noch nicht gemacht hast!)

Das wußte ich nicht, daß das Ergebnis-Array immer mit dem ersten ND-+ beginnen soll. Habe das geändert.

Zitat:Nur eine Sache ist merkwürdig! Ich habe mal die Indexe überprüft - wenn das Programm eine ND n->p findet merkt es sich nicht den ND sondern 3 Werte vorher und von p->n 2 Werte vorher. Mal sehen wie das zustande kommt! Ein wenig Arbeit muss mir schließlich auch noch bleiben.

Bei mir tritt das nicht auf, siehe Diagramm
   

Hab das VI abgeändert, und zwar so, daß es nun wirklich keiner mehr verstehen kann. Wenigstens diese Erläuterungen:[list=1]
[*]zu "Case1": Das Vi fängt immer mit den richtigen Werten an zu suchen. Wenn der erste Datenwert positiv ist (Wie in Deinen Daten), wird als erstes nach einem ND+- gesucht, und umgekehrt. (Wenn man dem "falschen" Wert beginnt, dann hat man beim Datenpunkt 0 fehlerhafterweise schon einen Durchgang).<>
[*]zu "Case2": Falls aber der erste gefundene Wert, wie in Deinen Daten, kein ND-+ ist, wird dieser nicht verwendet, d.h er kommt nicht in das Ergebnis-Array.<>
[*]Die Schleife "For1": Die Zahl N=17 hat den Hintergrund, daß ich - rein willkürlich - 17 Kursoren erzeugt habe.<>
[st]