Die Taktrückgewinnung wurde hier durch clockrec realisiert, die Funktion liefert die Taktverschiebung der Aufnahme. Mit Taktverschiebung ist hier der Abstand bis zum Anfang des ersten vollen Takts der Aufnahme gemeint, hier als dt bezeichnet. Ausgehend von dieser kann man den übrigen Teil der Aufnahme in Abschnitte mit der Länge eines Takts einteilen. Dabei bleibt allerdings auch am Ende der Aufnahme ein Rest übrig, welcher nicht verarbeitet werden kann und entweder weg fällt oder im darauf folgenden Aufnahmestück verarbeitet wird.

Wichtig ist, dass der Funktion mitgeteilt wird mit welcher Taktlänge und auf welcher Taktfrequenz gesendet wird. Das wird hier für eine bessere Übersichtlichkeit mit lokalen Variablen festgelegt, kann sonst aber auch über globale Variablen geschehen. Nun wird über einen digitalen Butterworth-Filter das Signal mit der Taktfrequenz aus der Aufnahme herausgefiltert. Dabei sind a und b Variablen, die von der Filter-Funktion benötigt werden, welche dann das gefilterte Signal yf zurück liefert. Jetzt wird die Hüllkurve um den positiven Bereich des gefilterten Signals gebildet. Das geschieht mithilfe der Hilbert-Transformation. Das so erhaltene Signal yh ist ein binäres Signal, welches die zwei Zustände 'An' und 'Aus' annehmen kann und dem gesendeten Taktsignal entspricht. Daraus wollen wir jetzt die Zeitpunkte der Taktgrenzen ermitteln. Dazu muss eine Schranke s festgelegt werden, die darstellt ab wann ein Signal als 'An' gilt. Diese berechnet sich hier als ein Fünftel des Maximalwerts der Hüllkurve. Nun werden die Indizes ermittelt, bei welchen diese Schranke (mit einer Toleranz tol) durchschritten wird. Das sind die Stellen wo ein Takt endet und der nächste anfängt. Wenn man diese Werte in idx ganzzahlig durch die Taktlänge in Samples teilt erhält man die Taktverschiebung. Aus diesen vielen Werten in dt wird jetzt noch der Mittelwert gebildet und man erhält die gesuchte Taktverschiebung, die für die Decodierung benötigt wird.

⇒ ⇒ ⇒

def clockrec(y):
	'''liefert die Verschiebung des Taktes'''
	tl = 0.05
	f_takt = 12600
 
	b, a = signal.butter(4, [(f_takt-100)*2./RATE, (f_takt+100)*2./RATE], btype="bandpass") # Parameter für Butterworth-Filter
	yf = signal.lfilter(b,a,y) # Taktfrequenzfilter
	yh = np.abs(signal.hilbert(yf)) # Hilbert-Transformation
 
	ma = max(y)
	tol = 0.002 * ma
	s = ma/5
	idx = np.where(np.abs(y-s)<tol)
	dt = np.fmod(idx, (tl*RATE))
	dt = np.mean(dt)
 
	return dt