Distribuzione di probabilità per un'onda sinusoidale rumorosa


12

Sto cercando di calcolare analiticamente una distribuzione di probabilità dei punti di campionamento da una funzione oscillante quando c'è un errore di misurazione. Ho già calcolato la distribuzione di probabilità per la parte "senza rumore" (lo metterò alla fine), ma non riesco a capire come includere "rumore".

Stima numerica

Per essere più chiari, immagina che ci sia una funzione cui scegli casualmente i punti durante un singolo ciclo; se raccogli i punti in un istogramma otterrai qualcosa relativo alla distribuzione.y(x)=sin(x)

Senza rumore

Ad esempio, ecco il e l'istogramma corrispondentesin(x)

inserisci qui la descrizione dell'immagine

Con il rumore

Ora, se c'è qualche errore di misurazione, cambierà la forma dell'istogramma (e quindi penso alla distribuzione sottostante). Per esempio

inserisci qui la descrizione dell'immagine

Calcolo analitico

Quindi speriamo di averti convinto che ci sia qualche differenza tra i due, ora scriverò come ho calcolato il caso "senza rumore":

Senza rumore

y(x)=sin(x)

Quindi se i tempi in cui campioniamo sono distribuiti uniformemente, allora la distribuzione di probabilità per deve soddisfare:y

P(y)dy=dx2π

allora da allora

dxdy=ddy(arcsin(y))=11y2

e così

P(y)=12π1y2

che con un'adeguata normalizzazione si adatta all'istogramma generato nel caso "nessun rumore".

Con il rumore

y(x)

Risposte:


10

Dipende da come è strutturato il processo del rumore.

Y

XixYi|Xi=xisin(xi)Yg

ϵiN(0,σ2)f(ϵ)=12πσexp(ϵ22σ2)fg

fY+ϵ(z)=(fg)(z)=fY(y)fϵ(zy)dy=fY(zw)fϵ(w)dw

inserisci qui la descrizione dell'immagine

(questa convoluzione è stata fatta numericamente; non so quanto sia trattabile questo integrale in questo esempio, perché non l'ho provato.)


Cose meravigliose, mi mancava l'idea di "convoluzione", poiché i tuoi numeri mostrano che questo è perfetto. Devo solo tentare l'integrazione. Grazie
Greg

2
Potresti trovarlo intrattabile, ma di solito non è difficile approssimare numericamente il risultato.
Glen_b

0

Penso che l'espressione derivata per P (x) sia disattivata di un fattore due. Il tempo di campionamento distribuito uniformemente equivale a distribuire uniformemente la fase nell'intervallo -pi, pi. La funzione trigonometrica distribuisce probabilità nell'intervallo y {-1,1}. L'integrazione di P (y) in questo intervallo deve essere = 1, non 2 come ottenuta usando il proprio integrando sopra. Penso che P (y) = 1 / (pi Sqrt (1-y ^ 2))


Potrebbe essere, è per questo che ho affermato "con un'adeguata normalizzazione" poiché ero troppo pigro per pensarci in quel momento. Grazie.
Greg,
Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.