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
|
// One sample delay
// usage: (delay1 signal)
(define delay1 (
(- (rzero $1 -1) $1)
)
)
// First-order IIR lowpass filter
// usage: (lp signal cutoff)
(define lp (
(let #coef (exp (/ (* (neg twopi) $2) sr)))
rpole (* $1 (- 1 #coef)) #coef
)
)
// First-order IIR highpass filter
// usage: (hp signal cutoff)
(define hp (
(let #1 (exp (/ (* (neg twopi) $2) sr)))
(let #2 (* (+ #1 1) 0.5))
rpole (rzero (* $1 #2) #2) #1
)
)
// First-order lowshelf filter
// usage: (lshelf signal cutoff gain)
(define lshelf (
(let #piOverSr (/ pi sr))
(let #gain (sqrt (pow 10 (* $3 0.05))))
(let #freq (* $2 #piOverSr))
rpole (rzero $1 (- 1 (* #freq #gain))) (- 1 (/ #freq #gain))
)
)
// First-order peaking filter
// usage: (peak signal cf res)
(define peak (
(let #p (sqrt (max 0 (min $3 0.99))))
(let #angle (/ (* twopi $2) sr))
(cpole (complex $1 0) (complex (* #p (cos #angle)) (* #p (sin #angle))))
)
)
// Notch filter
// usage: (notch signal cf res)
(define notch (
(let #ny (/ sr 2))
(let #cf (max 0 (min $2 #ny)))
(let #res (+ (* (max 0 (min $3 1)) 0.5) 0.5))
(let #a (+ (- 0.5 (/ #cf sr)) 0.5))
(let #a1 (cos (* #a twopi)))
(let #a2 (cos (* (- #a 0.25) twopi)))
(let #a3 (* #a2 -1))
(let #b1 (* #a1 #res))
(let #b2 (* #a3 #res))
(let #b3 (* #a2 #res))
(cpole
(cpole
(czero
(czero
(complex $1 0)
(complex #a1 #a2))
(complex #a1 #a3))
(complex #b1 #b3))
(complex #b1 #b2))
)
)
// Bandpass filter
// usage: (bp signal cf q)
(define bp (
(let #bw (/ $2 $3))
(let #b2 (pow e (* (/ (neg twopi) sr) #bw)))
(let #b1 (* (/ (* -4 #b2) (+ 1 #b2)) (cos (/ (* twopi $2) sr))))
(let #a1 (- 1 (sqrt #b2)))
(let #y2 (
(let #y1 (- (- (- (* #a1 $1) (* #a1 (delay (delay $1)))) (* #b1 #y1)) (* #b2 #y2)))
(delay #y1)
)
)
#y1
)
)
// DCblock filter
// usage (dcblock signal)
(define dcblock (
(rpole (rzero $1 1) 0.995)
)
)
|
