diff options
| author | wukong <wukong@longaeva> | 2018-11-04 10:59:44 -0800 |
|---|---|---|
| committer | wukong <wukong@longaeva> | 2018-11-04 10:59:44 -0800 |
| commit | e14342e827e2b42e43c006df90c7ad99e5124b3c (patch) | |
| tree | c73e930f3b111bb659957df1509e7d9f5a3c8493 /gaussian.awk | |
| parent | 72fdb25210c579beaabc35cae7ec803436887f20 (diff) | |
added gaussian.awk based on hamming script;
added pwr(x,p) function to sterling_approx;
minor clean up all around;
Diffstat (limited to 'gaussian.awk')
| -rw-r--r-- | gaussian.awk | 26 |
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
+ }
+ }
+
+}
|