python 圆形分布

import math

a = [月中角度]
b = [每月例数]
n =总例数
i=0
fsin_all = 0
fcos_all = 0
for aa,bb in zip(a,b):
    bb = bb
    i+=1
    value = math.sin(aa/180.0*math.pi)*bb
    fcos = math.cos(aa/180.0*math.pi)*bb
    fsin_all +=value
    fcos_all +=fcos
    print('month',i,':  yuezhongjiaodu=',aa,'lishu=',bb,'fsin=',value,'fcos=',fcos)

print('fsin_all=',fsin_all,'fcos_all=',fcos_all)
print('y_mean=',fsin_all/n,'xmean=',fcos_all/n)
y=fsin_all/n
x=fcos_all/n
print('arfa=',math.atan(y/x)*180/math.pi+180)
arfa=math.atan(y/x)*180/math.pi+180
print('tianshu=',arfa*365/360.0)

r = math.sqrt(x**2+y**2)
print('r=',r)
s = 122.9548*math.sqrt(-math.log(r))
print('s=',s,'r=',r)

print('r0.05=',math.sqrt(3/(n+0.28)))

注意：122.95479628483989=180/math.pi*math.sqrt(2)*math.sqrt(math.log(10))

import tensorflow as tf

with tf.Session() as sess:
    b=0
    a = sess.run(tf.random_normal([28], mean=4.62, stddev=2.58))
    print('tf.random_normal([100000], mean=9.67, stddev=2.795)')
    print(a)