ntnet.m

来自「这是一个数格网论的MATLAB程序源码」· M 代码 · 共 15 行

pp=[3,7,11,11,13,17,17,19,23,23,29];
s=2;
gama=2*cos([1:s]*2*pi/pp(s)); gama=gama-floor(gama);
%gama=sqrt(pp(1:s));
%h=[1  13];
n=0;
while (n<1000) ,
n=n+50; 
gp=[1:n]'*gama-fix([1:n]'*gama);
plot(gp(:,1),gp(:,2),'.'); pause;
     plot(sqrt(-2*log(gp(:,1))).*cos(2*pi*gp(:,2)),...
        sqrt(-2*log(gp(:,1))).*sin(2*pi*gp(:,2)),'.');
     axis([-4 4 -4 4]); grid; pause;
end