Làm cách nào để chọn ngẫu nhiên điểm k từ N điểm trong MATLAB?

Question

tôi sử dụng mã này để tạo ra và vẽ N điểm:

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

Làm thế nào để chọn k điểm (với xác suất p=0.25) ra khỏi N điểm, sau đó tô màu các phần k điểm màu đỏ và để lại các điểm khác như *.

Answer 1

Có hai cách tiếp cận bạn có thể thực hiện. Giải pháp đầu tiên là chọn ngẫu nhiên các giá trị k từ các giá trị N, đảm bảo rằng bạn luôn luôn có k điểm được chọn. Giải pháp thứ hai là chọn giá trị ngẫu nhiên với mỗi giá trị có xác suất trung bình p được chọn, điều này có thể dẫn đến ít nhất là 0 hoặc nhiều nhất là N được chọn ngẫu nhiên.

• Chọn k từ N giá trị:

  Bạn có thể sử dụng chức năng RANDPERM để tạo ra một hoán vị ngẫu nhiên của các số nguyên 1 qua N, sau đó chọn k giá trị đầu tiên trong danh sách hoán vị và replot họ màu đỏ:

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

• giá trị Đón khách với một xác suất trung bình p:

  Bạn có thể sử dụng RAND chức năng để chọn một giá trị ngẫu nhiên 0-1 cho mỗi N giá trị của bạn, sau đó chọn những người có một giá trị ngẫu nhiên nhỏ hơn hoặc bằng với bạn xác suất trung bình p và replot chúng như đỏ:

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

Answer 2

Từ những gì tôi hiểu, đối với mỗi người trong số các điểm ngẫu nhiên N bạn muốn lật một đồng xu để quyết định xem có nên Sel vv hay không (nơi đồng xu có một xác suất p=0.25 thành công!)

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 

này kết thúc lên được tương đương với chỉ tạo p*N điểm ngẫu nhiên ở nơi đầu tiên (ít nhất bạn tiếp cận con số này như N phát triển lớn hơn) ...

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

bạn có thể kiểm tra rằng tuyên bố cuối cùng cho các giá trị khác nhau của N:

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