Add octave scripts for examples into track list.
[Faustine.git] / interpretor / preprocessor / faust-0.9.47mr3 / examples / faust-stk / NLFeks.dsp
1 declare name "Nonlinear EKS";
2 declare author "Julius Smith and Romain Michon";
3 declare version "1.0";
4 declare license "STK-4.3";
5 declare copyright "Julius Smith";
6 declare reference "http://ccrma.stanford.edu/~jos/pasp/vegf.html";
7 // -> Virtual\_Electric\_Guitars\_Faust.html";
8
9 import("music.lib"); // Define SR, delay
10 import("instrument.lib");
11 import("effect.lib"); // stereopanner
12
13 //==================== GUI SPECIFICATION ================
14
15 // standard MIDI voice parameters:
16 // NOTE: The labels MUST be "freq", "gain", and "gate" for faust2pd
17 freq = nentry("freq", 440, 20, 7040, 1); // Hz
18 gain = nentry("gain", 1, 0, 10, 0.01); // 0 to 1
19 gate = button("gate"); // 0 or 1
20
21 // Additional parameters (MIDI "controllers"):
22
23 // Pick angle in [0,0.9]:
24 pickangle = 0.9 * hslider("pick_angle",0,0,0.9,0.1);
25
26 // Normalized pick-position in [0,0.5]:
27 beta = hslider("pick_position [midi: ctrl 0x81]", 0.13, 0.02, 0.5, 0.01);
28 // MIDI Control 0x81 often "highpass filter frequency"
29
30 // String decay time in seconds:
31 t60 = hslider("decaytime_T60", 4, 0, 10, 0.01); // -60db decay time (sec)
32
33 // Normalized brightness in [0,1]:
34 B = hslider("brightness [midi:ctrl 0x74]", 0.5, 0, 1, 0.01);// 0-1
35 // MIDI Controller 0x74 is often "brightness"
36 // (or VCF lowpass cutoff freq)
37
38 // Dynamic level specified as dB level desired at Nyquist limit:
39 L = hslider("dynamic_level", -10, -60, 0, 1) : db2linear;
40 // Note: A lively clavier is obtained by tying L to gain (MIDI velocity).
41
42 //Nonlinear filter parameters
43 typeModulation = nentry("v:Nonlinear Filter/typeMod",0,0,4,1);
44 nonLinearity = hslider("Nonlinearity",0,0,1,0.01) : smooth(0.999);
45 frequencyMod = hslider("freqMod",220,20,1000,0.1) : smooth(0.999);
46
47 //==================== SIGNAL PROCESSING ================
48
49 //----------------------- noiseburst -------------------------
50 // White noise burst (adapted from Faust's karplus.dsp example)
51 // Requires music.lib (for noise)
52 noiseburst(gate,P) = noise : *(gate : trigger(P))
53 with {
54 diffgtz(x) = (x-x') > 0;
55 decay(n,x) = x - (x>0)/n;
56 release(n) = + ~ decay(n);
57 trigger(n) = diffgtz : release(n) : > (0.0);
58 };
59
60 nlfOrder = 6;
61 P = SR/freq ; // fundamental period in samples
62 Pmax = 4096; // maximum P (for delay-line allocation)
63
64 ppdel = beta*P; // pick position delay
65 pickposfilter = ffcombfilter(Pmax,ppdel,-1); // defined in filter.lib
66
67 excitation = noiseburst(gate,P) : *(gain); // defined in signal.lib
68
69 rho = pow(0.001,1.0/(freq*t60)); // multiplies loop-gain
70
71 // Original EKS damping filter:
72 b1 = 0.5*B; b0 = 1.0-b1; // S and 1-S
73 dampingfilter1(x) = rho * ((b0 * x) + (b1 * x'));
74
75 // Linear phase FIR3 damping filter:
76 h0 = (1.0 + B)/2; h1 = (1.0 - B)/4;
77 dampingfilter2(x) = rho * (h0 * x' + h1*(x+x''));
78
79 loopfilter = dampingfilter2; // or dampingfilter1
80
81 filtered_excitation = excitation : smooth(pickangle)
82 : pickposfilter : levelfilter(L,freq); // see filter.lib
83
84 //nonlinear allpass filter (nonLinearModulator is declared in instrument.lib)
85 NLFM = nonLinearModulator(nonLinearity,1,freq,typeModulation,frequencyMod,nlfOrder);
86
87 //declared in instrument.lib
88 stereo = stereoizer(P);
89
90 stringloop = (+ : fdelay4(Pmax, P-2)) ~ (loopfilter : NLFM);
91
92 process = filtered_excitation : stringloop : stereo : instrReverb;