Una matrice asimmetrica di Bernoulli soddisfa il PIR?


9

Definire una matrice di rilevamento con con probabilità e con probabilità . Fa soddisfare la proprietà isometria ristretta ?A A i j = 0 p A i j = 1 / n×NAAij=0p 1-pAAij=1/n1pA

Per riferimento, al caso simmetrico viene data risposta nel seguente documento:

RG Baraniuk, MA Davenport, RA DeVore e MB Wakin, "Una semplice prova della proprietà isometria limitata per matrici casuali," Approssimazione costruttiva, 28 (3) pp. 253-263, dicembre 2008. ( pdf )


Questo può essere un puntatore: ieeexplore.ieee.org/xpls/abs_all.jsp?arnumber=5512379 (purtroppo è paywalled e non ne ho trovato una copia OA). Non conosco il documento in dettaglio, ma quello che posso vedere da una rapida occhiata è che non considerano un caso generale come chiedi; considerano p = 1/2. Inoltre, non so quanto siano approfonditi sul PIR di tali matrici.
Thomas Arildsen,

Questo potrebbe anche essere un suggerimento: rauhut.ins.uni-bonn.de/RauhutSlidesLinz.pdf (pagina 98). Sfortunatamente, sembra che quelle che lui chiama Bernoulli variabili casuali siano casuali +/- 1 - non 0/1 (io chiamerei questi Rademacher).
Thomas Arildsen,

2
Consentitemi di ripetere l'essenza di un commento che ho fatto sullo stesso post (ora eliminato) su stats.SE : Aiuterebbe a rendere questa domanda più precisa e indicare a cosa esattamente siete interessati e cosa state lottando per adattarvi. Il commento di @Thomas è rilevante; inoltre non sappiamo a quale grado (cioè ordine) di scarsità siete interessati. Anche se consideriamo le funzioni di Rademacher, la risposta è chiaramente no in alcun senso uniforme (in ), poiché è (o, sufficientemente vicino ) in modo che vi sia (un'alta probabilità di) una sottostruttura essendo tutte. (cont.)p 1pp1
cardinale il

2
pn(0,1)np pAij=(1p)/npp/n(1p)

AAA

Risposte:


1

Come altri hanno affermato nei commenti, la risposta è "No". La media diversa da zero della matrice impone che un vettore medio diverso da zero (diciamo tutti) avrà un guadagno sostanzialmente più elevato di un vettore casuale con media zero (diciamo uniformemente casuale + 1, -1).

Si consideri la norma quadrata di A volte si prevede che un vettore costante y sia n * (p * N) ^ 2. (iterazione delle aspettative)

La norma quadrata di A volte un vettore x disegnato uniformemente da (-1, + 1) dovrebbe essere n * (p * N). (calcolabile per somma delle varianze della distribuzione binomiale)

Le norme di xey sono le stesse, ma l'aspettativa delle norme trasformate differisce di un fattore di p * N - divergente con l'aumentare delle dimensioni.

Ecco il codice matlab per aiutare a dimostrare.

n=2000;
N=1000;
p=.9;
A=double(rand(n,N)<p); 
x=sign(randn(N,1)); 
y=ones(N,1);
Ex_normSqAx = n*(N*p);  % E[ squared norm of A times random signs ]
Ex_normSqAy = n*(N*p)^2; % E[ squared norm of A times constant vector ]
normSqAx = norm(A*x)^2;
normSqAy = norm(A*y)^2;
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.