my_rand.m

来自「用于实现非均匀采样FFT（NDFT）」· M 代码 · 共 20 行

function x=my_rand(M,h)

x=[];
while(length(x)<M)
  X=rand(M^2,1);
  for j=1:length(X)
    is_ok=1;	
    for l=1:length(x)  
      if((abs(X(j)-x(l))<h)||(abs(X(j)+1-x(l))<h)||(abs(X(j)-1-x(l))<h))
	is_ok=0;
      end	
    end
    if(is_ok)
      x=[x;X(j)];
    end
  end
end

x=sort(x(1:M))-0.5;