Wavelet-Domain-Gauß-Prozesse: Was ist die Kovarianz?

20

Ich habe Maraun et al . Gelesen, "Nichtstationäre Gauß-Prozesse im Wavelet-Bereich: Synthese, Abschätzung und signifikante Tests" (2007), die eine Klasse von nichtstationären GPs definieren, die durch Multiplikatoren im Wavelet-Bereich spezifiziert werden können. Eine Realisierung eines solchen GP ist: wobei weißes Rauschen ist, die kontinuierliche Wavelet-Transformation in Bezug auf Wavelet , ist der Multiplikator (ein bisschen wie ein Fourier-Koeffizient) mit der Skala und der Zeit , und ist die inverse Wavelet-Transformation mit dem Rekonstruktions-Wavelet .& eegr; ( t ) W g g m ( b , a ) a b M h h

s(t)=Mhm(b,a)Wgη(t),
η(t)Wggm(b,a)abMhh

Ein zentrales Ergebnis der Arbeit ist, dass, wenn sich die Multiplikatoren nur langsam ändern, die Realisierung selbst nur "schwach" von der tatsächlichen Wahl von und abhängt . Somit spezifiziert den Prozess. Sie erstellen einige aussagekräftige Tests, um anhand von Erkenntnissen auf die Wavelet-Multiplikatoren schließen zu können.g h m ( b , a )m(b,a)ghm(b,a)

Zwei Fragen:

1. Wie bewerten wir die Standard - GP Wahrscheinlichkeit , das ist ?p(D)=N(0,K)

Ich würde vermuten, dass wir effektiv eine Änderung der Koordinaten vornehmen, sodass wobei die Wavelets und die (diagonale?) Matrix der Wavelet-Koeffizienten . Sie verwenden jedoch eine nicht-orthonormale CWT, sodass ich nicht weiß, ob dies korrekt ist.W M m ( a , b )K1=WTM1WWMm(a,b)

2. Wie kann diese Wavelet-Domain-GP mit einer Real-Space-GP in Beziehung gesetzt werden ? Können wir konkret aus einen Realraumkern (nicht stationär) berechnen ?m ( a , b )km(a,b)

Zum Vergleich: Der Kern eines stationären Gaußschen Prozesses ist das Fourier-Dual seiner spektralen Dichte (Bochner-Theorem, siehe Rasmussen, Kapitel 4) - ein einfacher Weg, um zwischen einem Realraum-GP und einem Frequenzraum-GP zu wechseln. Hier frage ich, ob es eine solche Beziehung im Wavelet-Bereich gibt.

cgreen
quelle
Hast du damit irgendwohin gekommen? Ich bin nicht sicher, ob die Änderung der Variablen richtig ist, da dies widerspricht, wenn man sagt, dass ist nennt man den reproduzierenden kernel? Kg,h(bb/a,a/a)=Wg,h(bb/a)
TTC

Antworten:

0

Der Fahrvorgang, weißes Rauschen η (t), ist unabhängig von der Wahl der Basis. In einem CWT (im Gegensatz zum DWT-Springen in Oktaven) gibt es eine gewisse Redundanz, schmale Wellenbänder überlappen sich. Das "Merkmal", das auf Signifikanz getestet wird, ist eine Varianz (Leistung), die in einer engen Frequenz über eine kurze Zeit beobachtet wird. Dies hängt natürlich mathematisch von dem gewählten Wavelet ab, ist aber nicht sehr viel - eine schmalere Bandbreite kann sich langsamer ändernde Merkmale mit größerer Empfindlichkeit erkennen, eine größere Bandbreite reagiert schneller, hat jedoch einen rauschintensiveren Hintergrund und ist weniger spezifisch.

  1. Da dies den Wavelet-Raum misst, der über die Dauer des Wavelets integriert ist, ist die von Ihnen geschriebene Transformation für jeden "Zeitpunkt" gültig. Im Allgemeinen benötigt man Phaseninformationen, um die CWT zu invertieren. Marauns Test ist im Wesentlichen Chi-Quadrat in Kraft.

  2. Nein. Maraun hängt vom Signal-zu-Rauschen-Verhältnis in einem Frequenzband über einen Zeitbereich ab. Dies kann viele verschiedene Realisierungen im Rauschraum haben und ist phasenunabhängig. Es ist empfindlich gegenüber einem AR (1) -Signal im Wavelet-Bereich bei einer bestimmten Frequenz, dh einer zeitlich anhaltenden Oszillation. Die CWT-Domäne tendiert dazu, eine isolierte Spitze im Breitbandrauschen zu unterdrücken.

James Prichard
quelle