计算awk中每行的标准差

# 计算每行数据的标准差
{
    # 将每行的数据存储到数组中
    for (i = 1; i <= NF; i++) {
        data[NR, i] = $i
    }
    # 更新每列的和
    for (i = 1; i <= NF; i++) {
        sum[i] += $i
    }
}
END {
    # 计算每列的平均值
    for (i = 1; i <= NF; i++) {
        avg[i] = sum[i] / NR
    }
    # 计算每列的标准差
    for (i = 1; i <= NF; i++) {
        for (j = 1; j <= NR; j++) {
            diff = data[j, i] - avg[i]
            squareSum[i] += diff * diff
        }
        stdDev[i] = sqrt(squareSum[i] / NR)
    }
    # 输出每列的标准差
    for (i = 1; i <= NF; i++) {
        printf("第 %d 列的标准差为：%f\n", i, stdDev[i])
    }
}

awk -f calculate_stddev.awk input.txt