1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
|
// https://youtu.be/tOx5UI8GGns
$fn = 8;
function bezier_3(p0, p1, p2, t) =
(( 1 - t) ^ 2) * p0
+ 2 * (1 - t) * t * p1
+ (t ^ 2) * p2;
function bezier_4(p0, p1, p2, p3, t) =
(( 1 - t) ^ 3) * p0
+ 3 * ((1 - t) ^ 2) * t * p1
+ 3 * (1 - t) * (t ^ 2) * p2
+ (t ^ 3) * p3;
function sum(array, index = 0) =
index < len(array) - 1 ?
array[index] + sum(array, index + 1) :
array[index] ;
function factorial(x) =
x > 1 ?
x * factorial(x - 1) :
1 ;
function binomial_coefficient(degree, index) =
factorial(degree) / (factorial(index) * factorial(degree - index)) ;
function bezier_p0(point, degree, index, t) =
binomial_coefficient(degree, index) * ((1 - t) ^ (degree - index)) * (t ^ index) * point;
function bezier_p(points, t) =
[for (i=[1:len(points[0])]) sum(bezier_p0(points[i], len(points) - 1, i, t))];
function bezier_p_debug(points, t) =
[for (i=[1:len(points[0])]) sum(bezier_p0(points[i], len(points) - 1, i, t))];
p0 = [0, 0];
p1 = [10, 0];
pb = [1, 5];
pn = [0, 10];
points = $fn;
points_list = [for (i=[0:points]) bezier_3(p0, p1, pn, i/points)];
//rotate_extrude(angle = 360)
polygon(points_list);
points_list_2 = [for (i=[0:points]) bezier_4(p0, p1, pb, pn, i/points)];
translate([20, 0, 0])
//rotate_extrude(angle = 360)
polygon(points_list_2);
p_set_1 = [p0, p1, pn];
p_set_2 = [p0, p1, pb, pn];
points_list_p1 = [for (i=[0:$fn]) bezier_p(p_set_1, i/$fn)];
points_list_p2 = [for (i=[0:$fn]) bezier_p(p_set_2, i/$fn)];
points_list_p2d = [for (i=[0:$fn]) bezier_p_debug(p_set_2, i/$fn)];
for (i=[0:$fn]) {
echo("i", i, bezier_p(p_set_2, i/$fn));
echo("i", i, bezier_p_debug(p_set_2, i/$fn));
for (j=[1:len(p_set_2[0])]) {
echo("j", j, sum(bezier_p0(p_set_2[j], len(p_set_2[j]), j, i/$fn)));
echo("j", j, bezier_p0(p_set_2[j], len(p_set_2[j]), j, i/$fn));
}
}
//echo(points_list);
//echo(points_list_p1);
echo(points_list_2);
echo(points_list_p2d);
echo(points_list_p2);
translate([0, 20, 0])
//rotate_extrude(angle = 360)
polygon(points_list_p1);
translate([20, 20, 0])
//rotate_extrude(angle = 360)
polygon(points_list_p2d);
|
