shell插值：awk一行代码实现线性插值功能

awk一行代码实现线性插值功能

输入文件要求，第一列的数值为待插值数列，每一行必须大于等于上一行的数值。下面版本的输入文件共三列，后两列为插值数列对应的函数值。

#!/bin/bash
awk 'BEGIN{no=0}NR!=1{ 
    no=no+1;
    a[no]=$1; 
    b[no]=$2; 
    c[no]=$3;
    }END{
    print a[1],b[1],c[1]
    for (i=a[1]+0.01; i<a[no]; i+=0.01) {
        jless=0;
        jmore=0;
        for (ele=1; ele<=no; ele++) {
            if (a[ele]<=i) jless++
            if (a[ele]>=i) jmore++
        }; 
        if ((jless+jmore)==no)
            {j=jless
            w1=(a[j+1]-i)/(a[j+1]-a[j]); 
            w2=1.0-w1;} 
        else
            {j=jless; w1=1; w2=0}
        print i,w1*b[j]+w2*b[j+1],w1*c[j]+w2*c[j+1] 
    }
    print a[no],b[no],c[no]
}'  input.dat > output.dat