Python(三)——python实现斐波那契网格采样

python实现斐波那契网格采样

1.公式

具体参考我的另一篇博客：https://blog.csdn.net/m0_45866718/article/details/109461637

2.Code

作为python新手，没用解析式写

from mpl_toolkits import mplot3d
import numpy as np
import math
import matplotlib.pyplot as plt
%matplotlib inline
ax=plt.axes(projection="3d")

N=401
z=[]
x=[] 
y=[]
phi=(math.sqrt(5)+1)/2-1

for i in range(1,N+1):
    if(True):
        z.append((2*i-1)/N-1)
zline=np.array(z)


for i in range(1,N+1):
    if(True):
        x.append(math.sqrt(1-zline[i-1]*zline[i-1])*math.cos(2*math.pi*i*phi))
xline=np.array(x)

for i in range(1,N+1):
    if(True):
        y.append(math.sqrt(1-zline[i-1]*zline[i-1])*math.sin(2*math.pi*i*phi))
yline=np.array(y)

ax.scatter3D(2*xline,2*yline,2*zline,s=3)
plt.savefig("D:\\samples\\斐波那契.png")

3.结果