求一组数的偏差和累加

wpd0508

已知文本1中有1 3 3 5 4 2 这组数 文本2存放着他们的平均数3 然后通过SHELL脚本如何实现求出各数的偏差和累加。X1=$1-3 ,X2=X1+$2-3......

yestreenstars

酱紫？

1. awk 'BEGIN{getline t < "b"}{for(i=0;i++<NF;)printf "X%d=%d\n",i,X[i]=X[i-1]+$i-t}' a
   
2. X1=-2
   
3. X2=-2
   
4. X3=-2
   
5. X4=0
   
6. X5=1
   
7. X6=0

复制代码

wpd0508

大大厉害，谢了，顺便问一下其中的a,b 是指的文件名吗？
回复 2# yestreenstars


   

wpd0508

如果我想要继续对得到的数据进行处理，比如说对以上得到的数组 [ Yi ]  以L个为单元进行分组， 即第一组为 [  Y1 Y2 Y3 ... YL].    对这Ｌ个数字进行最小二乘法拟合，得 Y=a1 x + a0, (这里的x为数组对应的下标数，即1～ L).
计算偏差量  E = sqrt( sum (Yi - (a1*xi +a0) )^2)这就很麻烦了是不是?能帮忙解决一下吗？
回复 2# yestreenstars


   

yestreenstars

回复 3# wpd0508
是的~{:2_172:}

   

yestreenstars

回复 4# wpd0508
没看懂~{:2_169:}

   

q1208c

这么复杂的计算, shell可能处理不了.

就算是能处理, 你以后维护的成本也会非常高.

建议使用 perl 或 python 吧.

yestreenstars

回复 7# q1208c
你是来抬杠的~{:2_172:}

   

wpd0508

但是我已经得到了YI这一组数，a1和a0是随便选的2个数，只需要进行最后一步的运算E = sqrt( sum (Yi - (a1*xi +a0) )^2)比如说E = sqrt( sum (Yi - (1*xi +2) )^2) 这也不能实现吗，真的很需要解决这个问题，求大大们点拨一下。

wpd0508

a0 = （∑Yi) / n - a1（∑Xi) / n （式1-
a1 = [n∑Xi Yi - （∑Xi ∑Yi)] / [n∑Xi2 - （∑Xi)2 )]   我错了，a0 和a1居然是2个值，我去他妹的。