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
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
|
/*
Non-linear Dynamic Sound Generators
Lance Putnam 2004
lance@uwalumni.com
This is a set of iterative functions and differential equations that
are known to exhibit chaotic behavior. Internal calculations are
done with 64-bit words to ensure greater accuracy.
The name of the function is followed by one of N, L, or C. These
represent the interpolation method used between function iterations.
N -> None
L -> Linear
C -> Cubic
*/
ChaosGen : UGen {
}
// General Quadratic Map
QuadN : ChaosGen {
const <equation="x1 = a*x0^2 + b*x0 + c";
*ar { arg freq=22050, a=1, b= -1, c= -0.75, xi=0, mul=1, add=0;
^this.multiNew('audio', freq, a, b, c, xi).madd(mul, add)
}
}
QuadL : QuadN {}
QuadC : QuadN {}
// Cusp Map
CuspN : ChaosGen {
const <equation="x1 = a - b*sqrt(|x0|)";
*ar { arg freq=22050, a=1, b=1.9, xi=0, mul=1, add=0;
^this.multiNew('audio', freq, a, b, xi).madd(mul, add)
}
}
CuspL : CuspN {}
// Gingerbreadman Map
GbmanN : ChaosGen {
const <equation="x1 = 1 - y0 + |x0|\ny1 = x0";
*ar { arg freq=22050, xi=1.2, yi=2.1, mul=1, add=0;
^this.multiNew('audio', freq, xi, yi).madd(mul, add)
}
}
GbmanL : GbmanN {}
// Henon Map
HenonN : ChaosGen {
const <equation="x2 = 1 - a*(x1^2) + b*x0";
*ar { arg freq=22050, a=1.4, b=0.3, x0=0, x1=0, mul=1.0, add=0.0;
^this.multiNew('audio', freq, a, b, x0, x1).madd(mul, add)
}
}
HenonL : HenonN {}
HenonC : HenonN {}
// Latoocarfian
LatoocarfianN : ChaosGen {
const <equation="x1 = sin(b*y0) + c*sin(b*x0)\ny1 = sin(a*x0) + d*sin(a*y0)";
*ar { arg freq=22050, a=1, b=3, c=0.5, d=0.5, xi=0.5, yi=0.5, mul=1.0, add=0.0;
^this.multiNew('audio', freq, a, b, c, d, xi, yi).madd(mul, add)
}
}
LatoocarfianL : LatoocarfianN {}
LatoocarfianC : LatoocarfianN {}
// Linear Congruential
LinCongN : ChaosGen {
const <equation="x1 = ((a * x0) + c) % m";
*ar { arg freq=22050, a=1.1, c=0.13, m=1.0, xi=0, mul=1.0, add=0.0;
^this.multiNew('audio', freq, a, c, m, xi).madd(mul, add)
}
}
LinCongL : LinCongN {}
LinCongC : LinCongN {}
// Standard Map
StandardN : ChaosGen {
const <equation="x1 = (x0 + y1) % 2pi\ny1 = (y0 + k*sin(x0)) % 2pi";
*ar { arg freq=22050, k=1.0, xi=0.5, yi=0, mul=1.0, add=0.0;
^this.multiNew('audio', freq, k, xi, yi).madd(mul, add)
}
}
StandardL : StandardN {}
// Feedback Sine with Linear Congruential Phase Indexing
FBSineN : ChaosGen {
const <equation="x1 = sin(im*y0 + fb*x0)\ny1 = (a*y0 + c) % 2pi";
*ar { arg freq=22050, im=1, fb=0.1, a=1.1, c=0.5, xi=0.1, yi=0.1, mul=1, add=0;
^this.multiNew('audio',freq,im,fb,a,c,xi,yi).madd(mul, add)
}
}
FBSineL : FBSineN {}
FBSineC : FBSineN {}
// ODEs
// 'h' is integration time-step
// Lorenz Attractor
LorenzL : ChaosGen {
const <equation="x' = s*(y-x)\ny' = x*(r-z)-y\nz' = x*y - b*z";
*ar { arg freq=22050, s=10, r=28, b=2.667, h=0.05, xi=0.1, yi=0, zi=0, mul=1.0, add=0.0;
^this.multiNew('audio', freq, s, r, b, h, xi, yi, zi).madd(mul, add)
}
}
|
