From: kub <derkub@gmail.com>
Date: Wed, 26 Mar 2025 06:23:55 +0000 (+0100)
Subject: core, make dc filter reacting faster
X-Git-Tag: v2.04~12
X-Git-Url: https://notaz.gp2x.de/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=7b631dbe47856b35cc9e1a704110a98f37620abc;p=picodrive.git

core, make dc filter reacting faster

behaves closer to real hw according to measurements by notaz
---

diff --git a/pico/sound/mix.c b/pico/sound/mix.c
index 1a4785d1..430e6bc7 100644
--- a/pico/sound/mix.c
+++ b/pico/sound/mix.c
@@ -36,22 +36,22 @@ static struct iir {
 
 
 // exponential moving average combined DC filter and lowpass filter
-// y0[n] = (x[n]-y0[n-1])*alpha+y0[n-1], y1[n] = (y0[n] - y1[n-1])*(1-1/8192)
+// y0[n] = y0[n-1]+(x[n]-y0[n-1])*alpha, y1[n] = y[n-1]+(y0[n]-y1[n-1])*(1/512)
 static inline int filter_band(struct iir *fi2, int x)
 {
 	// low pass. alpha is Q8 to avoid loss by 32 bit overflow.
 //	fi2->y[0] += ((x<<(QB-8)) - (fi2->y[0]>>8)) * fi2->alpha;
 	fi2->y[0] += (x - (fi2->y[0]>>QB)) * fi2->alpha;
-	// DC filter. for alpha=1-1/8192 cutoff ~1HZ, for 1-1/1024 ~7Hz
-	fi2->y[1] += (fi2->y[0] - fi2->y[1]) >> QB;
+	// DC filter. for alpha=1/8192 cutoff ~1HZ, for 1/512 ~14Hz
+	fi2->y[1] += (fi2->y[0] - fi2->y[1]) >> 9;
 	return (fi2->y[0] - fi2->y[1]) >> QB;
 }
 
 // exponential moving average filter for DC filtering
-// y[n] = (x[n]-y[n-1])*(1-1/8192) (corner approx. 1Hz, gain 1)
+// y[n]= y[n-1] + (x[n]-y[n-1])*(1/512) (corner approx. 14Hz, gain 1)
 static inline int filter_exp(struct iir *fi2, int x)
 {
-	fi2->y[1] += ((x << QB) - fi2->y[1]) >> QB;
+	fi2->y[1] += ((x << QB) - fi2->y[1]) >> 9;
 	return x - (fi2->y[1] >> QB);
 }