Come seleziono a caso k punti da N punti in MATLAB?

Question

Io uso questo codice per creare e tracciare N punti:

N=input('No. of Nodes:'); 
data = rand(N,2) % Randomly generated n no. of nodes 
x = data(:,1); 
y = data(:,2); 
plot(x,y,'*'); 

Come selezionare k punti (con probabilità p=0.25) di N punti, poi il colore quei k punti rosso e lasciare gli altri punti come *.

Answer 1

Ci sono due approcci che è possibile adottare. La prima soluzione consiste nel selezionare in modo casuale i valori k dai valori N, i quali assicureranno che l'utente abbia sempre puntipunti scelti. La seconda soluzione consiste nel selezionare i valori in modo casuale con ognuno con una probabilità media di scegliere p, che potrebbe comportare un minimo di 0 o fino a N scelti a caso.

• Picking k da N valori:

  È possibile utilizzare la funzione RANDPERM per creare una permutazione casuale di interi 1 attraverso N, poi scegliere i primi k valori nell'elenco permutato e li rappresenterà nuovamente come rosso:

  index = randperm(N); 
  plot(x(index(1:k)),y(index(1:k)),'r*'); 
  

• valori Prelievo con una probabilità media p:

  È possibile utilizzare la funzione di RAND di scegliere un valore casuale 0-1 per ciascuno dei tuoi N valori, quindi scegliere quelli con un valore casuale minore o uguale al tuo probabilità media p e li rappresenterà nuovamente il rosso:

  index = (rand(N,1) <= p); 
  plot(x(index),y(index),'r*'); 
  

Answer 2

da quello che ho capito, per ciascuno del punto casuale N si vuole lanciare una moneta per decidere se a SEL ect o no (in cui la moneta ha un p=0.25 probabilità di successo!)

data = rand(N,2);    %# generate random points 
index = (rand(N,1) <= p);  %# roll coins to pick with prob p 
data(~index, :) = [];   %# keep only selected points 

Questo finisce per essere equivalente a generare solo p*N punti casuali, in primo luogo (almeno ci si avvicina a questo numero come N cresce più grande) ...

data = rand(p*N, 2);   %# directly generate p*N number of points 

si può verificare che l'ultima affermazione per diversi valori di N:

fprintf('1st = %d \n', p*N) 
fprintf('2nd = %d \n', sum(rand(N,1) <= p))