Wir dachtenHier Wie erstelle ich ein Hamming-64-Überlappungsfenster?

h = hamming(64);
h2 = hamming(38);
h = conv(h, h2);

Wir überlegen nun, wie Sie diese Fensterfunktion auf die resultierenden Variablen der Wigner-Ville-Verteilungsfunktion von Auger et al. In anwenden könnenZeit-Frequenz-Toolbox. Die Funktiontfrwv.m hat keinen Parameter für die Fensterfunktion.

Wir haben also diese Variablen

[B,T,F] = tfrwv(data, 1:length(data), length(data));

Hier ist eine Antwort auf ein verwandtes Problem, aber nicht ganz dasselbe. Man sagt, dass man die Fensterfunktion auf die Ergebnisse anwendet

Einfach Punkt für Punkt multiplizieren

Die Abmessungen vonh sind 101x1 doppelt, währendT undF 5001 x 1 doppelt. Daher scheint eine Extrapolation des Fenstervektors erforderlich zu sein, wenn Punkt für Punkt multipliziert wird.

Noch eine ErklärungHier

Ungefähr zur Hälfte des zweiten Codeblocks wende ich eine Fensterfunktion auf ein gepuffertes Signal an. Dies ist effektiv eine Vektormultiplikation der Fensterfunktion mit jedem gepufferten Block von Zeitreihendaten. Ich benutze einfach einen hinterhältigen diagonalen Matrixtrick, um es effizient zu machen.

Wie können Sie eine Fensterfunktion auf die Variablen anwenden?B, T, undF?