From e14342e827e2b42e43c006df90c7ad99e5124b3c Mon Sep 17 00:00:00 2001
From: wukong <wukong@longaeva>
Date: Sun, 4 Nov 2018 10:59:44 -0800
Subject: added gaussian.awk based on hamming script; added pwr(x,p) function
 to sterling_approx; minor clean up all around;

---
 gaussian.awk | 26 ++++++++++++++++++++++++++
 1 file changed, 26 insertions(+)
 create mode 100644 gaussian.awk

(limited to 'gaussian.awk')

diff --git a/gaussian.awk b/gaussian.awk
new file mode 100644
index 0000000..ef81b68
--- /dev/null
+++ b/gaussian.awk
@@ -0,0 +1,26 @@
+#!/usr/bin/awk -f
+
+### gaussian.awk
+# generate a Gaussian window
+# https://en.wikipedia.org/wiki/Window_function
+
+BEGIN {
+
+    N = ARGV[1]
+    M = 0.5*(N - 1)
+
+    # sigma <= 0.5
+    sigma = 0.4
+
+    for (n=0; n<N; n++) {
+        if (N > 1 && M > 0) {
+            e[n] = (n - M)/(sigma*M)
+            w[n] = exp(-0.5*e[n]*e[n])
+            print n, w[n]/M
+        }
+        else {
+            print n, 1.0
+        }
+    }
+
+}
-- 
cgit v1.2.3