conv.awk


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

#!/usr/bin/awk -f

### conv.awk
# [PoC] linear convolution (with hardcoded IR window).

BEGIN {
    ARGV[1] ? X = ARGV[1] : X = 0
    ARGV[2] ? OFMT = "%." ARGV[2] "g" : OFMT = "%g"
    input_size = split(X, X_arr)
    #H = "1.0 1.0 1.0"      # rect
    #H = "0.25 0.50 0.25"   # von Hann
    H = "0.23 0.54 0.23"   # Hamming
    window_size = split(H, H_arr)
    output_size = (input_size + window_size - 1)
    for (n=1; n <= output_size; n++) {
        Y_arr[n] = 0
        for (m=1; m <= input_size; m++) {
            if (n <= window_size) {
                Y_arr[n] = Y_arr[n] + H_arr[n-m+1]*X_arr[m]
                continue
            }
            if ((n > window_size) && (n <= input_size)) {
                Y_arr[n] = Y_arr[n] + H_arr[n-m+1]*X_arr[m]
                continue
            }
            if ((n > window_size) && (n > input_size)) {
                Y_arr[n] = Y_arr[n] + H_arr[n-m+1]*X_arr[m]
                continue
            }
            else {
                Y_arr[n] = Y_arr[n] + 0
                continue
            }
        }
        print Y_arr[n]
    }
}