Berechnungsproblem mit VBA

Hallo, zusammen!
Ich habe ein Problem mit VBA, kann jemand hier vielleicht mir helfen?
Um eindeutig zu sein, nehme ich hier ein Beispiel mit nur 10 Messdaten darin (tatsächlich gibt es 2000 Messdaten).
Meine Aufgabe ist Kxx zu berechnen, der der Durchschnitt von Summe von X(ti)X(ti+τ) ist.
Abtastzeit ti...........Messdaten X(ti)
t1.................................X(t1)
t2.................................X(t2)
t3.................................X(t3)
t4.................................X(t4)
t5.................................X(t5)
t6.................................X(t6)
t7.................................X(t7)
t8.................................X(t8)
t9.................................X(t9)
t10................................X(t10)
So läuft meine Rechnung:
1.τ = 0
Kxx = [X(t1)*X(t1) + X(t2)*X(t2) + X(t3)*X(t3) ... ... + X(t10)*X(t10)]
2.τ = 1
Kxx = [X(t1)*X(t2) + X(t2)*X(t3) + X(t3)*X(t4) + X(t4)*X(t5) + X(t5)*X(t6) + X(t6)*X(t7) + X(t7)*X(t8) + X(t8)*X(t9) + X(t9)*X(t10)]
3.τ = 2
Kxx = [X(t1)*X(t3) + X(t2)*X(t4) + X(t3)*X(t5) + X(t4)*X(t6) + X(t5)*X(t7) + X(t6)*X(t8) + X(t7)*X(t9) + X(t8)*X(t10)]
4.τ = 3
Kxx = [X(t1)*X(t4) + X(t2)*X(t5) + X(t3)*X(t6) + X(t4)*X(t7) + X(t5)*X(t8) + X(t6)*X(t9) + X(t7)*X(t10)]
5.τ = 4
Kxx = [X(t1)*X(t5) + X(t2)*X(t6) + X(t3)*X(t7) + X(t4)*X(t8) + X(t5)*X(t9) + X(t6)*X(t10) ]
Der max. Wert von τ ist halb von 10 gleich 5.
Da mein Kenntnis in VBA ganz wenig ist, kann ich dieses Problem nicht alleine schaffen, ich brauche eure Hilf, vielen Dank!

Beiträge zum Forumthread

Beiträge zu diesem Forumthread

Betreff

Datum

Anwender

Hallo, @caoli,
es war wohl etwas zu spät bei mir gewesen, als ichj meinen Beitrag verfasste, denn der die Ergebnisse des Summenproduktansatzes liefern nur unter bestimmten Vorberechnungen der Daten (gemittelte Werte der Diferenzen zum Mittelwert der Daten) die Autokorrelationsfunktion, mein entsprechnder Hinweis hatte also in der Antwort nichts zu suchen.
Auch die Indexvariblen im Ansatz lassen sich "überscihtlicher" ensprechnden Deiner Benutzung einsetzen, obwohl die vorgestellte Indexberechnung gleichwertig ist.
Beispieldatei zur Berechnung von Kxx

