数学图形(2.6)Satellit curve
这曲线有点像鼓,绕在球上两头是开口的.
#http://www.mathcurve.com/courbes3d/satellite/satellite.shtml vertices = 12000 t = from 0 to (40*PI) r = 10 k = rand2(0.5, 10) a = rand2(PI*0.1, PI*1.9) x = r*(cos(a)*cos(t)*cos(k*t) - sin(t)*sin(k*t)) y = r*(cos(a)*sin(t)*cos(k*t) + cos(t)*sin(k*t)) z = r*sin(a)*cos(k*t)
(1)当上面代码中的k == 1时
#http://www.mathcurve.com/courbes3d/satellite/satellite.shtml vertices = 12000 t = from 0 to (2*PI) r = 10 k = 1 a = rand2(PI*0.1, PI*1.9) x = r*(cos(a)*cos(t)*cos(k*t) - sin(t)*sin(k*t)) y = r*(cos(a)*sin(t)*cos(k*t) + cos(t)*sin(k*t)) z = r*sin(a)*cos(k*t)
生成一个帖在球上的伯努利双纽线
再将代码中的a = rand2(PI*0.1, PI*1.9)改为一个输入维度数据
vertices = D1:360 D2:100 u = from 0 to (PI) D1 v = from 0 to (2*PI) D2 r = 10 k = 1 x = r*(cos(v)*cos(u)*cos(k*u) - sin(u)*sin(k*u)) y = r*(cos(v)*sin(u)*cos(k*u) + cos(u)*sin(k*u)) z = r*sin(v)*cos(k*t)
这时生成一个曲面:
在这个曲面上,可以显示任意一个a值下生成的曲线.
(2)当上面代码中的k == 1时
2000 t = from 0 to (4*PI) r = 10 k = 0.5 a = rand2(PI*0.1, PI*1.9) x = r*(cos(a)*cos(t)*cos(k*t) - sin(t)*sin(k*t)) y = r*(cos(a)*sin(t)*cos(k*t) + cos(t)*sin(k*t)) z = r*sin(a)*cos(k*t)
再将代码中的a = rand2(PI*0.1, PI*1.9)改为一个输入维度数据
vertices = D1:360 D2:100 u = from 0 to (2*PI) D1 v = from 0 to (2*PI) D2 r = 10 k = 0.5 x = r*(cos(v)*cos(u)*cos(k*u) - sin(u)*sin(k*u)) y = r*(cos(v)*sin(u)*cos(k*u) + cos(u)*sin(k*u)) z = r*sin(v)*cos(k*t)
这时生成一个曲面: