added gaussian.awk based on hamming script;

added pwr(x,p) function to sterling_approx;
minor clean up all around;

1 files changed, 26 insertions, 0 deletions

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

+        }

+    }

+

+}