Option Explicit
Const WSDATEN = "Tabelle1"   'Tabellenblattname der Datenreihe UND Ausgabe
Const DATACOLUMN = 2         'Spaltennummer der Datenreihe
Const DATAHEADERROWS = 1     'Anzahl der Kopfzeilen vor Datenbeginn
Const RESULTCOLUMN = 4       'Spaltennummer der Ausgabespalten (drei aufeinanderfolgende Spalten)
Sub CalcKxx()
Dim ws As Worksheet
Dim x() As Double, acf() As Double, maxI As Long, i As Long, k As Long
Set ws = Worksheets(WSDATEN)
'Ermitteln der Datenanzahl der Datenreihe
maxI = ws.Cells(ws.Rows.Count, DATACOLUMN).End(xlUp).Row - DATAHEADERROWS
ReDim x(maxI)
ReDim Kxx(maxI / 2 + 1) 'Es werden die Ergebnisse in einem Datenfeld gespeichert
With ws
'Einlesen der Daten
For i = 1 To maxI
x(i) = .Cells(DATAHEADERROWS + i, DATACOLUMN)
Next i
'Inhalte der Ausgabe-/Ergebnisspalten löschen
.Columns(RESULTCOLUMN).Clear
.Columns(RESULTCOLUMN + 1).Clear
.Columns(RESULTCOLUMN + 2).Clear
'Überschrift über Ausgabe-/Ergebnisspalten
.Cells(DATAHEADERROWS, RESULTCOLUMN) = "lag k"
.Cells(DATAHEADERROWS, RESULTCOLUMN + 1) = "Kxx(i)"
.Cells(DATAHEADERROWS, RESULTCOLUMN + 2) = "Kxx(i)/Kxx(0)"
For k = 0 To maxI / 2 - 1
For i = 1 To maxI - k
Kxx(k) = Kxx(k) + x(i) * x(i + k)
Next i
.Cells(k + 1 + DATAHEADERROWS, RESULTCOLUMN) = k
.Cells(k + 1 + DATAHEADERROWS, RESULTCOLUMN + 1) = Kxx(k)
.Cells(k + 1 + DATAHEADERROWS, RESULTCOLUMN + 2) = Kxx(k) / Kxx(0) 'Programmzeile muß verändert werden, wenn Kxx nicht in einem Datenfeld gespeichrt wird
Next k
End With
Set ws = Nothing
End Sub

Auf die Abspeicherung der Ergebnisse für Kxx in den einzelnen Scheifendurchgängen der Laufvariablien k, kann versichtet werden, wenn sie für die weitere Auswertung nicht benötigt werden. Damit könnte dann das Datenfeöd Kxx() entfallen
Gruß,
Uwe

For k = 0 To maxI / 2 - 1 ' Setze Speicherstelle der Variablen dKxx auf Null, ' zur Berechnung mit dem Verschub k dKxx = 0 For i = 1 To maxI - k 'rechte Seite der Zuweisung von rechts nach links betrachtet: ' - Bilde Teilprodukt x(i)*x(i+k) und ' - addiere das Ergebnis zu der bisher ' in dieser Schleife (i=1 To maxI-k) ermittelten Summe ' (Inhalt der Speicherstelle dKxx) 'linke Seite der Zuweisung: ' - Schreibe das neue Ergebnis an die Speicherstelle dKxx dKxx = dKxx + x(i) * x(i + k) 'wenn (i<= maxI-k), führe dei Anweisung 'der FOR-Schleife für das nächste i aus Next i ' Hier wird nun einmal die Division zur Mittelwertbildung ' für die Verschiebung k ausgeführt dKxx = dKxx / (maxI-k) ' hier wird zentral der Faktor (1/n) berücksichtigt 'Hier wird das Ergebnis geschrieben .Cells(k + 1 + DATAHEADERROWS, RESULTCOLUMN) = k .Cells(k + 1 + DATAHEADERROWS, RESULTCOLUMN + 1) = dKxx 'Optinal kann aber auch das Ergebnis in einem Feld zwischen gespeichert werden Kxx(k)=dKxx Next k

Beliebteste Forumthreads (12 Monate)

	A	B	C	D	E	F	G	H
1	ti	Xi		r	Kxx	DKxx	Kxx	DKxx
2		1,1		0	2551,89	141,7716667	2551,89	141,7716667
3		2,2		1	2344,98	137,94	2344,98	137,94
4		3,3		2	2139,28	133,705	2139,28	133,705
5		4,4		3	1936	129,0666667	1936	129,0666667
6		5,5		4	1736,35	124,025	1736,35	124,025
7		6,6		5	1541,54	118,58	1541,54	118,58
8		7,7		6	1352,78	112,7316667	1352,78	112,7316667
9		8,8		7	1171,28	106,48	1171,28	106,48
10		9,9		8	998,25	99,825	998,25	99,825
11		11		9	834,9	92,76666667	834,9	92,76666667
12		12,1
13		13,2
14		14,3
15		15,4
16		16,5
17		17,6
18		18,7
19		19,8

8 Beiträge zum Forumthread Beiträge zu diesem Forumthread

Beliebteste Forumthreads (12 Monate)

Beliebteste Forumthreads (12 Monate)

8
Beiträge zum Forumthread

Beiträge zu diesem Forumthread